一类求解反应扩散方程的Newton波形松弛方法

对反应扩散方程提出一种新型的Newton波形松弛方法,并给出此方法的误差估计式.通过与传统的波形松弛方法比较,这种Newton波形松弛方法有更快的收敛性,且收敛速度不随网格加密而减慢.这种方法可以保持传统波形松弛方法可并行的特点.最后通过数值算例验证这种方法的有效性.     

基于非结构自适应网格的二维Euler方程数值求解方法研究

提出了一种基于非结构自适应网格的二维Euler方程的数值解法.采用有限体积法进行空间离散,通量计算采用Jamson中心格式,使得它适用于任意多边形计算单元.为了得到定常解,采用一种显式的四步Runge-Kutta迭代方法对时间进行积分.根据流场参数的变化梯度确定加密边,由加密准则进行自适应网格剖分,然后得到分布合理的加密过后的网格.求解二维Euler方程,对NACA0012翼型进行了数值模拟,通过对自适应前后的数值解的对比,说明所建立的方法是正确的.     

二维三温热传导方程的并行自适应多重网格算法求解

3-T heat conduct equation including electron, ion and photon (radiation) temperatures can be used to approximately describe the energy broadcast across multimedia for radial flow dynamics and discover the energy swapping among photon,electron and ion. Owing to the strong nonlinear diffusion coefficients and energy swapping coefficients and strong discontinuous coefficients across media interfaces,this equation is difficult to be solved with high numerical resolution. Based on the parallel adaptive multigrid software framework UG on 2-D unstructured grid, this paper successfully solved such equation with high resolution by combining the finite volume implicit discretization scheme and parallel adaptive multigrid algorithm, and gained much significant results.     

二维空间中分数阶扩散方程的变网格有限元方法

针对二维空间分数阶偏微分方程,给出了一个变网格全离散有限元格式,并得到了相应最优误差估计.其主要思想是对空间变量采用有限元离散,对时间交量采用差分,但不同时刻的有限元网格可以不同.这对于没计相应的自适应算法是十分有益的.     

二维奇异摄动对流扩散方程的自适应移动网格算法

针对二维奇异摄动对流扩散方程,在任意网格下给出了经典的迎风有限差分格式.利用二元多项式插值技术,推导出一阶最大范数的后验误差估计,并以此设计了一个自适应网格生成算法.数值实验表明本文构造的自适应移动网格算法是有效的.     

四阶分数阶扩散波动方程的两网格混合元快速算法

本文研究四阶分数阶扩散波动方程模型的基于新混合元方法的快速两网格算法.讨论该方法的稳定性,推导三个未知函数的$L^2$模意义下的最优误差估计.最后通过数值例子验证两网格混合元算法的高效性和理论结果的正确性.     

时间分数次扩散方程反演源项问题的迭代正则化方法

时间分数次扩散方程中反演源项问题是一类经典不适定问题.本文构造了一种新的迭代格式作为正则化方法,给出了先验和后验参数选取下相应的收敛性分析.数值算例验证该方法的有效性.     

半线性Sobolev方程的H~1-Galerkin混合有限元方法

利用H~1-Galerkin混合有限元方法研究了一维半线性Sobolev方程,得到了半离散解的最优阶误差估计,优点是不需验证LBB相容性条件.     

一类求解抛物型方程侧边值问题的最优滤波方法

该文考虑一类特殊的抛物型方程侧边值问题,即一类含有对流项的非标准逆热传导问题. 给定在x=1处的温度测量值来确定区间(0,1)上的未知解u(x, t). 这是一类不适定问题,即问题的解(如果解存在)不连续依赖于数据.为了求解这一问题, 必须采用某些正则化技巧. 该文给出了一种最优滤波方法, 使得问题的真实解和近似解之间的误差估计达到了Hölder型最优. 同时还证明了问题的解在x=0处的收敛性.     

广义神经传播方程的一种修正混合有限元方法的误差分析

利用修正的H~1-Galerkin混合有限元方法研究了广义神经传播方程,论证了其半离散解的存在唯一性,得到了半离散解的最优阶误差估计,该方法的优点是不需验证LBB相容性条件.     

AN OSCILLATION CRITERION FOR A DISCRETE ELLIPTIC EQUATION

ＡＮＯＳＣＩＬＬＡＴＩＯＮＣＲＩＴＥＲＩＯＮＦＯＲＡＤＩＳＣＲＥＴＥＥＬＬＩＰＴＩＣＥＱＵＡＴＩＯＮ￥ＳｕｉＳｕｎＣｈｅｎｇ（ＴｓｉｎｇＨｕａＵｎｉｖｅｒｓｉｔｙ，ＴＡＩＷＡＮ，３００４３）Ａｂｓｔｒａｃｔ：Ａｎｏｓｃｉｌｌａｔｉｏｎｃｒｔｅｒｉｏｎ...     

具有散度形式的拟线性椭圆型方程解的振动性质

利用偏Riccati变换和Young不等式技巧,得到了具有散度形式的拟线性椭圆型方程振动的一些新的充分准则。特别地,推广二阶常微分方程振动的Fite-Wintner定理与Kamenev定理到椭圆型方程。     

FORCED OSCILLATION FOR GENERAL DIFFERENTIAL EQUATIONS

New oscillation criteria for general differential equations of the form x(n)(t) Pn-1(t)x(n-1)(t) … p1(t)x'(t) p0(t)x(t) q1(t)xμ(t)=q2(t)xλ(t) e(t), whereλ,μare the ratios of positive odd integers, 0 <μ< 1 andλ> 1 are established.     

数值积分对完全非线性抛物型积分微分方程有限元方法的影响*

抛物型积分微分方程在带有记忆性的热传导 ,扩散 ,生物力学等实际问题中有广泛的应用 本文考虑如下模型 :ｃ(ｘ ,ｕ) ｔ = · {ａ(ｘ ,ｕ) ｕ + ∫ｔ0 ｂ(ｘ ,ｔ,τ ,ｕ(ｘ ,τ) ) ｕ(ｘ ,τ)ｄτ}+ ｆ(ｘ ,ｔ ,ｕ)　　　　　　　　　　　　　　　　 (ｘ ,ｔ) ∈Ω× [0 ,Ｔ]ｕ(ｘ ,ｔ) =0 ,　　 (ｘ ,ｔ)∈ Ω× [0 ,Ｔ]ｕ(ｘ ,0 ) =ｕ0 (ｘ) ,　　ｘ ∈Ω ,其中Ω Ｒｎ 为多角型区域 , Ω为边界 不用数值积分 ,对这类问题有限元方法研究已有很多工作[5 -8,11] ,导出了最优Ｌ2 及Ｈ1模估计 众所周知 ,解偏微分方程的有限元方法最终…     

求二维半线性椭圆问题混合元解的二重网格法

1引言求解偏微分方程的混合元方法在实际问题(如油藏和地下水的数值模拟)中有广泛的应用.与标准的有限元相比,混合元方法可以同时求解两个变量即压力和流速的近似,而且能保持问题的局部守恒性,从而得到了更好的结果(如Arnold).现已有许多文献使用混合元法求解线性的(如),拟线性的(如Milner)和非线性的(如Park)二阶椭圆问题.     

Riccati技巧与半线性椭圆型方程的振动性质

对二阶半线性椭圆型方程 sum from i,j=1 to N Di[A_(ij)(x)D_jy]+q(x)f(y)=0其中q(x)在外区域Ω R~N上变号。本文建立了一些新的振动性定理,所得结果推广和改进Kamenev,Philos和Li等人的振动准则。     

二阶椭圆问题的混合元方向交替法

本文研究了一般的二阶椭圆问题的混合元方向交替法，给出了两种迭代格式即Uzawa格式和Arrow－Hurwitz格式，并就系数是常数的情形给出了谱分析．本文的结果对油藏模拟有一定的理论和实际意义．     

OSCILLATION CRITERIA FOR NONLINEAR SECOND ORDER ELLIPTIC DIFFERENTIAL EQUATIONS

1.IntroductionOscillationtheoryforellipticdifferentialequationswithvariablecoefficielltshasbeenextensivelydevelopedinrecentyearsbymanyauthors(seee-g-,Allegretto[1,2l,Bugir['],Fiedler[5]'KitamuraandKusano['],Kural4]KusanoandNaitol1o]'NaitoandYoshidal12]'NoussairandSwanson[l3'14]'Swansonl15'l6]andthereferencescitedtherein).Inthispaper,weareconcernedwiththeoscillatorybehaviorofsolutionsofsecondorderellipticdifferelltialequatiollswithalternqtingcoefficients.Asusual,pointsinn-dimensionalEucli…     

二阶椭圆问题的混合元方法

１．引 言 令 是有界区域,边界 充分光滑．Ｓｏｂｏｌｅｖ空间 是熟知的．引入Ｑ＝ Ｈ（ｄｉｖ;Ω）,Ｕ＝ Ｈ１（Ω）,内积和范数记为而 是 的半范．令 ,其范数为 ． 考虑如下二阶椭圆问模型题：由问题（０．１）的位移有限元解通过求导的方法来求ｐ的近似解,会带来额外的舍入误差．应用Ｂａｂｕｓｋａ－Ｂｒｅｚｚｉ混合元法［２］则可得到ｐ足够精度的逼近解．但是,该方法要求离散Ｋ－椭圆性和Ｉｎｆ－Ｓｕｐ不等式同时成立,使得混合元的构造或自由度的选取变得相当复杂［２,１２－１４］．通过“增补”办法,能够克服Ｋ－椭圆性…     

SOME n-RECTANGLE NONCONFORMING ELEMENTS FOR FOURTH ORDER ELLIPTIC EQUATIONS

In this paper, three n-rectangle nonconforming elements are proposed with n ≥ 3. They are the extensions of well-known Morley element, Adini element and Bogner-Fox-Schmit element in two spatial dimensions to any higher dimensions respectively. These elements are all proved to be convergent for a model biharmonic equation in n dimensions.