Existence and zero‐electron‐mass limit of strong solutions to the stationary compressible Navier‐Stokes‐Poisson equation with large external force

In this study, we consider a viscous compressible model of plasma and semiconductors, which is expressed as a compressible Navier‐Stokes‐Poisson equation. We prove that there exists a strong solution to the boundary value problem of the steady compressible Navier‐Stokes‐Poisson equation with large external forces in bounded domain, provided that the ratio of the electron/ions mass is appropriately small. Moreover, the zero‐electron‐mass limit of the strong solutions is rigorously verified. The main idea in the proof is to split the original equation into 4 parts, a system of stationary incompressible Navier‐Stokes equations with large forces, a system of stationary compressible Navier‐Stokes equations with small forces, coupled with 2 Poisson equations. Based on the known results about linear incompressible Navier‐Stokes equation, linear compressible Navier‐Stokes, linear transport, and Poisson equations, we try to establish uniform in the ratio of the electron/ions mass a priori estimates. Further, using Schauder fixed point theorem, we can show the existence of a strong solution to the boundary value problem of the steady compressible Navier‐Stokes‐Poisson equation with large external forces. At the same time, from the uniform a priori estimates, we present the zero‐electron‐mass limit of the strong solutions, which converge to the solutions of the corresponding incompressible Navier‐Stokes‐Poisson equations.     

Increasing accuracy and efficiency in FE computations of the Leray‐Deconvolution model

This article develops, analyzes, and tests a finite element method for approximating solutions to the Leray‐deconvolution regularization of the Navier‐Stokes equations. The scheme combines three ideas to create an accurate and effective algorithm: the use of an incompressible filter, a linearization that decouples the velocity‐pressure system from the filtering and deconvolution operations, and a stabilization that works well with the linearization. A rigorous and complete numerical analysis of the scheme is given, and numerical experiments are presented that show clear advantages of the scheme. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 720–736, 2012     

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     

Uniform regularity and vanishing viscosity limit for the chemotaxis‐Navier‐Stokes system in a 3D bounded domain

We investigate the uniform regularity and vanishing viscosity limit for the incompressible chemotaxis‐Navier‐Stokes system with Navier boundary condition for velocity field and Neumann boundary condition for cell density and chemical concentration in a 3D bounded domain. It is shown that there exists a unique strong solution of the incompressible chemotaxis‐Navier‐Stokes system in a finite time interval, which is independent of the viscosity coefficient. Moreover, this solution is uniformly bounded in a conormal Sobolev space, which allows us to take the vanishing viscosity limit to obtain the incompressible chemotaxis‐Euler system.     

A Two‐Parameter Stabilized Finite Element Method for Incompressible Flows

Based on two‐grid discretizations, a two‐parameter stabilized finite element method for the steady incompressible Navier–Stokes equations at high Reynolds numbers is presented and studied. In this method, a stabilized Navier–Stokes problem is first solved on a coarse grid, and then a correction is calculated on a fine grid by solving a stabilized linear problem. The stabilization term for the nonlinear Navier–Stokes equations on the coarse grid is based on an elliptic projection, which projects higher‐order finite element interpolants of the velocity into a lower‐order finite element interpolation space. For the linear problem on the fine grid, either the same stabilization approach (with a different stabilization parameter) as that for the coarse grid problem or a completely different stabilization approach could be employed. Error bounds for the discrete solutions are estimated. Algorithmic parameter scalings of the method are also derived. The theoretical results show that, with suitable scalings of the algorithmic parameters, this method can yield an optimal convergence rate. Numerical results are provided to verify the theoretical predictions and demonstrate the effectiveness of the proposed method. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 425–444, 2017     

A parallel two-level finite element variational multiscale method for the Navier–Stokes equations

A combination method of two-grid discretization approach with a recent finite element variational multiscale algorithm for simulation of the incompressible Navier–Stokes equations is proposed and analyzed. The method consists of a global small-scale nonlinear Navier–Stokes problem on a coarse grid and local linearized residual problems in overlapped fine grid subdomains, where the numerical form of the Navier–Stokes equations on the coarse grid is stabilized by a stabilization term based on two local Gauss integrations at element level and defined by the difference between a consistent and an under-integrated matrix involving the gradient of velocity. By the technical tool of local a priori estimate for the finite element solution, error bounds of the discrete solution are estimated. Algorithmic parameter scalings are derived. Numerical tests are also given to verify the theoretical predictions and demonstrate the effectiveness of the method.     

A stabilization algorithm of the Navier–Stokes equations based on algebraic Bernoulli equation

In this paper, we consider the stabilization of the nonstationary incompressible Navier–Stokes equations around a stationary solution by a boundary linear feedback control. The feedback operator is obtained from the solution of the algebraic Bernoulli equation associated with the penalized linearized Navier–Stokes equations around an unstable stationary solution and is used to locally stabilize the original nonlinear equations. We give the explicit factorized form of the stabilizing solution of the algebraic Bernoulli equation. The numerical effectiveness of this approach is demonstrated by stabilizing the vortex shedding behind a circular obstacle. Copyright © 2011 John Wiley & Sons, Ltd.     

Application of a posteriori error estimates for Navier‐Stokes equations to fluid flow in a channel with corners

We investigate a posteriori error estimates for the Stokes and Navier‐Stokes equations for incompressible fluid problem on 2D polygonal domains. We have determined accurately the constant that appears in the estimates. We apply these results to the technique of adaptive mesh refinement ‐ we solve an incompressible flow problem in a domain with corners that cause singularities in the solution.     

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.     

Stability and convergence of efficient Navier‐Stokes solvers via a commutator estimate

For strong solutions of the incompressible Navier‐Stokes equations in bounded domains with velocity specified at the boundary, we establish the unconditional stability and convergence of discretization schemes that decouple the updates of pressure and velocity through explicit time stepping for pressure. These schemes require no solution of stationary Stokes systems, nor any compatibility between velocity and pressure spaces to ensure an inf‐sup condition, and are representative of a class of highly efficient computational methods that have recently emerged. The proofs are simple, based upon a new, sharp estimate for the commutator of the Laplacian and Helmholtz projection operators. This allows us to treat an unconstrained formulation of the Navier‐Stokes equations as a perturbed diffusion equation. © 2007 Wiley Periodicals, Inc.     

First‐order system least squares for the Oseen equations

Following earlier work for Stokes equations, a least squares functional is developed for two‐ and three‐dimensional Oseen equations. By introducing a velocity flux variable and associated curl and trace equations, ellipticity is established in an appropriate product norm. The form of Oseen equations examined here is obtained by linearizing the incompressible Navier–Stokes equations. An algorithm is presented for approximately solving steady‐state, incompressible Navier–Stokes equations with a nested iteration‐Newton‐FOSLS‐AMG iterative scheme, which involves solving a sequence of Oseen equations. Some numerical results for Kovasznay flow are provided. Copyright © 2006 John Wiley & Sons, Ltd.     

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     

Shape reconstruction for the time‐dependent Navier‐Stokes flow

This article deals with the shape reconstruction of a bounded domain with a viscous incompressible fluid driven by the time‐dependent Navier‐Stokes equations. For the approximate solution of the ill‐posed and nonlinear problem we propose a regularized Newton method. A theoretical foundation for the Newton method is given by establishing the differentiability of the initial boundary value problem with respect to the interior boundary curve in the sense of the domain derivative. Numerical examples indicate the feasibility of our method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

On the stabilizing effect of convection in three‐dimensional incompressible flows

We investigate the stabilizing effect of convection in three‐dimensional incompressible Euler and Navier‐Stokes equations. The convection term is the main source of nonlinearity for these equations. It is often considered destabilizing although it conserves energy due to the incompressibility condition. In this paper, we show that the convection term together with the incompressibility condition actually has a surprising stabilizing effect. We demonstrate this by constructing a new three‐dimensional model that is derived for axisymmetric flows with swirl using a set of new variables. This model preserves almost all the properties of the full three‐dimensional Euler or Navier‐Stokes equations except for the convection term, which is neglected in our model. If we added the convection term back to our model, we would recover the full Navier‐Stokes equations. We will present numerical evidence that seems to support that the three‐dimensional model may develop a potential finite time singularity. We will also analyze the mechanism that leads to these singular events in the new three‐dimensional model and how the convection term in the full Euler and Navier‐Stokes equations destroys such a mechanism, thus preventing the singularity from forming in a finite time. © 2008 Wiley Periodicals, Inc.     

The inviscid limit of the incompressible anisotropic Navier–Stokes equations with the non-slip boundary condition

In this paper, we study the asymptotic behavior for the incompressible anisotropic Navier–Stokes equations with the non-slip boundary condition in a half space of ${\mathbb{R}^3}$ when the vertical viscosity goes to zero. Firstly, by multi-scale analysis, we formally deduce an asymptotic expansion of the solution to the problem with respect to the vertical viscosity, which shows that the boundary layer appears in the tangential velocity field and satisfies a nonlinear parabolic–elliptic coupled system. Also from the expansion, it is observed that away from the boundary the solution of the anisotropic Navier–Stokes equations formally converges to a solution of a degenerate incompressible Navier–Stokes equation. Secondly, we study the well-posedness of the problems for the boundary layer equations and then rigorously justify the asymptotic expansion by using the energy method. We obtain the convergence results of the vanishing vertical viscosity limit, that is, the solution to the incompressible anisotropic Navier–Stokes equations tends to the solution to degenerate incompressible Navier–Stokes equations away from the boundary, while near the boundary, it tends to the boundary layer profile, in both the energy space and the L ^∞ space.     

On the Inviscid Limit Problem of the Vorticity Equations for Viscous Incompressible Flows in the Half‐Plane

We consider the Navier‐Stokes equations for viscous incompressible flows in the half‐plane under the no‐slip boundary condition. By using the vorticity formulation we prove the local‐in‐time convergence of the Navier‐Stokes flows to the Euler flows outside a boundary layer and to the Prandtl flows in the boundary layer in the inviscid limit when the initial vorticity is located away from the boundary. © 2014 Wiley Periodicals, Inc.     

Analysis of the Yosida method for the incompressible Navier–Stokes equations

The Yosida method was introduced in (Quarteroni et al., to appear) for the numerical approximation of the incompressible unsteady Navier–Stokes equations. From the algebraic viewpoint, it can be regarded as an inexact factorization of the matrix arising from the space and time discretization of the problem. However, its differential interpretation resides on an elliptic stabilization of the continuity equation through the Yosida regularization of the Laplacian (see (Brezis, 1983, Ciarlet and Lions, 1991)). The motivation of this method as well as an extensive numerical validation were given in (Quarteroni et al., to appear).In this paper we carry out the analysis of this scheme. In particular, we consider a first-order time advancing unsplit method. In the case of the Stokes problem, we prove unconditional stability and moreover that the splitting error introduced by the Yosida scheme does not affect the overall accuracy of the solution, which remains linear with respect to the time step. Some numerical experiments, for both the Stokes and Navier–Stokes equations, are presented in order to substantiate our theoretical results.     

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     

Numerical analysis of a Picard multilevel stabilization of mixed finite volume method for the 2D/3D incompressible flow with large data

In this article, we develop a branch of nonsingular solutions of a Picard multilevel stabilization of mixed finite volume method for the 2D/3D stationary Navier‐Stokes equations without relying on the unique solution condition. The method presented consists of capturing almost all information of initial problem (the nonlinear problems) on the coarsest mesh and then performs one Picard defect correction (the linear problems) on each subsequent mesh based on previous information thus only solving one large linear systems. What is more, the method presented can results in a better coefficient matrix in the model presented with small viscosity. Theoretical results show that the method presented is derived with the convergence rate of the same order as the corresponding finite volume method/finite element method solving the stationary Navier‐Stokes equations on a fine mesh. Therefore, the method presented is definitely more efficient than the standard finite volume method/finite element method. Finally, numerical experiments clearly show the efficiency of the method presented for solving the stationary Navier‐Stokes equations.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 34: 30–50, 2018     

Global large solutions to incompressible Navier‐Stokes equations with gravity

In this per, we consider a special class of initial data for the three‐dimensional incompressible Navier–Stokes equations with gravity. We show that, under such conditions, the incompressible Navier‐Stokes equations with gravity are globally well posed, and the velocity minus gravity term has finite energy. The important features of the initial data is that the velocity fields minus gravity term are almost parallel to the corresponding vorticity fields in a very large space domain. Copyright © 2014 John Wiley & Sons, Ltd.