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     

A two‐level method in space and time for the Navier‐Stokes equations

A two‐level method in space and time for the time‐dependent Navier‐Stokes equations is considered in this article. The approximate solution u_M∈H_M is decomposed into the large eddy component v∈H_m(m < M) and the small eddy component w∈H. We obtain the large eddy component v by solving a standard Galerkin equation in a coarse‐level subspace H_m with a time step length k, whereas the small eddy component w is derived by solving a linear equation in an orthogonal complement subspace H with a time step length pk, where p is a positive integer. The analysis shows that our two‐level scheme has long‐time stability and can reach the same accuracy as the standard Galerkin method in fine‐level subspace H_M for an appropriate configuration of p and m. Moreover, some numerical examples are provided to complement our theoretical analysis. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Discretization of the Navier‐Stokes equations with slip boundary condition

We propose and analyze a two‐level method of discretizing the nonlinear Navier‐Stokes equations with slip boundary condition. The slip boundary condition is appropriate for problems that involve free boundaries, flows past chemically reacting walls, and other examples where the usual no‐slip condition u = 0 is not valid. The two‐level algorithm consists of solving a small nonlinear system of equations on the coarse mesh and then using that solution to solve a larger linear system on the fine mesh. The two‐level method exploits the quadratic nonlinearity in the Navier‐Stokes equations. Our error estimates show that it has optimal order accuracy, provided that the best approximation to the true solution in the velocity and pressure spaces is bounded above by the data. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 26–42, 2001     

A two‐grid stabilization method for solving the steady‐state Navier‐Stokes equations

We formulate a subgrid eddy viscosity method for solving the steady‐state incompressible flow problem. The eddy viscosity does not act on the large flow structures. Optimal error estimates are obtained for velocity and pressure. The numerical illustrations agree completely with the theoretical results. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

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.     

A defect correction method for the time‐dependent Navier‐Stokes equations

A method for solving the time dependent Navier‐Stokes equations, aiming at higher Reynolds' number, is presented. The direct numerical simulation of flows with high Reynolds' number is computationally expensive. The method presented is unconditionally stable, computationally cheap, and gives an accurate approximation to the quantities sought. In the defect step, the artificial viscosity parameter is added to the inverse Reynolds number as a stability factor, and the system is antidiffused in the correction step. Stability of the method is proven, and the error estimations for velocity and pressure are derived for the one‐ and two‐step defect‐correction methods. The spacial error is O(h) for the one‐step defect‐correction method, and O(h²) for the two‐step method, where h is the diameter of the mesh. The method is compared to an alternative approach, and both methods are applied to a singularly perturbed convection–diffusion problem. The numerical results are given, which demonstrate the advantage (stability, no oscillations) of the method presented. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

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.     

Semi‐implicit schemes for transient Navier–Stokes equations and eddy viscosity models

This study presents two computational schemes for the numerical approximation of solutions to eddy viscosity models as well as transient Navier–Stokes equations. The eddy viscosity model is one example of a class of Large Eddy Simulation models, which are used to simulate turbulent flow. The first approximation scheme is a first order single step method that treats the nonlinear term using a semi‐implicit discretization. The second scheme employs a two step approach that applies a Crank–Nicolson method for the nonlinear term while also retaining the semi‐implicit treatment used in the first scheme. A finite element approximation is used in the spatial discretization of the partial differential equations. The convergence analysis for both schemes is discussed in detail, and numerical results are given for two test problems one of which is the two dimensional flow around a cylinder. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A two‐level finite element method for the stationary Navier‐Stokes equations based on a stabilized local projection

In this article, we propose a two‐level finite element method to analyze the approximate solutions of the stationary Navier‐Stokes equations based on a stabilized local projection. The local projection allows to circumvent the Babuska‐Brezzi condition by using equal‐order finite element pairs. The local projection can be used to stabilize high equal‐order finite element pairs. The proposed method combines the local projection stabilization method and the two‐level method under the assumption of the uniqueness condition. The two‐level method consists of solving a nonlinear equation on the coarse mesh and solving a linear equation on fine mesh. The nonlinear equation is solved by the one‐step Newtonian iteration method. In the rest of this article, we show the error analysis of the lowest equal‐order finite element pair and provide convergence rate of approximate solutions. Furthermore, the numerical illustrations coincide with the theoretical analysis expectations. From the view of computational time, the results show that the two‐level method is effective to solve the stationary Navier‐Stokes equations. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

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.     

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     

Remark on stability of rarefaction waves to the one‐dimensional compressible Navier–Stokes equations with density‐dependent viscosity coefficient

In this paper, we study one‐dimensional compressible isentropic Navier–Stokes equations with density‐dependent viscosity. We can obtain the asymptotic stability of rarefaction waves for the compressible isentropic Navier–Stokes equations when the power of viscosity coefficient , which enlarge the range of α in the article [Jiu Q, Wang Y, Xin ZP, Communication in Partial Differential Equations 2011; 36: 602‐634]. Copyright © 2013 John Wiley & Sons, Ltd.     

Global C∞‐solutions to 1D compressible Navier–Stokes equations with density‐dependent viscosity

In this paper, we consider the existence of global smooth solutions to 1D compressible isentropic Navier–Stokes equations with density‐dependent viscosity and free boundaries. The initial density ρ₀∈W^1,2n is bounded below away from zero and the initial velocity u₀∈L²ⁿ. The viscosity coefficient µ is proportional to ρ^θ with 0<θ?1, where ρis the density. The existence and uniqueness of global solutions in Hⁱ([0,1])(i = 1,2,4) have been established in (J. Math. Phys. 2009; 50 :023101; Meth. Appl. Anal. 2005; 12 :239–252; J. Differ. Equations 2008; 245:3956–3973; Commun. Pure Appl. Anal. 2008; 7 :373–381). By mathematical induction method, we will establish the existence of global smooth solutions to 1D compressible isentropic Navier–Stokes equations with density‐dependent viscosity and free boundaries when the initial data ρ₀ and u₀ are smooth. Copyright © 2011 John Wiley & Sons, Ltd.     

Two‐grid variational multiscale algorithms for the stationary incompressible Navier‐Stokes equations with friction boundary conditions

Two‐grid variational multiscale (VMS) algorithms for the incompressible Navier‐Stokes equations with friction boundary conditions are presented in this article. First, one‐grid VMS algorithm is used to solve this problem and some error estimates are derived. Then, two‐grid VMS algorithms are proposed and analyzed. The algorithms consist of nonlinear problem on coarse grid and linearized problem (Stokes problem or Oseen problem) on fine grid. Moreover, the stability and convergence of the present algorithms are established. Finally, Numerical results are shown to confirm the theoretical analysis. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 546–569, 2017     

Two‐dimensional stationary Navier–Stokes equations with 4‐cyclic symmetry

This paper is concerned with the stationary Navier–Stokes equation in the whole plane and in the two–dimensional exterior domain invariant under the action of the cyclic group of order 4, and gives a condition on the potentials yielding the external force, and on the boundary value, sufficient for the unique existence of a small solution equivariant with respect to the aforementioned cyclic group.     

Semi‐invariant solutions of the Navier‐Stokes equations

We present new exact solutions and reduced differential systems of the Navier‐Stokes equations of incompressible viscous fluid flow. We apply the method of semi‐invariant manifolds, introduced earlier as a modification of the Lie invariance method. We show that many known solutions of the Navier‐Stokes equations are, in fact, semi‐invariant and that the reduced differential systems we derive using semi‐invariant manifolds generalize previously obtained results that used ad hoc methods. Many of our semi‐invariant solutions solve decoupled systems in triangular form that are effectively linear. We also obtain several new reductions of Navier‐Stokes to a single nonlinear partial differential equation. In some cases, we can solve reduced systems and generate new analytic solutions of the Navier‐Stokes equations or find their approximations, and physical interpretation.     

A Liouville theorem for the compressible Navier‐Stokes equations

This paper is concerned with 3‐dimensional steady compressible Navier‐Stokes equations. A Liouville‐type theorem is proved when some suitable conditions are satisfied.     

Second‐order Galerkin‐Lagrange method for the Navier‐Stokes equations (retracted article)

It has come to the attention of the editors and publisher that an article published in Numerical Methods and Partial Differential Equations, “Second‐order Galerkin‐Lagrange method for the Navier‐Stokes equations,” by Mohamed Bensaada, Driss Esselaoui, and Pierre Saramito, Numer Methods Partial Differential Eq 21(6) (2005), 1099–1121 included large portions that were copied from the following paper without proper citation: “Convergence and nonlinear stability of the Lagrange‐Galerkin method for the Navier‐Stokes equations,” Endre Suli, Numerische Mathematik, Vol. 53, No. 4, pp. 459–486 (July, 1988). We have retracted the paper and apologize to Dr. Suli Numer Methods Partial Differential Eq (2007)23(1)211 .     

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.