“…in Ω, ∂ t is the time derivative in the Sobolev-Bochner sense, ∆ is the distributional Laplacian, B denotes the canonical embedding of L 2 (0, T ; L 2 (D)) into L 2 (0, T ; L 2 (Ω)) obtained from an extension by zero, and •, • denotes the dual pairing in H 1 0 (Ω). Then, using [3, Theorem 1.13, Equation (1.70)], [8,Theorem 2.3], and the lemma of Aubin-Lions, see [31,Theorem 10.12], it is easy to check that the variational inequality in (10) possesses a well-defined, weak-to-weak continuous solution map G : L 2 (0, T ; L 2 (D)) → H 1 (0, T ; L 2 (Ω)), u → y. (Note that, in order to apply [3,Theorem 1.13], one has to define the function Φ appearing in this theorem as in [3,Equation (4.9)]. )…”