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.     

SINGLE-LAYER FINITE ELEMENT APPROXIMATION OF 3D NAVIER–STOKES EQUATIONS WITH FREE BOUNDARIES FOR SHALLOW WATER FLOWS

In order to simulate flows in the shallow water limit, the full incompressible Navier–Stokes equations with free boundaries are solved using a single layer of finite elements. This implies a polynomial approximation of the velocity profile in the vertical direction, which in turn distorts the wave speed. This fact is verified by numerical results: the wave speed depends on the vertical discretization. When at least two layers of finite elements are used, the boundary layer at the bottom can be simulated and the correct solution for the shallow water limit is recovered. Then this algorithm is applied to the prediction of Tsunami event.     

A FINITE ELEMENT MODEL FOR THE SIMULATION OF DENSE GAS DISPERSION IN THE ATMOSPHERE

A state-of-the-art model is developed for the simulation of the dispersion of hazardous toxic or flammable gases heavier than air in the atmosphere. The model depends on solving the Reynolds-averaged Navier–Stokes equations in addition to the energy equation and a species concentration equation for the contaminant gas. Turbulence closure is achieved by using a buoyancy-extended version of the standard k– ϵ two-equation model. The buoyancy extension is introduced to account for the anisotropic turbulent viscosity resulting from the strong stratification introduced by the dense gas clouds. The spatial discretization is achieved via the Galerkin finite element method, while the solution is advanced in time using the forward Euler method. A special element layer is introduced in the near-ground region to bridge the gap between the solid wall and the main solution domain where the turbulence model can be applied. This special element layer eliminates the need to apply the wall function in the standard way where any oscillations in the pressure field could contaminate the velocity solution. The model was tested against the Burro-8 field trial and could predict the experiment satisfactorily to within the experimental uncertainties of the reported results. © by 1997 John Wiley & Sons, Ltd.     

AN IMPROVED NAVIER–STOKES PROCEDURE FOR ANALYSIS OF TWO-DIMENSIONAL AEROFOILS IN LAMINAR AND TURBULENT FLOWS

An accurate and robust Navier–Stokes procedure to predict the complex flow about an aerofoil has been developed. Much improvement over existing methods is achieved in various aspects of the solution procedure. The computational grid generated by conformal mapping, which is not only orthogonal but aligned with the inviscid streamlines, keeps the equations simple and minimizes the error due to false diffusion. Formal second-order accuracy is ensured by employing the QUICK scheme for the convective derivatives in the full Navier–Stokes and turbulence transport equations. To treat the separated region properly and to better resolve the flow field in the wake, the two-layer k–ε turbulence model is incorporated. The onset of transition is triggered in a unique fashion to warrant the smooth transition to turbulent flow. Sample calculations for various aerofoil sections show that the prediction is improved substantially over those by existing methods. The details of the flow extending to the wake, such as the surface pressure distribution, C_Lmax, the velocity fields and the Reynolds stress profiles, are found to be in excellent agreement with the data. © 1997 John Wiley & Sons, Ltd.     

自适应有限元方法及其工程应用

自适应有限元方法是一种能通过自适应分析自动调整算法以改进求解过程的数值方法．它以常规有限元法为基础,以误差估计和自适应网格改进技术为核心,是一种效率高、可靠性高的计算方法．文中简要介绍并综述了自适应有限元方法的重要发展及应用情况．并对其发展前景作了概要的预测      

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.     

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.     

多尺度材料模拟的自适应有限元方法研究

研究了材料模拟中一类新型耦合多尺度的自适应有限元方法. 采用微观分子动力学耦合宏观有限元的桥尺度方法来模拟材料破坏的前期行为,其中宏观有限元计算推广到了一般非结构三角形网格. 材料破坏形成后,停止微观尺度的计算,它的进一步发展和演化通过一个宏观模型来描述,采用自适应有限元方法来求解这一宏观模型. 其中,后验误差估计的基础是变分多尺度理论,即自适应网格加密是基于粗尺度上残差分布和细尺度上单元Green's函数. 计算中采用了破坏准则来模拟材料的断裂. 数值实验表明了方法的有效性.      

一种有限元模型修正中的参数选择方法

有限元模型修正面临的首要问题是待修正参数的选择.目前主要是基于灵敏度分析的结果进行参数选择.本文从结构固有频率分析的能量法出发,论证了某阶固有频率主要受到同阶振型中发生弹性变形部位刚度的影响,据此提出一种待修正参数的选择方法,以振型中弹性变形为依据,选择弹性变形较大部位的参数作为待修正参数.用一个梁模型仿真验证了该方法的正确性,并将其应用于系杆拱桥的有限元模型修正中,修正后有限元计算结果与试验结果的最大误差缩小至3.3%.     

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.     

AN EXPLICIT FINITE VOLUME SPATIAL MARCHING METHOD FOR REDUCED NAVIER–STOKES EQUATIONS

This paper develops a spatial marching method for high-speed flows based on a finite volume approach. The method employs the reduced Navier– Stokes equations and a pressure splitting in the streamwise direction based on the Vigneron strategy. For marching from an upstream station to one downstream the modified five-level Runge–Kutta integration scheme due to Jameson and Schmidt is used. In addition, for shock handling and for good convergence properties the method employs a matrix form of the artificial dissipation terms, which has been shown to improve the accuracy of predictions. To achieve a fast rate of convergence, a local time-stepping concept is used. The method retains the time derivative in the governing equations and the solution at every spatial station is obtained in an iterative manner. The developed method is validated against two test cases: (a) supersonic flow past a flat plate; and (b) hypersonic flow past a compression corner involving a strong viscous–inviscid interaction. The computed wall pressure and wall heat transfer coefficients exhibit good general agreement with previous computations by other investigators and with experiments.     

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 ADAPTIVE FINITE ELEMENT SCHEME FOR HYDRAULIC FRACTURING WITH PROPPANT TRANSPORT

A mathematical model and adaptive finite element scheme are developed for describing the distribution of proppant in a propagating hydraulic fracture. The governing equation for proppant concentration is derived by applying the conservation law of mass to the proppant and to the proppant-laden fluid. Shah's empirical equation, which relates the proppant concentration and the indices of the non-Newtonian fluid, is used to describe the proppant-laden fluid. The proppant distribution inside a hydraulic fracture can then be obtained by solving the proppant concentration equation together with the governing equations of fluid and elasticity for a hydraulic fracturing. A novel moving grid scheme is developed that combines grid point insertion with redistribution. Four examples corresponding to different in situ stress distributions are computed to demonstrate the scheme. © 1997 by John Wiley & Sons, Ltd.     

AN ADAPTIVE UNSTRUCTURED TRI-TREE ITERATIVE SOLVER FOR MIXED FINITE ELEMENT FORMULATION OF THE STOKES EQUATIONS

An iterative adaptive equation solver for solving the implicit Stokes equations simultaneously with tri-tree grid generation is developed. The tri-tree grid generator builds a hierarchical grid structure which is mapped to a finite element grid at each hierarchical level. For each hierarchical finite element grid the Stokes equations are solved. The approximate solution at each level is projected onto the next finer grid and used as a start vector for the iterative equation solver at the finer level. When the finest grid is reached, the equation solver is iterated until a tolerated solution is reached. In order to reduce the overall work, the element matrices are integrated analytically beforehand. The efficiency and behaviour of the present adaptive method are compared with those of the previously developed iterative equation solver which is preconditioned by incomplete LU factorization with coupled node fill-in. The efficiency of the incomplete coupled node fill-in preconditioner is shown to be largely dependent on the global node numbering. The preconditioner is therefore tested for the natural node ordering of the tri-tree grid generator and for different ways of sorting the nodes.     

NUMERICAL SIMULATION OF UNSTEADY VISCOUS FLOWS USING AN IMPLICIT PROJECTION METHOD

In this paper an implicit projection method for the solution of the two-dimensional, time-dependent, incompressible Navier– Stokes equations is presented. The basic principle of this method is that the evaluation of the time evolution is split into intermediate steps. The computational method is based on the approximate factorization technique. The coupled approach is used to link the equations of motion and the turbulence model equations. The standard k-ϵ turbulence model is used. The current methodology, which has been tested extensively for steady problems, is now applied for the numerical simulation of unsteady flows. Several cases were tested, such as plane or axisymmetric channels, a backward-facing step, a square cavity and an axisymmetric stenosis.     

An adaptive finite element method for forced convection

This paper presents an adaptive finite element method to solve forced convective heat transfer. Solutions are obtained in primitive variables using a high-order finite element approximation on unstructured grids. Two general-purpose error estimators are developed to analyse finite element solutions and to determine the characteristics of an improved mesh which is adaptively regenerated by the advancing front method. The adaptive methodology is validated on a problem with a known analytical solution. The methodology is then applied to heat transfer predictions for two cases of practical interest. Predictions of the Nusselt number compare well with measurements and constitute an improvement over previous results. © 1997 John Wiley & Sons, Ltd.     

比奥固结有限元中参数正交敏感性分析

用邓肯-张模型求解比奥平面固结有限元涉及到的参数较多,本文通过正交试验研究参数共同随机变化对软土地基固结沉降的影响,发现其中的k、c、G、φ和F对整个固结沉降过程的影响是高度显著的; Pc在固结的初期影响是显著的;在其他参数变异较大时,Rf对固结沉降的影响相对较小,常渗透系数Kz在固结初期影响显著。本文结果对模型参数的取值、反演、随机分析以及固结沉降的估计具有参考价值。     

NUMERICAL SOLUTIONS OF OPTIMAL CONTROL FOR THERMALLY CONVECTIVE FLOWS

We study the numerical solution of optimal control problems associated with two-dimensional viscous incompressible thermally convective flows. Although the techniques apply to more general settings, the presentation is confined to the objectives of minimizing the vorticity in the steady state case and tracking the velocity field in the non-stationary case with boundary temperature controls. In the steady state case we develop a systematic way to use the Lagrange multiplier rules to derive an optimality system of equations from which an optimal solution can be computed; finite element methods are used to find approximate solutions for the optimality system of equations. In the time-dependent case a piecewise-in-time optimal control approach is proposed and the fully discrete approximation algorithm for solving the piecewise optimal control problem is defined. Numerical results are presented for both the steady state and time-dependent optimal control problems. © 1997 John Wiley & Sons, Ltd.     

HIGH-ORDER-ACCURATE SCHEMES FOR INCOMPRESSIBLE VISCOUS FLOW

We present new finite difference schemes for the incompressible Navier–Stokes equations. The schemes are based on two spatial differencing methods; one is fourth-order-accurate and the other is sixth-order accurate. The temporal differencing is based on backward differencing formulae. The schemes use non-staggered grids and satisfy regularity estimates, guaranteeing smoothness of the solutions. The schemes are computationally efficient. Computational results demonstrating the accuracy are presented. © 1997 by John Wiley & Sons, Ltd.     

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.