High-order accurate iterative solution of the Navier–Stokes equations for incompressible flows

An iterative solution scheme for the incompressible Navier–Stokes equations is presented. It is split into inner and outer iteration cycles, such that the momentum and continuity equations are satisfied within prescribed accuracy. The spatial discretization is based on high-order finite differences which makes it well suited for massively parallel computers. This is demonstrated in a scaling test. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Large time behavior for the incompressible Navier–Stokes flows in 2D exterior domains

The incompressible Navier–Stokes flows in a 2D exterior domain are considered, for which the associated total net force to the boundary may not vanish. The decay properties are shown for the first and second derivatives of the Navier–Stokes flows in L ¹ and weighted spaces, respectively, which improve Theorem 1.2 in Bae and Jin (J Funct Anal 240:508–529, 2006).     

Analysis of the Yosida method for the incompressible Navier–Stokes equations

The Yosida method was introduced in (Quarteroni et al., to appear) for the numerical approximation of the incompressible unsteady Navier–Stokes equations. From the algebraic viewpoint, it can be regarded as an inexact factorization of the matrix arising from the space and time discretization of the problem. However, its differential interpretation resides on an elliptic stabilization of the continuity equation through the Yosida regularization of the Laplacian (see (Brezis, 1983, Ciarlet and Lions, 1991)). The motivation of this method as well as an extensive numerical validation were given in (Quarteroni et al., to appear).In this paper we carry out the analysis of this scheme. In particular, we consider a first-order time advancing unsplit method. In the case of the Stokes problem, we prove unconditional stability and moreover that the splitting error introduced by the Yosida scheme does not affect the overall accuracy of the solution, which remains linear with respect to the time step. Some numerical experiments, for both the Stokes and Navier–Stokes equations, are presented in order to substantiate our theoretical results.     

On the well-posedness of 3-D inhomogeneous incompressible Navier–Stokes equations with variable viscosity

In this paper, we mainly study the well-posedness for the 3-D inhomogeneous incompressible Navier–Stokes equations with variable viscosity. With some smallness assumption on the BMO-norm of the initial density, we first get the local well-posedness of (1.1) in the critical Besov spaces. Moreover, if the viscosity coefficient is a constant, we can extend this local solution to be a global one. Our theorem implies that we have successfully extended the integrability index p of the initial velocity which has been obtained by Abidi, Gui and Zhang in [3], Burtea in [8] and Zhai and Yin in [32] to approach the ideal one i.e. $1<p<6$. The main novelty of this work is to apply the CRW theorem obtained by Coifman, Rochberg, Weiss in [11] to get a new a priori estimate for an elliptic equation with variable coefficients. The uniqueness of the solution also relies on a Lagrangian approach as in [16], [17], [18].     

Modified fully discretized projection method for the incompressible Navier–Stokes equations

The stability and convergence of a second-order fully discretized projection method for the incompressible Navier–Stokes equations is studied. In order to update the pressure field faster, modified fully discretized projection methods are proposed. It results in a nearly second-order method. This method sacrifices a little of accuracy, but it requires much less computations at each time step. It is very appropriate for actual computations. The comparison with other methods for the driven-cavity problem is presented.     

Algorithm for solving the Navier–Stokes equations for the modeling of creeping flows

We study a coupled algorithm for solving the two-dimensional Navier–Stokes equations in the stream function–vorticity variables. The algorithm is based on a finite-difference scheme in which the inertial terms in the vortex transport equation are taken from the lower time layer and the dissipative terms, from the upper time layer. In the linear approximation, we study the stability of this algorithm and use test computations to show its advantages when modeling flows at moderate Reynolds numbers.     

Strong solutions to the inhomogeneous Navier–Stokes–BGK system

In this paper, we are concerned with the local-in-time well-posedness of a fluid-kinetic model in which the BGK model with density dependent collision frequency is coupled with the inhomogeneous Navier–Stokes equation through drag forces. To the best knowledge of authors, this is the first result on the existence of local-in-time smooth solution for particle–fluid model with nonlinear inter-particle operator for which the existence of time can be prolonged as the size of initial data gets smaller.     

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.     

Parallel spectral-element direction splitting method for incompressible Navier–Stokes equations

An efficient parallel algorithm for the time dependent incompressible Navier–Stokes equations is developed in this paper. The time discretization is based on a direction splitting method which only requires solving a sequence of one-dimensional Poisson type equations at each time step. Then, a spectral-element method is used to approximate these one-dimensional problems. A Schur-complement approach is used to decouple the computation of interface nodes from that of interior nodes, allowing an efficient parallel implementation. The unconditional stability of the full discretized scheme is rigorously proved for the two-dimensional case. Numerical results are presented to show that this algorithm retains the same order of accuracy as a usual spectral-element projection type schemes but it is much more efficient, particularly on massively parallel computers.     

High order time discretisation methods for incompressible Navier–Stokes equations

In this paper we solve incompressible Navier–Stokes equations with higher order time integration schemes such as Radau-IIA and Lobatto-IIIC methods. We compare the results with 4th order DIRK and ROW methods to show the advantages of higher order time integration methods. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The incompressible Navier–Stokes limit of the Boltzmann equation for hard cutoff potentials

The present paper proves that all limit points of sequences of renormalized solutions of the Boltzmann equation in the limit of small, asymptotically equivalent Mach and Knudsen numbers are governed by Leray solutions of the Navier–Stokes equations. This convergence result holds for hard cutoff potentials in the sense of H. Grad, and therefore completes earlier results by the same authors [Invent. Math. 155 (2004) 81–161] for Maxwell molecules.     

On the Navier–Stokes flows for heat-conducting fluids with mixed boundary conditions

We consider a coupled model for steady flows of viscous incompressible heat-conducting fluids with temperature dependent material coefficients in a fixed three-dimensional open cylindrical channel. We introduce the Banach spaces X and Y to be the space of possible solutions of this problem and the space of its data, respectively. We show that the corresponding operator of the problem acting between X and Y is Fréchet differentiable. Applying the local diffeomorphism theorem we get the local solvability results for a variational formulation.     

Convergence analysis of an implicit fractional-step method for the incompressible Navier–Stokes equations

In this paper, an implicit fractional-step method for numerical solutions of the incompressible Navier–Stokes equations is studied. The time advancement is decomposed into a sequence of two steps, and the first step can be seen as a linear elliptic problem; on the other hand, the second step has the structure of the Stokes problem. The two problems satisfy the full homogeneous Dirichlet boundary conditions on the velocity. At the same time, we introduce a diffusion term −θΔu in all steps of the schemes. It allows to calculate by the large time step and enhance numerical stability by choosing the proper parameter values of θ. The convergence analysis and error estimates for the intermediate velocities, the end-of step velocities and the pressure solution are derived. Finally, numerical experiments show that the feasibility and effectiveness of this method.     

A parallel implementation of the modified augmented Lagrangian preconditioner for the incompressible Navier–Stokes equations

We describe a parallel implementation of a block triangular preconditioner based on the modified augmented Lagrangian approach to the steady incompressible Navier–Stokes equations. The equations are linearized by Picard iteration and discretized with various finite element and finite difference schemes on two- and three-dimensional domains. We report strong scalability results for up to 64 cores.