三角翼上分离及涡流的数值模拟

综述了三角翼(包括双三角翼,边条-三角翼,近距耦合鸭翼)上分离及涡流问题的数值模拟进展,介绍了用Euler方程、N-S方程的数值方法模拟不同三角翼,在不同攻角、不同来流M数等多种条件下的复杂涡流形态和流动结构,研究了控制方程、网格、湍流模型和计算方法等对计算结果的影响。      

Phase separation in impacting wyes and tees

An experimental and analytical study of phase separation for impacting two-phase flows in branching conduits has been performed. The resulting analytical model is applicable to impacting flows in wyes and tees for various inlet flow regimes. This model is based on a dividing-streamline concept, and it assumes that there is a “zone of influence” for each of the two phases which is bounded by the conduit wall and the dividing streamlines. Good agreement was obtained between model predictions and the experimental results from this study and others.     

Eulerian prediction of the fluid/particle correlated motion in turbulent two-phase flows

An approximate equation governing the turbulent fluid velocity encountered along discrete particle path is used to derive the fluid/particle turbulent moments required for dispersed two-phase flows modelling. Then, closure model predictions are compared with results obtained from large-eddy simulation of particle fluctuating motion in forced isotropic fluid turbulence.     

An algorithm for rotating axisymmetric flows: model,validation and industrial applications

The study of axisymmetric flows is of interest not only from an academic point of view, due to the existence of exact solutions of Navier–Stokes equations, but also from an industrial point of view, since these kind of flows are frequently found in several applications. In the present work the development and implementation of a finite element algorithm to solve Navier–Stokes equations with axisymmetric geometry and boundary conditions is presented. Such algorithm allows the simulation of flows with tangential velocity, including free surface flows, for both laminar and turbulent conditions. Pseudo‐concentration technique is used to model the free surface (or the interface between two fluids) and the k–ε model is employed to take into account turbulent effects. The finite element model is validated by comparisons with analytical solutions of Navier–Stokes equations and experimental measurements. Two different industrial applications are presented. Copyright © 2005 John Wiley & Sons, Ltd.     

A practical implementation of turbulence models for the computation of three-dimensional separated flows

An upwind MUSCL-type implicit scheme for the three-dimensional Navier-Stokes equations is presented and details on the implementation for three-dimensional flows of a ‘diagonal’ upwind implicit operator are developed. Turbulence models for separated flows are also described with an emphasis on the numerical specificities of the Johnson-King non-equilibrium model. Good predictions of separated two- and three-dimensional flows are demonstrated.     

An improved large eddy simulation of two-phase flows in a pump impeller

An improved large eddy simulation using a dynamic second-order sub-grid-scale (SGS) stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the current second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body-fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbulent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results. The project supported by the National Natural Science Foundation of China (50779069 and 90510007), the Start-up Scientific Research Foundation of China Agricultural University (2006021) and the Beijing Natural Science Foundation (3071002).     

NUMERICAL PREDICTION OF TWO-PHASE FLOW IN BUBBLE COLUMNS

A numerical model is described for the prediction of turbulent continuum equations for two-phase gas–liquid flows in bubble columns. The mathematical formulation is based on the solution of each phase. The two-phase model incorporates interfacial models of momentum transfer to account for the effects of virtual mass, lift, drag and pressure discontinuities at the gas–liquid interface. Turbulence is represented by means of a two-equation k–ϵ model modified to account for bubble-induced turbulence production. The numerical discretization is based on a staggered finite-volume approach, and the coupled equations are solved in a segregated manner using the IPSA method. The model is implemented generally in the multipurpose PHOENICS computer code, although the present appllications are restricted to two-dimensional flows. The model is applied to simulate two bubble column geometries and the predictions are compared with the measured circulation patterns and void fraction distributions.     

An unstructured finite element model for incompressible two-phase flow based on a monolithic conservative level set method

We present a robust numerical method for solving incompressible, immiscible two-phase flows. The method extends both a monolithic phase conservative level set method with embedded redistancing and a semi-implicit high-order projection scheme for variable-density flows. The level set method can be initialized conveniently via a simple phase indicator field instead of a signed distance function (SDF). To process the indicator field into a SDF, we propose a new partial differential equation-based redistancing method. We also improve the monolithic level set scheme to provide more accuracy and robustness in full two-phase flow simulations. Specifically, we perform an extra step to ensure convergence to the signed distance level set function and simplify other aspects of the original scheme. Lastly, we introduce consistent artificial viscosity to stabilize the momentum equations in the context of the projection scheme. This stabilization is algebraic, has no tunable parameters and is suitable for unstructured meshes and arbitrary refinement levels. The overall methodology includes few numerical tuning parameters; however, for the wide range of problems that we solve, we identify only one parameter that strongly affects performance of the computational model and provide a value that provides accurate results across all the benchmarks presented. This methodology results in a robust, accurate, and efficient two-phase flow model, which is mass- and volume-conserving on unstructured meshes and has low user input requirements, making it attractive for real-world applications.     

Numerical simulation of the vertical structure of discontinuous flows

A numerical method to solve the Reynolds‐averaged Navier–Stokes equations with the presence of discontinuities is outlined and discussed. The pressure is decomposed into the sum of a hydrostatic component and a hydrodynamic component. The numerical technique is based upon the classical staggered grids and semi‐implicit finite difference methods applied for quasi‐ and non‐hydrostatic flows. The advection terms in the momentum equations are approximated in order to conserve mass and momentum following the principles recently developed for the numerical simulation of shallow water flows with large gradients. Conservation of these properties is the most important aspect to represent near local discontinuities in the solution, following from sharp bottom gradients or hydraulic jumps. The model is applied to reproduce the flow over a step where a hydraulic jump forms downstream. The hydrostatic pressure assumption fails to represent this type of flow mainly because of the pressure deviation from the hydrostatic values downstream the step. Fairly accurate results are obtained from the numerical model compared with experimental data. Deviation from the data is found to be inherent to the standard k–ε model implemented. Copyright © 2001 John Wiley & Sons, Ltd.     

An implicit Lagrangian lattice Boltzmann method for the compressible flows

In this paper, we propose a new Lagrangian lattice Boltzmann method (LBM) for simulating the compressible flows. The new scheme simulates fluid flows based on the displacement distribution functions. The compressible flows, such as shock waves and contact discontinuities are modelled by using Lagrangian LBM. In this model, we select the element in the Lagrangian coordinate to satisfy the basic fluid laws. This model is a simpler version than the corresponding Eulerian coordinates, because the convection term of the Euler equations disappears. The numerical simulations conform to classical results. Copyright © 2006 John Wiley & Sons, Ltd.     

The setting up of the equations for two-phase flows by the kinetic theory method

The fundamental equations for two-phase flows are deduced from the Boltzmann's equation. The collision terms are treated with a method similar to what is used in the classical kinetic theory for handling the transport properties of dense gases. It is shown that collision pressure and collision thermal flux exist in gas-particle flows in addition to the general partial pressure and partial thermal flux. Their physical natures are quite different from those of the general partial pressure and partial thermal flux. The applicability of the binary collision assumption and the molecular chaos assumption to gas-particle flows is also discussed. Finally, the equations for two-phase flows obtained by the method of the kinetic theory are compared with those obtained by average continuum models and by the model of particle clouds. The results from the kinetic theory show clearly the physical significance of various parameters and clarify some confusing concepts. Institute of Mechanics, Academia Sinica     

高雷诺数流动的基本控制方程体系和扩散跑物化Navier-Stokes方程组的意义和用途

在计算机发达的时代, 高雷诺($Re$)数绕流计算中有无必要使用简化NS方程组, 本文讨论这个问题. 主要内容如下: (1)高$Re$数绕流包含3种基本流动: 所有方向对流占优流动、所有方向对流扩散竞争流动和部分方向对流占优部分方向对流扩散竞争流动(简称干扰剪切流动), 3个基本流动的特征彼此不同且在流场中所占领域大小彼此相差悬殊, NS方程区域很小,它们的最简单控制方程组Euler、Navier-Stokes (NS)和扩散抛物化(DP) NS方程组的数学性质彼此不同, 因此利用Euler-DPNS-NS方程组体系分析计算高$Re$数绕流流动就是一个合乎逻辑的选择, 该法与利用单一NS方程组的常用方法可以彼此检验和补充. (2)流体之间以及流体与外界的动量、能量和质量交换, 流态从层流到湍流的演化主要发生在干扰剪切流动中, 干扰剪切流及其最简单控制方程------DPNS方程组具有基础意义; DPNS方程组笔者在1967年已提出. (3)诸简化NS方程组: DPNS、抛物化(P)NS、薄层(TL)NS、黏性层(VL)NS方程组的发展、相互关系, 它们的历史贡献和今后的用途; 它们的数学性质均为扩散抛物型, 但它们包含的黏性项彼此有所不同; 从流体力学角度来看, 它们中只有DPNS方程组能够准确描述干扰剪切流动. 提出把诸简化NS方程组统一为DPNS方程组的建议. (4)干扰剪切流------DPNS方程组与无干扰剪切流------边界层方程组之间的关系以及进一步研究干扰剪切流的意义.      

一维溃坝洪水波的高精度数值模拟

将ENO(Essentially Non-Oscillatory)格式和Runge-Kutta时间离散的思想应用于一维Saint-Venant方程组的求解,数值模拟溃坝洪水,得出了水位和流速的沿程分布。经与理论解比较可见,数值解在间断波附近没有出现数值振荡,水位和流速大小均符合较好,表明ENO格式是一类新的高精度无振荡差分格式,采用ENO格式所建立的高分辨率模型能够很好地模拟溃坝波的演进过程。     

APPLICATION OF FLUX-VECTOR-SPLITTING SCHEME TO A DILUTE GAS–PARTICLE JPL NOZZLE FLOW

A time-dependent numerical algorithm is developed for the two-fluid model Euler of TLNS (thin layer Navier–Stokes) equations. The analysis is based on a MUSCL (monotone upstream central scheme for conservation laws)-type flux-vector-splitting scheme with the multi-level technique. This algorithm is applied to investigate JPL (Jet Propulsion Laboratory) nozzle flow. Calculated results for both one- and two-phase flows are given to show the accuracy, the computational efficiency and the particle influence on the flow field.     

A PRESSURE-BASED ALGORITHM FOR HIGH-SPEED TURBOMACHINERY FLOWS

The steady state Navier–Stokes equations are solved in transonic flows using an elliptic formulation. A segregated solution algorithm is established in which the pressure correction equation is utilized to enforce the divergence-free mass flux constraint. The momentum equations are solved in terms of the primitive variables, while the pressure correction field is used to update both the convecting mass flux components and the pressure itself. The velocity components are deduced from the corrected mass fluxes on the basis of an upwind-biased density, which is a mechanism capable of overcoming the ellipticity of the system of equations, in the transonic flow regime. An incomplete LU decomposition is used for the solution of the transport-type equations and a globally minimized residual method resolves the pressure correction equation. Turbulence is resolved through the k–ε model. Dealing with turbomachinery applications, results are presented in two-dimensional compressor and turbine cascades under design and off-design conditions. © 1997 John Wiley & Sons, Ltd.     

A two-scale second-moment turbulence closure based on weighted spectrum integration

A two-scale second-moment turbulence closure has been derived based on the weighted integration of the dynamic equation for the covariance spectrum. The goal is to close the Reynolds stress equations with two additional scalar equations that provide separately the scales of the spectral energy transfer and of the turbulence energy dissipation rate. Such a model should provide better prediction of nonequilibrium turbulent flows. The derivation consists of analytical integration of the wave-number-weighted covariance spectrum using a model of the spectral equations with an assumed simple representation of the shape of the energy spectrum. The resulting closure consists of a set of three tensorial equations, one for the Reynolds stress and two for length scale tensors, the latter representing the energy containing- and dissipative eddies respectively. The trace of the two tensor-scale equations leads to a set of two scalar scale parameters. In the equilibrium limit, the model reduces to the standard second-moment single-scale closure. The approach makes it also possible to derive the scale equations in a more systematic manner as compared with the common single-scale and other multi-scale models. The performance of the model in capturing the scale dynamics is illustrated by predictions of several generic homogeneous and inhomogeneous unsteady flows, demonstrating the expected response of the two scale equations. PACS 03.50.De, 04.20-q, 42.65-k     

Computation of 3D vortex flows past a flat plate at incidence through a variational approach of the full steady euler equations

A variational method for solving directly the full steady Euler equations is presented. This method is based on both Newton's linearization and a least squares formulation. The validity of the Euler model and boundary conditions to capture the vortex sheet is discussed. A finite element approximation of the groups of conservative variables is described and results are given for 3D subsonic flows as well as supersonic flows past a flat plate at high angle of attack.     

Transition of Stratified Flow in Horizontal or Slightly Inclined Pipe with Elliptical Cross Section

The study concerns the separated two-phase gas-liquid flow in a horizontal or slightly inclined pipe. The cross section of the cylindrical pipe is elliptical. After deriving the partial differential equations, the analytical prediction of the transition between stratified/non-stratified flows is obtained by using a linear stability analysis. The stratified flow, in a pipe having one of its cross sections axes disposed horizontally, is more stable when this axis is the great major.     

Advances in the study of separated flows

This paper consists of three parts. The first deals with the separation conditions for three dimensional steady viscous separated flows, in which the behaviour of separated flow described by Navier-Stokes equations and boundary layer equations is studied. The second part involves an application of differential topology to qualitative analysis of flow fields. Here the distribution rule of singular points on the separation line is studied. The last part discusses the numerical method solving Navier-Stokes equations for separated flows. The obtained computational results are analysed by the above mentioned theories and methods.