1. 前置知识
足够光滑的向量场 u , v : Ω → R N u,v:\Omega\to \mathbb{R}^N u,v:Ω→RN ,其中
∇ v = ( ∂ i v ) i = 1 N , ∣ ∇ v ∣ : = ( ∑ i = 1 N ∣ ∂ i v ∣ 2 ) 1 2 , ∇ u ⋅ ∇ v = ∑ i = 1 N ∂ i u ∂ i v \nabla v = (\partial_i v)_{i=1}^N,|\nabla v| : = (\mathop{\sum}\limits_{i=1}^N |\partial_i v|^2)^{\frac{1}{2}}, \nabla u \cdot \nabla v = \mathop{\sum}\limits_{i=1}^N\partial_iu\partial_iv ∇v=(∂iv)i=1N,∣∇v∣:=(i=1∑N∣∂iv∣2)21,∇u⋅∇v=i=1∑N∂iu∂iv
其中 ∣ ⋅ ∣ , ⋅ |\cdot| , \cdot ∣⋅∣,⋅ 分别表示2范数和欧式内积
Ω \Omega Ω 是 R N \mathbb{R}^N RN 中的区域, ν = ( ν i ) i = 1 N \nu = (\nu_i)_{i=1}^N ν=(νi)i=1N 表示沿着 Γ = ∂ Ω \Gamma = \partial \Omega Γ=∂Ω (就是边界)的单位外法向量场。(很像梯度)
迹: v v v 是 Ω ‾ → R \overline \Omega \to \mathbb{R} Ω→R 上的连续函数, Ω \Omega Ω 是 R \mathbb{R} R 的开子集,则在边界 Γ : = ∂ Ω \Gamma :=\partial \Omega Γ:=∂Ω 上的迹为连续函数 t r v : Γ → R tr~v:\Gamma\to R tr v:Γ→R ,有 ( t r v ) ( x ) = v ( x ) , x ∈ Γ (tr~v)(x) = v(x),x\in \Gamma (tr v)(x)=v(x),x∈Γ
所以讲一个函数的迹就是限制在边界上的映射。
1.1. 格林公式
- 对任意 u ∈ H 2 ( Ω ) u\in H^2(\Omega) u∈H2(Ω) ,令
Δ u = ∑ i = 1 N ∂ i i u ∈ L 2 ( Ω ) \Delta u = \mathop{\sum}\limits_{i=1}^N \partial_{ii}u\in L^2(\Omega) Δu=i=1∑N∂iiu∈L2(Ω)
(看上去这好像就是个散度?)
∂ ν u = ∑ i = 1 N ν i ∂ i u ∈ L 2 ( Γ ) \partial_\nu u = \mathop{\sum}\limits_{i=1}^N\nu_i\partial_i u\in L^2(\Gamma) ∂νu=i=1∑Nνi∂iu∈L2(Γ)
(用单位法向量和函数的迹做欧式内积吗,得到迹在法向量方向的投影)
其中 ∂ i u ∈ L 2 ( Γ ) \partial_i u\in L^2(\Gamma) ∂iu∈L2(Γ) 表示函数 ∂ i u ∈ H 1 ( Ω ) \partial_i u\in H^1(\Omega) ∂iu∈H1(Ω) 在 Γ \Gamma Γ 上的迹。
此时下Green公式成立:
∫ Ω ∇ u ⋅ ∇ v d x = − ∫ Ω ( Δ u ) v d x + ∫ Γ ( ∂ ν u ) v d Γ \int_\Omega \nabla u \cdot \nabla v dx = - \int_\Omega (\Delta u)vdx + \int_\Gamma (\partial_\nu u) vd\Gamma ∫Ω∇u⋅∇vdx=−∫Ω(Δu)vdx+∫Γ(∂νu)vdΓ
- 给定函数 a ∈ C 2 ( Ω ‾ ) a\in \mathcal{C}^2(\overline{\Omega}) a∈C2(Ω) 以及 u ∈ H 1 ( Ω ) u\in H^1(\Omega) u∈H1(Ω) ,则函数 a u au au 属于空间 H 1 ( Ω ) H^1(\Omega) H1(Ω) ,且下格林公式成立
∫ Ω a u ∂ j v d x = − ∫ Ω ( ∂ j ( a u ) ) v d x + ∫ Γ a u v ν j d Γ \int_\Omega au\partial_j vdx = -\int_\Omega(\partial_j(au))vdx + \int_\Gamma auv\nu_j d\Gamma ∫Ωau∂jvdx=−∫Ω(∂j(au))vdx+∫ΓauvνjdΓ
1.2. 算子
Δ = ∑ i = 1 N ∂ i i \Delta = \mathop{\sum}\limits_{i=1}^N \partial_{ii} Δ=i=1∑N∂ii 作用在定义于 Ω \Omega Ω 内的函数上,称为Laplace算子
Δ u \Delta u Δu 称为 u 的 Laplace 算符。
∂ ν = ∑ i = 1 N ν i ∂ i \partial_\nu = \mathop{\sum}\limits_{i=1}^N\nu_i \partial_i ∂ν=i=1∑Nνi∂i 作用于定义在边界 Γ \Gamma Γ 上的函数上,称为外法向微分算子, ∂ ν u \partial_\nu u ∂νu 称为 u u u 的外法向导数。
2. 边值问题
给定
f
∈
L
2
(
Ω
)
f\in L^2(\Omega)
f∈L2(Ω) ,令
V
=
H
0
1
(
Ω
)
V = H_0^1(\Omega)
V=H01(Ω)
a ( u , v ) = ∫ Ω ( ∇ u ⋅ ∇ v + c u v ) d x , u , v ∈ V a(u,v) = \int_\Omega(\nabla u\cdot\nabla v + cuv)dx, u,v\in V a(u,v)=∫Ω(∇u⋅∇v+cuv)dx,u,v∈V
l ( v ) = ∫ Ω f v d x , v ∈ V l(v)=\int_\Omega fvdx,v\in V l(v)=∫Ωfvdx,v∈V
则有唯一的函数 u ∈ H 0 1 ( Ω ) u\in H_0^1(\Omega) u∈H01(Ω) 使以下泛函极小化
J ( v ) : = 1 2 a ( v , v ) − l ( v ) = 1 2 ∫ Ω ( ∣ ∇ v ∣ 2 + c v 2 ) d x − ∫ Ω f v d x , v ∈ H 0 1 ( Ω ) J(v):=\frac{1}{2}a(v,v)-l(v) =\frac{1}{2}\int_\Omega (|\nabla v|^2 + cv^2)dx - \int_\Omega fvdx,~ v\in H_0^1(\Omega) J(v):=21a(v,v)−l(v)=21∫Ω(∣∇v∣2+cv2)dx−∫Ωfvdx, v∈H01(Ω)
即满足变分方程 a ( u , v ) = l ( v ) , v ∈ H 0 1 ( Ω ) a(u,v) = l(v),v\in H_0^1(\Omega) a(u,v)=l(v),v∈H01(Ω)
一个 f f f 对应一个 u u u 嘛,所以可以定义线性映射来表示这种关系
f ∈ L 2 ( Ω ) → u ∈ H 0 1 ( Ω ) f\in L^2(\Omega) \to u \in H_0^1(\Omega) f∈L2(Ω)→u∈H01(Ω)
这个映射是连续的,而且函数 u u u 满足边值问题:
− Δ u + c u = f -\Delta u +cu = f −Δu+cu=f 在 Ω \Omega Ω 上,且 u = 0 u=0 u=0 在 Γ \Gamma Γ 上
别的就不证明了,最后这个边值问题来自于变分方程。
∫ Ω ( ∇ u ⋅ ∇ v + c u v ) d x = ∫ Ω f v d x \int_\Omega(\nabla u\cdot\nabla v + cuv)dx = \int_\Omega fvdx ∫Ω(∇u⋅∇v+cuv)dx=∫Ωfvdx
即为
∫ Ω [ ( ∇ u ⋅ ∇ v + c u v ) − f v ] d x = 0 \int_\Omega[(\nabla u\cdot\nabla v + cuv) - fv]dx = 0 ∫Ω[(∇u⋅∇v+cuv)−fv]dx=0
可得
− Δ u + c u − f = 0 -\Delta u + cu -f =0 −Δu+cu−f=0
(边界条件 u=0,第一项用格林公式替换)
是其次的
更一般的有
− Δ u + c u = f 在 Ω 中 , u = u 0 在 Γ 上 -\Delta u + cu = f 在 \Omega 中, u=u_0 在 \Gamma 上 −Δu+cu=f在Ω中,u=u0在Γ上
称为偏微分算子
L
:
v
→
L
v
=
−
Δ
v
+
c
v
\mathcal{L}:v\to \mathcal{L}v = -\Delta v + cv
L:v→Lv=−Δv+cv 的 Dirichlet问题。
(这是非齐次的)
特殊情况:
Possion方程: − Δ u = f -\Delta u = f −Δu=f
Laplace方程: − Δ u = 0 -\Delta u = 0 −Δu=0
(Laplace 方程的解是调和函数)