A PARALLEL BLOCK-STRUCTURED MULTIGRID METHOD FOR THE PREDICTION OF INCOMPRESSIBLE FLOWS

In this paper a parallel multigrid finite volume solver for the prediction of steady and unsteady flows in complex geometries is presented. For the handling of the complexity of the geometry and for the parallelization a unified approach connected with the concept of block-structured grids is employed. The parallel implementation is based on grid partitioning with automatic load balancing and follows the message-passing concept, ensuring a high degree of portability. A high numerical efficiency is obtained by a non-linear multigrid method with a pressure correction scheme as smoother. By a number of numerical experiments on various parallel computers the method is investigated with respect to its numerical and parallel efficiency. The results illustrate that the high performance of the underlying sequential multigrid algorithm can largely be retained in the parallel implementation and that the proposed method is well suited for solving complex flow problems on parallel computers with high efficiency.     

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

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

ADAPTIVE SOLUTIONS FOR UNSTEADY LAMINAR FLOWS ON UNSTRUCTURED GRIDS

An adaptive finite volume method for the simulation of time-dependent, viscous flow is presented. The Navier–Stokes equations are discretized by central schemes on unstructured grids and solved by an explicit Runge–Kutta method. The essential topics of the present study are a new concept for a local Runge–Kutta time-stepping scheme, called multisequence Runge–Kutta, which reduces the severe stability restriction in unsteady problems, a common grid generation and adaptation procedure and the application of dynamic grids for capturing moving flow structures. Results are presented for laminar, separated flow around an aerofoil with a flap.     

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

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

A CHARACTERISTIC-LIKE METHOD FOR THERMALLY EXPANDABLE FLOW ON UNSTRUCTURED TRIANGULAR GRIDS

Finite-difference-like discretizations are developed for the time-dependent Navier–Stokes equations and the thermal energy equation on Delaunay triangulations of the flow domain. The flow is assumed to be thermally expandable; that is, the density varies only with temperature. A characteristic- like (CL) method is used to discretize the temporal and convective terms. The dual variable method reduces the size of the discrete system by a factor of five.     

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.     

COVOLUME-DUAL VARIABLE METHOD FOR THERMALLY EXPANDABLE FLOW ON UNSTRUCTURED TRIANGULAR GRIDS

Finite difference like discretizations are developed for the time dependent Navier-Stokes equations and the thermal energy equation. The flow is assumed to be thermally expandable, that is, the density varies only with temperature. A new pointwise first order upwind scheme for convection is presented which is of nonnegative type. Also presented are new approaches to reconstructing the velocity vector field from the covolume primitive variables. The resulting difference equations reproduce linear flow fields.     

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方法,并且耦合了基于混合网格的多重网格方法.利用该方法数值模拟了不同雷诺数下的静止圆柱、振荡圆柱的绕流,得到了与实验和他人计算一致的结果.     

A HYBRID SPECTRAL/FINITE VOLUME METHOD FOR VISCOUS FLOWS

Abstract

This paper presents a hybrid spectral/finite volume method for steady-state compressible viscous flows. The method is evaluated for accuracy via test cases for various Mach numbers. The domain is divided into a viscous region and an inviscid region. The viscous region uses the full Navier-Stokes equations, while the inviscid region employs the Euler equations. A high order Chebyshev collocation spectral method is developed for the viscous region to resolve boundary layers. This method avoids the dense grids needed by finite-volume methods to resolve the viscous areas. A low order finite-volume method based on a Lax-Wendroff type scheme is employed for the inviscid region. A special interface formulation is developed for coupling the spectral with the finite-volume method. Comparisons with analytic results as well as convergence histories are presented.     

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.     

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.     

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

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

A VORTICITY–STREAMFUNCTION FORMULATION FOR STEADY INCOMPRESSIBLE TWO-DIMENSIONAL FLOWS

A vorticity–streamfunction formulation for incompressible planar viscous flows is presented. The standard kinematic field equations are discretized using centred finite difference schemes and solved in a coupled way via a Newton-like linearization scheme. The linearized system of partial differential equations is handled through the restarting linear GMRES algorithm, preconditioned by means of an incomplete LU approximate factorization. The proposed solution technique constitutes a fast and robust algorithm for treating laminar flows at high Reynolds numbers. The pressure field is obtained at a subsequent step by solving a convection– diffusion equation in terms of the stagnation pressure, which presents certain advantages compared with the widely used static pressure Poisson equation. Results are shown for a wide variety of applications including internal and external flows.     

A PRESSURE CORRECTION METHOD FOR UNSTRUCTURED MESHES WITH ARBITRARY CONTROL VOLUMES

A pressure correction procedure for general unstructured meshes is presented. It is a cell-centred, collocated finite volume method and the pressure–velocity coupling is treated using SIMPLEC. The cells can have an arbitrary number of grid points (cell vertices). In the present study the number of faces on the cells varies between three and six. The discretized equations are solved using either a symmetric Gauss–Seidel solver or a conjugate gradient solver with a preconditioner. The method is applied to three two-dimensional test cases in which the flow is incompressible and laminar. The extension to three dimensions as well as to turbulent flow using transport models is straightforward. It can also be extended to handle compressible flow.     

AN EULERIAN FINITE ELEMENT METHOD FOR TIME-DEPENDENT FREE SURFACE PROBLEMS IN HYDRODYNAMICS

A numerical method based on the finite element method is presented for simulating the two-dimensional transient motion of a viscous liquid with free surfaces. For ease of numerical treatment of the free surface expressed by a multiple-valued function, the marker particle method is employed. Numerous virtual particles are spread over all regions occupied by liquid. They move about on a fixed finite element mesh with the liquid velocity at their positions. These particles contribute nothing to the dynamics of the liquid and only serve as markers of liquid regions. The velocity field within liquid regions is calculated by solving the Navier– Stokes equations and the equation of continuity by the finite element method based on quadrilateral elements. A detailed discussion is given of the methodological problems arising in the implementation of the marker particle method on an unstructured finite element mesh and of the solutions to these problems. The proposed method is demonstrated on three sample problems: the broken dam problem, the impact of a falling liquid drop on a still liquid and the entry of a rigid block into water. Good agreement has been obtained in the comparison of the present numerical results with available experimental data.     

局部荷载作用下网架结构的拟夹层板分析法

把正交正放类网架结构简化为构造上正交异性的夹层板,采用考虑剪切变形的具有三个广义位移的平板弯曲理论来分析,给出了网架结构在局部荷载作用下的内力、位移解析计算公式,并用有限元方法进行了验证．     

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.     

THE SIMULATION OF SOME MODEL VISCOELASTIC EXTENSIONAL FLOWS

The numerical simulation of three model viscoelastic extensional flows is considered: sink flow, model draw-down and conical section draw-down. A transient finite element scheme with a pressure correction method is employed to analyse the numerical treatment of such problems for Oldroyd- Band Phan-Thien/Tanner constitutive models. Both decoupled and coupled formulations are compared for these highly convective flows and effective mechanisms are proposed for removing numerical oscillations in the temporally developing solution. In pure viscoelastic extensional flow from an initial stress-free state, the maximum stress level attained decreases with increase in material relaxation time. When this is followed by stress relaxation, as in conical section draw-down, increasing the relaxation time inhibits stress decay.