A new method for accelerating the rate of convergence of the simple-like algorithm

A pressure correction formula is proposed for the SIMPLE-like algorithm in order to improve the rate of the convergence when solving laminar Navier–Stokes equations when there is rapidly varying pressure. Based on global mass conservation, a line average pressure correction is derived by integration of the momentum equation for approximate one-dimensional flow. The use of this formula with the SIMPLE-like algorithm can rapidly build up the pressure distribution in the region where the pressure undergoes a very large change, which normally causes the rate of convergence of the SIMPLE or the SIMPLEC schemes to be slow. In order to illustrate the technique, the performances of SIMPLE and of SIMPLEC with the average pressur correction are investigated for axisymmetric flow past and through a sampler. A comparison of these two techniques shows that the average pressure correction proposed in this paper significantly accelerates the rate of convergence.     

A robust numerical method for flow through a pipe driven by an oscillating pressure gradient

The problem of periodic flow of an incompressible fluid through a pipe, which is driven by an oscillating pressure gradient (e.g. a reciprocating piston), is investigated in the case of a large Reynolds number. This process is described by a singularly perturbed parabolic equation with a periodic right‐hand side, where the singular perturbation parameter is the viscosity ν. The periodic solution of this problem is a solution of the Navier–Stokes equations with cylindrical symmetry. We are interested in constructing a parameter‐robust numerical method for this problem, i.e. a numerical method generating numerical approximations that converge uniformly with respect to the parameter ν and require a bounded time, independent of the value of ν, for their computation. Our method comprises a standard monotone discretization of the problem on non‐standard piecewise uniform meshes condensing in a neighbourhood of the boundary layer. The transition point between segments of the mesh with different step sizes is chosen in accordance with the behaviour of the analytic solution in the boundary layer region. In this paper we construct the numerical method and discuss the results of extensive numerical experiments, which show experimentally that the method is parameter‐robust. Copyright © 2006 John Wiley & Sons, Ltd.     

A Petrov-Galerkin method for the numerical solution of the Bradshaw-Ferriss-Atwell turbulence model

The Bradshaw-Ferriss-Atwell model for 2D constant property turbulent boundary layers is shown to be ill-posed with respect to numerical solution. It is shown that a simple modification to the model equations results in a well-posed system which is hyperbolic in nature. For this modified system a numerical algorithm is constructed by discretizing in space using the Petrov-Galerkin technique (of which the standard Galerkin method is a special case) and stepping in the timelike direction with the trapezoidal (Crank-Nicolson) rule. The algorithm is applied to a selection of test problems. It is found that the solutions produced by the standard Galerkin method exhibit oscillations. It is further shown that these oscillations may be eliminated by employing the Petrov-Galerkin method with the free parameters set to simple functions of the eigenvalues of the modified system.     

An implicit numerical algorithm for solving non‐hydrostatic free‐surface flow problems

Details are given of the development of a two‐dimensional vertical numerical model for simulating unsteady free‐surface flows, using a non‐hydrostatic pressure distribution. In this model, the Reynolds equations and the kinematic free‐surface boundary condition are solved simultaneously, so that the water surface elevation can be integrated into the solution and solved for, together with the velocity and pressure fields. An efficient numerical algorithm has been developed, deploying implicit parameters similar to those used in the Crank–Nicholson method, and generating a block tri‐diagonal algebraic system of equations. The model has been applied to simulate a range of unsteady flow problems involving relatively strong vertical accelerations. The results show that the numerical algorithm described is able to produce accurate predictions and is also easy to apply. Copyright © 2001 John Wiley & Sons, Ltd.     

Direct numerical simulation of a turbulent plane Couette-Poiseuille flow with zero-mean shear

Direct numerical simulations of a turbulent Couette-Poiseuille flow with zero-mean-shear at the moving wall (SL-flow) is performed to examine flow features compared to those for a turbulent pure Poiseuille flow (P-flow). Profiles of the streamwise mean velocity, indicator function and ratio of production to dissipation show that the logarithmic region is significantly elongated for the SL-flow compared to that for the P-flow at a similar Reynolds number. In addition, the magnitudes of the Reynolds stresses are found to be larger in both inner and outer layers for the SL-flow than those for the P-flow. The spanwise spectra of the production term in the turbulent kinetic energy equation are examined to provide a structural basis for explaining the statistical behaviors. In addition, because the growth of the energy-containing motions extends to the outer layer further for the SL-flow due to the presence of a positive mean shear throughout the entire wall layer, the self-similar behavior of the energy balance between the production and transport terms with respect to the self-similar wavenumber is found far from the wall. We also find the increase in the number of uniform momentum zones in the SL-flow, revealing the hierarchical distribution of the energy-containing eddies which are composed of multiple uniform momentum zones. These coherent motions lead to the elongation of the logarithmic region for the SL-flow. Finally, investigation of the turbulent energy transfer process in a spectral domain for the SL-flow demonstrates importance of outer layer very-long structures, and these structures attribute to the energy transport process in an entire flow field.     

A decoupling numerical method for fluid flow

A first-order non-conforming numerical methodology, Separation method, for fluid flow problems with a 3-point exponential interpolation scheme has been developed. The flow problem is decoupled into multiple one-dimensional subproblems and assembled to form the solutions. A fully staggered grid and a conservational domain centred at the node of interest make the decoupling scheme first-order-accurate. The discretization of each one-dimensional subproblem is based on a 3-point interpolation function and a conservational domain centred at the node of interest. The proposed scheme gives a guaranteed first-order accuracy. It is shown that the traditional upwind (or exponentially weighted upstream) scheme is less than first-order-accurate. The pressure is decoupled from the velocity field using the pressure correction method of SIMPLE. Thomas algorithm (tri-diagonal solver) is used to solve the algebraic equations iteratively. The numerical advantage of the proposed scheme is tested for laminar fluid flows in a torus and in a square-driven cavity. The convergence rates are compared with the traditional schemes for the square-driven cavity problem. Good behaviour of the proposed scheme is ascertained.     

An asymptotic iteration method for the numerical analysis of near-critical free-surface flows

Satisfying the boundary conditions at the free surface may impose severe difficulties to the computation of turbulent open-channel flows with finite-volume or finite-element methods, in particular, when the flow conditions are nearly critical. It is proposed to apply an iteration procedure that is based on an asymptotic expansion for large Reynolds numbers and Froude numbers close to the critical value 1.The iteration procedure starts by prescribing a first approximation for the free surface as it is obtained from solving an ODE that has been derived previously by means of an asymptotic expansion (Grillhofer and Schneider, 2003). The numerical solution of the full equations of motion then gives a surface pressure distribution that differs from the constant value required by the dynamic boundary condition. To determine a correction to the elevation of the free surface we next solve an ODE that is obtained from the asymptotic analysis of the flow with a prescribed pressure disturbance at the free surface. The full equations of motion are then solved for the corrected surface, and the procedure is repeated until criteria of accuracy for surface elevation and surface pressure, respectively, are satisfied.The method is applied to an undular hydraulic jump as a test case.     

A numerical study of flow through wavy-walled channels

A numerical procedure is developed for the analysis of flow in a channel whose walls describe a travelling wave motion. Following a perturbation method, the primitive variables are expanded in a series with the wall amplitude as the perturbation parameter. The boundary conditions are applied at the mean surface of the channel and the first-order perturbation quantities are calculated using the pseudospectral collocation method. Although limited by the linear analysis, the present approach is not restricted by the Reynolds number of the flow and the wave number and frequency of the wavy-walled channel. Using the computed wall shear stresses, the positions of flow separation and reattachment are determined. The variations in velocity and pressure with frequency of excitation are also presented. © 1998 John Wiley & Sons, Ltd.     

A scalable parallel algorithm for the direct numerical simulation of three-dimensional incompressible particulate flow

Particulate flow is of great importance from both the scientific and engineering points of view. Owing to the complexity of particle-flow interactions, direct numerical simulations (DNS) of inertial particulate flow with finite-size particles have been limited to a very small number of particles, while the industrial applications involve larger numbers with many orders of magnitude. This article presents a parallel implementation of a fictitious domain method for the DNS of particulate flows. The method is thoroughly tested and its parallel performance on distributed memory clusters is evaluated on a large-scale problem. Finally, we present the results for the separation of 21,336 particles of two different densities in a viscous fluid. Although there is still a significant gap between DNS and the industrial applications, the present algorithm allows to simulate significantly large number of particles so that a meaningful statistical analysis can be performed. This will help in the development of new closure relations for the averaged models of multiphase flows.     

AN IMPLICIT ALGORITHM OF THIN LAYER EQUATIONS IN VISCOUS，TRANSONIC，TWO－PHASE NOZZLE FLOW

ＡＮＩＭＰＬＩＣＴＡＬＧＯＲＩＴＨＭＯＦＴＨＩＮＬＡＹＥＲＥＱＵＡＴＩＯＮＳｉＮＶＩＳＣＯＵＳ,ＴＲＡＮＳＯＮＩＣ,ＴＷＯ－ＰＨＡＳＥＮＯＺＺＬＥＦＬＯＷＨｅＨｏｎｇ－ｑｉｎｇ（何洪庆）ＨｏｕＸｉａｏ（侯晓）ＣａｉＴｉ－ｍｉｎ（蔡体敏）ＷｕＸｉｎｇ－...     

基于萤火虫群优化的空间非合作目标相对导航粒子滤波算法

研究了空间非合作目标相对导航算法,针对标准粒子滤波的重采样过程导致的粒子贫化现象及其造成的相对导航精度下降问题,分析了萤火虫优化算法的运行机制,提出一种基于萤火虫智能优化算法的改进粒子滤波算法。改进算法通过优化粒子滤波的重采样过程,使粒子群智能的向高似然区域移动,同时在低似然区域也合理保留了部分粒子,保证了粒子的多样性,提高了样本的整体质量。仿真结果表明,改进算法导航精度较标准算法提高了39.35%,达到稳定精度所需粒子数较少,有效抑制了粒子贫化问题。     

A pressure correction method for fluid–particle interaction flow: Direct‐forcing method and sedimentation flow

A direct‐forcing pressure correction method is developed to simulate fluid–particle interaction problems. In this paper, the sedimentation flow is investigated. This method uses a pressure correction method to solve incompressible flow fields. A direct‐forcing method is introduced to capture the particle motions. It is found that the direct‐forcing method can also be served as a wall‐boundary condition. By applying Gauss's divergence theorem, the formulas for computing the hydrodynamic force and torque acting on the particle from flows are derived from the volume integral of the particle instead of the particle surface. The order of accuracy of the present method is demonstrated by the errors of velocity, pressure, and wall stress. To demonstrate the efficiency and capability of the present method, sedimentations of many spherical particles in an enclosure are simulated. Copyright © 2010 John Wiley & Sons, Ltd.     

A Helmholtz pressure equation method for the calculation of unsteady incompressibl eviscous flows

A time-implicit numerical method for solving unsteady incompressible viscous flow problems is introduced. The method is based on introducing intermediate compressibility into a projection scheme to obtain a Helmholtz equation for a pressure-type variable. The intermediate compressibility increases the diagonal dominance of the discretized pressure equation so that the Helmholtz pressure equation is relatively easy to solve numerically. The Helmholtz pressure equation provides an iterative method for satisfying the continuity equation for time-implicit Navier–Stokes algorithms. An iterative scheme is used to simultaneously satisfy, within a given tolerance, the velocity divergence-free condition and momentum equations at each time step. Collocated primitive variables on a non-staggered finite difference mesh are used. The method is applied to an unsteady Taylor problem and unsteady laminar flow past a circular cylinder.     

剧变截面圆管内渗流的数值计算方法

对于剧变截面圆管的渗流问题写出不可压缩渗流的基本方程组,对直接求解原始变量(速度和压力)的数值计算方法作出改进。先由非主流方向的运动方程计算压力,后由主流方向的运动方程计算主流方向的速度分量,再由连续性方程计算非主流方向的速度分量。这样可以避免在一般的求解原始变量方法中由连续性方程计算压力时出现的困难和麻烦。根据本方法和剧变截面圆管的特点,采用半交错不等距非正交贴体混合网格系。本文详细写出差分方程和迭代计算公式,对剧变截面圆管内的渗流算例进行数值计算。本方法的优点是简单和实用,在工程上具有较大的应用价值。     

A Taylor–Galerkin-based algorithm for viscous incompressible flow

In this paper the development and behaviour of a new finite element algorithm for viscous incompressible flow is presented. The stability and background theory are discussed and the numerical performance is considered for some benchmark problems. The Taylor–Galerkin approach naturally leads to a time-stepping algorithm which is shown to perform well for a wide range of Reynolds numbers (1 ? Re ? 400).

1 A conventional definition for Re is assumed.

Various modifications to the algorithm are investigated, particularly with respect to their effects on stability and accuracy.     

Characteristics of the critical heat flux for downward flow in a vertical tube at low flow rate and low pressure conditions

The characteristics of the critical heat flux (CHF) for downward flow were studied experimentally with an Inconel 600 circular tube test section in a water test loop at low-flow rate (0 200 kg/m²s) and low-pressure (0.1 0.7 MPa) conditions. The attention was given to the effects of upstream conditions—upper plenum and inlet throttling. Two totally different kinds of CHF behaviors were observed. It seems appropriate to interpret them as flooding-type CHF and dryout in annular flow. The CHF in downward flow may vary from extremely unstable flow CHF as low as near the flooding CHF value to stable flow CHF as high as that of upflow, depending on the upstream conditions of the test section. The CHF correlation by Mishima and that by Weber were proposed for the presentation of the lower and upper limits of the CHF for downward flow in a vertical tube at low-flow rate and low-pressure conditions.     

A concise method for determining a valve flow coefficient of a valve under compressible gas flow

This work presents a novel method of determining a valve flow coefficient for a valve and transient mass flow rate of compressible gas discharged from a reservoir. The proposed method consists of a set of equations to express the physical phenomena and measurement equipment to measure the indispensable data used in the above equations. Regardless of the kind of valve, the valve flow coefficient can be obtained efficiently and feasibly. The results of this study indicate not only that the valve flow characteristics of the diaphragm valve significantly differ from those of the ball valve, but also that the C_v flow equation conventionally used is no longer valid for the diaphragm valve. The valve flow coefficient of the ball valve determined by the proposed method is about 48 and the representative one proposed by the ANSI/ISA is about 56. In addition, the mass flow rate of gas flow through a valve under transient process can be estimated without using a flow meter. Moreover, the cumulated masses discharged predicted by the method proposed herein are consistent well with those of the experimental results. The deviations are smaller than 6%.     

Investigation of scale effect for the computation of turbulent flow around a circular cylinder

In order to investigate the scale effect of turbulent flow around a circular cylinder, two similarity numbers （criteria） based on turbulent kinetic and dissipation rates associ- ated with the fluctuation characteristics of turbulence wake are deduced by analyzing the Reynolds averaged NavierStokes equations （RANS）. The RNG k-s models and finite volume method are used to solve the governing equations and the second-order implicit time and upwind space discretization algorithms are used to discrete the governing equations. A numerical computation of flow parameters around a two-dimensional circular cylinder with Reynolds numbers ranging from 102 to l07 is accomplished and the result indicates that the fluctuation of turbulence flow along the center line in the wake of circular cylinder can never be changed with increasing Reynolds numbers when Re ≥ 3 × 10^6. This conclusion is useful for controlling the scale of numerical calculations and for applying model test data to engineering practice.     

A finite volume method for two-dimensional transonic potential flow through turbomachinery blade rows

To predict inviscid transonic flow through turbomachinery blade rows, the exact transonic potential flow equation is solved on a mesh constructed from small area elements. A transformation is introduced through which distorted squares of the physical plane are mapped into computational squares. Two sets of overlapping elements are used; while the thermodynamic properties are calculated at the primary element centres, the flux balance is established on the secondary elements. For transonic flows an artificial compressibility term (upwind density gradient) is added to density in order to produce the desired directional bias in the hyperbolic region. while the entropy does not increase across mass conservative shock jump regions. Comparisons withexperiments and with other numerical and analytical solutions for various turbomachinery configurations show that this approach is comparatively accurate, reliable, and fast.