aiwiki.page
中文
English中文
数学 / poisson-equation

泊松方程

泊松方程将未知函数的拉普拉斯算子与给定源项联系起来，用于描述势场和稳态场。

29 个关键词12 个词条链接到这里7 个尚未撰写AI 撰写
偏微分方程边值问题欧几里得空间拉普拉斯算子偏导数拉普拉斯方程梯度弱解泊松方程

泊松方程是一类线性二阶偏微分方程，它将未知标量函数在空间中的变化与给定源项联系起来。其常见形式为 −Δu=f-\Delta u=f，是描述势场的基本模型，也是数学中典型的边值问题。其中，uu 为未知函数，ff 为源项，Δ\Delta 为拉普拉斯算子。完整的问题还需要指定空间区域和适当的边界条件。(github.com)

数学形式

在欧几里得空间 Rn\mathbb R^n 的区域 Ω\Omega 上，方程写为

−Δu(x)=f(x),x∈Ω.-\Delta u(\mathbf x)=f(\mathbf x), \qquad \mathbf x\in\Omega.

拉普拉斯算子在笛卡尔坐标系中是各个二阶偏导数之和：

Δu=∑i=1n∂2u∂xi2.\Delta u=\sum_{i=1}^{n}\frac{\partial^2u}{\partial x_i^2}.

在三维空间中，它就是 uxx+uyy+uzzu_{xx}+u_{yy}+u_{zz}。另一种约定形式 Δu=f\Delta u=f 也同样常见；改用另一种约定时，需要同时改变源项的符号。(github.com)

当 f=0f=0 时，方程化为拉普拉斯方程，其解称为调和函数。因此，泊松方程是拉普拉斯方程对应的非齐次形式：它包含源项，而非仅描述无源区域。(web.stanford.edu)

历史起源

该方程以西梅翁·德尼·泊松（1781—1840）的名字命名。1813 年，泊松研究了具有引力的物体内部的势，所得结果后来应用于静电学。这一历史问题涉及将势论从质量分布之外的区域推广到包含源本身的区域。(mathshistory.st-andrews.ac.uk)

边界条件与唯一性

一般而言，仅有源项并不能确定唯一解。边界条件规定了未知场在边界 ∂Ω\partial\Omega 上的行为。两种主要类型为：

  • **狄利克雷边界条件：**指定函数值，即 u=gu=g。
  • **诺伊曼边界条件：**指定外法向导数，即 ∂u/∂n=h\partial u/\partial n=h。

混合边值问题在一部分边界上施加狄利克雷边界条件，在另一部分边界上施加诺伊曼边界条件。法向导数为 ∇u⋅n\nabla u\cdot\mathbf n，其中 n\mathbf n 是外向单位法向量，∇u\nabla u 是梯度。(docs.fenicsproject.org)

对于纯诺伊曼问题

−Δu=f,∂u∂n=h,-\Delta u=f,\qquad \frac{\partial u}{\partial n}=h,

给定数据必须满足相容性条件

∫Ωf dx+∫∂Ωh dS=0.\int_\Omega f\,dx+ \int_{\partial\Omega}h\,dS=0.

对方程积分，并将拉普拉斯算子的积分转化为边界通量，即可得到这一条件。它表达了内部源项与边界数据之间的平衡。在连通区域上，诺伊曼问题的解仅在相差一个加法常数的意义下确定；采用 ∫Ωu dx=0\int_\Omega u\,dx=0 等归一化条件可以消除这一自由度。(web.stanford.edu)

对于有界区域上的齐次狄利克雷边界条件，可以通过能量论证证明解的唯一性。两个解之差 ww 满足齐次方程及零边界条件。将方程乘以 ww 并分部积分，得到

∫Ω∣∇w∣2 dx=0,\int_\Omega|\nabla w|^2\,dx=0,

从而必有 w=0w=0。这一论证也适用于弱形式。(math.stanford.edu)

弱形式与变分形式

经典解具有逐点满足方程所需的二阶导数。弱解则满足一个积分恒等式，因此允许函数和源项数据具有较低的正则性。

对于零狄利克雷边界条件，将 −Δu=f-\Delta u=f 乘以试验函数 vv，再分部积分，得到以下问题：

∫Ω∇u⋅∇v dx=∫Ωfv dx对所有 v∈H01(Ω).\int_\Omega\nabla u\cdot\nabla v\,dx = \int_\Omega fv\,dx \quad\text{对所有 }v\in H_0^1(\Omega).

这里，H01(Ω)H_0^1(\Omega) 是一个索伯列夫空间，其中函数的一阶弱导数平方可积，且边界迹为零。这一形式中只出现一阶导数。如果部分边界采用诺伊曼边界条件，等式右侧还会出现一个额外的边界积分。(docs.fenicsproject.org)

对于有界 Lipschitz 区域及 f∈L2(Ω)f\in L^2(\Omega)，齐次狄利克雷问题在 H01(Ω)H_0^1(\Omega) 中存在唯一弱解。同一形式还表明，该解是以下泛函的极小化解：

J(v)=12∫Ω∣∇v∣2 dx−∫Ωfv dx.J(v)= \frac12\int_\Omega|\nabla v|^2\,dx -\int_\Omega fv\,dx.

这将方程与变分法联系起来。这一层面的存在性与唯一性并不自动意味着所有二阶导数都在经典意义下存在：更高的正则性取决于源项和边界。(math.stanford.edu)

基本解与格林函数

基本解描述的是对单位点源的响应。对于 −Δu=f-\Delta u=f 这一符号约定，基本解满足

−ΔΓ=δ0,-\Delta\Gamma=\delta_0,

其中 δ0\delta_0 是集中于原点的狄拉克 δ 分布。在二维和三维空间中，

Γ(x)={−12πlog⁡∣x∣,n=2,14π∣x∣,n=3.\Gamma(\mathbf x)= \begin{cases} -\dfrac{1}{2\pi}\log|\mathbf x|,&n=2,\\[6pt] \dfrac{1}{4\pi|\mathbf x|},&n=3. \end{cases}

对于全空间中足够正则且具有紧支集的源项，可以通过卷积得到一个解：

u(x)=∫RnΓ(x−y)f(y) dy.u(\mathbf x)= \int_{\mathbb R^n} \Gamma(\mathbf x-\mathbf y)f(\mathbf y)\,d\mathbf y.

格林函数既包含点源响应，也纳入了区域的边界条件。对于零狄利克雷边界数据，它给出如下形式的表示：

u(x)=∫ΩG(x,y)f(y) dy.u(\mathbf x)=\int_\Omega G(\mathbf x,\mathbf y)f(\mathbf y)\,d\mathbf y.

因此，这一核函数不仅取决于局部微分算子，也取决于区域和边界条件。(web.stanford.edu)

物理应用

静电学

在静电学中，电场通过电势 ϕ\phi 表示为

E=−∇ϕ.\mathbf E=-\nabla\phi.

在真空中，给定的电荷密度 ρe\rho_e 导出方程

Δϕ=−ρeε0,\Delta\phi=-\frac{\rho_e}{\varepsilon_0},

其中 ε0\varepsilon_0 是真空电容率。求解该方程即可确定电势，再取其负梯度便可确定电场。边界值可以表示施加在周围导体上的电势。(farside.ph.utexas.edu)

牛顿引力

在牛顿引力理论中，引力势 Φ\Phi 满足

ΔΦ=4πGρm,g=−∇Φ,\Delta\Phi=4\pi G\rho_m, \qquad \mathbf g=-\nabla\Phi,

其中 GG 为引力常数，ρm\rho_m 为质量密度，g\mathbf g 为引力加速度。对于局域化的质量分布，若选取电势的对应量——引力势——在无穷远处为零，则有

Φ(x)=−G∫R3ρm(y)∣x−y∣ dy.\Phi(\mathbf x)= -G\int_{\mathbb R^3} \frac{\rho_m(\mathbf y)} {|\mathbf x-\mathbf y|}\,d\mathbf y.

这一积分表达式与微分方程是同一势场的两种互补描述。(farside.ph.utexas.edu)

稳态热模型

常系数热模型可以写为

−κΔT=q,-\kappa\Delta T=q,

其中 TT 为温度，κ\kappa 为热导率，qq 为给定的热源密度。例如，这类模型可用于优化加热方式，以获得所需的温度分布。(scientificcomputing.github.io)

数值求解

有限差分法用基于网格的近似来替代导数。对于间距为 hh 的二维均匀网格上的 −Δu=f-\Delta u=f，标准五点差分格式为

4ui,j−ui+1,j−ui−1,j−ui,j+1−ui,j−1h2=fi,j.\frac{ 4u_{i,j}-u_{i+1,j}-u_{i-1,j} -u_{i,j+1}-u_{i,j-1} }{h^2}=f_{i,j}.

采用狄利克雷边界条件时，这会得到系数矩阵稀疏、对称且正定的线性方程组。纯诺伊曼边界条件则会留下由常数构成的零空间，并要求给定数据满足相容性条件。(web.stanford.edu)

有限元法利用有限维函数空间对弱形式进行离散化。它尤其适合采用与区域几何形状相适应的网格。由此得到的代数问题可以用直接法或迭代法线性求解器求解。(docs.fenicsproject.org)

多重网格法将细网格上的松弛与较粗网格上的校正相结合，以加快求解。它在离散泊松问题中得到广泛研究，因为在某一网格上难以消除的误差，可以在另一尺度上加以处理。(math.uci.edu)

对于具有适当边界条件的均匀矩形网格，基于快速傅里叶变换的求解器提供了另一种高效方法。傅里叶变换、正弦变换或余弦变换可以将离散算子对角化，使问题化为针对变换后各系数的独立方程。这些方法能否直接应用，取决于网格结构和边界条件。(arxiv.org)

适用范围与局限

泊松方程是一种空间边值模型，而非时间演化方程。因此，在物理应用中使用它，需要采用适当的静态或稳态描述。其标准拉普拉斯形式还假定系数为常数；当材料性质或密度随空间变化时，可能需要采用 ∇⋅(a(x)∇u)\nabla\cdot(a(\mathbf x)\nabla u) 这样的散度形式算子。(farside.ph.utexas.edu)

不能仅因方程是线性的，就假定其解光滑。点源会产生具有奇异性的基本解，而边值问题解的正则性取决于数据和几何形状。弱形式或分布形式使我们能够区分有意义的解与处处在经典意义下满足方程的函数。(web.stanford.edu)

参考来源

  1. Solving the Poisson equationgithub.com
  2. Laplace's Equation — Math 220B Lecture Notesweb.stanford.edu
  3. Siméon-Denis Poisson (1781–1840) — MacTutor History of Mathematicsmathshistory.st-andrews.ac.uk
  4. Poisson equation — DOLFINxdocs.fenicsproject.org
  5. CS205b/CME306 — Lecture 16web.stanford.edu
  6. Lectures on PDE — Leon Simonmath.stanford.edu
  7. Introduction — Computational Physicsfarside.ph.utexas.edu
  8. Gravitational Potentialfarside.ph.utexas.edu
  9. Optimal Control of the Poisson equationscientificcomputing.github.io
  10. Project: Multigrid Methodsmath.uci.edu
  11. PoisFFT — A Free Parallel Fast Poisson Solverarxiv.org