泛函 Φ(u) 驻点所满足的偏微分方程与自然边界条件
题目
问题13:写出泛函
Φ(u)=∬D(ux2−2uxuy+2uy2)dxdy+∫Γ (u2+2g(x)u)dx\Phi(u)=\iint_{D}\left(u_{x}^{2}-2u_{x}u_{y}+2u_{y}^{2}\right)dxdy+\int_{\Gamma }\!\left(u^{2}+2g(x)u\right)dxΦ(u)=∬D(ux2−2uxuy+2uy2)dxdy+∫Γ(u2+2g(x)u)dx
的驻点所满足的方程和边界条件,其中 D={(x,y) :0<x<2, 0<y<1}D=\{(x,y)\colon 0<x<2,\,0<y<1\}D={(x,y):0<x<2,0<y<1},Γ={(x,y) :y=0,0<x<1}\Gamma=\{(x,y)\colon y=0,0<x<1\}Γ={(x,y):y=0,0<x<1},且满足边界条件
u∣y=1=u∣x=0=u∣x=2=u∣y=0, 1<x<2=0.u|_{y=1}=u|_{x=0}=u|_{x=2}=u|_{y=0,\,1<x<2}=0.u∣y=1=u∣x=0=u∣x=2=u∣y=0,1<x<2=0.
解答
为了求泛函 Φ(u)\Phi(u)Φ(u) 的驻点,使用变分法。计算一阶变分 δΦ\delta \PhiδΦ 并令其为零,得到欧拉-拉格朗日方程和自然边界条件。
欧拉-拉格朗日方程
泛函中的二重积分部分为 ∬DF(ux,uy) dxdy\iint_D F(u_x, u_y) \, dx dy∬DF(ux,uy)dxdy,其中 F=ux2−2uxuy+2uy2F = u_x^2 - 2 u_x u_y + 2 u_y^2F=ux2−2uxuy+2uy2。欧拉-拉格朗日方程为:
∂∂x(∂F∂ux)+∂∂y(∂F∂uy)=0
\frac{\partial}{\partial x} \left( \frac{\partial F}{\partial u_x} \right) + \frac{\partial}{\partial y} \left( \frac{\partial F}{\partial u_y} \right) = 0
∂x∂(∂ux∂F)+∂y∂(∂uy∂F)=0
计算偏导数:
∂F∂ux=2ux−2uy,∂F∂uy=−2ux+4uy
\frac{\partial F}{\partial u_x} = 2 u_x - 2 u_y, \quad \frac{\partial F}{\partial u_y} = -2 u_x + 4 u_y
∂ux∂F=2ux−2uy,∂uy∂F=−2ux+4uy
代入方程:
∂∂x(2ux−2uy)+∂∂y(−2ux+4uy)=0
\frac{\partial}{\partial x} (2 u_x - 2 u_y) + \frac{\partial}{\partial y} (-2 u_x + 4 u_y) = 0
∂x∂(2ux−2uy)+∂y∂(−2ux+4uy)=0
计算二阶导数:
2uxx−2uyx−2uxy+4uyy=0
2 u_{xx} - 2 u_{yx} - 2 u_{xy} + 4 u_{yy} = 0
2uxx−2uyx−2uxy+4uyy=0
假设 uuu 光滑,uxy=uyxu_{xy} = u_{yx}uxy=uyx,化简得:
2uxx−4uxy+4uyy=0
2 u_{xx} - 4 u_{xy} + 4 u_{yy} = 0
2uxx−4uxy+4uyy=0
除以2:
uxx−2uxy+2uyy=0
u_{xx} - 2 u_{xy} + 2 u_{yy} = 0
uxx−2uxy+2uyy=0
因此,在区域 DDD 内部,驻点满足方程:
uxx−2uxy+2uyy=0
u_{xx} - 2 u_{xy} + 2 u_{yy} = 0
uxx−2uxy+2uyy=0
边界条件
给定的边界条件为:
- u∣y=1=0u|_{y=1} = 0u∣y=1=0(上边界)
- u∣x=0=0u|_{x=0} = 0u∣x=0=0(左边界)
- u∣x=2=0u|_{x=2} = 0u∣x=2=0(右边界)
- u∣y=0,1<x<2=0u|_{y=0, 1<x<2} = 0u∣y=0,1<x<2=0(下边界部分)
在边界 Γ\GammaΓ(即 y=0,0<x<1y=0, 0<x<1y=0,0<x<1)上,由于有线积分项,自然边界条件通过变分得到。总变分中的边界积分和线积分项导致:
∫x=01η[2ux−4uy+2u+2g]dx=0
\int_{x=0}^{1} \eta \left[ 2 u_x - 4 u_y + 2 u + 2 g \right] dx = 0
∫x=01η[2ux−4uy+2u+2g]dx=0
由于 η\etaη 任意,必有:
2ux−4uy+2u+2g=0
2 u_x - 4 u_y + 2 u + 2 g = 0
2ux−4uy+2u+2g=0
除以2:
ux−2uy+u+g=0
u_x - 2 u_y + u + g = 0
ux−2uy+u+g=0
因此,在 Γ\GammaΓ 上,自然边界条件为:
ux−2uy+u+g=0ony=0,0<x<1
u_x - 2 u_y + u + g = 0 \quad \text{on} \quad y=0, 0<x<1
ux−2uy+u+g=0ony=0,0<x<1
总结
驻点满足的方程和边界条件如下:
- 方程:在区域 DDD 内部,
uxx−2uxy+2uyy=0 u_{xx} - 2 u_{xy} + 2 u_{yy} = 0 uxx−2uxy+2uyy=0 - 边界条件:
- u=0u = 0u=0 on y=1y=1y=1(上边界)
- u=0u = 0u=0 on x=0x=0x=0(左边界)
- u=0u = 0u=0 on x=2x=2x=2(右边界)
- u=0u = 0u=0 on y=0,1<x<2y=0, 1<x<2y=0,1<x<2(下边界部分)
- ux−2uy+u+g=0u_x - 2 u_y + u + g = 0ux−2uy+u+g=0 on y=0,0<x<1y=0, 0<x<1y=0,0<x<1(自然边界条件)
此即泛函 Φ(u)\Phi(u)Φ(u) 的驻点所满足的方程和边界条件。