An efficient explicit full-discrete scheme for strong approximation of stochastic Allen–Cahn equation

In Becker and Jentzen (2019) and Becker et al. (2017), an explicit temporal semi-discretization scheme and a space–time full-discretization scheme were, respectively, introduced and analyzed for the additive noise-driven stochastic Allen–Cahn type equations, with strong convergence rates recovered. The present work aims to propose a different explicit full-discrete scheme to numerically solve the stochastic Allen–Cahn equation with cubic nonlinearity, perturbed by additive space–time white noise. The approximation is easily implementable, performing the spatial discretization by a spectral Galerkin method and the temporal discretization by a kind of nonlinearity-tamed accelerated exponential integrator scheme. Error bounds in a strong sense are analyzed for both the spatial semi-discretization and the spatio-temporal full discretization, with convergence rates in both space and time explicitly identified. It turns out that the obtained convergence rate of the new scheme is, in the temporal direction, twice as high as existing ones in the literature. Numerical results are finally reported to confirm the previous theoretical findings.     

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

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

Fast and Robust Sixth Order Multigrid Computation for 3D Convection Diffusion Equation

We present a sixth-order explicit compact finite difference scheme to solve the three-dimensional (3D) convection-diffusion equation. We first use a multiscale multigrid method to solve the linear systems arising from a 19-point fourth-order discretization scheme to compute the fourth-order solutions on both a coarse grid and a fine grid. Then an operator-based interpolation scheme combined with an extrapolation technique is used to approximate the sixth-order accurate solution on the fine grid. Since the multigrid method using a standard point relaxation smoother may fail to achieve the optimal grid-independent convergence rate for solving convection-diffusion equations with a high Reynolds number, we implement the plane relaxation smoother in the multigrid solver to achieve better grid independency. Supporting numerical results are presented to demonstrate the efficiency and accuracy of the sixth-order compact (SOC) scheme, compared with the previously published fourth-order compact (FOC) scheme.     

A formally fourth-order accurate compact scheme for 3D Poisson equation in cylindrical coordinates

In this paper, we extend our previous work (M.-C. Lai, A simple compact fourth-order Poisson solver on polar geometry, J. Comput. Phys. 182 (2002) 337–345) to 3D cases. More precisely, we present a spectral/finite difference scheme for Poisson equation in cylindrical coordinates. The scheme relies on the truncated Fourier series expansion, where the partial differential equations of Fourier coefficients are solved by a formally fourth-order accurate compact difference discretization. Here the formal fourth-order accuracy means that the scheme is exactly fourth-order accurate while the poles are excluded and is third-order accurate otherwise. Despite the degradation of one order of accuracy due to the presence of poles, the scheme handles the poles naturally; thus, no pole condition is needed. The resulting linear system is then solved by the Bi-CGSTAB method with the preconditioner arising from the second-order discretization which shows the scalability with the problem size.     

Convergence of the Swartztrauber-Sweet Method for Parabolic Problems on a Disk with Singularities

In this article, we study the convergence analysis for the initial and boundary value problem of parabolic equations on a disk with singular solutions. It is assumed that the exact solution performs singular properties that its derivatives go to infinity at the boundary of the disk. We propose a fully implicit time-stepping numerical scheme. A stretching polynomial-like function with a parameter is used to construct a local grid refinement. Over the nonuniform partition, we combine the Swartztrauber-Sweet scheme and the backward Euler method in spatial and temporal discretization, respectively. We carry out convergence analysis and analyze the effects of the parameter. It is shown that our numerical scheme is of first order accuracy for temporal discretization and of almost second order accuracy for spatial discretization. Numerical experiments are performed to illustrate our analysis results and show that there exists an optimal value for the parameter to obtain a best approximate solution.     

A compact difference scheme combined with extrapolation techniques for solving a class of neutral delay parabolic differential equations

In this work, we present an implicit compact difference scheme for solving a class of neutral delay parabolic differential equations (NDPDEs). The unique solvability and unconditional stability of the scheme are proved. The temporal accuracy of the scheme is improved by using different Richardson extrapolation techniques for linear and nonlinear problems, and fourth-order accuracy in both temporal and spatial dimensions is obtained. Finally, numerical experiments are conducted to verify the accuracy and efficiency of the algorithms.     

A fourth-order compact ADI method for solving two-dimensional unsteady convection–diffusion problems

In this article, an exponential high-order compact (EHOC) alternating direction implicit (ADI) method, in which the Crank–Nicolson scheme is used for the time discretization and an exponential fourth-order compact difference formula for the steady-state 1D convection–diffusion problem is used for the spatial discretization, is presented for the solution of the unsteady 2D convection–diffusion problems. The method is temporally second-order accurate and spatially fourth order accurate, which requires only a regular five-point 2D stencil similar to that in the standard second-order methods. The resulting EHOC ADI scheme in each ADI solution step corresponds to a strictly diagonally dominant tridiagonal matrix equation which can be inverted by simple tridiagonal Gaussian decomposition and may also be solved by application of the one-dimensional tridiagonal Thomas algorithm with a considerable saving in computing time. The unconditionally stable character of the method was verified by means of the discrete Fourier (or von Neumann) analysis. Numerical examples are given to demonstrate the performance of the method proposed and to compare mostly it with the high order ADI method of Karaa and Zhang and the spatial third-order compact scheme of Note and Tan.     

求解带刚性源项标量双曲型守恒律方程的保有界WCNS格式

带刚性源项的双曲守恒律方程是很多物理问题,特别是化学反应流的数学模型.本文考虑带刚性源项的标量双曲型守恒律方程,通过时空分离的方式,发展了一类保有界的WCNS格式.对于空间离散,我们将参数化的通量限制器推广到WCNS框架,使得方程对流项离散后满足极值原理.对于时间离散,我们将半离散的WCNS改写成指数形式,采用三阶修正...     

基于Maxwell模型的线性粘弹性波传导问题的半离散/全离散杂交应力有限元法

1引言设ΩСR²为一个凸的有界开集,边界为ЭΩ;T为一个正常数.我们考虑如下基于Maxwell模型的二维粘弹性固体介质波传导问题。     

Second Order Unconditionally Stable and Convergent Linearized Scheme for a Fluid-Fluid Interaction Model

In this paper, a fully discrete finite element scheme with second-order temporal accuracy is proposed for a fluid-fluid interaction model, which consists of two Navier-Stokes equations coupled by a linear interface condition. The proposed fully discrete scheme is a combination of a mixed finite element approximation for spatial discretization, the second-order backward differentiation formula for temporal discretization, the second-order Gear's extrapolation approach for the interface terms and extrapolated treatments in linearization for the nonlinear terms. Moreover, the unconditional stability is established by rigorous analysis and error estimate for the fully discrete scheme is also derived. Finally, some numerical experiments are carried out to verify the theoretical results and illustrate the accuracy and efficiency of the proposed scheme.     

A fourth-order method of the convection-diffusion equations with Neumann boundary conditions

In this paper, we have developed a fourth-order compact finite difference scheme for solving the convection-diffusion equation with Neumann boundary conditions. Firstly, we apply the compact finite difference scheme of fourth-order to discrete spatial derivatives at the interior points. Then, we present a new compact finite difference scheme for the boundary points, which is also fourth-order accurate. Finally, we use a Padé approximation method for the resulting linear system of ordinary differential equations. The presented scheme has fifth-order accuracy in the time direction and fourth-order accuracy in the space direction. It is shown through analysis that the scheme is unconditionally stable. Numerical results show that the compact finite difference scheme gives an efficient method for solving the convection-diffusion equations with Neumann boundary conditions.     

A high-order compact difference scheme for 2D Laplace and Poisson equations in non-uniform grid systems

In this study, a high-order compact scheme for 2D Laplace and Poisson equations under a non-uniform grid setting is developed. Based on the optimal difference method, a nine-point compact difference scheme is generated. Difference coefficients at each grid point and source term are derived. This is accomplished through the consideration of compatibility between the partial differential equation and its difference discretization. Theoretically, the proposed scheme has third- to fourth-order accuracy; its fourth-order accuracy is achieved under uniform grid settings. Two examples are provided to examine performance of the proposed scheme. Compared with the traditional five-point difference scheme, the proposed scheme can produce more accurate results with faster convergence. Another reference scheme with the same nine-point grid stencil is derived based on the five-point scheme. The two nine-point schemes have the same coefficients for each grid points; however, their coefficients for the source term are different. The overall accuracy level of the solution resulting from the proposed scheme is higher than that of the nine-point reference scheme. It is also indicated that the smoothness of grids has significant effects on accuracy and convergence of the solutions; efforts in optimizing the grid configuration and allocation can improve solution accuracy and efficiency. Consequently, with the proposed method, solution under the non-uniform grid setting with appropriate grid allocation would be more accurate than that under the uniform-grid manipulation, with the same number of grid points.     

Approximation of Invariant Measure for Damped Stochastic Nonlinear Schrödinger Equation via an Ergodic Numerical Scheme

In order to inherit numerically the ergodicity of the damped stochastic nonlinear Schrödinger equation with additive noise, we propose a fully discrete scheme, whose spatial direction is based on spectral Galerkin method and temporal direction is based on a modification of the implicit Euler scheme. We not only prove the unique ergodicity of the numerical solutions of both spatial semi-discretization and full discretization, but also present error estimations on invariant measures, which gives order 2 in spatial direction and order \({\frac 12}\) in temporal direction under certain hypotheses.     

带跳随机波动率模型的高阶ADI分裂格式

针对带跳随机波动率模型满足的偏积分微分方程,提出一种新的高阶交替方向隐式（ADI）有限差分格式,该模型是一个具有混合导数和非常数系数的对流扩散型初边值问题.我们将不同的高阶空间离散与时间步ADI分裂格式相结合,得到了一种空间四阶精度、时间二阶精度的有效方法,并采用Fourier方法分析了高阶ADI格式的稳定性.最后,通过对欧式看跌期权定价模型进行数值实验证实了数值方法的高阶收敛性.     

Numerical results for the Klein‐Gordon equation in de Sitter spacetime

We consider the initial value problem for the Klein‐Gordon equation in de Sitter spacetime. We use the central difference scheme on the temporal discretization. We also discretize the spatial variable using the finite element method with implicit and the Crank‐Nicolson schemes for the numerical solution of the initial value problem. In order to show the accuracy for the results of the solutions, we also examine the finite difference methods. We observe that the numerical results obtained by using these methods are compatible.     

Two-Grid Algorithm of $H^{1}$-Galerkin Mixed Finite Element Methods for Semilinear Parabolic Integro-Differential Equations

In this paper, we present a two-grid discretization scheme for semilinear parabolic integro-differential equations by $H^{1}$-Galerkin mixed finite element methods. We use the lowest order Raviart-Thomas mixed finite elements and continuous linear finite element for spatial discretization, and backward Euler scheme for temporal discretization. Firstly, a priori error estimates and some superclose properties are derived. Secondly, a two-grid scheme is presented and its convergence is discussed. In the proposed two-grid scheme, the solution of the nonlinear system on a fine grid is reduced to the solution of the nonlinear system on a much coarser grid and the solution of two symmetric and positive definite linear algebraic equations on the fine grid and the resulting solution still maintains optimal accuracy. Finally, a numerical experiment is implemented to verify theoretical results of the proposed scheme. The theoretical and numerical results show that the two-grid method achieves the same convergence property as the one-grid method with the choice $h=H^2$.     

A higher-order hybrid spline difference method on adaptive mesh for solving singularly perturbed parabolic reaction–diffusion problems with robin-boundary conditions

We propose a hybrid numerical scheme to discretize a class of singularly perturbed parabolic reaction–diffusion problems with robin-boundary conditions on an equidistributed grid. The hybrid difference scheme is developed by using a modified backward difference scheme in time, a combination of the cubic spline and exponential spline difference scheme in space. The proposed scheme uses a cubic spline difference scheme for the discretization of robin-boundary conditions. For the time discretization of the problem, we use the standard uniform mesh while a layer adapted equidistributed grid is generated for the spatial discretization. By equidistributing a curvature-based monitor function, the spatial adaptive grid is able to capture the presence of parabolic boundary layers without using any prior information about the solution. Parameter uniform error estimates are derived to illustrate an optimal convergence of first-order in time and second-order in space for the proposed discretization. The accuracy of the proposed scheme is confirmed by the numerical experiments that underpin the theoretical analysis.     

A numerical technique based on B-spline for a class of time-fractional diffusion equation

This paper presents an efficient numerical technique for solving a class of time-fractional diffusion equation. The time-fractional derivative is described in the Caputo form. The L1 scheme is used for discretization of Caputo fractional derivative and a collocation approach based on sextic B-spline basis function is employed for discretization of space variable. The unconditional stability of the fully-discrete scheme is analyzed. Two numerical examples are considered to demonstrate the accuracy and applicability of our scheme. The proposed scheme is shown to be sixth order accuracy with respect to space variable and (2 − α)-th order accuracy with respect to time variable, where α is the order of temporal fractional derivative. The numerical results obtained are compared with other existing numerical methods to justify the advantage of present method. The CPU time for the proposed scheme is provided.     

Compact Difference Scheme with High Accuracy for One-Dimensional Unsteady Quasi-Linear Biharmonic Problem of Second Kind: Application to Physical Problems

The present paper uses a new two-level implicit difference formula for the numerical study of one-dimensional unsteady biharmonic equation with appropriate initial and boundary conditions. The proposed difference scheme is second-order accurate in time and third-order accurate in space on non-uniform grid and in case of uniform mesh, it is of order two in time and four in space. The approximate solutions are computed without using any transformation and linearization. The simplicity of the proposed scheme lies in its three-point spatial discretization that yields block tri-diagonal matrix structure without the use of any fictitious nodes for handling the boundary conditions. The proposed scheme is directly applicable to singular problems, which is the main utility of our work. The method is shown to be unconditionally stable for model linear problem for uniform mesh. The efficacy of the proposed approach has been tested on several physical problems, including the complex fourth-order nonlinear equations like Kuramoto–Sivashinsky equation and extended Fisher–Kolmogorov equation, where comparison is done with some earlier work. It is clear from numerical experiments that the obtained results are not only in good agreement with the exact solutions but also competent with the solutions derived in earlier research studies.