Stabilized multiscale finite element method for the stationary Navier-Stokes equations

In the paper, a stabilized multiscale finite element method for the stationary incompressible Navier-Stokes equations is considered. The method is a Petrov-Galerkin approach based on the multiscale enrichment of the standard polynomial space enriched with the unusual bubble functions which no longer vanish on every element boundary for the velocity space. The stability of the P₁-P₀ triangular element (or the Q₁-P₀ quadrilateral element) is established. And the optimal error estimates of the stabilized multiscale finite element method for the stationary Navier-Stokes equations are obtained.     

A fully discrete stabilized finite element method for the time-dependent Navier-Stokes equations

In this article, we consider a fully discrete stabilized finite element method based on two local Gauss integrations for the two-dimensional time-dependent Navier-Stokes equations. It focuses on the lowest equal-order velocity-pressure pairs. Unlike the other stabilized method, the present approach does not require specification of a stabilization parameter or calculation of higher-order derivatives, and always leads to a symmetric linear system. The Euler semi-implicit scheme is used for the time discretization. It is shown that the proposed fully discrete stabilized finite element method results in the optimal order bounds for the velocity and pressure.     

Stability of a finite element method for 3D exterior stationary Navier-Stokes flows

We consider numerical approximations of stationary incompressible Navier-Stokes flows in 3D exterior domains, with nonzero velocity at infinity. It is shown that a P1-P1 stabilized finite element method proposed by C. Rebollo: A term by term stabilization algorithm for finite element solution of incompressible flow problems, Numer. Math. 79 (1998), 283–319, is stable when applied to a Navier-Stokes flow in a truncated exterior domain with a pointwise boundary condition on the artificial boundary.     

Convergence and stability of two-level penalty mixed finite element method for stationary Navier-Stokes equations

The two-level penalty mixed finite element method for the stationary Navier-Stokes equations based on Taylor-Hood element is considered in this paper. Two algorithms are proposed and analyzed. Moreover, the optimal stability analysis and error estimate for these two algorithms are provided. Finally, the numerical tests confirm the theoretical results of the presented algorithms.     

Two-level multiscale finite element methods for the steady navier-stokes problem

In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.     

Optimal error estimate of the penalty finite element method for the time-dependent Navier-Stokes equations

A fully discrete penalty finite element method is presented for the two-dimensional time-dependent Navier-Stokes equations. The time discretization of the penalty Navier-Stokes equations is based on the backward Euler scheme; the spatial discretization of the time discretized penalty Navier-Stokes equations is based on a finite element space pair which satisfies some approximate assumption. An optimal error estimate of the numerical velocity and pressure is provided for the fully discrete penalty finite element method when the parameters and are sufficiently small.

     

A new defect correction method for the Navier-Stokes equations at high Reynolds numbers

In this paper, a new defect correction method for the Navier-Stokes equations is presented. With solving an artificial viscosity stabilized nonlinear problem in the defect step, and correcting the residual by linearized equations in the correction step for a few steps, this combination is particularly efficient for the Navier-Stokes equations at high Reynolds numbers. In both the defect and correction steps, we use the Oseen iterative scheme to solve the discrete nonlinear equations. Furthermore, the stability and convergence of this new method are deduced, which are better than that of the classical ones. Finally, some numerical experiments are performed to verify the theoretical predictions and show the efficiency of the new combination.     

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.     

The Crouzeix-Raviart type nonconforming finite element method for the nonstationary Navier-Stokes equations on anisotropic meshes

This paper is devoted to study the Crouzeix-Raviart （C-R） type nonconforming linear triangular finite element method （FEM） for the nonstationary Navier-Stokes equations on anisotropic meshes. By intro- ducing auxiliary finite element spaces, the error estimates for the velocity in the L2-norm and energy norm, as well as for the pressure in the L2-norm are derived.     

Interior maximum norm estimates for finite element discretizations of the Stokes equations

Interior estimates are proved in the L ^∞ norm for stable finite element discretizations of the Stokes equations on translation invariant meshes. These estimates yield information about the quality of the finite element solution in subdomains a positive distance from the boundary. While they have been established for second-order elliptic problems, these interior, or local, maximum norm estimates for the Stokes equations are new. By applying finite differenciation methods on a translation invariant mesh, we obtain optimal convergence rates in the mesh size h in the maximum norm. These results can be used for analyzing superconvergence in finite element methods for the Stokes equations.     

A priori estimates for two multiscale finite element methods using multiple global fields to wave equations

We consider a scalar wave equation with nonseparable spatial scales. If the solution of the wave equation smoothly depends on some global fields, then we can utilize the global fields to construct multiscale finite element basis functions. We present two finite element approaches using the global fields: partition of unity method and mixed multiscale finite element method. We derive a priori error estimates for the two approaches and theoretically investigate the relation between the smoothness of the global fields and convergence rates of the approximations for the wave equation. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011     

Dimension splitting method for the three dimensional rotating Navier-Stokes equations

In this paper, we propose a dimensional splitting method for the three dimensional (3D) rotating Navier-Stokes equations. Assume that the domain is a channel bounded by two surfaces and is decomposed by a series of surfaces ■i into several sub-domains, which are called the layers of the flow. Every interface i between two sub-domains shares the same geometry. After establishing a semi-geodesic coordinate (S-coordinate) system based on ■i , Navier-Stoke equations in this coordinate can be expressed as the sum of two operators, of which one is called the membrane operator defined on the tangent space on ■i , another one is called the bending operator taking value in the normal space on ■i . Then the derivatives of velocity with respect to the normal direction of the surface are approximated by the Euler central difference, and an approximate form of Navier-Stokes equations on the surface ■i is obtained, which is called the two-dimensional three-component (2D-3C) Navier-Stokes equations on a two dimensional manifold. Solving these equations by alternate iteration, an approximate solution to the original 3D Navier-Stokes equations is obtained. In addition, the proof of the existence of solutions to 2D-3C Navier-Stokes equations is provided, and some approximate methods for solving 2D-3C Navier-Stokes equations are presented.     

Analysis of a sequential regularization method for the unsteady Navier-Stokes equations

The incompressibility constraint makes Navier-Stokes equations difficult. A reformulation to a better posed problem is needed before solving it numerically. The sequential regularization method (SRM) is a reformulation which combines the penalty method with a stabilization method in the context of constrained dynamical systems and has the benefit of both methods. In the paper, we study the existence and uniqueness for the solution of the SRM and provide a simple proof of the convergence of the solution of the SRM to the solution of the Navier-Stokes equations. We also give error estimates for the time discretized SRM formulation.

     

Adaptive finite element methods for Boussinesq equations

An a posteriori error analysis for Boussinesq equations is derived in this article. Then we compare this new estimate with a previous one developed for a regularized version of Boussinesq equations in a previous work. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 214–236, 2000     

Application of the least-squares spectral element method using Chebyshev polynomials to solve the incompressible Navier-Stokes equations

In this paper the extension of the Legendre least-squares spectral element formulation to Chebyshev polynomials will be explained. The new method will be applied to the incompressible Navier-Stokes equations and numerical results, obtained for the lid-driven cavity flow at Reynolds numbers varying between 1000 and 7500, will be compared with the commonly used benchmark results. The new results reveal that the least-squares spectral element formulations based on the Legendre and Chebyshev Gauss-Lobatto Lagrange interpolating polynomials are equally accurate.     

Two-level pressure projection finite element methods for Navier-Stokes equations with nonlinear slip boundary conditions

The two-level pressure projection stabilized finite element methods for Navier-Stokes equations with nonlinear slip boundary conditions are investigated in this paper, whose variational formulation is the Navier-Stokes type variational inequality problem of the second kind. Based on the P₁-P₁ triangular element and using the pressure projection stabilized finite element method, we solve a small Navier-Stokes type variational inequality problem on the coarse mesh with mesh size H and solve a large Stokes type variational inequality problem for simple iteration or a large Oseen type variational inequality problem for Oseen iteration on the fine mesh with mesh size h. The error analysis obtained in this paper shows that if h=O(H²), the two-level stabilized methods have the same convergence orders as the usual one-level stabilized finite element methods, which is only solving a large Navier-Stokes type variational inequality problem on the fine mesh. Finally, numerical results are given to verify the theoretical analysis.     

Unified analysis for stabilized methods of low-order mixed finite elements for stationary Navier-Stokes equations

A unified analysis is presented for the stabilized methods including the pres-sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed fi...     

Error estimates for an operator-splitting method for Navier-Stokes equations: Second-order schemes

We present in this paper an error analysis of a fractional-step method for the approximation of the unsteady incompressible Navier-Stokes equations. Under mild regularity assumptions on the continuous solution, we obtain second-order error estimates in the time step size, both for velocity and pressure. Numerical results in agreement with the error analysis are also presented.     

A two-stage least-squares finite element method for the stress-pressure-displacement elasticity equations

A new stress-pressure-displacement formulation for the planar elasticity equations is proposed by introducing the auxiliary variables, stresses, and pressure. The resulting first-order system involves a nonnegative parameter that measures the material compressibility for the elastic body. A two-stage least-squares finite element procedure is introduced for approximating the solution to this system with appropriate boundary conditions. It is shown that the two-stage least-squares scheme is stable and, with respect to the order of approximation for smooth exact solutions, the rates of convergence of the approximations for all the unknowns are optimal both in the H¹-norm and in the L²-norm. Numerical experiments with various values of the parameter are examined, which demonstrate the theoretical estimates. Among other things, computational results indicate that the behavior of convergence is uniform in the nonnegative parameter. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 297–315, 1998