Least‐squares mixed finite element methods for the RLW equations

A least‐squares mixed finite element (LSMFE) schemes are formulated to solve the 1D regularized long wave (RLW) equations and the convergence is discussed. The L² error estimates of LSMFE methods for RLW equations under the standard regularity assumption on the finite element partition are given.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Superconvergence of discontinuous Galerkin finite element method for the stationary Navier‐Stokes equations

This article focuses on discontinuous Galerkin method for the two‐ or three‐dimensional stationary incompressible Navier‐Stokes equations. The velocity field is approximated by discontinuous locally solenoidal finite element, and the pressure is approximated by the standard conforming finite element. Then, superconvergence of nonconforming finite element approximations is applied by using least‐squares surface fitting for the stationary Navier‐Stokes equations. The method ameliorates the two noticeable disadvantages about the given finite element pair. Finally, the superconvergence result is provided under some regular assumptions. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 421–436, 2007     

Stabilized mixed finite element methods for the Navier‐Stokes equations with damping

In this paper, the stabilized mixed finite element methods are presented for the Navier‐Stokes equations with damping. The existence and uniqueness of the weak solutions are proven by use of the Brouwer fixed‐point theorem. Then, optimal error estimates for the H¹‐norm and L²‐norm of the velocity and the L²‐norm of the pressure are derived. Moreover, on the basis of the optimal L²‐norm error estimate of the velocity, a stabilized two‐step method is proposed, which is more efficient than the usual stabilized methods. Finally, two numerical examples are implemented to confirm the theoretical analysis.     

Least‐squares spectral collocation method for the Stokes equations

First‐order system least‐squares spectral collocation methods are presented for the Stokes equations by adopting the first‐order system and modifying the least‐squares functionals in 2 . Then homogeneous Legendre and Chebyshev (continuous and discrete) functionals are shown to be elliptic and continuous with respect to appropriate product weighted norms. The spectral convergence is analyzed for the proposed methods with some numerical experiments. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 128–139, 2004     

Study of the mixed finite volume method for Stokes and Navier‐Stokes equations

We present finite volume schemes for Stokes and Navier‐Stokes equations. These schemes are based on the mixed finite volume introduced in (Droniou and Eymard, Numer Math 105 (2006), 35‐71), and can be applied to any type of grid (without “orthogonality” assumptions as for classical finite volume methods) and in any space dimension. We present numerical results on some irregular grids, and we prove, for both Stokes and Navier‐Stokes equations, the convergence of the scheme toward a solution of the continuous problem. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

An augmented stress‐based mixed finite element method for the steady state Navier‐Stokes equations with nonlinear viscosity

A new stress‐based mixed variational formulation for the stationary Navier‐Stokes equations with constant density and variable viscosity depending on the magnitude of the strain tensor, is proposed and analyzed in this work. Our approach is a natural extension of a technique applied in a recent paper by some of the authors to the same boundary value problem but with a viscosity that depends nonlinearly on the gradient of velocity instead of the strain tensor. In this case, and besides remarking that the strain‐dependence for the viscosity yields a more physically relevant model, we notice that to handle this nonlinearity we now need to incorporate not only the strain itself but also the vorticity as auxiliary unknowns. Furthermore, similarly as in that previous work, and aiming to deal with a suitable space for the velocity, the variational formulation is augmented with Galerkin‐type terms arising from the constitutive and equilibrium equations, the relations defining the two additional unknowns, and the Dirichlet boundary condition. In this way, and as the resulting augmented scheme can be rewritten as a fixed‐point operator equation, the classical Schauder and Banach theorems together with monotone operators theory are applied to derive the well‐posedness of the continuous and associated discrete schemes. In particular, we show that arbitrary finite element subspaces can be utilized for the latter, and then we derive optimal a priori error estimates along with the corresponding rates of convergence. Next, a reliable and efficient residual‐based a posteriori error estimator on arbitrary polygonal and polyhedral regions is proposed. The main tools used include Raviart‐Thomas and Clément interpolation operators, inverse and discrete inequalities, and the localization technique based on triangle‐bubble and edge‐bubble functions. Finally, several numerical essays illustrating the good performance of the method, confirming the reliability and efficiency of the a posteriori error estimator, and showing the desired behavior of the adaptive algorithm, are reported. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1692–1725, 2017     

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     

A superconvergent nonconforming mixed finite element method for the Navier–Stokes equations

The superconvergence for a nonconforming mixed finite element approximation of the Navier–Stokes equations is analyzed in this article. The velocity field is approximated by the constrained nonconforming rotated Q₁ (CNRQ₁) element, and the pressure is approximated by the piecewise constant functions. Under some regularity assumptions, the superconvergence estimates for both the velocity in broken H¹‐norm and the pressure in L²‐norm are obtained. Some numerical examples are presented to demonstrate our theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 646–660, 2016     

Mixed finite element methods for incompressible flow: Stationary Stokes equations

In this article, we develop and analyze a mixed finite element method for the Stokes equations. Our mixed method is based on the pseudostress‐velocity formulation. The pseudostress is approximated by the Raviart‐Thomas (RT) element of index k ≥ 0 and the velocity by piecewise discontinuous polynomials of degree k. It is shown that this pair of finite elements is stable and yields quasi‐optimal accuracy. The indefinite system of linear equations resulting from the discretization is decoupled by the penalty method. The penalized pseudostress system is solved by the H(div) type of multigrid method and the velocity is then calculated explicitly. Alternative preconditioning approaches that do not involve penalizing the system are also discussed. Finally, numerical experiments are presented. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Superconvergence of a 3D finite element method for stationary Stokes and Navier‐Stokes problems

For the Poisson equation on rectangular and brick meshes it is well known that the piecewise linear conforming finite element solution approximates the interpolant to a higher order than the solution itself. In this article, this type of supercloseness property is established for a special interpolant of the Q₂ ? P element applied to the 3D stationary Stokes and Navier‐Stokes problem, respectively. Moreover, applying a Q₃ ? P postprocessing technique, we can also state a superconvergence property for the discretization error of the postprocessed discrete solution to the solution itself. Finally, we show that inhomogeneous boundary values can be approximated by the Lagrange Q₂‐interpolation without influencing the superconvergence property. Numerical experiments verify the predicted convergence rates. Moreover, a cost‐benefit analysis between the two third‐order methods, the post‐processed Q₂ ? P discretization, and the Q₃ ? P discretization is carried out. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

An analysis for the compressible Stokes equations by first‐order system of least‐squares finite element method

This article applies the first‐order system least‐squares (fosls) finite element method developed by Cai, Manteuffel and McCormick to the compressible Stokes equations. By introducing a new dependent velocity flux variable, we recast the compressible Stokes equations as a first‐order system. Then it is shown that the ellipticity and continuity hold for the least‐squares functionals employing the mixture of H^?1 and L², so that the fosls finite element methods yield best approximations for the velocity flux and velocity. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17:689–699, 2001     

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.     

Low‐order stabilized finite element methods for the unsteady Stokes/Navier‐Stokes equations with friction boundary conditions

In this paper, we consider low‐order stabilized finite element methods for the unsteady Stokes/Navier‐Stokes equations with friction boundary conditions. The time discretization is based on the Euler implicit scheme, and the spatial discretization is based on the low‐order element (P₁−P₁ or P₁−P₀) for the approximation of the velocity and pressure. Moreover, some error estimates for the numerical solution of fully discrete stabilized finite element scheme are obtained. Finally, numerical experiments are performed to confirm our theoretical results.     

Stability of a colocated finite volume scheme for the Navier‐Stokes equations

The aim of this article is to describe a colocated finite volume approximation of the incompressible Navier‐Stokes equation and study its stability. One of the advantages of colocated finite volume space discretizations over staggered space discretizations is that all the variables share the same location; hence, the possibility to more easily use complex geometries and hierarchical decompositions of the unknowns. The time discretization used in the scheme studied here is a projection method. First, we give the full discretization of the incompressible Navier‐Stokes equations, then, we state the stability result and prove it following the methods of Marion and Temam. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Least‐squares spectral method for velocity‐vorticity‐pressure form of the Stokes equations

The aim of this article is to present and analyze first‐order system least‐squares spectral method for the Stokes equations in two‐dimensional spaces. The Stokes equations are transformed into a first‐order system of equations by introducing vorticity as a new variable. The least‐squares functional is then defined by summing up the ‐ and ‐norms of the residual equations. The ‐norm in the least‐squares functional is replaced by suitable operator. Continuous and discrete homogeneous least‐squares functionals are shown to be equivalent to ‐norm of velocity and ‐norm of vorticity and pressure for spectral Galerkin and pseudospectral method. The spectral convergence of the proposed methods are given and the theory is validated by numerical experiment. Mass conservation is also briefly investigated. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 661–680, 2016     

Iterative methods in penalty finite element discretizations for the Steady Navier–Stokes equations

Three penalty finite element methods are designed to solve numerically the steady Navier–Stokes equations, where the Stokes, Newton, and Oseen iteration methods are used, respectively. Moreover, the stability analysis and error estimate for these nine algorithms are provided. Finally, the numerical tests confirm the theoretical results of the presented algorithms. Meanwhile, the numerical investigations are provided to show that the proposed methods are efficient for solving the steady Navier–Stokes equations with the different viscosity. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 74‐94, 2014     

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     

Normal and tangential continuous elements for least‐squares mixed finite element methods

We consider a finite element discretization of the primal first‐order least‐squares mixed formulation of the second‐order elliptic problem. The unknown variables are displacement and flux, which are approximated by equal‐order elements of the usual continuous element and the normal continuous element, respectively. We show that the error bounds for all variables are optimal. In addition, a field‐based least‐squares finite element method is proposed for the 3D‐magnetostatic problem, where both magnetic field and magnetic flux are taken as two independent variables which are approximated by the tangential continuous and the normal continuous elements, respectively. Coerciveness and optimal error bounds are obtained. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004.     

A multilevel finite element method in space‐time for the Navier‐Stokes problem

A multilevel finite element method in space‐time for the two‐dimensional nonstationary Navier‐Stokes problem is considered. The method is a multi‐scale method in which the fully nonlinear Navier‐Stokes problem is only solved on a single coarsest space‐time mesh; subsequent approximations are generated on a succession of refined space‐time meshes by solving a linearized Navier‐Stokes problem about the solution on the previous level. The a priori estimates and error analysis are also presented for the J‐level finite element method. We demonstrate theoretically that for an appropriate choice of space and time mesh widths: h_j ～ h, k_j ～ k, j = 2, …, J, the J‐level finite element method in space‐time provides the same accuracy as the one‐level method in space‐time in which the fully nonlinear Navier‐Stokes problem is solved on a final finest space‐time mesh. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Existence of global weak solutions for Navier‐Stokes‐Poisson equations with quantum effect and convergence to incompressible Navier‐Stokes equations

In this paper, we consider a three dimensional quantum Navier‐Stokes‐Poisson equations. Existence of global weak solutions is obtained, and convergence toward the classical solution of the incompressible Navier‐Stokes equation is rigorously proven for well prepared initial data. Furthermore, the associated convergence rates are also obtained. Copyright © 2014 John Wiley & Sons, Ltd.