A non‐standard finite difference scheme for an advection‐diffusion‐reaction equation

In this note, a non‐standard finite difference (NSFD) scheme is proposed for an advection‐diffusion‐reaction equation with nonlinear reaction term. We first study the diffusion‐free case of this equation, that is, an advection‐reaction equation. Two exact finite difference schemes are constructed for the advection‐reaction equation by the method of characteristics. As these exact schemes are complicated and are not convenient to use, an NSFD scheme is derived from the exact scheme. Then, the NSFD scheme for the advection‐reaction equation is combined with a finite difference space‐approximation of the diffusion term to provide a NSFD scheme for the advection‐diffusion‐reaction equation. This new scheme could preserve the fixed points, the positivity, and the boundedness of the solution of the original equation. Numerical experiments verify the validity of our analytical results. Copyright © 2014 JohnWiley & Sons, Ltd.     

Application of Exp‐function method to EW‐Burgers equation

In this article, Exp‐function method is used to obtain an exact solution of the equal‐width wave‐Burgers equation (EW‐Burgers). The method is straightforward and concise, and its applications are promising. It is shown that Exp‐function method, with the help of symbolic computation, provides a very effective and powerful mathematical tool for solving EW‐Burgers equation. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

On a boundary value problem for a p‐Dirac equation

The p‐Laplace equation is a nonlinear generalization of the Laplace equation. This generalization is often used as a model problem for special types of nonlinearities. The p‐Laplace equation can be seen as a bridge between very general nonlinear equations and the linear Laplace equation. The aim of this paper is to solve the p‐Laplace equation for 1 < p < 2 and to find strong solutions. The idea is to apply a hypercomplex integral operator and spatial function theoretic methods to transform the p‐Laplace equation into the p‐Dirac equation. This equation will be solved iteratively by using a fixed‐point theorem. Applying operator‐theoretical methods for the p‐Dirac equation and p‐Laplace equation, the existence and uniqueness of solutions in certain Sobolev spaces will be proved. Copyright © 2016 John Wiley & Sons, Ltd.     

Initial conditions for the cylindrical Korteweg‐de Vries equation

In this paper, we find suitable initial conditions for the cylindrical Korteweg‐de Vries equation by first solving exactly the initial‐value problem for localized solutions of the underlying axisymmetric linear long‐wave equation. The far‐field limit of the solution of this linear problem then provides, through matching, an initial condition for the cylindrical Korteweg‐de Vries equation. This initial condition is associated only with the leading wave front of the far‐field limit of the linear solution. The main motivation is to resolve the discrepancy between the exact mass conservation law, and the “mass” conservation law for the cylindrical Korteweg‐de Vries equation. The outcome is that in the linear initial‐value problem all the mass is carried behind the wave front, and then the “mass” in the initial condition for the cylindrical Korteweg‐de Vries equation is zero. Hence, the evolving solution in the cylindrical Korteweg‐de Vries equation has zero “mass.” This situation arises because, unlike the well‐known unidirectional Korteweg‐de Vries equation, the solution of the initial‐value problem for the axisymmetric linear long‐wave problem contains both outgoing and ingoing waves, but in the cylindrical geometry, the latter are reflected at the origin into outgoing waves, and eventually the total outgoing solution is a combination of these and those initially generated.     

Well‐posedness of time‐space fractional stochastic evolution equations driven by α‐stable noise

This paper is devoted to the well‐posedness for time‐space fractional Ginzburg‐Landau equation and time‐space fractional Navier‐Stokes equations by α‐stable noise. The spatial regularity and the temporal regularity of the nonlocal stochastic convolution are firstly established, and then the existence and uniqueness of the global mild solution are obtained by the Banach fixed point theorem and Mittag‐Leffler functions, respectively. Numerical simulations for time‐space fractional Ginzburg‐Landau equation are provided to verify the analysis results.     

Fourth‐order compact scheme for the one‐dimensional sine‐Gordon equation

Finite difference scheme to the generalized one‐dimensional sine‐Gordon equation is considered in this paper. After approximating the second order derivative in the space variable by the compact finite difference, we transform the sine‐Gordon equation into an initial‐value problem of a second‐order ordinary differential equation. Then Padé approximant is used to approximate the time derivatives. The resulting fully discrete nonlinear finite‐difference equation is solved by a predictor‐corrector scheme. Both Dirichlet and Neumann boundary conditions are considered in our proposed algorithm. Stability analysis and error estimate are given for homogeneous Dirichlet boundary value problems using energy method. Numerical results are given to verify the condition for stability and convergence and to examine the accuracy and efficiency of the proposed algorithm. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Efficient local structure‐preserving schemes for the RLW‐type equation

The multisymplectic schemes have been used in numerical simulations for the RLW‐type equation successfully. They well preserve the local geometric property, but not other local conservation laws. In this article, we propose three novel efficient local structure‐preserving schemes for the RLW‐type equation, which preserve the local energy exactly on any time‐space region and can produce richer information of the original problem. The schemes will be mass‐ and energy‐preserving as the equation is imposed on appropriate boundary conditions. Numerical experiments are presented to verify the efficiency and invariant‐preserving property of the schemes. Comparisons with the existing nonconservative schemes are made to show the behavior of the energy affects the behavior of the solution.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1678–1691, 2017     

Complexiton solutions to the Hirota‐Satsuma‐Ito equation

The Hirota bilinear method is a powerful tool for solving nonlinear evolution equations. Together with the linear superposition principle, it can be used to find a special class of explicit solutions that correspond to complex eigenvalues of associated characteristic problems. These solutions are known as complexiton solutions or simply complexitons. In this article, we study complexiton solutions of the the Hirota‐Satsuma‐Ito equation which is a (2 + 1)‐dimensional extension of the Hirota‐Satsuma shallow water wave equation known to describe propagation of unidirectional shallow water waves. We first construct hyperbolic function solutions and consequently derive the so‐called complexitons via the Hirota bilinear method and the linear superposition principle. In particular, we find nonsingular complexiton solutions to the Hirota‐Satsuma‐Ito equation. Finally, we give some illustrative examples and a few concluding remarks.     

A new kind of nonlocal symmetry for the μ‐Camassa‐Holm equation with linear dispersion

The μ‐Camassa‐Holm equation with linear dispersion is a completely integrable model. In this paper, it is shown that this equation has quadratic pseudo‐potentials that allow us to construct pseudo‐potential–type nonlocal symmetries. As an application, we obtain its recursion operator by using this kind of nonlocal symmetry, and we construct a Darboux transformation for the μ‐Camassa‐Holm equation.     

The boundary integral equation approach for numerical solution of the one‐dimensional Sine‐Gordon equation

This article describes a numerical method based on the boundary integral equation and dual reciprocity method for solving the one‐dimensional Sine‐Gordon (SG) equation. The time derivative is approximated by the time‐stepping method and a predictor–corrector scheme is employed to deal with the nonlinearity which appears in the problem. Numerical results are presented for some problems to demonstrate the usefulness and accuracy of this approach. In addition, the conservation of energy in SG equation is investigated. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

D'Alembert‐type solution of the Cauchy problem for the Boussinesq‐Klein‐Gordon equation

In this paper, we construct a weakly‐nonlinear d'Alembert‐type solution of the Cauchy problem for the Boussinesq‐Klein‐Gordon (BKG) equation. Similarly to our earlier work based on the use of spatial Fourier series, we consider the problem in the class of periodic functions on an interval of finite length (including the case of localized solutions on a large interval), and work with the nonlinear partial differential equation with variable coefficients describing the deviation from the oscillating mean value. Unlike our earlier paper, here we develop a novel multiple‐scales procedure involving fast characteristic variables and two slow time scales and averaging with respect to the spatial variable at a constant value of one or another characteristic variable, which allows us to construct an explicit and compact d'Alembert‐type solution of the nonlinear problem in terms of solutions of two Ostrovsky equations emerging at the leading order and describing the right‐ and left‐propagating waves. Validity of the constructed solution in the case when only the first initial condition for the BKG equation may have nonzero mean value follows from our earlier results, and is illustrated numerically for a number of instructive examples, both for periodic solutions on a finite interval, and localized solutions on a large interval. We also outline an extension of the procedure to the general case, when both initial conditions may have nonzero mean values. Importantly, in all cases, the initial conditions for the leading‐order Ostrovsky equations by construction have zero mean, while initial conditions for the BKG equation may have nonzero mean values.     

N solutions for a derivative nonlinear Schrödinger‐type equation via Riemann‐Hilbert approach

The inverse scattering transform for the derivative nonlinear Schrödinger‐type equation is studied via the Riemann‐Hilbert approach. In the direct scattering process, the spectral analysis of the Lax pair is performed, from which a Riemann‐Hilbert problem is established for the derivative nonlinear Schrödinger‐type equation. In the inverse scattering process, N‐soliton solutions of the derivative nonlinear Schrödinger‐type equation are obtained by solving Riemann‐Hilbert problems corresponding to the reflectionless cases. Moreover, the dynamics of the exact solutions are discussed.     

A multigrid compact finite difference method for solving the one‐dimensional nonlinear sine‐Gordon equation

The aim of this paper is to propose a multigrid method to obtain the numerical solution of the one‐dimensional nonlinear sine‐Gordon equation. The finite difference equations at all interior grid points form a large sparse linear system, which needs to be solved efficiently. The solution cost of this sparse linear system usually dominates the total cost of solving the discretized partial differential equation. The proposed method is based on applying a compact finite difference scheme of fourth‐order for discretizing the spatial derivative and the standard second‐order central finite difference method for the time derivative. The proposed method uses the Richardson extrapolation method in time variable. The obtained system has been solved by V‐cycle multigrid (VMG) method, where the VMG method is used for solving the large sparse linear systems. The numerical examples show the efficiency of this algorithm for solving the one‐dimensional sine‐Gordon equation. Copyright © 2014 John Wiley & Sons, Ltd.     

Spatial behaviour of solutions of the dual‐phase‐lag heat equation

In this paper we study the spatial behaviour of solutions of some problems for the dual‐phase‐lag heat equation on a semi‐infinite cylinder. The theory of dual‐phase‐lag heat conduction leads to a hyperbolic partial differential equation with a third derivative with respect to time. First, we investigate the spatial evolution of solutions of an initial boundary‐value problem with zero boundary conditions on the lateral surface of the cylinder. Under a boundedness restriction on the initial data, an energy estimate is obtained. An upper bound for the amplitude term in this estimate in terms of the initial and boundary data is also established. For the case of zero initial conditions, a more explicit estimate is obtained which shows that solutions decay exponentially along certain spatial‐time lines. A class of non‐standard problems is also considered for which the temperature and its first two time derivatives at a fixed time T are assumed proportional to their initial values. Copyright © 2004 John Wiley & Sons, Ltd.     

A Jacobi–Davidson method for two‐real‐parameter nonlinear eigenvalue problems arising from delay‐differential equations

The critical delays of a delay‐differential equation can be computed by solving a nonlinear two‐parameter eigenvalue problem. The solution of this two‐parameter problem can be translated to solving a quadratic eigenvalue problem of squared dimension. We present a structure preserving QR‐type method for solving such quadratic eigenvalue problem that only computes real‐valued critical delays; that is, complex critical delays, which have no physical meaning, are discarded. For large‐scale problems, we propose new correction equations for a Newton‐type or Jacobi–Davidson style method, which also forces real‐valued critical delays. We present three different equations: one real‐valued equation using a direct linear system solver, one complex valued equation using a direct linear system solver, and one Jacobi–Davidson style correction equation that is suitable for an iterative linear system solver. We show numerical examples for large‐scale problems arising from PDEs. Copyright © 2012 John Wiley & Sons, Ltd.     

Asymptotic numerical analysis of the diffusion‐discrete absorption equation

The diffusion‐discrete absorption (DDA) equation is considered. This equation contains the standard diffusion term and the discrete sorption expressed by a sum of large number of δ‐functions with the support at a non‐uniform mesh multiplied by the unknown function (concentration). The main result of the paper is the homogenization (continualization) of this equation when the small parameter is the characteristic step h of the mesh. The error estimates are proved for the difference of the exact solution of the DDA equation and the solution of the homogenized differential equation. Copyright © 2012 John Wiley & Sons, Ltd.     

Symmetry analysis of the nonlinear two‐dimensional Klein–Gordon equation with a time‐varying delay

The group analysis method is applied to the two‐dimensional nonlinear Klein–Gordon equation with time‐varying delay. Determining equations for equations with a time‐varying delay are derived. A complete group classification of the studied equation with respect to the function involved into the equation is obtained. All admitted Lie algebras are classified. By using the classifications, representations of all invariant solutions are found. Copyright © 2017 John Wiley & Sons, Ltd.     

Blow‐up profiles of solutions for the exponential reaction‐diffusion equation

We consider the blow‐up of solutions for a semilinear reaction‐diffusion equation with exponential reaction term. It is known that certain solutions that can be continued beyond the blow‐up time possess a non‐constant self‐similar blow‐up profile. Our aim is to find the final time blow‐up profile for such solutions. The proof is based on general ideas using semigroup estimates. The same approach works also for the power nonlinearity. Copyright © 2011 John Wiley & Sons, Ltd.     

Finite difference discretization of the Benjamin‐Bona‐Mahony‐Burgers equation

Numerical solutions of the Benjamin‐Bona‐Mahony‐Burgers equation in one space dimension are considered using Crank‐Nicolson‐type finite difference method. Existence of solutions is shown by using the Brower's fixed point theorem. The stability and uniqueness of the corresponding methods are proved by the means of the discrete energy method. The convergence in L^∞‐norm of the difference solution is obtained. A conservative difference scheme is presented for the Benjamin‐Bona‐Mahony equation. Some numerical experiments have been conducted in order to validate the theoretical results.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007