Optimal error estimate of a conformal Fourier pseudo‐spectral method for the damped nonlinear Schrödinger equation

In this article, a Fourier pseudospectral method, which preserves the conforal conservation la, is proposed for solving the damped nonlinear Schrödinger equation. Based on the energy method and the semi‐norm equivalence between the Fourier pseudospectral method and the finite difference method, a priori estimate for the new method is established, which shows that the proposed method is unconditionally convergent with order of in the discrete ‐norm, where is the time step and is the number of collocation points used in the spectral method. Some numerical results are addressed to confirm our theoretical analysis.     

Fourth‐order compact and energy conservative scheme for solving nonlinear Klein‐Gordon equation

In this article, a fourth‐order compact and conservative scheme is proposed for solving the nonlinear Klein‐Gordon equation. The equation is discretized using the integral method with variational limit in space and the multidimensional extended Runge‐Kutta‐Nyström (ERKN) method in time. The conservation law of the space semidiscrete energy is proved. The proposed scheme is stable in the discrete maximum norm with respect to the initial value. The optimal convergent rate is obtained at the order of in the discrete ‐norm. Numerical results show that the integral method with variational limit gives an efficient fourth‐order compact scheme and has smaller error, higher convergence order and better energy conservation for solving the nonlinear Klein‐Gordon equation compared with other methods under the same condition. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1283–1304, 2017     

Numerical solutions to a BBM‐Burgers model with a nonlocal viscous term

In this paper, we numerically investigate the BBM‐Burgers equation with a nonlocal viscous term (1) where is the Riemann‐Liouville half derivative. In particular, we implement different numerical schemes to approximate the solution and its asymptotical behavior. Also, we compare our numerical results with those given in 2013, 2014 for similar models.     

Fast solution algorithms for exponentially tempered fractional diffusion equations

In this article, a fast‐iterative method and a fast‐direct method is proposed for solving one‐dimensional and two‐dimensional tempered fractional diffusion equations with constant coefficients. The proposed iterative method is accelerated by circulant preconditioning which is shown to converge superlinearly while the proposed direct method is based on circulant and skew‐circulant representation for Toeplitz matrix inversion. In one‐dimensional case, the operation cost of the proposed methods are both shown to be with memory requirement in each time step, where is the number of spatial nodes. With the alternating direction implicit method, it is proven that the proposed fast solution algorithms can be extended to handle two‐dimensional tempered fractional diffusion equations with operation cost and memory requirement in each time step, where the number of spatial nodes in ‐direction and ‐direction both equal to . Numerical examples are provided to illustrate the effectiveness and efficiency of the proposed methods.     

Numerical simulation and solutions of the two‐component second order KdV evolutionarysystem

In this study, with the aid of Wolfram Mathematica 11, the modified exp ‐expansion function method is used in constructing some new analytical solutions with novel structure such as the trigonometric and hyperbolic function solutions to the well‐known nonlinear evolutionary equation, namely; the two‐component second order KdV evolutionary system. Second, the finite forward difference method is used in analyzing the numerical behavior of this equation. We consider equation (6.5) and (6.6) for the numerical analysis. We examine the stability of the two‐component second order KdV evolutionary system with the finite forward difference method by using the Fourier‐Von Neumann analysis. We check the accuracy of the finite forward difference method with the help of and norm error. We present the comparison between the exact and numerical solutions of the two‐component second order KdV evolutionary system obtained in this article which and support with graphics plot. We observed that the modified exp ‐expansion function method is a powerful approach for finding abundant solutions to various nonlinear models and also finite forward difference method is efficient for examining numerical behavior of different nonlinear models.     

Long‐time stability and asymptotic analysis of the IFE method for the multilayer porous wall model

In this article, we study the long‐time stability and asymptotic behavior of the immersed finite element (IFE) method for the multilayer porous wall model for the drug‐eluting stents. First, with the IFE method for the spatial descretization, and the implicit Euler scheme for the temporal discretization, respectively, we deduce the global stability of fully discrete solution. Then, we investigate the asymptotic behavior of the discrete scheme which reveals that the multilayer porous wall model converges to the corresponding elliptic equation if approaches to a steady‐state in both and norms as . Finally, some numerical experiments are given to verify the theoretical predictions.     

Error analysis of mixed finite element method for Poisson‐Nernst‐Planck system

To improve the convergence rate in L² norm from suboptimal to optimal for both electrostatic potential and ionic concentrations in Poisson‐Nernst‐Planck (PNP) system, we propose the mixed finite element method in this article to discretize the electrostatic potential equation, and still use the standard finite element method to discretize the time‐dependent ionic concentrations equations. Optimal error estimates in norm for the electrostatic potential, and in and norms for the ionic concentrations are attained. As a by‐product, the electric field can also achieve a higher approximation order in contrast with the standard finite element method for PNP system. Numerical experiments are performed to validate the theoretical results.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1924–1948, 2017     

Analysis of expanded mixed finite element methods for the generalized forchheimer flows of slightly compressible fluids

The nonlinear Forchheimer equations are used to describe the dynamics of fluid flows in porous media when Darcy's law is not applicable. In this article, we consider the generalized Forchheimer flows for slightly compressible fluids, and then study the expanded mixed finite element method applied to the initial boundary value problem for the resulting degenerate parabolic equation for pressure. The bounds for the solutions, time derivative, and gradient of solutions are established. Utilizing the monotonicity properties of Forchheimer equation and boundedness of solutions, a priori error estimates for solution are obtained in ‐norm, ‐norm as well as for its gradient in ‐norm for all . Optimal ‐error estimates are shown for solutions under some additional regularity assumptions. Numerical results using the lowest order Raviart–Thomas mixed element confirm the theoretical analysis regarding convergence rates. Published 2015. Numer Methods Partial Differential Eq 32: 60–85, 2016     

A new nonconforming finite element with a conforming finite element approximation for a coupled continuum pipe‐flow/Darcy model in Karst aquifers

Numerical method is considered for a coupled continuum pipe‐flow/Darcy model describing flow in porous media with an embedded conduit pipe. A new nonconforming element is constructed to solve the Darcy equation on porous matrix. The existence and uniqueness of the approximation solution are deduced. Optimal error estimates are obtained in and norms. Some numerical examples show the accuracy and efficiency of the presented method. With the same number of nodal‐points and the same amount of computation, the results using the new nonconforming element are much better than those by both conforming element and Wilson nonconforming element on the same mesh. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 778–798, 2016     

A fast discontinuous finite element discretization for the space‐time fractional diffusion‐wave equation

In this article, we study fast discontinuous Galerkin finite element methods to solve a space‐time fractional diffusion‐wave equation. We introduce a piecewise‐constant discontinuous finite element method for solving this problem and derive optimal error estimates. Importantly, a fast solution technique to accelerate Toeplitz matrix‐vector multiplications which arise from discontinuous Galerkin finite element discretization is developed. This fast solution technique is based on fast Fourier transform and it depends on the special structure of coefficient matrices. In each temporal step, it helps to reduce the computational work from required by the traditional methods to log , where is the size of the coefficient matrices (number of spatial grid points). Moreover, the applicability and accuracy of the method are verified by numerical experiments including both continuous and discontinuous examples to support our theoretical analysis.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2043–2061, 2017     

Finite analytic numerical method for three‐dimensional quasi‐laplace equation with conductivity in tensor form

The finite analytic numerical method for 3D quasi‐Laplace equation with conductivity in full tensor form is constructed in this article. For cubic grid system, the gradient of the potential variable will diverge when tending to the common edge joining the four grids with different conductivities. However, the potential gradient along the tangential direction is of limited value. As a consequence, the 3D quasi‐Laplace equations will behave as a quasi‐2D one. An approximate analytical solution of the 3D quasi‐Laplace equation can be found around the common edge, which is expressed as a combination of a power‐law function and a linear function. With the help of this approximate analytical solution, a 3D finite analytical numerical scheme is then constructed. Numerical examples show that the proposed numerical scheme can provide rather accurate solutions only with or subdivisions. More important, the convergent speed of the numerical scheme is independent of the conductivity heterogeneity. In contrast, when using the traditional numerical schemes, typically such as the MPFA method, the refinement ratio for the grid cell needs to increase dramatically to get an accurate result for the strong heterogeneous case.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1475–1492, 2017     

Numerical solution of one‐dimensional Burgers' equation

In this article, an iterative method for the approximate solution of a class of Burgers' equation is obtained in reproducing kernel space . It is proved the approximation converges uniformly to the exact solution u(x, t) for any initial function under trivial conditions, the derivatives of are also convergent to the derivatives of u(x, t), and the approximate solution is the best approximation under the system © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1251–1264, 2015     

On the Euler implicit/explicit iterative scheme for the stationary Oldroyd fluid

In this article, we consider the stationary Oldroyd fluid equations from the large time behavior research of the nonstationary equations. Thus, to obtain its numerical solution, we first solve the nonstationary Oldroyd fluid equations via the Euler implicit/explicit finite element method with the integral term discretized by the right‐hand rectangle rule, then increase the total time (i.e., number of time steps) to approximate the solution of the original stationary equations. Under a new uniqueness condition (A2), we prove the exponential stability of the solution pair for the stationary equations and the almost unconditional stability of the numerical method. Furthermore, we also obtain the uniform optimal and error estimates in time integral . Finally, several numerical experiments are provided to verify our theoretical results.     

Two unconditionally stable and convergent difference schemes with the extrapolation method for the one‐dimensional distributed‐order differential equations

The Grünwald formula is used to solve the one‐dimensional distributed‐order differential equations. Two difference schemes are derived. It is proved that the schemes are unconditionally stable and convergent with the convergence orders and in maximum norm, respectively, where and are step sizes in time, space and distributed order. The extrapolation method is applied to improve the approximate accuracy to the orders and respectively. An illustrative numerical example is given to confirm the theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 591–615, 2016     

Asymptotic Analysis and Optimal Error estimates for Benjamin‐Bona‐Mahony‐Burgers' Type Equations

In this article, stabilization result for the Benjamin‐Bona‐Mahony‐Burgers' (BBM‐B) equation, that is, convergence of unsteady solution to steady state solution is established under the assumption that a linearized steady state eigenvalue problem has a minimal positive eigenvalue. Based on appropriate conditions on the forcing function, exponential decay estimates in , and ‐norms are derived, which are valid uniformly with respect to the coefficient of dispersion as it tends to zero. It is, further, observed that the decay rate for the BBM‐B equation is smaller than that of the decay rate for the Burgers equation. Then, a semidiscrete Galerkin method for spatial direction keeping time variable continuous is considered and stabilization results are discussed for the semidiscrete problem. Moreover, optimal error estimates in ‐norms preserving exponential decay property are established using the steady state error estimates. For a complete discrete scheme, a backward Euler method is applied for the time discretization and stabilization results are again proved for the fully discrete problem. Subsequently, numerical experiments are conducted, which verify our theoretical results. The article is finally concluded with a brief discussion on an extension to a multidimensional nonlinear Sobolev equation with Burgers' type nonlinearity.     

Error estimates of fully discrete finite element solutions for the 2D Cahn–Hilliard equation with infinite time horizon

In this article, we deal with a rigorous error analysis for the finite element solutions of the two‐dimensional Cahn–Hilliard equation with infinite time. The error estimates with respect to are proven for the fully discrete conforming piecewise linear element solution under Assumption (A1) on the initial value and Assumption (A2) on the discrete spectrum estimate in the finite element space. The analysis is based on sharp a‐priori estimates for the solutions, particularly reflecting their behavior as . Numerical experiments are carried out to support the theoretical analysis and demonstrate the efficiency of the fully discrete mixed finite element methods. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 742–762, 2017     

Unconditional superconvergence analysis of an H1‐galerkin mixed finite element method for nonlinear Sobolev equations

An efficient H¹‐Galerkin mixed finite element method (MFEM) is presented with and zero order Raviart‐Thomas elements for the nonlinear Sobolev equations. On one hand, the existence and uniqueness of the solutions of the semidiscrete approximation scheme are proved and the super close results of order for the original variable u in a broken H¹ norm and the auxiliary variable in norm are deduced without the boundedness of the numerical solution in ‐norm. Conversely, a linearized Crank‐Nicolson fully discrete scheme with the unconditional super close property is also developed through a new approach, while previous literature always require certain time step conditions (see the references below). Finally, a numerical experiment is included to illustrate the feasibility of the proposed method. Here h is the subdivision parameter and τ is the time step.     

A block‐centered finite difference method for fractional Cattaneo equation

In this article, a block‐centered finite difference method for fractional Cattaneo equation is introduced and analyzed. The unconditional stability and the global convergence of the scheme are proved rigorously. Some a priori estimates of discrete norm with optimal order of convergence both for pressure and velocity are established on nonuniform rectangular grids. Moreover, the applicability and accuracy of the scheme are demonstrated by numerical experiments to support our theoretical analysis.     

High‐order implicit staggered‐grid finite differences methods for the acoustic wave equation

Motivated by the idea that staggered‐grid methods give a greater stability and give energy conservation, this article presents a new family of high‐order implicit staggered‐grid finite difference methods with any order of accuracy to approximate partial differential equations involving second‐order derivatives. In particular, we numerically analyze our new methods for the solution of the one‐dimensional acoustic wave equation. The implicit formulation is based on the plane wave theory and the Taylor series expansion and only involves the solution of tridiagonal matrix equations resulting in an attractive method with higher order of accuracy but nearly the same computation cost as those of explicit formulation. The order of accuracy of the proposal staggered formulas are similar to the methods with conventional grids for a ‐point operator: the explicit formula is th‐order and the implicit formula is th‐order; however, the results demonstrate that new staggered methods are superior in terms of stability properties to the classical methods in the context of solving wave equations.     

Analysis of second order and unconditionally stable BDF2‐AB2 method for the Navier‐Stokes equations with nonlinear time relaxation

In this study, we first consider a second order time stepping finite element BDF2‐AB2 method for the Navier‐Stokes equations (NSE). We prove that the method is unconditionally stable and accurate. Second, we consider a nonlinear time relaxation model which consists of adding a term “” to the Navier‐Stokes Equations with the algorithm depends on BDF2‐AB2 method. We prove that this method is unconditionally stable, too. We applied the BDF2‐AB2 method to several numeral experiments including flow around the cylinder. We have also applied BDF2‐AB2 method with nonlinear time relaxation to some problems. It is observed that when the equilibrium errors are high, applying BDF2‐AB2 with nonlinear time relaxation method to the problem yields lower equilibrium errors.