A parareal waveform relaxation algorithm for semi-linear parabolic partial differential equations

We report a new parallel iterative algorithm for semi-linear parabolic partial differential equations (PDEs) by combining a kind of waveform relaxation (WR) techniques into the classical parareal algorithm. The parallelism can be simultaneously exploited by WR and parareal in different directions. We provide sharp error estimations for the new algorithm on bounded time domain and on unbounded time domain, respectively. The iterations of the parareal and the WR are balanced to optimize the performance of the algorithm. Furthermore, the speedup and the parallel efficiency of the new approach are analyzed. Numerical experiments are carried out to verify the effectiveness of the theoretic work.     

Locally linearized fractional step methods for nonlinear parabolic problems

This work deals with the efficient numerical solution of a class of nonlinear time-dependent reaction-diffusion equations. Via the method of lines approach, we first perform the spatial discretization of the original problem by applying a mimetic finite difference scheme. The system of ordinary differential equations arising from that process is then integrated in time with a linearly implicit fractional step method. For that purpose, we locally decompose the discrete nonlinear diffusion operator using suitable Taylor expansions and a domain decomposition splitting technique. The totally discrete scheme considers implicit time integrations for the linear terms while explicitly handling the nonlinear ones. As a result, the original problem is reduced to the solution of several linear systems per time step which can be trivially decomposed into a set of uncoupled parallelizable linear subsystems. The convergence of the proposed methods is illustrated by numerical experiments.     

Computing periodic solutions of linear differential-algebraic equations by waveform relaxation

We propose an algorithm, which is based on the waveform relaxation (WR) approach, to compute the periodic solutions of a linear system described by differential-algebraic equations. For this kind of two-point boundary problems, we derive an analytic expression of the spectral set for the periodic WR operator. We show that the periodic WR algorithm is convergent if the supremum value of the spectral radii for a series of matrices derived from the system is less than 1. Numerical examples, where discrete waveforms are computed with a backward-difference formula, further illustrate the correctness of the theoretical work in this paper.

     

New conservative schemes with discrete variational derivatives for nonlinear wave equations

New conservative finite difference schemes for certain classes of nonlinear wave equations are proposed. The key tool there is “discrete variational derivative”, by which discrete conservation property is realized. A similar approach for the target equations was recently proposed by Furihata, but in this paper a different approach is explored, where the target equations are first transformed to the equivalent system representations which are more natural forms to see conservation properties. Applications for the nonlinear Klein–Gordon equation and the so-called “good” Boussinesq equation are presented. Numerical examples reveal the good performance of the new schemes.     

A parareal approach of semi‐linear parabolic equations based on general waveform relaxation

We present a parareal approach of semi‐linear parabolic equations based on general waveform relaxation (WR) at the partial differential equation (PDE) level. An algorithm for initial‐boundary value problem and two algorithms for time‐periodic boundary value problem are constructed. The convergence analysis of three algorithms are provided. The results show that the algorithm for initial‐boundary value problem is superlinearly convergent while both algorithms for the time‐periodic boundary value problem linearly converge to the exact solutions at most. Numerical experiments show that the parareal algorithms based on general WR at the PDE level, compared with the parareal algorithm based on the classical WR at the ordinary differential equations (ODEs) level (the PDEs is discretized into ODEs), require much fewer number of iterations to converge.     

First-Order and Second-Order,Chaos-Free,Finite Difference Schemes for Fisher Equation

A new class of finite difference schemes is constructed for Fisher partial differential equation i.e. the reaction-diffusion equation with stiff source term: $au(1-u)$. These schemes have the properties that they reduce to high fidelity algorithms in the diffusion-free case namely in which the numerical solutions preserve the properties in the exact solutions for arbitrary time step-size and reaction coefficient α＞0 and all nonphysical spurious solutions including bifurcations and chaos that normally appear in the standard discrete models of Fisher partial differential equation will not occur. The implicit schemes so developed obtain the numerical solutions by solving a single linear algebraic system at each step. The boundness and asymptotic behaviour of numerical solutions obtained by all these schemes are given. The approach constructing the above schemes can be extended to reaction-diffusion equations with other stiff source terms.     

Fourier-Laplace analysis of the multigrid waveform relaxation method for hyperbolic equations

The multigrid waveform relaxation (WR) algorithm has been fairly studied and implemented for parabolic equations. It has been found that the performance of the multigrid WR method for a parabolic equation is practically the same as that of multigrid iteration for the associated steady state elliptic equation. However, the properties of the multigrid WR method for hyperbolic problems are relatively unknown. This paper studies the multigrid acceleration to the WR iteration for hyperbolic problems, with a focus on the convergence comparison between the multigrid WR iteration and the multigrid iteration for the corresponding steady state equations. Using a Fourier-Laplace analysis in two case studies, it is found that the multigrid performance on hyperbolic problems no longer shares the close resemblance in convergence factors between the WR iteration for parabolic equations and the iteration for the associated steady state equations.     

Error and extrapolation of a compact LOD method for parabolic differential equations

This paper is concerned with a compact locally one-dimensional (LOD) finite difference method for solving two-dimensional nonhomogeneous parabolic differential equations. An explicit error estimate for the finite difference solution is given in the discrete infinity norm. It is shown that the method has the accuracy of the second-order in time and the fourth-order in space with respect to the discrete infinity norm. A Richardson extrapolation algorithm is developed to make the final computed solution fourth-order accurate in both time and space when the time step equals the spatial mesh size. Numerical results demonstrate the accuracy and the high efficiency of the extrapolation algorithm.     

Finite Difference Method for Reaction-Diffusion Equation with Nonlocal Boundary Conditions

In this paper, we present a numerical approach to a class of nonlinear reaction-diffusion equations with nonlocal Robin type boundary conditions by finite difference methods. A second-order accurate difference scheme is derived by the method of reduction of order. Moreover, we prove that the scheme is uniquely solvable and convergent with the convergence rate of order two in a discrete L2-norm. A simple numerical example is given to illustrate the efficiency of the proposed method.     

The Waveform Relaxation method for systems of differential/algebraic equations

This paper reports efforts towards establishing a parallel numerical algorithm known as Waveform Relaxation (WR) for simulating large systems of differential/algebraic equations. The WR algorithm was established as a relaxation based iterative method for the numerical integration of systems of ODEs over a finite time interval. In the WR approach, the system is broken into subsystems which are solved independently, with each subsystem using the previous iterate waveform as “guesses” about the behavior of the state variables in other subsystems. Waveforms are then exchanged between subsystems, and the subsystems are then resolved repeatedly with this improved information about the other subsystems until convergence is achieved.

In this paper, a WR algorithm is introduced for the simulation of generalized high-index DAE systems. As with ODEs, DAE systems often exhibit a multirate behavior in which the states vary as differing speeds. This can be exploited by partitioning the system into subsystems as in the WR for ODEs. One additional benefit of partitioning the DAE system into subsystems is that some of the resulting subsystems may be of lower index and, therefore, do not suffer from the numerical complications that high-index systems do. These lower index subsystems may therefore be solved by less specialized simulations. This increases the efficiency of the simulation since only a portion of the problem must be solved with specially tailored code. In addition, this paper established solvability requirements and convergence theorems for varying index DAE systems for WR simulation.     

A Linearized Spectral-Galerkin Method for Three-Dimensional Riesz-Like Space Fractional Nonlinear Coupled Reaction-Diffusion Equations

In this paper, we establish a novel fractional model arising in the chemical reaction and develop an efficient spectral method for the three-dimensional Riesz-like space fractional nonlinear coupled reaction-diffusion equations. Based on the backward difference method for time stepping and the Legendre-Galerkin spectral method for space discretization, we construct a fully discrete numerical scheme which leads to a linear algebraic system. Then a direct method based on the matrix diagonalization approach is proposed to solve the linear algebraic system, where the cost of the algorithm is of a small multiple of $N^4$ ($N$ is the polynomial degree in each spatial coordinate) flops for each time level. In addition, the stability and convergence analysis are rigorously established. We obtain the optimal error estimate in space, and the results also show that the fully discrete scheme is unconditionally stable and convergent of order one in time. Furthermore, numerical experiments are presented to confirm the theoretical claims. As the applications of the proposed method, the fractional Gray-Scott model is solved to capture the pattern formation with an analysis of the properties of the fractional powers.     

一类非线性反应-扩散方程有限差分格式的稳定性研究

In the article,the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established.Then the new function space is introduced and the stability problem for the finite difference scheme is discussed by means of variational approximation method in this function space.The approach used is of a simple characteristic in gaining the stability condition of the scheme.     

Traveling waves connecting equilibrium and periodic orbit for reaction-diffusion equations with time delay and nonlocal response

A class of reaction-diffusion equations with time delay and nonlocal response is considered. Assuming that the corresponding reaction equations have heteroclinic orbits connecting an equilibrium point and a periodic solution, we show the existence of traveling wave solutions of large wave speed joining an equilibrium point and a periodic solution for reaction-diffusion equations. Our approach is based on a transformation of the differential equations to integral equations in a Banach space and the rigorous analysis of the property for a corresponding linear operator. Our approach eventually reduces a singular perturbation problem to a regular perturbation problem. The existence of traveling wave solution therefore is obtained by the application of Liapunov-Schmidt method and the Implicit Function Theorem.     

带有多项式非线性项的高维反应扩散方程有限差分格式的长时间行为

本文考虑了一类带有多项式非线性项的高维反应扩散方程．建立了一个全离散的有限差分格式,并证明了差分解的存在唯一性．分析了由差分格式生成的离散系统的动力性质,在对差分解先验估计的基础上得到了离散动力系统的整体吸引子的存在性．最后证明了差分格式的长时间稳定性和收敛性．     

Convergence of approximate attractors for a fully discrete system for reaction-diffusion equations

The reaction-diffusion equations are approximated by a fully discrete system: a Legendre-Galerkin approximation for the space variables and a semi-implicit scheme for the time integration. The stability and the convergence of the fully discrete system are established. It is also shown that, under a restriction on the space dimension and the growth rate of the nonlinear term, the approximate attractors of the discrete finite dimensional dynamical systems converge to the attractor of the original infinite dimensional dynamical systems. An error estimate of optimal order is derived as well without any further regularity assumption.     

Global asymptotic stability of stochastic BAM neural networks with distributed delays and reaction-diffusion terms

This paper is concerned with global asymptotic stability of a class of reaction-diffusion stochastic Bi-directional Associative Memory (BAM) neural networks with discrete and distributed delays. Based on suitable assumptions, we apply the linear matrix inequality (LMI) method to propose some new sufficient stability conditions for reaction-diffusion stochastic BAM neural networks with discrete and distributed delays. The obtained results are easy to check and improve upon the existing stability results. An example is also given to demonstrate the effectiveness of the obtained results.     

非线性分数阶反应扩散方程组的间断时空有限元方法

利用时间间断空间连续的时空有限元方法构造了空间分数阶反应扩散方程组的可以逐时间层求解的全离散格式.在时间离散区间上,采用Radau积分公式,将插值理论与有限元理论相结合,给出了全离散格式解的存在唯一性结果,并证明了所给格式是无条件稳定的,进而详细给出最优阶L~∞(L~2)模误差估计过程.最后用数值算例验证了理论分析的正确性.     

Waveform Relaxation Methods for Lie-Group Equations

In this paper, we derive and analyse waveform relaxation (WR) methods for solving differential equations evolving on a Lie-group. We present both continuous-time and discrete-time WR methods and study their convergence properties. In the discrete-time case, the novel methods are constructed by combining WR methods with Runge-Kutta-Munthe-Kaas (RK-MK) methods. The obtained methods have both advantages of WR methods and RK-MK methods, which simplify the computation by decoupling strategy and preserve the numerical solution of Lie-group equations on a manifold. Three numerical experiments are given to illustrate the feasibility of the new WR methods.     

Singular functional differential equations of neutral type in Banach spaces

The well-posedness of a large class of singular partial differential equations of neutral type is discussed. Here the term singularity means that the difference operator of such equations is nonatomic at zero. This fact offers many difficulties in applying the usual methods of perturbation theory and Laplace transform technique and thus makes the study interesting. Our approach is new and it is based on functional analysis of semigroup of operators in an essential way, and allows us to introduce a new concept of solutions for such equations. Finally, we study the well-posedness of a singular reaction-diffusion equation of neutral type in weighted Lebesgue's spaces.     

Existence and nonexistence of monotone traveling waves for the delayed Fisher equation

In this paper a new approach based on a shooting method in a half line coupled with the technique of upper-lower solution pair is used to study the existence and nonexistence of monotone waves for one form of the delayed Fisher equation that does not have the quasimonotonicity property. A necessary and sufficient condition is provided. This new method can be extended to investigate many other nonlocal and non-monotone delayed reaction-diffusion equations.