拓扑01 拉普拉斯方程 若u 是调和函数,那么u 在球心处的取值为其在球面上取值的平均。于是我们可以立即得出以下结论:任意一个调

来源: marketreflections 2011-06-26 07:51:47 [] [博客] [旧帖] [给我悄悄话] 本文已被阅读: 次 (47694 bytes)

拉普拉斯方程

维基百科,自由的百科全书
跳转到: 导航, 搜索
微积分学
e = \lim_{n\to\infty} \left(1+\frac{1}{n}\right)^n  
函数 · 导数 · 微分 · 积分

拉普拉斯方程(Laplace's equation),又名调和方程位势方程,是一种偏微分方程。因为由法国数学家拉普拉斯首先提出而得名。求解拉普拉斯方程是电磁学天文学流体力学等领域经常遇到的一类重要的数学问题,因为这种方程以势函数的形式描写了电场引力场流场等物理对象(一般统称为“保守场”或“有势场”)的性质。

目录

[隐藏]

[编辑] 定義

三维情况下,拉普拉斯方程可由下面的形式描述,问题归结为求解对实自变量xyz二阶可微的实函数φ :

 {\partial^2 \varphi\over \partial x^2 } + {\partial^2 \varphi\over \partial y^2 } + {\partial^2 \varphi\over \partial z^2 } = 0.

上面的方程常常简写作:

\nabla^2 \varphi = 0

\operatorname{div}\,\operatorname{grad}\,\varphi = 0,

其中div表示矢量场散度(结果是一个标量场),grad表示标量场的梯度(结果是一个矢量场),或者简写作:

\Delta \varphi = 0

其中Δ称为拉普拉斯算子.

拉普拉斯方程的解称为调和函数

如果等号右边是一个给定的函数f(x, y, z),即:

\Delta \varphi = f

则该方程称为泊松方程。 拉普拉斯方程和泊松方程是最简单的椭圆型偏微分方程。偏微分算子\nabla^2Δ(可以在任意维空间中定义这样的算子)称为拉普拉斯算子,英文是 Laplace operator 或简称作 Laplacian

拉普拉斯方程的狄利克雷问题可归结为求解在区域D内定义的函数φ,使得\varphiD的边界上等于某给定的函数。为方便叙述,以下采用拉普拉斯算子应用的其中一个例子——热传导问题作为背景进行介绍:固定区域边界上的温度(是边界上各点位置坐标的函数),直到区域内部热传导使温度分布达到稳定,这个温度分布场就是相应的狄利克雷问题的解。

拉普拉斯方程的诺伊曼边界条件不直接给出区域D边界处的温度函数φ本身,而是φ沿D的边界法向的导数。从物理的角度看,这种边界条件给出的是矢量场的势分布在区域边界处的已知效果(对热传导问题而言,这种效果便是边界热流密度)。

拉普拉斯方程的解称为调和函数,此函数在方程成立的区域内是解析的。任意两个函数,如果它们都满足拉普拉斯方程(或任意线性微分方程),这两个函数之和(或任意形式的线性组合)同样满足前述方程。这种非常有用的性质称为叠加原理。可以根据该原理将复杂问题的已知简单特解组合起来,构造适用面更广的通解

[编辑] 二维拉普拉斯方程

狄利克雷边界条件(u(r=2)=0、u(r=4)=4sin(5*θ))下的环形拉普拉斯方程(r=2、R=4)图形

两个自变量的拉普拉斯方程具有以下形式:

\varphi_{xx} + \varphi_{yy} = 0.\,

[编辑] 解析函数

解析函数的实部和虚部均满足拉普拉斯方程。换言之,若z = x + iy,并且

f(z) = u(x,y) + iv(x,y),\,

那么f(z)是解析函数的充要条件是u(x,y),v(x,y)可微,且满足下列柯西-黎曼方程

u_x = v_y, \quad v_x = -u_y.\,

上述方程继续求导就得到

u_{yy} = (-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_{xy} = \psi_{yx},\,

所以可以通过一个线积分来定义ψ。可积条件和斯托克斯定理的满足说明线积分的结果与积分经过的具体路径无关,仅由起点和终点决定。于是,我们便通过复变函数方法得到了φ和ψ这一对拉普拉斯方程的解。这样的解称为一对共轭调和函数。这种构造解的方法只在局部(复变函数f(z))的解析域内)有效,或者说,构造函数的积分路径不能围绕有f(z)的奇点。譬如,在极坐标平面(r,θ)上定义函数

\varphi = \log r, \,

那么相应的解析函数为

f(z) = \log z = \log r + i\theta. \,

在这里需要注意的是,极角θ 仅在不包含原点的区域内才是单值的。

拉普拉斯方程与解析函数之间的紧密联系说明拉普拉斯方程的任何解都无穷阶可导(这是解析函数的一个性质),因此可以展开成幂级数形式,至少在不包含奇点的圆域内是如此。这与波动方程的解形成鲜明对照,后者包含任意函数,其中一些的可微分阶数是很小的。

幂级数和傅里叶级数之间存在着密切的关系。如果我们将函数f 在复平面上以原点为中心,R 为半径的圆域内展开成幂级数,即

f(z) = \sum_{n=0}^\infty c_n z^n,\,

将每一项系数适当地分离出实部和虚部

c_n = a_n + i b_n.\,

那么

f(z) = \sum_{n=0}^\infty \left[ a_n r^n \cos n \theta - b_n \sin n \theta\right] + i \sum_{n=1}^\infty \left[ a_n \sin n\theta + b_n \cos n \theta\right],\,

这便是f 的傅里叶级数。

[编辑] 在流场中的应用

uv 分别为满足定常不可压缩无旋条件的流体速度场的xy 方向分量(这里仅考虑二维流场),那么不可压缩条件为:

u_x + v_y=0,\,

无旋条件为:

v_x - u_y =0. \,

若定义一个标量函数ψ,使其微分满足:

d \psi = -v\, dx + u\, dy,\,

那么不可压缩条件便是上述微分式的可积条件。积分的结果函数ψ称为流函数,因为它在同一条流线上各点的值是相同的。ψ的一阶偏导为:

\psi_x = -v, \quad \psi_y=u, \,

无旋条件即令 ψ 满足拉普拉斯方程。ψ的共轭调和函数\varphi称为速度势。 柯西-黎曼方程要求

\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_{xx} + \varphi_{yy} = -\rho,\,

这是一个泊松方程

[编辑] 三维拉普拉斯方程

[编辑] 基本解

拉普拉斯方程的基本解满足

 \nabla \cdot \nabla u = u_{xx} + u_{yy} + u_{zz} = -\delta(x-x',y-y',z-z'), \,

其中的三维δ函数代表位于 (x',\, y', \, z')的一个点源。 由基本解的定义,若对u 作用拉普拉斯算子,再把结果在包含点源的任意体积内积分,那么

 \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{1}{4\pi r^2},\,

所以

 u = \frac{1}{4\pi r}.\,

经过类似的推导同样可求得二维形式的解

 u = \frac{-\log r}{2\pi}. \,

[编辑] 格林函数

格林函数是一种不但满足前述基本解的定义,而且在体积域V 的边界S 上还满足一定的边界条件的基本解。譬如, G(x,y,z;x',y',z')\,可以满足

 \nabla \cdot \nabla G = -\delta(x-x',y-y',z-z') \quad \hbox{in} \quad V, \,
 G = 0 \quad \hbox{if} \quad (x,y,z) \quad \hbox{on} \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. \,

unGn 分别代表两个函数在边界S 上的法向导数。考虑到uG 满足的条件,可将上式化简为

 u(x',y',z') =  \iiint_V G f \, dV - \iint_S G_n g \, dS. \,

所以格林函数描述了量fg(x',y',z')点函数值的影响。格林函数在半径为a 的球面内的点上得值可以通过镜像法求得(Sommerfeld, 1949):距球心ρ的源点P 的通过球面的“反射镜像”P' 距球心

 \rho' = \frac{a^2}{\rho}. \,

需要注意的是,如果P 在球内,那么P' 将在球外。于是可得格林函数为

 \frac{1}{4 \pi R} - \frac{a}{4 \pi \rho R'}, \,

式中R 表示距源点P 的距离,R' 表示距镜像点P' 的距离。从格林函数上面的表示式可以推出泊松积分公式。设ρ、θ和φ为源点P 的三个球坐标分量。此处θ按照物理学界的通用标准定义为坐标矢径与竖直轴(z 轴)的夹角(与欧洲习惯相同,与美国习惯不同)。于是球面内拉普拉斯方程的解为:

 u(P) = \frac{1}{4\pi} a^3\left( 1 - \frac{\rho^2}{a^2} \right) \iint \frac{g(\theta',\varphi') \sin \theta' \, d\theta' \, d\varphi'}{(a^2 + \rho^2 - 2 a \rho \cos \Theta)^{3/2} },  \,

式中

 \cos \Theta = \cos \theta \cos \theta' + \sin\theta \sin\theta'\cos(\theta -\theta'). \,

这个公式的一个显见的结论是:若u 是调和函数,那么u 在球心处的取值为其在球面上取值的平均。于是我们可以立即得出以下结论:任意一个调和函数(只要不是常函数)的最大值必然不会在其定义域的内部点取得。

[编辑] 参见

[编辑] 外部链接

[编辑] 参考文献

  • 严镇军编,《数学物理方程》,第二版,中国科学技术大学出版社,合肥,2002,ISBN 7-312-00799-6/O·177
  • L.C. Evans, Partial Differential Equations, American Mathematical Society, Providence, 1998. ISBN 0-8218-0772-2
  • I. G. Petrovsky, Partial Differential Equations, W. B. Saunders Co., Philadelphia, 1967.
  • A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC Press, Boca Raton, 2002. ISBN 1-58488-299-9
  • A. Sommerfeld, Partial Differential Equations in Physics, Academic Press, New York, 1949.
  • Pijush K.Kundu, Fluid Mechanics, Academic Press, 2002.
请您先登陆,再发跟帖!

发现Adblock插件

如要继续浏览
请支持本站 请务必在本站关闭Adblock

关闭Adblock后 请点击

请参考如何关闭Adblock

安装Adblock plus用户请点击浏览器图标
选择“Disable on www.wenxuecity.com”

安装Adblock用户请点击图标
选择“don't run on pages on this domain”