[step:Pass the mollified identity to the limit in the duality pairing]
As $\varepsilon \to 0$, Bochner convolution gives
\begin{align*}
u_\varepsilon \to u \quad \text{in } L^2((a,b);V)
\end{align*}
and
\begin{align*}
f_\varepsilon \to u' \quad \text{in } L^2((a,b);V^*).
\end{align*}
The continuity of the embedding $V \hookrightarrow H$ also gives
\begin{align*}
u_\varepsilon \to u \quad \text{in } L^2((a,b);H).
\end{align*}
We compare the duality terms by writing
\begin{align*}
f_\varepsilon(u_\varepsilon)-u'(u)=f_\varepsilon(u_\varepsilon-u)+(f_\varepsilon-u')(u).
\end{align*}
Using the defining operator norm on $V^*$ and Cauchy-Schwarz on the [measure space](/page/Measure%20Space) $((a,b),\mathcal{B}((a,b)),\mathcal{L}^1)$,
\begin{align*}
\int_a^b |f_\varepsilon(\tau)(u_\varepsilon(\tau)-u(\tau))|\,d\mathcal{L}^1(\tau) \le \|f_\varepsilon\|_{L^2((a,b);V^*)}\|u_\varepsilon-u\|_{L^2((a,b);V)}.
\end{align*}
The right-hand side tends to $0$ because $(f_\varepsilon)$ is bounded in $L^2((a,b);V^*)$ and $u_\varepsilon \to u$ in $L^2((a,b);V)$. Similarly,
\begin{align*}
\int_a^b |(f_\varepsilon(\tau)-u'(\tau))(u(\tau))|\,d\mathcal{L}^1(\tau) \le \|f_\varepsilon-u'\|_{L^2((a,b);V^*)}\|u\|_{L^2((a,b);V)},
\end{align*}
which tends to $0$. Therefore
\begin{align*}
f_\varepsilon(u_\varepsilon)\to u'(u)
\end{align*}
in $L^1((a,b))$.
Moreover,
\begin{align*}
\int_a^b \left|\|u_\varepsilon(\tau)\|_H^2-\|u(\tau)\|_H^2\right|\,d\mathcal{L}^1(\tau) \le \left(\|u_\varepsilon\|_{L^2((a,b);H)}+\|u\|_{L^2((a,b);H)}\right)\|u_\varepsilon-u\|_{L^2((a,b);H)},
\end{align*}
so $\|u_\varepsilon\|_H^2\to \|u\|_H^2$ in $L^1((a,b))$. Let $\psi \in C_c^\infty((a,b))$. For each $\varepsilon$, the differential identity
\begin{align*}
\frac{d}{dt}\|u_\varepsilon(t)\|_H^2=2f_\varepsilon(t)(u_\varepsilon(t))
\end{align*}
holds on $(a,b)$. Testing this identity against $\psi$, integrating the left-hand side by parts on $(a,b)$ with respect to $\mathcal{L}^1$, and using the compact support of $\psi$ gives
\begin{align*}
-\int_a^b \|u_\varepsilon(\tau)\|_H^2\psi'(\tau)\,d\mathcal{L}^1(\tau)=2\int_a^b f_\varepsilon(\tau)(u_\varepsilon(\tau))\psi(\tau)\,d\mathcal{L}^1(\tau).
\end{align*}
Since $\psi$ and $\psi'$ are bounded and compactly supported in $(a,b)$, convergence in $L^1((a,b))$ implies convergence after multiplication by $\psi$ and $\psi'$, respectively, by the elementary estimate
\begin{align*}
\|\psi F\|_{L^1((a,b))}\le \|\psi\|_\infty\|F\|_{L^1((a,b))}.
\end{align*}
Passing to the limit in the two tested integrals therefore gives
\begin{align*}
-\int_a^b \|u(\tau)\|_H^2\psi'(\tau)\,d\mathcal{L}^1(\tau)=2\int_a^b u'(\tau)(u(\tau))\psi(\tau)\,d\mathcal{L}^1(\tau).
\end{align*}
Thus, in the sense of distributions on $(a,b)$,
\begin{align*}
\frac{d}{dt}\|u(t)\|_H^2=2u'(t)(u(t)).
\end{align*} Since $a$ and $b$ were arbitrary interior endpoints, the identity holds distributionally on $(0,T)$. The function
\begin{align*}
q: (0,T) \to \mathbb{R}
\end{align*}
defined by
\begin{align*}
q(t)=2u'(t)(u(t))
\end{align*}
belongs to $L^1(0,T)$, because
\begin{align*}
\int_0^{\mathsf T} |q(t)|\,d\mathcal{L}^1(t) \le 2\|u'\|_{L^2((0,T);V^*)}\|u\|_{L^2((0,T);V)}.
\end{align*}
Hence $\|u(t)\|_H^2$ has an absolutely continuous representative $g:[0,T]\to \mathbb{R}$ satisfying
\begin{align*}
g(t)-g(s)=2\int_s^t u'(\tau)(u(\tau))\,d\mathcal{L}^1(\tau)
\end{align*}
for all $0\le s\le t\le T$.
[/step]