Global error estimation and control in linearly-implicit parallel two-step peer W-methods

The class of linearly-implicit parallel two-step peer W-methods has been designed recently for efficient numerical solutions of stiff ordinary differential equations. Those schemes allow for parallelism across the method, that is an important feature for implementation on modern computational devices. Most importantly, all stage values of those methods possess the same properties in terms of stability and accuracy of numerical integration. This property results in the fact that no order reduction occurs when they are applied to very stiff problems. In this paper, we develop parallel local and global error estimation schemes that allow the numerical solution to be computed for a user-supplied accuracy requirement in automatic mode. An algorithm of such global error control and other technical particulars are also discussed here. Numerical examples confirm efficiency of the presented error estimation and stepsize control algorithm on a number of test problems with known exact solutions, including nonstiff, stiff, very stiff and large-scale differential equations. A comparison with the well-known stiff solver RODAS is also shown.     

Global error estimation and control in linearly-implicit parallel two-step peer W-methods

The class of linearly-implicit parallel two-step peer W-methods has been designed recently for efficient numerical solutions of stiff ordinary differential equations. Those schemes allow for parallelism across the method, that is an important feature for implementation on modern computational devices. Most importantly, all stage values of those methods possess the same properties in terms of stability and accuracy of numerical integration. This property results in the fact that no order reduction occurs when they are applied to very stiff problems. In this paper, we develop parallel local and global error estimation schemes that allow the numerical solution to be computed for a user-supplied accuracy requirement in automatic mode. An algorithm of such global error control and other technical particulars are also discussed here. Numerical examples confirm efficiency of the presented error estimation and stepsize control algorithm on a number of test problems with known exact solutions, including nonstiff, stiff, very stiff and large-scale differential equations. A comparison with the well-known stiff solver RODAS is also shown.     

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

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

A best finite-difference scheme for the fisher equation

A new class of finite-difference schemes is constructed for the Fisher partial differential equation. These schemes are constructed according to the nonstandard modeling rules formulated by Mickens. They have the property that, in the appropriate limits, the discrete models obtained are either “exact” or “best” finite-difference schemes for corresponding differential equation. Consequently, the elementary numerical instabilities will not occur. © 1994 John Wiley & Sons, Inc.     

用格子Boltzmann方法模拟非线性热传导方程

对具有非线性源项和非线性扩散项的热传导方程建立格子Boltzmann求解模型.在演化方程中增加了两个关于源项分布函数的微分算子,对演化方程实施Chapman-Enskog展开.通过对演化方程的进一步改进,恢复出具有高阶截断误差的宏观方程.对不同参数选取下的非线性热传导方程进行了数值模拟,数值解与精确解吻合得很好.该模型也可以用于同类型的其他偏微分方程的数值计算中.     

Numerical studies of a nonlinear heat equation with square root reaction term

Interest in calculating numerical solutions of a highly nonlinear parabolic partial differential equation with fractional power diffusion and dissipative terms motivated our investigation of a heat equation having a square root nonlinear reaction term. The original equation occurs in the study of plasma behavior in fusion physics. We begin by examining the numerical behavior of the ordinary differential equation obtained by dropping the diffusion term. The results from this simpler case are then used to construct nonstandard finite difference schemes for the partial differential equation. A variety of numerical results are obtained and analyzed, along with a comparison to the numerics of both standard and several nonstandard schemes. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Qualitative behavior of numerical traveling solutions for reaction–diffusion equations with memory

In this article the qualitative properties of numerical traveling wave solutions for integro- differential equations, which generalize the well known Fisher equation are studied. The integro-differential equation is replaced by an equivalent hyperbolic equation which allows us to characterize the numerical velocity of traveling wave solutions. Numerical results are presented.     

Stability Analysis and Application of the Exponential Time Differencing Schemes

Exponential time differencing schemes are time integration methods that can be efficiently combined with spatial spectral approximations to provide very high resolution to the smooth solutions of some linear and nonlinear partial differential equations. We study in this paper the stability properties of some exponential time differencing schemes. We also present their application to the numerical solution of the scalar Allen-Cahn equation in two and three dimensional spaces.     

Numerical simulation for a nonlinear partial differential equation with variable coefficients by means of the discrete variational derivative method

Partial differential equations with possibly discontinuous coefficients play an important part in engineering, physics and ecology. In this paper, we will study nonlinear partial differential equations with variable coefficients arising from population models. Generally speaking, it is difficult to analyze the behavior of nonlinear partial differential equations; therefore, we usually rely on the numerical approximation. Currently, there is an increasing interest in designing numerical schemes that preserve energy properties for differential equations. We will design the numerical schemes that preserve discrete energy property and show numerical experiments for a nonlinear partial differential equation with variable coefficients.     

Symmetry preserving numerical schemes for partial differential equations and their numerical tests

The method of equivariant moving frames is used to construct symmetry preserving finite difference schemes of partial differential equations invariant under finite-dimensional symmetry groups. Invariant numerical schemes for a heat equation with logarithmic source and the spherical Burgers' equation are obtained. Numerical tests show how invariant schemes can be more accurate than standard discretizations.     

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.     

Finite Element Simulations with Adaptively Moving Mesh for the Reaction Diffusion System

A moving mesh method is proposed for solving reaction-diffusion equations. The finite element method is used to solving the partial different equation system, and an efficient scheme is applied to implement mesh moving. In the practical calculations, the moving mesh step and the problem equation solver are performed alternatively. Serveral numerical examples are presented, including the Gray-Scott, the Activator-Inhibitor and a case with a growing domain. It is illustrated numerically that the moving mesh method costs much lower, compared with the numerical schemes on a fixed mesh. Even in the case of complex pattern dynamics described by the reaction-diffusion systems, the adapted meshes can capture the details successfully.     

Numerical solution of nonlinear Jaulent–Miodek and Whitham–Broer–Kaup equations

In this work we investigate the numerical solution of Jaulent–Miodek (JM) and Whitham–Broer–Kaup (WBK) equations. The proposed numerical schemes are based on the fourth-order time-stepping schemes in combination with discrete Fourier transform. We discretize the original partial differential equations (PDEs) with discrete Fourier transform in space and obtain a system of ordinary differential equations (ODEs) in Fourier space which will be solved with fourth order time-stepping methods. After transforming the equations to a system of ODEs, the linear operator in JM equation is diagonal but in WBK equation is not diagonal. However for WBK equation we can also implement the methods such as diagonal case which reduces the CPU time. Comparing numerical solutions with analytical solutions demonstrates that those methods are accurate and readily implemented.     

Adaptive decomposition finite difference methods for solving singular problems—A review

Decomposition, or splitting, finite difference methods have been playing an important role in the numerical solution of nonsingular differential equation problems due to their remarkable efficiency, simplicity, and flexibility in computations as compared with their peers. Although the numerical strategy is still in its infancy for solving singular differential equation problems arising from many applications, explorations of the next generation decomposition schemes associated with various kinds of adaptations can be found in many recent publications. The novel approaches have been proven to be highly effective and reliable in operations. In this article, we will focus on some of the latest developments in the area. Key comments and discussion will be devoted to two particularly interesting issues in the research, that is, direct solutions of degenerate singular reaction-diffusion equations and nonlinear sine-Gordon wave equations. Numerical experiments with simulated demonstrations will be given.     

On the asymptotic positivity of solutions for the extended Fisher–Kolmogorov equation with nonlinear diffusion

The objective of this paper aims to prove positivity of solutions for a semilinear dissipative partial differential equation with non‐linear diffusion. The equation is a generalized model of the well‐known Fisher–Kolmogorov equation and represents a class of dissipative partial differential equations containing differential operators of higher order than the Laplacian. It arises in a variety of meaningful physical situations including gas flows, diffusion of an electron–ion plasma and the dynamics of biological populations whose mobility is density dependent. In all these situations, the solutions of the equation must be positive functions. Copyright © 2002 John Wiley & Sons, Ltd.     

On discreteness of the Hopf equation

The principle aim of this essay is to illustrate how different phenomena is captured by different discretizations of the Hopf equation and general hyperbolic conservation laws. This includes dispersive schemes, shock capturing schemes as well as schemes for computing multi-valued solutions of the underlying equation. We introduce some model equations which describe the behavior of the discrete equation more accurate than the original equation. These model equations can either be conveniently discretized for producing novel numerical schemes or further analyzed to enrich the theory of nonlinear partial differential equations.     

Implicit Difference Methods for Degenerate Hyperbolic Equations of Second Order

This paper is a sequel to [2]. A two parameter family of explicit and implicit schemes is constructed for the numerical solution of the degenerate hyperbolic equations of second order. We prove the existence and the uniqueness of the solutions of these schemes. Furthermore, we prove that these schemes are stable for the initial values and that the numerical solution is convergent to the unique generalized solution of the partial differential equation.     

The Pathwise Numerical Approximation of Stationary Solutions of Semilinear Stochastic Evolution Equations

Under a one-sided dissipative Lipschitz condition on its drift, a stochastic evolution equation with additive noise of the reaction-diffusion type is shown to have a unique stochastic stationary solution which pathwise attracts all other solutions. A similar situation holds for each Galerkin approximation and each implicit Euler scheme applied to these Galerkin approximations. Moreover, the stationary solution of the Euler scheme converges pathwise to that of the Galerkin system as the stepsize tends to zero and the stationary solutions of the Galerkin systems converge pathwise to that of the evolution equation as the dimension increases. The analysis is carried out on random partial and ordinary differential equations obtained from their stochastic counterparts by subtraction of appropriate Ornstein-Uhlenbeck stationary solutions.     

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.     

一类时间分数阶偏微分方程的解

考虑一类时间分数阶偏微分方程,该方程包含几种特殊情况:时间分数阶扩散方程、时间分数阶反应-扩散方程、时间分数阶对流-扩散方程以及它们各自相对应的整数阶偏微分方程． 通过Laplace-Fourier变换及其逆变换,该方程在空间全平面和半平面内的基本解可以求出,但其表达式则是通过适当的变形来求．另外,对于有限域上的初边值问题,则可由Sine(Cosine)-Laplace变换导出该方程的一种级数形式的解,并通过两个数值例子来说明该方法的有效性．