Solitary wave solutions of the modified equal width wave equation

In this paper we use a linearized numerical scheme based on finite difference method to obtain solitary wave solutions of the one-dimensional modified equal width (MEW) equation. Two test problems including the motion of a single solitary wave and the interaction of two solitary waves are solved to demonstrate the efficiency of the proposed numerical scheme. The obtained results show that the proposed scheme is an accurate and efficient numerical technique in the case of small space and time steps. A stability analysis of the scheme is also investigated.     

Finite difference discretization of the Rosenau-RLW equation

In this paper, numerical solutions of the Rosenau-RLW equation are considered using Crank–Nicolson type finite difference method. Existence of the numerical solutions is derived by Brouwer fixed point theorem. A priori bound and the error estimates as well as conservation of discrete mass and discrete energy for the finite difference solutions are discussed. Unconditional stability, second-order convergence and uniqueness of the scheme are proved using discrete energy method. Some numerical experiments have been conducted in order to verify the theoretical results.     

A compact finite difference method for the solution of the generalized Burgers–Fisher equation

In this article, numerical solutions of the generalized Burgers–Fisher equation are obtained using a compact finite difference method with minimal computational effort. To verify this, a combination of a sixth‐order compact finite difference scheme in space and a low‐storage third‐order total variation diminishing Runge–Kutta scheme in time have been used. The computed results with the use of this technique have been compared with the exact solution to show the accuracy of it. The approximate solutions to the equation have been computed without transforming the equation and without using linearization. Comparisons indicate that there is a very good agreement between the numerical solutions and the exact solutions in terms of accuracy. The present method is seen to be a very good alternative to some existing techniques for realistic problems. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

On the convergence of a conservative numerical scheme for the usual Rosenau-RLW equation

In this paper, we study the initial-boundary value problem of the usual Rosenau-RLW equation by finite difference method. We design a conservative numerical scheme which preserves the original conservative properties for the equation. The scheme is three-level and linear-implicit. The unique solvability of numerical solutions has been shown. Priori estimate and second order convergence of the finite difference approximate solutions are discussed by discrete energy method. Numerical results demonstrate that the scheme is efficient and accurate.     

Numerical solutions of random mean square Fisher-KPP models with advection

This paper deals with the construction of numerical stable solutions of random mean square Fisher-Kolmogorov-Petrosky-Piskunov (Fisher-KPP) models with advection. The construction of the numerical scheme is performed in two stages. Firstly, a semidiscretization technique transforms the original continuous problem into a nonlinear inhomogeneous system of random differential equations. Then, by extending to the random framework, the ideas of the exponential time differencing method, a full vector discretization of the problem addresses to a random vector difference scheme. A sample approach of the random vector difference scheme, the use of properties of Metzler matrices and the logarithmic norm allow the proof of stability of the numerical solutions in the mean square sense. In spite of the computational complexity, the results are illustrated by comparing the results with a test problem where the exact solution is known.     

A spline collocation approach for a generalized wave equation subject to non-local conservation condition

In the present paper, a collocation finite element approach based on cubic splines is presented for the numerical solution of a generalized wave equation subject to non-local conservation condition. The efficiency, accuracy and stability of the method are assessed by applying it to a number of test problems. The results are compared with the existing closed-form solutions; the scheme demonstrates that the numerical outcomes are reliable and quite accurate when contrasted with the analytical solutions and an existing numerical method.     

ID-WAVELETS METHOD FOR HAMMERSTEIN INTEGRAL EQUATIONS

1.IntroductionInthispaperwewillconsiderthenumericalsolutionsofthenon--linearintegralequationsofHammersteintype:wheref,kandgaregivenfunctionandyistheunknown.TherehasbeenmuchinterestinthisproblemsinceHammersteinintegralequations,whichcamefromtheelectromagneticfluiddynamics,yieldsstrongphysicalbackground.Moreover,theFredholmintegralequationsofsecondkindarethespecialcaseoftheHammersteinintegralequations.In[6,p.700]thestandardcollocationmethodisappliedtoobtaintheapproximationsolutionofEq.(1).Int…     

Analysis of new conservative difference scheme for two-dimensional Rosenau-RLW equation

In this article, numerical solution for the Rosenau-RLW equation in 2D is considered and a conservative Crank–Nicolson finite difference scheme is proposed. Existence of the numerical solutions for the difference scheme has been shown by Browder fixed point theorem. A priori bound and uniqueness as well as conservation of discrete mass and discrete energy for the finite difference solutions are discussed. Unconditional stability and a second-order accuracy on both space and time of the difference scheme are proved. Numerical experiments are given to support our theoretical results.     

A spline collocation approach for the numerical solution of a generalized nonlinear Klein-Gordon equation

In this paper, a finite element collocation approach using cubic B-splines is employed for the numerical solution of a generalized form of the nonlinear Klein-Gordon equation. The efficiency of the method is tested on a number of examples that represent special cases of the extended equation including the sine-Gordon equation. The numerical results are compared with existing numerical and analytic solutions and the outcomes confirm that the scheme yields accurate and reliable results even when few nodes are used at the time levels.     

The method of fundamental solutions with eigenfunction expansion method for nonhomogeneous diffusion equation

In this article we describe a numerical method to solve a nonhomogeneous diffusion equation with arbitrary geometry by combining the method of fundamental solutions (MFS), the method of particular solutions (MPS), and the eigenfunction expansion method (EEM). This forms a meshless numerical scheme of the MFS‐MPS‐EEM model to solve nonhomogeneous diffusion equations with time‐independent source terms and boundary conditions for any time and any shape. Nonhomogeneous diffusion equation with complex domain can be separated into a Poisson equation and a homogeneous diffusion equation using this model. The Poisson equation is solved by the MFS‐MPS model, in which the compactly supported radial basis functions are adopted for the MPS. On the other hand, utilizing the EEM the diffusion equation is first translated to a Helmholtz equation, which is then solved by the MFS together with the technique of the singular value decomposition (SVD). Since the present meshless method does not need mesh generation, nodal connectivity, or numerical integration, the computational effort and memory storage required are minimal as compared with other numerical schemes. Test results for two 2D diffusion problems show good comparability with the analytical solutions. The proposed algorithm is then extended to solve a problem with irregular domain and the results compare very well with solutions of a finite element scheme. Therefore, the present scheme has been proved to be very promising as a meshfree numerical method to solve nonhomogeneous diffusion equations with time‐independent source terms of any time frame, and for any arbitrary geometry. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

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.     

二维抛物型方程的一族高精度分支稳定显格式

1 引言在渗流、扩散、热传导等领域中经常会遇到求解二维抛物型方程的初边值问题■(1)其中,φ,f1,f2,f3,f4为已知的光滑函数,a>0为热扩散项系数.对问题(1)的求解,有限差分法是解决此类问题的常用方法,常见的差分格式有古典显式格式与Crank-Nicolson格式[1-2],古典显式格式稳定性条件为r≤1/4,局部截断误差为O (Δt+Δx2).     

A meshless method for numerical solution of a linear hyperbolic equation with variable coefficients in two space dimensions

A meshless method is proposed for the numerical solution of the two space dimensional linear hyperbolic equation subject to appropriate initial and Dirichlet boundary conditions. The new developed scheme uses collocation points and approximates the solution employing thin plate splines radial basis functions. Numerical results are obtained for various cases involving variable, singular and constant coefficients, and are compared with analytical solutions to confirm the good accuracy of the presented scheme. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Multiscale asymptotic expansions and numerical algorithms for the wave equations of second order with rapidly oscillating coefficients

In this paper, we discuss the multiscale analysis and numerical algorithms for the wave equations of second order with rapidly oscillating coefficients. The formal multiscale asymptotic expansions of the solutions for these problems in four specific cases are presented. Higher order corrector methods are constructed and associated explicit convergence rates are obtained in some cases. A multiscale numerical method and a symplectic geometric scheme are introduced. Finally, some numerical results and unsolved problems are presented, and these numerical results support strongly the convergence theorem of this paper.     

Study of a Generalized Nonlinear Euler-Poisson-Darboux System: Numerical and Bessel Based Solutions

In this paper a nonlinear Euler-Poisson-Darboux system is considered. In a first part, we proved the genericity of the hypergeometric functions in the development of exact solutions for such a systemin some special cases leading to Bessel type differential equations. Next, a finite difference scheme in two-dimensional case has been developed. The continuous system is transformed into an algebraic quasi linear discrete one leading to generalized Lyapunov-Sylvester operators. The discrete algebraic system is proved to be uniquely solvable, stable and convergent based on Lyapunov criterion of stability and Lax-Richtmyer equivalence theorem for the convergence. A numerical example has been provided at the end to illustrate the efficiency of the numerical scheme developed in section 3. The present method is thus proved to be more accurate than existing ones and lead to faster algorithms.     

用线方法求解带混合边值条件的抛物方程

研制了分别用显式Euler法、隐式Euler法、Crank-Nicolson格式(梯形方法)求解带第一、第二及混合边值条件的抛物问题的应用软件,通过求解若干抛物问题对该软件作了测试,获得了预期的数值结果,讨论了时间和空间步长的变化对格式计算结果的影响,得到了三种方法的稳定性、收敛精度和计算量.     

High-order linear compact conservative method for the nonlinear Schrödinger equation coupled with the nonlinear Klein–Gordon equation

In this paper, we design a linear-compact conservative numerical scheme which preserves the original conservative properties to solve the Klein–Gordon–Schrödinger equation. The proposed scheme is based on using the finite difference method. The scheme is three-level and linear-implicit. Priori estimate and the convergence of the finite difference approximate solutions are discussed by the discrete energy method. Numerical results demonstrate that the present scheme is conservative, efficient and of high accuracy.     

非定常Burgers方程基于POD方法的全二阶差分格式降阶外推算法

利用特征投影分解方法和奇值分解技术对非定常Burgers 方程的经典全二阶有限差分格式做降阶处理, 给出一种时间和空间变量都是二阶精度的降阶差分格式, 并给出这种降阶全二阶精度差分解的误差估计和外推算法的实现, 最后用数值例子说明数值结果与理论结果是相吻合的.     

Stability preserving non-standard finite difference scheme for a harvesting Leslie–Gower predator–prey model

A non-standard finite difference scheme for a harvesting Leslie–Gower equations is constructed. It is shown that the obtained difference system has the same dynamics as the original continuous system, such as positivity of solutions, equilibria and their local stability properties, irrespective of the size of numerical time step. To illustrate the analytical results, we present some numerical simulations.     

Numerical solution of the nonlinear Klein–Gordon equation using radial basis functions

The nonlinear Klein–Gordon equation is used to model many nonlinear phenomena. In this paper, we propose a numerical scheme to solve the one-dimensional nonlinear Klein–Gordon equation with quadratic and cubic nonlinearity. Our scheme uses the collocation points and approximates the solution using Thin Plate Splines (TPS) radial basis functions (RBF). The implementation of the method is simple as finite difference methods. The results of numerical experiments are presented, and are compared with analytical solutions to confirm the good accuracy of the presented scheme.