Let

- $d\in\mathbb N$
- $\lambda^d$ denote the Lebesgue measure on $\mathbb R^d$
- $\Lambda\subseteq\mathbb R^d$ be open
- $\mathcal V:=\left\{v\in C_c^\infty(\Lambda)^d:\nabla\cdot v=0\right\}$, $$V:=\overline{\mathcal V}^{\left\|\;\cdot\;\right\|_{H^1(\Lambda,\:\mathbb R^d)}}$$ and $$H:=\overline{\mathcal V}^{\left\|\;\cdot\;\right\|_{L^2(\Lambda,\:\mathbb R^d)}}=\overline V^{\left\|\;\cdot\;\right\|_{L^2(\Lambda,\:\mathbb R^d)}}\tag 1$$
- $W^1((0,T),V):=\left\{u\in L^1_{\text{loc}}((0,T),V):u\text{ is weakly differentiable}\right\}$

I want to reformulate $$\left\{\begin{array}{rll}\displaystyle\frac{\partial u}{\partial t}-\nu\Delta u+\left(u\cdot\nabla\right)u+\frac1\rho\nabla p&=&f&&\text{in }[0,T]\times\Lambda\\ u(0,\;\cdot\;)&=&u_0&&\text{in }\Lambda\end{array}\right.\tag 2$$ as a semilinear evolution equation on $H$.

It's easy to see that $$(A_0u)v:=\sum_{i=1}^d\langle\nabla u_i,\nabla v_i\rangle_{L^2(\Lambda,\:\mathbb R^d)}\;\;\;\text{for }u,v\in V$$ is a bounded linear operator from $V$ to $V'$. Moreover, if $d\le 4$ and $\Lambda$ is bounded, then $$b(u,v,w):=\int_\Lambda(u\cdot\nabla)v\cdot w\:{\rm d}\lambda^d\;\;\;\text{for }u,v,w\in H_0^1(\Lambda,\mathbb R^d)$$ is a well-defined bounded trilinear form on $H_0^1(\Lambda,\mathbb R^d)$ and hence $$B(u,v):=\left.b(u,v,\;\cdot\;)\right|_V\;\;\;\text{for }u,v\in V$$ is a bounded bilinear operator from $V\times V$ to $V'$. Now, let $u$ be a classical solution of $(2)$ and $$\tilde u(t):=u(t,\;\cdot\;)\;\;\;\text{for }t\in(0,T)\;.$$ A straightforward calculation yields $\tilde u\in W^1((0,T),V)$ with $$\tilde u'(t)+\nu A_0\tilde u(t)+B(\tilde u(t),\tilde u(t))=0\;\;\;\text{for all }t\in[0,T].\tag 3$$

Note that $(3)$ is an equation in $V'$ and hence any solution $\tilde u\in W^1((0,T),V)$ of $(3)$ corresponds to the notion of a

weaksolution of $(2)$. Since I'm interested in amildsolution of $(2)$, I somehow need to reformulate $(2)$ as an equation in $H$.

It's well-known that $$A:=\left.A_0\right|_{\mathcal D(A)}$$ with $D(A):=V\cap H^2(\Lambda,\mathbb R^d)$ is a densely-defined linear operator on $H$. Moreover, $-A$ is the generator of a contraction $C_0$-semigroup $S$ on $H$.

Now, I've often seen that people call $$v(t)=S(t)v_0+\int_0^tS(t-s)B(v(s),v(s))\:{\rm ds}\;\;\;\text{for all }t\in[0,T]\tag 4$$ the

mild equationcorresponding to $(2)$.

However, I expect that a solution $v$ of $(4)$ is $H$-valued! As stated above, if $v(t)\in\mathcal D(A)$, then $A_0v(t)\in H$ for all $t\in[0,T]$, but the $B(v(t),v(t))$ are still $V'$-valued.

Is this simply ruled out by $V'\cong V$ or has $(4)$ to be understood in a special sense?

**In order to prevent any confusion**: My problem with $(4)$ is that $A_0v(t)+B(v(t),v(t))\not\in H$, cause $B(v(t),v(t))\not\in H$, even when $v(t)\in\mathcal D(A)$, for any $t\in[0,T]$.