Least‐squares mixed finite element methods for the incompressible Navier‐Stokes equations

Least‐squares mixed finite element schemes are formulated to solve the evolutionary Navier‐Stokes equations and the convergence is analyzed. We recast the Navier‐Stokes equations as a first‐order system by introducing a vorticity flux variable, and show that a least‐squares principle based on L² norms applied to this system yields optimal discretization error estimates. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 441–453, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10015     

Optimal estimates on stabilized finite volume methods for the incompressible Navier–Stokes model in three dimensions

Optimal estimates on stabilized finite volume methods for the three dimensional Navier–Stokes model are investigated and developed in this paper. Based on the global existence theorem [23], we first prove the global bound for the velocity in the H¹‐norm in time of a solution for suitably small data, and uniqueness of a suitably small solution by contradiction. Then, a full set of estimates is then obtained by some classical Galerkin techniques based on the relationship between finite element methods and finite volume methods approximated by the lower order finite elements for the three dimensional Navier–Stokes model.     

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 L2 finite element approximation for the incompressible Navier–Stokes equations

This paper utilizes the Picard method and Newton's method to linearize the stationary incompressible Navier–Stokes equations and then uses an LL* approach, which is a least-squares finite element method applied to the dual problem of the corresponding linear system. The LL* approach provides an L²-approximation to a given problem, which is not typically available with conventional finite element methods for nonlinear second-order partial differential equations. We first show that the proposed combination of linearization scheme and LL* approach provides an L²-approximation to the stationary incompressible Navier–Stokes equations. The validity of L²-approximation is proven through the analysis of the weak problem corresponding to the linearized Navier–Stokes equations. Then, the convergence is analyzed, and numerical results are presented.     

On Two Iterative Least-Squares Finite Element Schemes for the Incompressible Navier–Stokes Problem

This paper is mainly devoted to a comparative study of two iterative least-squares finite element schemes for solving the stationary incompressible Navier–Stokes equations with velocity boundary condition. Introducing vorticity as an additional unknown variable, we recast the Navier–Stokes problem into a first-order quasilinear velocity–vorticity–pressure system. Two Picard-type iterative least-squares finite element schemes are proposed to approximate the solution to the nonlinear first-order problem. In each iteration, we adopt the usual L ² least-squares scheme or a weighted L ² least-squares scheme to solve the corresponding Oseen problem and provide error estimates. We concentrate on two-dimensional model problems using continuous piecewise polynomial finite elements on uniform meshes for both iterative least-squares schemes. Numerical evidences show that the iterative L ² least-squares scheme is somewhat suitable for low Reynolds number flow problems, whereas for flows with relatively higher Reynolds numbers the iterative weighted L ² least-squares scheme seems to be better than the iterative L ² least-squares scheme. Numerical simulations of the two-dimensional driven cavity flow are presented to demonstrate the effectiveness of the iterative least-squares finite element approach.     

On the stability at all times of linearly extrapolated BDF2 timestepping for multiphysics incompressible flow problems

We prove long‐time stability of linearly extrapolated BDF2 (BDF2LE) timestepping methods, together with finite element spatial discretizations, for incompressible Navier‐Stokes equations (NSE) and related multiphysics problems. For the NSE, Boussinesq, and magnetohydrodynamics schemes, we prove unconditional long time L² stability, provided external forces (and sources) are uniformly bounded in time. We also provide numerical experiments to compare stability of BDF2LE to linearly extrapolated Crank‐Nicolson scheme for NSE, and find that BDF2LE has better stability properties, particularly for smaller viscosity values. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 999–1017, 2017     

Convergence analysis and computational testing of the finite element discretization of the Navier–Stokes alpha model

This report performs a complete analysis of convergence and rates of convergence of finite element approximations of the Navier–Stokes‐α (NS‐α) regularization of the NSE, under a zero‐divergence constraint on the velocity, to the true solution of the NSE. Convergence of the discrete NS‐α approximate velocity to the true Navier–Stokes velocity is proved and rates of convergence derived, under no‐slip boundary conditions. Generalization of the results herein to periodic boundary conditions is evident. Two‐dimensional experiments are performed, verifying convergence and predicted rates of convergence. It is shown that the NS‐α‐FE solutions converge at the theoretical limit of O(h²) when choosing α = h, in the H¹ norm. Furthermore, in the case of flow over a step the NS‐α model is shown to resolve vortex separation in the recirculation zone. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

The discrete duality finite volume method for a class of quasi‐Newtonian Stokes flows

In this paper, we propose a discrete duality finite volume (DDFV) scheme for the incompressible quasi‐Newtonian Stokes equation. The DDFV method is based on the use of discrete differential operators which satisfy some duality properties analogous to their continuous counterparts in a discrete sense. The DDFV method has a great ability to handle general geometries and meshes. In addition, every component of the velocity gradient can be reconstructed directly, which makes it suitable to deal with the nonlinear terms in the quasi‐Newtonian Stokes equation. We prove that the proposed DDFV scheme is uniquely solvable and of first‐order convergence in the discrete L²‐norms for the velocity, the strain rate tensor, and the pressure, respectively. Ample numerical tests are provided to highlight the performance of the proposed DDFV scheme and to validate the theoretical error analysis, in particular on locally refined nonconforming and polygonal meshes.     

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     

Stabilized DDFV schemes for stokes problem with variable viscosity on general 2D meshes

“Discrete Duality Finite Volume” schemes (DDFV for short) on general meshes are studied here for Stokes problems with variable viscosity with Dirichlet boundary conditions. The aim of this work is to analyze the well‐posedness of the scheme and its convergence properties. The DDFV method requires a staggered scheme, the discrete unknowns, the components of the velocity and the pressure, are located on different nodes. The scheme is stabilized using a finite volume analogue to Brezzi‐Pitkäranta techniques. This scheme is proved to be well‐posed on general meshes and to be first order convergent in a discrete H¹ ‐norm and a discrete L² ‐norm for respectively the velocity and the pressure. Finally, numerical experiments confirm the theoretical prediction, in particular on locally refined non conformal meshes. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1666–1706, 2011     

Non‐homogeneous Navier–Stokes systems with order‐parameter‐dependent stresses

We consider the Navier–Stokes system with variable density and variable viscosity coupled to a transport equation for an order‐parameter c. Moreover, an extra stress depending on c and ?c, which describes surface tension like effects, is included in the Navier–Stokes system. Such a system arises, e.g. for certain models of granular flows and as a diffuse interface model for a two‐phase flow of viscous incompressible fluids. The so‐called density‐dependent Navier–Stokes system is also a special case of our system. We prove short‐time existence of strong solution in L^q‐Sobolev spaces with q>d. We consider the case of a bounded domain and an asymptotically flat layer with a combination of a Dirichlet boundary condition and a free surface boundary condition. The result is based on a maximal regularity result for the linearized system. Copyright © 2010 John Wiley & Sons, Ltd.     

Finite volume element approximations of nonlocal reactive flows in porous media

In this article, we study finite volume element approximations for two‐dimensional parabolic integro‐differential equations, arising in the modeling of nonlocal reactive flows in porous media. These types of flows are also called NonFickian flows and exhibit mixing length growth. For simplicity, we consider only linear finite volume element methods, although higher‐order volume elements can be considered as well under this framework. It is proved that the finite volume element approximations derived are convergent with optimal order in H¹‐ and L²‐norm and are superconvergent in a discrete H¹‐norm. By examining the relationship between finite volume element and finite element approximations, we prove convergence in L^∞‐ and W^1,∞‐norms. These results are also new for finite volume element methods for elliptic and parabolic equations. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 285–311, 2000     

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     

Analysis of newton multilevel stabilized finite volume method for the three‐dimensional stationary Navier‐Stokes equations

This article proposes and analyzes a multilevel stabilized finite volume method(FVM) for the three‐dimensional stationary Navier–Stokes equations approximated by the lowest equal‐order finite element pairs. The method combines the new stabilized FVM with the multilevel discretization under the assumption of the uniqueness condition. The multilevel stabilized FVM consists of solving the nonlinear problem on the coarsest mesh and then performs one Newton correction step on each subsequent mesh thus only solving one large linear systems. The error analysis shows that the multilevel‐stabilized FVM provides an approximate solution with the convergence rate of the same order as the usual stabilized finite element solution solving the stationary Navier–Stokes equations on a fine mesh for an appropriate choice of mesh widths: h_j ～ h_j‐1², j = 1,…,J. Therefore, the multilevel stabilized FVM is more efficient than the standard one‐level‐stabilized FVM. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Periodic solutions of the Navier–Stokes equations in a perturbed half‐space and an aperture domain

We shall construct a periodic strong solution of the Navier–Stokes equations for some periodic external force in a perturbed half‐space and an aperture domain of the dimension n?3. Our proof is based on L^p–L^q estimates of the Stokes semigroup. We apply L^p–L^q estimates to the integral equation which is transformed from the original equation. As a result, we obtain the existence and uniqueness of periodic strong solutions. Copyright © 2005 John Wiley & Sons, Ltd.     

Existence and exponential stability in Lr‐spaces of stationary Navier–Stokes flows with prescribed flux in infinite cylindrical domains

We prove existence, uniqueness and exponential stability of stationary Navier–Stokes flows with prescribed flux in an unbounded cylinder of ?ⁿ,n?3, with several exits to infinity provided the total flux and external force are sufficiently small. The proofs are based on analytic semigroup theory, perturbation theory and L^r ? L^q‐estimates of a perturbation of the Stokes operator in L^q‐spaces. Copyright © 2006 John Wiley & Sons, Ltd.     

Error estimates for two‐level penalty finite volume method for the stationary Navier–Stokes equations

Two‐level penalty finite volume method for the stationary Navier–Stokes equations based on the P₁ ? P₀ element is considered in this paper. The method involves solving one small penalty Navier–Stokes problem on a coarse mesh with mesh size H = ?^{1 / 4}h^{1 / 2}, a large penalty Stokes problem on a fine mesh with mesh size h, where 0 < ? < 1 is a penalty parameter. The method we study provides an approximate solution with the convergence rate of same order as the penalty finite volume solution (u_?h,p_?h), which involves solving one large penalty Navier–Stokes problem on a fine mesh with the same mesh size h. However, our method can save a large amount of computational time. Copyright © 2013 John Wiley & Sons, Ltd.     

Global Regularity of the 3D Axi-Symmetric Navier–Stokes Equations with Anisotropic Data

In this paper, we study the 3D axisymmetric Navier–Stokes equations with swirl. We prove the global regularity of the 3D Navier–Stokes equations for a family of large anisotropic initial data. Moreover, we obtain a global bound of the solution in terms of its initial data in some L ^p norm. Our results also reveal some interesting dynamic growth behavior of the solution due to the interaction between the angular velocity and the angular vorticity fields.     

Strong solutions to 3D compressible magnetohydrodynamic equations with Navier‐slip condition

We consider the short time strong solutions to the compressible magnetohydrodynamic equations with initial vacuum, in which the velocity field satisfies the Navier‐slip condition. The Navier‐slip condition differs in many aspects from no‐slip conditions, and it has attracted considerable attention in nanoscale and microscale flows research. Inspired by Kato and Lax's idea, we use the Lax–Milgram theorem and contraction mapping argument to prove local existence. Moreover, under the Navier‐slip condition, we establish a criterion for possible breakdown of such solutions at finite time in terms of the temporal integral of L^∞ norm of the deformation tensor D(u). Copyright © 2015 John Wiley & Sons, Ltd.     

Finite volume method based on stabilized finite elements for the nonstationary Navier–Stokes problem

A finite volume method based on stabilized finite element for the two‐dimensional nonstationary Navier–Stokes equations is investigated in this work. As in stabilized finite element method, macroelement condition is introduced for constructing the local stabilized formulation of the nonstationary Navier–Stokes equations. Moreover, for P₁ ? P₀ element, the H¹ error estimate of optimal order for finite volume solution (u_h,p_h) is analyzed. And, a uniform H¹ error estimate of optimal order for finite volume solution (u_h, p_h) is also obtained if the uniqueness condition is satisfied. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007