Methods for the numerical solution of the Benjamin‐Bona‐Mahony‐Burgers equation

L^∞‐error estimates for finite element for Galerkin solutions for the Benjamin‐Bona‐Mahony‐Burgers (BBMB) equation are considered. A priori bound and the semidiscrete Galerkin scheme are studied using appropriate projections. For fully discrete Galerkin schemes, we consider the backward Euler method and analyze the corresponding error estimates. For a second order accuracy in time, we propose a three‐level backward method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Space‐time spectral collocation method for the one‐dimensional sine‐Gordon equation

In this article, we introduce a new space‐time spectral collocation method for solving the one‐dimensional sine‐Gordon equation. We apply a spectral collocation method for discretizing spatial derivatives, and then use the spectral collocation method for the time integration of the resulting nonlinear second‐order system of ordinary differential equations (ODE). Our formulation has high‐order accurate in both space and time. Optimal a priori error bounds are derived in the L²‐norm for the semidiscrete formulation. Numerical experiments show that our formulation have exponential rates of convergence in both space and time. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 670–690, 2015     

A fifth order semidiscrete mKdV equation

In this paper,aiming to get more insight on the relation between the higher order semidiscrete mKdV equations and higher order mKdV equations,we construct a fifth order semidiscrete mKdV equation from the three known semidiscrete mKdV fluxes.We not only give its Lax pairs,Darboux transformation,explicit solutions and infinitely many conservation laws,but also show that their continuous limits yield the corresponding results for the fifth order mKdV equation.We thus conclude that the fifth order discrete mKdV equation is extremely an useful discrete scheme for the fifth order mKdV equation.     

Analysis and computational method based on quadratic B‐spline FEM for the Rosenau–Burgers equation

L_∞‐error estimates for B‐spline Galerkin finite element solution of the Rosenau–Burgers equation are considered. The semidiscrete B‐spline Galerkin scheme is studied using appropriate projections. For fully discrete B‐spline Galerkin scheme, we consider the Crank–Nicolson method and analyze the corresponding error estimates in time. Numerical experiments are given to demonstrate validity and order of accuracy of the proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 877–895, 2016     

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     

Entropy dissipative one‐leg multistep time approximations of nonlinear diffusive equations

New one‐leg multistep time discretizations of nonlinear evolution equations are investigated. The main features of the scheme are the preservation of the non‐negativity and the entropy dissipation structure of the diffusive equations. The key ideas are to combine Dahlquist's G‐stability theory with entropy dissipation methods and to introduce a nonlinear transformation of variables, which provides a quadratic structure in the equations. It is shown that G‐stability of the one‐leg scheme is sufficient to derive discrete entropy dissipation estimates. The general result is applied to a cross‐diffusion system from population dynamics and a nonlinear fourth‐order quantum diffusion model, for which the existence of semidiscrete weak solutions is proved. Under some assumptions on the operator of the evolution equation, the second‐order convergence of solutions is shown. Moreover, some numerical experiments for the population model are presented, which underline the theoretical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1119–1149, 2015     

Hopf bifurcation in a reaction–diffusion model with Degn–Harrison reaction scheme

In this paper, we consider a chemical reaction–diffusion model with Degn–Harrison reaction scheme under homogeneous Neumann boundary conditions. The existence of Hopf bifurcation to ordinary differential equation (ODE) and partial differential equation (PDE) models are derived, respectively. Furthermore, by using the center manifold theory and the normal form method, we establish the bifurcation direction and stability of periodic solutions. Finally, some numerical simulations are shown to support the analytical results, and to reveal new phenomenon on the Hopf bifurcation.     

An explicit spectral collocation method using nonpolynomial basis functions for the time‐dependent Schrödinger equation

We propose a spectral collocation method for the numerical solution of the time‐dependent Schrödinger equation, where the newly developed nonpolynomial functions in a previous study are used as basis functions. Equipped with the new basis functions, various boundary conditions can be imposed exactly. The preferable semi‐implicit time marching schemes are employed for temporal discretization. Moreover, the new basis functions build in a free parameter λ intrinsically, which can be chosen properly so that the semi‐implicit scheme collapses to an explicit scheme. The method is further applied to linear Schrödinger equation set in unbounded domain. The transparent boundary conditions are constructed for time semidiscrete scheme of the linear Schrödinger equation. We employ spectral collocation method using the new basis functions for the spatial discretization, which allows for the exact imposition of the transparent boundary conditions. Comprehensive numerical tests both in bounded and unbounded domain are performed to demonstrate the attractive features of the proposed method.     

Similar constructing method for solving nonhomogeneous mixed boundary value problem of n‐interval composite ode and its application

This paper studies a simple method—Similar Constructing Method (SCM)—for constructing the exact solutions of the nonhomogeneous mixed boundary value problem for sets of n‐interval composite second‐order ordinary differential equation (ODE) with variable coefficient. Then this paper proves the correctness of the solution obtained by SCM. After that, this paper has done simulation experiment. This section uses the SCM to solve the nonhomogeneous boundary value problem of three‐interval composite Bessel equation. Solutions are presented in graphical form for various parameter values, and the influence of parameters on the solution is analyzed. The example shows that using SCM to solve the class of nonhomogeneous mixed boundary value problems of n‐interval composite second‐order linear ODE is easy, convenient, and effective.     

A reduced‐order extrapolated natural boundary element method based on POD for the 2D hyperbolic equation in unbounded domain

In this article, we primarily focuses to study the order‐reduction for the classical natural boundary element (NBE) method for the two‐dimensional (2D) hyperbolic equation in unbounded domain. To this end, we first build a semi‐discretized format about time for the hyperbolic equation and discuss the existence, stability, and convergence of the time semi‐discretized solutions. We then establish the classical fully discretized NBE format from the time semi‐discretized one and analyze the existence, stability, and convergence of the classical NBE solutions. Next, using proper orthogonal decomposition method, we build a reduced‐order extrapolated NBE (ROENBE) format containing very few unknowns but having adequately high accuracy, and we also discuss the existence, stability, and convergence of the ROENBE solutions. Finally, we use some numerical examples to show that the ROENBE method is far superior to the classical NBE one. It shows that the ROENBE method is reliable and effective for solving the 2D hyperbolic equation with the unbounded domain.     

Explicit Peakon solutions to a family of wave-breaking equations

The singular traveling wave solutions of a general 4-parameter family equation which unifies the Camass-Holm equation, the Degasperis-Procesi equation and the Novikov equation are investigated in this paper. At first, we obtain the explicit peakon solutions for one of its specific case that $a=(p+2)c$, $b=(p+1)c$ and $c=1$, which is referred to a generalized Camassa-Holm-Novikov (CHN) equation, by reducing it to a second-order ordinary differential equation (ODE) and solving its associated first-order integrable ODE. By observing the characteristics of peakon solutions to the CHN equation, we construct the peakon solutions for the general 4-parameter breaking wave equation. It reveals that singularities of the peakon solutions come up only when the solutions attain singular points of the equation, which might be a universal principal for all singular traveling wave solutions for wave breaking equations.     

Une méthode de décomposition en temps avec des schémas dʼintégration réversible pour la résolution de système dʼéquations différentielles ordinaires

We propose a time domain decomposition method that breaks the sequentiality of the integration scheme for systems of ODE. Under the condition of differentiability of the flow, we transform the initial value problem into a well-posed boundary values problem using the symmetrization of the interval of time integration and time-reversible integration scheme. For systems of linear ODE, we explicitly construct the block tridiagonal system satisfied by the solutions at the time sub-intervals extremities. We then propose an iterative algorithm of Schwarz type for updating the interfaces conditions which can extend the method to systems of nonlinear ODE.     

四阶抛物偏微分方程的H1-Galerkin混合元方法及数值模拟

到目前为止, H¹-Galerkin 混合有限元方法研究的问题仅局限于二阶发展方程. 然而对于高阶发展方程, 特别是重要的四阶发展方程问题的研究却没有出现. 本文首次提出四阶发展方程的H¹-Galerkin 混合有限元方法, 为了给出理论分析的需要, 我们考虑四阶抛物型发展方程. 通过引进三个适当的中间辅助变量, 形成四个一阶方程组成的方程组系统, 提出四阶抛物型方程的H¹-Galerkin 混合有限元方法. 得到了一维情形下的半离散和全离散格式的最优收敛阶误差估计和多维情形的半离散格式误差估计, 并采用迭代方法证明了全离散格式的稳定性. 最后, 通过数值例子验证了提出算法的可行性. 在一维情况下我们能够同时得到未知纯量函数、一阶导数、负二阶导数和负三阶导数的最优逼近解, 这一点是以往混合元方法所不能得到的.     

On the numerical approach of the enthalpy method for the Stefan problem

In this article an error bound is derived for a piecewise linear finite element approximation of an enthalpy formulation of the Stefan problem; we have analyzed a semidiscrete Galerkin approximation and completely discrete scheme based on the backward Euler method and a linearized scheme is given and its convergence is also proved. A second‐order error estimates are derived for the Crank‐Nicolson Galerkin method. In the second part, a new class of finite difference schemes is proposed. Our approach is to introduce a new variable and transform the given equation into an equivalent system of equations. Then, we prove that the difference scheme is second order convergent. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

Modified ‐expansion method for finding exact wave solutions of the coupled Klein–Gordon–Schrödinger equation

The coupled Klein–Gordon–Schrödinger equation is reduced to a nonlinear ordinary differential equation (ODE) by using Lie classical symmetries, and various solutions of the nonlinear ODE are obtained by the modified ‐expansion method proposed recently. With the aid of solutions of the nonlinear ODE, more explicit traveling wave solutions of the coupled Klein–Gordon–Schrödinger equation are found out. The traveling wave solutions are expressed by the hyperbolic functions, trigonometric functions, and rational functions. Copyright © 2012 John Wiley & Sons, Ltd.     

Lie symmetry analysis to Fisher's equation with time fractional order

The aim of this letter is to apply the Lie group analysis method to the Fisher''s equation with time fractional order. We considered the symmetry analysis, explicit solutions to the time fractional Fisher''s(TFF) equations with Riemann-Liouville (R-L) derivative. The time fractional Fisher''s is reduced to respective nonlinear ordinary differential equation(ODE) of fractional order. We solve the reduced fractional ODE using an explicit power series method.     

Solving 2D time‐fractional diffusion equations by a pseudospectral method and Mittag‐Leffler function evaluation

Two‐dimensional time‐fractional diffusion equations with given initial condition and homogeneous Dirichlet boundary conditions in a bounded domain are considered. A semidiscrete approximation scheme based on the pseudospectral method to the time‐fractional diffusion equation leads to a system of ordinary fractional differential equations. To preserve the high accuracy of the spectral approximation, an approach based on the evaluation of the Mittag‐Leffler function on matrix arguments is used for the integration along the time variable. Some examples along with numerical experiments illustrate the effectiveness of the proposed approach. Copyright © 2016 John Wiley & Sons, Ltd.     

Convergence analysis of virtual element methods for semilinear parabolic problems on polygonal meshes

In this article, we discuss and analyze new conforming virtual element methods (VEMs) for the approximation of semilinear parabolic problems on convex polygonal meshes in two spatial dimension. The spatial discretization is based on polynomial and suitable nonpolynomial functions, and a Euler backward scheme is employed for time discretization. The discrete formulation of both the proposed schemes—semidiscrete and fully discrete (with time discretization) is discussed in detail, and the unique solvability of the resulted schemes is discussed. A priori error estimates for the proposed schemes (semidiscrete and fully discrete) in H¹‐ and L²‐norms are derived under the assumption that the source term f is Lipschitz continuous. Some numerical experiments are conducted to illustrate the performance of the proposed scheme and to confirm the theoretical convergence rates.     

Fast evaluation and high accuracy finite element approximation for the time fractional subdiffusion equation

In this article, an efficient algorithm for the evaluation of the Caputo fractional derivative and the superconvergence property of fully discrete finite element approximation for the time fractional subdiffusion equation are considered. First, the space semidiscrete finite element approximation scheme for the constant coefficient problem is derived and supercloseness result is proved. The time discretization is based on the L1‐type formula, whereas the space discretization is done using, the fully discrete scheme is developed. Under some regularity assumptions, the superconvergence estimate is proposed and analyzed. Then, extension to the case of variable coefficients is also discussed. To reduce the computational cost, the fast evaluation scheme of the Caputo fractional derivative to solve the fractional diffusion equations is designed. Finally, numerical experiments are presented to support the theoretical results.     

Generalized Lattice Heisenberg Magnet Model and Its Quasideterminant Soliton Solutions

We consider a Darboux transformation of a generalized lattice (or semidiscrete) Heisenberg magnet (GLHM) model. We define a Darboux transformation on solutions of the Lax pair and on solutions of the spin evolution equation of the GLHM model. The solutions are expressed in terms of quasideterminants. We give a general expression for K-soliton solutions in terms of quasideterminants. Finally, we obtain one- and two-soliton solutions of the GLHM model using quasideterminant properties.