Bivariate spline method for numerical solution of steady state Navier–Stokes equations over polygons in stream function formulation

We use the bivariate spline finite elements to numerically solve the steady state Navier–Stokes equations. The bivariate spline finite element space we use in this article is the space of splines of smoothness r and degree 3r over triangulated quadrangulations. The stream function formulation for the steady state Navier–Stokes equations is employed. Galerkin's method is applied to the resulting nonlinear fourth‐order equation, 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 H²(Ω) of the nonlinear fourth‐order problem and give an estimate of how fast the numerical solution converges to the weak solution. The Galerkin method with C¹ cubic splines is implemented in MATLAB. Our numerical experiments show that the method is effective and efficient. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 147–183, 2000     

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 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     

Study of the mixed finite volume method for Stokes and Navier‐Stokes equations

We present finite volume schemes for Stokes and Navier‐Stokes equations. These schemes are based on the mixed finite volume introduced in (Droniou and Eymard, Numer Math 105 (2006), 35‐71), and can be applied to any type of grid (without “orthogonality” assumptions as for classical finite volume methods) and in any space dimension. We present numerical results on some irregular grids, and we prove, for both Stokes and Navier‐Stokes equations, the convergence of the scheme toward a solution of the continuous problem. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A block preconditioning technique for the streamfunction‐vorticity formulation of the Navier‐Stokes equations

Iterative methods of Krylov‐subspace type can be very effective solvers for matrix systems resulting from partial differential equations if appropriate preconditioning is employed. We describe and test block preconditioners based on a Schur complement approximation which uses a multigrid method for finite element approximations of the linearized incompressible Navier‐Stokes equations in streamfunction and vorticity formulation. By using a Picard iteration, we use this technology to solve fully nonlinear Navier‐Stokes problems. The solvers which result scale very well with problem parameters. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

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     

Clifford analysis and the Navier‐Stokes equations over unbounded domains

We describe a finite difference scheme based on a discrete function theory for the Navier‐Stokes equations over unbounded domains which has only low regularity demands on the solution for convergence. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Global mild solutions of Navier‐Stokes equations

We establish a global well‐posedness of mild solutions to the three‐dimensional, incompressible Navier‐Stokes equations if the initial data are in the space ${\cal{X}}^{-1}$ defined by \input amssym ${\cal{X}}^{‐1} = \{f \in {\cal{D}}^\prime(R^3): \int_{{\Bbb{R}}^3}|\xi|^{‐1}|\widehat{f}|d\xi < \infty\}$ and if the norms of the initial data in ${\cal{X}}^{-1}$ are bounded exactly by the viscosity coefficient μ. © 2010 Wiley Periodicals, Inc.     

Uniqueness of the weak solution to the fractional anisotropic Navier‐Stokes equations

In this work, we demonstrate uniqueness of the weak solution to the fractional anisotropic Navier‐Stokes system with only horizontal dissipation.     

Block factorized preconditioners for high‐order accurate in time approximation of the Navier‐Stokes equations

Computationally efficient solution methods for the unsteady Navier‐Stokes incompressible equations are mandatory in real applications of fluid dynamics. A typical strategy to reduce the computational cost is to split the original problem into subproblems involving the separate computation of velocity and pressure. The splitting can be carried out either at a differential level, like in the Chorin‐Temam scheme, or in an algebraic fashion, like in the algebraic reinterpretation of the Chorin‐Temam method, or in the Yosida scheme (see 1 and 19 ). These fractional step schemes indeed provide effective methods of solution when dealing with first order accurate time discretizations. Their extension to high order time discretization schemes is not trivial. To this end, in the present work we focus our attention on the adoption of inexact algebraic factorizations as preconditioners of the original problem. We investigate their properties and show that some particular choices of the approximate factorization lead to very effective schemes. In particular, we prove that performing a small number of preconditioned iterations is enough to obtain a time accurate solution, irrespective of the dimension of the system at hand. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 487–510, 2003     

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.     

A multilevel finite element method in space‐time for the Navier‐Stokes problem

A multilevel finite element method in space‐time for the two‐dimensional nonstationary Navier‐Stokes problem is considered. The method is a multi‐scale method in which the fully nonlinear Navier‐Stokes problem is only solved on a single coarsest space‐time mesh; subsequent approximations are generated on a succession of refined space‐time meshes by solving a linearized Navier‐Stokes problem about the solution on the previous level. The a priori estimates and error analysis are also presented for the J‐level finite element method. We demonstrate theoretically that for an appropriate choice of space and time mesh widths: h_j ～ h, k_j ～ k, j = 2, …, J, the J‐level finite element method in space‐time provides the same accuracy as the one‐level method in space‐time in which the fully nonlinear Navier‐Stokes problem is solved on a final finest space‐time mesh. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

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.     

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     

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 estimate and regularity for the compressible Navier‐Stokes equations by Newton's method

The finite element discretization error estimate and H¹ regularity are shown for the solution generated by Newton's method to the stationary compressible Navier‐Stokes equations by interpreting Newton's method as an equivalent iterative method. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 511–524, 2003     

Numerical simulations for two‐dimensional stochastic incompressible Navier‐Stokes equations

This article mainly concerns modeling the stochastic input and its propagation in incompressible Navier‐Stokes(N‐S) flow simulations. The stochastic input is represented spectrally by employing orthogonal polynomial functionals from the Askey scheme as trial basis to represent the random space. A standard Galerkin projection is applied in the random dimension to derive the equations in the weak form. The resulting set of deterministic equations is then solved with standard methods to obtain the mean solution. In this article, the main method employs the Hermite polynomial as the basis in random space. Numerical examples are given and the error analysis is demonstrated for a model problem. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010