非线性对流扩散方程的逆风型差分格式

§1 对流扩散方程可以用来描述水中和大气中污染物质的分布、流体流动和流体中的传热等,以上均有对流扩散的特征。数值求解对流扩散方程很重要,近年来工作不少,如[1,2]。我们考虑二维非线性对流扩散方程的初边值问题     

双曲型守恒律组的一类差分格式及其熵条件

其中u(x,t)=(u_1(x,t),…,u_m(x,t))~T,f(u(x,t))=(f_1(u(x,t)),…,f_m(u(x,t)))~T,f的Jacobian记为 A(u)=?f(u)/?u,具有m个实特征值λ_1(u)≤λ_2(u)≤… ≤λ_m(u)以及完备的古特征向量系{γ_k(u)}_k~m=1.对区域R~+={(x,t)|x∈(-∞,+∞),t∈     

四阶奇异摄动边值问题在自适应网格上的一致收敛分析

we study a difference scheme for the fourth-order singular pertur-bation differential equation on the Bakhvalov-Shishkin grid by Green‘‘s function.The method is shown to be uniformly convergent with respect to the perturbation parameter,of order N^-2 in the maxmum norm on Bakhvalov-Shishkin meshes.Numerical results support our theoretical results.     

二维抛物型方程的高精度多重网格解法

提出了数值求解二维抛物型方程的一种新的高精度加权平均紧隐格式，利用Fourier分析方法证明了该格式是无条件稳定的，为了克服传统迭代法在求解隐格式是收敛速度慢的缺陷，利用了多重网格加速技术，大大加快了迭代收敛速度，提高了求解效率，数值实验结果验证了方法的精确性和可靠性。     

连续平均场Ising模型的显式差分格式

将离散的平均场Ｉｓｉｎｇ模型转换为连续的模型，给出了显式差分格式。利用有界延拓法，证明了差分解的收敛性与稳定性。最后通过数值例子说明了此格式的可信性。     

非线性Schr■dinger方程的一种数值模拟方法

该文通过对非线性Schr■dinger方程增加耗散项,提出了一种新的三层线性差分格式.证明了该格式满足连续方程所具有的两个守恒量及收敛性和稳定性.通过数值例子与已知格式进行比较,结果表明该格式计算简单且具有较高精度.     

一种新的求解双曲型守恒律的三阶中心差分格式

本文提出了一种求解双曲型守恒律新的三阶中心差分格式，主要是引入了一种推广的三阶重构，并证明了这种重构在网格边界无振荡．所提的格式保持了中心差分格式简单的优点，不需用Riemann解算器，避免了进行特征解耦．数值试验结果表明本文格式是高精度、高分辨率的。     

平面热传导方程Douglas交替方向隐格式的稳定性与收敛性

1 引言 解抛物方程的交替方向差分方法产生于本世纪五十年代,由于其稳定性好、易编制计算程序,几十年来在科学与工程计算中得到了广泛的应用,最早的交替方向差分格式由Peace- man,D.W.和Rachford,H.H.以及Douglas,J,Jr.和Rachford,H.H.给出,Douglas,J.Jr.于1962年提出了后人称之为Douglas格式的一类交替方向隐格式,其特点是有二阶截断误差而且适用于三维情形.下面简要给出其格式。 考虑平面方形区域上热传导方程 (1.1) (1.2) (1.3)其中Ω=(0,1)×(0,1),为其边界,t为时间变量,以及u_o为已知函数. 令x与y方向网格步长均为h,x_i=ih,y_j=jh,N=1/h为正整数;时间步长为△t,t~n=n△t,记w_(ij)~n=w(x_i,y_j,t~n)以及     

高精度TVD有限差分格式的自由梯度条件

1自由流条件守恒性气动方程组的高精度TVD有限差分格式一般通过使用通量限制器加入非线性的人工粘性来构成。在曲线坐标系中，通量限制器比较的两个因子与各自所在单元的界面雅可比转换行列式相关。为了消除畸变网格单元界面雅可比转换行列式的变化对通量限制器比较结果的影响，本文提出了高精度TVD格式的自由梯度性质，并分别就Osher和HartenTVD格式讨论了同时满足自由流和自由梯度的度量张量计算方法，使高精度TVD格式在畸变网格上得以实现。任意曲线坐标系中守恒型二维欧拉方程表示为其中，E、F为笛卡尔坐标系中的通量；J为坐标变…     

一类抛物-椭圆耦合方程组混合初边值问题的二阶收敛差分格式(Ⅱ)

2.差分格式的可能性和收敛性我们证明[1]中建立的差分格式(6.1-6.17)是唯一可解且二阶收敛的.记(3.1-3.10)的解为{ψ ,p,q},(4.1-4.13)的解为{ψ,p,q;u_1,u_2,v_1,v_2,w_1,w_2}.假设(3)的系数满足如下条件:当|ε_l|≤ε_0,1≤l≤4时,使得