Robin型非重叠区域分解法的收敛性-LSEC
服从OldroydB型微分模型的粘弹性流体问题的数值解法
一ｎ ｒ ｌ ｒｏ ｏ ｎ ｒ ｒ ｅ
Ｎｏ — ｏ ｅｌ ｐ ｉｇ Ｄｏ ｉ ｃ ｍｐ ｓｔ ｎ Ｐ ｏ ｅ ｕ ｎ ｖ ｒａ ｐ ｎ ｍａｎ Ｄｅ ｏ ｏ ｉ ｏ ｒ ｃ ｄ ｒ ｉ ｅ
Ｍｉｅ ｎｔ ｅ ｎ ｘ ｄ Ｆｉｉ Ｅｌｍｅ ｔ ｅ
１ 引言
７ ４
数学理论与应用
片多项式去逼近 ， 并且对于附加应力张量的对流项使用了 Ｌｓｉ — ａｉ ｔ ｅ ｎ Ｒ ｖ ｒ方法［ ，］ ａｔ ａ ４５ ．
本 文考虑 服从 ＯｄｏｄＢ型微分 模型 的粘 弹性 问题． ｌｒｙ 其模 型 的一般 形式 （ ｌｒｙ Ｏｄｏｄ问题 ） 为
‘ ｆ
，
从 上面 的定 义可 知 ， 限元 空 间 有
≤Ｋ ｌＩ ｌｌ≤Ｋ ～ｌｌ Ｖ Ｍ ． ｈ ｌ ｌ ｈ ｌｌ ∈．ｙ ， 。 ｘ，
中的 函数 在 厂上 有一 个 良定 义 的跳跃 ［ ］：
［］ ）＝ ）一 （ 一 ，Ｖｘ∈Ｆ， （ （ ｖｘ ）
其 中 ｖｘ ）ｉ ｉ （ （ ｌ ｍｖ ｘ＋Ａ，） ｔ１． ，
其 中 是 Ｒ 中边 界为 ａ 的开 区域 ， 有分段 一致 光滑 Ｌｐｃｉ 边界 ．模 型 （ ．） 关 于 区域 力 具 ｉｓｈｔ ｚ １ １是
一
周期的． １ １ 中， 在（ ． ） 是流速矢量 ， ， 是给定的已知函数 ； 。 × ＋ 是双线性 ｔ ） ｇ： Ｒ 一
映射： 。 ）＝ （）一∞ 一ａｄｕｔ ｇ（ ， （） （（） ｒ＋ａ （） — ｔ ｕ＋ ｄｕ ）＋ 三 ｒ （ Ｖ
第３ ２卷 第 ３期 ２１ ０ ２年 ９月
数 学 理 论 与 应 用
ＭＡＴ ＨＥＭＡＴ Ｃ Ｉ ＡＬ ＴＨＥ０ＲＹ ＡＮＤ ＰＵＣＡ１ ＯＮＳ ＡＰ １
第五节收敛准则
第五节收敛准则一、收敛准则我们希望，当单元的划分逐渐加密的时候，位移、应变和应力能收敛到精确值，而且收敛得越快越好。
这样就要求所选择的位移模式满足某些条件：1. 位移模式必须包含单元的刚体位移和单元的常应变。
——该条件是收敛准则的必要条件，称为完备条件。
满足该条件的单元称为完备单元。
当节点位移是由某个刚体位移所引起时，弹性体内不会有应变。
这样，位移模式就不但要具有描述单元本身形变的能力，而且还要具有描述由于其它单元形变而通过节点位移引起单元刚体位移的能力。
每个单元的应变一般总是包含着两个部分：一部分是与该单元中各点的位置坐标有关的（即所谓各点的变应变）；另一部分是与位置坐标无关的，即所谓的常应变。
从物理意义上看，当单元尺寸无限缩小时，每个单元中的应变应趋于常数。
除非我们的位移模式包含着这些常应变，否则就没有可能收敛与正确解。
例题6 试证明3节点三角形单元是完备单元。
证明：3节点三角形单元的位移模式123546u x y v x y αααααα⎧⎪⎨⎪⎩=++=++ （3-1）1. 单元刚体运动时，有 0x y x y εεγ=== ，即0,0,0du dv du dv dx dy dy dx==+= 代入式（3-1），得i 图3.18 52630,0,0αααα==+= 则有55533313115342222u y y y y v x ααααααααααααα⎧⎪⎪⎨⎪⎪⎩-+-=+=-+=--=+ （a ）设三角形单元沿x 轴的刚体平移u 0和y 轴的刚体平移v 0，三角形单元绕z 轴作刚体转动w 0角度。
由于w 0引起单元内任意点A(x,y)的位移0000:sin :cos x u r y y v r xωθωωθω=-=-'==' 则A(x,y)点的总位移为000000u u u u y v v v v x ωω⎧⎪⎨⎪⎩=+=-'=+=+' （b ）比较（a ）式和（b ）式，得5310400,,2u v ααααω-===这说明位移模式包含单元的刚体位移。
有序Banach空间分数阶Robin边值问题的正解
有序Banach空间分数阶Robin边值问题的正解李小龙;张丽丽【摘要】讨论了有序Banach空间E中Riemann-Liouville分数阶Robin边值问题:-Dα0+u(t) =f(t,u(t)),0≤t≤ 1,u(0) =u'(1) =θ正解的存在性,其中1＜α≤2,f:[0,1]×P→P连续,P为E中的正元锥.利用非紧性测度的估计技巧及凝聚映射的不动点指数理论获得了该边值问题正解的存在性结果.【期刊名称】《宁夏大学学报（自然科学版）》【年(卷),期】2019(040)002【总页数】5页(P111-115)【关键词】分数阶微分方程;Robin边值问题;正解;凝聚映射;不动点指数【作者】李小龙;张丽丽【作者单位】陇东学院数学与统计学院,甘肃庆阳 745000;陇东学院数学与统计学院,甘肃庆阳 745000【正文语种】中文【中图分类】O175.15分数阶微分方程在流体力学、流变学、生物系统的电传导、各种回路、黏弹性力学、分数控制系统与分数控制器、电分析化学、神经的分数模型以及分数回归模型等领域有广泛应用,特别是在与分形维有关的物理与工程问题中.近年来许多学者应用相关的不动点定理与上下解的单调迭代技巧研究了分数阶边值问题的正解及其多个正解的存在性[1—5],但在一般的无穷维Banach空间中对该类问题的研究还比较少,并且讨论的是分数阶Dirichlet边值问题解的存在性,本文研究了分数阶Robin边值问题正解的存在性.在研究无穷维Banach空间中分数阶微分方程边值问题时,非线性项f的连续性保证不了解的存在性,为了对相应的积分算子应用凝聚映射的拓扑度理论及相关的不动点定理,需要给f附加以下条件:紧型条件或者是耗散型条件,而耗散型条件适用于特殊情形;又Banach空间的微分方程与普通微分方程的最大差异是把微分方程转换为与之等价的积分方程后,相应的积分算子不再具有紧性.为了对该积分算子应用凝聚映射的不动点定理,通常需要给非线性项f附加一些非紧性测度条件.本文使用了下列非紧性测度条件:(H0) 对任意R>0,f(I×PR)有界,且存在常数使得对任意t∈I,D⊂PR,有α(f(t,D))≤Lα(D),其中在研究Banach空间中微分方程的正解时,很多文献(如文献[6])都要求f在有界集上一致连续.文中所利用的新的非紧性测度估计技巧[7]只需要f连续.本文将在一般的有序Banach空间E中讨论非线性分数阶Robin边值问题:(1)正解的存在性,其中是标准的Riemann-Liouville分数阶导数,f:[0,1]×P→P连续,P 为E中的正元锥.1 预备知识设E为有序Banach空间,其正元锥P为正规锥,正规常数为N,记I=[0,1].设C(I,E)为定义于I取值于E的全体连续函数按范数构成的Banach空间.记C(I,P)={u∈C(I,E)|u(t)∈P,t∈I},则C(I,P)为C(I,E)中的正规锥.正规常数亦为N,以下使用的C(I,E)中的半序“≤”由C(I,P)引出.定义1[1] 设α>0,函数f:(0,+∞)→R的α阶Riemann-Liouville积分为其中Γ(·)为Gamma函数.定义2[1] 设α>0,函数f:(0,+∞)→R的α阶Riemann-Liouville导数为其中Γ(·)为Gamma函数,n=[α]+1.由Riemann-Liouville型微分的定义可得下列结论.引理1[1] 设α>0,u∈C(0,1)∩L(0,1)是分数阶微分方程的解,则u(t)具有形式u(t)=c1tα-1+c2tα-2+…+cNtα-N,ci∈R, i=1,2,…,N,其中N是大于或等于α的最小正整数.引理2[1] 假设u∈C(0,1)∩L(0,1)有α(α>0)阶导数属于C(0,1)∩L(0,1),则ci∈R, i=1,2,…,N,其中N是大于或等于α的最小正整数.引理3 设1<α≤2,则对任意h∈C(I,E),Banach空间E中的线性分数阶边值问题(2)存在唯一解u(t)=G(t,s)h(s)dsTh(t),(3)其中证明由文献[1]知的解为再由u(0)=u′(1)=θ,可得于是下面证明唯一性.设u1,u2∈C(I,E)为方程(2)的两个解,则对任意φ∈E*,r(t)=φ(u1(t)-u2(t))为纯量线性方程的解,由文献[3]知r(t)=0,由φ∈E*的任意性知u1(t)-u2(t)≡θ,即u1(t)≡u2(t)于I,因此方程(2)的解唯一.引理4 由(3)式知算子T:C(I,E)→C(I,E)满足证明由(3)式知故从而显然算子T:C(I,E)→C(I,E)为正的线性连续算子,T有相应于第一特征值λ1的正特征函数u*,即λ1Tu*=u*.文中E与C(I,E)中有界集的Kuratiwski非紧性测度均由α(·)表示.对B⊂C(I,E), 记B(t)={u(t)|u∈B}⊂E,t∈I.引理5[8] 设B⊂C(I,E)为等度连续的有界函数族,则α(B(t))在I上连续,且引理6[9] 设B={un}⊂C(I,E)为可列集,若存在ψ∈L1(I)使得则α(B(t))在I上可积,且引理7[7] 设D⊂E有界, 则存在D的可列子集D0,使得α(D)≤2α(D0).定义算子Q:C(I,P)→C(I,P)如下:(Qu)(t)=G(t,s)f(s,u(s))ds,(4)则Q:C(I,P)→C(I,P)连续,且方程(1)的解等价于积分算子Q的不动点.引理8 设f:I×P→P满足假设(H0),则由(4)式定义的算子Q:C(I,P)→C(I,P)为凝聚映射.证明由(4)式易证Q把C(I,P)中的有界集映为有界的等度连续集.任取非相对紧的有界集B⊂C(I,P),下面证明α(Q(B))<α(B). 令则对∀t∈I,B(t)⊂PR,设为假设(H0)中的非紧性测度系数. 由引理7知,存在可列集B1={un}⊂B,使得α(Q(B))≤2α(Q(B1)).故对任意t∈I,由引理6及假设(H0)可得α(Q(B1(t)))=2α({G(t,s)f(s,un(s))|n=1,2,…})ds=2G(t,s)α(f(s,B1(s)))ds≤2LG(t,s)α(B1(s))ds≤2LG(t,s)dsα(B1)≤因为Q(B1)等度连续,由引理5知α(Q(B1))=于是有α(Q(B))≤2α(Q(B1))≤因此Q:C(I,P)→C(I,P)为凝聚映射.取C(I,P)的子锥:K={u∈C(I,P)|u(t)≥θ,∀t∈I},容易证明Q(C(I,P))⊂K,从而当f:I×P→P时, Q:K→K为凝聚映射,方程(1)的正解等价于Q在K中的不动点.本文将用凝聚映射的不动点指数理论寻找Q的不动点.引理9[10] 设E为Banach空间, K为E中的锥, Ω⊂E为有界开集,为凝聚映射.若Q 满足u≠λQu,∀u∈K∩∂Ω,0<λ≤1,则不动点指数i(Q,K∩Ω,K)=1.引理10[11] 设E为Banach空间, K为E中的锥, Ω⊂E为有界开集,为凝聚映射.若存在v0∈K,v0≠θ,使得Q满足u-Qu≠μv0,∀u∈K∩∂Ω,μ≥0,则不动点指数i(Q,K∩Ω,K)=0.2 主要结果及其证明定理1 设E为有序Banach空间,其正元锥P为正规锥,f:I×P→P连续,满足条件(H0).若f满足下列条件之一:(H1) ① 存在ε∈(0,(α-1)Γ(α+1))及δ>0,使得当x∈Pδ时f(t,x)≤εx;② 存在η>λ1及h0∈C(I,P),使得当x∈P时f(t,x)≥ηx-h0(t).(H2) ① 存在ε>λ1及δ>0,使得当x∈Pδ时f(t,x)≥εx;② 存在η∈(0,(α-1)Γ(α+1))及h0∈C(I,P),使得当x∈P时f(t,x)≤ηx+h0(t).则边值问题(1)至少存在一个正解.证明由上面的论述知,只需证明由(4)式定义的凝聚映射Q:K→K存在非零的不动点.取0<r<R<∞,记以下分2种情形分别证明当r充分小R充分大时Q在上存在不动点.情形1 f满足假设(H1). 取0<r<δ,其中δ为假设(H1)中的常数, 证明Q满足引理9中的条件:u≠λQu, ∀u∈K∩∂Ωr, 0<λ≤1.(5)反设(5)式不成立,则存在u0∈K∩∂Ωr及0<λ0≤1,使得u0=λ0Qu0.根据Q的定义及条件(H1)中① 得u0(t)=λ0Qu0(t)≤G(t,s)f(s,u0(s))ds≤εG(t,s)u0(s)ds=εTu0(t).累次使用上式,可得u0(t)≤εTu0(t)≤…≤εnTnu0(t),∀t∈I,n∈N.由锥K的正规性和引理4知其中N为正规常数,故这与矛盾.于是(5)式成立,再由引理9知i(Q,K∩Ωr,K)=1.(6)下面证明当R充分大时u-Qu≠τu*，∀u∈K∩∂ΩR, τ≥0.(7)反设(7)式不成立，则存在u0∈K∩∂ΩR及τ0≥0,使得u0-Qu0=τ0u*,则u0=Qu0+τ0u*,根据算子Q的定义及条件(H1)中② 得u0=Qu0+τ0u*=G(t,s)f(s,u0(s))ds+τ0u*≥ηG(t,s)u0(s)ds-G(t,s)h0(s)ds+τ0u*=ηTu0-G(t,s)h0(s)ds+τ0u*.从而有(ηT-I)u0≤G(t,s)h0(s)ds-τ0u*≤G(t,s)h0(s)ds.又由η>λ1知(ηT-I)为正算子,故逆算子(ηT-I)-1存在, 由锥K的正规性得(8)取则(7)式成立,由引理10知i(Q,K∩ΩR,K)=0,从而根据不动点指数理论的区域可加性,由该式结合(6)式可得i(Q,K∩ΩR,K)-i(Q,K∩Ωr,K)=-1.由可解性知Q在中至少存在一个不动点,该不动点即为方程(1)的正解.情形2 f满足假设(H2). 取0<r<δ,证明u-Qu≠τu*, ∀u∈K∩∂Ωr, τ≥0.(9)反设(9)式不成立,则存在u0∈K∩∂Ωr及τ0≥0,使得u0-Qu0=τ0u*,从而u0=Qu0+τ0u*.令τ*=sup{τ|u0≥τu*},即0<τ0<τ*<+∞,且u0≥τ*u*.又由T的正性知λ1Tu0≥τ*λ1Tu*=τ*u*.由条件(H2)中① 可得u0=Qu0+τ0u*=G(t,s)f(s,u0(s))ds+τ0u*≥εG(t,s)u0(s)ds+τ0u*=εTu0+τ0u*≥λ1Tu0+τ0u*≥(τ*+τ0)u*.这与τ*的定义矛盾.故根据引理10知i(Q,K∩Ωr,K)=0.(10)再证当R充分大时u≠λQu, ∀u∈K∩∂ΩR, 0<λ≤1.(11)假设存在u0∈K及0<λ0≤1,使得u0=λ0Qu0.从而由条件(H2)中② 得u0=λ0Qu0≤G(t,s)f(s,u0(s))ds≤G(t,s)(ηu0(s)+h0(s))ds=ηTu0+G(t,s)h0(s)ds，即(I-ηT)u0≤G(t,s)h0(s)ds,又所以由微扰定理知I-ηT存在有界逆算子(I-ηT)-1,且从而由锥K的正规性得取ε0>0,使得η+ε0<(α-1)Γ(α+1),则有而收敛,即级数收敛,令则有取则(11)式成立,由引理9知i(Q,K∩ΩR,K)=1.于是,由该式结合(10)式可得i(Q,K∩ΩR,K)-i(Q,K∩Ωr,K)=1.由可解性知Q在中至少存在一个不动点,该不动点即为方程(1)的正解. 参考文献：【相关文献】[1] BAI Zhanbing, LÜ Haishen.Positive solutions for boundary value problem of nonlinear fractionl differential equation[J].J Math Anal Appl,2005,311(2):495-505.[2] BAI Zhanbing.On positive solutions of a nonlocal fractional boundary valueproblem[J].Nonlinear Anal,2010,72(2):916-924.[3] JIANG Daqing, YUAN Chengjun.The positive properties of the Green function for Dirichlet-type boundary value problems of nonlinear fractional differential equations and its application[J].Nonlinear Anal,2010,72(2):710-719.[4] WANG Yingqing, LIU Lishan, WU Yonghong. Positive solutions for a nonlocal fractional differential equation[J].Nonlinear Anal,2011,74(11):3599-3605.[5] LIN Legang, LIU Xiping, FANG Haiqin.Method of upper and lower solutions for fractional differential equations[J].Electronic J Differential Equations,2012,100:1-13. [6] 王永庆,刘永山.Banach空间中分数阶微分方程m点边值问题的正解[J].数学物理学:A,2012,32(1):246-256.[7] 李永祥.抽象半线性发展方程初值问题解的存在性[J].数学学报,2005,48(6):1103-1108.[8] 郭大钧,孙经先.抽象空间常微分方程[M].济南:山东科学技术出版社,1989:188-222.[9] HEINZ H R. On the behaviour of measure of noncompactness with respect to differentiation and integration of vector-valued functions[J].NonlinearAnal,1983,7(12):1351-1371.[10] 郭大钧.非线性泛函分析[M].济南:山东科学技术出版社,1985:234-353.[11] 余庆余.半序Banach空间中凝聚映射及其正不动点[J].兰州大学学报:自然科学版,1979,15(3):1-5.。
并行矢量有限元法分析新型波导侧面馈电天线
并行矢量有限元法分析新型波导侧面馈电天线彭文峰;宛汀;郭继承【摘要】采用并行矢量有限元区域分解法对一种新型波导侧面馈电天线进行有效的分析。
计算中采用矢量棱边元消除传统节点有限元存在的伪解问题。
当所分析的天线规模较大时，在普通单台计算机上采用传统有限元方法分析会面临着内存不足和效率不高的问题，引入一种非重叠型矢量有限元区域分解法有效地克服了这一问题。
该方法将原始大型求解区域划分成一系列小的子区域单独求解，具有很高的可并行性，从而大大缓解了内存需求，提高了计算效率。
通过对一种新型波导侧面馈电天线的分析，验证了这种方法的精确性和有效性。
%A parallel vector finite element domain decomposition method is employed to analyze a new-style antenna with waveguide side feed. In this method, edge-based vector elements are employed to avoid the spurious solutions in traditional method of node-based elements. However, on single personal computer, the finite element solution of large-scale problems faces the problems of low operating efficiency and insufficient physical memory. A non-overlapping domain decomposition method is introduced to overcome these problems. The original large domain is divided into a series of small sub-domains, which can be analyzed independently. Hence, this method can easily be parallel implemented to improve efficiency and relieve memory requirement. Numerical experiments of new-style antenna with waveguide side feed demonstrate the accuracy and efficiency of this method.【期刊名称】《电波科学学报》【年(卷),期】2012(027)006【总页数】7页(P1110-1115,1128)【关键词】矢量有限元;非重叠型区域分解法;并行技术;波导侧面馈电天线【作者】彭文峰;宛汀;郭继承【作者单位】南京理工大学通信工程系,江苏南京210094;南京理工大学通信工程系,江苏南京210094;南京理工大学通信工程系,江苏南京210094【正文语种】中文【中图分类】TN925引言天线在现代通信领域中具有广阔的应用前景，因此，天线的分析与设计引起了越来越多的关注。
《概率论四种收敛性》PPT课件
第三章 3.1四种收敛性
1
主要内容
车贝晓夫不等式 几乎处处收敛 依概率收敛 依分布收敛 r-阶收敛
X
E(X)
2)
D(X)
2
2 22
1 2
23
证明：已知Xi (i 1, 2, , n)相互独立，且方差有限
证明lim
P
n
1 n
n i 1
Xi
1 n
n i 1
E( Xi )
1
证明：设随机变量Z
X
1 n
n i 1
Xi ,
1 n
1
n
1n
E(Z)
E(X )
E( n
i 1
Xi )
n
E(
的期望的偏差不小于 的概率的估计式 .
如取 3
P{|
X
E( X ) |
3 }
2 9 2
0.111
可见，对任给的分布，只要期望和方差 2存在，
则 r.v X取值偏离E(X)超过 3 的概率小于0.111 .
6
车贝晓夫不等式的用途：
（1）证明大数定律；（2）估计事件的概率。
车贝晓夫不等式只利用随机变量的数学期望及方差就可对的 概率分布进行估计。
( x E( X ))2dF( x)
xE( X )
2dF( x) xE( X )
2P X E(X)
从而P(
X
E( X )
)
D( X )
关于奇摄动robin边值问题的几个定理
关于奇摄动robin边值问题的几个定理随着科学技术的发展，奇摄动robin边值问题也受到了广泛的关注，并成为研究者们需要解决的一个重要问题。
该问题涉及了一些重要的数学定理，其中主要涉及到几个定理，其中最为重要的有Liouville定理，Caccioppoli定理和Rellich-Kondrachov定理，它们在解决奇摄动robin边值问题中均扮演重要角色。
首先，我们介绍Liouville定理，又称Liouville-Neumann定理。
它是一个把有限区域外部源的能量从内部传至外部的关系，其主要的表达式为：V(x)*u(x) = S(x)其中V(x)是robin边值中的一个常数，S(x)表示区域内部的源，u(x)表示u(x)的梯度；此外，当V(x)=0时，公式约化为：u(x) = 0这个定理可以有效地处理奇摄动robin边值问题，它实质上是在一个紧张的区域内求解某些不定方程的问题。
其次，我们来讨论Caccioppoli定理。
它的核心概念是利用一个所谓的Caccioppoli方程来描述传热方程的解，即：α2u2 +2u2 +2u2 = 0其中α，β，γ都是常数，其中α表示温度梯度，β表示声速梯度，γ表示吸收率。
由于Caccioppoli定理可以非常有效地求解不定方程，因此它被广泛用于奇摄动robin边值问题。
最后，我们来谈谈Rellich-Kondrachov定理。
它是一种利用函数间隙和函数梯度来描述某一单元的解的定理。
其主要表达式为：u(x) =u(x)其中λ是一个常数，它表示某一单元内的解的空间变形系数。
通过利用Rellich-Kondrachov定理，人们可以更有效地求解奇摄动robin边值问题。
综上所述，Liouville定理，Caccioppoli定理和Rellich-Kondrachov定理是研究奇摄动robin边值问题的重要理论基础，在解决问题时可以极大地提高计算效率，有助于我们进一步了解该问题。
解带Robin边界条件的变分不等式的区域分解算法
解带Robin边界条件的变分不等式的区域分解算法
曾金平;陈高洁
【期刊名称】《系统仿真学报》
【年(卷),期】2007(19)17
【摘要】针对一类带Robin边界条件的椭圆型变分不等式问题,构造基于Robin 内边界传输条件的非重叠加性区域分解算法,并建立了算法的收敛性。
这类区域分解算法广泛应用于求解偏微分方程边值问题并取得了一系列收敛性结果。
数值结果表明,基于Robin边界传输条件的区域分解法可通过调节内边界传输条件中的Robin参数,来加快算法的收敛速度。
【总页数】3页(P3949-3950)
【关键词】Robin条件;区域分解算法;变分不等式;收敛性
【作者】曾金平;陈高洁
【作者单位】东莞理工学院软件学院;湖南大学数学与计量经济学院
【正文语种】中文
【中图分类】O211
【相关文献】
1.大地电磁非重叠型区域分解算法子域边界条件比较 [J], 李丹
2.解Stokes问题的区域分解算法 [J], 顾金生;胡显承
3.解四阶椭圆问题的区域分解算法——多子域重叠情形 [J],
4.解椭圆变分不等式问题的区域分裂与异步并行算法 [J], 芮洪兴
5.带约束广义变分不等式问题的一般分解算法的收敛性分析(英文) [J], 鲁其辉;朱道立
因版权原因，仅展示原文概要，查看原文内容请购买。
有限元边界积分结合撕裂对接法分析电磁散射
（３）
根据变分原理可知，内部区域中电场的泛函满足
ｔ？Ｆ（Ｅ）＝Ｏ，其中
舢，＝壶Ⅲ，吉ｃＶ㈣·（Ｖ Ｘ Ｅ）ｄｖ～
Ⅲ，砍ｒＥ．ＥｄＶ＋。Ｅ矗ｄＳ
（４）
式中，Ｊ。＝蠢×Ｈ＝痞×槲，将Ｅ和Ｊ。分别用四面体棱边基 函数以及表面ＲＷＧ（Ｒａｏ－Ｗｉｈｏｎ－Ｇｌｉｓｓｏｎ）基函数来展开
Ｅ一∑Ｗ。Ｅ。，Ｊｓ＝∑Ｊ。厂。
了硬件资源不足的问题．提高了求解问题的效率。该方法在处理内部未知量大和内部材料复杂的电大尺寸问题
方面具有较大优势。给出的数值算例充分证明了该方法的可行性和有效性。
关键词：有限元边界积分；有限元撕裂对接法；并行技术；电磁散射
中图分类号：ＴＭ １５
文献标志码：Ａ
ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１００１—５０６Ｘ．２０１０．０９．１５
万方数据
第９期
宛汀等：有限元边界积分结合撕裂对接法分析电磁散射
的产生口…。其中边界积分部分还能引入多层快速多极子 技术来提高其求解效率。目前，有限元边界积分方法已经 被广泛地应用于电磁散射、辐射、传输以及电磁兼容等问题 的分析中，是一种功能强大的数值分析方法。
计算电磁学领域分析问题的效率不仅取决于好的数值 分析方法，还受到硬件资源的限制。随着科技进步和工程 应用要求的提高，计算电磁学领域分析问题的电尺寸越来 越大，求解问题规模的增大给计算机的运行速度和内存容 量带来了巨大挑战。单台ＰＣ机已经不能满足大规模计算 的需求，ＦＥＢＩ方法也同样面临着这个问题。采用多机并联 的并行机群来扩展计算机的内存并提高计算效率，是解决 这一问题的有效方法［７］。区域分解法（ｄｏｍａｉｎ ｄｅｃｏｍｐｏｓｉ－ ｔｉｏｎ ｍｅｔｈｏｄ，ＤＤＭ）是一种极其适合于进行大规模数值分 析及并行计算的方法。它采用“分而治之”的思想．把原始 的待求区域划分成一系列较小的子区域来求解，缩减了计 算规模，降低了内存需求。同时该方法具有很高的并行度， 这一特性使得人们可以很方便地将并行算法引入其 中¨１…。本文将有限元区域分解法中的撕裂对接法（ｆｉｎｉｔｅ ｅｌｅｍｅｎｔ ｔｅａｒｉｎｇ ａｎｄ ｉｎｔｅｒｃｏｎｎｅｃｔｉｎｇ，ＦＥＴＩ）引入到ＦＥＢＩ方 法的求解中，对于ＦＥＭ部分的处理采用ＦＥＴＩ方法，极易 在并行平台上实现；对于ＢＩ部分的处理采用并行多层快速 多极子技术（ｍｕｌｔｉｌｅｖｅｌ ｆａｓｔ ｍｕｈｉｐｏｌｅ ｍｅｔｈｏｄ，ＭＩ。ＦＭＭ）。 整套算法的高度并行性使得这种方法能够分析电大尺寸问 题，尤其适合于分析内部未知量大且材料复杂的问题。文 中给出的算例的数值分析结果充分证明了该方法的正确性 和有效性。
热传导方程有限差分法的MATLAB实现
万方数据万方数据万方数据万方数据热传导方程有限差分法的MATLAB实现作者：史策作者单位：西安建筑科技大学,理学院,陕西,西安,710055刊名：咸阳师范学院学报英文刊名：JOURNAL OF XIANYANG NORMAL UNIVERSITY年，卷(期)：2009，24(4)被引用次数：0次1.曹钢,王桂珍,任晓荣.一维热传导方程的基本解[J].山东轻工业学院学报,2005,19(4):76-80.2.万正苏,方春华,张再云.关于热传导方程有限差分区域分解并行算法精度的注记[J].湖南理工学院学报(自然科学版),2007,20(3):12-14.3.StephenJ.Chapman.MATLAB编程[M].邢树军,郑碧波,译.北京:科学出版社,2008.4.田兵.用MATLAB解偏微分方程[J].阴山学刊,2006,20(4):12-13.5.王飞,裴永祥.有限差分方法的MATLAB编程[J].新疆师范大学学报(自然科学版),2003,22(4):21-27.6.王宝红.热传导方程的可视化探讨[J].忻州师范学院学报,2008,24(2):31-36.7.李先枝.热传导方程差分解法的最佳网格[J].河南大学学报(自然科学版),2004,34(3):16-18.8.赵德奎,刘勇.MATLAB在有限差分数值计算中的应用[J].四川理工学院学报,2005,18(4):61-64.9.谢焕田,吴艳.拉普拉斯有限差分法的MATLAB实现[J].四川理工学院学报,2008,21(3):1-2.10.南京大学数学系计算数学专业.偏微分方程数值解法[M].北京:科学出版社,1979.1.学位论文申卫东热传导方程有限差分区域分解算法研究2003区域分解算法是在并行机上求解偏微分方程数值解的一种较自然的方法.该方法先将偏微分方程求解区域划分为若干个子区域,然后在各个子区域并行求解.全文共五章.第一章为引言,简要介绍了热传导方程并行算法的概况及该文所讨论的基本内容.在第二章,我们在内边界点为等距分划的多子区域条件下,得到Dawson等人关于求解热传导方程区域分解算法差分解的误差估计.在第三章,我们以Saul'yev非对称格式作内边界处理,发展了新的区域分解算法,得到了差分解的先验误差估计,并与Dawson等人的算法作了比较.给出了关于算法计算精度的数值结果.在第四章,我们发展了一些新技术,在子区域的边界处采用小时间步长古典显式格式求解,构造了新的区域分解算法,得到了差分解的先验误差估计.给出了关于算法计算精度的数值结果.在第五章,我们在二维热传导方程求解上扩充了Dawson等人的区域分解算法.给出了关于算法计算精度的数值结果.第六章为该研究工作的主要结论.2.期刊论文张守慧.王文洽.ZHANG Shou-hui.WANG Wen-qia热传导方程有限差分逼近的数学Stencil及其新型迭代格式-山东大学学报（理学版）2006,41(6)将Stencil应用于偏微分方程有限元差分逼近过程,以两类差分格式为基础建立了求解热传导方程的两种新型迭代算法.此两种算法与经典的Jacobi方法同样具有并行的性质,但比Jacobi方法收敛快.给出的算例说明方法的适用性.3.期刊论文吕桂霞.马富明.Lü Guixia.Ma Fuming二维热传导方程有限差分区域分解算法-数值计算与计算机应用2006,27(2)本文讨论了一类数值求解二维热传导方程的并行差分格式.在这个算法中,通过引进内界点将求解区域分裂成若干子区域.在子区域间内界点上采用非对称格式计算,一旦这些点的值被计算出来,各子区域间的计算可完全并行.本文得到了稳定性条件和最大模误差估计.它表明我们的格式有令人满意的稳定性,并且有着较高的收敛阶.4.学位论文田源地下煤火三维数理模型正演数值模拟2006本文首先给出了几个地下煤火随空间、温度变化的动态和稳态热数学物理模型及其简化模型。
偏微分方程的Robin型非重叠区域分解方法的开题报告
偏微分方程的Robin型非重叠区域分解方法的开题报告一、研究背景及意义偏微分方程是自然界中许多现象的数学描述，如物理过程，化学反应，经济变化等。
由于其广泛应用，使得人们对数值解法的求解效率提出了更高的要求。
而分块法作为一种分治法，具有分解区域的优点，能够有效提高数值解的计算速度。
在偏微分方程的求解中，非重叠区域分解方法是一种常用的分块法，并且该方法在计算过程中考虑了边界条件的影响，使得其求解结果更加准确。
因此，深入研究非重叠区域分解方法在偏微分方程求解中的应用，对于提高数值解的计算速度和精度有着重要的意义。
本文将研究偏微分方程的Robin型非重叠区域分解方法及其在数值求解中的应用。
Robin型边界条件常常出现在物理过程中，如热传导，扩散，传感器等问题，因此研究该问题对于物理学，化学和工程学领域有着重要的意义。
二、研究内容及方法本文将探究偏微分方程的Robin型非重叠区域分解方法，该方法采用分解区域来提高计算效率和准确性。
具体研究内容包括：1. 研究偏微分方程Robin型边界条件的特点及其在数值求解中的应用。
2. 探讨非重叠区域分解方法在偏微分方程数值求解中的基本原理。
3. 针对Robin型偏微分方程，建立相应的非重叠区域分解方法数值模型，并对不同情况下的Robin型边界条件进行数值求解。
4. 对数值求解结果进行分析和比较，在典型问题中验证非重叠区域分解方法在Robin型偏微分方程求解中的有效性和优越性。
本文将采用数学理论分析和程序实现相结合的方法，运用计算机编程语言和数学仿真软件进行偏微分方程Robin型非重叠区域分解方法的求解和结果分析，以期达到深入研究该问题的目的。
三、预期成果及意义通过对偏微分方程Robin型非重叠区域分解方法的研究，我们可以得到以下预期成果：1. 研究分析偏微分方程的Robin型边界条件及其在物理问题中的应用特点，为相关学科领域提供理论支持。
2. 建立基于Matlab等编程语言的非重叠区域分解方法数值模型，提高数值计算求解效率和精度。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Robin型非重叠区域分解法的收敛性
秦理真
摘要
本文对Robin型非重叠区域分解法的收敛性作了系统的研究与总结.此方法在本文中被用来求解二阶线性椭圆方程.
自H.A.Schwarz提出以来,区域分解方法方法得到了极大的发展,现已成为偏微分方程求解的最有效的方法之一.本文所研究的Robin型非重叠区域分解法最初由P.L.Lions首次提出.作为一种非重叠型区域分解法,此方法以Robin条件为交界面传递条件,因此被称为Robin型方法.近年来,此方法已被推广而应用到其它类型的方程上去.
Robin型方法在子区域之间传递的是Robin条件.它由迭代解在子区域交界面上的迹及关于交界面的外法向导数线性组合而得.此组合含有一个松弛参数λ.大量数值分析表明,Robin型方法的收敛率对于λ十分敏感.自Lions提出此方法以来,此方法的收敛率及最优参数的选取一直是令人关注的问题.本文将对此问题作系统的研究,并用数值实验来验证理论分析的结果.
对于连续问题的Robin型方法,本文遵循Deng的方法给出一个与Deng略有不同的证明.且本文还将给出两子区域情形的反例,来证明连续问题的Robin型方法在任何范数下都不可能几何收敛.
对于有限元离散问题的Robin型方法,本文改进了以往文献中的收敛率分析.当方程中低阶项系数b∈[O(H−2),O(h−2)]时,取λ=O(h−1/2b1/4),收敛率为1−O(h1/2b1/4).当b≥O(h−2)时,取λ=O(bh),收敛率为1−δ,其中δ为不依赖于h, H和b的常数.当环绕数N=1,且b∈[0,O(H−2)]时,取λ=O(h−1/2H−1/2),收敛率为1−O(h1/2H−1/2).根据我们对两子区域情形所作的细致的理论分析,可知以上结果不可改进.另外,当N>1,且b∈[0,O(H−2)]时,我们给出了此方法的两个上界1−O(h1/2H1/2b1/2)和1−O((C0)N h1/2H−1/2).以上h为有限元网格尺寸,H 为子区域尺寸,N为区域分解的环绕数,本文将给出它们的定义.需要指出的是,当N>1,且b=0时,Robin型方法的收敛率分析很困难.本文所给的此情形下的收敛率估计是这方面的第一个结果.
为了对收敛率的下界进行估计,本文系统地对D－N算子和离散D－N算子的作了谱分析.与以前这方面的工作相比,本文的工作更为细致而又不失一般性.
本文系统地介绍了Robin型方法的收敛性分析所采取的三种常用方法:直接求解法,能量估计法和改进的能量估计法.此三种方法的特点将在适当的应用中得以体现.
关键词:有限元,Robin型区域分解法,二阶椭圆问题,几何收敛,收敛率。