线性部分:古典解1-极值定理及其应用
1.定义
对于二阶的线性偏微分算子,一般有以下两种形式:散度型形式(divergence)
和非散度型:
这两种形式各有各的特点,其中散度型利于分部积分,非散度型利用证明古典解的极值原理.
在这里研究\(\Omega\)是\(\mathbb{R}^n\)中有界开集的情况.
在研究
这类方程时,我们一般会要求系数满足如下条件:
以及经常出现的如下条件:
[椭圆] :如果存在一个常数\(\theta>0\)使得:
\[\sum_{i, j=1}^n a^{i j}(\boldsymbol{x}) \xi_i \xi_j \geqslant \theta|\boldsymbol{\xi}|^2 \text { a.e. } \boldsymbol{x} \in \Omega,\forall \xi \in \mathbb{R}^n \]那么我们称偏微分算子\(L\)是\(\color{red}{严格椭圆}\)如果存在\(\lambda>0,\Lambda\ge 0\)使的:
\[\Lambda|\boldsymbol{\xi}|^2\ge \sum_{i, j=1}^n a^{i j}(\boldsymbol{x}) \xi_i \xi_j \geqslant \lambda|\boldsymbol{\xi}|^2 \text { a.e. } \boldsymbol{x} \in \Omega,\forall \xi \in \mathbb{R}^n \]则称\(L\)是一致椭圆的.
类似我们学习的三大方程一样,我们也研究椭圆型算子的初边值问题:
2.古典解的极值原理
下边我们证明非散度型算子的极值原理.
同时约定\(a^{ij},b^i,c\)是连续的.满足一致椭圆条件.
弱极值原理
[\(c=0\)的弱极值原理]
设 \(u \in C^2(\Omega) \cap C(\bar{\Omega})\) 且在 \(\Omega\) 上 \(c(\boldsymbol{x}) \equiv 0\).(1) 如果 \(L u \leqslant 0\) 在 \(\Omega\) 上成立, 则 \(\max\limits_{\bar{\Omega}} u=\max\limits _{\partial \Omega} u\);
(2) 如果 \(L u \geqslant 0\) 在 \(\Omega\) 上成立, 则 \(\min\limits _{\bar{\Omega}} u=\min\limits _{\partial \Omega} u\).
证明:\(\color{blue}{第一步}\).先考虑\(Lu < 0\)的情况,我们假设存在\(x_0\in \Omega\),使得其达到\(\bar{\Omega}\)上的最大值.断言\(x\in \partial \Omega\).不妨设\(x_0\)位于内点,因此根据古典的极值理论,我们知道\(Du(x_0)=0\).因此我们只需要对:
进行判断即可.我们知道Hessi矩阵\(D^2u(x_0)\le 0\).由于\(A\)是正定的,因此对于\(x_0\)点,存在正交矩阵\(O\)使得:
我们记:\(y=x+O(x-x_0)\),因此:
直接计算:
又因为同时有\(x=x_0+O'(y-x_0)\),因此:
综上:
因此\(Lu\ge 0\).这与\(Lu<0\)是矛盾的.因此断言成立,\(x\in \partial \Omega\).
\(\color{blue}{第二步:}\)我们讨论\(Lu\le 0\)的情况,我们对\(\varepsilon\)进行扰动,记:
我们计算\(u^{\varepsilon}\)的\(Lu\).直接计算有:
取\(\lambda\)充分大,那么上边的式子就可以小于0,因此根据第一步我们知道:
同时我们有:
令\(\varepsilon\to 0\)我们就可以得到结论.
\(\color{blue}{第三步:}\)对于(2)的情况,只需要将\(u\)换为\(-u\)然后利用(1)的结论就行了.
下边我们讨论\(c(x)\ge 0\)的弱极值原理,同时我们记\(u^+=\max\{u,0\},u^-=\min\{u,0\}\).
[\(c\ge 0\)的弱极值原理]
设 \(u \in C^2(\Omega) \cap C(\bar{\Omega})\) 且在 \(\Omega\) 上 \(c(\boldsymbol{x}) \ge 0\).(1) 如果 \(L u \leqslant 0\) 在 \(\Omega\) 上成立, 则 \(\max\limits _{\bar{\Omega}} u\le \max\limits _{\partial \Omega} u^+\);
(2) 如果 \(L u \geqslant 0\) 在 \(\Omega\) 上成立, 则 \(\min\limits _{\bar{\Omega}} u\ge \min\limits _{\partial \Omega} u^-\).
证明:仅对(1)进行证明.
第一步: 同样还是先证\(Lu<0\)的情况,我们设\(M=\max_{\bar{\Omega}}u\).如果\(M=0\)那么结论显然,因此我们证明\(M>0\)的情况.假设不然,那么必然在内点取得非负极大值,同上我们有:
因此矛盾.
第二步: 对\(Lu\le 0\),我们还是对\(u\)作扰动,令
同理计算\(Lu^{\varepsilon}\),取\(\lambda\)充分大,我们有:
按\(c=0\)的弱极值原理分析,我们就知道结论是成立的.
注意判断\(a^{11}>\theta\)
强最大值定理.
现在我们建立强最大值原理,同样的我们需要引入Hofp引理.
[Hofp引理]
设 \(u \in C^2(\Omega) \cap C(\bar{\Omega})\) 且在 \(\Omega\) 上 \(c(\boldsymbol{x}) \equiv 0\).且:\[Lu\le 0,x\in \Omega \]若存在一个点\(x_0\in \partial \Omega\),使得:
\[u(x_0)>u(x),\forall x\in \Omega \]以及\(\Omega\)满足内部球条件,那么:
(i):
\[\frac{\partial u}{\partial \nu}|_{x_0}>0 \](ii)如果\(c\ge 0\),且有\(u(x_0)\ge 0\).则(i)的结论同样成立.
证明:考虑\(c\ge 0\),不妨假设\(U=B(0,r)\),证明想法是类似的,引入一个\(v(x)\),使得:\(f(x)=u(x)+\varepsilon v(x)\)在\(|x|=r\)上达到最大值,此时必然有
进而得到
下边开始详细的证明,取
直接计算\(Lv\),就可以得到:
其中\(A=(a^{ij}),B=(b^i)\),我们取\(R=B(0,r)-B(0,r/2)\),因此当\(\lambda\)取得充分大时就有\(Lv\le 0\).
现在令\(f(x)=f(x)=u(x)+\varepsilon v(x)-u(x_0)\),可以看出\(Lf\le 0\),由于连续性因此因此\(u\)在\(|x|=r/2\)有最大值根据\(u(x_0)>u(x)\)可知,取\(\varepsilon\)充分小,那么\(f\)在\(|x|=r/2\)上小于0,而在\(|x|=r\)上显然大于等于\(u(x_0)\),根据数学分析的知识我们知道:
因此:
故命题得证.
利用Hofp引理我们证明强极大值定理.
[\(c=0\)的强极大值定理]
设 \(u \in C^2(\Omega) \cap C(\bar{\Omega})\) 且在 \(\Omega\) 上 \(c(\boldsymbol{x}) \equiv 0\).且\(\Omega\)是有界连通集,则有:(1) 如果在 \(\Omega\) 内有 \(L u \leqslant 0\) 且 \(u\) 在 \(\Omega\) 的内点取到 \(\bar{\Omega}\) 上的最大值, 则在 \(\Omega\) 内 \(u \equiv\) 常数;
(2) 如果在 \(\Omega\) 内有 \(L u \geqslant 0\) 且 \(u\) 在 \(\Omega\) 的内点取到 \(\bar{\Omega}\) 上的最小值, 则在 \(\Omega\) 内 \(u \equiv\) 常数.
我们记\(M=\max_{\bar{\Omega}}\),并且记\(C:=\{x\in \Omega|u(x)=M\}\),由连续性我们知道\(C\)是闭集,现在我们证明他还是开集.现在记:
因此\(\Omega=C\cup V\).\(V\)是开集,现在取\(y\in V\),使得\(dist(y,C) < dist(y,\partial \Omega)\),作以\(y\)球球心的内点位于\(V\)的最大球\(B\),因此存在\(x_0\in C\)使得\(x_0\in \partial B\),满足内部球条件,且满足Hofp引理条件,但是\(x_0\)是内点,因此\(Du=0\),但是根据Hofp引理\(\frac{\partial u}{\partial \nu} > 0\),因此矛盾,所以\(V=\varnothing\),所以\(\Omega=C\),即\(u\equiv M\).命题得证.
[\(c\ge 0\)的强极大值定理]
设 \(u \in C^2(\Omega) \cap C(\bar{\Omega})\), 在 \(\Omega\) 内 \(c \geqslant 0\), \(\Omega\) 是有界的连通开集.(1) 如果在 \(\Omega\) 内有 \(L u \leqslant 0\) 且 \(u\) 在 \(\Omega\) 的内点达到 \(u\) 在 \(\bar{\Omega}\) 上的非负最大值, 则 \(u \equiv M\);
(2) 如果在 \(\Omega\) 内有 \(L u \geqslant 0\) 且 \(u\) 在 \(\Omega\) 的内点达到 \(u\) 在 \(\bar{\Omega}\) 上的非正最小值, 则 \(u \equiv M\).
应用\(c\ge 0\)的Hofp引理即可.
3.先验估计:Drichlet边值条件
本节我们建立古典解的先验估计.
仍然设\(\Omega\)是\(\mathbb{R}^n\)中的有界区域,以及非散度型椭圆偏微分算子:
其中\(a_{ij},b_i,c\in C(\bar{\Omega})\),且:
以及:
[先验估计1] :设\(u\in C^2(\Omega)\cap C(\bar{\Omega})\)是方程:
\[\begin{cases}L u=f & \text { in } \Omega \\ u=\varphi & \text { on } \partial \Omega\end{cases} \]的古典解,其中\(f\in C(\bar{\Omega}),\varphi\in C(\partial \Omega)\)如果\(c\le 0\),则:
\[|u(x)| \leq \max _{\partial \Omega}|\varphi|+C(\lambda,\Lambda,\mathrm{diam}\Omega) \max _{\Omega}|f| \]
证明: 现在我们想要利用最大值原理来证明该定理,我们的想法是找到一个函数\(w\),使得:
如此一来我们就可以利用最大值原理,得到\(w+u\)的最大值,我们相信如果对\(w\)取特别的值,就可以得到上述的先验估计.令\(F=\max_{ \Omega}|f|,\Phi=\max_{\partial \Omega}|\varphi|\),为了能够让上述的构造更加方便,我们直接令:
由于\(\Omega\)是有界的,因此我们不妨设\(\Omega\subset \{0 < x_1 < d\}\)中,令\(w=\Phi=\left(e^{\alpha d}-e^{\alpha x_1}\right)F,\alpha\)待定.直接计算:
注意到\(\alpha^2\lambda+b_1\alpha\)是关于\(\alpha\)的二次方程,因此\(\alpha\)充分大时就可以使得其大于1.
故这里的\(w\)满足我们想要的构造.利用最大值原理我们就知道:
因此:
于是定理得证.
4.先验估计:Robin和Nueman边值条件
下边考虑Robin和Nueman边值问题的先验估计.
[先验估计2] :设\(u\in C^2(\Omega)\cap C^1(\bar{\Omega})\)是方程:
\[\begin{cases} Lu=f,x\in \Omega\\ \frac{\partial u}{\partial \boldsymbol{n}}+\alpha(x)u=\varphi,x\in \partial\Omega \end{cases} \]其中\(\boldsymbol{n}\)是\(\partial\Omega\)的单位外法方向.如果\(c(x)\le 0x\in \Omega,\alpha(x)\ge \alpha_0>0,x\in \partial\Omega\),则:
\[|u(x)| \leq C(\lambda,\Lambda,\alpha_0,\mathrm{diam}\Omega)\left\{\max _{\partial \Omega}|\varphi|+\max _{\Omega}|f|\right\},x\in \Omega \]
证明: Case1:我们先考虑\(c(x)\le -c_0<0\)的情形.我们证明:
证明想法是类似的,构造\(v\)使得其和满足Hofp最大值定理的条件.令\(v=\frac{1}{c_0}F+\frac{1}{\alpha_0}\Phi\pm u\).则我们有:
如果\(v\)在\(\bar{\Omega}\)中负的最小值,则其在边界上可达,根据Hofp极值定理,我们知道存在\(x_0\)使得\(\frac{\partial v}{\partial n}(x_0)\le 0\),这意味着:
矛盾.因此\(v\ge 0,x\in \Omega\).故我们得到了:
Case2:考虑一般的情况,\(c\le 0\).我们考虑辅助函数\(u(x)=z(x)w(x)\),其中\(z\)是一个正的待定的函数,我们直接计算\(Lu\).就得到了:
其中\(B_i=\frac{1}{2}(a_{ij}+a_{ji})D_jz+b_i\),现在我们希望\(z>0,x\in \bar{\Omega}\),使得:
如此一来我们就可以利用\(c(x)<-c_0\)的情况.为此我们仍然设\(\Omega\subset \{0<x_1<d\}\),并令\(z(x)=A+e^{\beta d}-e^{\beta x_1}\)其中\(A,\beta>0\)待定.计算可得:
这里我们选择\(\beta\)充分大使得\(\beta^2a_{11}+\beta b_1\ge 1\),则我们有:
这里我们选择\(A\)充分大使得上式成立.
然后利用\(c(x)\le -c_0\)的情况即可得证.、
参考文献
- 韩青, 林芳华. Elliptic Partial Differential Equations[M]. AMS, 2011
- Evans. Partial differential equations[M]. 2010.
- 王术. Sobolev 空间和偏微分方程引论 [M]. 科学出版社, 2009