A PRESSURE-CORRECTION METHOD FOR THE SOLUTION OF INCOMPRESSIBLE VISCOUS FLOWS ON UNSTRUCTURED GRIDS

An unstructured grid, finite volume method is presented for the solution of two-dimensional viscous, incompressible flow. The method is based on the pressure-correction concept implemented on a semi-staggered grid. The computational procedure can handle cells of arbitrary shape, although solutions presented herein have been obtained only with meshes of triangular and quadrilateral cells. The discretization of the momentum equations is effected on dual cells surrounding the vertices of primary cells, while the pressure-correction equation applies to the primary-cell centroids and represents the conservation of mass across the primary cells. A special interpolation scheme s used to suppress pressure and velocity oscillations in cases where the semi-staggered arrangement does not ensure a sufficiently strong coupling between pressure and velocity to avoid such oscillations. Computational results presented for several viscous flows are shown to be in good agreement with analytical and experimental data reported in the open literature.     

THE PARALLEL BLOCK ADAPTIVE MULTIGRID METHOD FOR THE IMPLICIT SOLUTION OF THE EULER EQUATIONS

A method capable of solving very fast and robust complex non-linear systems of equations is presented. The block adaptive multigrid (BAM) method combines mesh adaptive techniques with multigrid and domain decomposition methods. The overall method is based on the FAS multigrid, but instead of using global grids, locally enriched subgrids are also employed in regions where excessive solution errors are encountered. The final mesh is a composite grid with uniform rectangular subgrids of various mesh densities. The regions where finer grid resolution is necessary are detected using an estimation of the solution error by comparing solutions between grid levels. Furthermore, an alternative domain decomposition strategy has been developed to take advantage of parallel computing machines. The proposed method has been applied to an implicit upwind Euler code (EuFlex) for the solution of complex transonic flows around aerofoils. The efficiency and robustness of the BAM method are demonstrated for two popular inviscid test cases. Up to 19-fold acceleration with respect to the single-grid solution has been achieved, but a further twofold speed-up is possible on four-processor parallel computers.     

ADAPTIVE PARALLEL MULTIGRID SOLUTION OF 2D INCOMPRESSIBLE NAVIER–STOKES EQUATIONS

In this paper an adaptive parallel multigrid method and an application example for the 2D incompressible Navier–Stokes equations are described. The strategy of the adaptivity in the sense of local grid refinement in the multigrid context is the multilevel adaptive technique (MLAT) suggested by Brandt. The parallelization of this method on scalable parallel systems is based on the portable communication library CLIC and the message-passing standards: PARMACS, PVM and MPI. The specific problem considered in this work is a two-dimensional hole pressure problem in which a Poiseuille channel flow is disturbed by a cavity on one side of the channel. Near geometric singularities a very fine grid is needed for obtaining an accurate solution of the pressure value. Two important issues of the efficiency of adaptive parallel multigrid algorithms, namely the data redistribution strategy and the refinement criterion, are discussed here. For approximate dynamic load balancing, new data in the adaptive steps are redistributed into distributed memories in different processors of the parallel system by block remapping. Among several refinement criteria tested in this work, the most suitable one for the specific problem is that based on finite-element residuals from the point of view of self-adaptivity and computational efficiency, since it is a kind of error indicator and can stop refinement algorithms in a natural way for a given tolerance. Comparisons between different global grids without and with local refinement have shown the advantages of the self-adaptive technique, as this can save computer memory and speed up the computing time several times without impairing the numerical accuracy. © 1997 By John Wiley & Sons, Ltd. Int. J. Numer. Methods Fluids 24, 875–892, 1997.     

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

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

不可压粘流N－S方程的边界积分解法

对原变量的Ｎ－Ｓ方程进行一阶时间离散，采用共轭梯度法解除压强－速度的耦合．对所得的一系列Ｌａｐｌａｃｅ方程、Ｐｏｓｓｉｏｎ方程和Ｈｅｌｍｈｏｔｚ方程均进行边界积分法求解，首次得到了粘性Ｎ－Ｓ方程的边界积分表示式．圆柱的定常、非定常尾迹计算结果表明了本文方法的有效性．     

不可压粘流N－S方程的边界积分解法

对原变量的Ｎ－Ｓ方程进行一阶时间离散，采用共轭梯度法解除压强－速度的耦合．对所得的一系列Ｌａｐｌａｃｅ方程、Ｐｏｓｓｉｏｎ方程和Ｈｅｌｍｈｏｔｚ方程均进行边界积分法求解，首次得到了粘性Ｎ－Ｓ方程的边界积分表示式．圆柱的定常、非定常尾迹计算结果表明了本文方法的有效性．     

用BGK格式计算不可压缩流场

简要介绍BGK方法的基本思想,数值方法及其对不可压缩流场数值计算的推广。然后详细介绍将BGK方法应用于平面方框流的结果,将其结果与Ghia的结果进行了详细的比较。本方法还用来模拟后台阶流动,所得结果与相应的实验数据进行比较。以上比较表明本方法的有效性,同时也确定在计算不可压缩流体时参数的正确选择。     

NUMERICAL SIMULATION OF INCOMPRESSIBLE FLOWS WITH MOVING INTERFACES

A numerical model has been developed for the 2D simulation of free surface flows or, more generally speaking, moving interface ones. The bulk fluids on both sides of the interface are taken into account in simulating the incompressible laminar flow state. In the case of heat transfer the whole system, i.e. walls as well as possible obstacles, is considered. This model is based on finite element analysis with an Eulerian approach and an unstructured fixed mesh. A special technique to localize the interface allows its temporal evolution through this mesh. Several numerical examples are presented to demonstrate the capabilities of the model. © 1997 by John Wiley & Sons, Ltd.     

MULTIGRID CONVERGENCE ACCELERATION FOR TURBULENT SUPERSONIC FLOWS

A multigrid convergence acceleration technique has been developed for solving both the Navier–Stokes and turbulence transport equations. For turbulence closure a low-Reynolds-number q–ω turbulence model is employed. To enable convergence, the stiff non-linear turbulent source terms have to be treated in a special way. Further modifications to standard multigrid methods are necessary for the resolution of shock waves in supersonic flows. An implicit LU algorithm is used for numerical time integration. Several ramped duct test cases are presented to demonstrate the improvements in performance of the numerical scheme. Cases with strong shock waves and separation are included. It is shown to be very effective to treat fluid and turbulence equations with the multigrid method. A comparison with experimental data demonstrates the accuracy of the q–ω turbulence closure for the simulation of supersonic flows. © 1997 by John Wiley & Sons, Ltd. Int. j. numer. methods fluids 24: 1019–1035, 1997.     

PEGASE: A NAVIER–STOKES SOLVER FOR DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE FLOWS

A hybrid conservative finite difference/finite element scheme is proposed for the solution of the unsteady incompressible Navier–Stokes equations. Using velocity–pressure variables on a non-staggeredgrid system, the solution is obtained with a projection method basedon the resolution of a pressure Poisson equation. The new proposed scheme is derived from the finite element spatial discretization using the Galerkin method with piecewise bilinear polynomial basis functions defined on quadrilateral elements. It is applied to the pressure gradient term and to the non-linear convection term as in the so-called group finite element method. It ensures strong coupling between spatial directions, inhibiting the development of oscillations during long-term computations, as demonstrated by the validation studies. Two- and three-dimensional unsteady separated flows with open boundaries have been simulated with the proposed method using Cartesian uniform mesh grids. Several examples of calculations on the backward-facing step configuration are reported and the results obtained are compared with those given by other methods. © 1997 by John Wiley & Sons, Ltd. Int. j. numer. methods fluids 24: 833–861, 1997.     

两种浸入式边界方法的比较

介绍两类不同的浸入式边界方法及其对它的改进. 然后采用均匀矩形交错网格和压力校正投影法,对不可压流场中的二维圆柱绕流进行了数值求解并对比了两类方法的精度.计算分析表明,连续显力法具有构造简单,适用性强的优点. 但离散隐力法在物面边界精度上要优于前者. 改进后,在二阶精度的离散格式下物面边界精度较低的显示力源法的精度有一定提高,同时发现,加密网格以提高数值精度的方法对于连续显力法并不总是有效.而同样格式下,离散隐力法具有更高精度,其中预测-校正离散隐力法可以在此基础上获得更小的计算误差和更快的收敛速度. 数值解与文献已有的数值和实验结果吻合得很好,表明边界算法及其程序是可靠和有效的.      

AN ARTIFICIAL BOUNDARY CONDITION FOR TWO-DIMENSIONAL INCOMPRESSIBLE VISCOUS FLOWS USING THE METHOD OF LINES

We design an artificial boundary condition for the steady incompressible Navier–Stokes equations in streamfunction–vorticity formulation in a flat channel with slip boundary conditions on the wall. The new boundary condition is derived from the Oseen equations and the method of lines. A numerical experiment for the non-linear Navier–Stokes equations is presented. The artificial boundary condition is compared with Dirichlet and Neumann boundary conditions for the flow past a rectangular cylinder in a flat channel. The numerical results show that our boundary condition is more accurate.     

A SEGREGATED SOLUTION ALGORITHM FOR INCOMPRESSIBLE FLOWS IN GENERAL CO-ORDINATES

To analyse an incompressible Navier–Stokes flow problem in a boundary- fitted curvilinear co-ordinate system is definitely not a trivial task. In the primitive variable formulation, choices between working variables and their storage points have to be made judiciously. The present work engages contravariant velocity components and scalar pressure which stagger each other in the mesh to prevent even–odd pressure oscillations from emerging. Now that smoothness of the pressure field is attainable, the remaining task is to ensure a discrete divergence-free velocity field for an incompressible flow simulation. Aside from the flux discretizations, the indispensable metric tensors, Jacobian and Christoffel symbols in the transformed equations should be approximated with care. The guiding idea is to get the property of geometric identity pertaining to these grid-sensitive discretizations. In addition, how to maintain the revertible one-to-one equivalence at the discrete level between primitive and contravariant velocities is another theme in the present staggered formulation. A semi-implicit segregated solution algorithm felicitous for a large-scale flow simulation was utilized to solve the entire set of basic equations iteratively. Also of note is that the present segregated solution algorithm has the virtue of requiring no user-specified relaxation parameters for speeding up the satisfaction of incompressibility in an optimal sense. Three benchmark problems, including an analytic problem, were investigated to justify the capability of the present formulation in handling problems with complex geometry. The test cases considered and the results obtained herein make a useful contribution in solving problems subsuming cells with arbitrary shapes in a boundary-fitted grid system.     

基于动态混合网格的不可压非定常流计算方法

鱼类、昆虫等运动速度较低,对它们的数值模拟需要解决不可压问题.虚拟压缩方法通过在连续性方程中加入压强对虚拟时间的偏导数,从而把压力场和速度场耦合起来,解决了不可压缩流的计算问题.基于动态混合网格技术,利用双时间步方法耦合虚拟压缩方法来解决非定常不可压缩流的计算问题.为了加快每一虚拟时间步内的收敛速度,子迭代采用了高效的块LU-SGS方法,并且耦合了基于混合网格的多重网格方法.利用该方法数值模拟了不同雷诺数下的静止圆柱、振荡圆柱的绕流,得到了与实验和他人计算一致的结果.     

亚网格尺度稳定化有限元求解不可压黏性流动

从亚网格尺度稳定化方法的基本原理出发, 提出了适合时间推进求解非定常Navier-Stokes方程获得定常解的SGS稳定化方法. 基于一定程度的近似和简化, 获得了与时间步长相关的稳定化参数, 从而排除了传统SGS稳定化方法在求解高Re数、小时间步长问题时所引发的数值不稳定性. 把SGS稳定化方法应用于求解不可压湍流, 结合标准k-\varepsilon湍流模型和壁面函数法估计湍流黏性系数, 详细讨论了壁面函数法的实施、湍流输运方程的求解和保证湍流变量非负性的限制策略, 发展了时间推进求解不可压湍流的分离式算法. 二维外掠后台阶层流和湍流计算结果表明,该方法求解不可压黏性流动是可行的, 并且具有稳定性好、计算精度高的特点.      

A MULTIBLOCK IMPLEMENTATION OF A NON-ORTHOGONAL,COLLOCATED FINITE VOLUME ALGORITHM FOR COMPLEX TURBULENT FLOWS

A multiblock algorithm for general 2D and 3D turbulent flows is introduced and applied to three cases: a compressor cascade passage, a two-element high-lift aerofoil and a round-to-square transition duct. The method is a generalization of a single-block scheme which is based on a non-orthogonal, fully collocated finite volume framework, applicable to incompressible and compressible flows and incorporating a range of turbulence transport models, including second-moment closure. The multiblock implementation is essentially block-unstructured, each block having its own local co-ordinate system unrelated to those of its neighbours. Any one block may interface with more than one neighbour along any one block face. Interblock communication is handled by connectivity matrices and effected via a two-cell overlap region along block boundaries in which ‘halo data’ reside. The algorithm and the associated data communication are explained in detail, and its effectiveness is verified, with particular reference to improved numerical resolution and parallel computing.     

轴对称高速碰撞问题的拉格朗日无结构网格有限体积法的并行计算

给出计算轴对称高速碰撞问题的拉格朗日无结构三角形网格有限体积法的并行格式，并给出以小巨型机ＡｌｉａｎｔＦＸ／４０为目标计算机的算例数值模拟结果和效率分析     

数值求解不可压缩流动的投影方法研究进展

在过去的20多年中, 投影方法通过速度和压力的解耦计算, 获得了比全耦合方法更高的计算效率, 这个显著优点使之得以广泛应用. 目前, 在计算非定常不可压缩流动的原始变量形式的数值方法中, 投影方法得到了越来越广泛的应用. 本文根据投影方法的构造思路,将众多的投影方法分成了3类, 即: Helmholtz-Hodge分解类投影方法、算子分裂类投影方法和局部连续投影方法, 并详细的介绍了3类投影方法的发展历程和求解步骤. 从投影方法的求解过程不难发现, 通过速度和压力的解耦计算, 提高了投影方法的计算效率, 但同时也给投影方法的时间精度分析带来了困难, 并长期成为大家争论的焦点. 普遍认为, 速度的时间精度比较容易达到高阶, 但是压力一般来说只有一阶精度. 但通过对3类投影的对比分析后, 我们认为, 局部连续投影方法将有助于澄清目前投影方法存在的相关争议, 并使得发展高阶精度的投影方法在理论上和技术上成为可能.      

三维非结构聚合多重网格法数值模拟研究

在三维非结构网格上应用聚合式多重网格技术来加速Euler方程的收敛过程．自行设计了一种高效率的网格聚合方法．采用四重三维非结构网格，在每一层网格上采用有限体积法进行计算．通过对M6翼型的数值求解验证了多重网格加速收敛的高效性．