  We show that the notions of weak solution to the total variation flow based
on the Anzellotti pairing and the variational inequality coincide under some
restrictions on the boundary data. The key ingredient in the argument is a
duality result for the total variation functional, which is based on an
approximation of the total variation by area-type functionals.
