双项时间分数阶慢扩散方程的一类高效差分方法

反常扩散既是一个重要的物理课题,也是工程中普遍涉及的一个现实问题.针对双项时间分数阶慢扩散方程,本文结合古典显式格式和古典隐式格式,提出了显-隐(Explicit-Implicit,E-I)差分方法和隐-显(Implicit-Explicit,I-E)差分方法.分析证明E-I格式解和I-E格式解的存在唯一性,稳定性和收敛性.理论分析和数值试验结果均表明E-I和I-E差分方法无条件稳定,具有空间2阶精度、时间2-α阶精度.在计算精度一致的要求下,E-I和I-E差分方法相较于经典隐式差分方法具有省时性,证实了E-I差分方法和I-E差分方法求解双项时间分数阶慢扩散方程是高效可行的.     

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

构造了求解一维抛物型方程的一族高精度隐式差分格式.首先,推导了抛物型方程解的一阶偏导数在特殊节点处的一个差分近似式,利用该差分近似式和二阶中心差商近似式用待定系数法构造了一族隐式差分格式,通过选取适当的参数使格式具有高阶截断误差;然后,利用Fourier分析法证明了当r大于1/6时,差分格式是稳定的.最后,通过数值试验将差分格式的解与具有同样精度的其它差分格式的解和精确解进行了比较,并比较了差分格式与经典差分格式的计算效率.结果说明了差分格式的有效性.     

时间分数阶慢扩散方程的一类有效差分方法

对时间分数阶慢扩散方程提出一类数值差分方法:显-隐(Explicit-Implicit, E-I)和隐-显(Implicit-Explicit, I-E)差分方法.它是将古典显式格式与古典隐式格式相结合构造出的一类有效差分格式.理论证明了格式解的存在唯一性,用傅里叶方法证明了格式的稳定性和收敛性.数值试验验证了理论分析,表明E-I格式和I-E格式在具有良好的精度且无条件稳定的情况下,计算速度比隐式格式提高了75%.从而用此格式解决分数阶慢扩散方程是可行的.     

解四阶抛物型方程高精度紧致差分格式

针对四阶抛物型方程周期初值问题,提出了一个两层隐式差分格式和一个三层隐式差分格式.它们的局部截断误差分别为O((Δt)2+(Δx)4)和O((Δt)2+(Δt)(Δx)2+(Δx)4),其中Δt,Δx分别为时间步长和空间步长.误差分析和数值实验均表明,本文构造的差分格式比经典的Crank-Nicolson格式和Saul’ev构造的差分格式精度更高.从精度及稳定性方面考虑,本文构造的格式也比文[5]的显式格式要好.     

时间分数阶反应-扩散方程混合差分格式的并行计算方法

分数阶反应-扩散方程有深刻的物理和工程背景,其数值方法的研究具有重要的科学意义和应用价值.文中提出时间分数阶反应-扩散方程混合差分格式的并行计算方法,构造了一类交替分段显-隐格式(alternative segment explicit-implicit,ASE-I)和交替分段隐-显格式(alternative segment implicit-explicit,ASI-E),这类并行差分格式是基于Saul'yev非对称格式与古典显式差分格式和古典隐式差分格式的有效组合.理论分析格式解的存在唯一性,无条件稳定性和收敛性.数值试验验证了理论分析,表明ASE-I格式和ASI-E格式具有理想的计算精度和明显的并行计算性质,证实了这类并行差分方法求解时间分数阶反应-扩散方程是有效的.     

Kuramoto-Sivashinsky方程的线性化L_∞-差分方法

利用有限差分方法研究Kuramoto-Sivashinsky方程初边值问题的数值解.首先,给出了二阶线性化隐式差分格式,该格式在每一时间层均为线性方程组.其次,给出差分格式的守恒性和数值解的有界性.第三,证明差分格式在最大模意义下的收敛性.最后,通过数值算例验证差分格式的收敛阶,并数值模拟方程的混沌解.     

广义sine-Gordon方程高精度隐式紧致差分方法

本文研究一类二维非线性的广义sine-Gordon(简称SG)方程的有限差分格式.首先构造三层时间的紧致交替方向隐式差分格式,并用能量分析法证明格式具有二阶时间精度和四阶空间精度.然后应用改进的Richardson外推算法将时间精度提高到四阶.最后,数值算例证实改进后的算法在空间和时间上均达到四阶精度.     

求解一维定常对流扩散反应方程的混合型紧致差分格式

借助显式紧致格式和隐式紧致格式的思想,基于截断误差余项修正,并结合原方程本身,构造出了一种求解一维定常对流扩散反应方程的高精度混合型紧致差分格式.格式仅用到三个点上的未知函数值及一阶导数值,而一阶导数值利用四阶Pade格式进行计算,格式整体具有四阶精度.数值实验结果验证了格式的精确性和可靠性.     

解弥散问题的高精度差分格式

用待定系数法 ,对弥散方程构造了一个二层六阶精度的差分格式 ,给出了稳定条件 .用该格式可以直接从初始条件出发逐层求解 ,也可以在使用三层差分格式时 ,用来求第一层的数值解 u1j.     

解对流方程的加耗散项的差分格式

解对流方程的大多数常见的显式差分格式 ,其稳定性条件是苛刻的 .这一困难可由在常规的显式差分格式中引入耗散项而得到克服 .基于此 ,我们导出一类新的无条件稳定的两层的半显式差分格式及若干具有高稳定性的显式格式 .它们包含了若干已知的具有高稳定性的显式格式 .     

Numeric solution of advection–diffusion equations by a discrete time random walk scheme

Explicit numerical finite difference schemes for partial differential equations are well known to be easy to implement but they are particularly problematic for solving equations whose solutions admit shocks, blowups, and discontinuities. Here we present an explicit numerical scheme for solving nonlinear advection–diffusion equations admitting shock solutions that is both easy to implement and stable. The numerical scheme is obtained by considering the continuum limit of a discrete time and space stochastic process for nonlinear advection–diffusion. The stochastic process is well posed and this guarantees the stability of the scheme. Several examples are provided to highlight the importance of the formulation of the stochastic process in obtaining a stable and accurate numerical scheme.     

An inverse problem of finding a source parameter in a semilinear parabolic equation

An inverse problem concerning diffusion equation with source control parameter is considered. Several finite-difference schemes are presented for identifying the control parameter. These schemes are based on the classical forward time centred space (FTCS) explicit formula, and the 5-point FTCS explicit method and the classical backward time centred space (BTCS) implicit scheme, and the Crank–Nicolson implicit method. The classical FTCS explicit formula and the 5-point FTCS explicit technique are economical to use, are second-order accurate, but have bounded range of stability. The classical BTCS implicit scheme and the Crank–Nicolson implicit method are unconditionally stable, but these schemes use more central processor (CPU) times than the explicit finite difference mehods. The basis of analysis of the finite difference equations considered here is the modified equivalent partial differential equation approach, developed from the 1974 work of Warming and Hyett. This allows direct and simple comparison of the errors associated with the equations as well as providing a means to develop more accurate finite difference schemes. The results of a numerical experiment are presented, and the accuracy and CPU time needed for this inverse problem are discussed.     

Comparison of higher-order accurate schemes for solving the two-dimensional unsteady Burgers' equation

This paper is devoted to the testing and comparison of numerical solutions obtained from higher-order accurate finite difference schemes for the two-dimensional Burgers' equation having moderate to severe internal gradients. The fourth-order accurate two-point compact scheme, and the fourth-order accurate Du Fort Frankel scheme are derived. The numerical stability and convergence are presented. The cases of shock waves of severe gradient are solved and checked with the fourth-order accurate Du Fort Frankel scheme solutions. The present study shows that the fourth-order two-point compact scheme is highly stable and efficient in comparison with the fourth-order accurate Du Fort Frankel scheme.     

High-order accurate implicit running schemes

An approach to the construction of high-order accurate implicit predictor-corrector schemes is proposed. The accuracy is improved by choosing a special time integration step for computing numerical fluxes through cell interfaces by using an unconditionally stable implicit scheme. For smooth solutions of advection equations with constant coefficients, the scheme is second-order accurate. Implicit difference schemes for multidimensional advection equations are constructed on the basis of Godunov’s method with splitting over spatial variables as applied to the computation of “large” values at an intermediate layer. The numerical solutions obtained for advection equations and the radiative transfer equations in a vacuum are compared with their exact solutions. The comparison results confirm that the approach is efficient and that the accuracy of the implicit predictor-corrector schemes is improved.     

Two difference schemes for the numerical solution of Maxwell’s equations as applied to extremely and super low frequency signal propagation in the Earth-ionosphere waveguide

Two explicit two-time-level difference schemes for the numerical solution of Maxwell’s equations are proposed to simulate propagation of small-amplitude extremely and super low frequency electromagnetic signals (200 Hz and lower) in the Earth-ionosphere waveguide with allowance for the tensor conductivity of the ionosphere. Both schemes rely on a new approach to time approximation, specifically, on Maxwell’s equations represented in integral form with respect to time. The spatial derivatives in both schemes are approximated to fourth-order accuracy. The first scheme uses field equations and is second-order accurate in time. The second scheme uses potential equations and is fourth-order accurate in time. Comparative test computations show that the schemes have a number of important advantages over those based on finite-difference approximations of time derivatives. Additionally, the potential scheme is shown to possess better properties than the field scheme.     

Spectral problems with mixed Dirichlet–Neumann boundary conditions: Isospectrality and beyond

In this study, an implicit semi-discrete higher order compact (HOC) scheme, with an averaged time discretization, has been presented for the numerical solution of unsteady two-dimensional (2D) Schrödinger equation. The scheme is second order accurate in time and fourth order accurate in space. The results of numerical experiments are presented, and are compared with analytical solutions and well established numerical results of some other finite difference schemes. In all cases, the present scheme produces highly accurate results with much better computational efficiency.     

任意精度的三点紧致显格式及其在CFD中的应用

通过在泰勒级数展开中运用逐阶迭代的方法,推导出了空间任意精度的三点紧致显格式的表达式,又由Fourier分析法得到了格式的数值弥散和耗散特性．与以往的高精度紧致差分格式不同,提出的格式不用隐式求解代数方程组并且可以达到任意精度．通过方波问题和顶盖方腔流的算例表明,格式在稀疏网格下可以得到很高的精度,不仅能节省计算量,而且易于编程,有很高的计算效率．     

二维热传导方程的三层显式差分格式

对二维热传导方程构造了一个稳定的三层显式差分格式求其数值解,其背景源于高维热力学反问题迭代算法中对正问题小计算量算法的需求。首先建立一个含参数的一般差分格式去逼近微分方程,并得到了最优截断误差。然后导出了参数应满足的条件以保证差分格式的稳定性。最后给出了数值的例子并和其它算法进行比较,说明了格式在精度上的有效性和计算量上的优越性。     

求解对流扩散方程的一致四阶紧致格式

目前,许多高精度差分格式,由于未成功地构造与其精度匹配的稳定的边界格式,不得不采用低精度的边界格式.本文针对对流扩散方程证明了存在一致四阶紧致格式,它的边界点的计算格式和内点的计算格式的截断误差主项保持一致,给出了具体内点和边界格式;并分析了此半离散格式的渐近稳定性.数值结果表明该格式是四阶精度;在对流占优情况下,本文边界格式的数值结果比四阶精度的显式差分格式的的数值结果的数值振荡小,取得了不错的效果,理论结果得到了数值验证;驱动方腔数值结果显示,本文对N-S方程的离散格式具有很好的可靠性,适合对复杂流体流动的数值模拟和研究.     

On exact dimensional splitting for a multidimensional scalar quasilinear hyperbolic conservation law

A dimensional splitting scheme is applied to a multidimensional scalar homogeneous quasilinear hyperbolic equation (conservation law). It is proved that the splitting error is zero. The proof is presented for the above partial differential equation in an arbitrary number of dimensions. A numerical example is given that illustrates the proved accuracy of the splitting scheme. In the example, the grid convergence of split (locally one-dimensional) compact and bicompact difference schemes and unsplit bicompact schemes combined with high-order accurate time-stepping schemes (namely, Runge–Kutta methods of order 3, 4, and 5) is analyzed. The errors of the numerical solutions produced by these schemes are compared. It is shown that the orders of convergence of the split schemes remain high, which agrees with the conclusion that the splitting error is zero.