数学 偏微分方程

偏微分方程:第 5 次作业

整理热方程的最大值原理、唯一性与能量估计,以及典型区域 Green 函数练习。

本页目录10 节

20251126 #

红圈数字 12

3-16 #

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) 是定解问题

{ult2ulx2=0,(x,t)Ql,ult=0=0,0xl,ul(0,t)=g1(t), ul(l,t)=0,0tT\begin{cases} \dfrac{\partial u_l}{\partial t}-\dfrac{\partial^2 u_l}{\partial x^2}=0, & (x,t)\in Q^l,\\ u_l|_{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)\geqslant 0. 证明 ul(x,t)u_l(x,t) 关于 ll 是递增的, 即对于 l1<l2l_1<l_2,

ul1(x,t)ul2(x,t),(x,t)Ql1,u_{l_1}(x,t)\leqslant u_{l_2}(x,t),\quad (x,t)\in Q^{l_1},

并给出物理解释.

答案

对每个 l>0l>0, ulu_l 满足

{utuxx=0,0<x<l, 0<tT,ul(x,0)=0,0xl,ul(0,t)=g1(t)0,ul(l,t)=0,0tT. \begin{cases} u_t-u_{xx}=0,& 0<x<l,\ 0<t\leqslant T,\\[2mm] u_l(x,0)=0,& 0\leqslant x\leqslant l,\\[2mm] u_l(0,t)=g_1(t)\ge0,\quad u_l(l,t)=0,& 0\leqslant t\leqslant T. \end{cases}

v=ulv=-u_l, 则 vtvxx=(utuxx)=0v_t-v_{xx}=-(u_t-u_{xx})=0, 且在抛物边界上

v(x,0)=ul(x,0)=0,v(0,t)=g1(t)0,v(l,t)=ul(l,t)=0. v(x,0)=-u_l(x,0)=0,\quad v(0,t)=-g_1(t)\le0,\quad v(l,t)=-u_l(l,t)=0.

因此 vv 在边界上最大值为 00, 由最大值原理知 v0v\le0, 即 ul0u_l\ge0Ql\overline{Q^l}.

0<l1<l20<l_1<l_2, 设对应解为 ul1,ul2u_{l_1},u_{l_2}. 在区域 Ql1={0<x<l1,0<tT}Q^{l_1}=\{0<x<l_1,0<t\leqslant T\} 上考虑 w(x,t):=ul2(x,t)ul1(x,t).w(x,t):=u_{l_2}(x,t)-u_{l_1}(x,t). 两者都满足 utuxx=0u_t-u_{xx}=0, 故 wtwxx=0w_t-w_{xx}=0 ww 在抛物边界上的值:

w(x,0)=ul2(x,0)ul1(x,0)=0,0xl1, w(x,0)=u_{l_2}(x,0)-u_{l_1}(x,0)=0,\quad 0\leqslant x\leqslant l_1, w(0,t)=ul2(0,t)ul1(0,t)=g1(t)g1(t)=0, w(0,t)=u_{l_2}(0,t)-u_{l_1}(0,t)=g_1(t)-g_1(t)=0, w(l1,t)=ul2(l1,t)ul1(l1,t)=ul2(l1,t)0=ul2(l1,t)0, w(l_1,t)=u_{l_2}(l_1,t)-u_{l_1}(l_1,t)=u_{l_2}(l_1,t)-0=u_{l_2}(l_1,t)\ge0,

于是 ww 在抛物边界上满足 w0w\ge0, 且在内部满足 wtwxx=0w_t-w_{xx}=0. 再由最大值原理得 w(x,t)0,(x,t)Ql1,w(x,t)\ge0,\quad (x,t)\in \overline{Q^{l_1}},ul1(x,t)ul2(x,t),(x,t)Ql1.u_{l_1}(x,t)\le u_{l_2}(x,t),\quad (x,t)\in Q^{l_1}. 这就证明了 ulu_l 关于 ll 的递增性.

物理解释:

可以把问题看作一根长度为 ll 的均匀细杆, 其初温为 00, 左端 x=0x=0 的温度被控制为非负的给定函数 g1(t)g_1(t), 右端 x=lx=l 一直保持在 00^\circ. 当把杆从长度 l1l_1 延长到 l2>l1l_2>l_1 时, 对同一内部点 x(0,l1)x\in(0,l_1) 来说, 右端这个“冷端”被移得更远, 热量更不容易通过右端散失掉, 因此该点的温度在任何时刻都不会比短杆情况更低, 即

ul1(x,t)ul2(x,t). u_{l_1}(x,t)\le u_{l_2}(x,t).

3-18 #

uC(Q)C2,1(Q)u\in C(\overline{Q})\cap C^{2,1}(Q) 且满足

Luuta2uxx+c(x,t)u0,(x,t)Q,Lu\equiv u_t-a^2u_{xx}+c(x,t)u\leqslant0,\quad (x,t)\in Q,

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

答案

Q=(0,l)×(0,T]Q=(0,l)\times(0,T], a2>0a^2>0, 其抛物边界为

Γ=([0,l]×{0})({0,l}×[0,T]).\Gamma=([0,l]\times\{0\})\cup(\{0,l\}\times[0,T]).

注意顶面 t=Tt=T 不属于抛物边界。若函数在 0<x<l0<x<l, 0<tT0<t\leqslant T 处达到全局最大值, 则空间二阶导数不大于零, 时间导数满足 vt0v_t\geqslant0;在 t=Tt=T 处使用左侧时间导数即可。

N>supQcN>\sup_Q c, 对任意 ε>0\varepsilon>0 定义严格扰动

vε(x,t)=u(x,t)+εeNt.v_\varepsilon(x,t)=u(x,t)+\varepsilon e^{-Nt}.

直接求导得到

Lvε=Lu+ε(cN)eNt<0.Lv_\varepsilon =Lu+\varepsilon(c-N)e^{-Nt}<0.

假如 vεv_\varepsilonQQ 中达到正的全局最大值, 则在该点

(vε)t0,(vε)xx0,cvε0,(v_\varepsilon)_t\geqslant0,\qquad (v_\varepsilon)_{xx}\leqslant0,\qquad c\,v_\varepsilon\geqslant0,

从而 Lvε0Lv_\varepsilon\geqslant0, 与严格不等式矛盾。因此 vεv_\varepsilon 的正最大值只能在 Γ\Gamma 上取得。

M=maxQu0M=\max_{\overline Q}u\geqslant0。在 uu 的最大值点有 vε>0v_\varepsilon>0, 故 vεv_\varepsilon 确有正的最大值, 并且

MmaxQvε=maxΓvεmaxΓu+ε.M\leqslant\max_{\overline Q}v_\varepsilon =\max_\Gamma v_\varepsilon \leqslant\max_\Gamma u+\varepsilon.

ε0\varepsilon\downarrow0, 得 MmaxΓuM\leqslant\max_\Gamma u。反向不等式由 ΓQ\Gamma\subset\overline Q 直接成立, 所以

maxQu=maxΓu.\max_{\overline Q}u=\max_\Gamma u.

这个论证也包含 M=0M=0 的情形。

3-21 #

证明半无界问题

{uta2uxx=f(x,t),0<x<,t>0,ut=0=φ(x),0x<,ux=0=μ(t),t0.\begin{cases} u_t-a^2u_{xx}=f(x,t),& 0<x<\infty,t>0,\\ u|_{t=0}=\varphi(x), & 0\leqslant x<\infty,\\ u|_{x=0}=\mu(t), & t\geqslant 0. \end{cases}

的有界解是唯一的.

答案

u1,u2u_1,u_2 是两个有界经典解, 在初始面和边界上连续。令 w=u1u2w=u_1-u_2, 则

{wta2wxx=0,x>0, t>0,w(x,0)=0,x0,w(0,t)=0,t0.\begin{cases} w_t-a^2w_{xx}=0, & x>0,\ t>0,\\ w(x,0)=0, & x\geqslant0,\\ w(0,t)=0, & t\geqslant0. \end{cases}

固定任意 T>0T>0。由有界性, 存在 MT<M_T<\infty 使得 w(x,t)MT|w(x,t)|\leqslant M_T 对所有 x0x\geqslant0, 0tT0\leqslant t\leqslant T 成立。事实上, 每个有限时间条带上的一致有界性已经足够。

ε>0\varepsilon>0 构造屏障

bε(x,t)=ε(1+x2+(2a2+1)t),(bε)ta2(bε)xx=ε>0.b_\varepsilon(x,t)=\varepsilon\bigl(1+x^2+(2a^2+1)t\bigr), \qquad (b_\varepsilon)_t-a^2(b_\varepsilon)_{xx}=\varepsilon>0.

在有限柱体 QR=(0,R)×(0,T]Q_R=(0,R)\times(0,T] 中分别考虑 v±=±wbεv_\pm=\pm w-b_\varepsilon。它们满足

(v±)ta2(v±)xx=ε<0.(v_\pm)_t-a^2(v_\pm)_{xx}=-\varepsilon<0.

t=0t=0x=0x=0 上, v±0v_\pm\leqslant0。只要取 RR 足够大, 使得 ε(1+R2)MT\varepsilon(1+R^2)\geqslant M_T, 就有

v±(R,t)MTε(1+R2)0,0tT.v_\pm(R,t)\leqslant M_T-\varepsilon(1+R^2)\leqslant0, \qquad 0\leqslant t\leqslant T.

由有限柱体上的最大值原理, v±0v_\pm\leqslant0QR\overline{Q_R}。对于任意固定的 (x,t)(x,t), 取还满足 R>xR>x 的上述 RR, 得到

w(x,t)ε(1+x2+(2a2+1)t).|w(x,t)|\leqslant\varepsilon\bigl(1+x^2+(2a^2+1)t\bigr).

ε0\varepsilon\downarrow0, 可得 w(x,t)=0w(x,t)=0。由于 TT 任意, w0w\equiv0, 即 u1=u2u_1=u_2。整个比较只在有限区间进行, 不需要假设 wxw_x 在无穷远衰减。

3-22 #

u(x,t)C2,1(Q)u(x,t)\in C^{2,1}(\overline{Q}) 是问题

{utuxx=f,(x,t)Q,u(x,0)=φ(x),0xl,u(0,t)=u(l,t)=0,0tT\begin{cases} u_{t}-u_{xx}=f, & (x,t)\in Q,\\ u(x,0)=\varphi(x), & 0\leqslant x\leqslant l,\\ u(0,t)=u(l,t)=0, & 0\leqslant t\leqslant T \end{cases}

的解, 证明 uu 满足以下估计

sup0tT0l(ux)2dx+0T0l(ut)2dxdtM[0l(φ(x))2dx+0T0lf2(x,t)dxdt],\begin{aligned} & \sup\limits_{0\leqslant t\leqslant T} \int_0^l (u_x)^2\text{d} x + \int_0^T\int_0^l (u_t)^2\text{d} x\text{d} t\\ \leqslant & M\left[\int_0^l (\varphi'(x))^2\text{d} x + \int_0^T\int_0^l f^2(x,t)\text{d} x\text{d} t\right], \end{aligned}

其中 MM 只依赖于 T,lT,l.

答案

对方程两边乘以 utu_t, 在 (0,l)(0,l) 上积分, 得 0lut2dx0luxxutdx=0lfutdx.\int_0^l u_t^2\,dx-\int_0^l u_{xx}u_t\,dx=\int_0^l f\,u_t\,dx. 对第二项分部积分:

0luxxutdx=[uxut]0l+0luxutxdx. -\int_0^l u_{xx}u_t\,dx = -\big[u_x u_t\big]_{0}^{l}+\int_0^l u_x u_{tx}\,dx.

由于边界条件 u(0,t)=u(l,t)=0u(0,t)=u(l,t)=0 与时间无关, 对 tt 求导得 ut(0,t)=ut(l,t)=0,u_t(0,t)=u_t(l,t)=0, 故边界项 [uxut]0l=0[u_x u_t]_{0}^{l}=0, 于是

0luxxutdx=0luxuxtdx=12ddt0lux2dx. -\int_0^l u_{xx}u_t\,dx=\int_0^l u_x u_{xt}\,dx =\frac12\frac{d}{dt}\int_0^l u_x^2\,dx.

代回原式:

0lut2dx+12ddt0lux2dx=0lfutdx. \int_0^l u_t^2\,dx + \frac12\frac{d}{dt}\int_0^l u_x^2\,dx = \int_0^l f\,u_t\,dx.

应用 Cauchy—Schwarz 与 Young 不等式:

0lfutdx120lut2dx+120lf2dx, \left|\int_0^l f\,u_t\,dx\right| \le \frac12\int_0^l u_t^2\,dx + \frac12\int_0^l f^2\,dx,

从而

12ddt0lux2dx+0lut2dx120lut2dx+120lf2dx, \frac12\frac{d}{dt}\int_0^l u_x^2\,dx + \int_0^l u_t^2\,dx \le \frac12\int_0^l u_t^2\,dx + \frac12\int_0^l f^2\,dx,

ddt0lux2dx+0lut2dx0lf2dx. \frac{d}{dt}\int_0^l u_x^2\,dx + \int_0^l u_t^2\,dx \le \int_0^l f^2\,dx.

tt00 积分到任意 τ(0,T]\tau\in(0,T], 得到

0lux2(x,τ)dx0lux2(x,0)dx+0τ0lut2dxdt0τ0lf2dxdt. \int_0^l u_x^2(x,\tau)\,dx - \int_0^l u_x^2(x,0)\,dx + \int_0^\tau\int_0^l u_t^2\,dx\,dt \le \int_0^\tau\int_0^l f^2\,dx\,dt.

ux(x,0)=φ(x)u_x(x,0)=\varphi'(x), 并放宽上限 τT\tau\leqslant T, 得

0lux2(x,τ)dx+0τ0lut2dxdt0lφ(x)2dx+0T0lf2dxdt. \int_0^l u_x^2(x,\tau)\,dx + \int_0^\tau\int_0^l u_t^2\,dx\,dt \le \int_0^l \varphi'(x)^2\,dx + \int_0^T\int_0^l f^2\,dx\,dt.

A=0lφ(x)2dx+0T0lf2dxdt.A=\int_0^l\varphi'(x)^2\,dx+\int_0^T\int_0^l f^2\,dx\,dt.

由上式分别舍去两个非负项, 可得

sup0tT0lux2(x,t)dxA,0T0lut2dxdtA.\sup_{0\leqslant t\leqslant T}\int_0^l u_x^2(x,t)\,dx\leqslant A, \qquad \int_0^T\int_0^l u_t^2\,dx\,dt\leqslant A.

将两个估计相加, 得到

sup0tT0lux2(x,t)dx+0T0lut2(x,t)dxdtM[0lφ(x)2dx+0T0lf2dxdt], \sup_{0\leqslant t\leqslant T}\int_0^l u_x^2(x,t)\,dx +\int_0^T\int_0^l u_t^2(x,t)\,dx\,dt \le M\left[\int_0^l\varphi'(x)^2\,dx + \int_0^T\int_0^l f^2\,dx\,dt\right],

其中可取 M=2M=2, 与初值、源项无关, 也不依赖 T,lT,l。这里不能把随 τ\tau 变化的积分上限直接换成 TT 后仍取同一个上界。

4-17 #

求边值问题

{Δu=f(x,y),(x,y)Ω,uΩ=φ\begin{cases} -\Delta u =f(x,y), & (x,y)\in \Omega,\\ u|_{\partial\Omega}=\varphi & \end{cases}

的 Green 函数, 其中

  • (1) Ω\Omega 是上半平面;
  • (2) Ω\Omega 是第一象限;
答案
  • (1)
G(x,y;ξ1,ξ2)=12πln(xξ1)2+(yξ2)2(xξ1)2+(y+ξ2)2,y>0. G(x,y;\xi_1,\xi_2) =-\frac{1}{2\pi} \ln\frac{\sqrt{(x-\xi_1)^2+(y-\xi_2)^2}} {\sqrt{(x-\xi_1)^2+(y+\xi_2)^2}}, \qquad y>0.
  • (2)
G(x,y;ξ1,ξ2)=12π[ln(xξ1)2+(yξ2)2ln(xξ1)2+(y+ξ2)2ln(x+ξ1)2+(yξ2)2+ln(x+ξ1)2+(y+ξ2)2]. \begin{aligned} G(x,y;\xi_1,\xi_2) &=-\frac{1}{2\pi}\big[ \ln\sqrt{(x-\xi_1)^2+(y-\xi_2)^2} -\ln\sqrt{(x-\xi_1)^2+(y+\xi_2)^2}\\ &-\ln\sqrt{(x+\xi_1)^2+(y-\xi_2)^2} +\ln\sqrt{(x+\xi_1)^2+(y+\xi_2)^2} \big]. \end{aligned}

4-19 #

B+(R)={(x,y)x2+y2<R2,y>0}B^+(R)=\{(x,y)|x^2+y^2<R^2,y>0\}, 求定解问题

{Δu=f(x,y),(x,y)B+(R),uB+(R){y>0}=φ(x,y),uyy=0=ψ(x,0),RxR\begin{cases} -\Delta u = f(x,y), & (x,y)\in B^+(R),\\ u|_{\partial B^+(R)\cap\{y>0\}}=\varphi(x,y), & \\ u_y|_{y=0}=\psi(x,0), & -R\leqslant x\leqslant R \end{cases}

的 Green 函数. 如果 uC1(B+(R))C2(B+(R))u\in C^1(\overline{B^+}(R))\cap C^2(B^+(R)) 是上述问题的解, 试给出解的表达式.

答案

D=B+(R)D=B^+(R), 并将边界分成开圆弧与开直径

ΓD=B(R){y>0},ΓN={(x,0):x<R}.\Gamma_D=\partial B(R)\cap\{y>0\},\qquad \Gamma_N=\{(x,0):|x|<R\}.

两端点不影响边界积分。以下固定采用 ΔG=δ-\Delta G=\delta 的约定, nn 始终表示区域的外法向。

圆盘的 Dirichlet Green 函数 #

用复数记号 z=x+iyz=x+iy, ζ=ξ+iη\zeta=\xi+i\eta, 定义

GR(z,ζ)=12πlnR(zζ)R2zζ,z,ζ<R,zζ.G_R(z,\zeta) =-\frac{1}{2\pi}\ln\left| \frac{R(z-\zeta)}{R^2-z\overline\zeta} \right|, \qquad |z|,|\zeta|<R,\quad z\ne\zeta.

分母在圆盘内不为零, 所以对数修正项关于 zz 调和, 而奇异部分为 12πlnzζ-\frac1{2\pi}\ln|z-\zeta|。因此

ΔzGR(z,ζ)=δζ.-\Delta_zG_R(z,\zeta)=\delta_\zeta.

z=R|z|=R 时,

R2zζ2=R2zζ2,|R^2-z\overline\zeta|^2=R^2|z-\zeta|^2,

GR=0G_R=0 于圆周。这个表达式也适用于 ζ=0\zeta=0, 此时 GR(z,0)=12πln(z/R)G_R(z,0)=-\frac1{2\pi}\ln(|z|/R)

半圆盘的混合 Green 函数 #

设极点 ζ=ξ+iηD\zeta=\xi+i\eta\in D, 其镜像为 ζ\overline\zeta。取同号镜像

G+(z,ζ)=GR(z,ζ)+GR(z,ζ).G^+(z,\zeta)=G_R(z,\zeta)+G_R(z,\overline\zeta).

镜像极点位于下半圆盘, 因而第二项在 DD 内调和, 并有

ΔzG+(z,ζ)=δζ于 D.-\Delta_zG^+(z,\zeta)=\delta_\zeta\quad\text{于 }D.

圆弧上两项均为零, 所以 G+ΓD=0G^+|_{\Gamma_D}=0。由

GR(z,ζ)=GR(z,ζ)G_R(\overline z,\zeta)=G_R(z,\overline\zeta)

可知 G+(x+iy,ζ)G^+(x+iy,\zeta) 关于 yy 为偶函数。因此在直径上

G+yy=0=0,G+nΓN=0.\left.\frac{\partial G^+}{\partial y}\right|_{y=0}=0, \qquad \left.\frac{\partial G^+}{\partial n}\right|_{\Gamma_N}=0.

这就得到所求的 Green 函数。由显式公式还可验证 G+(z,ζ)=G+(ζ,z)G^+(z,\zeta)=G^+(\zeta,z)

解的表示 #

记求值点 P=(ξ,η)DP=(\xi,\eta)\in D, 积分变量 Q=(x,y)Q=(x,y)。在题设经典解下, 对去掉极点附近小圆盘的区域应用 Green 第二恒等式, 再令小圆半径趋于零, 得

u(P)=DG+(P,Q)f(Q)dQ+D(G+(P,Q)un(Q)u(Q)G+nQ(P,Q))dsQ.\begin{aligned} u(P)&=\int_DG^+(P,Q)f(Q)\,dQ\\ &\quad+\int_{\partial D} \left(G^+(P,Q)\frac{\partial u}{\partial n}(Q) -u(Q)\frac{\partial G^+}{\partial n_Q}(P,Q)\right)\,ds_Q. \end{aligned}

若仅由题设的内部 C2C^2 与闭域 C1C^1 正则性出发, 体积分可先按内部截断区域的极限理解;当 ff 连续到边界时, 它就是通常的面积积分。

ΓD\Gamma_Du=φu=\varphi, G+=0G^+=0;在 ΓN\Gamma_Nn=(0,1)n=(0,-1), 故

un=uy=ψ,G+n=0.\frac{\partial u}{\partial n}=-u_y=-\psi, \qquad \frac{\partial G^+}{\partial n}=0.

代入得到

u(P)=DG+(P,Q)f(Q)dQΓDφ(Q)G+nQ(P,Q)dsQRRG+((ξ,η);(x,0))ψ(x,0)dx.\begin{aligned} u(P)&=\int_DG^+(P,Q)f(Q)\,dQ -\int_{\Gamma_D}\varphi(Q)\frac{\partial G^+}{\partial n_Q}(P,Q)\,ds_Q\\ &\quad-\int_{-R}^{R}G^+\bigl((\xi,\eta);(x,0)\bigr)\psi(x,0)\,dx. \end{aligned}

两个边界项都取负号:圆弧项来自 unG+-u\,\partial_nG^+, 直径项来自 nu=ψ\partial_nu=-\psi

讨论

评论

正在加载评论…

输入关键词开始搜索。