自然边界归化方法在求解抛物方程初边值外区域问题中的应用.docxVIP

自然边界归化方法在求解抛物方程初边值外区域问题中的应用 0 自然边界二元耦合算法 冯康和余德浩率先并开发了这种自然边界元方法。与经典的边界分布方法相比，它具有独特的优势：不仅易于实现，而且数值稳定性好，适合将开元方法基于同一变分原理，并与开元方法进行自然结合。目前对椭圆型问题而言,自然边界元方法已形成相当完整的理论,并已成功地应用于区域分解。但是,这一有效而重要的方法一直未能在发展型问题(如热传导方程,波动方程等)中得到实现。近期作者寻找到一种有效方法,已成功地实现了抛物型方程的自然边界归化,以此为基础可作一些相关问题的研究。 本文以二维热传导方程第一初边值问题为例提出一种新的耦合算法。我们首先给出圆外区域上的自然边界归化,然后以此为基础,研究一种自然边界元与有限元的耦合算法。由于人工边界的选取,将原无界区域化为一个有界子区域及一个规则的无界区域(圆外区域)。有界区域上的子问题用有限元法去求解,而无界区域上的问题则用自然边界元法去处理。由于用自然边界元法处理圆外区域,所以计算量大为减少。最后我们给出数值例子,数值结果表明此法是有效的。 1 时间步长的离散化 设Ω为R2中具有适当光滑边界Γ:=?Ω的单连通有界区域,Ωc:=R2\ˉΩ.对任意固定的正实数T,记J:=(0,T].我们考虑如下初边值问题 {ut-Δu=f(x,t),(x,t)∈Ωc×Ju(x,t)=g(x,t),(x,t)∈Γ×Ju(x,0)=ψ(x),x∈Ωc(1) 其中,ut表示未知函数u(x,t)关于时间t的导数,f(x,t),g(x,t),ψ(x)为满足适当条件的已知函数。我们还假设函数u(x,t)在无穷远处有界。但在此我们不需要在无穷远处的辐射条件。 首先,我们将问题(1)对时间进行离散化。设τ为时间步长,并记tk=k·τ,uk(x)=u(x,tk),zk(x)=ut(x,tk).我们取(1)式中第一式如下离散化形式 zk(x)-Δuk(x)=fk(x)(2)uk(x)=uk-1(x)+τ{(1-β)zk-1(x)+βzk(x)}(3) 其中,β∈(0,1],k=1,2,…,[T/τ].记 λ∶=(√τβ)-1??uk=uk-1+τ(1-β)zk-1??fk∶=-?uk-τβfk 经过简单的代数运算,由方程(2)与(3)可得 Δuk-λ2uk=λ2?fk?x∈Ωc(4)zk=λ2(uk-?uk)(5) 在对时间离散化(1)中的控制方程的同时,对(1)中的边界条件也取离散化形式。结合在无穷远处u(x,t)的有界性,我们可将问题(1)归化为如下离散问题 1)求已知量?uk与?fk(k=1?2???[Τ/τ]). ?uk=uk-1+τ(1-β)zk-1(6)?fk=-?uk-τβfk(7) 2)求解问题 {Δuk-λ2uk=λ2?fk?uk=gk?|uk|＜+∞?x∈Ωcx∈Γ|x|→+∞(8) 3)更新zk的值 zk=λ2(uk-?uk)(9) 由上述可知,问题(1)对时间离散化后,对每一时间步长τ得到关于时间步长的离散化问题(8),而(8)式形式上为椭圆问题,可采用分层方式求解。对每一时间层,先用(6)～(7)式求出?uk与?fk,然后重点求解问题(8)得uk之后,为了求下一时间层问题利用(9)式来更新zk之值,依次重复上述过程,便可得到每一时刻tk时uk之值。 2 基于自然积分码的kc 我们引入一条嵌入Ωc中且半径为R的圆周ΓR,ΓR将Ωc分为两部分:由Γ与ΓR所围成的有界环形区域Ω1;ΓR的圆外无界区域Ω2.由文可知,无界区域Ω2归化的Poisson积分公式与自然积分方程分别为 uk(r,θ)=12π∫2π0[Κ0(λr)Κ0(λR)+2+∞∑n=1Κn(λr)Κn(λR)cosn(θ-θ′)]?uk(R?θ′)dθ′+F(λ?R；?fk?r?θ)?r＞R(10)?uk(R?θ)?n+G(λ?R；?fk?θ)=λ2π∫2π0ˉΚn(λ?R；θ-θ′)?uk(R,θ′)dθ′(11) 其中 上述式中的In(x)与Kn(x)分别为第一类、第二类变型的Bessel函数。自然积分算子Kcλ由下式定义 Kcλv=λ2π∫2π0ˉΚn(λ?R；θ-θ′)?v(R?θ′)dθ′(22) 问题(8)可归结为求解如下问题 {Δuk-λ2uk=λ2?fkuk=gk(x)??uk(R?θ)?n+G(λ?R；?fk?θ)=Kcλ?uk?x∈Ω1x∈ΓΓR上(23) 不失一般性,我们以下假定g(x,t)≡0.记 Η10Γ(Ω1)∶={u∶u∈Η1(Ω1)?u|Γ=0} 问题(23)等价于如下变分问题 {求uk∈Η10Γ(Ω1)?便得对任意vk∈Η10Γ(Ω1)满足D1(uk?vk)+?D2(uk?vk)=?f(vk)(24) 其中 对p0,引入Sobolev空间Hp(ΓR) Ηp