The last property measures the distance between the continuous and the
discrete solutions through the discrete trilinear form.
The proof differs from that of [1]},
since here we are interested in the analysis of the time-dependent case.
For this reason, we also use some techniques from [2]}.
Further, we require some extra regularity on the exact velocity solution \(\mathbf {u}\)  to ().
