一类强色散DGH方程的多辛普雷斯曼格式

DGH方程作为一类重要的非线性水波方程有着许多广泛的应用前景.基于Hamilton系统的多辛理论研究了一类强色散DGH方程的数值解法,利用多辛普雷斯曼方法构造了一种典型的半隐式的多辛格式.分析了该格式的局部能量和动量守恒律误差,并给出了数值算例.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.     

New analysis and application of fractional order Schrödinger equation using with Atangana–Batogna numerical scheme

In this work, an analytical approximation to the solution of Schrodinger equation has been provided. The fractional derivative used in this equation is the Caputo derivative. The existence and uniqueness conditions of solutions for the proposed model are derived based on the power law. While solving the fractional order Schrodinger equation, Atangana–Batogna numerical method is presented for fractional order equation. We obtain an efficient recurrence relation for solving these kinds of equations. To illustrate the usefulness of the numerical scheme, the numerical simulations are presented. The results show that the numerical scheme is very effective and simple.     

求解非线性互补问题的微分方程方法

本文构造了一种求解非线性互补问题的微分方程方法.在一定条件下,证明了微分方程系统的平衡点是非线性互补问题的解并且基于一般微分方程系统的数值积分建立了一个数值算法.在适当的条件下,证明了此算法产生的序列解是收敛的.本文最后给出了数值结果,该结果表明了此微分方程方法的有效性.     

Numerical Treatment of Homogeneous Semi-Markov Processes in Transient Case–a Straightforward Approach

This paper presents the numerical solution of the process evolution equation of a homogeneous semi-Markov process (HSMP) with a general quadrature method. Furthermore, results that justify this approach proving that the numerical solution tends to the evolution equation of the continuous time HSMP are given. The results obtained generalize classical results on integral equation numerical solutions applying them to particular kinds of integral equation systems. A method for obtaining the discrete time HSMP is shown by applying a very particular quadrature formula for the discretization. Following that, the problem of obtaining the continuous time HSMP from the discrete one is considered. In addition, the discrete time HSMP in matrix form is presented and the fact that the solution of the evolution equation of this process always exists is proved. Afterwards, an algorithm for solving the discrete time HSMP is given. Finally, a simple application of the HSMP is given for a real data social security example.     

Numerical treatment of the spherically symmetric solutions of a generalized Fisher-Kolmogorov-Petrovsky-Piscounov equation

In the present work, the connection of the generalized Fisher-KPP equation to physical and biological fields is noted. Radially symmetric solutions to the generalized Fisher-KPP equation are considered, and analytical results for the positivity and asymptotic stability of solutions to the corresponding time-independent elliptic differential equation are quoted. An energy analysis of the generalized theory is carried out with further physical applications in mind, and a numerical method that consistently approximates the energy of the system and its rate of change is presented. The method is thoroughly tested against analytical and numerical results on the classical Fisher-KPP equation, the Heaviside equation, and the generalized Fisher-KPP equation with logistic nonlinearity and Heaviside initial profile, obtaining as a result that our method is highly stable and accurate, even in the presence of discontinuities. As an application, we establish numerically that, under the presence of suitable initial conditions, there exists a threshold for the relaxation time with the property that solutions to the problems considered are nonnegative if and only if the relaxation time is below a critical value. An analytical prediction is provided for the Heaviside equation, against which we verify the validity of our computational code, and numerical approximations are provided for several generalized Fisher-KPP problems.     

Sobolev方程的全离散有限体积元格式及数值模拟

本文研究二维Sobolev方程的有限体积元方法, 给出一种全离散化有限体积元格式及其有限体积元解的误差估计,并用数值例子说明数值计算的结果与理论结果是相吻合的, 进一步说明了有限体积元方法比其他数值方法更优越.     

一类非线性延迟微分方程数值解的振动性分析

本文考虑一类非线性延迟微分方程-带有单调造血率的造血模型数值解的振动性.通过研究特征方程根的情况得到数值解振动的条件并且讨论了非振动的数值解的一些性质.为了更有力的说明我们的结果,最后给出了相应的算例.     

Lattice Boltzmann model for two-dimensional generalized sine-Gordon equation

The nonlinear sine-Gordon equation arises in various problems in science and engineering. In this paper, we propose a numerical model based on lattice Boltmann method to obtain the numerical solutions of two-dimensional generalized sine-Gordon equation, including damped and undamped sine-Gordon equation. By choosing properly the conservation condition between the macroscopic quantity $u_t$ and the distribution functions and applying the Chapman-Enskog expansion, the governing equation is recovered correctly from the lattice Boltzmann equation. Moreover, the local equilibrium distribution function is obtained. The numerical results of the first three examples agree well with the analytic solutions, which indicates the lattice Boltzmann model is satisfactory and efficient. Numerical solutions for cases involving the most known from the bibliography line and ring solitons are given. Numerical experiments also show that the present scheme has a good long-time numerical behavior for the generalized sine-Gordon equation. Moreover, the model can also be applied to other two-dimensional nonlinear wave equations, such as nonlinear hyperbolic telegraph equation and Klein-Gordon equation.     

Horizontal water flow in unsaturated porous media using a fractional integral method with an adaptive time step

Predicting the horizontal groundwater flow in unsaturated porous media is a challenge in many areas of science and engineering. The governing equation associated with this phenomenon is a nonlinear partial differential equation known as the Richards equation. However, the numerical results obtained using this equation can differ substantially from the experimental results. In order to overcome this difficulty, a new version of the Richards equation was proposed recently that considers a time derivative of fractional order. In this study, we present a numerical method for solving this fractional Richards equation. Our method comprises an adaptive time marching scheme that uses Picard iterations to solve the corresponding nonlinear equations. A computational code was implemented for the proposed method using the Scilab programming language. We performed numerical simulations of the anomalous diffusion of water in a white siliceous brick and showed that the numerical results were consistent with the available experimental data.     

Implicit numerical approximation scheme for the fractional Fokker-Planck equation

In this paper, we consider an anomalous subdiffusion process, governed by fractional Fokker-Planck equation. An effective numerical method for approximating Fokker-Planck equation in a bounded domain is presented. The stability and convergence of the numerical method are analyzed. Some numerical examples are presented to show the application of the present technique. The numerical results exhibit the good performance of our theoretical analysis.     

A scheme for stable numerical differentiation

A method for stable numerical differentiation of noisy data is proposed. The method requires solving a Volterra integral equation of the second kind. This equation is solved analytically. In the examples considered its solution is computed analytically. Some numerical results of its application are presented. These examples show that the proposed method for stable numerical differentiation is numerically more efficient than some other methods, in particular, than variational regularization.     

Numerical simulation for the variable-order Galilei invariant advection diffusion equation with a nonlinear source term

In this paper, we consider the variable-order Galilei advection diffusion equation with a nonlinear source term. A numerical scheme with first order temporal accuracy and second order spatial accuracy is developed to simulate the equation. The stability and convergence of the numerical scheme are analyzed. Besides, another numerical scheme for improving temporal accuracy is also developed. Finally, some numerical examples are given and the results demonstrate the effectiveness of theoretical analysis.     

A Multi-Symplectic Scheme for RLW Equation

The Hamiltonian and multi-symplectic formulations for RLW equation are considered in this paper. A new twelve-point difference scheme which is equivalent to multi-symplectic Preissmann integrator is derived based on the multi-symplectic formulation of RLW equation. And the numerical experiments on solitary waves are also given. Comparing the numerical results for RLW equation with those for KdV equation, the inelastic behavior of RLW equation is shown.     

带强耗散项的拟线性波方程的数值解

研究了一类拟线性波方程的数值解．构造了带强耗散项的拟线性波方程的三级差分格式,并证明其收敛性,估计了差分解的误差．最后给出数值例子．     

An implicit numerical method of a new time distributed-order and two-sided space-fractional advection-dispersion equation

Distributed-order differential equations have recently been investigated for complex dynamical systems, which have been used to describe some important physical phenomena. In this paper, a new time distributed-order and two-sided space-fractional advection-dispersion equation is considered. Firstly, we transform the time distributed-order fractional equation into a multi-term time-space fractional partial differential equation by applying numerical integration. Then an implicit numerical method is constructed to solve the multi-term fractional equation. The uniqueness, stability and convergence of the implicit numerical method are proved. Some numerical results are presented to demonstrate the effectiveness of the method. The method and techniques can be extended to other time distributed-order and space-fractional partial differential equations.     

Modelling dispersion in laminar and turbulent flows in an open channel based on centre manifolds using 1D-IRBFN method

Centre manifold method is an accurate approach for analytically constructing an advection–diffusion equation (and even more accurate equations involving higher-order derivatives) for the depth-averaged concentration of substances in channels. This paper presents a direct numerical verification of this method with examples of the dispersion in laminar and turbulent flows in an open channel with a smooth bottom. The one-dimensional integrated radial basis function network (1D-IRBFN) method is used as a numerical approach to obtain a numerical solution for the original two-dimensional (2-D) advection–diffusion equation. The 2-D solution is depth-averaged and compared with the solution of the 1-D equation derived using the centre manifolds. The numerical results show that the 2-D and 1-D solutions are in good agreement both for the laminar flow and turbulent flow. The maximum depth-averaged concentrations for the 1-D and 2-D models gradually converge to each other, with their velocities becoming practically equal. The obtained numerical results also demonstrate that the longitudinal diffusion can be neglected compared to the advection.     

Polynomial and rational wave solutions of Kudryashov-Sinelshchikov equation and numerical simulations for its dynamic motions

Polynomial and rational wave solutions of Kudryashov-Sinelshchikov equation and numerical simulations for its dynamic motions are investigated. Conservation flows of the dynamic motion are obtained utilizing multiplier approach. Using the unified method, a collection of exact solitary and soliton solutions of Kudryashov-Sinelshchikov equation is presented. Collocation finite element method based on quintic B-spline functions is implemented to the equation to evidence the accuracy of the proposed method by test problems. Stability analysis of the numerical scheme is studied by employing von Neumann theory. The obtained analytical and numerical results are in good agreement.     

A Crank-Nicolson scheme for the Landau-Lifshitz equation without damping

An accurate and efficient numerical approach, based on a finite difference method with Crank-Nicolson time stepping, is proposed for the Landau-Lifshitz equation without damping. The phenomenological Landau-Lifshitz equation describes the dynamics of ferromagnetism. The Crank-Nicolson method is very popular in the numerical schemes for parabolic equations since it is second-order accurate in time. Although widely used, the method does not always produce accurate results when it is applied to the Landau-Lifshitz equation. The objective of this article is to enumerate the problems and then to propose an accurate and robust numerical solution algorithm. A discrete scheme and a numerical solution algorithm for the Landau-Lifshitz equation are described. A nonlinear multigrid method is used for handling the nonlinearities of the resulting discrete system of equations at each time step. We show numerically that the proposed scheme has a second-order convergence in space and time.     

Numerical solution of a heated subsidence mound problem in a porous medium

Two different numerical models are constructed to solve a two dimensional subsidence mound problem heated along the moving wet/dry interface. One numerical model is based on cartesian coordinates while the other is based on polar coordinates. In both approaches coordinate transformations are used that render the interface stationary. The problem involves a system of three coupled equations; an elliptic equation for a stream function, a parabolic equation for the temperature and a non-linear equation for the boundary location. Good agreement is found between the results of both methods. Graphic results are presented for the decay of a subsidence mound for different values of the various parameters in the model problem.     

一维Allen-Cahn方程紧差分格式的离散最大化原则和能量稳定性研究

本文主要研究相场模拟中的Allen-Cahn模型,考虑一维Allen-Cahn方程紧差分方法的数值逼近.建立具有O(∫τ2+h4)精度的全离散紧差分格式,证明在合理的步长比和时间步长的约束下,其数值解满足离散最大化原则,在此基础上,研究了全离散格式的能量稳定性.最后给出数值算例.