Stabilized multiphysics finite element method with Crank–Nicolson scheme for a poroelasticity model

In the paper, we propose a stabilized multiphysics finite element method with Crank–Nicolson scheme for a poroelasticity model. The method can eliminate the locking phenomenon and reveal the multi‐physical process. The lowest equal order finite element pair is used to reduce the computational cost. Furthermore, the method needs no constraint condition Δt = O(h²) and achieves optimal convergent order. Numerical tests are provided to illustrate the optimal accuracy and good performance in eliminating locking phenomenon of the method.     

The stability and convergence of two linearized finite difference schemes for the nonlinear epitaxial growth model

The numerical simulation of the dynamics of the molecular beam epitaxy (MBE) growth is considered in this article. The governing equation is a nonlinear evolutionary equation that is of linear fourth order derivative term and nonlinear second order derivative term in space. The main purpose of this work is to construct and analyze two linearized finite difference schemes for solving the MBE model. The linearized backward Euler difference scheme and the linearized Crank‐Nicolson difference scheme are derived. The unique solvability, unconditional stability and convergence are proved. The linearized Euler scheme is convergent with the convergence order of O(τ + h²) and linearized Crank‐Nicolson scheme is convergent with the convergence order of O(τ² + h²) in discrete L²‐norm, respectively. Numerical stability with respect to the initial conditions is also obtained for both schemes. Numerical experiments are carried out to demonstrate the theoretical analysis. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Error estimates of expanded mixed methods for optimal control problems governed by hyperbolic integro‐differential equations

In this article, we investigate the L^∞(L²) ‐error estimates of the semidiscrete expanded mixed finite element methods for quadratic optimal control problems governed by hyperbolic integrodifferential equations. The state and the costate are discretized by the order k Raviart‐Thomas mixed finite element spaces, and the control is approximated by piecewise polynomials of order k(k ≥ 0). We derive error estimates for both the state and the control approximation. Numerical experiments are presented to test the theoretical results. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A Petrov‐Galerkin method with quadrature for semi‐linear second‐order hyperbolic problems

We propose and analyze a Crank–Nicolson quadrature Petrov–Galerkin (CNQPG) ‐spline method for solving semi‐linear second‐order hyperbolic initial‐boundary value problems. We prove second‐order convergence in time and optimal order H² norm convergence in space for the CNQPG scheme that requires only linear algebraic solvers. We demonstrate numerically optimal order H^k, k = 0,1,2, norm convergence of the scheme for some test problems with smooth and nonsmooth nonlinearities. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A Priori Error Estimates for Finite Element Methods for H (2,1)-Elliptic Equations

The convergence of finite element methods for linear elliptic boundary value problems of second and forth order is well understood. In this article, we introduce finite element approximations of some linear semi-elliptic boundary value problem of mixed order on a two-dimensional rectangular domain Q. The equation is of second order in one direction and forth order in the other and appears in the optimal control of parabolic partial differential equations if one eliminates the control and the state (or the adjoint state) in the first order optimality conditions. We establish a regularity result and estimate for the finite element error of conforming approximations of this equation. The finite elements in use have a tensor product structure, in one dimension we use linear, quadratic or cubic Lagrange elements in the other dimension cubic Hermite elements. For these elements, we prove the error bound O(h ² + τ ^k ) in the energy norm and O((h ² + τ ^k )(h ² + τ)) in the L ²(Q)-norm.     

Error estimates and superconvergence of semidiscrete mixed methods for optimal control problems governed by hyperbolic equations

In this paper, we investigate the L ^??(L ₂)-error estimates and superconvergence of the semidiscrete mixed finite elementmethods for quadratic optimal control problems governed by linear hyperbolic equations. The state and the co-state are discretized by the order k Raviart-Thomas mixed finite element spaces and the control is approximated by piecewise polynomials of order k(k ?? 0). We derive error estimates for approximation of both state and control. Moreover, we present the superconvergence analysis for mixed finite element approximation of the optimal control problems.     

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.     

Superconvergence of H 1-Galerkin mixed finite element methods for parabolic problems

In this article, we study the semidiscrete H ¹-Galerkin mixed finite element method for parabolic problems over rectangular partitions. The well-known optimal order error estimate in the L ²-norm for the flux is of order 𝒪(h ^k+1) (SIAM J. Numer. Anal. 35 (2), (1998), pp. 712–727), where k ≥ 1 is the order of the approximating polynomials employed in the Raviart–Thomas element. We derive a superconvergence estimate of order 𝒪(h ^k+3) between the H ¹-Galerkin mixed finite element approximation and an appropriately defined local projection of the flux variable when k ≥ 1. A the new approximate solution for the flux with superconvergence of order 𝒪(h ^k+3) is realized via a postprocessing technique using local projection methods.     

A second order nonconforming rectangular finite element method for approximating Maxwell’s equations

The main objective of this paper is to present a new rectangular nonconforming finite element scheme with the second order convergence behavior for approximation of Maxwell’s equations.Then the corresponding optimal error estimates are derived.The difficulty in construction of this finite element scheme is how to choose a compatible pair of degrees of freedom and shape function space so as to make the consistency error due to the nonconformity of the element being of order O(h 3 ) ,properly one order higher than that of its interpolation error O(h 2 ) in the broken energy norm,where h is the subdivision parameter tending to zero.     

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     

Quintic trigonometric spline based numerical scheme for nonlinear modified Burgers' equation

This paper presents a numerical method based on quintic trigonometric B‐splines for solving modified Burgers' equation (MBE). Here, the MBE is first discretized in time by Crank–Nicolson scheme and the resulting scheme is solved by quintic trigonometric B‐splines. The proposed method tackles nonlinearity by using a linearization process known as quasilinearization. A rigorous analysis of the stability and convergence of the proposed method are carried out, which proves that the method is unconditionally stable and has order of convergence O(h⁴ + k²). Numerical results presented are very much in accordance with the exact solution, which is established by the negligible values of L₂ and L_∞ errors. Computational efficiency of the scheme is proved by small values of CPU time. The method furnishes results better than those obtained by using most of the existing methods for solving MBE.     

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     

Spline methods for the solution of hyperbolic equation with variable coefficients

In this study, we developed the methods based on nonpolynomial cubic spline for numerical solution of second‐order nonhomogeneous hyperbolic partial differential equation. Using nonpolynomial cubic spline in space and finite difference in time directions, we obtained the implicit three level methods of O(k² + h²) and O(k² + h⁴). The proposed methods are applicable to the problems having singularity at x = 0, too. Stability analysis of the presented methods have been carried out. The presented methods are applied to the nonhomogeneous examples of different types. Numerical comparison with Mohanty's method (Mohanty, Appl Math Comput, 165 (2005), 229–236) shows the superiority of our presented schemes. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

An optimal nonlinear Galerkin method with mixed finite elements for the steady Navier‐Stokes equations

An optimal nonlinear Galerkin method with mixed finite elements is developed for solving the two‐dimensional steady incompressible Navier‐Stokes equations. This method is based on two finite element spaces X_H and X_h for the approximation of velocity, defined on a coarse grid with grid size H and a fine grid with grid size h ? H, respectively, and a finite element space M_h for the approximation of pressure. We prove that the difference in appropriate norms between the solutions of the nonlinear Galerkin method and a classical Galerkin method is of the order of H⁵. If we choose H = O(h^2/5), these two methods have a convergence rate of the same order. We numerically demonstrate that the optimal nonlinear Galerkin method is efficient and can save a large amount of computational time. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 762–775, 2003.     

High‐order numerical methods for one‐dimensional parabolic singularly perturbed problems with regular layers

In this work we construct and analyze some finite difference schemes used to solve a class of time‐dependent one‐dimensional convection‐diffusion problems, which present only regular layers in their solution. We use the implicit Euler or the Crank‐Nicolson method to discretize the time variable and a HODIE finite difference scheme, defined on a piecewise uniform Shishkin mesh, to discretize the spatial variable. In both cases we prove that the numerical method is uniformly convergent with respect to the diffusion parameter, having order near two in space and order one or 3/2, depending on the method used, in time. We show some numerical examples which illustrate the theoretical results, in the case of using the Euler implicit method, and give better numerical behaviour than that predicted theoretically, showing order two in time and order N^?2log²N in space, if the Crank‐Nicolson scheme is used to discretize the time variable. Finally, we construct a numerical algorithm by combining a third order A‐stable SDIRK with two stages and a third‐order HODIE difference scheme, showing its uniformly convergent behavior, reaching order three, up to a logarithmic factor. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Convergence analysis of a linearized Crank–Nicolson scheme for the two‐dimensional complex Ginzburg–Landau equation

A linearized Crank–Nicolson‐type scheme is proposed for the two‐dimensional complex Ginzburg–Landau equation. The scheme is proved to be unconditionally convergent in the L² ‐norm by the discrete energy method. The convergence order is \begin{align*}\mathcal{O}(\tau^2+h_1^2+h^2_2)\end{align*}, where τ is the temporal grid size and h₁,h₂ are spatial grid sizes in the x ‐ and y ‐directions, respectively. A numerical example is presented to support the theoretical result. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A quadrature finite element method for semilinear second‐order hyperbolic problems

In this work we propose and analyze a fully discrete modified Crank–Nicolson finite element (CNFE) method with quadrature for solving semilinear second‐order hyperbolic initial‐boundary value problems. We prove optimal‐order convergence in both time and space for the quadrature‐modified CNFE scheme that does not require nonlinear algebraic solvers. Finally, we demonstrate numerically the order of convergence of our scheme for some test problems. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

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     

Approximations of linear control problems with bang-bang solutions

We analyse the Euler discretization to a class of linear optimal control problems. First we show convergence of order h for the discrete approximation of the adjoint solution and the switching function, where h is the mesh size. Under the additional assumption that the optimal control has bang-bang structure we show that the discrete and the exact controls coincide except on a set of measure O(h). As a consequence, the discrete optimal control approximates the optimal control with order 1 w.r.t. the L ¹-norm and with order 1/2 w.r.t. the L ²-norm. An essential assumption is that the slopes of the switching function at its zeros are bounded away from zero which is in fact an inverse stability condition for these zeros. We also discuss higher order approximation methods based on the approximation of the adjoint solution and the switching function. Several numerical examples underline the results.     

The theta‐Galerkin finite element method for coupled systems resulting from microsensor thermistor problems

This paper is devoted to the analysis of a linearized theta‐Galerkin finite element method for the time‐dependent coupled systems resulting from microsensor thermistor problems. Hereby, we focus on time discretization based on θ‐time stepping scheme with including the standard Crank‐Nicolson ( ) and the shifted Crank‐Nicolson ( , where δ is the time‐step) schemes. The semidiscrete formulation in space is presented and optimal error bounds in L²‐norm and the energy norm are established. For the fully discrete system, the optimal error estimates are derived for the standard Crank‐Nicolson, the shifted Crank‐Nicolson, and the general case where with k=0,1 . Finally, numerical simulations that validate the theoretical findings are exhibited.