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

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

Lagrangian finite element analysis applied to viscous free surface fluid flow

A new Lagrangian finite element formulation is presented for time-dependent incompressible free surface fluid flow problems described by the Navier-Stokes equations. The partial differential equations describing the continuum motion of the fluid are discretized using a Galerkin procedure in conjunction with the finite element approximation. Triangular finite elements are used to represent the dependent variables of the problem. An effective time integration procedure is introduced and provides a viable computational method for solving problems with equality of representation of the pressure and velocity fields. Its success has been attributed to the strict enforcement of the continuity constraint at every stage of the iterative process. The capabilities of the analysis procedure and the computer programs are demonstrated through the solution of several problems in viscous free surface fluid flow. Comparisons of results are presented with previous theoretical, numerical and experimental results.     

谱元方法求解环形空间内自然对流换热

提出了将谱元方法应用到极坐标系下，利用极坐标系下的谱元方法求解环形空间内自然对流问题。具体求解了原始变量速度和压力的不可压缩Navier-Stokes方程和能量方程，通过在时间方向采用时间分裂方法和空间采用谱元方法对方程进行离散求解，取得了与基准解较一致的计算结果。     

边界元法分析粘性液体小幅晃动

本文采用了边界元法对容器中粘性、不可压缩液体小幅晃动进行数值分析。在频域内考虑二维线性化Navier-Stokes方程,以问题的物理变量作为数值分析的未知函数,并推导了该问题分析的边界积分方程。自由面上的动力学条件为法向正应力和切向剪应力为零,这两个条件本身是线性的,避免了采用无粘势理论边界条件的非线性,固壁面上采用流体质点与固壁质点速度相等的条件,在数值计算过程中,结合有限差分法对边界条件进行了处理,由此建立了问题的一个边界元数值求解过程。     

Using a segregated finite element scheme to solve the incompressible Navier-Stokes equations

In this paper, a segregated finite element scheme for the solution of the incompressible Navier-Stokes equations is proposed which is simpler in form than previously reported formulations. A pressure correction equation is derived from the momentum and continuity equations, and equal-order interpolation is used for both the velocity components and pressure. Algorithms such as this have been known to lead to checkerboard pressure oscillations; however, the pressure correction equation of this scheme should not produce these oscillations. The method is applied to several laminar flow situations, and details of the methods used to achieve converged solutions are given.     

A hybrid vertex-centered finite volume/element method for viscous incompressible flows on non-staggered unstructured meshes

This paper proposes a hybrid vertex-centered finite volume/finite element method for solution of the two dimensional (2D) incompressible Navier-Stokes equations on unstructured grids.An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling.The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by joining the centroid of cells sharing the common vertex.For the temporal integration of the momentum equations,an implicit second-order scheme is utilized to enhance the computational stability and eliminate the time step limit due to the diffusion term.The momentum equations are discretized by the vertex-centered finite volume method (FVM) and the pressure Poisson equation is solved by the Galerkin finite element method (FEM).The momentum interpolation is used to damp out the spurious pressure wiggles.The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both velocity and pressure.The classic test cases,the lid-driven cavity flow,the skew cavity flow and the backward-facing step flow,show that numerical results are in good agreement with the published benchmark solutions.     

A streamline diffusion nonconforming finite element method for the time-dependent linearized Navier-Stokes equations

A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.     

Finite element solution of the incompressible Navier-Stokes equations by a Helmholtz velocity decomposition

Finite element solution methods for the incompressible Navier-Stokes equations in primitive variables form are presented. To provide the necessary coupling and enhance stability, a dissipation in the form of a pressure Laplacian is introduced into the continuity equation. The recasting of the problem in terms of pressure and an auxiliary velocity demonstrates how the error introduced by the pressure dissipation can be totally eliminated while retaining its stabilizing properties. The method can also be formally interpreted as a Helmholtz decomposition of the velocity vector. The governing equations are discretized by a Galerkin weighted residual method and, because of the modification to the continuity equation, equal interpolations for all the unknowns are permitted. Newton linearization is used and at each iteration the linear algebraic system is solved by a direct solver. Convergence of the algorithm is shown to be very rapid. Results are presented for two-dimensional flows in various geometries.     

Low order nonconforming mixed finite element method for nonstationary incompressible Navier-Stokes equations

This paper studies a low order mixed finite element method (FEM) for nonstationary incompressible Navier-Stokes equations. The velocity and pressure are approximated by the nonconforming constrained Q ₁ ^rot element and the piecewise constant, respectively. The superconvergent error estimates of the velocity in the broken H ¹-norm and the pressure in the L ²-norm are obtained respectively when the exact solutions are reasonably smooth. A numerical experiment is carried out to confirm the theoretical results.     

A Galerkin finite element algorithm based on third-order Runge-Kutta temporal discretization along the uniform streamline for unsteady incompressible flows

In this paper, for two-dimensional unsteady incompressible flow, the Navier-Stokes equations without convection term are derived by the coordinate transformation along the streamline characteristic. The third-order Runge-Kutta method along the streamline is introduced to discrete the alternative Navier-Stokes equations in time, and spacial discretization is carried out by the Galerkin method, and then, the third-order accuracy finite element method is obtained. Meanwhile, the streamline velocity is uniformly approximated by initial velocity in each time step in order to reduce update frequency of total element matrix and improve calculation efficiency. Finally, some classic unsteady flow examples are calculated and analyzed by different calculation methods, which further demonstrate that the present method has more advantages in stability, permissible time step, dissipation, computational cost, and accuracy. The code can be downloaded at https://doi.org/10.13140/RG.2.2.27706.44484 .     

A high-order characteristics/finite element method for the incompressible Navier-Stokes equations

In this paper we consider a discretization of the incompressible Navier-Stokes equations involving a second-order time scheme based on the characteristics method and a spatial discretization of finite element type. Theoretical and numerical analyses are detailed and we obtain stability results abnd optimal eror estimates on the velocity and pressure under a time step restriction less stringent than the standard Courant-Freidrichs-Levy condition. Finally, some numerical results obtained wiht the code N3S are shown which justify the interest of this scheme and its advantages with respect to an analogous first-order time scheme. © 1997 John Wiley & Sons, Ltd.     

不可压Navier-Stokes方程基于误差估算的网格自适应解法

首先导出了广义Stokes方程Petrov—Galerkin有限元数值解的当地事后误差估算公式;以非连续二阶鼓包（bump）函数空间为速度、压强误差的近似空间,该估算基于求解当地单元上的广义Stokes问题。然后,证明了误差估算值与精确误差之间的等价性。最后,将误差估算方法应用于Navier—Stokes环境,以进行不可压粘流计算中的网格自适应处理。数值实验中成功地捕获了多强度物理现象,验证了本文所发展的方法。     

不可压黏性流问题的无网格法分析

不可压缩黏性流问题一般采用Navier-Stokes方程来描述，基于加权残值法，推导了问题的无网格伽辽金法（EFGM）离散Navier-Stokes方程，在时间域上采用分步方法计算，速度和压力由相互独立的方程以解耦的形式求解，并采用同阶移动最小二乘近似，在每一时间步中，对压力解和速度解采用了Newton-Raphson迭代法进行修正，最后将所得到的方法应用到剪切驱动空腔流问题中，验证了方法的有效性，且解的精度高、稳定性好。     

Numerical computation of an impinging jet with uniform wall suction

The paper presents numerical predictions of a turbulent axisymmetric jet impinging onto a porous plate, based on a finite volume method of solving the Navier-Stokes equations for an incompressible air jet with the K–ε turbulence model. The velocity and pressure terms of the momentum equations are solved by the SIMPLE (semi-implicit method for pressure-linked equation) method. In this study, non-uniform staggered grids are used. The parameters of interest include the nozzle-to-wall distance and the suction velocity. The results of the present calculations are compared with available data reported in the literature. It is found that suction effects reduce the boundary layer thickness and increase the velocity gradient near the wall.     

Viscous incompressible flow over the surface of a rotary disk

Viscous incompressible film flow over the surface of an impermeable rotary disk is studied. An exact self-similar solution of the complete Navier-Stokes system of equations is obtained and the velocity and pressure fields together with the radial profiles of the fluid film are determined. A physical interpretation of the results obtained is given.Volgograd. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 39–43, November–December, 1995.     

A new method for computing solenoidal vector fields on arbitrary regions

We develop a new method for the efficient calculation of solenoidal vector fields on general regions. The method takes advantage of fast direct methods and uses boundary integral equations to satisfy boundary conditions. For the latter we give an effective scheme for computing far-field boundary influences (based on discrete charges). Examples and numerical results are given. The method is applicable to incompressible Navier-Stokes calculations.     

隐格式近似不可压四面体4节点大变形单元

由于在处理体积自锁方面的优势,近似不可压问题的大变形求解多采用六面体单元/网格,但对于复杂工程问题,由于网格剖分上的限制,往往更需要一种可以很好解决体积自锁的四面体单元。Bonet和Burton的平均节点压力4节点四面体单元是为数不多能够较好处理体积自锁问题的四面体单元之一,但是该单元目前主要用于显式计算。利用单元平均压力对位移增量的精确方向导数,得到了严格的一致切线阵,保证了Newton-Raphson迭代的二阶收敛,从而使得该单元可以用于隐式计算。该单元的压力平均计算会耦合相邻单元的节点自由度,从而增加切线刚度阵的非零带宽,但不增加自由度总数。分别采用线性六面体选择缩减积分单元、标准线性四面体单元和本文的单元计算了3个近似不可压的典型算例。算例表明,本文推导的单元可以有效克服体积自锁,达到与常用六面体单元相近的效果,使得四面体网格可以方便地用于不可压问题的大变形隐式求解。     

Friction and heat transfer in an incompressible fluid near a plane stagnation point

In the neighborhood of a plane stagnation point, the flow and heat transfer of an incompressible fluid are studied. In the inner flow region, the velocity and pressure fields are described by the complete Navier-Stokes equations, and the temperature field is described by the complete energy equation. In the outer flow region, a two-term asymptotic solution of the corresponding equations is obtained. The problem is reduced to the numerical solution of ordinary differential equations. Numerical results are discussed.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 52–65, July–August, 1996.     

靠近排列的两个交错方柱三维绕流的数值模拟

通过用时间分裂算法求解Navier－Stokes方程,对中等Reynolds数下靠近排列的两个交错方柱三维绕流进行了数值模拟,其中,中间速度场用四阶Adams格式计算,压力场通过结合近似因子分解方法AF1与稳定的双共轭梯度方法Bi-CGSTAB进行迭代求解．数值模拟发现当两个方柱靠得较近时,有互相吸引趋势,而且上游方柱的Strouhal数较大．方柱的交错排列方式对绕流影响明显．计算结果与实验定性吻合,而且比用MAC－AF1方法计算的结果好．     

Multifluid flows with weak and strong discontinuous interfaces using an elemental enriched space

In a previous paper, the authors presented an elemental enriched space to be used in a finite‐element framework (EFEM) capable of reproducing kinks and jumps in an unknown function using a fixed mesh in which the jumps and kinks do not coincide with the interelement boundaries. In this previous publication, only scalar transport problems were solved (thermal problems). In the present work, these ideas are generalized to vectorial unknowns, in particular, the incompressible Navier‐Stokes equations for multifluid flows presenting internal moving interfaces. The advantage of the EFEM compared with global enrichment is the significant reduction in computing time when the internal interface is moving. In the EFEM, the matrix to be solved at each time step has not only the same amount of degrees of freedom (DOFs) but also the same connectivity between the DOFs. This frozen matrix graph enormously improves the efficiency of the solver. Another characteristic of the elemental enriched space presented here is that it allows a linear variation of the jump, thus improving the convergence rate, compared with other enriched spaces that have a constant variation of the jump. Furthermore, the implementation in any existing finite‐element code is extremely easy with the version presented here because the new shape functions are based on the usual finite‐element method shape functions for triangles or tetrahedrals, and once the internal DOFs are statically condensed, the resulting elements have exactly the same number of unknowns as the nonenriched finite elements.