高阶抛物型方程的具有高稳定性的显式与半显式差分格式

高阶抛物型方程的具有高稳定性的显式与半显式差分格式曾文平（华侨大学数学系，泉州３６２０１１）１引言１９６０年Ｃａｖｅｂ在文［１］中，讨论了如下的高阶抛物型方程混合问题提出了一类含权因子α（０≤α≤１）的两层差分格式（初边值条件处理同［１］下同，从略）...     

关于色散方程的一类二阶恒稳显格式

1　引　　言对于具有高阶空间导数的发展方程 ,其显格式因结构简单 ,易于计算 ,具有明显的计算优越性 ,但已有的绝大多数显格式的稳定性条件都十分苛刻 (见 [6 ] -[1 5] ) ,远不如一般隐格式 ,使其应用受到限制 .1 994年《计算物理》中关于“色散方程的一类具任意稳定性的显格式”一文 (见 [1 4 ] ) ,把色散方程显格式的稳定性条件提高到了可以任意选择的程度 ,但截断误差仅为 O(τ+h) .本文构造了新一类双参数显式差分格式 ,它是绝对稳定的 ,且其截断误差是 O(τ+h2 ) ,它结构简单 ,易于实现计算 ,利于实际应用 .我们用数值例子验证了理论…     

关于色散方程u_t=au_(xxx)一类显式差分格式的讨论

关于色散方程u_t=au_(xxx)差分格式的讨论,在[1]和[2]中,分别提出了中层为五点和六点的显式差分格式,其稳定区域分别为 0≤r≤0.7016和-0.0625 ≤r≤1.1851.本文针对这一问题,讨论中层为七点的一类差分格式的稳定性.[1]中格式是本文的特例,并且这类格式的最佳稳定区域为0≤r≤2.394,大约是[2]中稳定范围的二倍,[1]中稳定范围的三倍.     

解三抛物型微分方程的一族高精度显式差分格式

对求解三维抛物型微分方程利用待定参数法构造出截断误差为Ｏ的一族高精度的三层显式差分格式，并讨论了其稳定性。     

解三维抛物型方程的一个高精度显式差分格式

提出了求解三维抛物型方程的一个高精度显式差分格式.首先,推导了一个特殊节点处一阶偏导数(■u)/(■/t)的一个差分近似表达式,利用待定系数法构造了一个显式差分格式,通过选取适当的参数使格式的截断误差在空间层上达到了四阶精度和在时间层上达到了三阶精度.然后,利用Fourier分析法证明了当r1/6时,差分格式是稳定的.最后,通过数值试验比较了差分格式的解与精确解的区别,结果说明了差分格式的有效性.     

关于色散方程的具有高稳定性的显式差分格式

本文对色散方程u_t=au_(xxx)构造了显式差分格式J_4,其截断误差和稳定条件分别为O(τ+ h~2)和|r|≤4.0884,稳定性比[1]的结果|r|≤0.7016和[2]的结果|r|≤1.1851有很大改进,而且格式的形式也比[2]的格式简单得多.     

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

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

高阶Schrodinger方程的差分格式

高阶Ｓｃｈｒｏｄｉｎｇｅｒ方程的差分格式曾文平（泉州华侨大学应用数学系泉州３６２０ｉｆ）关键词：显式差分格式；稳定性分析；高阶Ｓｃｈｒｏｄｉｎｇｅｒ方程ＡＭＳ（１９９１）主题分类：６５Ｍ０６，６５Ｍ１２．高阶Ｓｃｈｒｓｄｉｎｇｅｒ方程在量子力学、非线...     

二维抛物型偏微分方程的绝对稳定显格式

提出了一个解二维抛物型偏微分方程初边值问题的绝对稳定分支三层显式差分格式，格式的局部截断误差阶为Ｏ（△ｔ２＋△ｘ２＋△ｙ２）．实算表明，格式的稳定性能与理论分析是一致的．     

一类广义的Kdv方程的显式差分解

本文考察一类具耗散的广义KdV方程组■ (gradφ(■))_x ■_(xxx)-α■_(xx) γ■ =■(x,t,■)的周期初值问题的显式差分格式．利用有界延拓法证明了该差分格式的收敛性与稳定性,并给出了算法和数值例子．     

五维热传导方程的一族两层显格式

构造了五维热传导方程的一族两层显格式,证明了当截断误差阶为O（τ＋h2）时,其稳定性条件为网比r=hτ2≤21,优于同类的其它显格式,当截断误差阶为O（τ2＋h2）时,可以得到一个简洁而实用的二阶精度的两层显格式.     

Weak second order S-ROCK methods for Stratonovich stochastic differential equations

It is well known that the numerical solution of stiff stochastic ordinary differential equations leads to a step size reduction when explicit methods are used. This has led to a plethora of implicit or semi-implicit methods with a wide variety of stability properties. However, for stiff stochastic problems in which the eigenvalues of a drift term lie near the negative real axis, such as those arising from stochastic partial differential equations, explicit methods with extended stability regions can be very effective. In the present paper our aim is to derive explicit Runge–Kutta schemes for non-commutative Stratonovich stochastic differential equations, which are of weak order two and which have large stability regions. This will be achieved by the use of a technique in Chebyshev methods for ordinary differential equations.     

带吸收边界条件的声波方程显式差分格式的稳定性分析

１．引言 在进行无界或半无界区域上各种波动方程的数值求解时,需引进入工边界以限制计算范围,在这些边界上应加相应的人工边界条件．这种边界条件应保证所求得的有界区域上的解很好地逼近原来无界区域上的解．对波动方程来说,就是在边界上人工反射应尽可能地小,使之对区域内部解的影响在允许的误差范围以内.因而它们被称为无反射边界条件或吸收边界条件．这种边界条件还应保证所形成的有界区域上的微分方程定解问题是适定的．这也是各种数值方法稳定的必要条件。 近二十多年来,已发展了声波方程的各种类型的吸收边界条件,其中以Ｃｌ…     

Incremental displacement-correction schemes for incompressible fluid-structure interaction

In this paper we introduce a class of incremental displacement-correction schemes for the explicit coupling of a thin-structure with an incompressible fluid. These methods enforce a specific Robin–Neumann explicit treatment of the interface coupling. We provide a general stability and convergence analysis that covers both the incremental and the non-incremental variants. Their stability properties are independent of the added-mass effect. The superior accuracy of the incremental schemes (with respect to the original non-incremental variant) is highlighted by the error estimates, and then confirmed in a benchmark by numerical experiments.     

Unconditional stability of parallel alternating difference schemes for semilinear parabolic systems

The general alternating schemes with intrinsic parallelism for semilinear parabolic systems are studied. First we prove the a priori estimates in the discrete H¹ space of the difference solution for these schemes. Then the existence of the difference solution for these schemes follows from the fixed point principle. Finally the unconditional stability of the general alternating schemes is proved. The alternating group explicit scheme, the alternating segment explicit–implicit scheme and the alternating segment Crank–Nicolson scheme are the special cases of the general alternating schemes.     

Linearly implicit methods for nonlinear parabolic equations

We construct and analyze combinations of rational implicit and explicit multistep methods for nonlinear parabolic equations. The resulting schemes are linearly implicit and include as particular cases implicit-explicit multistep schemes as well as the combination of implicit Runge-Kutta schemes and extrapolation. An optimal condition for the stability constant is derived under which the schemes are locally stable. We establish optimal order error estimates.

     

Explicit finite difference schemes with extended stability for advection equations

In this study an explicit central difference approximation of the generalized leap-frog type is applied to the one- and two-dimensional advection equations. The stability of the considered numerical schemes is investigated and the scheme with the largest stable time step is found. For the linear and nonlinear advection equations numerical experiments with different schemes from the considered class are performed in order to evaluate the practical stability of the designed schemes.     

Stable non-standard implicit finite difference schemes for non-linear heat transfer in a thin finite rod

In this study, first, three non-standard implicit finite difference schemes are proposed for solving the initial-boundary value problem involving a quartic non-linearity that arises in heat transfer involving conduction with thermal radiation. A thin finite rod exposed to radiating heat across its lateral surface into a medium of constant temperature and convection is ignored. Stability and consistency of the third scheme is proved. Numerical results are compared with non-standard explicit finite difference schemes that show fully stability of our third proposed scheme. Then, three non-standard implicit and three non-standard explicit finite difference schemes are proposed for solving the heat transfer problem with additional convection term. It is shown that in the second case when the model involves conduction, radiation and convection terms, the rod reaches steady state sooner. Numerical results for implicit and explicit schemes are compared and the effect of the convection term is discussed.     

三维热传导方程的一族两层显式格式

提出了一族三维热传导方程的两层显式差分格式,当截断误差阶为O(Δt+(Δx)2)时,稳定性条件为网格比r=Δt/(Δx)2=Δt/(Δy)2=Δt/(Δz)2≤1/2,优于其他显式差分格式。而当截断误差阶为O((Δt)2+(Δx)4)时,稳定性条件为r≤1/6,包含了已有的结果。