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     

An analysis for the compressible Stokes equations by first‐order system of least‐squares finite element method

This article applies the first‐order system least‐squares (fosls) finite element method developed by Cai, Manteuffel and McCormick to the compressible Stokes equations. By introducing a new dependent velocity flux variable, we recast the compressible Stokes equations as a first‐order system. Then it is shown that the ellipticity and continuity hold for the least‐squares functionals employing the mixture of H^?1 and L², so that the fosls finite element methods yield best approximations for the velocity flux and velocity. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17:689–699, 2001     

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.     

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     

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.     

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     

Vanishing viscosity limit of the Navier‐Stokes equations to the euler equations for compressible fluid flow

We establish the vanishing viscosity limit of the Navier‐Stokes equations to the isentropic Euler equations for one‐dimensional compressible fluid flow. For the Navier‐Stokes equations, there exist no natural invariant regions for the equations with the real physical viscosity term so that the uniform sup‐norm of solutions with respect to the physical viscosity coefficient may not be directly controllable. Furthermore, convex entropy‐entropy flux pairs may not produce signed entropy dissipation measures. To overcome these difficulties, we first develop uniform energy‐type estimates with respect to the viscosity coefficient for solutions of the Navier‐Stokes equations and establish the existence of measure‐valued solutions of the isentropic Euler equations generated by the Navier‐Stokes equations. Based on the uniform energy‐type estimates and the features of the isentropic Euler equations, we establish that the entropy dissipation measures of the solutions of the Navier‐Stokes equations for weak entropy‐entropy flux pairs, generated by compactly supported C² test functions, are confined in a compact set in H^?1, which leads to the existence of measure‐valued solutions that are confined by the Tartar‐Murat commutator relation. A careful characterization of the unbounded support of the measure‐valued solution confined by the commutator relation yields the reduction of the measurevalued solution to a Dirac mass, which leads to the convergence of solutions of the Navier‐Stokes equations to a finite‐energy entropy solution of the isentropic Euler equations with finite‐energy initial data, relative to the different end‐states at infinity. © 2010 Wiley Periodicals, Inc.     

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.     

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     

Technical note: A fourth‐order finite difference scheme for a system of a 2D nonlinear elliptic partial differential equations

In this article we present a fourth‐order finite difference scheme, for a system of two‐dimensional, second‐order, nonlinear elliptic partial differential equations with mixed spatial derivative terms, using 13‐point stencils with a uniform mesh size h on a square region R subject to Dirichlet boundary conditions. The scheme of order h⁴ is derived using the local solution of the system on a single stencil. The resulting system of algebraic equations can be solved by iterative methods. The difference scheme can be easily modified to obtain formulae for grid points near the boundary. Computational results are given to demonstrate the performance of the scheme on some problems including Navier‐Stokes equations. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 43–53, 2001     

On the decay and stability of global solutions to the 3‐D inhomogeneous Navier‐Stokes equations

In this paper, we investigate the large‐time decay and stability to any given global smooth solutions of the 3‐D incompressible inhomogeneous Navier‐Stokes equations. In particular, we prove that given any global smooth solution (a,u) of (1.2), the velocity field u decays to 0 with an explicit rate, which coincides with the L² norm decay for the weak solutions of the 3‐D classical Navier‐Stokes system [26,29] as t goes to ∞. Moreover, a small perturbation to the initial data of (a,u) still generates a unique global smooth solution to (1.2), and this solution keeps close to the reference solution (a,u) for t > 0. We should point out that the main results in this paper work for large solutions of (1.2). © 2010 Wiley Periodicals, Inc.     

Bivariate spline method for numerical solution of time evolution Navier‐Stokes equations over polygons in stream function formulation

We use a bivariate spline method to solve the time evolution Navier‐Stokes equations numerically. The bivariate splines we use in this article are in the spline space of smoothness r and degree 3r over triangulated quadrangulations. The stream function formulation for the Navier‐Stokes equations is employed. Galerkin's method is applied to discretize the space variables of the nonlinear fourth‐order equation, Crank‐Nicholson's method is applied to discretize the time variable, and Newton's iterative method is then used to solve the resulting nonlinear system. We show the existence and uniqueness of the weak solution in L₂(0, T; H²(Ω)) ∩ L_∞(0, T; H¹(Ω)) of the 2D nonlinear fourth‐order problem and give an estimate of how fast the numerical solution converges to the weak solution. The C¹ cubic splines are implemented in MATLAB for solving the Navier‐Stokes equations numerically. Our numerical experiments show that the method is effective and efficient. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 776–827, 2003.     

Orthogonal spline collocation methods for the stream function‐vorticity formulation of the Navier–Stokes equations

An orthogonal spline collocation (OSC) spatial discretization is proposed for the solution of the fully coupled stream function‐vorticity formulation of the Navier–Stokes equations in two dimensions. For the time‐stepping, a three‐level leapfrog scheme is employed. This method is algebraically linear, and, at each time step, gives rise to a system of linear equations of the form arising in the OSC approximation of the biharmonic Dirichlet problem and can be solved by a fast direct method. Error estimates in the H^l–norm in space, l = 1,2, are derived for the semi‐discrete method and the fully‐discrete leapfrog scheme which is also shown to be second order accurate in time. Numerical results are presented which confirm the theoretical analysis and exhibit superconvergence phenomena, which provide superconvergent approximations to the components of the velocity. © John Wiley & Sons, Inc. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

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     

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.     

Local existence and blow‐up criterion for the generalized Boussinesq equations in Besov spaces

In this paper, we consider the three‐dimensional generalized Boussinesq equations, a system of equations resulting from replacing the Laplacian ? Δ in the usual Boussinesq equations by a fractional Laplacian ( ? Δ)^α. We prove the local existence in time and obtain a regularity criterion of solution for the generalized Boussinesq equations by means of the Littlewood–Paley theory and Bony's paradifferential calculus. The results in this paper can be regarded as an extension to the Serrin‐type criteria for Navier–Stokes equations and magnetohydrodynamics equations, respectively. Copyright © 2012 John Wiley & Sons, Ltd.     

Global Small Solutions to an MHD‐Type System: The Three‐Dimensional Case

In this paper, we consider the global well‐posedness of a three‐dimensional incompressible MHD type system with smooth initial data that is close to some nontrivial steady state. It is a coupled system between the Navier‐Stokes equations and a free transport equation with a universal nonlinear coupling structure. The main difficulty of the proof lies in exploring the dissipative mechanism of the system due to the fact that there is a free transport equation of ? in the coupled equations and only the horizontal derivatives of ? is dissipative with respect to time. To achieve this, we first employ anisotropic Littlewood‐Paley analysis to establish the key L¹(? ⁺ ; Lip(?³)) estimate to the third component of the velocity field. Then we prove the global well‐posedness to this system by the energy method, which depends crucially on the divergence‐free condition of the velocity field. © 2014 Wiley Periodicals, Inc.     

Stability of small‐amplitude shock profiles of symmetric hyperbolic‐parabolic systems

Combining pointwise Green's function bounds obtained in a companion paper [36] with earlier, spectral stability results obtained in [16], we establish nonlinear orbital stability of small‐amplitude Lax‐type viscous shock profiles for the class of dissipative symmetric hyperbolic‐parabolic systems identified by Kawashima [20], notably including compressible Navier‐Stokes equations and the equations of magnetohydrodynamics, obtaining sharp rates of decay in L^p with respect to small L¹ ∩ H³ perturbations, 2 ≤ p ≤ ∞. Our analysis extends and somewhat refines the approach introduced in [35] to treat stability of relaxation profiles. © 2004 Wiley Periodicals, Inc.     

A remark on the decay of solutions to the 3‐D Navier–Stokes equations

In this paper we derive a decay rate of the L²‐norm of the solution to the 3‐D Navier–Stokes equations. Although the result which is proved by Fourier splitting method is well known, our method is new, concise and direct. Moreover, it turns out that the new method established here has a wide application on other equations. Copyright © 2007 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