计算高维带弱奇异核发展型方程的交替方向隐式欧拉方法

本文主要研究高维带弱奇异核的发展型方程的交替方向隐式(ADI)差分方法.向后欧拉(Euler)方法联立一阶卷积求积公式处理时间方向的离散,有限差分方法处理空间方向的离散,并进一步构造了ADI全离散差分格式.然后将二维问题延伸到三维问题,构造三维空间问题的ADI差分格式.基于离散能量法,详细证明了全离散格式的稳定性和误差分析.随后给出了2个数值算例,数值结果进一步验证了时间方向的收敛阶为一阶,空间方向的收敛阶为二阶,和理论分析结果一致.     

Positive numerical solution for a nonarbitrage liquidity model using nonstandard finite difference schemes

In this article, we construct a numerical method based on a nonstandard finite difference scheme to solve numerically a nonarbitrage liquidity model with observable parameters for derivatives. This nonlinear model considers that the parameters involved are observable from order book data. The proposed numerical method use a exact difference scheme in the linear convection‐reaction term, and the spatial derivative is approximated using a nonstandard finite difference scheme. It is shown that the proposed numerical scheme preserves the positivity as well as stability and consistence. To illustrate the accuracy of the method, the numerical results are compared with those produced by other methods. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 210‐221, 2014     

Cosine expansion-based differential quadrature method for numerical solution of the KdV equation

Numerical solution of the Korteweg–de Vries equation is obtained using space-splitting technique and the differential quadrature method based on cosine expansion (CDQM). The details of the CDQM and its implementation to the KdV equation are given. Three test problems are studied to demonstrate the accuracy and efficiency of the proposed method. Accuracy and efficiency are discussed by computing the numerical conserved laws and L₂, L_∞ error norms.     

A higher-order moment method of the lattice Boltzmann model for the conservation law equation

In this paper, we proposed a higher-order moment method in the lattice Boltzmann model for the conservation law equation. In contrast to the lattice Bhatnagar–Gross–Krook (BGK) model, the higher-order moment method has a wide flexibility to select equilibrium distribution function. This method is based on so-called a series of partial differential equations obtained by using multi-scale technique and Chapman–Enskog expansion. According to Hirt’s heuristic stability theory, the stability of the scheme can be controlled by modulating some special moments to design the third-order dispersion term and the fourth-order dissipation term. As results, the conservation law equation is recovered with higher-order truncation error. The numerical examples show the higher-order moment method can be used to raise the accuracy of the truncation error of the lattice Boltzmann scheme for the conservation law equation.     

Computation of the attractive force of an ellipsoid

A numerical method for computing the attractive force of an ellipsoid is proposed that does not involve separating subdomains with singularities. The sought function is represented as a triple integral such as the inner integral of the kernel can be evaluated analytically with the kernel treated as a weight function. The inner integral is approximated by a quadrature for the product of functions, of which one has an integrable singularity. As a result, the integrand obtained before the second integration has only a weak logarithmic singularity. The subsequent change of variables yields an integrand without singularities. Based on this approach, at each stage of integral evaluation with respect to a single variable, quadrature formulas are derived that do not have singularities at integration nodes and do not take large values at these nodes. For numerical experiments, a rather complicated test function is constructed that is the exact attractive force of an ellipsoid of revolution with an elliptic density distribution.     

A combined mixed finite element ADI scheme for solving Richards' equation with mixed derivatives on irregular grids

We propose and analyze an efficient numerical method for solving semilinear parabolic problems with mixed derivative terms on non-rectangular domains. The spatial semidiscretization process is based on an expanded mixed finite element scheme which, combined with suitable quadrature rules, is converted into a cell-centered finite difference scheme. This choice preserves the asymptotic accuracy and local conservation of mass of the method, while substantially reducing the computational cost of the totally discrete scheme. To obtain it, an alternating direction implicit scheme is used for the integration in time. The resulting numerical algorithm involves sets of uncoupled tridiagonal systems which can be solved in parallel. We set out some theoretical results of unconditional convergence (of second order in space and first order in time) for our method. Finally, a numerical experiment is shown in order to illustrate the theoretical results.     

Mixed Laguerre–Legendre Spectral Method for Incompressible Flow in an Infinite Strip

This paper concerns the mixed Laguerre–Legendre spectral approximation and its application to numerical simulation of incompressible flow in an infinite strip. Some approximation results in weighted Sobolev spaces are given. A Laguerre–Legendre spectral scheme for the stream function form of Navier–Stokes equations is constructed. The stability and the convergence of the proposed scheme are proved. The numerical experiments show the high accuracy of this method. The main techniques used in this paper are also applicable to other nonlinear partial differential equations in an infinite strip.     

Numerical solution of nonlinear two-dimensional Fredholm integral equations of the second kind using Gauss product quadrature rules

The Gauss product quadrature rules and collocation method are applied to reduce the second-kind nonlinear two-dimensional Fredholm integral equations (FIE) to a nonlinear system of equations. The convergence of the proposed numerical method is proved under certain conditions on the kernel of the integral equation. An iterative method for approximating the solution of the obtained nonlinear system is provided and its convergence is proved. Also, some numerical examples are presented to show the efficiency and accuracy of the proposed method.     

Chaotic bursting lag synchronization of Hindmarsh–Rose system via a single controller

Chaotic bursting lag synchronization of Hindmarsh–Rose system is investigated. Two lag synchronization schemes with only a single controller are proposed to synchronize Hindmarsh–Rose chaotic system via back stepping method. Especially in the second scheme, only one state variable is contained in the controller. Based on Lyapunov stability theory, the sufficient conditions for synchronization are obtained analytically in both cases. Finally, numerical simulations are provided to show the effectiveness of the developed methods.     

A fast numerical solution of the 3D nonlinear tempered fractional integrodifferential equation

In this paper, we investigate the numerical solution of the three-dimensional (3D) nonlinear tempered fractional integrodifferential equation which is subject to the initial and boundary conditions. The backward Euler (BE) method in association with the first-order convolution quadrature rule is employed to discretize this equation for time, and the Galerkin finite element method is applied for space, which is combined with an alternating direction implicit (ADI) algorithm, in order to reduce the computational cost for solving the three-dimensional nonlocal problem. Then a fully discrete BE ADI Galerkin finite element scheme can be obtained by linearizing the non-linear term. Thereafter we prove a positive-type lemma, from which the stability and convergence of the proposed numerical scheme are derived based on the energy method. Numerical experiments are performed to verify the effectiveness of the proposed approach.