作者Emcc (mass, momentum, energy)
看板Physics
标题Re: [问题] 请教 Green's Function 用在 PDE 的概念?
时间Thu Jan 28 22:32:34 2010
: 解 PDE 除了 eigenfunction 的方法, 还可以用 Green's Function
: 2
: ex. ▽ u = Q st. u = f (r) , δu/δn = g (r)
^^^^^^^^^
Poisson's eqn.
这个方程式 除了你提的inhomogeneous 这项特性之外
还有一个重点 就是它是"线性"的 (
linear)
Poisson's eqn. 代表的物理意义
就是一个有input/output的扩散系统
inhomogeneous term, Q(r) 就是系统外来的input
所要求解的u(r) 就是系统对这个input Q(r) 所产生的"响应"(response)
就我熟悉的热传学来举例
Q(r)就是一个热源, heat source
u(r)就是物质因为热源 所形成的温度场T
Laplacian operator ▽^2
在物理上代表一个很重要的特性
就是"扩散"
这家伙很有意思
它喜欢"齐头式的平等"
也就是 你高 它把你压下来
你矮 它把你拉上去
至於为什麽 温度场的扩散 ▽^2[u(r)] 会等於热源Q(r)
中间那个"等号" 为何能成立?
那是因为物理性质的 "守恒定律"
以热传的例子来说
这个方程式的意义为 "能量的守恒"
同样的Poisson's eqn. 在电磁学里
Q(r)为电荷密度
u(r)为 electric potential field
方程式代表 "电位的扩散"
(电磁的部分小弟没有很熟 还请指正)
现在回到PDE的 Green's function解法
如果今天对一物理系统
施以一个"脉冲"(impulse input) 在一点r0
即 δ(r-r0)
则满足Poisson's eqn.的系统对这个 impulse所产生的响应
就是它的Green's function
"系统对某一点 r0 的input 在另一点 r 所造成的响应" 为G(r,r0) r≠r0
满足governing eqn.
=> ▽^2[G(r,r0)] = δ(r-r0)
重点来了
如果impulse input不只发生在一点
而是数个点
那麽系统会如何响应?
"因为Poisson's eqn.是线性的"
线性的好处是可以叠加(superposition)
所以对所有impulse input的点作叠加 ∫ [Q(r0).δ(r-r0)] d(r0) = Q(r)
其响应为 ∫ [Q(r0).G(r,r0)] dr0 = u(r)
当然
Poisson's eqn. 只是最基本的Green's function方法求解DE的例子
其他还有更复杂的PDE或是ODE
但是Green's function的基本精神就是这样
"线性 所以可叠加"
另外有BC存在
系统的响应也会受BC影响
那就是你写的後面那项了
: 解法就是视 nonhomogeneous term Q 与 boundary conditions 为许多离散点
: 而这些离散点影响着 u
: 离散点( Free Space Green's Function )定义为:
: 2
: ▽ G ( r,r0 ) = delta ( r-r0 )
: 使用 Green's second identity 可以导出
: u (r) =
: δG δu
: ∫ G(r,r0).Q(r0).dV(r0) + ∫ [ u(r1).——(r,r1) - G(r,r1)——.(r1) ] dS(r1)
: V S δn δn
: source boundary conditions
: 而且可求出 G(r) = -1/4π(r-r0), for r =/= r0
: ----------------------------------------------------------------------------
: 问题来了, 我们可以观察出 u 受到 Q 和 BCs 两者的影响
: 可是:
: 2
: 1) 为何 ▽ G 会 model 成 delta function ,而非 ▽G 或 G
: 有实际的物理系统可以辅助理解吗?
: 可以解释 G 和 Q 交互作用, u 和 G'作用, G 和 u'负作用, 三者合起来影响 u 吗?
: 2) G 如果是正比 (1/r), 用一维的例子来看, G' 长怎样呢?
: G''又会是delta function 吗? 原点外的其他点的二阶微分不是零啊?
: 可以用一维空间举个三者叠加成 u 的例子吗?
--
The World Series champions of 2006 is
St.Louis Cardinals !!!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.247.229
1F:→ kuromu:好详细 但是我也想问一下 边界项代表的意义..... 01/28 23:19
2F:推 cobeycobey:推.好详细 01/28 23:56
3F:推 WINDHEAD:恩,只有当▽^2[u(r)]=0才是"齐头式的平等喔" 01/29 08:50
4F:推 qna:边界的贡献是说 边界外面也可能有电荷 他也会贡献内部电位阿 01/29 10:39
5F:推 LeeSeDol:推 01/29 15:10
6F:→ kuromu:但是如何从公式上解释边界外电荷如何影响内部电位? 01/29 17:39
7F:→ kuromu:请问一下 为什麽 01/29 19:28
8F:→ kuromu: δG δu 01/29 19:28
9F:→ kuromu:∫ [ u.—— - G.—— ]dS 中的 01/29 19:28
10F:→ kuromu: S δn δn 01/29 19:28
11F:→ kuromu:u可解释为surface dipole density 01/29 19:28
12F:→ kuromu: δu 01/29 19:28
13F:→ kuromu:而 ——可解释为surface charge density ? 01/29 19:28
14F:→ kuromu: δn 01/29 19:28
15F:推 nevinyrrals:专业推 01/30 01:57
16F:→ nevinyrrals:其实电磁学中的电荷守衡应该是divJ+dρ/dt=0 01/30 02:04
17F:→ nevinyrrals:Poisson's eqn.指的是在有电荷源之下电位的扩散情形 01/30 02:06
18F:→ nevinyrrals:不过我想其基本精神是相同的 01/30 02:07
19F:推 LeeSeDol:原po指的应该是"无源场中一封闭曲面的净通量恒为零" 01/30 02:11
20F:→ Emcc:感谢nevinyrrals指正 01/30 02:30
※ 编辑: Emcc 来自: 140.112.247.229 (01/30 02:34)