Error analysis of fully discrete velocity-correction methods for incompressible flows

A fully discrete version of the velocity-correction method, proposed by Guermond and Shen (2003) for the time-dependent Navier-Stokes equations, is introduced and analyzed. It is shown that, when accounting for space discretization, additional consistency terms, which vanish when space is not discretized, have to be added to establish stability and optimal convergence. Error estimates are derived for both the standard version and the rotational version of the method. These error estimates are consistent with those by Guermond and Shen (2003) as far as time discretiztion is concerned and are optimal in space for finite elements satisfying the inf-sup condition.

     

Extrapolation of the Nédélec element for the Maxwell equations by the mixed finite element method

In this paper, we use the integral-identity argument to obtain asymptotic error expansions for the mixed finite element approximation of the Maxwell equations on a rectangular mesh. The extrapolation method is applied to improve the accuracy of the approximation via an interpolation postprocessing technique. With the extrapolation, the approximation accuracy can be improved from O(h) to O(h ⁴) in the L ²-norm. Illustrative numerical results are given to demonstrate the higher order accuracy of the extrapolation method. This research was supported by the National Natural Science Foundation of China (No.10471103), Social Science Foundation of the Ministry of Education of China (06JA630047), Tianjin Natural Science Foundation (07JCYBJC14300).     

A partitioned second‐order method for magnetohydrodynamic flows at small magnetic reynolds numbers

This article aims to study the partitioned method for magnetohydrodynamic flows at small magnetic Reynolds numbers. We design a partitioned second‐order method and show that this method is stable under a time step () restrict condition. Our method can decouple the magnetohydrodynamic equations so that we can solve two relatively simple subproblems separately at each time step, which is computationally economic. A complete theoretical analysis of error estimates is also given. Finally, we present numerical experiments to support our theory.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1966–1986, 2017     

Error estimate of a fully discrete defect correction finite element method for unsteady incompressible Magnetohydrodynamics equations

In this study, a fully discrete defect correction finite element method for the unsteady incompressible Magnetohydrodynamics (MHD) equations, which is leaded by combining the Back Euler time discretization with the two-step defect correction in space, is presented. It is a continuous work of our formal paper [Math Method Appl Sci. 2017. DOI:10.1002/mma.4296]. The defect correction method is an iterative improvement technique for increasing the accuracy of a numerical solution without applying a grid refinement. Firstly, the nonlinear MHD equation is solved with an artificial viscosity term. Then, the numerical solutions are improved on the same grid by a linearized defect-correction technique. Then, we introduce the numerical analysis including stability analysis and error analysis. The numerical analysis proves that our method is stable and has an optimal convergence rate. Some numerical results [see Math Method Appl Sci. 2017. DOI:10.1002/mma.4296] show that this method is highly efficient for the unsteady incompressible MHD problems.     

Error correction iterative method for the stationary incompressible MHD flow

A new iterative finite element method for solving the stationary incompressible magnetohydrodynamics (MHD) equations is derived in this paper. The method consists of two steps at each iteration step, we need first to solve the MHD equations by the Oseen-type iterative scheme, and then an error correction strategy is applied to control the error arising from the linearization of the nonlinear MHD equations. The new method not only maintains the advantage of the standard Oseen-type scheme but also possesses a rapid rate of convergence. It is proved that the convergence rate of the proposed method is increased greatly under the uniqueness condition. The uniform stability and convergence of the new scheme are analyzed. Ample numerical experiments are performed to validate the accuracy and the efficiency of the new numerical scheme.     

Error analysis of a fully discrete finite element variational multiscale method for time‐dependent incompressible Navier–Stokes equations

A finite element variational multiscale method based on two local Gauss integrations is applied to solve numerically the time‐dependent incompressible Navier–Stokes equations. A significant feature of the method is that the definition of the stabilization term is derived via two local Guass integrations at element level, making it more efficient than the usual projection‐based variational multiscale methods. It is computationally cheap and gives an accurate approximation to the quantities sought. Based on backward Euler and Crank–Nicolson schemes for temporal discretization, we derive error bounds of the fully discrete solution which are first and second order in time, respectively. Numerical tests are also given to verify the theoretical predictions and demonstrate the effectiveness of the proposed method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Error estimates for a finite volume element method for parabolic equations in convex polygonal domains

We analyze the spatially semidiscrete piecewise linear finite volume element method for parabolic equations in a convex polygonal domain in the plane. Our approach is based on the properties of the standard finite element Ritz projection and also of the elliptic projection defined by the bilinear form associated with the variational formulation of the finite volume element method. Because the domain is polygonal, special attention has to be paid to the limited regularity of the exact solution. We give sufficient conditions in terms of data that yield optimal order error estimates in L₂ and H 1 . The convergence rate in the L_∞ norm is suboptimal, the same as in the corresponding finite element method, and almost optimal away from the corners. We also briefly consider the lumped mass modification and the backward Euler fully discrete method. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

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.     

Numerical analysis of two partitioned methods for uncoupling evolutionary MHD flows

Magnetohydrodynamics (MHD) studies the dynamics of electrically conducting fluids, involving Navier–Stokes (NSE) equations in fluid dynamics and Maxwell equations in eletromagnetism. The physical processes of fluid flows and electricity and magnetism are quite different and numerical simulations of each subprocess can require different meshes, time steps, and methods. In most terrestrial applications, MHD flows occur at low‐magnetic Reynold numbers. We introduce two partitioned methods to solve evolutionary MHD equations in such cases. The methods we study allow us at each time step to call NSE and Maxwell codes separately, each possibly optimized for the subproblem's respective physics. Complete error analysis and computational tests supporting the theory are given.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1083–1102, 2014     

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     

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     

A non-overlapping domain decomposition for low-frequency time-harmonic Maxwell’s equations in unbounded domains

In this paper, we are concerned with a non-overlapping domain decomposition method for solving the low-frequency time-harmonic Maxwell’s equations in unbounded domains. This method can be viewed as a coupling of finite elements and boundary elements in unbounded domains, which are decomposed into two subdomains with a spherical artificial boundary. We first introduce a discretization for the coupled variational problem by combining Nédélec edge elements of the lowest order and curvilinear elements. Then we design a D-N alternating method for solving the discrete problem. In the method, one needs only to solve the finite element problem (in a bounded domain) and calculate some boundary integrations, instead of solving a boundary integral equation. It will be shown that such iterative algorithm converges with a rate independent of the mesh size. The work of Qiya Hu was supported by Natural Science Foundation of China G10371129.     

Streamline diffusion finite element method for stationary incompressible magnetohydrodynamics

In this article, a streamline diffusion finite element method is proposed and analyzed for stationary incompressible magnetohydrodynamics (MHD) equations. This method is stable for any combinations of velocity, pressure, and magnet finite element spaces, without requiring Ladyzenskaja‐Babu?ka‐Brezzi (LBB) condition. The well‐posedness and convergence (at optimal error rate) of this scheme are proved in terms of some conditions. Two numerical experiments are illustrated to validate our theoretical analysis and show the streamline diffusion finite element approach is effective for solving the MHD problems. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1877–1901, 2014     

Richardson extrapolation of iterated discrete projection methods for eigenvalue approximation

In this paper, the eigenvalue approximation of a compact integral operator with a smooth kernel is discussed. We propose asymptotic error expansions of the iterated discrete Galerkin and iterated discrete collocation methods, and asymptotic error expansion of approximate eigenvalues. We then apply Richardson extrapolation to obtain higher order super-convergence of eigenvalue approximations. Numerical examples are presented to illustrate the theoretical estimate.     

Finite element error analysis of a variational multiscale method for the Navier-Stokes equations

The paper presents finite element error estimates of a variational multiscale method (VMS) for the incompressible Navier–Stokes equations. The constants in these estimates do not depend on the Reynolds number but on a reduced Reynolds number or on the mesh size of a coarse mesh. This work is partially supported by NSF grants DMS9972622, DMS20207627 and INT9814115.     

Improved stability and error analysis for a class of local projection stabilizations applied to the Oseen problem

We consider a class of local projection stabilizations with projection spaces defined on (possibly) overlapping sets applied to the Oseen problem. We prove that the underlying bilinear form satisfies an inf–sup condition with respect to a stronger norm than coercivity suggests. A modification of the stabilization of the convection allows an optimal estimation of the consistency error. A priori estimates in the stronger norm and in the L² norm for the pressure are established. Discontinuous pressure approximations are included in the analysis. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Error analysis of discontinuous Galerkin finite element method for optimal control problem governed by the transport equation

This article discusses a priori and a posteriori error estimates of discontinuous Galerkin finite element method for optimal control problem governed by the transport equation. We use variational discretization concept to discretize the control variable and discontinuous piecewise linear finite elements to approximate the state and costate variable. Based on the error estimates of discontinuous Galerkin finite element method for the transport equation, we get a priori and a posteriori error estimates for the transport equation optimal control problem. Finally, two numerical experiments are carried out to confirm the theoretical analysis.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1493–1512, 2017     

A projection technique for incompressible flow in the meshless finite volume particle method

The finite volume particle method is a meshless discretization technique, which generalizes the classical finite volume method by using smooth, overlapping and moving test functions applied in the weak formulation of the conservation law. The method was originally developed for hyperbolic conservation laws so that the compressible Euler equations particularly apply. In the present work we analyze the discretization error and enforce consistency by a new set of geometrical quantities. Furthermore, we introduce a discrete Laplace operator for the scheme in order to extend the method to second order partial differential equations. Finally, we transfer Chorins projection technique to the finite volume particle method in order to obtain a meshless scheme for incompressible flow. AMS subject classification 65M99, 68U20, 76B99, 76M12, 76M25, 76M28     

一种非线性对流扩散方程的特征有限元方法的误差估计

给出求解一种二维非线性对流扩散方程组的 Grank-Nicolson型特征有限元方法 ,并给出该方法的 H1模最优阶误差估计 .     

A posteriori error analysis for locally conservative mixed methods

In this work we present a theoretical analysis for a residual-type error estimator for locally conservative mixed methods. This estimator was first introduced by Braess and Verfürth for the Raviart-Thomas mixed finite element method working in mesh-dependent norms. We improve and extend their results to cover any locally conservative mixed method under minimal assumptions, in particular, avoiding the saturation assumption made by Braess and Verfürth. Our analysis also takes into account discontinuous coefficients with possibly large jumps across interelement boundaries. The main results are applied to the nonconforming finite element method and the interior penalty discontinuous Galerkin method as well as the mixed finite element method.