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     

Finite element analysis of parabolic integro-differential equations of Kirchhoff type

The aim of this paper is to study parabolic integro-differential equations of Kirchhoff type. We prove the existence and uniqueness of the solution for this problem via Galerkin method. Semidiscrete formulation for this problem is presented using conforming finite element method. As a consequence of the Ritz–Volterra projection, we derive error estimates for both semidiscrete solution and its time derivative. To find the numerical solution of this class of equations, we develop two different types of numerical schemes, which are based on backward Euler–Galerkin method and Crank–Nicolson–Galerkin method. A priori bounds and convergence estimates in spatial as well as temporal direction of the proposed schemes are established. Finally, we conclude this work by implementing some numerical experiments to confirm our theoretical results.     

Numerical approximation and error control for a thermoelastic contact problem

A numerical method using finite elements for the spatial discretization and the Crank–Nicolson scheme for the time stepping is applied to a partial differential equation problem involving thermoelastic contact. The Crank–Nicolson scheme is interpreted as a low order continuous Galerkin method. By exploiting the variational framework inherent in this approach, an a posteriori error estimate is derived. This estimate gives a bound on the approximation error that depends on computable quantities such as the mesh parameters, time step and numerical solution. In this paper, the a posteriori estimate is used to develop a time step refinement strategy. Several computational examples are included that demonstrate the performance of the method and validity of the estimate.     

A linearized Crank–Nicolson Galerkin FEMs for the nonlinear fractional Ginzburg–Landau equation

In this paper, we study a linearized Crank–Nicolson Galerkin finite element method for solving the nonlinear fractional Ginzburg–Landau equation. The boundedness, existence and uniqueness of the numerical solution are studied in details. Then we prove that the optimal error estimates hold unconditionally, in the sense that no restriction on the size of the time step in terms of the spatial mesh size needs to be assumed. Finally, numerical tests are investigated to support our theoretical analysis.     

The Crank–Nicolson–Galerkin finite element method for a nonlocal parabolic equation with moving boundaries

The aim of this article is to establish the convergence and error bounds for the fully discrete solutions of a class of nonlinear equations of reaction–diffusion nonlocal type with moving boundaries, using a linearized Crank–Nicolson–Galerkin finite element method with polynomial approximations of any degree. A coordinate transformation which fixes the boundaries is used. Some numerical tests to compare our Matlab code with some existing moving finite element methods are investigated. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1515–1533, 2015     

A fully Galerkin method for the damped generalized regularized long‐wave (DGRLW) equation

In this article, a fully discrete Galerkin scheme based on a nonlinear Crank–Nicolson method to approximate the solution of the DGRLW equation is constructed. Some a priori bounds are proved as well as error estimates. Then, a linearized modification scheme by an extrapolation method is discussed. The two schemes are time second order convergence. The last part is devoted to some numerical results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A finite element approach for analysis and computational modelling of coupled reaction diffusion models

In this article, we establish the existence and uniqueness of solutions to the coupled reaction–diffusion models using Banach fixed point theorem. The Galerkin finite element method is used for the approximation of solutions, and an a priori error estimate is derived for such approximations. A scheme is proposed by combining the Crank–Nicolson and the predictor–corrector methods for the time discretization. Some numerical examples are considered to illustrate the accuracy and efficiency of the proposed scheme. It is found that the scheme is second‐order convergent. In addition, nonuniform grids are used in some cases to enhance the accuracy of the scheme.     

Stability and convergence of two‐grid Crank‐Nicolson extrapolation scheme for the time‐dependent natural convection equations

In this paper, we consider the Crank‐Nicolson extrapolation scheme for the 2D/3D unsteady natural convection problem. Our numerical scheme includes the implicit Crank‐Nicolson scheme for linear terms and the recursive linear method for nonlinear terms. Standard Galerkin finite element method is used to approximate the spatial discretization. Stability and optimal error estimates are provided for the numerical solutions. Furthermore, a fully discrete two‐grid Crank‐Nicolson extrapolation scheme is developed, the corresponding stability and convergence results are derived for the approximate solutions. Comparison from aspects of the theoretical results and computational efficiency, the two‐grid Crank‐Nicolson extrapolation scheme has the same order as the one grid method for velocity and temperature in H¹‐norm and for pressure in L²‐norm. However, the two‐grid scheme involves much less work than one grid method. Finally, some numerical examples are provided to verify the established theoretical results and illustrate the performances of the developed numerical schemes.     

Numerical study using explicit multistep Galerkin finite element method for the MRLW equation

In this article, an explicit multistep Galerkin finite element method for the modified regularized long wave equation is studied. The discretization of this equation in space is by linear finite elements, and the time discretization is based on explicit multistep schemes. Stability analysis and error estimates of our numerical scheme are derived. Numerical experiments indicate the validation of the scheme by L₂– and L_∞– error norms and three invariants of motion.4 © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1875–1889, 2015     

Interior penalty discontinuous Galerkin methods with implicit time‐integration techniques for nonlinear parabolic equations

We prove existence and numerical stability of numerical solutions of three fully discrete interior penalty discontinuous Galerkin methods for solving nonlinear parabolic equations. Under some appropriate regularity conditions, we give the l²(H¹) and l^∞(L²) error estimates of the fully discrete symmetric interior penalty discontinuous Galerkin–scheme with the implicit θ ‐schemes in time, which include backward Euler and Crank–Nicolson finite difference approximations. Our estimates are optimal with respect to the mesh size h. The theoretical results are confirmed by some numerical experiments. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Multilevel augmentation methods for solving the Burgers' equation

In this article, multilevel augmentation method (MAM) for solving the Burgers' equation is developed. The Crank–Nicolson–Galerkin scheme of the Burgers' equation results in nonlinear algebraic systems at each time step, the computational cost for solving these nonlinear systems is huge. The MAM allows us to solve the nonlinear system at a fixed initial lower level and then compensate the error by solving a linear system at the higher level. We prove that the method has the same optimal convergence order as the projection method, while reducing the computational complexity greatly. Finally, numerical experiments are presented to confirm the theoretical analysis and illustrate the efficiency of the proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1665–1691, 2015     

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 finite element method for the Sivashinsky equation

The Sivashinsky equation is a nonlinear evolutionary equation of fourth order in space. In this paper we have analyzed a semidiscrete finite element method and completely discrete scheme based on the backward Euler method and Crank–Nicolson–Galerkin scheme. A linearized backward Euler method have been developed and error bounds are derived for an L² projection.     

Adaptive operator splitting finite element method for Allen–Cahn equation

In this paper, a new numerical method is proposed and analyzed for the Allen–Cahn (AC) equation. We divide the AC equation into linear section and nonlinear section based on the idea of operator splitting. For the linear part, it is discretized by using the Crank–Nicolson scheme and solved by finite element method. The nonlinear part is solved accurately. In addition, a posteriori error estimator of AC equation is constructed in adaptive computation based on superconvergent cluster recovery. According to the proposed a posteriori error estimator, we design an adaptive algorithm for the AC equation. Numerical examples are also presented to illustrate the effectiveness of our adaptive procedure.     

Algorithms for numerical solution of the modified equal width wave equation using collocation method

Quintic B-spline collocation algorithms for numerical solution of the modified equal width wave (MEW) equation have been proposed. The algorithms are based on Crank–Nicolson formulation for time integration and quintic B-spline functions for space integration. Quintic B-spline collocation method over the finite intervals is also applied to the time split MEW equation and space split MEW equation. Results for the three algorithms are compared by studying the propagation of the solitary wave, interaction of the solitary waves, wave generation and birth of solitons.     

Crank–Nicolson Finite Element Scheme and Modified Reduced-Order Scheme for Fractional Sobolev Equation

Abstract

In this paper, a Crank–Nicolson finite difference/finite element method is considered to obtain the numerical solution for a time fractional Sobolev equation. Firstly, the classical finite element method is presented. Stability and error estimation for the fully discrete scheme are rigorously established. However, the amount of calculation and computing time are too large due to many degrees of freedom of classical finite element scheme and nonlocality of fractional differential operator. And then the modified reduced-order finite element scheme with low dimensions and sufficiently high accuracy, which is based on proper orthogonal decomposition technique, is provided. Stability and convergence for the reduced-order scheme are also studied. At last, numerical examples show that the results of numerical computation are consistent with previous theoretical conclusions.     

Numerical results for the Klein‐Gordon equation in de Sitter spacetime

We consider the initial value problem for the Klein‐Gordon equation in de Sitter spacetime. We use the central difference scheme on the temporal discretization. We also discretize the spatial variable using the finite element method with implicit and the Crank‐Nicolson schemes for the numerical solution of the initial value problem. In order to show the accuracy for the results of the solutions, we also examine the finite difference methods. We observe that the numerical results obtained by using these methods are compatible.     

A Crank–Nicolson collocation spectral method for the two‐dimensional viscoelastic wave equation

In this paper, we first establish the Crank–Nicolson collocation spectral (CNCS) method for two‐dimensional (2D) viscoelastic wave equation by means of the Chebyshev polynomials. And then, we analyze the existence, uniqueness, stability, and convergence of the CNCS solutions. Finally, we use some numerical experiments to verify the correctness of theoretical analysis. This implies that the CNCS model is very effective for solving the 2D viscoelastic wave equations.     

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     

Energy stability and convergence of the scalar auxiliary variable Fourier-spectral method for the viscous Cahn–Hilliard equation

In this paper, we develop two linear and unconditionally energy stable Fourier-spectral schemes for solving viscous Cahn–Hilliard equation based on the recently scalar auxiliary variable approach. The temporal discretizations are built upon the first-order Euler method and second-order Crank–Nicolson method, respectively. We carry out the energy stability and error analysis rigorously. Various classical numerical experiments are performed to validate the efficiency and accuracy of the proposed schemes.