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     

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.     

A posteriori error estimators for the steady incompressible Navier–Stokes equations

A residual-based a posteriori error estimator for finite element discretizations of the steady incompressible Navier–Stokes equations in the primitive variable formulation is discussed. Though the estimator is similar to existing ones, an alternate derivation is presented, involving an abstract estimate that may prove of some intrinsic value. The estimator is particularized to Hood–Taylor and modified Hood–Taylor finite elements and proved to be a global upper bound (up to a positive multiplicative constant) of the true error. Numerical examples are provided to illustrate the performance of the resulting mesh adaptation process. © 1997 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 561–574, 1997     

An iterative solver for the Oseen problem and numerical solution of incompressible Navier–Stokes equations

Incompressible unsteady Navier–Stokes equations in pressure–velocity variables are considered. By use of the implicit and semi‐implicit schemes presented the resulting system of linear equations can be solved by a robust and efficient iterative method. This iterative solver is constructed for the system of linearized Navier–Stokes equations. The Schur complement technique is used. We present a new approach of building a non‐symmetric preconditioner to solve a non‐symmetric problem of convection–diffusion and saddle‐point type. It is shown that handling the differential equations properly results in constructing efficient solvers for the corresponding finite linear algebra systems. The method has good performance for various ranges of viscosity and can be used both for 2D and 3D problems. The analysis of the method is still partly heuristic, however, the mathematically rigorous results are proved for certain cases. The proof is based on energy estimates and basic properties of the underlying partial differential equations. Numerical results are provided. Additionally, a multigrid method for the auxiliary convection–diffusion problem is briefly discussed. Copyright © 1999 John Wiley & Sons, Ltd.     

Coupling method for the exterior stationary Navier–Stokes equations

We represent a new numerical method to solve the stationary Navier–Stokes equations in an unbounded domain. This technique consists in coupling the boundary integral and finite element methods. Moreover, we derive the variational formulation and well-posedness of the coupling method and provide the convergence result for the approximate solution. © 1993 John Wiley & Sons, Inc.     

Two-level Galerkin–Lagrange multipliers method for the stationary Navier–Stokes equations

In this paper, we combine the Galerkin–Lagrange multiplier (GLM) method with the two-level method to solve the stationary Navier–Stokes equations in order to avoid the time-consuming process and the construction of zero-divergence elements. Different quadrilateral partitions are used for approximating the velocity and the pressure. Then some error estimates are obtained and some numerical results of the GLM method and the two-level GLM method are given. The results show that the two-level method based on the GLM method is more efficient than the GLM method under the convergence rate of same order.     

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     

On the asymptotic stability of steady solutions of the Navier–Stokes equations in unbounded domains

We consider the problem of the asymptotic behaviour in the L²‐norm of solutions of the Navier–Stokes equations. We consider perturbations to the rest state and to stationary motions. In both cases we study the initial‐boundary value problem in unbounded domains with non‐compact boundary. In particular, we deal with domains with varying and possibly divergent exits to infinity and aperture domains. Copyright © 2007 John Wiley & Sons, Ltd.     

A two-grid method based on Newton iteration for the Navier–Stokes equations

In this paper, we consider a two-grid method for resolving the nonlinearity in finite element approximations of the equilibrium Navier–Stokes equations. We prove the convergence rate of the approximation obtained by this method. The two-grid method involves solving one small, nonlinear coarse mesh system and two linear problems on the fine mesh which have the same stiffness matrix with only different right-hand side. The algorithm we study produces an approximate solution with the optimal asymptotic in h and accuracy for any Reynolds number. Numerical example is given to show the convergence of the method.     

A superconvergent nonconforming mixed finite element method for the Navier–Stokes equations

The superconvergence for a nonconforming mixed finite element approximation of the Navier–Stokes equations is analyzed in this article. The velocity field is approximated by the constrained nonconforming rotated Q₁ (CNRQ₁) element, and the pressure is approximated by the piecewise constant functions. Under some regularity assumptions, the superconvergence estimates for both the velocity in broken H¹‐norm and the pressure in L²‐norm are obtained. Some numerical examples are presented to demonstrate our theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 646–660, 2016     

A nearly optimal preconditioner for the Navier–Stokes equations

We present a preconditioner for the linearized Navier–Stokes equations which is based on the combination of a fast transform approximation of an advection diffusion problem together with the recently introduced ‘BFBT^T’ preconditioner of Elman (SIAM Journal of Scientific Computing, 1999; 20 :1299–1316). The resulting preconditioner when combined with an appropriate Krylov subspace iteration method yields the solution in a number of iterations which appears to be independent of the Reynolds number provided a mesh Péclet number restriction holds, and depends only mildly on the mesh size. The preconditioner is particularly appropriate for problems involving a primary flow direction. Copyright © 2001 John Wiley & Sons, Ltd.     

Stabilized spectral element approximation for the Navier–Stokes equations

The conforming spectral element methods are applied to solve the linearized Navier–Stokes equations by the help of stabilization techniques like those applied for finite elements. The stability and convergence analysis is carried out and essential numerical results are presented demonstrating the high accuracy of the method as well as its robustness. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 115–141, 1998     

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.     

Energy stable numerical scheme for the viscous Cahn–Hilliard–Navier–Stokes equations with moving contact line

In the present article we study the numerics of the viscous Cahn–Hilliard–Navier–Stokes model, endowed with dynamic boundary conditions which allow us to take into account the interaction between the fluids interface and the moving walls of the physical domain. In what follows, we propose an energy stable temporal scheme for the problem and we prove the stability and the unconditional solvability of the discretization proposed. We also propose a fully discrete scheme for which we prove the stability and the unconditional solvability. Numerical simulations are presented to illustrate the theoretical results.     

A fully discrete nonlinear Galerkin method for the 3D Navier–Stokes equations

The purpose of this paper is twofold: (i) We show that the Fourier‐based Nonlinear Galerkin Method (NLGM) constructs suitable weak solutions to the periodic Navier–Stokes equations in three space dimensions provided the large scale/small scale cutoff is appropriately chosen. (ii) If smoothness is assumed, NLGM always outperforms the Galerkin method by a factor equal to 1 in the convergence order of the H ¹‐norm for the velocity and the L²‐norm for the pressure. This is a purely linear superconvergence effect resulting from standard elliptic regularity and holds independently of the nature of the boundary conditions (whether periodicity or no‐slip BC is enforced). © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A sparse grid method for the Navier–Stokes equations based on hyperbolic cross

A sparse grid method for the time‐dependent Navier–Stokes equations based on hyperbolic cross approximation is considered in this article. Subsequent truncation of the associated series expansion results in a sparse grid discretization. Stability and convergence of the fully discrete sparse grid method are established. Finally, the numerical experiment is presented to show the effectiveness of this sparse grid method. Copyright © 2013 John Wiley & Sons, Ltd.     

Compressible Navier–Stokes equations with vacuum state in the case of general pressure law

In this paper, we consider the one‐dimensional compressible isentropic Navier–Stokes equations with a general ‘pressure law’ and the density‐dependent viscosity coefficient when the density connects to vacuum continuously. Precisely, the viscosity coefficient µ is proportional to ρ^θ and 0<θ<1, where ρ is the density. And the pressure P = P(ρ) is a general ‘pressure law’. The global existence and the uniqueness of weak solutions is proved, and a decay result for the pressure as t→ + ∞ is given. It is also proved that no vacuum states and no concentration states develop, and the free boundary do not expand to infinite. Copyright © 2006 John Wiley & Sons, Ltd.     

Parabolic Dirac operators and the Navier–Stokes equations over time‐varying domains

We consider parabolic Dirac operators which do not involve fractional derivatives and use them to show the solvability of the in‐stationary Navier–Stokes equations over time‐varying domains. Copyright © 2005 John Wiley & Sons, Ltd.     

Logarithmically improved regularity criteria for Navier–Stokes and related equations

We use an interpolation inequality on Besov spaces to show some logarithmically improved regularity criteria for Navier–Stokes equations, the harmonic heat flow, the Landau–Lifshitz equations, and the Landau–Lifshitz–Maxwell system. Copyright © 2009 John Wiley & Sons, Ltd.