On error estimates of some higher order projection and penalty-projection methods for Navier-Stokes equations

Summary This paper is a continuation of our previous work [10] on projection methods. We study first existing higher order projection schemes in the semidiscretized form for the Navier-Stokes equations. One error analysis suggests that the precision of these schemes is most likely plagued by the inconsistent Neumann boundary condition satisfied by the pressure approximations. We then propose a penalty-projection scheme for which we obtain improved error estimates.This work is partially supported by NSF grant MS-8802596.     

Fourier spectral projection method and nonlinear convergence analysis for Navier-Stokes equations

In this paper, we investigate the convergence rate of the Fourier spectral projection methods for the periodic problem of n-dimensional Navier-Stokes equations. Based on some alternative formulations of the Navier-Stokes equations and the related projection methods, the error estimates are carried out by a global nonlinear error analysis. It simplifies the analysis, relaxes the restriction on the time step size, weakens the regularity requirements on the genuine solution, and leads to some improved convergence results. A new correction technique is proposed for improving the accuracy of the numerical pressure.     

Navier-Stokes方程的局部投影稳定化方法

将Matthies,Skrzypacz和Tubiska的思想从线性的Oseen方程拓展到了非线性的Navier-Stokes方程,针对不可压缩的定常Navier-Stokes方程,提出了一种局部投影稳定化有限元方法．该方法既克服了对流占优,又绕开了inf-sup条件的限制．给出的局部投影空间既可以定义在两种不同网格上,又可以定义在相同网格上．与其他两级方法相比,定义在同一网格空间上的局部投影稳定化格式更紧凑．在同一网格上,除了给出需要bubble函数来增强的逼近空间外,还特别考虑了两种不需要用bubble函数来增强的新的空间．基于一种特殊的插值技巧,给出了稳定性分析和误差估计．最后,还列举了两个数值算例,进一步验证了理论结果的正确性．     

关于非定常不可压Navier-Stokes方程的时间高精度隐式差分方法

The incompressible Navier-Stokes equations,upon spatial discretization,become a system of differential algebraic equations,formally of index2.But due to the special forms of the discrete gradient and disrete divergence,its index can be regarded as 1.Thus,in this paper,a systematic approach following the ODE theory and methods is presented for the construction of high-order time-accurate implicit schemes for the incompressible Navier-Stokes equations,with projection methods for efficiency of numerical solution.The 3rd order 3-step BDF with componentconsistent pressure-correction projection method is a first attempt in this direction;the related iterative solution of the auxiliary velocyty,the boundary conditions and the stability of the algorithm are discussed.Results of numerical tests on the incompressible Navier-Stokes equations with an exact solution are presented,confirming the accureacy,stability and component-consistency of the proposed method.     

Efficient splitting schemes for magneto-hydrodynamic equations

We present in this paper several efficient numerical schemes for the magneto-hydrodynamic (MHD) equations. These semi-discretized (in time) schemes are based on the standard and rotational pressure-correction schemes for the Navier-Stokes equations and do not involve a projection step for the magnetic field. We show that these schemes are unconditionally energy stable, present an effective algorithm for their fully discrete versions and carry out demonstrative numerical experiments.     

A new error correction method for the stationary Navier-Stokes equations based on two local Gauss integrations

A new error correction method for the stationary Navier-Stokes equations based on two local Gauss integrations is presented. Applying the orthogonal projection technique, we introduce two local Gauss integrations as a stabilizing term in the error correction method, and derive a new error correction method. In both the coarse solution computation step and the error computation step, a locally stabilizing term based on two local Gauss integrations is introduced. The stability and convergence of the new error correction algorithm are established. Numerical examples are also presented to verify the theoretical analysis and demonstrate the efficiency of the proposed method.     

Estimating the number of asymptotic degrees of Freedom for Nonlinear Dissipative Systems

We show that the long-time behavior of the projection of the exact solutions to the Navier-Stokes equations and other dissipative evolution equations on the finite-dimensional space of interpolant polynomials determines the long-time behavior of the solution itself provided that the spatial mesh is fine enough. We also provide an explicit estimate on the size of the mesh. Moreover, we show that if the evolution equation has an inertial manifold, then the dynamics of the evolution equation is equivalent to the dynamics of the projection of the solutions on the finite-dimensional space spanned by the approximating polynomials. Our results suggest that certain numerical schemes may capture the essential dynamics of the underlying evolution equation.

     

Derivation of projection methods from formal integration of the Navier-Stokes equations

Projection methods constitute a class of numerical methods for solving the incompressible Navier-Stokes equations. These methods operate using a two-step procedure in which the zero-divergence constraint on the velocity is first relaxed while the velocity evolves, then after a certain period of time the resulting velocity field is projected onto a divergence-free subspace. Although these methods can be quite efficient, there have been certain concerns regarding their formulation. In this paper we show how a formal integration of the Navier-Stokes equations leads to a new and general procedure for the derivation of projection methods. By following this procedure, we show how each of three practical projection methods approximates a system of equations that is equivalent to the Navier-Stokes equations. We also show how the auxiliary boundary conditions required in projection methods are related to the physical boundary conditions. These results should allay the concerns regarding the legitimacy of projection methods, and may assist in their future development.     

Two-level pressure projection finite element methods for Navier-Stokes equations with nonlinear slip boundary conditions

The two-level pressure projection stabilized finite element methods for Navier-Stokes equations with nonlinear slip boundary conditions are investigated in this paper, whose variational formulation is the Navier-Stokes type variational inequality problem of the second kind. Based on the P₁-P₁ triangular element and using the pressure projection stabilized finite element method, we solve a small Navier-Stokes type variational inequality problem on the coarse mesh with mesh size H and solve a large Stokes type variational inequality problem for simple iteration or a large Oseen type variational inequality problem for Oseen iteration on the fine mesh with mesh size h. The error analysis obtained in this paper shows that if h=O(H²), the two-level stabilized methods have the same convergence orders as the usual one-level stabilized finite element methods, which is only solving a large Navier-Stokes type variational inequality problem on the fine mesh. Finally, numerical results are given to verify the theoretical analysis.     

ON DISCRETE PROJECTION AND NUMERICAL BOUNDARY CONDITIONS FOR THE NUMERICAL SOLUTION OF THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

1. IntroductionLet us consider the unsteady incompressible Navier--Stokes equations (INSE)on a two--dimensional rectangular region fl with boundary 0fl. Here w = (u, v)" is tl1e velocityvector, p is the pressure, and f a known vector function of x) y, and…     

Finite difference scheme based on proper orthogonal decomposition for the nonstationary Navier-Stokes equations

The proper orthogonal decomposition(POD)and the singular value decomposition(SVD) are used to study the finite difference scheme(FDS)for the nonstationary Navier-Stokes equations. Ensembles of data are compiled from the transient solutions computed from the discrete equation system derived from the FDS for the nonstationary Navier-Stokes equations.The optimal orthogonal bases are reconstructed by the elements of the ensemble with POD and SVD.Combining the above procedures with a Galerkin projection approach yields a new optimizing FDS model with lower dimensions and a high accuracy for the nonstationary Navier-Stokes equations.The errors between POD approximate solutions and FDS solutions are analyzed.It is shown by considering the results obtained for numerical simulations of cavity flows that the error between POD approximate solution and FDS solution is consistent with theoretical results.Moreover,it is also shown that this validates the feasibility and efficiency of POD method.     

A posteriori error estimators for the Stokes equations

Summary We present two a posteriori error estimators for the mini-element discretization of the Stokes equations. One is based on a suitable evaluation of the residual of the finite element solution. The other one is based on the solution of suitable local Stokes problems involving the residual of the finite element solution. Both estimators are globally upper and locally lower bounds for the error of the finite element discretization. Numerical examples show their efficiency both in estimating the error and in controlling an automatic, self-adaptive mesh-refinement process. The methods presented here can easily be generalized to the Navier-Stokes equations and to other discretization schemes.This work was accomplished at the Universität Heidelberg with the support of the Deutsche Forschungsgemeinschaft     

A NEW STABILIZED FINITE ELEMENT METHOD FOR SOLVING TRANSIENT NAVIER-STOKES EQUATIONS WITH HIGH REYNOLDS NUMBER

In this paper, we present a new stabilized finite element method for transient Navier-Stokes equations with high Reynolds number based on the projection of the velocity and pressure. We use Taylor-Hood elements and the equal order elements in space and second order difference in time to get the fully discrete scheme. The scheme is proven to possess the absolute stability and the optimal error estimates. Numerical experiments show that our method is effective for transient Navier-Stokes equations with high Reynolds number and the results are in good agreement with the value of subgrid-scale eddy viscosity methods, Petro-Galerkin finite element method and streamline diffusion method.     

Spectral-finite element method for solving three-dimensional unsteady Navier-Stokes equations

This paper is devoted to the combined Fourier spectral and finite element approximations of three-dimensional, semi-periodic, unsteady Navier-Stokes equations. Fourier spectral method and finite element method are employed in the periodic and non-periodic directions respectively. A class of fully discrete schemes are constructed with artificial compression. Strict error estimations are proved. The analysis shows also that the classical two-dimensional velocity-pressure elements can be readily extended to solving such three-dimensional semi-periodic problems, provided they satisfy the two-dimensional “inf-suf” condition.     

Finite element approximation of incompressible Navier-Stokes equations with slip boundary condition II

Summary We consider mixed finite element approximations of the stationary, incompressible Navier-Stokes equations with slip boundary condition simultaneously approximating the velocity, pressure, and normal stress component. The stability of the schemes is achieved by adding suitable, consistent penalty terms corresponding to the normal stress component and to the pressure. A new method of proving the stability of the discretizations allows, us to obtain optimal error estimates for the velocity, pressure, and normal stress component in natural norms without using duality arguments and without imposing uniformity conditions on the finite element partition. The schemes can easily be implemented into existing finite element codes for the Navier-Stokes equations with standard Dirichlet boundary conditions.     

非定常Navier-Stokes方程的质量集中各向异性非协调有限元逼近

该文将一个低阶Crouzeix-Raviart型非协调三角形元应用到非定常Navier-Stokes方程,给出了其质量集中有限元逼近格式.在不需要传统Ritz-Volterra投影下,通过引入两个辅助有限元空间对边界进行估计的技巧,在各向异性网格下导出了速度的L~2模和能量模及压力的L~2模的误差估计.     

非定常 N-S 方程的非线性 Galerkin 算法及其误差估计

本文给出了二维非定常N-S方程的三种数值格式,其中空间变量用谱非线性Galerkin算法进行离散,时间变量用有限差分离散,并研究了这些格式数值解的逼近精度．最后,给出了部分数值计算结果．     

The second order projection method in time for the time-dependent natural convection problem

We consider the second-order projection schemes for the time-dependent natural convection problem. By the projection method, the natural convection problem is decoupled into two linear subproblems, and each subproblem is solved more easily than the original one. The error analysis is accomplished by interpreting the second-order time discretization of a perturbed system which approximates the time-dependent natural convection problem, and the rigorous error analysis of the projection schemes is presented. Our main results of the second order projection schemes for the time-dependent natural convection problem are that the convergence for the velocity and temperature are strongly second order in time while that for the pressure is strongly first order in time.     

Analysis and convergence of the MAC scheme. II. Navier-Stokes equations

The MAC discretization scheme for the incompressible Navier-Stokes equations is interpreted as a covolume approximation to the equations. Using some results from earlier papers dealing with covolume error estimates for div-curl equation systems, and under certain conditions on the data and the solutions of the Navier-Stokes equations, we obtain first-order error estimates for both the vorticity and the pressure.

     

非线性sine-Gordon方程的各向异性线性元高精度分析新模式

在各向异性网格下,针对一类非线性sine-Gordon方程提出了线性三角形元新的高精度分析模式.基于该元的积分恒等式结果,导出了插值与Riesz投影之间的误差估计,再借助于插值后处理技术得到了在半离散和全离散格式下单独利用插值或Riesz投影所无法得到的超逼近和超收敛结果.最后,对一些常见的单元作了进一步探讨.