非线性抛物方程耦合的离散化及其误差分析

0　引　　言对于线性抛物型初边值问题的有限元与边界元耦合法,我们已作过研究(可参见[2],[3]).然而,许多实际问题,如流体、对流扩散、热辐射及热传输等涉及到非线性问题,因此研究非线性问题的数值方法显得尤为重要.近年来,G.N.Gatica与G.C.Hsiao已将有限元(FEM)与边界元(BEM)耦合法拓广到非线性椭圆问题(如[5],[6]),但如何应用FEM与BEM耦合法来处理非线性抛物型初边值问题,就作者所知迄今为止尚属空白.这里我们试图对此进行研究.设Ω为R2中一有界单连通区域,其边界为Γ:=Ω.Ωc:=R2\Ω为闭区域Ω的补区域,I:=(0,T].我们考虑如…     

BOUNDARY INTEGRAL FORMULAS FOR ELASTIC PLANE PROBLEM OF EXTERIOR CIRCULAR DOMAIN

After the stress function and the normal derivative on the boundary for the plane problem of exterior circular domain are expanded into Laurent series, comparing them with the Laurent series of the complex stress function and making use of some formulas in Fourier series and the convolutions, the boundary integral formula of the stress function is derived further. Then the stress function can be obtained directly by the integration of the stress function and its normal derivative on the boundary. Some examples are given. It shows that the boundary integral formula of the stress function is convenient to be used for solving the elastic plane problem of exterior circular domain.     

BOUNDARY INTEGRAL FORMULA OF ELASTIC PROBLEMS IN CIRCLE PLANE

IntroductionConcerningtheelasticplaneprobleminaunitcircle ,ZhengShenzhouandZhengXueliangdevelopedaboundaryintegralformulaofthestressfunction[1]:Φ(r,θ) =-( 1 -r2 ) 24π ∫2π0ν( φ)1 -2rcos(θ-φ) r2 dφ　　 12π∫2π011 -2rcos(θ-ω) r2 dω∫2π0μ( φ)1 -cos(ω-φ) dφ　　 1 -r22π∫2π0μ( φ)1 -2rcos(θ -φ) r2 dφ　　 ( 0 ≤r <1 ) ,( 1 )whereμ(θ) =Φ(r,θ) ｜r=1,ν(θ) = Φ n r=1= Φ r r=1.Intheformula ( 1 )theseconditemisastrongsingularintegral,itshouldbeunderstoodasanintegra…     

波动方程基于自然边界归化的区域分解算法

提出了无界区域波动方程的区域分解算法,基于自然边界归化,分别研究了重叠型与非重叠型区域分解算法,首先将控制方程对时间进行离散化,得到关于时间步长离散化格式,对每一时间步长给出了Dirichlet-Neumann和Schwartz交替算法,对Schwartz交替算法,给出了算法的收敛性,对圆外区域研究了压缩因子,并给出了数值例子。