数学 旧 .com 迁移

偏微分方程:热方程 / 极值原理

从旧 .com 全量搬运的历史内容,来源路径:/math/课程/偏微分方程/chapters/热方程/极值原理/

迁移来源

极值原理 #

弱极值原理 #

Q={(x,t):0<x<l,0<tT}Q=\{(x,t):0<x<l,0<t\leqslant T\}.

tip

uC2,1(Q)C(Q)u\in C^{2,1}(Q)\cap C(\overline{Q}) 且满足 Lu=f0Lu=f\leqslant 0, 则 uuQ\overline{Q} 上的最大值必在 QQ 的抛物边界 Γ\Gamma 上达到, 即

maxQu(x,t)=maxΓu(x,t)\max\limits_{\overline{Q}} u(x,t)=\max\limits_{\Gamma}u(x,t)
note

f<0f<0, 用反证法容易得到矛盾.

f0f\leqslant 0, 取 v(x,t)=u(x,t)εtv(x,t)=u(x,t)-\varepsilon t, 则 Lv=Luε=f(x,t)ε<0Lv=Lu-\varepsilon=f(x,t)-\varepsilon<0, 于是

maxQu(x,t)=maxQ(v(x,t)+εt)maxQv(x,t)+εT\max\limits_{\overline{Q}}u(x,t)=\max\limits_{\overline{Q}}(v(x,t)+\varepsilon t)\leqslant \max\limits_{\overline{Q}} v(x,t)+\varepsilon T
tip

Lu=f0Lu=f\geqslant 0, 则最小值在边界取到

minQu=minΓu\min\limits_{\overline{Q}}u=\min\limits_{\Gamma}u
tip

u,vC2,1(Q)C(Q)u,v\in C^{2,1}(Q)\cap C(\overline{Q}), 且有 LuLvLu\leqslant Lv, uΓvΓu\big|_\Gamma\leqslant v\big|_\Gamma, 则在 Q\overline{Q}u(x,t)v(x,t)u(x,t)\leqslant v(x,t).

第一边值问题解的最大模估计 #

question

考虑第一边值问题

{Lu=uta2uxx=f,(x,tQ)ut=0=φ(x),(0xl)ux=0=g1(t),ux=l=g2(t),0tT\begin{equation} \begin{cases} Lu=u_t-a^2u_{xx}=f, & (x,t\in Q)\\ u\big|_{t=0}=\varphi(x), & (0\leqslant x\leqslant l)\\ u\big|_{x=0}=g_1(t), u\big|_{x=l}=g_2(t), & 0\leqslant t\leqslant T \end{cases} \end{equation}
tip

u(x,t)C2,1(Q)C(Q)u(x,t)\in C^{2,1}(Q)\cap C(\overline{Q}) 是问题 \eqref{热方程第一边值最大模} 的解, 则

maxQuFT+B,\max\limits_{\overline{Q}}|u|\leqslant FT+B,

其中

F=supQfB=max{sup[0,l]φ,sup[0,T]g1,sup[0,T]g2}\begin{aligned} F=\sup\limits_{Q}|f|\\ B=\max\left\lbrace\sup\limits_{[0,l]}|\varphi|, \sup\limits_{[0,T]}|g_1|,\sup\limits_{[0,T]}|g_2|\right\rbrace \end{aligned}
note

考虑辅助函数 w(x,t)=Ft+B±u(x,t)w(x,t)=Ft+B\pm u(x,t), 则有 Lw=F±f0Lw=F\pm f\geqslant 0.

且在边界上 wΓB+±u0w\big|_\Gamma\geqslant B+\pm u\geqslant 0. 由弱极值原理, 在 QQw(x,t)0w(x,t)\geqslant 0.

tip

问题 \eqref{热方程第一边值最大模} 在 C2,1(Q)C(Q)C^{2,1}(Q)\cap C(\overline{Q}) 中的解是唯一的.

tip

问题 \eqref{热方程第一边值最大模} 在 C2,1(Q)C(Q)C^{2,1}(Q)\cap C(\overline{Q}) 中的解连续依赖于非齐次项 ff, 初值 φ\varphi 与边值 g1,g2g_1,g_2.

第二, 三边值问题解的最大模估计 #

question

第二, 三边值问题可以统一写成

{Lu=uta2uxx=f,(x,tQ)ut=0=φ(x),(0xl)[ux+α(t)u]x=0=g1(t),0tT[ux+β(t)u]x=0=g2(t),0tT\begin{equation} \begin{cases} Lu=u_t-a^2u_{xx}=f, & (x,t\in Q)\\ u\big|_{t=0}=\varphi(x), & (0\leqslant x\leqslant l)\\ \left[-\frac{\partial u}{\partial x}+\alpha(t)u\right]_{x=0}=g_1(t), & 0\leqslant t\leqslant T \\ \left[\frac{\partial u}{\partial x}+\beta(t)u\right]_{x=0}=g_2(t), & 0\leqslant t\leqslant T \end{cases} \end{equation}

其中 α(t),β(t)0\alpha(t),\beta(t)\geqslant 0, 当 α(t)β(t)0\alpha(t)\equiv\beta(t)\equiv 0 时就是第二边值问题.

tip

u(x,t)C2,1(Q)C1,0(Q)u(x,t)\in C^{2,1}(Q)\cap C^{1,0}(\overline{ Q}) 满足

{Lu=uta2uxx0(x,t)Qut=00,0xl[ux+α(t)u]x=00,0tT[ux+β(t)u]x=l0,0tT\begin{cases} Lu=u_t-a^2u_{xx}\geqslant 0 & (x,t)\in Q\\ u\big|_{t=0}\geqslant 0, & 0\leqslant x\leqslant l\\ \left[-\frac{\partial u}{\partial x}+\alpha(t) u\right]_{x=0}\geqslant 0, & 0\leqslant t\leqslant T\\ \left[\frac{\partial u}{\partial x}+\beta(t)u\right]_{x=l}\geqslant 0, & 0\leqslant t\leqslant T \end{cases}

则在 Q\overline{Q}u(x,t)0u(x,t)\geqslant 0.

note

类似弱极值原理证明, 先设后两个条件是大于. 证明没有负的最小值.

由弱极值原理, 最小值必在边界取得, 设 P0P_0 是负的最小值, 那么由

uxP00,α(t)u(P0)0-\frac{\partial u}{\partial x}\big|_{P_0}\leqslant 0, \alpha(t)u(P_0)\leqslant 0

从而与假设矛盾. 故没有负的最小值.

再取辅助函数

v(x,t)=u(x,t)+ε[2a2t+(xl2)2]v(x,t)=u(x,t)+\varepsilon\left[2a^2t+\left(x-\frac l2\right)^2\right]

从而 v(x,t)0v(x,t)\geqslant 0, 故 u(x,t)ε(2a2T+l2/4)u(x,t)\geqslant-\varepsilon(2a^2T+l^2/4) 再由 ε\varepsilon 的任意性得到结论.

tip

u(x,t)C2,1(Q)C1,0(Q)u(x,t)\in C^{2,1}(Q)\cap C^{1,0}(\overline{Q}) 是问题 \eqref{热方程第三边值最大模} 的解, 则

u(x,t)C(F+B)|u(x,t)|\leqslant C(F+B)

其中常数 CC 只依赖于 a,l,Ta,l,T, 而

F=supQf,B=max{max[0,T]g1,max[0,T]g2,max[0,l]φ}\begin{aligned} F=\sup\limits_{Q}|f|,\\ B=\max\left\lbrace\max\limits_{[0,T]}|g_1|,\max\limits_{[0,T]}|g_2|,\max\limits_{[0,l]}|\varphi|\right\rbrace \end{aligned}
note

考虑辅助函数

w(x,t)=Ft+Bz(x,t)±u(x,t)w(x,t)=Ft+Bz(x,t)\pm u(x,t)

, 其中

z(x,t)=1+1l[2a2t+(xl2)2]z(x,t)=1+\frac 1l\left[2a^2t+\left(x-\frac l2\right)^2\right]

w(x,t)w(x,t) 满足引理条件, 从而

u(x,t)FT+Bz(x,t)FT+(1+2a2Tl+l4)B|u(x,t)|\leqslant FT+Bz(x,t)\leqslant FT+\left(1+\frac{2a^2T}{l}+\frac{l}{4}\right)B

C=max{T,1+2a2Tl+l4}C=\max\left\lbrace T,1+\dfrac{2a^2T}{l}+\dfrac{l}{4}\right\rbrace, 即可.

初值问题的最大模估计 #

question

在区域 Q={(x,t):xR,0<tT}Q=\{(x,t):x\in\mathbb{R},0<t\leqslant T\} 上考虑初值问题

{uta2uxx=f(x,t),(x,t)Qu(x,0)=φ(x),xR\begin{equation} \begin{cases} u_t-a^2u_{xx}=f(x,t), & (x,t)\in Q\\ u(x,0)=\varphi (x), & x\in \mathbb{R} \end{cases} \end{equation}
tip

u(x,t)C2,1(Q)C(Q)u(x,t)\in C^{2,1}(Q)\cap C(\overline{ Q}) 是初值问题 \eqref{热方程初值问题最大模} 的有界解, 则有估计

supQu(x,t)TsupQf(x,t)+sup(,)φ(x)\sup\limits_{Q}|u(x,t)|\leqslant T\sup\limits_{Q}|f(x,t)|+\sup\limits_{(-\infty,\infty)}|\varphi(x)|
note

对任意 L>0L>0, 考虑区域 QL={(x,t):x<L,0<tT}Q_L=\{(x,t):|x|<L,0<t\leqslant T\}.

F=supQf,ϕ=sup(,)φ.F=\sup\limits_{Q}|f|, \phi=\sup\limits_{(-\infty,\infty)}|\varphi|.supQu=K\sup\limits_{Q}|u|=K.

考虑辅助函数 w(x,t)=Ft+ϕ+vL(x,t)±u(x,t)w(x,t)=Ft+\phi+v_L(x,t)\pm u(x,t), 其中 vL(x,t)=KL2(x2+2a2t)v_L(x,t)=\frac{K}{L^2}(x^2+2a^2t).

容易验证边界值均非负, 故有极值原理知在 QLQ_Lw(x,t)0w(x,t)\geqslant 0. 故对任意的 (x0,t0)Q(x_0,t_0)\in Q, 总存在 LL, 可得 w(x0,t0)0w(x_0,t_0)\geqslant 0. 从而

u(x0,t0)Ft0+ϕ+KL2(x02+2a2t0)|u(x_0,t_0)|\leqslant Ft_0+\phi+\frac{K}{L^2}(x_0^2+2a^2t_0)

再令 LL\to \infty 则有

u(x0,t0)Ft0+ϕ|u(x_0,t_0)|\leqslant Ft_0+\phi

于是定理得证.

边值问题解的能量模估计 #

question

QT={(x,t):0<x<l,0<tT}Q_T=\{(x,t):0<x<l,0<t\leqslant T\}, 在 QTQ_T 上考虑混合问题

{uta2uxx=f,(x,t)QTu(x,0)=φ(x),0xl,u(0,t)=u(l,t)=0,0tT.\begin{equation} \begin{cases} u_t-a^2u_{xx}=f, & (x,t)\in Q_T\\ u(x,0)=\varphi(x), & 0\leqslant x\leqslant l,\\ u(0,t)=u(l,t)=0, & 0\leqslant t\leqslant T. \end{cases} \end{equation}
tip

uC2,1(QT)C1,0(QT)u\in C^{2,1}(Q_T)\cap C^{1,0}(\overline{Q_T}) 是问题 \eqref{边值问题解的能量模估计} 的解, 则有估计

sup0tT0lu2(x,t)dx+2a20T0l(ux)2dxdtM(0lφ2(x)dx+0T0lf2(x,t)dxdt)\sup\limits_{0\leqslant t\leqslant T}\int_0^l u^2(x,t)\text{d} x+2a^2\int_0^T\int_0^l \left(\frac{\partial u}{\partial x}\right)^2\text{d} x\text{d} t\leqslant M\left(\int_0^l \varphi^2(x)\text{d} x+\int_0^T\int_0^l f^2(x,t)\text{d} x\text{d} t\right)

其中 MM 只与 TT 有关.

question

题目 #

uC2,1(Q),utC2,1(Q)u\in C^{2,1}(\overline{Q}), u_t\in C^{2,1}(Q) 且满足以下定解问题

{utuxx=f(x,t),(x,t)Qut=0=φ(x)0xlux=0=ux=l=0,0tT\begin{cases} u_t-u_{xx}=f(x,t), & (x,t)\in Q\\ u\big|_{t=0} =\varphi(x) & 0\leqslant x\leqslant l\\ u\big|_{x=0}=u\big|_{x=l}=0, & 0\leqslant t\leqslant T \end{cases}

则有估计

maxQut(x,t)C[fC1(Q)]+φC[0,1]\max\limits_{\overline{ Q}}|u_t(x,t)|\leqslant C[\Vert f \Vert_{C^1(\overline{Q})}]+\Vert \varphi'' \Vert_{C[0,1]}

其中 CC 仅依赖于 TT.

note

v=utv=u_t, 则 vv 满足

{vtvxx=ft(x,t),(x,t)Qvt=0=f(x,t)+uxxt=0=f(x,0)+φ(x)0xlvx=0=vx=l=0,0tT\begin{cases} v_t-v_{xx}=f_t(x,t), & (x,t)\in Q\\ v\big|_{t=0}=f(x,t)+u_{xx}\big|_{t=0}=f(x,0)+\varphi''(x) & 0\leqslant x\leqslant l\\ v\big|_{x=0}=v\big|_{x=l}=0, & 0\leqslant t\leqslant T \end{cases}

从而由最大模原理知

supQv(x,t)TsupQft+sup[0,l]ε(x)+f(x,0)\sup\limits_{\overline{Q}}|v(x,t)|\leqslant T\sup\limits_{Q}|f_t|+\sup\limits_{[0,l]}|\varepsilon '' (x)+f(x,0)|

再进行放缩 supf(x,0)supf(x,t)\sup|f(x,0)|\leqslant\sup|f(x,t)|, 又 supft+supf=fC1\sup|f_t|+\sup|f|=\Vert f \Vert_{C^1}. 故取 C=T+1C=T+1 即可得到估计.

题目 #

uC1,0(Q)C2,1(Q)u\in C^{1,0}(\overline{Q})\cap C^{2,1}(Q) 且满足定解问题

{utuxx=0,(x,t)Qut=0=φ(x)0xlux=0=ux=l=0,0tT\begin{cases} u_t-u_{xx}=0, & (x,t)\in Q\\ u\big|_{t=0} =\varphi(x) & 0\leqslant x\leqslant l\\ u\big|_{x=0}=u\big|_{x=l}=0, & 0\leqslant t\leqslant T \end{cases}

证明

  • (1)
max(0,T)u(0,t)x,max(0,T)u(l,t)xC\max\limits_{(0,T)}\left|\frac{\partial u(0,t)}{\partial x}\right|,\max\limits_{(0,T)}\left|\frac{\partial u(l,t)}{\partial x}\right|\leqslant C

其中 CC 仅依赖于 φC1\Vert \varphi \Vert_{C^1}.

  • (2) 又设 uxC2,1(Q)u_x\in C^{2,1}(Q)
maxQu(x,t)xC~\max\limits_{\overline{Q}}\left|\frac{\partial u(x,t)}{\partial x}\right|\leqslant \widetilde{C}

其中 C~\widetilde{C} 也仅依赖于 φC1\Vert \varphi \Vert_{C^1}.

note

(1) 考虑 u(0,t)x=limx0u(x,t)u(0,t)x=limx0u(x,t)x\frac{\partial u(0,t)}{\partial x}=\lim\limits_{x\to 0}\frac{u(x,t)-u(0,t)}{x}=\lim\limits_{x\to 0}\frac{u(x,t)}{x}, 故只需证 u(x,t)Cx|u(x,t)|\leqslant Cx 即可.

则考虑辅助函数 w(x,t)=Cx±u(x,t)w(x,t)=Cx\pm u(x,t), 使用弱极值原理, 边界非负需使用拉格朗日中值定理, 结合 φ(0)=0\varphi(0)=0, 得到 Cx±φ(x)0Cx\pm \varphi(x)\geqslant 0. 从而 w(x,t)0w(x,t)\geqslant 0.

x=lx=l 则考虑 w(x,t)=C(lx)±u(x,t)w(x,t)=C(l-x)\pm u(x,t).

(2) 设 v=uxv=u_x, 则 vtvxx=fx=0v_t-v_{xx}=f_x=0, v(x,0)=φ(x)v(x,0)=\varphi'(x). 则由第一边值问题的最大模估计

maxQFT+B\max\limits_{\overline{Q}}\leqslant FT+B

其中 F=supf=0,B=supφ(x)φC1F=\sup|f|=0, B=\sup|\varphi'(x)|\leqslant\Vert \varphi \Vert_{C^1}.

题目 #

u,uxC2,1(Q)C1,0(Q)u,u_x\in C^{2,1}(Q)\cap C^{1,0}(\overline{Q}), uu 满足第三边值问题

{Lu=utuxx=f(x,t)(x,t)Qut=0=φ(x),0xl[ux+αu]x=0=g1(t),0tT[ux+βu]x=l=g2(t),0tT\begin{cases} Lu=u_t-u_{xx}=f(x,t) & (x,t)\in Q\\ u|_{t=0}=\varphi(x), & 0\leqslant x\leqslant l\\ [-\frac{\partial u}{\partial x}+\alpha u]_{x=0}=g_1(t), & 0\leqslant t\leqslant T\\ [\frac{\partial u}{\partial x}+\beta u]_{x=l}=g_2(t), & 0\leqslant t\leqslant T \end{cases}

其中 α,β0\alpha,\beta\geqslant 0, 给出 maxux\max|\frac{\partial u}{\partial x}| 的估计.

note

v=uxv=u_x, 首先由 uu 的第三边值问题可以给出估计

uC(F+B)|u|\leqslant C(F+B)

, F=supf,B=max{supφ,supg1,supg2}F=\sup|f|, B=\max\{\sup|\varphi|,\sup|g_1|,\sup|g_2|\}, C=max{T,1+2a2Tl+l4}C=\max\{T,1+\frac{2a^2T}{l}+\frac l4\}.

则考虑 vv 满足第一边值问题

{Lv=vtvxx=fx,(x,t)Qv(x,0)=φ(x),0xlv(0,t)=αu(0,t)g1(t),0tTv(l,t)=βu(l,t)g2(t),0tT\begin{cases} Lv=v_t-v_{xx}=f_x, & (x,t)\in Q\\ v(x,0)=\varphi'(x), & 0\leqslant x\leqslant l\\ v(0,t)=\alpha u(0,t)-g_1(t), & 0\leqslant t\leqslant T\\ v(l,t)=-\beta u(l,t)-g_2(t), & 0\leqslant t\leqslant T \end{cases}

从而有估计

vFT+B|v|\leqslant FT+B

其中 F=supfx,B=max{supφ,supu(0,t)g1(t),supu(l,t)g2(t)}F=\sup|f_x|, B=\max\{\sup|\varphi'|,\sup|u(0,t)-g_1(t)|,\sup|u(l,t)-g_2(t)|\}, 由 uu 的估计可得 BB 存在.

题目 #

Ql={0<x<l,0<tT}Q^l=\{0<x<l,0<t\leqslant T\}, 设 ulC2,1(Ql)C(Ql)u_l\in C^{2,1}(Q^l)\cap C(\overline{Q^l}) 是定解问题

{Lu=utuxx=0,(x,t)Qlult=0=0,0xlul(0,t)=g1(t),ul(l,t)=0,0tT\begin{cases} Lu=u_t-u_{xx}=0, & (x,t)\in Q^l\\ u_l\big|_{t=0}=0, & 0\leqslant x\leqslant l\\ u_l(0,t)=g_1(t), u_l(l,t)=0, & 0\leqslant t\leqslant T \end{cases}

的解, 其中 g1(t)0g_1(t)\geqslant0 证明: l1<l2l_1<l_2ul1(x,t)ul2(x,t),(x,t)Ql1u_{l_1}(x,t)\leqslant u_{l_2}(x,t),\quad (x,t)\in Q^{l_1}.

note

先证明 ulu_l 非负, 直接使用弱极值原理即可得到.

于是定义 w(x,t)=ul2(x,t)ul1(x,t)w(x,t)=u_{l_2}(x,t)-u_{l_1}(x,t), 在 Ql1Q^{l_1} 上.

那么再用一次弱极值原理, 对于 w(x,t)w(x,t) 的边界 w(0,t)=0,w(l1,t)=ul2(l1,t)0w(0,t)=0, w(l_1,t)=u_{l_2}(l_1,t)\geqslant 0, 从而 w(x,t)w(x,t) 非负. 证毕.

题目 #

uu 满足

Lu=uta2uxx+c(x,t)u0Lu=u_t-a^2u_{xx}+c(x,t)u\leqslant 0

其中 c(x,t)c(x,t) 有界, 且 c(x,t)0c(x,t)\geqslant 0. 证明如果 uuQ\overline{Q} 上存在非负最大值, 则必在边界上达到非负最大值.

note

按弱极值原理证明过程走一遍即可, 先考虑严格小于 0.

再考虑辅助函数 v(x,t)=u(x,t)+εeNtv(x,t)=u(x,t)+\varepsilon e^{-Nt}, 变为严格小于 0. 其中 N>c(x,t)N>|c(x,t)|, 由 c(x,t)c(x,t) 有界, 保证存在这样的 N. {{< /admonition >}}

讨论

评论

正在加载评论...