A convergent FEM-DG method for the compressible Navier–Stokes equations

This paper presents a new numerical method for the compressible Navier–Stokes equations governing the flow of an ideal isentropic gas. To approximate the continuity equation, the method utilizes a discontinuous Galerkin discretization on piecewise constants and a basic upwind flux. For the momentum equation, the method is a new combined discontinuous Galerkin and finite element method approximating the velocity in the Crouzeix–Raviart finite element space. While the diffusion operator is discretized in a standard fashion, the convection and time-derivative are discretized using discontinuous Galerkin on the element average velocity and a Lax–Friedrich type flux. Our main result is convergence of the method to a global weak solution as discretization parameters go to zero. The convergence analysis constitutes a numerical version of the existence analysis of Lions and Feireisl.     

A Lipschitz semigroup approach to two-dimensional Navier–Stokes equations

Motivated by the notion of integral solutions in the sense of Bénilan, a dissipative condition combined with a tangential condition is introduced. A generation problem of semigroups of Lipschitz operators is discussed under such a dissipative condition and the result is applied to the initial–boundary value problem for the Navier–Stokes equation in two-dimensional space.     

A first Integral of Navier­Stokes equations for two-dimensional flows

As wellknown, Bernoulli's equation is obtained as the first integral of Euler's equations in the absence of vorticity. Even in case of non-vanishing vorticity, a first integral from Euler's equations is obtained by using the so called Clebsch transformation [1] for inviscid flows. In contrast to this, a generalisation of this procedure towards viscous flows has not been established so far. In the present paper a first integral of Navier-Stokes equations is constructed in the case of two-dimensional flow by making use of an alternative representation of the fields in terms of complex coordinates and introducing a potential representation for the pressure. The associated boundary conditions are also considered. The first integral is a suitable tool for the development of new analytical methods and numerical codes in fluid dynamics. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Boundary partial regularity for the high dimensional Navier–Stokes equations

We consider suitable weak solutions of the incompressible Navier–Stokes equations in two cases: the 4D time-dependent case and the 6D stationary case. We prove that up to the boundary, the two-dimensional Hausdorff measure of the set of singular points is equal to zero in both cases.     

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.     

Analysis of hybrid discontinuous Galerkin methods for linearized Navier–Stokes equations

In this article, we introduce and analyze arbitrary-order, locally conservative hybrid discontinuous Galkerin methods for linearized Navier–Stokes equations. The unknowns of the global system are reduced to trace variables on the skeleton of a triangulation and the average of pressure on each cell via embedded static condensation. We prove that the lifting operator associated with trace variables is injective for any polynomial degree. This generalizes the result in (Y. Jeon and E.-J. Park, Numerische Mathematik 123 [2013], no. 1, pp. 97–119), where quadratic and cubic rectangular elements are analyzed. Moreover, optimal error estimates in the energy norm are obtained by introducing nonstandard projection operators for the hybrid DG method. Several numerical results are presented to show the performance of the algorithm and to validate the theory developed in the article.     

A parallel two-level finite element variational multiscale method for the Navier–Stokes equations

A combination method of two-grid discretization approach with a recent finite element variational multiscale algorithm for simulation of the incompressible Navier–Stokes equations is proposed and analyzed. The method consists of a global small-scale nonlinear Navier–Stokes problem on a coarse grid and local linearized residual problems in overlapped fine grid subdomains, where the numerical form of the Navier–Stokes equations on the coarse grid is stabilized by a stabilization term based on two local Gauss integrations at element level and defined by the difference between a consistent and an under-integrated matrix involving the gradient of velocity. By the technical tool of local a priori estimate for the finite element solution, error bounds of the discrete solution are estimated. Algorithmic parameter scalings are derived. Numerical tests are also given to verify the theoretical predictions and demonstrate the effectiveness of the method.     

A mixed problem for the steady Navier–Stokes equations

We consider a mixed boundary problem for the Navier–Stokes equations in a bounded Lipschitz two-dimensional domain: we assign a Dirichlet condition on the curve portion of the boundary and a slip zero condition on its straight portion. We prove that the problem has a solution provided the boundary datum and the body force belong to a Lebesgue’s space and to the Hardy space respectively.     

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.     

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.     

A defect-correction stabilized finite element method for Navier–Stokes equations with friction boundary conditions

In this paper, we consider a defect-correction stabilized finite element method for incompressible Navier–Stokes equations with friction boundary conditions whose variational formulation is the variational inequality problem of the second kind with Navier–Stokes operator. In the defect step, an artificial viscosity parameter σ is added to the Reynolds number as a stability factor, and the Oseen iterative scheme is applied in the correction step. $H^1\times L^2$ error estimations are derived for the one-step defect-correction stabilized finite element method. In the end, some numerical results are presented to verify the theoretical analysis.     

A stabilized implicit fractional-step method for the time-dependent Navier–Stokes equations using equal-order pairs

A stabilized implicit fractional-step method for numerical solutions of the time-dependent Navier–Stokes equations is presented in this paper. The time advancement is decomposed into a sequence of two steps: the first step has the structure of the linear elliptic problem; the second step can be seen as the generalized Stokes problem. The two problems satisfy the full homogeneous Dirichlet boundary conditions on the velocity. On the other hand, a locally stabilized term is added in the second step of the schemes. It allows one to enhance the numerical stability and efficiency by using the equal-order pairs. Convergence analysis and error estimates for the velocity and pressure of the schemes are established via the energy method. Some numerical experiments are also used to demonstrate the efficiency of this new method.     

Asymptotic behavior of stochastic two-dimensional Navier–Stokes equations with delays

The paper proves the L ²-exponential stability of weak solutions of two-dimensional stochastic Navier?CStokes equations in the presence of delays. The results extend some of the existing results.     

A logarithmically improved regularity criterion for the Navier–Stokes equations

In this note we prove a logarithmically improved regularity criterion in terms of the Besov space norm for the Navier–Stokes equations. The result shows that if a mild solution u satisfies ${\int_{0}^{T}\frac{\|u (t,\cdot)\|_{{\dot{B}}_{\infty,\infty}^{-r}}^{\frac{2}{1-r}}}{1+\ln(e+\| u(t,\cdot)\|_{H^{s}})}\text{d}t < \infty}$ for some 0?≤ r?<?1 and ${s\geq\frac{n}{2}-1}$ , then u is regular at t?=?T.     

A class of hemivariational inequalities for nonstationary Navier–Stokes equations

This paper is devoted to the study of a class of hemivariational inequalities for the time-dependent Navier–Stokes equations, including both boundary hemivariational inequalities and domain hemivariational inequalities. The hemivariational inequalities are analyzed in the framework of an abstract hemivariational inequality. Solution existence for the abstract hemivariational inequality is explored through a limiting procedure for a temporally semi-discrete scheme based on the backward Euler difference of the time derivative, known as the Rothe method. It is shown that solutions of the Rothe scheme exist, they contain a weakly convergent subsequence as the time step-size approaches zero, and any weak limit of the solution sequence is a solution of the abstract hemivariational inequality. It is further shown that under certain conditions, a solution of the abstract hemivariational inequality is unique and the solution of the abstract hemivariational inequality depends continuously on the problem data. The results on the abstract hemivariational inequality are applied to hemivariational inequalities associated with the time-dependent Navier–Stokes equations.     

A factorization method for numerical solution of the Navier–Stokes equations for a viscous incompressible liquid

Some implicit difference scheme of approximate factorization is proposed for numerical solution of the Navier–Stokes equations for an incompressible liquid in curvilinear coordinates. Testing of the algorithm is carried out on the solution of the problems concerning the Couette and Poiseuille flows; and the results are presented of numerical simulation of a flow between the rotating cylinders with covers.     

A second order accuracy for a full discretized time-dependent Navier–Stokes equations by a two-grid scheme

We study a second-order two-grid scheme fully discrete in time and space for solving the Navier–Stokes equations. The two-grid strategy consists in discretizing, in the first step, the fully non-linear problem, in space on a coarse grid with mesh-size H and time step Δt and, in the second step, in discretizing the linearized problem around the velocity u _H computed in the first step, in space on a fine grid with mesh-size h and the same time step. The two-grid method has been applied for an analysis of a first order fully-discrete in time and space algorithm and we extend the method to the second order algorithm. This strategy is motivated by the fact that under suitable assumptions, the contribution of u _H to the error in the non-linear term, is measured in the L ² norm in space and time, and thus has a higher-order than if it were measured in the H ¹ norm in space. We present the following results: if h ² = H ³ = (Δt)², then the global error of the two-grid algorithm is of the order of h ², the same as would have been obtained if the non-linear problem had been solved directly on the fine grid.