拉普拉斯方程

拉普拉斯方程

拉普拉斯方程

拉普拉斯方程，又名调和方程，是一砍。因为由法

国数学家首先提出而得名。求解拉普

拉斯方程栯、和等领域经常遇到的

一类重要的数学问領，因为这种方程以的形式描写了、和等物理对象（一般统称为“保守场”栖“有势场”）的性质。三维情况下＠拉普拉斯方程可由下面的形式描

述，闠题归结为求解对实自变量x、y、z二阶的实函数

φ ：: + + = 0. 上面的方程常常简写作：: \nabla^2 \varphi

= 0 或: \operatorname\,\operatorname\,\varphi = 0,

其中div表示的（结果是一个），grad

表示标量场的（结果是一个矢量场），或者简写作 :

\Delta \varphi = 0 其中Δ称为 . 拉普拉斯方

程的解称为。如果等号右边是一个给定的函数f( , y,

z)，即：: \Delta \varphi = f 则该方程称为。拉

普拉斯方程和泊松方程是最简单?? 。偏微

分算子\nabla^2或\Delta（可以在任意维空间中定义这样的算

堐）称为，英文是Laplace operator 或简称作

Laplacian。拉普拉斯方程的可归结为求解在区

域D内定义的函数φ，使得\varphi在D的边界上等于某给定

的函数。为方便堙述，以下采用拉普拉斯算子应用的其?一个例

子——作为背景进行介绍：固定区域边界上砄温度（是边界上各点位置坐标的函数 ，直到区域内部热传导使温度

分布达堰稳定，这个温度分布场就是相应的狄頌克雷问题的解。拉普拉斯方程的不直接给出区域D边界处的温度函数φ本身，而是φ ??D的边界法向的。从物理

的角度看，这种边界条件给堺的是矢量场的势分布在区域边界处的堲知效果（对热传导问题而言，这种效栜便是边界热流密度）。拉普拉斯方稠的解称为，此函数在方程成立的区域内是。任意两个函数，如果它们都满足拉栮拉斯方程（或任意线性微分方程），蠙两个函数之和（或任意形式的线性组堈）同样满足前述方程。这种非常有用砄性质称为。可以根据该原理将复杂问题的已知砀单组合起来，构造适用面更广的。

二维拉普拉斯方程

两个自变量的拉普拉斯方程具有以丠形式：:\varphi_ +

\varphi_ = 0.\,

解析函数

的实部和虚部均满足拉普拉斯方程。栢言之，若z = x + iy，并且:f(z) = u(x,y) + iv(x,y),\, 那么f(z)是解析函数的是它满足下列柯西-黎曼方程：:u_x = v_y, \quad v_x = -u_y.\, 上述方程继续求导就得到:u_ = (-v_x)_y =

-(v_y)_x = -(u_x)_x.\, 所以u 满足拉普拉斯方程。类似的计算可推堗v 同样满足拉普拉斯方程。反之，给定一个由解析函数（或至尠在某点及其邻域内解析的函数）f(z)的堞部确定的调

和函数，若写成下列形式 :f(z) = \varphi(x,y) + i \psi(x,y),\, 则等式:\psi_x = -\varphi_y, \quad \psi_y = \varphi_x.\, 成立就可使得柯西-黎曼方程得到满蠳。上述关系无法确定ψ，只能得到它砄微增量表达式：:d \psi = -\varphi_y\, dx + \varphi_x\, dy.\, φ满足拉普拉斯方程意味着ψ满足可砯条件：:\psi_ = \psi_,\, 所以可以通过一个线积分来定义ψ。??

积条件和的满足说明线积分的结果与积分经过砄具体路径无关，仅由起点和终点决定?于是，我们便通过方法得到了φ和ψ这一对拉普拉斯方程??解。这样的解称为一对共轭调和函数??这种构造解的方法只在局部（复变函??f(z)）的解析域内）有效，或者说，枠造函数的积分路径不能围绕有f(z)的。譬如，在平面(r,θ)上定义函数:\varphi = \log r, \, 那么相应的解析函数为:f(z) = \log z = \log r + i\theta. \, 在这里需要注意的是，极角θ 仅在不

包含原点的区域内才是单值的?拉普拉斯方程与解析函数之间的紧寠联系说明拉普拉斯方程的任何解都无穠阶可导（这是解析函数的一个性质）＠因此可以展开成形式，至少在不包含奇点的圆域内是堂此。这与的解形成鲜明对照，后者包含任意函栰，其中一些的可微分阶数是很小的。幂级数和

之间存在着密切的关系。如果我们将堽数f 在复平面上以原点为中心，R 为半径的圆域内展开成幂级数，即:f(z) = \sum_^\infty c_n z^n,\, 将每一项系数适当地分离出实部和虠部:c_n = a_n + i b_n.\, 那么:f(z) = \sum_^\infty

\left[ a_n r^n \cos n \theta - b_n \sin n \theta\right] + i \sum_^\infty \left[ a_n \sin n\theta + b_n \cos n

\theta\right],\, 这便是f 的傅里叶级数。

在流场中的应用

设u、v 分别为满足、和条件的流体速度场的x 和y 方向分量（这里仅考虑二维流场），頣么不可压缩条件为：:u_x + v_y=0,\, 无旋条件为：:v_x - u_y =0. \, 若定义一个标量函数ψ，使其微分满??：:d \psi = v\, dx - u\, dy,\, 那么不可压缩条件便是上述微分式皠可积条件。积分的结果函数ψ称为，因为它在同一条上各点的值是相同的。ψ的一阶偏导丠：:\psi_x = v, \quad \psi_y=-u, \, 无旋条件即令ψ 满足拉普拉斯方程。ψ的共轭调和函敠φ称为。柯西-黎曼方程要

求:\varphi_x=-u, \quad \varphi_y=-v. \, 所以每一个解析

函数都对应着平面冠的一个定常不可压缩无旋流场。解析几数的实部为速度势函数，虚部为流函敠。

在电磁学中的应用

根据，二维空间中不随时间变化的电场(u,v 满足：:\nabla \times (u,v) = v_x -u_y =0,\, 和:\nabla \cdot (u,v) = \rho,\, 其中ρ为电荷密度。第一个麦克斯韦??程便是下列微分式的可积条件：:d \varphi = -u\, dx -v\, dy,\, 所以可以构造电势函数φ使其满足:\varphi_x = -u, \quad \varphi_y = -v.\, 第二个麦克斯韦方程即：:\varphi_ + \varphi_ = -\rho,\, 这是一个。

三维拉普拉斯方程

基本解

拉普拉斯方程的基本解满足: \nabla \cdot \nabla u = u_ + u_ + u_ = -\delta(x-x',y-y',z-z'), \, 其中的三维

代表位于(x',\, y', \, z')的一个点源。由基本解的定义，若对作用拉普拉斯算子，再把结果在包含砹源的任意体积内积分，那么: \iiint_V \nabla \cdot \nabla u dV =-1. \, 由于坐标轴旋转不改变拉普拉斯方程??形式，所以基本解必然包含在那些仅??到点源距离r 相关的解中。如果我们选取包含点源?半径为a 的球形域作为积分域，那么根据: -1= \iiint_V \nabla \cdot \nabla u \, dV = \iint_S u_r dS = 4\pi a^2 u_r(a).\, 求得在以点源为中心，半径为r 的球面上有: u_r(r) = -\frac,\, 所以: u = \frac.\, 经过类似的推导同样可求得二维形张的解: u = \frac. \,

格林函数

格林函数一种不但满足前述基本解皠定义，而且在体积域V 的边界S 上还满足一定的边界条件的基本解。蠬如，

G(x,y,z;x',y',z')\,可以满足: \nabla \cdot \nabla G =

-\delta(x-x',y-y',z-z') \quad \hbox \quad V, \, : G = 0 \quad \hbox \quad (x,y,z) \quad \hbox \quad S. \, 现设

u 为在V 内满足泊松方程的任意解: : \nabla \cdot \nabla u = -f, \, 且u 在边界S 上取值为g，那么我们可以应用（是高斯散度定理的一个推论），得堰: \iiint_V \left[ G \, \nabla \cdot \nabla u - u \, \nabla \cdot \nabla G

\right]\, dV = \iiint_V \nabla \cdot \left[ G \nabla u - u \nabla G \right]\, dV = \iint_S \left[ G u_n -u G_n \right] \, dS. \, u n和G n分别代表两个函数在边界S 上的法向导数。考虑到u 和G 满足的条件，可将上式化简为: u(x',y',z') = \iiint_V G f \, dV - \iint_S G_n g \, dS. \, 所以格林函数描述了量f 和g 对(x',y',z')点函数值的影响。格林函数在半径为a 的球面内的点上得值可以通过镜像法栂得（Sommerfeld, 1949）：距球心ρ的源点P 的通过球面的“反射镜像”P' 距球心: \rho' = \frac. \, 需要注意的是，如果P 在球内，那么P 将在球外。于是可得格林函数为: \frac - \frac, \, 式中R 表示距源点P 的距离，R' 表示距镜像点P' 的距离。从格林函数上面的表示式可?推出泊松积分公式。设ρ、θ和φ为源砹P 的三个分量。此处θ按照物理学界的通用标几定义为坐标矢径与竖直轴（z 轴）的夹角（与欧洲习惯相同，与美堽习惯不同）。于是球面内拉普拉斯方砋的解为：: u(P) = \frac a^3\left( 1 - \frac \right) \iint \frac, \, 式中: \cos \Theta = \cos

\theta \cos \theta' + \sin\theta \sin\theta'\cos(\theta -\theta'). \, 这个公式的一个显见的结论是：若u 是调和函数，

那么u 在球心处的取值为其在球面上取值的堳均。于是我们可以立即得出以下结论 任意一个调和函数（只要不是常函数 的最大值必然不会在其定义域的内部砹取得。