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     

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     

Dissipativity of the linearly implicit Euler scheme for Navier‐Stokes equations with delay

In this article, we study the dissipativity of the linearly implicit Euler scheme for the 2D Navier‐Stokes equations with time delay volume forces (NSD). This scheme can be viewed as an application of the implicit Euler scheme to linearized NSD. Therefore, only a linear system is needed to solve at each time step. The main results we obtain are that this scheme is L² dissipative for any time step size and H¹ dissipative under a time‐step constraint. As a consequence, the existence of a numerical attractor of the discrete dynamical system is established. A by‐product of the dissipativity analysis of the linearly implicit Euler scheme for NSD is that the dissipativity of an implicit‐explicit scheme for the celebrated Navier‐Stokes equations that treats the volume forces term explicitly is obtained.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2114–2140, 2017     

Nonoverlapping domain decomposition characteristic finite differences for three‐dimensional convection‐diffusion equations

One domain decomposition method modified with characteristic differences is presented for non‐periodic three‐dimensional equations by multiply‐type quadratic interpolation and variant time‐step technique. This method consists of reduced‐scale, two‐dimensional computation on subdomain interface boundaries and fully implicit subdomain computation in parallel. A computational algorithm is outlined and an error estimate in discrete l²‐ norm is established by introducing new inner products and norms. Finally, numerical examples are given to illustrate the theoretical results, efficiency and parallelism of this method. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 17‐37, 2012     

Error estimates for a discontinuous Galerkin method with interior penalties applied to nonlinear Sobolev equations

A Discontinuous Galerkin method with interior penalties is presented for nonlinear Sobolev equations. A semi‐discrete and a family of fully‐discrete time approximate schemes are formulated. These schemes are symmetric. Hp‐version error estimates are analyzed for these schemes. For the semi‐discrete time scheme a priori L^∞(H¹) error estimate is derived and similarly, l^∞(H¹) and l²(H¹) for the fully‐discrete time schemes. These results indicate that spatial rates in H¹ and time truncation errors in L² are optimal. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A linearized Crank–Nicolson–Galerkin FEM for the time‐dependent Ginzburg–Landau equations under the temporal gauge

We propose a decoupled and linearized fully discrete finite element method (FEM) for the time‐dependent Ginzburg–Landau equations under the temporal gauge, where a Crank–Nicolson scheme is used for the time discretization. By carefully designing the time‐discretization scheme, we manage to prove the convergence rate , where τ is the time‐step size and r is the degree of the finite element space. Due to the degeneracy of the problem, the convergence rate in the spatial direction is one order lower than the optimal convergence rate of FEMs for parabolic equations. Numerical tests are provided to support our error analysis. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1279–1290, 2014     

Fully discrete interior penalty discontinuous Galerkin methods for nonlinear parabolic equations

In this article, we investigate interior penalty discontinuous Galerkin (IPDG) methods for solving a class of two‐dimensional nonlinear parabolic equations. For semi‐discrete IPDG schemes on a quasi‐uniform family of meshes, we obtain a priori bounds on solutions measured in the L² norm and in the broken Sobolev norm. The fully discrete IPDG schemes considered are based on the approximation by forward Euler difference in time and broken Sobolev space. Under a restriction related to the mesh size and time step, an hp ‐version of an a priori l^∞(L²) and l²(H¹) error estimate is derived and numerical experiments are presented.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 288–311, 2012     

Negative norm least‐squares methods for the velocity‐vorticity‐pressure Navier–Stokes equations

We develop and analyze a least‐squares finite element method for the steady state, incompressible Navier–Stokes equations, written as a first‐order system involving vorticity as new dependent variable. In contrast to standard L² least‐squares methods for this system, our approach utilizes discrete negative norms in the least‐squares functional. This allows us to devise efficient preconditioners for the discrete equations, and to establish optimal error estimates under relaxed regularity assumptions. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 237–256, 1999     

Error estimates of CFVE method for fully nonlinear convection‐dominated diffusion problems

Finite volume method and characteristics finite element method are two important methods for solving the partial differential equations. These two methods are combined in this paper to establish a fully discrete characteristics finite volume method for fully nonlinear convection‐dominated diffusion problems. Through detailed theoretical analysis, optimal order H¹ norm error estimates are obtained for this fully discrete scheme. Copyright © 2010 John Wiley & Sons, Ltd.     

A Crank‐Nicolson Petrov‐Galerkin method with quadrature for semi‐linear parabolic problems

We propose and analyze an application of a fully discrete C² spline quadrature Petrov‐Galerkin method for spatial discretization of semi‐linear parabolic initial‐boundary value problems on rectangular domains. We prove second order in time and optimal order H¹ norm convergence in space for the extrapolated Crank‐Nicolson quadrature Petrov‐Galerkin scheme. We demonstrate numerically both L² and H¹ norm optimal order convergence of the scheme even if the nonlinear source term is not smooth. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

Finite element approximation of a forward and backward anisotropic diffusion model in image denoising and form generalization

A new forward–backward anisotropic diffusion model is introduced. The two limit cases are the Perona‐Malik equation and the Total Variation flow model. A fully discrete finite element scheme is studied using C⁰‐piecewise linear elements in space and the backward Euler difference scheme in time. A priori estimates are proven. Numerical results in image denoising and form generalization are presented.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Approximation of the long‐term dynamics of the dynamical system generated by the multilayer quasigeostrophic equations of the ocean

In this article, we study the multilayer quasigeostrophic equations of the ocean. More precisely, we discretize these equations in time using the implicit Euler scheme and using the classical and uniform discrete Gronwall lemmas we prove that the approximate solution is uniformly bounded in H^?1, L² and H¹. Using the uniform stability of the scheme and the theory of the multivalued attractors, we then prove that the discrete attractors generated by the numerical scheme converge to the global attractor of the continuous system as the time‐step approaches zero. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1041–1065, 2016     

Positive numerical splitting method for the Hull and White 2D Black–Scholes equation

We consider the locally one‐dimensional backward Euler splitting method to solve numerically the Hull and White problem for pricing European options with stochastic volatility in the presence of a mixed derivative term. We prove the first‐order convergence of the time‐splitting. The parabolic equation degenerates on the boundary x = 0 and we apply a fitted finite volume scheme to the equation to resolve the degeneracy and derive the fully discrete problem as we also investigate the discrete maximum principle. Numerical experiments illustrate the efficiency of our difference scheme. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 822–846, 2015     

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     

Stability of the Crank–Nicolson–Adams–Bashforth scheme for the 2D Leray‐alpha model

We consider the stability of an efficient Crank–Nicolson–Adams–Bashforth method in time, finite element in space, discretization of the Leray‐α model. We prove finite‐time stability of the scheme in L², H¹, and H², as well as the long‐time L‐stability of the scheme under a Courant‐Freidrichs‐Lewy (CFL)‐type condition. Numerical experiments are given that are in agreement with the theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1155–1183, 2016     

A new fourth‐order numerical algorithm for a class of three‐dimensional nonlinear evolution equations

In this article, a new compact alternating direction implicit finite difference scheme is derived for solving a class of 3‐D nonlinear evolution equations. By the discrete energy method, it is shown that the new difference scheme has good stability and can attain second‐order accuracy in time and fourth‐order accuracy in space with respect to the discrete H¹ ‐norm. A Richardson extrapolation algorithm is applied to achieve fourth‐order accuracy in temporal dimension. Numerical experiments illustrate the accuracy and efficiency of the extrapolation algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Computational methods for delay parabolic and time‐fractional partial differential equations

This article is concerned with ?‐methods for delay parabolic partial differential equations. The methodology is extended to time‐fractional‐order parabolic partial differential equations in the sense of Caputo. The fully implicit scheme preserves delay‐independent asymptotic stability and the solution continuously depends on the time‐fractional order. Several numerical examples of interest are included to demonstrate the effectiveness of the method. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

An energy‐preserving Crank–Nicolson Galerkin method for Hamiltonian partial differential equations

A semidiscretization based method for solving Hamiltonian partial differential equations is proposed in this article. Our key idea consists of two approaches. First, the underlying equation is discretized in space via a selected finite element method and the Hamiltonian PDE can thus be casted to Hamiltonian ODEs based on the weak formulation of the system. Second, the resulting ordinary differential system is solved by an energy‐preserving integrator. The relay leads to a fully discretized and energy‐preserved scheme. This strategy is fully realized for solving a nonlinear Schrödinger equation through a combination of the Galerkin discretization in space and a Crank–Nicolson scheme in time. The order of convergence of our new method is if the discrete L²‐norm is employed. An error estimate is acquired and analyzed without grid ratio restrictions. Numerical examples are given to further illustrate the conservation and convergence of the energy‐preserving scheme constructed.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1485–1504, 2016     

New energy identities and super convergence analysis of the energy conserved splitting FDTD methods for 3D Maxwell's equations

The energy‐conserved splitting finite‐difference time‐domain (EC‐S‐FDTD) method has recently been proposed to solve the Maxwell equations with second order accuracy while numerically keep the L² energy conservation laws of the equations. In this paper, the EC‐S‐FDTD scheme for the 3D Maxwell equations is proved to be energy‐conserved and unconditionally stable in the discrete H¹ norm. The EC‐S‐FDTD scheme is of second‐order accuracy both in time step and spatial steps, which suggests the super‐convergence of this scheme in the discrete H¹ norm. And the divergence of the electric field of the EC‐S‐FDTD scheme in the discrete L² norm is second‐order accurate. Numerical experiments confirm our theoretical analysis. Copyright © 2012 John Wiley & Sons, Ltd.     

The Crank‐Nicolson/Adams‐Bashforth scheme for the time‐dependent Navier‐Stokes equations with nonsmooth initial data

In this article, we study the stability and convergence of the Crank‐Nicolson/Adams‐Bashforth scheme for the two‐dimensional nonstationary Navier‐Stokes equations with a nonsmooth initial data. A finite element method is applied for the spatial approximation of the velocity and pressure. The time discretization is based on the implicit Crank‐Nicolson scheme for the linear terms and the explicit Adams‐Bashforth scheme for the nonlinear term. Moreover, we prove that the scheme is almost unconditionally stable for a nonsmooth initial data u₀ with div u₀ = 0, i.e., the time step τ satisfies: τ ≤ C₀ if u₀ ∈ H¹ ∩ L^∞; τ |log h| ≤ C₀ if u₀ ∈ H¹ for the mesh size h and some positive constant C₀. Finally, we obtain some error estimates for the discrete velocity and pressure under the above stability condition. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 28: 155‐187, 2012