\section{Introduction}\label{intro}

We deal with Cauchy-Dirichlet problems for parabolic equations and systems  of the form
\begin{equation}\label{eqdirichlet}
\begin{cases}
u_t - {\rm div} (|\nabla u|^{p-2} \nabla u ) = f  & {\rm in}\,\, \Omega_T\\
 u =0  &
{\rm on}\,\, \partial \Omega \times (0, T) \,
\\ u(\cdot ,0)= \psi (\cdot) &
\hbox{in} \,\, \Omega\,.
\end{cases}
\end{equation}
\iffalse
{\color{yellow} or
\begin{equation}\label{eqneumann}
\begin{cases}
u_t - {\rm div} (|\nabla u|^{p-2} \nabla u ) = f  & {\rm in}\,\,\, \Omega \times (0, T)\\
 \frac{\partial u}{\partial \nu} =0  &
\hbox{on} \partial \Omega \times (0, T) \,
\\ u(x,0)= \psi (x) &
{\rm for}\,\,\, x\in \Omega
\end{cases}
\end{equation}}
\fi
 Here, {\color{black} $p>1$,} $\Omega$ is an open set  in ${{\R}^n}$, $n \geq 1$, with  finite Lebesgue measure $|\Omega|$,  and $T>0$. Moreover,
$$\Omega _T = \Omega \times (0, T),$$
 the functions $f: \Omega_T \to \mathbb R^N$ and $\psi : \Omega \to \mathbb R^N$, $N\geq 1$ are given, and   $u: \Omega \to \mathbb R^N$ is the unknown. According to usage, $\nabla u$ stands for the gradient of $u$ with respect to the space variables $x \in \Omega$, and $u_t$ for its  derivative in time $t \in (0,T)$.
\par We are concerned with global second-order  regularity properties, with respect to the variables $x$, of the solution $u$ to problem \eqref{eqdirichlet}. Our main results can be summarized as follows. Assume that
\begin{equation}\label{hpf}
f \in L^2(\Omega_T)\,,
\end{equation}
\begin{equation}\label{hppsi}
\psi \in W^{1,p}_0(\Omega)\,,
\end{equation}
and $\Omega$ satisfies suitable regularity conditions (if $n \geq 2$). Then
\begin{equation}\label{loose}
|\nabla u|^{p-2} \nabla u \in L^2((0,T); W^{1,2}(\Omega))\,,
\end{equation}
and the norm of $|\nabla u|^{p-2} \nabla u$ in $L^2((0,T); W^{1,2}(\Omega))$ is bounded by the norms of the data $f$ and $\psi$. Consequently, $u$ is actually a strong solution to problem \eqref{eqdirichlet}.
This
provides a    natural nonlinear counterpart of the classical $L^2((0,T); W^{2,2}(\Omega))$ regularity of solutions to Cauchy-Dirichlet problems for the heat equation \cite{Eid, Fr, LaSoUr, Li}.
\par The results of the present paper are  new even in the case of smooth domains $\Omega$.  They will however be established under minimal regularity assumptions on $\partial \Omega$. In particular, they hold in any convex bounded domain $\Omega$.
\par
Let us notice that, if $p\geq 2$, then assumptions \eqref{hpf} and \eqref{hppsi} ensure that $f$ and $\psi$ belong to  proper function spaces for a classical weak solution to problem \eqref{eqdirichlet} to be well defined. On the other hand, this is not guaranteed if $1<p<2$. Some specification is thus in order.
\par
When $N=1$, namely when dealing with a single equation, a generalized notion of solution can still be introduced to cover the whole range of exponents $p \in (1, \infty)$, and our regularity theory holds   for every such $p$. One kind of solution that fits the situation at hand can be defined as the limit of solutions to approximating problems involving smooth data \cite{BG, DallA, Pr}. Such a solution will be called approximable throughout this paper. Its existence, uniqueness and basic regularity under \eqref{hpf} and \eqref{hppsi} is established  in Theorem \ref{exist}.  Its second-order differentiability properties -- the central issue of our contribution --  are addressed  in Theorems \ref{secondcap} and \ref{secondconvex}.
\par When $N>1$, i.e. when systems are in question, we restrict our attention to the case $p\geq 2$. This is the subject of Theorems \ref{secondcapsys} and \ref{secondconvexsys}. A reason for this limitation on $p$ is the lack of a suitable  existence and uniqueness theory of approximable solutions for parabolic systems.
One obstacle for this gap is the failure of standard truncation methods for vector-valued functions. In fact, also the elliptic theory of approximable solutions  for systems is incomplete. This affects
 our approach, which makes critical use of parallel results for elliptic equations and systems  recently obtained in \cite{cmsecond} and \cite{cmsecondsys}.
\par To conclude this section, let us point out that, in spite of  huge developments of the regularity theory of nonlinear singular and degenarate parabolic problems, presented e.g. in the reference monographs \cite{Di, DiGiVe, Li, lions} and in recent papers including \cite{AKN, BaHa, BoDuMa, BoDuMi, DSSV, FS, KuM1, KuM2, Schw},
 information available in the literature about second-order regularity of solutions to nonlinear parabolic problems is still limited. A result in the spirit of \eqref{loose} can be found in  \cite{Beirao}, which yet  requires additional regularity of the datum $f$   and is restricted to values of $p$ smaller than and close to $2$.
More classical contributions instead concern differentiability properties of the nonlinear expression $|\nabla u|^{\frac{p-2}2} \nabla u$, which differes from
 that appearing in \eqref{loose}. Furthermore,  they tipically apply to local solutions, and request higher regularity of $f$.

\iffalse
the exponent $p$ appearing in \eqref{eqdirichlet} is allowed to range in the whole interval $(1,\infty)$.

Since assumptions \eqref{hpf} and \eqref{hppsi} do not ensure that $f$ and $\psi$ belong to the proper function spaces for a classical weak solution to problem \eqref{exist} to be well defined if $p<2$, a generalized notion of solution has to be adopted. One kind of solutions that fits this situation has been introduced in the literature as the limit of solutions to approximationg problems involving smooth data \cite{???}. Such a solution will be called approximable throughout this paper. Its existence, uniqueness and basic regularity under \eqref{hpf} and \eqref{hppsi} is established  in Theorem \ref{exist} before attacking the issue of second-order differentiability in Theorems \ref{secondcap} and \ref{secondconvex}.
\par When $N>1$, i.e. if systems are in questions, we restrict our attention to the case $p\geq 2$. This is the subject of Theorems \ref{secondcapsys} and \ref{secondconvexsys}. A reason for this limitation on $p$ is the lack of a theory of approximable solutions for parabolic systems,
one obstacle being the failure of standard truncation methods for vector-valued functions.
Critical use in our approach is made of parallel results for elliptic equations and systems that we have recently obtained in \cite{cmsecond} and \cite{cmsecondsys}.
\fi

\iffalse

$p\in (1, \infty)$.  Our assumptions on the data amount to equiring that
 $f \in L^2(\Omega\times (0, T))$ and $\psi \in W^{1,p}_0(\Omega)$.
\iffalse
{\color{yellow} Of course, the compatibility condition
\begin{equation}\label{0mean}
\int _\Omega f\, dx =0
\end{equation}
 has to be required when dealing with \eqref{eqneumann}. }
\fi
\\ The notation
$$\Omega _T = \Omega \times (0, T)$$
will also be adopted in what follows.

{\color{red} Comment on definitions of solutions depending on assumptions on $f$ and $\psi$. A standard definition and existence result for weak solution requires in addition that $f \in L^{p'} (0,T; W^{-1,p'}(\Omega))$ and $\psi \in L^2(\Omega)$.}

\fi

\section{Main results}\label{main}

Let $N \geq 1$ and let $1<p<\infty$. Assume that $f\in L^{p'}((0,T); W^{-1,p'}(\Omega))$ and $\psi \in L^2(\Omega)$, where $p'= \tfrac p{p-1}$, the H\"older's conjugate of $p$. Then there exists a unique weak solution $u$ to problem \eqref{eqdirichlet}, namely a function $u \in C([0,T]; L^2(\Omega))\cap L^p((0,T); W^{1,p}_0(\Omega))$  such that $u_t \in L^{p'}((0,T); W^{-1,p'}(\Omega))$, $u(\cdot , 0) = \psi (\cdot)$ and
\begin{equation}\label{lions}
\int _0^\tau\langle u_t, \phi\rangle \, dt + \int _0^\tau \int_\Omega |\nabla u|^{p-2} \nabla u \cdot \nabla \phi \,dx\, dt  = \int _0^\tau\langle f, \phi\rangle \, dt
\end{equation}
for every $\tau \in (0, T]$ and
 every  function $\phi \in L^p((0,T); W^{1,p}_0(\Omega))$. Here, $\langle \cdot , \cdot \rangle$ stands for the duality pairing between $W^{-1,p'}(\Omega)$ and $W^{1,p}_0(\Omega)$. This follows from classical monotonicity arguments -- see e.g. \cite[Chapter 2, Sections 1.1 and 1.5]{lions}.

\iffalse

\par
{\color{red}
A function $u \in C(0,T; L^2(\Omega))\cap L^p(0,T; W^{1,p}_0(\Omega))$ \todo{Has $C(0,T; L^2(\Omega))$ to be replaced by $L^\infty(0,T; L^2(\Omega))$? See \cite{FreSch}} is called a weak solution to problem \eqref{eqdirichlet} if
\begin{multline}\label{weakpar}
\int _\Omega u(x,\tau)\phi (x,\tau)\, dx  + \int _0^\tau\int _\Omega -u \phi _t\, dx \, dt + \int _0^\tau\int _\Omega |\nabla u|^{p-2} \nabla u \cdot \nabla \phi \,dx\, dt  \\
=  \int_\Omega \psi (x) \phi(x, 0)\, dx  + \int _0^\tau\int _\Omega f \phi\, dx\, dt
\end{multline}
for  every $\tau \in (0, T]$ and every function $\phi \in W^{1,2}(0,T; L^2(\Omega))\cap  L^p(0,T; W^{1,p}_0(\Omega))$. If $f\in L^{p'}(0,T; (W^{1,p}_0(\Omega))')$ and {\color{red} $\psi \in L^2(\Omega)$}, then there exists a weak solution to problem \eqref{eqdirichlet} -- see \color{blue} \cite{FreSch}\color{black}.
\\  \color{red} The assumption $\cap L^\infty(\Omega \times (0,T))$ should not be needed, since $W^{1,2}(0,T; L^2(\Omega)) \subset L^2(\Omega _T)$, and hence the last integral in \eqref{weakpar} is well defined for every $\phi$ as above. }

\fi
\par As mentioned above, if the assumptions that $f \in  L^{p'}((0,T); W^{-1,p'}(\Omega))$ and $\psi \in L^2(\Omega)$ are dropped, then the notion of weak solution to problem \eqref{eqdirichlet} does not apply  anymore. Still, in the spirit of  \cite{BG, DallA, Pr}, generalized solutions in the approximable sense can be defined whenever $f$ and $\psi$ are merely integrable, namely if  $f\in L^1(\Omega_T)$ and $\psi \in L^1(\Omega)$.
\\
A function $u \in C([0,T]; L^1(\Omega)) \cap L^1((0,T); W^{1,1}_0(\Omega))$   will be called an approximable solution to problem \eqref{eqdirichlet}  if there exist  sequences $\{f_k\}\subset C_0^\infty(\Omega_T)$  and $\{\psi_k\}\subset C^\infty_0(\Omega)$  such that
$$_k \to f \quad \hbox{in $L^1(\Omega_T)$,} \quad \quad \psi_k \to \psi \quad \hbox{in $L^1(\Omega)$,}$$
and the sequence $\{u_k\}$ of
 weak solutions to the problems
\begin{equation}\label{eqdirichletk}
\begin{cases}
(u_k)_t - {\rm div} (|\nabla u_k|^{p-2} \nabla u_k ) = f_k  & {\rm in}\,\,\Omega_T\\
 u_k=0  &
{\rm on}\,\,
 \partial \Omega \times (0, T) \,
\\ u_k(\cdot ,0)= \psi_k (\cdot) &
{\rm in}\,\, \Omega\,
\end{cases}
\end{equation}
satisfies
\begin{equation}\label{june0}
u_k \to u \quad \hbox{and} \quad \nabla u_k \to \nabla u \quad \hbox{a.e. in $\Omega_T$,}
\end{equation}
and
\begin{equation}\label{june0bis}
u_k \to u \quad \hbox{in $ C([0,T]; L^1(\Omega))$.}
\end{equation}

Our first result ensures that, if $N=1$, problem \eqref{eqdirichlet} actually admits a unique approximable solution for every $p \in (1, \infty)$ which, under assumptions \eqref{hpf} and \eqref{hppsi}, enjoys additional regularity properties. Its uniqueness easily implies  that such a solution agrees with the weak solution if $p\geq 2$. This follows, for instance, by an argument as in the proof of Theorem \ref{secondcapsys}.

\begin{theorem}\label{exist} {\rm {\bf [Existence, uniqueness and basic regularity of approximable solutions]}}
Let $N=1$, $1<p<\infty$  and  $T>0$. Assume that $\Omega$ is an open set with finite measure in ${{\R}^n}$,  \  $n \geq 1$. Let $f \in
L^2(\Omega_T )$ and $\psi \in W^{1,p}_0(\Omega)$. Then
there exists a unique  approximable solution $u$ to problem \eqref{eqdirichlet}.  Moreover,
\begin{equation}\label{solamain1}
 u\in  L^\infty((0,T); W^{1,p}_0(\Omega)) \quad \hbox{and}    \quad u_t\in L^2(\Omega_T)\,,
\end{equation}
and there exists a constant $C=C(n, p, |\Omega|)$ such that
\begin{align}\label{solamain}
\| u\|_{L^\infty((0,T); W^{1,p}_0(\Omega))}^{\frac p2}+ \|u_t\|_{L^2(\Omega_T)}\leq C \big(\|f\|_{L^2(\Omega_T)} + \|\nabla \psi\|_{L^p(\Omega)}^{\frac p2}\big)\,.
\end{align}
\end{theorem}

Let us now turn to the primary objective of this paper, namely the global square integrability of weak derivatives  in $x$ of the  nonlinear expression of the gradient $|\nabla u|^{p-2} \nabla u$.
\\ We begin by introducing a few notions to be used in prescribing the regularity of $\Omega$ when $n \geq 2$.
A minimal regularity condition on $\Omega$ for our result to hold involves   a local isocapacitary inequality for
 the integral of the weak curvatures of $\partial \Omega$.  Specifically,
$\Omega$ is assumed to be a bounded Lipschitz domain. Moreover,  the functions of    $(n-1)$ variables that locally describe the boundary of $\Omega$ are  required to be twice weakly differentiable and integrable on $\partial \Omega$ with respect to the $(n-1)$-dimensional Hausdorff measure ${\cal H}^{n-1}$. This will be denoted by  $\partial \Omega \in W^{2,1}$. In particular, the weak second fundamental form
$\mathcal B$  on $\partial \Omega$ belongs to $L^1(\partial \Omega)$
These assumptions are not yet sufficient,
 as  demonstrated, for instance, by examples avaliable in the literature in the stationary  case -- see Remark \ref{sharp} below.
A local smallness condition on the $L^1$-norm of $\mathcal B$ is also needed. Specifically, denote
 by
$|\mathcal B|$ the norm of $\mathcal B$, and set
\begin{equation}\label{defK}
\mathcal K_\Omega(r) =
\displaystyle
\sup _{
\begin{tiny}
 \begin{array}{c}{E\subset \partial \Omega \cap B_r(x)} \\
x\in \partial \Omega
 \end{array}
  \end{tiny}
} \frac{\int _E |\mathcal B|d{\cal H}^{n-1}}{{\rm cap}_{B_1(x)} (E)}\qquad \hbox{for $r\in (0,1)$}\,.
\end{equation}
Here, $B_r(x)$ stands for the ball centered at $x$, with radius $r$,  and the notation ${\rm cap}_{B_1(x)}(E)$ is adopted for the capacity of the set $E$ relative to the ball $B_1(x)$.
Then we require that
\begin{equation}\label{capcond}
\lim _{r\to 0^+} \mathcal K _\Omega(r) < c_1
\end{equation}
for a suitable constant $c_1=c_1(n, N, p,d_\Omega, L_\Omega)$,
where,
  $d_\Omega$ and $L_\Omega$  denote the diameter
  and the   Lipschitz constant of $\Omega$, respectively. Here, and in similar occurrences in what follows, the dependence of a constant on $d_\Omega$  and  $L_\Omega$ is understood just via an upper bound for them.
\\ Recall that
the capacity ${\rm cap}_{B_1(x)}(E)$  of a set $E$ relative to $B_1(x)$ is defined as
\begin{equation}\label{cap}
{\rm cap}_{B_1(x)}(E) = \inf \bigg\{\int _{B_1(x)} |\nabla v|^2\, dy : v\in C^{0,1}_0(B_1(x)),\, v\geq 1 \,\hbox{on}\, E\bigg\}\,,
\end{equation}
where  $C^{0,1}_0(B_1(x))$ denotes the space of  Lipschitz continuous, compactly supported functions in $B_1(x)$.
Let us notice that,
if $n \geq 3$, then the  capacity ${\rm cap}_{B_1(x)}$  is equivalent (up to multiplicative constants depending on $n$) to the standard capacity in the whole of ${{\R}^n}$.

A more transparent condition on $\partial \Omega$, though slightly stronger  than  \eqref{capcond}, for our regularity result to hold can be given in terms of an integrability property of $\mathcal B$. It involves  the space of weak type, also called Marcinkiewicz space, with respect to the measure ${\cal H}^{n-1}$ on $\partial \Omega$ defined as
\begin{equation}\label{Xspace}
X= \begin{cases} L^{n-1, \infty} & \quad \hbox{if $n \ge 3$,}
\\
L^{1, \infty} \log L  & \quad \hbox{if $n =2$.}
\end{cases}
\end{equation}
Here,  $L^{n-1, \infty}$ denotes the
weak Lebesgue space endowed with the norm
\begin{equation}\label{weakleb}
\|h\|_{L^{q, \infty} (\partial \Omega)} = \sup _{s \in (0, {\cal H}^{n-1}(\partial \Omega))} s ^{\frac 1q} h^{**}(s),
\end{equation}
and $L^{1, \infty}\log L$ the weak Zygmund space endowed with the norm
\begin{equation}\label{weaklog}
\|h\|_{L^{1, \infty} \log L (\partial \Omega)} = \sup _{s \in (0, {\cal H}^{n-1}(\partial \Omega))} s \log\big(1+ \tfrac{C}s\big) h^{**}(s),
\end{equation}
for any constant  $C>{\cal H}^{n-1}(\partial \Omega)$. Observe that  different constants $C$ result in equivalent norms in \eqref{weaklog}.
 Here, $h ^{**}(s)= \smallint _0^s h^* (r)\, dr$ for $s >0$, where $h^*$ denotes the decreasing rearrangement of a measurable function $h: \partial \Omega \to \mathbb R$ with respect to ${\cal H}^{n-1}$.
\\ The relevant integrability property on $\mathcal B$ amounts to requiring that $\mathcal B \in X(\partial \Omega)$, an assumption that will be denoted by $\partial \Omega \in W^2X$.  An additional  smallness condition is however again needed on local norms of $\mathcal B$ in $X(\partial \Omega)$, which takes the form
\begin{equation}\label{smalln}
 \lim _{r\to 0^+} \Big(\sup _{x \in \partial \Omega} \|\mathcal B \|_{X(\partial \Omega \cap B_r(x))}\Big) < c_2 \,,
\end{equation}
for a suitable constant $c_2=c_2(n, N, p, L_\Omega, d_\Omega)$.
\par Observe that condition \eqref{smalln} is certainly  fulfilled if either $n\geq 3$ and $\partial \Omega \in W^{2,{n-1}}$, or $n=2$ and  $\partial \Omega \in W^{2, q}$ for some $q>1$,  and hence, in particular, if $\partial \Omega \in C^2$. Let us however emphasize that
 condition \eqref{smalln} does not even entail that $\partial \Omega \in C^1$.
 \par The link bewteen assumptions \eqref{smalln} and \eqref{capcond} is provided by
 \cite[Lemmas 3.5 and 3.7]{cmsecondsys}.  Those results  ensure that, given any constant $c_1$, there exists a constant $c_2=c_2(n, d_\Omega, L_\Omega, c_1)$ such that if $\Omega$ fulfills  \eqref{smalln}, then it also satisfies \eqref{capcond}.

\iffalse

{\color{orange}
\begin{theorem}\label{secondconvex} {\rm {\bf [Estimates in convex domains]}}
Let $p\in (1,\infty)$ and let $T>0$. Assume that $\Omega$ is any bounded  convex  open set in ${{\R}^n}$, with $n \geq 2$, let    $f \in
L^2(\Omega )$ and $\psi \in W^{1,p}_0(\Omega)$.
Let
$u$ be the generalized  solution   to     problem \eqref{eqdirichlet}
{\color{yellow}, or   the Neumann problem
\eqref{eqneumann}. }
Then
\begin{equation}\label{secondconvex1} u \in L^\infty(0, T; W^{1,p}_0(\Omega)), \quad |\nabla u|^{p-2} \nabla u
\in L^2(0,T; W^{1,2}(\Omega)), \quad u_t \in L^2(\Omega \times (0,T))\,.
\end{equation}
Moreover, there exists a constant $C=C(n,p, |\Omega|, T)$ such that
\begin{align}\label{secondconvex2}
\|\nabla u\|_{L^\infty(0,T; L^p(\Omega))}^{\frac p2}+ \||\nabla u|^{p-2} \nabla u\|_{L^2(0,T; W^{1,2}(\Omega))} + \|u_t\|_{L^2(\Omega \times (0,T))}\leq C \big(\|f\|_{L^2(\Omega_T)} + \|\nabla \psi\|_{L^p(\Omega)}^{\frac p2}\big)\,.
\end{align}
\end{theorem}
}

\fi

\begin{theorem}\label{secondcap} {\rm {\bf [Second-order estimates under minimal boundary regularity]}}
Let $N=1$, $1<p <\infty$ and  $T>0$.
Assume that
 $\Omega$ is a bounded  Lipschitz  domain in ${{\R}^n}$, $n \geq 2$, with $\partial \Omega \in W^{2,1}$.
Let $f \in
L^2(\Omega_T )$ and $\psi \in W^{1,p}_0(\Omega)$, and let $u$ be the approximable solution to problem \eqref{eqdirichlet}.
There exists a constant $c_1=c_1(n,p, d_\Omega, L_\Omega)$ such that, if  condition \eqref{capcond} is fulfilled, then
\begin{equation}\label{secondcap1} |\nabla u|^{p-2} \nabla u
\in L^2((0,T); W^{1,2}(\Omega))\,.
\end{equation}
Moreover, there exists a constant $C=C(n,p, \Omega, T)$ such that
\begin{align}\label{secondcap2}
 \||\nabla u|^{p-2} \nabla u\|_{L^2((0,T); W^{1,2}(\Omega))} \leq C \big(\|f\|_{L^2(\Omega_T)} + \|\nabla \psi\|_{L^p(\Omega)}^{\frac p2}\big)\,.
\end{align}
In particular, there exists a constant  $c_2=c_2(n,p, d_\Omega, L_\Omega)$ such that properties \eqref{secondcap1} and \eqref{secondcap2} hold if $\partial \Omega \in W^2X$ and fulfills condition \eqref{smalln}.
\end{theorem}

The sharpness of Theorem \ref{secondcap}, both in the  assumptions on $\Omega$ and in the estimate for $u$, is pointed out in the following remarks.

\begin{remark}\label{twosided}  {\rm {\bf [Sharpness of estimates]}}
{\rm
The  regularity for the derivative in $t$ of $u$ and the derivatives in $x$ of $|\nabla u|^{p-2} \nabla u$, provided by Theorems \ref{exist} and \ref{secondcap}, is sharp. Indeed, assume, for instance, that $\psi =0$. One trivially has that
\begin{align}\label{twosided1}
\|f\|_{L^2(\Omega_T)}  & = \|u_t - {\rm div} (|\nabla u|^{p-2} \nabla u )  \|_{L^2(\Omega_T)}
 \leq \|u_t \|_{L^2(\Omega_T)} +  \|{\rm div} (|\nabla u|^{p-2} \nabla u )  \|_{L^2(\Omega_T)}
\\ \nonumber & \leq
 \|u_t \|_{L^2(\Omega_T)} + c \||\nabla u|^{p-2} \nabla u\|_{L^2((0,T); W^{1,2}(\Omega))}
\end{align}
for some constant $c=c(n,p)$. Thus,  if $\Omega$ is as in Theorem \ref{secondcap}, then the following two-sided estimate holds:
\begin{align}\label{twosided2}
c_1 \|f\|_{L^2(\Omega_T)} \leq    \|u_t \|_{L^2(\Omega_T)} +  \||\nabla u|^{p-2} \nabla u\|_{L^2((0,T); W^{1,2}(\Omega))}
  \leq c_2  \|f\|_{L^2(\Omega_T)}
\end{align}
for suitable positive constants $c_1$ and $c_2$.
\iffalse

couple of estimates \eqref{solamain} and \eqref{secondcap2} is sharp.
Indeed, assume, for instance, that $\psi =0$ and that $\Omega$ is as in Theorem \ref{secondcap}. Then by Theorems \ref{exist} and \ref{secondcap} we deduce that
\begin{align}\label{twosided10}
c_1\|f\|_{L^2(\Omega_T)}   = c_1\|u_t - {\rm div} (|\nabla u|^{p-2} \nabla u )  \|_{L^2(\Omega_T)}
  \leq c_1( \|u_t \|_{L^2(\Omega_T)} +  \|{\rm div} (|\nabla u|^{p-2} \nabla u )  \|_{L^2(\Omega_T)})
 \leq
 \|u_t \|_{L^2(\Omega_T)} +  \||\nabla u|^{p-2} \nabla u\|_{L^2((0,T); W^{1,2}(\Omega))} \leq c_2\|f\|_{L^2(\Omega_T)}
\end{align}
for suitable positive constants $c_1$ and $c_2$.

 and that the solution $u$ to problem \eqref{eqdirichlet} satisfies properties  \eqref{solamain1} and \eqref{secondcap1}.  Then
\begin{align}\label{twosided1}
\|f\|_{L^2(\Omega_T)} & = \|u_t - {\rm div} (|\nabla u|^{p-2} \nabla u )  \|_{L^2(\Omega_T)}
 \\ \nonumber & \leq \|u_t \|_{L^2(\Omega_T)} +  \|{\rm div} (|\nabla u|^{p-2} \nabla u )  \|_{L^2(\Omega_T)}
\\ \nonumber & \leq
 \|u_t \|_{L^2(\Omega_T)} + c \||\nabla u|^{p-2} \nabla u\|_{L^2((0,T); W^{1,2}(\Omega))}
\end{align}
for some constant $c=c(n,p)$. Thus,  if $\Omega$ is as in Theorem \ref{secondcap}, then the following two-sided estimate holds:
\begin{multline}\label{twosided2}
c_1 \|f\|_{L^2(\Omega_T)} \\ \leq    \|u_t \|_{L^2(\Omega_T)} +  \||\nabla u|^{p-2} \nabla u\|_{L^2((0,T); W^{1,2}(\Omega))}
 \\   \leq c_2  \|f\|_{L^2(\Omega_T)}
\end{multline}
for suitable positive constants $c_1$ and $c_2$.
\fi
}
\end{remark}

\begin{remark}\label{sharp}  {\rm {\bf [Sharpness of boundary regularity]}}
{\rm
Membership of $\partial \Omega$ in $W^2X$, namely the mere finiteness of  the limit in \eqref{capcond}, is not sufficient for the conclusions of Theorem \ref{secondcap} to hold.
Actually,  in \cite{Ma67, Ma73} a family of domains $\{\Omega_\beta\}$ with the following properties is exhibited: (i) the limit in \eqref{capcond}, with $\Omega = \Omega_\beta$, is finite and depends continuously on $\beta$; (ii)
the stationary solution in $\Omega _\beta$ to the heat equation,
with a suitable smooth right-hand side,   belongs to $W^{2,2}(\Omega_T)$ if and only if the relevant limit  does not exceed an explicit threshold.
The boundary of each domain $\Omega_\beta$ is smooth outside a small region, where it agrees with the graph of a function $\Theta_\beta$ depending on the variables $(x_1, \dots , x_{n-1})$ only through $x_1$ and having the form
\begin{equation}\label{theta}
\Theta_\beta (x_1, \dots , x_{n-1})= \beta |x_1|(\log |x_1|)^{-1}
\end{equation}
  for small $x_1$. The limit  in  \eqref{capcond}   is a multiple of  $\beta$, and the  solution $u$ turns out to belong to $W^{2,2}(\Omega)$ if and only if the constant $\beta$ is smaller than or equal to an explicit value, depending only on $n$.
\\
The sharpness of the alternate
assumption \eqref{smalln} in Theorem \ref{secondcap}   can be shown, for instance, when $n=3$ and   $p\in (\tfrac 32,2]$, by an
 example from \cite{KrolM}, which again applies to the stationary case.  In that paper,  open sets  $\Omega \subset \mathbb R^3$ are constructed such that  the limit in \eqref{smalln} is finite,  but too large,  and  the stationary solution $u$ to problem \eqref{eqdirichlet}, with  a smooth right-hand side,  is so irregular that $|\nabla u|^{p-2}\nabla u \notin W^{1,2}(\Omega_T)$.   Similarly, in $\mathbb R^2$ there exist open sets $\Omega$
for which the limit in  \eqref{smalln}   is finite but larger than some critical value,
and  where the stationary solution $u$ to  the heat equation  with  a smooth right-hand side fails to belong to $W^{2,2}(\Omega)$ \cite{Ma67}.
\\ Let us notice that, if $\Omega_\beta$ is a domain as above, then  $\partial \Omega_\beta \notin   W^2L^{2, \infty}$ if $n \geq 3$.  Hence, if $\beta $ is sufficiently small, the capacitary criterion \eqref{capcond}  of Theorem \ref{secondcap} applies to deduce properties \eqref{secondcap1} and \eqref{secondcap2}, whereas the integrability condition \eqref{smalln} does not.
 }
\end{remark}

\begin{remark}\label{1-dim}  {\rm {\bf [Case $n=1$]}}
{\rm An inspection of the proof will reveal that, if $n=1$, then Theorem \ref{secondcap} holds  if $\Omega$ is any bounded interval. The result is considerarbly easier in this case, since the divergence operator agrees with plain differentiation in dimension one. }
\end{remark}

\iffalse

Capacities can be dispensed in the description of the  boundary regularity for the conclusions of Theorem \ref{eqdirichlet} to hold, and can be replaced with a
condition involving only integrability properties of $\mathcal B$. The relevant condition reads
\begin{equation}\label{smalln}
 \lim _{r\to 0^+} \Big(\sup _{x \in \partial \Omega} \|\mathcal B \|_{X(\partial \Omega \cap B_r(x))}\Big) < c \,,
\end{equation}
for a suitable constant $c=c(n, N, p, L_\Omega, d_\Omega)$,
where
\begin{equation}\label{Xspace}
X= \begin{cases} L^{n-1, \infty} & \quad \hbox{if $n \ge 3$,}
\\
L^{1, \infty} \log L  & \quad \hbox{if $n =2$.}
\end{cases}
\end{equation}
The spaces appearing on the right-hand side of \eqref{Xspace} are Marcinkiewicz type spaces, also called weak type space. Precisely,  $L^{n-1, \infty}$ denotes the
weak Lebesgue space $L^{n-1}$, and  $L^{1, \infty} \log L$   the weak Zygmund space $L\log L$. Let us point out that the assumption $\partial \Omega \in W^2X$ and condition \eqref{smalln} do not even entail that $\partial \Omega \in C^1$. Here, the notation $\partial \Omega \in W^2X$ means that the functions which  locally  describe $\partial \Omega$ are twice weakly differentiable, with second derivatives in  $X$.
\\ The following corollary follows from Theorem \ref{secondcap}, via \cite[Lemmas 3.5 and 3.7]{cmsecondsys}, which ensure that assumption \eqref{smalln} implies \eqref{capcond}.

\begin{corollary}\label{secondlorentz} {\rm {\bf [Regularity under integrability conditions on curvatures]}}
Let $p\in (1,\infty)$ and let $T>0$. Let $\Omega$ be a bounded  Lipschitz  domain in ${{\R}^n}$, $n \geq 2$, such that $\partial \Omega \in W^2X$. Let $f \in
L^2(\Omega )$ and $\psi \in W^{1,p}_0(\Omega)$, and let $u$ be the approximable solution to problem \eqref{eqdirichlet}. There eixts a constant $c=c(n,p, d_\Omega, L_\Omega)$ such that, if $\Omega$ satisfies condition \eqref{smalln}, then
$ |\nabla u|^{p-2} \nabla u
\in L^2(0,T; W^{1,2}(\Omega))$ and inequality \eqref{secondcap2}  holds for some positive constant  $C=C(n,p, \Omega, T)$.
\end{corollary}

\fi

Under the assumption that $\Omega$ is  a bounded convex open set, the conclusions of Theorem \ref{secondcap} hold without any additional regularity condition on $\partial \Omega$.

\begin{theorem}\label{secondconvex} {\rm {\bf [Second-order estimates in convex domains]}}
Let $N=1$, $1 <p <\infty$ and $T>0$. Assume that $\Omega$ is a bounded  convex open set in ${{\R}^n}$, $n \geq 2$. Let $f \in
L^2(\Omega_T )$ and $\psi \in W^{1,p}_0(\Omega)$, and let $u$ be the approximable solution to problem \eqref{eqdirichlet}. Then
$ |\nabla u|^{p-2} \nabla u
\in L^2((0,T); W^{1,2}(\Omega))$ and inequality \eqref{secondcap2}  holds for some positive constant  $C=C(n,p, \Omega, T)$.
\end{theorem}

Our results for systems are stated in Theorems \ref{secondcapsys} and \ref{secondconvexsys} below. They parallel Theorems  \ref{secondcap} and \ref{secondconvex}. The difference here is that they hold for $p \geq 2$ and can hence be stated in terms of weak solutions to problem \eqref{eqdirichlet}. Of course, Remarks \ref{twosided}--\ref{1-dim} carry over to this case.

\begin{theorem}\label{secondcapsys} {\rm {\bf [Second-order estimates for systems under minimal boundary regularity]}}
Let $N\geq 1$, $2 \leq p < \infty$ and  $T>0$.
Assume that
 $\Omega$ is a bounded  Lipschitz  domain with $\partial \Omega \in W^{2,1}$.
Let $f \in
L^2(\Omega_T )$ and $\psi \in W^{1,p}_0(\Omega)$, and let $u$ be the weak solution to problem \eqref{eqdirichlet}.
There exists a constant $c_1=c_1(n, N, p, d_\Omega, L_\Omega)$ such that, if  condition \eqref{capcond} is fulfilled, then
\begin{equation}\label{secondcapsys1} |\nabla u|^{p-2} \nabla u
\in L^2((0,T); W^{1,2}(\Omega))\,.
\end{equation}
Moreover, there exists a constant $C=C(n,N, p, \Omega, T)$ such that
\begin{align}\label{secondcapsys2}
 \||\nabla u|^{p-2} \nabla u\|_{L^2((0,T); W^{1,2}(\Omega))} \leq C \big(\|f\|_{L^2(\Omega_T)} + \|\nabla \psi\|_{L^p(\Omega)}^{\frac p2}\big)\,.
\end{align}
In particular, there exists a constant  $c_2=c_2(n, N, p, d_\Omega, L_\Omega)$ such that properties \eqref{secondcapsys1} and \eqref{secondcapsys2} hold if $\Omega \in W^2X$ and fulfils condition \eqref{smalln}.
\end{theorem}

\begin{theorem}\label{secondconvexsys} {\rm {\bf [Second-order estimates for systems in convex domains]}}
Let $N\geq 1$, $2\leq p <\infty$ and $T>0$. Assume that $\Omega$ is a bounded  convex open set in ${{\R}^n}$, $n \geq 2$. Let $f \in
L^2(\Omega_T)$ and $\psi \in W^{1,p}_0(\Omega)$, and let $u$ be the weak solution to problem \eqref{eqdirichlet}. Then
$ |\nabla u|^{p-2} \nabla u
\in L^2((0,T); W^{1,2}(\Omega))$ and inequality \eqref{secondcapsys2}  holds for some positive constant  $C=C(n, N, p, \Omega, T)$.
\end{theorem}

\section{Proofs}\label{proofs}

\iffalse
Let $N \geq 1$ and $1<p<\infty$. Assume that $f\in L^{p'}((0,T); W^{-1,p'}(\Omega))$ and $\psi \in L^2(\Omega)$. Then there exists a unique weak solution $u$ to problem \eqref{eqdirichlet}, namely a function $u \in C([0,T]; L^2(\Omega))\cap L^p((0,T); W^{1,p}_0(\Omega))$  such that $u_t \in L^{p'}((0,T); W^{-1,p'}(\Omega))$, $u(\cdot , 0) = \psi (\cdot)$ and
\begin{equation}\label{lions}
\int _0^\tau\langle u_t, \phi\rangle \, dt + \int _0^\tau \int_\Omega |\nabla u|^{p-2} \nabla u \cdot \nabla \phi \,dx\, dt  = \int _0^\tau\langle f, \phi\rangle \, dt \quad \hbox{for $\tau \in [0, T]$,}
\end{equation}
for every $\tau \in (0, T]$ and
 every  function $\phi \in L^p((0,T); W^{1,p}_0(\Omega))$. Here, $\langle \cdot , \cdot \rangle$ stands for the duality pairing between $W^{-1,p'}(\Omega)$ and $W^{1,p}_0(\Omega)$. This follows from classical monotonicity arguments -- see e.g. \cite[Chapter 2, Sections ????]{lions}.

\par
{\color{red}
A function $u \in C(0,T; L^2(\Omega))\cap L^p(0,T; W^{1,p}_0(\Omega))$ \todo{Has $C(0,T; L^2(\Omega))$ to be replaced by $L^\infty(0,T; L^2(\Omega))$? See \cite{FreSch}} is called a weak solution to problem \eqref{eqdirichlet} if
\begin{multline}\label{weakpar}
\int _\Omega u(x,\tau)\phi (x,\tau)\, dx  + \int _0^\tau\int _\Omega -u \phi _t\, dx \, dt + \int _0^\tau\int _\Omega |\nabla u|^{p-2} \nabla u \cdot \nabla \phi \,dx\, dt  \\
=  \int_\Omega \psi (x) \phi(x, 0)\, dx  + \int _0^\tau\int _\Omega f \phi\, dx\, dt
\end{multline}
for  every $\tau \in (0, T]$ and every function $\phi \in W^{1,2}(0,T; L^2(\Omega))\cap  L^p(0,T; W^{1,p}_0(\Omega))$. If $f\in L^{p'}(0,T; (W^{1,p}_0(\Omega))')$ and {\color{red} $\psi \in L^2(\Omega)$}, then there exists a weak solution to problem \eqref{eqdirichlet} -- see \color{blue} \cite{FreSch}\color{black}.
\\  \color{red} The assumption $\cap L^\infty(\Omega \times (0,T))$ should not be needed, since $W^{1,2}(0,T; L^2(\Omega)) \subset L^2(\Omega _T)$, and hence the last integral in \eqref{weakpar} is well defined for every $\phi$ as above. }

\medskip
\par
A function $u \in C(0,T; L^1(\Omega)) \cap L^1(0,T; W^{1,1}_0(\Omega))$ on $\Omega_T$  is called a generalized solution to problem \eqref{eqdirichlet} (in the approximable sense) if there exist  sequences $\{f_k\}\subset C_0^\infty(\Omega_T)$  and $\{\psi_k\}\subset C^\infty_0(\Omega)$  such that
\begin{equation}\label{june1}
f_k \to f \qquad \hbox{in $L^1(\Omega_T)$,}
\end{equation}
\begin{equation}\label{june2}
\psi_k \to \psi \qquad \hbox{in $L^1(\Omega)$,}
\end{equation}
and the sequence $\{u_k\}$ of
 weak solutions to the problems
\begin{equation}\label{eqdirichletk}
\begin{cases}
(u_k)_t - {\rm div} (|\nabla u_k|^{p-2} \nabla u_k ) = f_k  & {\rm in}\,\,\, \Omega \times (0, T)\\
 u_k(x,t) =0  &
{\rm for}\,\,\,
(x,t)\in \partial \Omega \times (0, T) \,
\\ u_k(x,0)= \psi_k (x) &
{\rm for}\,\,\, x\in \Omega\,
\end{cases}
\end{equation}
satisfies \todo{Check if this is a suitable definition of approximable solution}
\begin{equation}\label{june0}
u_k \to u \quad \hbox{and} \quad \nabla u_k \to \nabla u \quad \hbox{a.e. in $\Omega_T$,}
\end{equation}
and
\begin{equation}\label{june0bis}
u_k \to u \quad \hbox{in $ C(0,T; L^1(\Omega))$.}
\end{equation}
\fi

\iffalse
\bigskip
\par
Let $g \in L^1(\Omega)$. A function $v \in \mathcal T^{1,p}_0(\Omega)$ is called an approximable solution to the Dirichlet problem
\begin{equation}\label{L1}
\begin{cases}
- {\rm div} (|\nabla v|^{p-2} \nabla v ) = g  & {\rm in}\,\,\, \Omega \\
 u =0  &
{\rm on}\,\,\,
\partial \Omega
\end{cases}
\end{equation}
if there exist sequences $\{g_k\} \subset C^\infty_0(\Omega)$ and $v_k \in W^{1,p}_0(\Omega)$ such that $g_k \to g$ in $L^1(\Omega)$, $v_k$ is the weak solution to problem
\begin{equation}\label{L1k}
\begin{cases}
- {\rm div} (|\nabla v_k|^{p-2} \nabla v_k ) = g_k  & {\rm in}\,\,\, \Omega \\
 u_k =0  &
{\rm on}\,\,\,
\partial \Omega\,,
\end{cases}
\end{equation}
and
\begin{equation}\label{L1conv}
v_k \to v \quad \hbox{and} \quad \nabla v_k \to \nabla v \quad \hbox{a.e. in $\Omega$.}
\end{equation}
Here, $\nabla v$ denotes the generalized gradient of $v$ in the sense of \eqref{??}.

\begin{proposition}\label{sola} Let $\Omega$ be a bounded open subset of ${{\R}^n}$, and let $1<p<\infty$. Assume that $g \in L^1(\Omega)$. Let $v$ be the approximable solution of the Dirichlet problem \eqref{L1}. Assume that the function $w \in W^{1,p}_0(\Omega)$ fulfills
\begin{equation}\label{sola1}
\int_\Omega |\nabla w|^{p-2}\nabla w \cdot \nabla \varphi\, dx = \int _\Omega g \,\varphi\, dx
\end{equation}
for every $\varphi \in W^{1,p}_0(\Omega) \cap L^\infty (\Omega)$. Then
\begin{equation}\label{sola2}
w=v\,.
\end{equation}
\end{proposition}
\par\noindent
{\bf Proof}. Since   $w \in W^{1,p}_0(\Omega)$ and $v_k \in W^{1,p}_0(\Omega)$, one has that $T_t(w-v_k) \in W^{1,p}_0(\Omega) \cap L^\infty (\Omega)$ for every $k\in \mathbb N$ and $t>0$. Making use of $T_t(w-v_k)$ as test function $\varphi$ in equation \eqref{sola1} and in the definition of weak solution to problem \eqref{L1k}, and subtracting the resultant equations tells us that
\begin{equation}\label{sola3}
\int_{\{|w-v_k|< t\}}(|\nabla w|^{p-2}\nabla w - |\nabla v_k|^{p-2}\nabla v_k)\cdot (\nabla w - \nabla v_k)\, dx = \int _\Omega (g -g_k)\, T_t(w-v_k) \, dx\,.
\end{equation}
Since $|T_t(w-v_k)|\leq t$ and $g_k \to g$ in $L^1(\Omega)$, the right-hand side of equation \eqref{sola3} tends to $0$ as $k \to \infty$. Thus, inasmuch as  the integrand on the left-hand side of \eqref{sola3} is nonnegative,  passing to the limit in \eqref{sola3} as $k \to \infty$ yields, by Fatou's lemma and equation \eqref{L1conv},
\begin{equation}\label{sola4}
\int_{\{|w-v|< t\}}(|\nabla w|^{p-2}\nabla w - |\nabla v|^{p-2}\nabla v)\cdot (\nabla w - \nabla v)\, dx =0
\end{equation}
for every $t>0$. The integrand in \eqref{sola4} in nonnegative, and vanishes only if $\nabla w = \nabla v$. Therefore, $\nabla w = \nabla v$ a.e. in the set $\{|w-v|< t\}$, for every $t>0$, whence
\begin{equation}\label{sola5}
\nabla w = \nabla v \qquad \hbox{a.e. in $\Omega$.}
\end{equation}
Now, the function $T_\tau (w-T_t(v)) \in W^{1,p}_0(\Omega)\subset W^{1,1}_0(\Omega)$ for every $t, \tau >0$. An application of the Sobolev inequality and the use of equation \eqref{sola5} then imply that
\begin{equation}\label{sola6}
\bigg( \int_{\Omega }|T_{\tau}(w-T_t(v)|^{n'}\, dx \bigg
)^{\frac{1}{n'}}
\le C
\bigg(
\int_{\{t<|u|<t+\tau\} }|\nabla w|\, dx+
\int_{\{t-\tau<|u|<t\} }|\nabla w|\, dx
\bigg)
\end{equation}
for some constant $C=C(n)$ and for every $t, \tau >0$. Since $|\nabla w| \in L^1(\Omega)$, for each $\tau >0$, the right-hand side of inequality \eqref{sola6} converges to $0$ as $t\to \infty$. Therefore,  passing to the limit  in \eqref{sola6} as $t\to \infty$ tells us that
$$\int_{\Omega }|T_{\tau}(w-v)|^{n'}\, dx =0$$
for every $\tau >0$. Hence, equation \eqref{sola2} follows. \qed

\fi
\iffalse
{\color{yellow}
The next result is a version, global in time, of Proposition \cite[]{FS}. It provides us with an $L^2(\Omega_T)$-estimate for $u_t$, and an $L^\infty(0, T; W^{1,p}(\Omega))$-estimate for the solution $u$ to problem \eqref{eqdirichlet}. The proof is analogous to (in fact, simpler than) that of Proposition \cite[]{FS}, and will be omitted.

\begin{proposition}\label{ut} Let $\Omega$ be a bounded open subset of ${{\R}^n}$, and let $T>0$.  Assume that $f \in C^\infty _0 (\Omega \times (0,T))$ and
$\psi \in C^\infty _0(\Omega)$, and let $u$ be the weak solution to problem \eqref{eqdirichlet}. Then
\begin{equation}\label{ut2} u \in L^\infty(0, T; W^{1,p}_0(\Omega)), \quad  u_t \in L^2(\Omega \times (0,T))\,
\end{equation}
and there exists a constant $C=C(??)$ such that
\begin{equation}\label{ut1}
\int _0^\tau\int _{\Omega} u_t^2 \, dx \,dt + \sup _{\tau \in (0, T)}\int _\Omega |\nabla u(x, \tau)|^p\, dx \leq C \bigg(\int_0^T \int_{\Omega} f^2\, dx\, dt + \int _\Omega |\psi|^p \, dx\bigg)\,.
\end{equation}
\end{proposition}
}
\fi

Our proof of Theorem \ref{exist} combines arguments introduced in \cite{BBGGPV} in the elliptic case and developed in \cite{Pr} for parabolic problems, with an estimate for $u_t$ and $\nabla u$ in the spirit of  \cite[Proposition 4.1]{FS}, and of earlier results from \cite[Chapter 1, Theorem 8.1]{lions} and \cite{BEKP}.

\medskip
\par\noindent
{\bf Proof  of Theorem \ref{exist}}.  We assume that $n \geq 2$, the case when $n=1$ being analogous, and even simpler.
Let $\{f_k\} \subset C^\infty_0(\Omega_T)$ and $\psi_k \in C^{\infty}_0(\Omega)$ be sequences such that
\begin{equation}\label{aug1}
f_k \to f \quad \hbox{in $L^2(\Omega_T)$} \quad \hbox{and} \quad \|f_k\|_{L^2(\Omega_T)} \leq 2  \|f\|_{L^2(\Omega_T)}\,,
\end{equation}
\begin{equation}\label{aug1bis}
\psi_k \to \psi \quad \hbox{in $W^{1,p}(\Omega)$} \quad \hbox{and} \quad \|\nabla \psi_k\|_{L^{p}(\Omega)} \leq 2  \|\nabla \psi\|_{L^{p}(\Omega)}\,,
\end{equation}
and let $\{u_k\}$ be the corresponding sequence of weak solutions to problems \eqref{eqdirichletk}.
A global in time version of \cite[Proposition 4.1]{FS}  tells us that  $u_k \in L^\infty((0, T); W^{1,p}(\Omega))$, $(u_k)_t \in L^2(\Omega_T)$ and
\begin{align}\label{solamaink}
\| u_k\|_{L^\infty((0,T); W^{1,p}_0(\Omega))}^{\frac p2}+ \|(u_k)_t\|_{L^2(\Omega_T)}& \leq C \big(\|f_k\|_{L^2(\Omega_T)} + \|\nabla \psi_k\|_{L^p(\Omega)}^{\frac p2}\big) \\ \nonumber & \leq C'  \big(\|f\|_{L^2(\Omega_T)} + \|\nabla \psi\|_{L^p(\Omega)}^{\frac p2}\big)  \,,
\end{align}
for suitable constants $C$ and $C'$ depending on $n$, $p$, $|\Omega|$ and $T$.
Inequality \eqref{solamaink} holds in the whole interval $[0,T]$, namely up to $t=0$, instead of just locally, as in the result of \cite{FS},  thanks to the present assumption  that $\psi$ belongs  to $W^{1,p}_0(\Omega)$.  The proof  is completely analogous to (in fact, simpler than) that of \cite[Proposition 4.1]{FS},  and will be omitted.
\\ One clearly has that $L^\infty((0,T); W^{1,p}_0(\Omega)) \to L^p(\Omega_T)$.
Thus, owing to \eqref{solamaink}, the sequence $\{u_k\}$ is bounded in the anisotropic Sobolev space $W^{1,(p,2)}(\Omega_T)$ defined as
$$W^{1,(p,2)}(\Omega_T)=\{v:\, v\in L^p(\Omega_T), \,|\nabla v|\in L^p(\Omega_T), v_t \in L^2(\Omega_T)\}\,,$$
and equipped with the norm
$$\|v\|_{W^{1,(p,2)}(\Omega_T)} = \|v\|_{L^p(\Omega_T)} + \|\nabla v\|_{L^p(\Omega_T)}+ \|v_t\|_{L^2(\Omega_T)}\,.$$
Since $p \in (1,\infty)$, the space $W^{1,(p,2)}(\Omega_T)$ is reflexive. Moreover, on setting
$q=\min\{p, 2\}$, one trivially has $W^{1,(p,2)}(\Omega_T)\to W^{1,q}(\Omega_T)$, and hence the sequence $\{u_k\}$ is also bounded in $W^{1,q}(\Omega_T)$.
Therefore, given any number  $r \in \big[1, \tfrac{q(n+1)}{n+1-q}\big)$,
there exists
a subsequence, still denoted by $\{u_k\}$, such that
\begin{equation}\label{aug2}
u_k \rightharpoonup u \quad \hbox{weakly in $W^{1, (p,2)}(\Omega_T)$} \quad \hbox{and} \quad u\to u \quad \hbox{in $L^r(\Omega_T)$.}
\end{equation}
\\ Our goal is now to show that $\{\nabla u_k\}$ is a Cauchy sequence in measure. Fix any $\lambda >0$ and $\varepsilon>0$. Given any $\theta , \sigma >0$, one has that
\begin{align}\label{aug3}
|\{|\nabla u_k -\nabla u_m|>\lambda\}| & \leq |\{|\nabla u_k|>\theta\}| +|\{|\nabla u_m|>\theta\}| + |\{|u_k -u_m|>\sigma\}|\\ \nonumber & \quad + |\{|\nabla u_k -\nabla u_m|>\lambda, |u_k -u_m|\leq\sigma, |\nabla u_k|\leq \theta, |\nabla u_m|\leq\theta \}|
\end{align}
for $k,m \in \mathbb N$. Since the sequence $\{u_k\}$ is bounded in $W^{1, (p,2)}(\Omega_T)$, there exists $\theta >0$ such that
\begin{align}\label{aug4}
|\{|\nabla u_k|>\theta\}| +|\{|\nabla u_m|>\theta\}|< \varepsilon
\end{align}
for every $k, m \in \mathbb N$. Define, for $\sigma >0$, the function $T_\sigma : \mathbb R \to \mathbb R$ as
$$T_\sigma (s) = \begin{cases} s & \quad \hbox{if $|s|< \sigma$,}
\\ \sigma \,{\rm sign} (s) & \quad \hbox{if $|s|\geq \sigma$.}
\end{cases}
$$
Clearly, $|T_\sigma (s)|\leq \sigma $ for $s \in \mathbb R$.
 Choose the test function $\phi = T_\sigma (u_k-u_m)$ in the weak formulation of problem \eqref{eqdirichletk}, namely in equation \eqref{lions} with $u$, $f$, $\psi$ replaced by $u_k$, $f_k$, $\psi_k$. Next, choose the same test function in the same problem with $k$ replaced by $m$, and subtract the resultant equations to obtain that
\begin{align}\label{aug5}
\int _0^T   & \langle (u_k - u_m)_t ,T_\sigma (u_k - u_m) \rangle \,dt \\ \nonumber & \quad +  \int _0^T \int_\Omega (|\nabla u_k|^{p-2}\nabla u_k - |\nabla u_m|^{p-2}\nabla u_m)\cdot \nabla (T_\sigma (u_k - u_m))\, dx\,dt \\ \nonumber &
= \int _0^T \int_\Omega  (f_k - f_m)T_\sigma (u_k - u_m)\, dx\,dt \,.
\end{align}
 Define
the function $\Lambda_\sigma : \mathbb R \to \mathbb R$ as $\Lambda _\sigma (s) = \int _0^sT_\sigma (r)\, dr$ for $s \in \mathbb R$. Observe that
\begin{equation}\label{sep1}
0\leq \Lambda _\sigma (s) \leq \sigma |s|  \quad \hbox{for $s\in \mathbb R$.}
\end{equation}
 Moreover, one has that
\begin{equation}\label{aug30}
\int _0^\tau \langle v_t , T_\sigma (v) \rangle \, dt  = \int_\Omega \Lambda _\sigma  (v(x,\tau)\, dx -   \int_\Omega \Lambda _\sigma  (v(x,0)\, dx\quad \hbox{if $\tau  \in(0,T)$,}
\end{equation}
provided that $v\in L^2(\Omega_T)$ and $v_t  \in L^2(\Omega_T)$ -- see e.g. by \cite[Section 2.1.1]{GaMa}.
From equations \eqref{aug5}, \eqref{sep1}    and  \eqref{aug30}, with $v=u_k -u_m$,  one can deduce that
\begin{align}\label{aug6}
\int _0^T &\int_{\{|u_k-u_m|<\sigma\}}(|\nabla u_k|^{p-2}\nabla u_k - |\nabla u_m|^{p-2}\nabla u_m)\cdot (\nabla u_k - \nabla u_m)\, dxdt \\ \nonumber & \leq
 \sigma \int _0^T \int_\Omega  |f_k - f_m| \, dxdt + \sigma \int _\Omega |\psi _k - \psi_m|\, dx \leq \sigma C (\|f\|_{L^2(\Omega_T)} + \|\nabla \psi\|_{L^p(\Omega)}) \,
\end{align}
for  some constant $C=C(n,p, |\Omega|)$ and every $k, m \in \mathbb N$. Define
$$\kappa = \min \{(|\xi|^{p-2}\xi - |\eta|^{p-2}\eta) \cdot (\xi -\eta): |\xi|\leq \theta, |\eta|\leq \theta, |\xi-\eta|\geq \lambda\}\,.$$
Since $\kappa >0$, from inequality \eqref{aug6} one infers that
\begin{multline}\label{aug7}
\kappa \, |\{|\nabla u_k -\nabla u_m|>\lambda, |u_k -u_m|\leq\sigma, |\nabla u_k|\leq \theta, |\nabla u_m|\leq\theta \}| \\ \leq \sigma C (\|f\|_{L^2(\Omega_T)} + \|\nabla \psi\|_{L^p(\Omega)})\,,
\end{multline}
whence
\begin{align}\label{aug8}
|\{|\nabla u_k -\nabla u_m|>\lambda, |u_k -u_m|\leq\sigma, |\nabla u_k|\leq \theta, |\nabla u_m|\leq\theta \}| < \varepsilon
\end{align}
for every $k, m \in \mathbb N$, provided that $\sigma $ is sufficiently small. On the other hand, by property \eqref{aug2}, $\{u_k\}$ is a Cauchy sequence in measure, and hence
\begin{align}\label{aug9}
|\{|u_k-u_m|>\sigma\}|< \varepsilon
\end{align}
if $k$ and $m$ are sufficiently large. Inequalities \eqref{aug3}, \eqref{aug4}, \eqref{aug8} and \eqref{aug9} tell us that $\{\nabla u_k\}$ is actually a Cauchy sequence in measure. As a consequence, there exists a subsequence, still indexed by $k$, such that
\begin{equation}\label{aug10}
\nabla u_k \to \nabla u \quad \hbox{a.e. in $\Omega_T$.}
\end{equation}
Hence, equation \eqref{june0} holds. Moreover, passing to the limit  in inequality \eqref{solamaink} as $k \to \infty$, yields inequality \eqref{solamain}.
\\ In order to establish property \eqref{june0bis}, fix any $\tau \in [0,T]$ and make use of the test function $\phi= T_1(u_k-u_m)\chi _{[0, \tau)}$ in the weak formulation of problem \eqref{eqdirichletk}  and in its analogue with $k$ replaced by $m$. Here $\chi_E$ stands for the characteristic function of the set $E$.  Subtracting the resultant equations yields
\begin{align}\label{aug11}
\int _0^\tau & \langle (u_k - u_m)_t , T_1(u_k - u_m)\rangle \, dt\\  \nonumber & \quad  +  \int _0^\tau \int_\Omega (|\nabla u_k|^{p-2}\nabla u_k - |\nabla u_m|^{p-2}\nabla u_m)\cdot \nabla (T_1 (u_k - u_m))\, dx\,dt \\ \nonumber &
= \int _0^\tau \int_\Omega  (f_k - f_m)T_1 (u_k - u_m)\, dx\,dt \,.
\end{align}
Since the integrand in the second integral on the left-hand side of inequality \eqref{aug11} is nonnegative, one  infers via equation \eqref{aug30}, with $v= u_k -u_m$,  and \eqref{sep1} that
\begin{align}\label{aug12}
\int_\Omega \Lambda_1(u_k(x,\tau) - u_m(x,\tau))\, dx & \leq  \int_\Omega \Lambda_1(\psi_k  - \psi_m )\, dx  +
\int _0^\tau \int_\Omega  (f_k - f_m)T_1 (u_k - u_m)\, dxdt
\\ \nonumber & \leq \|\psi_k - \psi _m\|_{L^1(\Omega)} +
T  \|f_k - f_m\|_{L^1(\Omega_T)}  \\ \nonumber &\leq C( \|\nabla \psi_k - \nabla \psi _m\|_{L^p(\Omega)} +
 \|f_k - f_m\|_{L^2(\Omega_T)}) \,,
\end{align}
for some constant $C=C(n,p,|\Omega|, T)$ and every $\tau \in [0,T]$. On the other hand,
\begin{align}\label{aug13}
& \int_\Omega  \Lambda_1 (u_k(x,\tau) - u_m(x,\tau))\, dx  \\  \nonumber & \geq \int _{\{|u_k-u_m|\leq 1\}} |u_k(x,\tau)-u_m(x,\tau)|^2\, dx + \frac 12\int _{\{|u_k-u_m|> 1\}} |u_k(x,\tau)-u_m(x,\tau)|\, dx
\\ \nonumber & \geq
\frac 1{|\Omega|}\bigg(\int _{\{|u_k-u_m|\leq 1\}} |u_k(x,\tau)-u_m(x,\tau)|\, dx\bigg)^2 + \frac 12\int _{\{|u_k-u_m|> 1\}} |u_k(x,\tau)-u_m(x,\tau)|\, dx
\end{align}
for every $\tau \in [0,T]$.
Inequalities \eqref{aug12} and \eqref{aug13} ensure that $u_k$ is a Cauchy sequence in $C([0,T]; L^1(\Omega))$, and hence that  \eqref{june0bis} holds (up to subsequences).
\\ Finally, as far as uniqueness is concerned, assume that $u$ and $\overline u$ are   approximable solutions to problem \eqref{eqdirichlet}.
Let $\{f_k\}$, $\{\psi_k\}$ and $\{u_k\}$ be sequences as in the definition of approximable solution for $u$, and $\{\overline f_k\}$, $\{\overline \psi_k\}$ and $\{\overline u_k\}$ sequences as in a parallel definition for $\overline u$. On choosing, for $\sigma >0$, the test function $\phi = T_\sigma (u_k - \overline u_k)$ in the definitions of weak solutions for $u_k$ and $\overline u_k$,  and subtracting the equations so obtained, one deduces, analogously to \eqref{aug6},
\begin{align}\label{aug14}
\int _0^T &\int_{\{|u_k-\overline u_k|<\sigma\}}(|\nabla u_k|^{p-2}\nabla u_k - |\nabla \overline u_k|^{p-2}\nabla \overline u_k)\cdot (\nabla u_k - \nabla \overline u_k)\, dx\,dt \\ \nonumber & \leq
 \sigma \int _0^T \int_\Omega  |f_k - \overline f_k| \, dxdt + \sigma \int _\Omega |\psi _k - \overline \psi_k|\, dx
\end{align}
for $k \in \mathbb N$. By our assumptions, the right-hand side of inequality \eqref{aug14} approaches $0$ as $k \to \infty$, and $u_k \to u$, $\nabla u_k \to \nabla u$, $\overline u_k \to \overline u$ and $\nabla \overline u_k \to \overline u$ a.e. in $\Omega_T$. Hence, inequality  \eqref{aug14} implies, via Fatou's lemma, that
\begin{align}\label{aug15}
\int _0^T &\int_{\{|u-\overline u|<\sigma\}}(|\nabla u|^{p-2}\nabla u - |\nabla \overline u|^{p-2}\nabla \overline u)\cdot (\nabla u - \nabla \overline u)\, dx\,dt =0
\end{align}
for every $\sigma >0$. The integrand in \eqref{aug15} is nonnegative, and vanishes if and only if $\nabla u = \nabla \overline u$. Hence, by the arbitrariness of $\sigma$, we have that $\nabla u = \nabla \overline u$ a.e. in $\Omega_T$. Inasmuch as $u$ and $\overline u \in L^1((0,T); W^{1,1}_0(\Omega))$, we conclude that $u=\overline u$. \qed

We are now ready to prove Theorem \ref{secondcap}.

\medskip
\par\noindent
{\bf Proof of Theorem \ref{secondcap}}.
Assume, for the time being,
that
\begin{equation}\label{smooth}
f \in C^\infty _0 (\Omega_T) \quad \hbox{and} \quad \psi \in C^\infty _0(\Omega)\,.
\end{equation}
Then, there exists a unique weak solution $u \in C([0,T]; L^2(\Omega))\cap L^p((0,T); W^{1,p}_0(\Omega))$, with $u_t \in L^{p'}((0,T); W^{-1,p'}(\Omega))$,  to problem \eqref{eqdirichlet} in the sense of \eqref{lions}.
By inequality \eqref{solamaink}, with $f_k$, $\psi_k$ and $u_k$ replaced with  $f$, $\psi$ and $u$, we also have that $u_t \in L^2(\Omega _T)$
  Choose any test function $\phi$ of the form $\phi (x, t) = \varphi (x)\rho(t)$ in \eqref{lions}, with $\varphi \in W^{1,p}_0(\Omega)\cap L^\infty(\Omega)$ and $\rho \in C^\infty_0(0,T)$. Equation
\eqref{lions}
then entails that
\begin{equation}\label{june2}
 \int _0^T \rho(t)\int _\Omega |\nabla u(x,t)|^{p-2}\nabla u(x,t) \cdot \nabla \varphi (x)\, dx \,dt = \int _0^T\rho (t)\int_\Omega \varphi (x) (f(x,t) - u_t(x,t))\, dx \, dt\,.
\end{equation}
Hence, by  the arbtrariness of the function $\rho$,
\begin{equation}\label{june3}
\int _\Omega |\nabla u(x,t)|^{p-2}\nabla u(x,t) \cdot \nabla \varphi (x)\, dx = \int_\Omega \varphi (x) (f(x,t) - u_t(x,t))\, dx \quad \hbox{for a.e. $t \in (0,T)$,}
\end{equation}
and for every $\varphi \in W^{1,p}_0(\Omega)\cap L^\infty(\Omega)$. Now, define,
 for each $t \in (0,T)$, the function $g^t: \Omega \to \mathbb R$ as
\begin{equation}\label{g}
g^t(x) =f(x,t) - u_t(x,t) \quad \hbox{for $x \in \Omega$,}
\end{equation}
and the function $w^t : \Omega \to \mathbb R$ as
\begin{equation}\label{w}
w^t(x) = u(x,t) \quad \hbox{for $x \in \Omega$.}
\end{equation}
Then $g^t \in L^2(\Omega)$ and $w^t \in W^{1,p}_0(\Omega)$ for a.e. $t \in (0, T)$. Moreover, equation \eqref{june3} reads
\begin{equation}\label{june3bis}
\int _\Omega |\nabla w^t|^{p-2}\nabla w^t \cdot \nabla \varphi \, dx = \int_\Omega \varphi \,  g^t \, dx
\end{equation}
for every $\varphi \in W^{1,p}_0(\Omega)\cap L^\infty(\Omega)$ and for a.e. $t \in (0, T)$. Fix any such $t$, and consider the (elliptic) Dirichlet problem
\begin{equation}\label{aux}
\begin{cases}
- {\rm div} (|\nabla v^t|^{p-2}\nabla v^t ) = g^t  & {\rm in}\,\,\, \Omega \\
 v^t =0  &
{\rm on}\,\,\,
\partial \Omega \,.
\end{cases}
\end{equation}
This problem has a generalized approximable solution $v^t$ in the following sense. For every $\sigma >0$, the function $T_\sigma (v^t) \in W^{1,p}_0(\Omega)$,  and there exists a measurable function $Z^t : \Omega \to {{\R}^n}$ such that
\begin{equation}\label{sep2} \nabla T_\sigma (v^t)  = Z^t \chi_{\{|v^t|<\sigma\}} \quad \hbox{a.e. in $\Omega$,}
 \end{equation}
 for every $\sigma >0$. Moreover, there exist sequences
$\{g_k^t\} \subset C^\infty_0(\Omega)$ and $\{v_k^t\} \subset W^{1,p}_0(\Omega)$ such that $g_k^t \to g$ in $L^1(\Omega)$,
 $v_k^t$ is the weak solution to problem
\begin{equation}\label{L1k}
\begin{cases}
- {\rm div} (|\nabla v_k^t|^{p-2} \nabla v_k^t ) = g_k^t  & {\rm in}\,\,\, \Omega \\
 v_k^t =0  &
{\rm on}\,\,\,
\partial \Omega\,,
\end{cases}
\end{equation}
and
\begin{equation}\label{L1conv}
v_k^t \to v^t \quad \hbox{and} \quad \nabla v_k \to Z^t \quad \hbox{a.e. in $\Omega$.}
\end{equation}
This follows e.g. from   \cite[Theorem 3.2]{cmapprox}. We claim that
\begin{equation}\label{june4}
w^t= v^t\,.
\end{equation}
In order to verify this claim, observe that $T_\sigma(w^t-v_k^t) \in W^{1,p}_0(\Omega) \cap L^\infty (\Omega)$ for every $k\in \mathbb N$ and $\sigma>0$,
since   $w^t \in W^{1,p}_0(\Omega)$ and $v_k^t \in W^{1,p}_0(\Omega)$. Making use of the test function $\varphi = T_\sigma(w^t-v_k^t)$   in equation \eqref{june3bis} and in the definition of weak solution to problem \eqref{L1k}, and subtracting the resultant equations tells us that
\begin{equation}\label{sola3}
\int_{\{|w^t-v_k^t|< \sigma\}}(|\nabla w^t|^{p-2}\nabla w^t - |\nabla v_k^t|^{p-2}\nabla v_k^t)\cdot (\nabla w^t - \nabla v_k^t)\, dx = \int _\Omega (g^t -g_k^t)\, T_\sigma(w^t-v_k^t) \, dx\,.
\end{equation}
Since $|T_\sigma(w^t-v_k^t)|\leq \sigma$ and $g_k^t \to g^t$ in $L^1(\Omega)$, the right-hand side of equation \eqref{sola3} tends to $0$ as $k \to \infty$. Thus, inasmuch as  the integrand on the left-hand side  is nonnegative,  passing to the limit in \eqref{sola3} as $k \to \infty$ yields, by Fatou's lemma and equation \eqref{L1conv},
\begin{equation}\label{sola4}
\int_{\{|w^t-v^t|< \sigma\}}(|\nabla w^t|^{p-2}\nabla w^t - |Z^t|^{p-2}Z^t)\cdot (\nabla w^t - Z^t)\, dx =0
\end{equation}
for every $\sigma>0$. The integrand in \eqref{sola4} in nonnegative, and vanishes if and only if $\nabla w^t = Z^t$. Therefore,  for every $\sigma>0$,  we have that $\nabla w^t = Z^t$ a.e. in the set $\{|w^t-v^t|< \sigma\}$, whence
\begin{equation}\label{sola5}
\nabla w^t = Z^t \qquad \hbox{a.e. in $\Omega$.}
\end{equation}
Now, the function $T_\sigma (w^t-T_\varrho(v^t)) \in W^{1,p}_0(\Omega)\subset W^{1,1}_0(\Omega)$ for every $\sigma, \varrho >0$. An application of the Sobolev inequality in $W^{1,1}_0(\Omega)$ and the use of equations \eqref{sep2} and  \eqref{sola5} imply that
\begin{equation}\label{sola6}
\bigg( \int_{\Omega }|T_{\sigma}(w^t-T_\varrho(v^t)|^{n'}\, dx \bigg
)^{n'}
\le C
\bigg(
\int_{\{\varrho<|w^t|<\varrho+\sigma\} }|\nabla w^t|\, dx+
\int_{\{\varrho-\sigma<|w^t|<\varrho\} }|\nabla w^t|\, dx
\bigg)
\end{equation}
for some constant $C=C(n)$ and for every $\sigma, \varrho >0$. Since $|\nabla w^t| \in L^1(\Omega)$, for each fixed $\sigma >0$ the right-hand side of inequality \eqref{sola6} converges to $0$ as $\varrho\to \infty$. Therefore,  passing to the limit  in \eqref{sola6} as $\varrho\to \infty$ tells us that
$$\int_{\Omega }|T_{\sigma}(w^t-v^t)|^{n'}\, dx =0$$
for every $\sigma >0$. Hence, equation \eqref{june4} follows, on letting $\sigma \to \infty$.
\iffalse

One has  that $u_t (\cdot, \tau)  \in L^2(\Omega )$ and $u(\cdot , t) \in  W^{1,p}_0(\Omega)$ for a.e. $t \in (0,T)$. Fix any such $t$, and define $g^t: \Omega \to \mathbb R$ as $g^t(x) =f(x,t) - u_t(x,t)$. Thus, assumptions of Proposition \ref{sola} are satisfied with $w(\cdot) = u(\cdot , t)$ and $g = g^t$.
Consider the approximable soluton $v^t$ to the Dirichlet problem
\begin{equation}\label{aux}
\begin{cases}
- {\rm div} (|\nabla v^t|^{p-2}\nabla v^t ) = g^t  & {\rm in}\,\,\, \Omega \\
 v^t =0  &
{\rm on}\,\,\,
\partial \Omega \,.
\end{cases}
\end{equation}
By Proposition \ref{sola},
\begin{equation}\label{june4}
u(\cdot , t)= v^t\,.
\end{equation}

\fi
An application of \cite[Theorem 2.1]{cmsecondsys} tells us that
\begin{equation}\label{aug25}
|\nabla v^t|^{p-2} \nabla v^t \in W^{1,2}(\Omega),\end{equation}
 and there exists a constant $C=C(n,p, \Omega)$ such that
\begin{equation}\label{aug26}
\||\nabla v^t|^{p-2} \nabla v^t\|_{W^{1,2}(\Omega)}^2 \leq C \|g^t\|_{L^2(\Omega)}.
\end{equation}
Hence, owing to equations \eqref{june4}, \eqref{g} and \eqref{w}
$$ |\nabla u(\cdot , t)|^{p-2} \nabla u(\cdot , t) \in W^{1,2}(\Omega)$$
for a.e. $t \in (0,T)$,
and
\begin{equation}\label{june5}
\||\nabla u(\cdot , t)|^{p-2} \nabla u(\cdot , t)\|_{W^{1,2}(\Omega)}^2 \leq C \int_\Omega f(x,t)^2 + u_t(x,t)^2\, dx
\end{equation}
for the same constant  $C=C(n,p, \Omega)$.
Integrating inequality \eqref{june5} with respect to $t$ over $(0,T)$ yields
\begin{equation}\label{june6}
\||\nabla u|^{p-2} \nabla u\|_{L^2((0,T); W^{1,2}(\Omega))}^2 \leq C \big(\|f\|_{L^2(\Omega_T)}^2 + \|u_t\|_{L^2(\Omega_T)}^2\big)\,.
\end{equation}
Coupling inequalities  \eqref{june6} and \eqref{solamain} tells us that $|\nabla u|^{p-2} \nabla u \in L^2((0,T); W^{1,2}(\Omega))$, and
\begin{align}\label{june7}
\||\nabla u|^{p-2} \nabla u\|_{L^2((0,T); W^{1,2}(\Omega))} \leq C \big(\|f\|_{L^2(\Omega_T)} + \|\nabla \psi\|_{L^p(\Omega)}^{\frac p2}\big)\,
\end{align}
constant  $C=C(n,p, \Omega)$.
\\
We have now to remove the additional assumptions  \eqref{smooth}. To this purpose, consider sequences  $\{f_k\}$, $\{\psi_k\}$ and $\{u_k\}$  as in the definition of approximable solution to problem \eqref{eqdirichlet}. Clearly, we may assume that $\|f_k\|_{L^2(\Omega_T)} \leq 2 \|f\|_{L^2(\Omega_T)}$ and $\|\nabla \psi _k\|_{L^p(\Omega)} \leq 2\|\nabla \psi \|_{L^p(\Omega)}$ for $k \in \mathbb N$.
\iffalse

such that
\begin{equation}\label{june1}
f_k \to f \qquad \hbox{in $L^2(\Omega_T)$,}
\end{equation}
\begin{equation}\label{june2}
\psi_k \to \psi \qquad \hbox{in $W^{1,p}(\Omega)\cap L^2(\Omega)$,}
\end{equation}
and the sequence $\{u_k\}$ of
 weak solutions to the problems
\begin{equation}\label{eqdirichletk}
\begin{cases}
(u_k)_t - {\rm div} (|\nabla u_k|^{p-2} \nabla u_k ) = f_k  & {\rm in}\,\,\, \Omega \times (0, T)\\
 u_k(x,t) =0  &
{\rm for}\,\,\,
(x,t)\in \partial \Omega \times (0, T) \,
\\ u_k(x,0)= \psi_k (x) &
{\rm for}\,\,\, x\in \Omega\,
\end{cases}
\end{equation}
is such that \todo{Check if this is a suitable definition of approximable solution}
\begin{equation}\label{june0}
u_k \to u \quad \hbox{and} \quad \nabla u_k \to Z \quad \hbox{a.e. in $\Omega_t$.}
\end{equation}

\fi
\iffalse

Definition \eqref{weakpar} of weak solution $u_k$ to problem \eqref{eqdirichletk} tells us that
\begin{multline}\label{weakapp}
\int _\Omega u_k(x,\tau)\phi (x,\tau)\, dx  + \int _0^\tau\int _\Omega -u_k \phi _t\, dx \, dt + \int _0^\tau\int _\Omega |\nabla u_k|^{p-2} \nabla u_k \cdot \nabla \phi \,dx\, dt  \\
=  \int_\Omega \psi _k(x) \phi(x, 0)\, dx  + \int _0^\tau\int _\Omega f_k \phi\, dx\, dt
\end{multline}
for  every $\tau \in (0, T]$ and every function $\phi \in W^{1,2}(0,T; L^2(\Omega))\cap  L^p(0,T; W^{1,p}_0(\Omega))$.
\\

\fi
By inequality    \eqref{june7}, applied with $u$ replaced by $u_k$, and our assumptions on the sequences $\{f_k\}$ and $\{\psi _k\}$, there exists a constant $C=C(n,p, \Omega)$  such that
\begin{align}\label{june8}

 \||\nabla u_k|^{p-2} \nabla u_k\|_{L^2((0,T); W^{1,2}(\Omega))}
  \leq C  \big(\|f\|_{L^2(\Omega_T)} + \|\nabla \psi\|_{L^p(\Omega)}^{\frac p2}\big)\,
\end{align}
for $k \in \mathbb N$.
\iffalse

By the boundedness of the first and third norm on the left-hand side of inequality \eqref{june8},
there exists a subsequence, still denoted by $\{u_k\}$, and a function $w\in W^{1,(p,2)}(\Omega_T)\cap L^2(0,T; W^{1,2}(\Omega))$ such that
\begin{equation}\label{june20}
u_k \to w \quad \hbox{in $L^r(\Omega_T)$,}
\end{equation}
where $r=???$,
\begin{equation}\label{june21}
u_k \rightharpoonup w \quad \hbox{in $W^{1,(p,2)}(\Omega_T)$.}
\end{equation}
Moreover, by the boundedness of the middle norm on the left-hand side of inequality \eqref{june8},

\fi
Therefore, there exist
 an ${{\R}^n}$--valued  function $U \in L^2(\Omega_T)$ and an $\mathbb R^{n\times n}$--valued function $V \in L^2(\Omega_T)$ such that
\begin{equation}\label{par25k}
|\nabla u_{k}|^{p-2} \nabla u_k\rightharpoonup U \quad \hbox{weakly in $L^2(\Omega_T)$}
\end{equation}
and
\begin{equation}\label{par26k}
\nabla (|\nabla u_{k}|^{p-2} \nabla u_{k})\rightharpoonup V \quad \hbox{weakly in $L^2(\Omega_T)$,}
\end{equation}
up to subsequences.
Thereby, $U \in L^2((0, T); W^{1,2}(\Omega))$ and
\begin{equation}\label{par38}
V= \nabla U\,.
\end{equation}
Owing to \eqref{june0} and \eqref{par25k},
$U = |\nabla u|^{p-2} \nabla u$, whence
\begin{equation}\label{june6k}
|\nabla u_k|^{p-2} \nabla u_k  \rightharpoonup |\nabla u|^{p-2} \nabla u \quad \hbox{weakly in $L^2((0,T); W^{1,2}(\Omega))$.}
\end{equation}
 Hence, inequality \eqref{secondcap2} follows,   via \eqref{june8}  and \eqref{june6k}.
\qed

\medskip
\par\noindent
{\bf Proof of Theorem \ref{secondconvex}}. The proof is completely analogous to that of Theorem \ref{secondcap}, save that   properties \eqref{aug25} and \eqref{aug26} now follow from \cite[Theorem 2.3]{cmsecond}, that holds for any bounded open convex set.
 \qed

\medskip
\par
We conclude with the proofs of Theorems \ref{secondcapsys} and \ref{secondconvexsys} for systems. They  require some minor variant with respect to those of the corresponding results for equations.

\medskip
\par\noindent
{\bf Proof of Theorem \ref{secondcapsys}}.
Since we are assuming that $p\geq 2$,
$$f \in L^2(\Omega_T) \to L^{p'}(\Omega_T) \to  L^{p'}((0,T); W^{-1,p'}(\Omega)) \quad \hbox{and} \quad \psi \in W^{1,p}_0(\Omega) \to L^2(\Omega)\,.$$
 Therefore, there exists a unique weak solution $u \in C([0,T]; L^2(\Omega))\cap L^p((0,T); W^{1,p}_0(\Omega))$, with $u_t \in L^{p'}((0,T); W^{-1,p'}(\Omega))$,  to problem \eqref{eqdirichlet}. Moreover, owing to estimate  \eqref{solamaink} (that also holds when $N>1$) with $f_k$, $\psi_k$ and $u_k$ replaced by  $f$, $\psi$ and $u$, we have that
 $u_t \in L^2(\Omega _T)$ and there exists a constant $C$ such that
\begin{align}\label{solaN}
\| u\|_{L^\infty((0,T); W^{1,p}_0(\Omega))}^{\frac p2}+ \|u_t\|_{L^2(\Omega_T)}& \leq C \big(\|f\|_{L^2(\Omega_T)} + \|\nabla \psi\|_{L^p(\Omega)}^{\frac p2}\big)  \,.
\end{align}
For  $t\in [0,T]$, define the functions $g^t: \Omega \to \mathbb R^N$ and $w^t : \Omega \to \mathbb R^N$ as
$$g^t(x) =f(x,t) - u_t(x,t) \quad  \hbox{and} \quad   w^t(x) = u(x,t) \quad \hbox{for $x \in \Omega$.}$$
The same argument as in the proof of  equation \eqref{june3bis} tells us that the function $w^t$ is the weak solution to the Dirichlet problem
\begin{equation}\label{auxN}
\begin{cases}
- {\rm div} (|\nabla w^t|^{p-2}\nabla w^t ) = g^t  & {\rm in}\,\,\, \Omega \\
 w^t =0  &
{\rm on}\,\,\,
\partial \Omega \,
\end{cases}
\end{equation}
for a.e. $t\in (0, T)$. Observe that $w^t$ is actually the weak solution to problem \eqref{auxN} in the standard sense, since equation \eqref{june2}, and hence \eqref{june3bis} now hold for every $\varphi \in W^{1,p}_0(\Omega)$. Indeed, inasmuch as $p \geq 2$, one has that $L^p((0,T); W^{1,p}_0(\Omega)) \to L^2(\Omega)$, and since $\rho \varphi \in L^p((0,T); W^{1,p}_0(\Omega))$ if $\rho \in C^\infty _0(0, T)$ and $\varphi \in W^{1,p}_0(\Omega)$, one has that
$$  \langle u_t , \rho\, \varphi \rangle = \int _0^T \rho (t) \int _\Omega u_t (x,t) \varphi (x,t)\, dxdt\,.$$
Next,  fix any  $t\in (0, T)$ for which \eqref{auxN} holds, and let
$\{g^t_k\} \subset C^\infty_0(\Omega)$ be a sequence such that $g^t_k \to g^t$ in $L^2(\Omega)$ and $\|g^t_k\|_{L^2(\Omega)} \leq 2 \|g^t\|_{L^2(\Omega)}$ for $k \in \mathbb N$. For each $k$, let $w^t_k$ be the weak solution to the problem
\begin{equation}\label{auxNk}
\begin{cases}
- {\rm div} (|\nabla w^t_k|^{p-2}\nabla w^t_k ) = g^t_k  & {\rm in}\,\,\, \Omega \\
 w^t_k =0  &
{\rm on}\,\,\,
\partial \Omega \,.
\end{cases}
\end{equation}
The use of $w^t_k$ as a test function in the weak formulation of problem \eqref{auxNk} and of the H\"older and  the Sobolev inequalities tells  us that
\begin{equation}\label{auxN2}
\|w^t_k\|_{L^2(\Omega)} \leq c \|g^t\|_{L^2(\Omega)}
\end{equation}
for $k \in\mathbb N$, for some constant $c=c(n,N,p, |\Omega|)$. On the other hand,
choosing  $w^t_k - w^t$ as a test function in the weak formulations of problems \eqref{auxN} and \eqref{auxNk}, subtracting the resultant equations, and observing that
$$c |\xi - \eta|^p \leq (|\xi|^{p-2} \xi - |\eta |^{p-2}\eta) \cdot (\xi - \eta) \quad \hbox{for $\xi, \eta \in \mathbb R^{N\times n}$,}$$
for some positive constant $c=c(n,N,p)$ yield
\begin{align}\label{auxN1}
c \int _\Omega |\nabla w^t_k - \nabla w^t|^p\, dx&  \leq
\int_\Omega(|\nabla w^t_k|^{p-2}\nabla w^t_k - |\nabla w^t|^{p-2}\nabla w^t)\cdot (\nabla w^t_k - \nabla w^t)\, dx \\ \nonumber & = \int _\Omega (g^t_k -g^t) \cdot (w^t_k-w^t) \, dx  \leq
\|g^t_k -g^t\|_{L^2(\Omega)} \|w^t_k-w^t\|_{L^2(\Omega)} \\ \nonumber &\leq  \|g^t_k -g^t\|_{L^2(\Omega)} (\|w^t_k\|_{L^2(\Omega)}+ \|w^t\|_{L^2(\Omega)})
\end{align}
for $k \in\mathbb N$. Hence, by inequality \eqref{auxN2}, $w^t_k \to w^t$ in $W^{1,p}_0(\Omega)$ as $k \to \infty$. As a conequence, for a.e. $t\in (0,T)$ the function $w^t$ is an approximable solution to problem \eqref{auxN}. An application of
  \cite[Theorem 2.1]{cmsecondsys} then entails that
$ |\nabla u(\cdot , t)|^{p-2} \nabla u(\cdot , t)= |\nabla w^t(\cdot)|^{p-2}\nabla w^t(\cdot) \in W^{1,2}(\Omega)$ and
\begin{equation}\label{june5N}
\||\nabla u(\cdot , t)|^{p-2} \nabla u(\cdot , t)\|_{W^{1,2}(\Omega)}^2 \leq C \int_\Omega f(x,t)^2 + u_t(x,t)^2\, dx \quad \hbox{for a.e. $t \in (0,T)$,}
\end{equation}
some constant  $C=C(n,N,p, \Omega)$. Equations \eqref{secondcapsys1} and \eqref{secondcapsys2} follow on
integrating inequality \eqref{june5N} with respect to $t$ over $(0,T)$, and exploiting inequality \eqref{solaN}. \qed

\medskip
\par\noindent
{\bf Proof of Theorem \ref{secondconvexsys}}. The proof is the same as that of Theorem \ref{secondcapsys},  inequality \eqref{june5N}  being  a consequence of \cite[Theorem 2.6]{cmsecondsys} under the assumption that $\Omega$ is a bounded convex set.
 \qed

\color{black}

\bigskip
\par\noindent
{\bf Acknowledgment}. This research was initiated during a visit of the first-named author at the Institut Mittag-Leffler in August 2017. He wishes to thank the Director and the Staff of the Institut for their support and hospitality.

\color{black}
\section*{Compliance with Ethical Standards

\smallskip
\par\noindent
{\bf Funding}. This research was partly funded by:
\\ (i) Research Project of the
Italian Ministry of University and Research (MIUR) Prin 2012 \lq\lq Elliptic and
parabolic partial differential equations: geometric aspects, related
inequalities, and applications"  (grant number 2012TC7588);
\\ (ii) GNAMPA   of the Italian INdAM - National Institute of High Mathematics (grant number not available);

\smallskip
\par\noindent
{\bf Conflict of Interest}. The authors declare that they have no conflict of interest.