II. Weak Solutions and Uniqueness
PDF Version
Weak Solutions
Let me remind you what I did in a previous lecture. I considered the problem
[latex]
\begin{align}
u_{t}-\text{div}\big( u^{\sigma} \nabla u\big) &=0,~x \in\R^{N},~t>0; \ \sigma > 0,\label{eq:inst-pme-1}\\
\int_{\R^{N}}u(x,t)\,dx &=1,~t > 0\label{eq:inst-pme-2}\\
u(x,0) = \delta(x)\label{eq:inst-pme-3}
\end{align}
[/latex]
and constructed instantaneous point-source solution, also called the Zeldovich-Kompaneets-Barenblatt (ZKB) solution.
\begin{equation} u_*(x,t)=t^{-\frac{N}{2+N\sigma}}\left[\frac{\sigma}{2(2+N\sigma)} \left( \eta_{0}^{2}-\frac{|x|^{2}}{t^{\frac{2}{2+N\sigma}}}\right)_{+} \right]^{\frac{1}{\sigma}}\end{equation}
Note that the nonlinear diffusion equation is invariant under translation of time and space coordinate, and the solution to the problem
\[ \begin{cases}
u_{t}=\text{div}\big(u^{\sigma }u\big),~&(x,t) \in \R^{N}\times \R_{+}\\
u(x,0)=u_*(x,T),~&x \in \R^{N}
\end{cases}\]
is exactly $u_*(x,t+T)$, which has compact support for all $t\geq 0$.
However, several important questions are left open. First of all, ZKB solution is not a classical solution: it is not even differentiable on the boundary of the support. First important question which we need to answer is the following: In what sense is $u_*$ actually a solution of the problem? Hence, we need to define the notion of the weak solution. Second important question we need to answer is the following: May be there is a smooth solution of the same problem with different properties, and ZKB solution is just physically irrelevant mathematical example. Main goal of this lecture is to answer these important questions.
For clarity, I will concentrate on the one particular problem – Dirichlet problem for the nonlinear diffusion equation. I am going to formulate the notion of the weak solution and prove its uniqueness. Let me note that the pioneering work on the mathematical theory of the nonlinear diffusion type equations is [1]. For complete review and modern treatment I refer to [2].
Consider the Dirichlet problem for the nonlinear diffusion equation:
\begin{align}
u_{t}&=\lapl u^{m},\quad \text{in}~Q_{T}=\Omega \times (0,T]\label{eq:pme}\\
u(x,0)&=u_{0}(x),\quad \text{in}~\Omega\label{eq:pme-iv}\\
u(x,t)&=0,\quad\text{on}~S_{T}=\partial \Omega \times (0,T]\label{eq:pme-bc}
\end{align}
where $m>1$, $\Omega \subset \R^{N}$, with $N \geq 1$, is a bounded domain with boundary $\bdy \Omega \in C^{2+\alpha}$, $\alpha \in (0,1)$, $u_{0}(x)\geq 0$;
Assuming temporarily that $u$ is a smooth solution, multiply \eqref{eq:pme} by a smooth test function (whose properties are to be determined) and integrate over $Q_{T}$ to find
\begin{align*}
0&=\iint_{Q_{T}} \left( u_{t}-\lapl u^{m}\right)\phi\,dx\,dt
\end{align*}
Integrating by parts with respect to $x,t$, it follows that
\begin{align*}
0&=\iint_{Q_{T}} \left[ \nabla u^{m}\cdot \nabla \phi – u \phi_{t}\right]\,dx\,dt
+\int_{\Omega}u(x,T)\phi(x,T)\,dx-\\
&\qquad – \int_{\Omega}u_{0}(x)\phi(x,0)\,dx – \int_{S_{T}}\D{u^{m}}{\nu}\phi\,dx\,dt
\end{align*}
where $\nu$ is an exterior normal at the boundary point $x \in \partial\Omega$.
To simplify further, choose $\phi\big\vert_{t=T}=\phi\big\vert_{S_{T}}=0$ to derive
\begin{equation} \iint_{Q_{T}} \left( \nabla u^{m} \cdot \nabla \phi – u \phi_{t}\right)\,dx\,dt = \int_{\Omega}u_{0}(x)\phi(x,0)\,dx\label{eq:pme-weak-int}\end{equation}
which holds for all $\phi \in C^{1}(\bar{Q}_{T})$ with $\phi(x,T)=0$, $\phi\big\vert_{S_{T}}=0$.
The integral identity clearly holds for any smooth or classical solution of \eqref{eq:pme}-\eqref{eq:pme-bc}. However, opposite is not necessarily true, and in particular ZKB solution satisfies the integral identity \eqref{eq:pme-weak-int}, but not the equation \eqref{eq:pme}. Hence, in order to give sense to ZKB type solutions, it makes sense to weaken the smoothness requirements on the solution by replacing PDE \eqref{eq:pme} with the integral identity \eqref{eq:pme-weak-int}. To exploit this idea, let us formulate minimal requirements on $u$ which are enough to hold \eqref{eq:pme-weak-int}. Hence, integral identity will dictate the relevant Sobolev space for solutions. First of all, we are going to understand integration in the Lebesgue’s sense, and all the derivatives are understood in a weak sense.
For the integral identity to make sense, we will need $u^{m}$ in the Sobolev space
\[ L_{2}\big( 0,T:H_{0}^{1}(\Omega)\big).\]
Let’s unpack this notation a bit.
Recall that $L_2(\Omega)$ is a Banach space of Lebesgue measurable functions $u$ with finite norm
\[ \norm{u}_{L_2(\Omega)}= \left( \int_{\Omega} |u|^{2}dx\right)^{1/2}\]
It is an Hilbert space with inner product
\[ (u,v)=\int_{\Omega}uvdx \]
and
\[ H^{1}(\Omega)=\bk{u:\Omega \to \R,~u,\nabla u \in L_{2}(\Omega)}\]
is a Banach space with norm
\[ \norm{u}_{H^{1}(\Omega)}=\left( \int_{\Omega} |u|^{2} + \norm{\nabla u}^{2}\right)^{1/2}\]
Then $H_{0}^{1}(\Omega)$ is the linear subspace of elements of $H^{1}(\Omega)$ which are “zero” on the boundary of $\Omega$. To be more rigorous, $H_{0}^{1}(\Omega)$ is the closure of $C_{0}^{\infty}(\Omega)$ with respect to the $H^{1}(\Omega)$ norm.
Now, we can define the Sobolev space $L_{2}\big( 0,T:H_{0}^{1}(\Omega)\big)$ comprising functions mapping time into Bnach spaces:
\[ L_{2}\big( 0,T:H_{0}^{1}(\Omega)\big)=\bk{u=u(t):[0,T]\to H_{0}^{1}(\Omega)}\]
that is, points in this space are maps from the interval $[0,T]$ to the Banach space $H_{0}^{1}(\Omega)$. This space has a norm
\[ \norm{u}_{L_{2}\big( 0,T:H_{0}^{1}(\Omega)\big)} = \left( \int_{0}^{T} \norm{u(\cdot,t)}_{H^{1}(\Omega)}^{2}\,dt\right)^{1/2}=\left( \int_{0}^{T}\int_{\Omega} |u|^{2} + \norm{\nabla u}^{2}\,dx\,dt\right)^{1/2}\]
Finally note that for the right-hand side of \eqref{eq:pme-weak-int} to make sense, it is enough to assume $u_0\in L_1(\Omega)$, space of Lebesgue integrable functions with the norm
\[ \norm{u}_{L_1(\Omega)}=\int_\Omega|u|dx \]
We can now formulate the notion of the weak solution of the problem \eqref{eq:pme}–\eqref{eq:pme-bc}:
Definition
We say that a function $u=u(x,t)$ is a weak solution of the Dirichlet problem for the Nonlinear Diffusion Equation, \eqref{eq:pme}–\eqref{eq:pme-bc} if
- $u \geq 0$
- $u^{m} \in L_{2}(0,T;H^{1}_{0}(\Omega))$,
- $u$ satisfies the integral identity
\[\iint_{Q_{T}}\left( \nabla u^{m}\cdot \nabla \phi-u \phi_{t}\right)\,dx\,dt = \int_{\Omega}u_{0}(x)\phi(x,0)\,dx\]
for any $\phi \in C^{1}(\overline{Q}_{T})$ satisfying $\phi(x,T)=\phi\big\vert_{S_{T}}=0$.
This definition gives sense to ZKB solution, and one can check that $u_*$ is a weak solution of the relevant Dirichlet problem (recall that $\nabla u_*^m$ is a continuous function). How to make sure that the ZKB solution is the physically relevant solution, and in particular, there is no smooth classical solution of the nonlinear diffusion equation for the same data? To justify this we need to prove the uniqueness of the weak solution.
Uniqueness
Consider two solutions $u_{1},u_{2}$ in the sense of Definition 1. Their difference satisfies the integral identity
\begin{equation} \begin{gathered}\iint_{Q_{T}}\left[ \nabla \left( u_{1}^{m}-u_{2}^{m}\right)\cdot \nabla \phi – \left( u_{1}-u_{2}\right)\phi_{t}\right]\,dx\,dt=0,\\ \forall \phi \in C^{1}(\bar{Q}_{T})~\text{s.t.}~\phi(x,T)=\phi\big\vert_{S_{T}}=0\end{gathered}\label{eq:pme-uniq-1}\end{equation}
Due to the density of the space of smooth functions $C^1(\bar{Q}_{T})$ in a Sobolev space $H^1(Q_T)$, this integral identity holds for all $\phi \in H^{1}(Q_{T})$ with $\phi(x,T)=\phi\big\vert_{S_{T}}=0$. Set
\[ \eta(x,t)=\begin{cases}
\int_t^T u_{1}^{m}(x,s)-u_{2}^{m}(x,s)~ds & 0<t< T\\
0,& t\geq T
\end{cases}\]
We calculate
\begin{align*}
\eta_{t}&=-\left( u_{1}^{m} – u_{2}^{m}\right) \in L_{2}(Q_{T}),\\
\nabla \eta &= \int_{t}^{T}\left( \nabla u_{1}^{m} – \nabla u_{2}^{m}\right)\,ds \in L_{2}(Q_{T})
\end{align*}
so $\eta$ is an appropriate test function. Substituting it to \eqref{eq:pme-uniq-1} gives
\begin{align*}
0
&=\iint_{Q_{T}}\left[ \nabla \left( u_{1}^{m}-u_{2}^{m}\right)\cdot \int_{t}^{T}\left( \nabla u_{1}^{m} – \nabla u_{2}^{m}\right)\,ds + \left( u_{1}-u_{2}\right)\left( u_{1}^{m} – u_{2}^{m}\right)\right]\,dx\,dt\\
&=\iint_{Q_{T}}\left[ -\frac{1}{2}\dD{}{t}\left( \int_{t}^{T}\left( \nabla u_{1}^{m} – \nabla u_{2}^{m}\right)\,ds\right)^{2} + \left( u_{1}-u_{2}\right)\left( u_{1}^{m} – u_{2}^{m}\right)\right]\,dx\,dt\\
&=\frac{1}{2}\int_{\Omega}\left( \int_{0}^{T}\left( \nabla u_{1}^{m} – \nabla u_{2}^{m}\right)\,ds\right)^{2}\,dx + \iint_{Q_{T}}\left( u_{1}-u_{2}\right)\left( u_{1}^{m} – u_{2}^{m}\right)\,dx\,dt
\end{align*}
Note that both terms on the right-hand side are nonnegative, and accordingly both must be equal to zero. Hence, integrands must be equal to zero almost everywhere (a.e.). It follows that $u_{1}=u_{2}$ a.e. in $Q_{T}$. Thus we proved the following result:
Theorem
There exists at most one solution of the Dirichlet problem if $u_0\in L_1(\Omega)$.
In particular, it follows
Corollary
There is no classical solution for ZKB datum i.e. initial data of the form $u_*(x,t_{0})$ for fixed $t_{0}$.
References
- O.A. Oleinik, A.S. Kalashnikov, Yui-Lin Czou, The Cauchy problem and boundary problems for equations of the type of non-stationary filtration (Russian), Izv. Akad. Nauk SSSR. Ser. Mat. 22 (1958), 667–704.
- J.L. Vazquez, The Porous Medium Equation. Mathematical Theory., Oxford Mathematical Monographs. The Clarendon Press, Oxford University Press, Oxford, 2007. xxii+624 pp.