A posteriori error estimations for mixed finite-element approximations to the Navier–Stokes equations

A posteriori estimates for mixed finite element discretizations of the Navier–Stokes equations are derived. We show that the task of estimating the error in the evolutionary Navier–Stokes equations can be reduced to the estimation of the error in a steady Stokes problem. As a consequence, any available procedure to estimate the error in a Stokes problem can be used to estimate the error in the nonlinear evolutionary problem. A practical procedure to estimate the error based on the so-called postprocessed approximation is also considered. Both the semidiscrete (in space) and the fully discrete cases are analyzed. Some numerical experiments are provided.     

Two-grid stabilized algorithms for the steady Navier–Stokes equations with damping

This paper presents and studies three two-grid stabilized quadratic equal-order finite element algorithms based on two local Gauss integrations for the steady Navier–Stokes equations with damping. In these algorithms, we first solve a stabilized nonlinear problem on a coarse grid, and then pass the coarse grid solution to a fine grid and solve a stabilized linear problem. Using some nonlinear analysis techniques, we analyze stability of the algorithms and derive optimal order error estimates of the approximate solutions. Theoretical and numerical results show that, when the algorithmic parameters are chosen appropriately, the accuracy of the approximate solutions computed by our two-grid stabilized algorithms is comparable to that of solving a fully stabilized nonlinear problem on the same fine grid; however, our two-grid algorithms save a large amount of CPU time than the one-grid stabilized algorithm.     

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.     

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     

Decoupled modified characteristics finite element method for the time dependent Navier–Stokes/Darcy problem

In this paper, a modified characteristics finite element method for the time dependent Navier–Stokes/Darcy problem with the Beavers–Joseph–Saffman interface condition is presented. In this method, the Navier–Stokes/Darcy equation is decoupled into two equations, one is the Navier–Stokes equation, the other is the Darcy equation, and the Navier–Stokes equation is solved by the modified characteristics finite element method. The theory analysis shows that this method has a good convergence property. In order to show the effect of our method, some numerical results was presented. The numerical results show that this method is highly efficient. Copyright © 2013 John Wiley & Sons, Ltd.     

On the semi-discrete stabilized finite volume method for the transient Navier–Stokes equations

A stabilized finite volume method for solving the transient Navier–Stokes equations is developed and studied in this paper. This method maintains conservation property associated with the Navier–Stokes equations. An error analysis based on the variational formulation of the corresponding finite volume method is first introduced to obtain optimal error estimates for velocity and pressure. This error analysis shows that the present stabilized finite volume method provides an approximate solution with the same convergence rate as that provided by the stabilized linear finite element method for the Navier–Stokes equations under the same regularity assumption on the exact solution and a slightly additional regularity on the source term. The stability and convergence results of the proposed method are also demonstrated by the numerical experiments presented.     

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.     

Two‐level Newton iterative method for the 2D/3D steady Navier‐Stokes equations

A combination method of the Newton iteration and two‐level finite element algorithm is applied for solving numerically the steady Navier‐Stokes equations under the strong uniqueness condition. This algorithm is motivated by applying the m Newton iterations for solving the Navier‐Stokes problem on a coarse grid and computing the Stokes problem on a fine grid. Then, the uniform stability and convergence with respect to ν of the two‐level Newton iterative solution are analyzed for the large m and small H and h << H. Finally, some numerical tests are made to demonstrate the effectiveness of the method. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

Effective Chorin–Temam algebraic splitting schemes for the steady Navier–stokes equations

This paper continues some recent work on the numerical solution of the steady incompressible Navier–Stokes equations. We present a new method, similar to the one presented in Rebholz et al., but with superior convergence and numerical properties. The method is efficient as it allows one to solve the same symmetric positive‐definite system for the pressure at each iteration, allowing for the simple preconditioning and the reuse of preconditioners. We also demonstrate how one can replace the Schur complement system with a diagonal matrix inversion while maintaining accuracy and convergence, at a small fraction of the numerical cost. Convergence is analyzed for Newton and Picard‐type algorithms, as well as for the Schur complement approximation.     

Decoupled characteristic stabilized finite element method for time‐dependent Navier–Stokes/Darcy model

In this article, we propose and analyze a new decoupled characteristic stabilized finite element method for the time‐dependent Navier–Stokes/Darcy model. The key idea lies in combining the characteristic method with the stabilized finite element method to solve the decoupled model by using the lowest‐order conforming finite element space. In this method, the original model is divided into two parts: one is the nonstationary Navier–Stokes equation, and the other one is the Darcy equation. To deal with the difficulty caused by the trilinear term with nonzero boundary condition, we use the characteristic method. Furthermore, as the lowest‐order finite element pair do not satisfy LBB (Ladyzhen‐Skaya‐Brezzi‐Babuska) condition, we adopt the stabilized technique to overcome this flaw. The stability of the numerical method is first proved, and the optimal error estimates are established. Finally, extensive numerical results are provided to justify the theoretical analysis.     

Improving the pressure accuracy in a projection scheme for incompressible fluids with variable viscosity

The incremental projection scheme and its enhanced version, the rotational projection scheme are powerful and commonly used approaches producing efficient numerical algorithms for solving the Navier–Stokes equations. However, the much improved rotational projection scheme cannot be used on models with non-homogeneous viscosity, imposing the use of the less accurate incremental projection. This paper presents a projection method for the Navier–Stokes equations for fluids having variable viscosity, giving a consistent pressure and increased accuracy in pressure when compared to the incremental projection. The accuracy of the method will be illustrated using a manufactured solution.     

Two-level defect-correction locally stabilized finite element method for the steady Navier–Stokes equations

This paper proposes a two-level defect-correction stabilized finite element method for the steady Navier–Stokes equations based on local Gauss integration. The method combines the two-level strategy with the defect-correction method under the assumption of the uniqueness condition. Both the simplified and the Newton scheme are proposed and analyzed. Moreover, the numerical illustrations agree completely with the theoretical expectations.     

Variational iteration method for solving three‐dimensional Navier–Stokes equations of flow between two stretchable disks

The similarity transform for the steady three‐dimensional Navier–Stokes equations of flow between two stretchable disks gives a system of nonlinear ordinary differential equations. In this article, the variational iteration method was used for solving these equations. The results have been compared with the numerical results. This article depicts that the VIM is an efficient and powerful method for solving nonlinear differential equations. This method is applicable to strongly and weakly nonlinear problems. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Large time behavior for the non‐isentropic Navier–Stokes–Maxwell system

In this paper, we are concerned with the system of the non‐isentropic compressible Navier–Stokes equations coupled with the Maxwell equations through the Lorentz force in three space dimensions. The global existence of solutions near constant steady states is established, and the time‐decay rates of perturbed solutions are obtained. The proof for existence is due to the classical energy method, and the investigation of large‐time behavior is based on the linearized analysis of the non‐isentropic Navier–Stokes–Poisson equations and the electromagnetic part for the linearized isentropic Navier–Stokes–Maxwell equations. In the meantime, the time‐decay rates obtained by Zhang, Li, and Zhu [J. Differential Equations, 250(2011), 866‐891] for the linearized non‐isentropic Navier–Stokes–Poisson equations are improved. Copyright © 2016 John Wiley & Sons, Ltd.     

Stability and convergence of optimum spectral non‐linear Galerkin methods

Our objective in this article is to present some numerical schemes for the approximation of the 2‐D Navier–Stokes equations with periodic boundary conditions, and to study the stability and convergence of the schemes. Spatial discretization can be performed by either the spectral Galerkin method or the optimum spectral non‐linear Galerkin method; time discretization is done by the Euler scheme and a two‐step scheme. Our results show that under the same convergence rate the optimum spectral non‐linear Galerkin method is superior to the usual Galerkin methods. Finally, numerical example is provided and supports our results. Copyright © 2001 John Wiley & Sons, Ltd.     

Preconditioning for boundary control problems in incompressible fluid dynamics

PDE‐constrained optimization problems arise in many physical applications, prominently in incompressible fluid dynamics. In recent research, efficient solvers for optimization problems governed by the Stokes and Navier–Stokes equations have been developed, which are mostly designed for distributed control. Our work closes a gap by showing the effectiveness of an appropriately modified preconditioner to the case of Stokes boundary control. We also discuss the applicability of an analogous preconditioner for Navier–Stokes boundary control and provide some numerical results.     

The coupling system of Navier–Stokes equations and elastic Navier–Lame equations in a blood vessel

In this paper, the blood flow problem is considered in a blood vessel, and a coupling system of Navier–Stokes equations and linear elastic equations, Navier–Lame equations, in a cylinder with cylindrical elastic shell is given as the governing equations of the problem. We provide two finite element models to simulating the three-dimensional Navier–Stokes equations in the cylinder while the asymptotic expansion method is used to solving the linearly elastic shell equations. Specifically, in order to discrete the Navier–Stokes equations, the dimensional splitting strategy is constructed under the cylinder coordinate system. The spectral method is adopted along the rotation direction while the finite element method is used along the other directions. By using the above strategy, we get a series of two-dimensional-three-components (2D-3C) fluid problems. By introduce the S-coordinate system in E³ and employ the thickness of blood vessel wall as the expanding parameter, the asymptotic expansion method can be established to approximate the solution of the 3D elastic problem. The interface contact conditions can be treated exactly based on the knowledge of tensor analysis. Finally, numerical test shows that our method is reasonable.     

Spectral method for vorticity‐stream function form of Navier–Stokes equations with slip boundary conditions

In this paper, we propose a spectral method for the vorticity‐stream function form of the Navier–Stokes equations with slip boundary conditions. The numerical solutions fulfill the incompressibility and the physical boundary conditions automatically. The stability and convergence of the proposed methods are proven. Numeric results demonstrate the efficiency of suggested algorithm. Copyright © 2011 John Wiley & Sons, Ltd.     

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.