Stability of Navier-Stokes Equation(Ⅰ)

ＳＴＡＢＩＬＩＴＹＯＦＮＡＶＩＥＲ－ＳＴＯＫＥＳＥＱＵＡＴＩＯＮ（Ⅰ）ＳｈｉＷｅｉ－ｈｕｉ（施惟慧）（ＳｈａｎｇｈａｉＵｎｉｖｅｒｓｉｔｙ，Ｓｈａｎｇｈａｉ）（ＲｅｃｅｉｖｅｄＤｅｃ．１０，１９９３；ＣｏｍｍｕｎｉｃａｔｅｄｂｙＣｈｉｅｎＷｅｉ－ｚａ...     

二维定常不可压缩粘性流动N-S方程的数值流形方法

将流形方法应用于定常不可压缩粘性流动N-S方程的直接数值求解,建立基于Galerkin加权余量法的N-S方程数值流形格式,有限覆盖系统采用混合覆盖形式,即速度分量取1阶和压力取0阶多项式覆盖函数,非线性流形方程组采用直接线性化交替迭代方法和Nowton-Raphson迭代方法进行求解.将混合覆盖的四节点矩形流形单元用于阶梯流和方腔驱动流动的数值算例,以较少单元获得的数值解与经典数值解十分吻合.数值实验证明,流形方法是求解定常不可压缩粘性流动N-S方程有效的高精度数值方法.     

A new method for obtaining the exact solutions of the Navier-Stokes (NS) equations and its applications

In this paper we propose a new method for obtaining the exact solutions of the Mavier-Stokes (NS) equations for incompressible viscous fluid in the light of the theory of simplified Navier-Stokes (SNS) equations developed by the first author^[1,2], Using the present method we can find some new exact solutions as well as the well-known exact solutions of the NS equations. In illustration of its applications, we give a variety of exact solutions of incompressible viscous fluid flows for which NS equations of fluid motion are written in Cartesian coordinates, or in cylindrical polar coordinates, or in spherical coordinates. The project supported by National Natural Science Foundation of China.     

A boundary integral equation method for Navier-Stokes equations—application to flow in annulus of eccentric cylinders

A novel Navier-Stokes solver based on the boundary integral equation method is presented. The solver can be used to obtain flow solutions in arbitrary 2D geometries with modest computational effort. The vorticity transport equation is modelled as a modified Helmholtz equation with the wave number dependent on the flow Reynolds number. The non-linear inertial terms partly manifest themselves as volume vorticity sources which are computed iteratively by tracking flow trajectories. The integral equation representations of the Helmholtz equation for vorticity and Poisson equation for streamfunction are solved directly for the unknown vorticity boundary conditions. Rapid computation of the flow and vorticity field in the volume at each iteration level is achieved by precomputing the influence coefficient matrices. The pressure field can be extracted from the converged streamfunction and vorticity fields. The solver is validated by considering flow in a converging channel (Hamel flow). The solver is then applied to flow in the annulus of eccentric cylinders. Results are presented for various Reynolds numbers and compared with the literature.     

Impact of singularity of Navier-Stokes equation upon atmospheric motion equations

Some conclusions about the smooth function classes stability for the basic system of equations of atmospheric motion and instability for Navier-Stokes equation are summarized. On the basis of this, by taking the basic system of equations of atmospheric motion via Boussinesq approximation as example to explain in detail that the instability about some simplified models of the basic system of equations for atmospheric motion is caused by the instability of Navier-Stokes equation, thereby, a principle to guarantee the stability of simplified equation is drawn in simplifying the basic system of equations.     

On pressure boundary conditions for the incompressible Navier-Stokes equations

The pressure is a somewhat mysterious quantity in incompressible flows. It is not a thermodynamic variable as there is no ‘equation of state’ for an incompressible fluid. It is in one sense a mathematical artefact—a Lagrange multiplier that constrains the velocity field to remain divergence-free; i.e., incompressible—yet its gradient is a relevant physical quantity: a force per unit volume. It propagates at infinite speed in order to keep the flow always and everywhere incompressible; i.e., it is always in equilibrium with a time-varying divergence-free velocity field. It is also often difficult and/or expensive to compute. While the pressure is perfectly well-defined (at least up to an arbitrary additive constant) by the governing equations describing the conservation of mass and momentum, it is (ironically) less so when more directly expressed in terms of a Poisson equation that is both derivable from the original conservation equations and used (or misused) to replace the mass conservation equation. This is because in this latter form it is also necessary to address directly the subject of pressure boundary conditions, whose proper specification is crucial (in many ways) and forms the basis of this work. Herein we show that the same principles of mass and momentum conservation, combined with a continuity argument, lead to the correct boundary conditions for the pressure Poisson equation: viz., a Neumann condition that is derived simply by applying the normal component of the momentum equation at the boundary. It usually follows, but is not so crucial, that the tangential momentum equation is also satisfied at the boundary.     

UNSTEADY／STEADY NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON ARTIFICIAL COMPRESSIBILITY

A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.     

COMPARISON OF STABILITY BETWEEN NAVIER-STOKES AND EULER EQUATIONS

The stability about Navier-Stokes equation and Euler equation was brought into comparison. And by taking their typical initial value problem for example, the reason of leading to the difference in stability between Navier-Stokes equation and Euler equation was also analyzed.     

Several exact solutions of the simplified Navier-Stokes equations (SNSE)

The Simplified Navier-Stokes equations (SNSE) and their exact solutions for the flow near a rotating disk and the flow in the vicinity of a stagnation point for both two- and three-dimensional flows are presented in this paper. The analysis shows that in the aforementioned cases the exact solutions of the inner-outer-layer-matched SNSE^[4] are completely consistent with those of the complete Navier-Stokes equations (CNSE) and that the exact velocity solutions of D-SNSE^[1,3] agree with those of CNSE, however, the exact pressure solutions of D-SNSE do not agree with those of CNSE. The maximum relative pressure errors between the exact solutions of D-SNSE and CNSE can be as high as a hundred per cent.     

Stability of Navier-Stokes Equation(Ⅱ)

ＳＴＡＢＩＬＩＴＹＯＦＮＡＶＩＥＲ－ＳＴＯＫＥＳＥＱＵＡＴＩＯＮ（ＩＩ）ＳｈｉＷｅｉ－ｈｕｉ（施惟慧）;ＦａｎｇＸｉａｏ－ｚｕｏ（方晓佐）（ＳｈａｎｇｈａｉＵｎｉｖｅｒｓｉｔｙ）,Ｓｈａｎｇｈａｉ（ＲｅｃｅｉｖｅｄＤｅｃ．１１,１９９３;Ｃｏｍｍｕｎ...     

High-order accurate and high-resolution upwind finite volume scheme for solving euler/reynolds-averaged Navier-Stokes equations

A high-order accurate explicit scheme is proposed for solving Euler/Reynolds-averaged Navier-Stokes equations for steady and unsteady flows, respectively. Baldwin-Lomax turbulence model is utilized to obtain the turbulent viscosity. For the explicit scheme, the Runge-Kutta time-stepping methods of third orders are used in time integration, and space discretization for the right-hand side (RHS) terms of semi-discrete equations is performed by third-order ENN scheme for inviscid terms and fourth-order compact difference for viscous terms. Numerical experiments suggest that the present scheme not only has a fairly rapid convergence rate, but also can generate a highly resolved approximation to numerical solution, even to unsteady problem. The project supported by the National Natural Science Foundation of China under Contract No. 59576007 and 19572038     

A deterministic vortex method for solving the Navier-Stokes equations

In this paper, a new 2-D vortex method is developed, which treats the vorticity diffusion in a deterministical way. The Laplacian operator, which describes vorticity diffusion, is approximated by a contour integral. The numerical results of two model problems show that this method has a good accuracy. A primary error estimation is given, and the self-adaptive vortex blob and the boundary conditions are discussed. The project supported by the National Natural Science Foundation of China     

90°弯管内流动的理论模型及流动特性的数值研究

从三维不可压缩雷诺时均Navier-Stokes方程出发,对90°弯曲管道内湍流流动进行数值模拟。网格划分采用六面体网格,湍流模型为RNGk-ε模型,在近壁区采用两层壁面模型进行修正,流场的计算结果与实验数据吻合较好。在此基础上,本文数值研究了来流方向对流场结构和流动特性的影响。得出在弯管流场中发生了分离现象,且随着来流侧滑角的增大,分离区范围增大。此外,随着来流从同一侧滑角变换至同一攻角时,横截面的二次流图像中也从具有两个对称主涡变成只具有一个主涡的现象。     

Numerical study of Navier-Stokes equations

The problem of a viscous incompressible fluid flow around a body of revolution at incidence, which is described by Navier-Stokes equations, is considered. For low Reynolds numbers, the solutions of these equations are smooth functions. A numerical algorithm without saturation is constructed, which responds to solution smoothness. The calculations are performed on grids consisting of 900 (10 × 10 × 9) and 700 (10 × 10 × 7) nodes. On the grid consisting of 900 nodes, a system of 3600 nonlinear equations is solved by a standard code. The pressures on the shaded side of the body of revolution are compared. It is found that a numerical study (on this grid) is feasible for problems with Re ≈ 1. For high Reynolds numbers, the number of grid nodes has to be increased. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 5, pp. 43–52, September–October, 2007.     

Dimensions of attractors for discretizations for Navier-Stokes equations

In this paper, we discretize the 2-D incompressible Navier-Stokes equations with the periodic boundary condition by the finite difference method. We prove that with a shift for discretization, the global solutions exist. After proving some discrete Sobolev inequalities in the sense of finite differences, we prove the existence of the global attractors of the discretized system, and we estimate the upper bounds for the Hausdorff and the fractal dimensions of the attractors. These bounds are indepent of the mesh sizes and are considerably close to those of the continuous version.     

TWO-GRID ERROR ESTIMATES FOR THE STREAM FUNCTION FORM OF NAVIER-STOKES EQUATIONS

ＩｎｔｒｏｄｕｃｔｉｏｎＷｅｃｏｎｓｉｄｅｒｔｗｏ＿ｇｒｉｄｍｅｔｈｏｄｆｏｒｔｈｅｓｔｒｅａｍｆｕｎｃｔｉｏｎｆｏｒｍｏｆｔｈｅｓｔａｔｉｏｎａｒｙＮａｖｉｅｒ＿Ｓｔｏｋｅｓｅｑｕａｔｉｏｎｓ.Ｔｈｅａｄｖａｎｔａｇｅｓｏｆｔｈｅｓｔｒｅａｍｆｕｎｃｔｉｏｎｆｏｒｍａｒｅｔｈａｔｔｈｅｉｎｃｏｍｐｒｅｓｓｉｂｉｌｉｔｙｃｏｎｄｉｔｉｏｎｉｓｓａｔｉｓｆｉｅｄａｕｔｏｍａｔｉｃａｌｌｙａｎｄｔｈｅｐｒｅｓｓｕｒｅｉｓｎｏｔｐｒｅｓｅｎｔｉｎｔｈｅｗｅａｋｆｏｒｍ .Ｔｈｅｍｅｔｈｏｄｉｓｂａｓｅｄ…     

Optimal Control in Navier-Stokes Equations

This paper presents a formulation for optimal control of a forced convection flow. The state equation that governs the forced convection flow can be expressed as the incompressible Navier-Stokes equations and energy equation. The optimal control can be formulated as finding a control force to minimize a performance function that is defined to evaluate a control object. The stabilized finite element method is used for the spatial discretization, while the Crank-Nicolson scheme is used for the temporal discretization. The Sakawa-Shindo method, which is an iterative procedure, is applied for minimizing the performance function.     

Residual a posteriori error estimate two-grid methods for the steady Navier-Stokes equation with stream function form

Residual based on a posteriori error estimates for conforming element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level method were derived. The posteriori error estimates contained additional terms in comparison to the error estimates for the solution obtained by the standard finite element method. The importance of these additional terms in the error estimates was investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than of convergence of discrete solution.     

TheC k instability of navier-Stokes equation appending polynomials of unknown functions

Applying the theory of stratification, the solution space structure about a class of deformed Navier-Stokes equation isdetermined. It is proved that such kind of equation has no C^k (k≥2) stable solution by the fact that the strate transversale is a null set. Foundation item: the National Natural Science Foundation of China (19971054), Shanghai Natural Science Foundation (99ZA14034) Biography: HE You-hua (1960-)     

On the stability of general Navier-Stokes type equation

Ｉｎｔｈｉｓｐａｐｅｒ,ｗｈａｔｗｅｄｉｓｃｕｓｂｅｌｏｎｇｓｔｏａｓｔａｂｉｌｉｔｙｐｒｏｂｌｅｍａｂｏｕｔｔｈｅｓｅｃｏｎｄｃｌａｓｇｅｎｅｒａｌＮａｖｉｅｒ＿Ｓｔｏｋｅｓｅｑｕａｔｉｏｎｓ．Ｉｔｔｏｕｃｈｅｓｓｏｍｅｉｍｐｏｒｔａｎｔｐｒｏｂｌ...