使用定点迭代解耦PDE系统


12

假设我有一个边值问题:

d2udx2+dvdx=f in Ω
dudx+d2vdx2=g in Ω
u=h in Ω

我的目标是将这个耦合问题的解决方案分解为一系列非耦合PDE。为了使系统解耦,我对一系列近似值应用定点迭代(uk,vk)以便

d2ukdx2+dvk1dx=f
duk1dx+d2vkdx2=g

从理论上讲,这将使我能够将两个方程作为纯椭圆PDE求解。但是,我从未见过以这种方式将定点迭代应用于PDE。我已经看到定点迭代应用于数字离散方程(有限差分法,有限元方法等),但从未直接应用于连续方程。

我这样做违反了任何公然的数学原理吗?这在数学上有效吗?我可以通过使用应用于连续变量问题而不是离散变量问题的定点迭代来将耦合PDE求解为非耦合PDE的序列吗?

在这一点上,我并不真正在乎使用这种方法是否可行,而是在理论上是否可行。任何反馈将不胜感激!


6
在双曲线PDE文献中,分数步和运算符拆分方法可以完成您在上面描述的操作。
Geoff Oxberry

您是说而不是 u kp k吗?(uk,vk)(uk,pk)
比尔·巴特

@BillBarth:是的!我只是更正了。
保罗

@GeoffOxberry:我发现运算符拆分的字符非常不同。
2014年

@Paul:我可以想到至少一个其他问题,其中“耦合PDE”通过定点迭代来解决(而不仅仅是公式化为定点问题):域分解,请参见例如Neumann–Dirichlet方法。(这里的区别是您有两个PDE,但它们位于不同的域上,并且仅通过接口进行耦合)。
2014年

Answers:


11

C(Ω)×C(Ω)

d2ukdx2+dvk1dx=fd2vkdx2+duk1dx=g
(加上边界条件)。

显然,如果此序列收敛,它将是您原始PDE集的解决方案。

xkxk+1u0v0

(ukvk)(u^kv^k)q(uk1vk1)(u^k1v^k1)
|q|<1(uk1,vk1)(u^k1,v^k1)

此逻辑在连续空间和离散空间中均有效。


3
不应该吗?|q|<1
保罗
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.