The flowfield characteristics of a precombustion chamber of coflowing jets with large velocity differences

A new type of precombustion chamber has been developed in which coflowing jets with large velocity differences are used. A finite difference calculation procedure based on the SIMPLE algorithm is developed for the numerical simulation of this flowfield. This paper presents a statement of the algorithm, boundary condition specifications, and the results of computational experiments. The computational results demonstrate the existence of the large recirculation regions. It displays significant three-dimensional characteristics. The effects of design parameters on the recirculation features are discussed briefly in the present study.     

Poiseuille-Bénard流的出口边界条件及其谱元法计算

研究二维矩形管道中底部加热的不可压缩Ｐｏｉｓｅｕｉｌｌｅ－Ｂｅｎａｒｄ流的谱元法数值计算问题．讨论各种不同的出口边界条件的处理及其对谱元法数值模拟的影响．通过干扰区、干扰幅度和计算时间的比较,确定比较理想的出口边界条件．     

Spectral element method for viscoelastic flows in a planar contraction channel

A new algorithm, which combines the spectral element method with elastic viscous splitting stress (EVSS) method, has been developed for viscoelastic fluid flows in a planar contraction channel. The system of spectral element approximations to the velocity, pressure, extra stress and the rate of deformation variables is solved by a preconditioned conjugate gradient method based on the Uzawa iteration procedure. The numerical approach is implemented on a planar four‐to‐one contraction channel for a fluid governed by an Oldroyd‐B constitutive equation. The behaviour of the Oldroyd‐B fluids in the contraction channel is investigated with various Weissenberg numbers. It is shown that numerical solutions obtained here agree well with experimental measurements and other numerical predictions. Copyright © 2003 John Wiley & Sons, Ltd.     

The characteristic‐based split (CBS) meshfree method for free surface flow problems in ALE formulation

A semi‐implicit characteristic‐based split (CBS) meshfree algorithm in the arbitrary Lagrangian Eulerian (ALE) framework is proposed for the numerical solution of incompressible free surface flow problem in the paper. The algorithm is the extension of general CBS method which was initially introduced in finite element framework, this is due to the fact that CBS method not only can enhance the stability, but also avoid LBB condition when equal order basis function is used to approximate velocity and pressure variables. Meanwhile, a simple way for node update and node speed calculation is developed which is used to capture the free surface exactly. The numerical solutions are compared with available analytical and numerical solutions, which shows that the proposed method has better ability to simulate the free surface incompressible flow problem. Copyright © 2009 John Wiley & Sons, Ltd.     

The structure of a three-dimensional boundary layer subjected to streamwise-varying spanwise-homogeneous pressure gradient

A spatially-evolving three-dimensional boundary layer, subjected to a streamwise-varying spanwise-homogeneous pressure gradient, equivalent to a body force, is investigated by way of direct numerical simulation. The pressure gradient, prescribed to change its sign half-way along the boundary layer, provokes strong skewing of the velocity vector, with a layer of nearly collateral flow forming close to the wall up to the position of maximum spanwise velocity. A wide range of flow-physical properties have been studied, with particular emphasis on the near-wall layer, including second-moments, major budget contributions and wall-normal two-point correlations of velocity fluctuations and their angles, relative to wall-shear fluctuations. The results illustrate the complexity caused by skewing, including a damping in turbulent mixing and a significant lag between strains and stresses. The study has been undertaken in the context of efforts to develop and test novel hybrid LES–RANS schemes for non-equilibrium near-wall flows, with an emphasis on three-dimensional near-wall straining. Fundamental flow-physical issues aside, the data derived should be of particular relevance to a priori studies of second-moment RANS closure and the development and validation of RANS-type near-wall approximations implemented in LES schemes for high-Reynolds-number complex flows.     

THREE-DIMENSIONAL VISCOUS FLOW THROUGH A ROTATING CHANNEL: A PSEUDOSPECTRAL MATRIX METHOD APPROACH

A Fourier–Chebyshev pseudospectral method is used for the numerical simulation of incompressible flows in a three-dimensio nal channel of square cross-section with rotation. Realistic, non-periodic boundary conditions that impose no-slip conditions in two directions (spanwis e and vertical directions) are used. The Navier–Stokes equations are integrated in time using a fractional step method. The Poisson equations for pressure and the Helmholtz equation for velocity are solved using a matrix diagonalization (eigenfunction decomposition) method, through which we are able to reduce a three-dimensional matrix problem to a simple algebraic vector equation. This results in signficant savings in computer storage requirement, particularly for large-scale computations. Verification of the numerical algorithm and code is carried out by comparing with a limiting case of an exact steady state solution for a one-dimensional channel flow and also with a two-dimensional rotating channel case. Two-cell and four-cell two-dimensional flow patterns are observed in the numerical experiment. It is found that the four-cell flow pattern is stable to symmetri cal disturbances but unstable to asymmetrical disturbances.     

Unsteady PTV velocity field past an airfoil solved with DNS: Part 1. Algorithm of hybrid simulation and hybrid velocity field at <Emphasis Type="Italic">Re</Emphasis> ≈ 10<Superscript>3</Superscript>

We develop a hybrid unsteady-flow simulation technique combining direct numerical simulation (DNS) and particle tracking velocimetry (PTV) and demonstrate its capabilities by investigating flows past an airfoil. We rectify instantaneous PTV velocity fields in a least-squares sense so that they satisfy the equation of continuity, and feed them to the DNS by equating the computational time step with the frame rate of the time-resolved PTV system. As a result, we can reconstruct unsteady velocity fields that satisfy the governing equations based on experimental data, with the resolution comparable to numerical simulation. In addition, unsteady pressure distribution can be solved simultaneously. In this study, particle velocities are acquired on a laser-light sheet in a water tunnel, and unsteady flow fields are reconstructed with the hybrid algorithm solving the incompressible Navier–Stokes equations in two dimensions. By performing the hybrid simulation, we investigate nominally two-dimensional flows past the NACA0012 airfoil at low Reynolds numbers. In part 1, we introduce the algorithm of the proposed technique and discuss the characteristics of hybrid velocity fields. In particular, we focus on a vortex shedding phenomenon under a deep stall condition (α = 15°) at Reynolds numbers of Re = 1000 and 1300, and compare the hybrid velocity fields with those computed with two-dimensional DNS. In part 2, the extension to higher Reynolds numbers is considered. The accuracy of the hybrid simulation is evaluated by comparing with independent experimental results at various angles of attack and Reynolds numbers up to Re = 10⁴. The capabilities of the hybrid simulation are also compared with two-dimensional unsteady Reynolds-Averaged Navier–Stokes (URANS) solutions in part 2. In the first part of these twin papers, we demonstrate that the hybrid velocity field approaches the PTV velocity field over time. We find that intensive alternate vortex shedding past the airfoil, which is predicted by the two-dimensional DNS, is substantially suppressed in the hybrid simulation and the resultant flow field is similar to the PTV velocity field, which is projection of the three-dimensional velocity field on the streamwise plane. We attempt to identify the motion that originates three-dimensional flow patterns by highlighting the deviation of the PTV velocity field from the two-dimensional governing equations at each snapshot. The results indicate that the intensive spots of the deviation appear in the regions in which three-dimensional instabilities are induced in the shear layer separated from the pressure side.     

贯流风机流场的实验和数值研究

贯流风机流场的分析研究是对其进行结构调整、性能优化的前提。本文运用烟迹和闪光照相的方法对贯流风机的内部流动进行了流场显示研究，得到了瞬态流场照片，并通过多点总压和静压测量得到了速度分布。本文还对贯流风机的流场进行了数值模拟。结果表明：数值计算与实验结果符合得较好，这两者的结合是研究贯流风机流动和改进其性能的有效手段。本文最后就贯流风机的噪声进行了初步研究，得出了一些有价值的结果。     

Gas-kinetic numerical method for solving mesoscopic velocity distribution function equation

A gas-kinetic numerical method for directly solving the mesoscopic velocity distribution function equation is presented and applied to the study of three-dimensional complex flows and micro-channel flows covering various flow regimes. The unified velocity distribution function equation describing gas transport phenomena from rarefied transition to continuum flow regimes can be presented on the basis of the kinetic Boltzmann–Shakhov model equation. The gas-kinetic finite-difference schemes for the velocity distribution function are constructed by developing a discrete velocity ordinate method of gas kinetic theory and an unsteady time-splitting technique from computational fluid dynamics. Gas-kinetic boundary conditions and numerical modeling can be established by directly manipulating on the mesoscopic velocity distribution function. A new Gauss-type discrete velocity numerical integration method can be developed and adopted to attack complex flows with different Mach numbers. HPF parallel strategy suitable for the gas-kinetic numerical method is investigated and adopted to solve three-dimensional complex problems. High Mach number flows around three-dimensional bodies are computed preliminarily with massive scale parallel. It is noteworthy and of practical importance that the HPF parallel algorithm for solving three-dimensional complex problems can be effectively developed to cover various flow regimes. On the other hand, the gas-kinetic numerical method is extended and used to study micro-channel gas flows including the classical Couette flow, the Poiseuille- channel flow and pressure-driven gas flows in two-dimensional short micro-channels. The numerical experience shows that the gas-kinetic algorithm may be a powerful tool in the numerical simulation of micro-scale gas flows occuring in the Micro-Electro-Mechanical System (MEMS). The project supported by the National Natural Science Foundation of China (90205009 and 10321002), and the National Parallel Computing Center in Beijing. The English text was polished by Yunming Chen.     

Stability of a numerical algorithm for gas bubble modelling

A free-surface-tracking algorithm based on the SOLA-VOF method is analysed for numerical stability when modelling gas bubble evolution in a fluid. It is shown that an instability can arise from the fact that the bubble pressure varies with its volume. A time step stability criterion is introduced which is a function of the natural oscillation period but does not depend on the mesh size. This dependence suggests that the instability is likely to arise in the case of a general motion of a bubble, especially if break-up occurs. The effect is shown using linear Fourier analysis of the discretized equation for radial bubble oscillation and demonstrated numerically using a CFD code FLOW-3D. One- and three-dimensional situations are considered: a bubble in a fluid bounded by two concentric surfaces and a bubble floating in a fluid chamber with and without gravity. In cases where no analytical solution is available, a numerical method for the stability time step limit calculation is suggested based on finding the natural oscillation frequency. The nature of the instability suggests that it can be a feature of any numerical algorithm which models transient fluid flow with a free surface.     

NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE FLOW BY A NEW FORMULATION

A new mathematical formulation, called the pseudovorticity–velocity formulation, of the three-dimensional incompressible Navier–Stokes equations is presented as an alternative to the vorticity–velocity approach. For the model lid-driven cavity flow problem in two and three dimensions, combined with an explicit mixed spectral /finite different numerical scheme the proposed formulation is found to be efficient and very accurate as compared with the results available in the literature. In particular, the simulation results demonstrate an attractive feature of the present formulation compared with the vorticity–velocity approach, namely that the divergence-free condition of the velocity field can always be achieved on a non-staggered mesh.     

Examining a modified algorithm of smoothed particle hydrodynamics for a high velocity perforation of an aluminum beam

The wide spreading of utilizing of smoothed particle hydrodynamics (SPH) for numerical studies of the complex and high rate deformations of continuums, led the current study to gain a more reliable simulation by employing a modified compressible smoothed particle hydrodynamics (MCSPH) algorithm which could be a more accurate and stable technique in high tension regions, in despite of incompressible standard SPH. The main feature of the modified compressible SPH algorithm relies on a three steps solution procedure to calculate the pressure gradient, the deviatoric stress tensor, and the body forces separately. This algorithm is free of any artificial viscosity in its formulations, as well as welcoming to compressible effects which permits the pressure shock waves in high rate plastic deformation. To examine the accuracy of the algorithm, a benchmark problem of colliding rubber cylinders was simulated first and then a high velocity perforation process of an aluminum beam struck by a rigid projectile was simulated in various projectile speeds, and the failure response of the beam in each case was accompanied by crack propagation process. The prominent capability of the utilized MCSPH can be more illustrated when it was used in simulation of thickness crack propagation a tiny crack paths and defragmentation which can be encountered as a not easy numerical case study. The adequate assurance has been more fortified when the results were compared to those reported from a Finite Element method study.     

An algebraic factorisation scheme for spectral element solution of incompressible flow and scalar transport

A spectral element algorithm for solution of the unsteady incompressible Navier–Stokes and scalar (species/heat) transport equations is developed using the algebraic factorisation scheme. The new algorithm utilises Nth order Gauss–Lobatto–Legendre points for velocity and the scalar, while (N-2)th order Gauss–Legendre points are used for pressure. As a result, the algorithm does not require inter-element continuity for pressure and pressure boundary conditions on solid surfaces. Implementations of the algorithm are performed for conforming and non-conforming grids. The latter is accomplished using both the point-wise matching and integral projection methods, and applied for grids with both polynomial and geometric non-conformities. Code validation cases include the unsteady scalar convection equation, and Kovasznay flow in two- and three-dimensional domains. Using cases with analytical solutions, the algorithm is shown to achieve spectral accuracy in space and second-order accuracy in time. The results for the Boussinesq approximation for buoyancy-driven flows, and the species mixing in a continuous flow micro-mixer are also included as examples of applications that require long-time integration of the scalar transport equations.     

A numerical method for accelerating the rate of convergence of the SIMPLE-like algorithm for flow through a thin filter

The control volume, finite difference method and the k-? tubulence model are employed in a numerical simulation of the turbulent fluid flow both outside and inside a blunt cylindrical sampler which houses a paper filter in its chamber. The presence of a paper filter, which has a very large resistance, results in a large pressure drop across the filter and this causes difficulties in making the SIMPLE or the SIMPLEC scheme converge. In order to improve the rate of convergence of the SIMPLE-like algorithm when the resistance of the filter is very large, an average pressure correction formula is proposed. Based on global mass conservation, a line average pressure correction for the paper filter is derived using a modified Darcy law for a porous medium. A combination of this formula and the SIMPLE-like algorithm can rapidly build up the pressure drop across the filter and hence dramatically improve the rate of convergence of the iterative scheme. Comparisons of the convergence histories and the numerical results for the fluid flow when using SIMPLE and SIMPLEC with the average pressure correction method show that the average pressure correction method for dealing with the paper filter significantly accelerates the rate of convergence of the iterative scheme.     

基于Boltzmann模型方程的气体运动论统一算法研究

模型方程出发,研究确立含流态控制参数可描述不同流域气体流动特征的气体分子速度分布函数方程; 研究发展气体运动论离散速度坐标法, 借助非定常时间分裂数值计算方法和NND差分格式, 结合DSMC方法关于分子运动与碰撞去耦技术, 发展直接求解速度分布函数的气体运动论耦合迭代数值格式; 研制可用于物理空间各点宏观流动取矩的离散速度数值积分方法, 由此提出一套能有效模拟稀薄流到连续流不同流域气体流动问题统一算法. 通过对不同Knudsen数下一维激波内流动、二维圆柱、三维球体绕流数值计算表明, 计算结果与有关实验数据及其它途径研究结果(如DSMC模拟值、N-S数值解)吻合较好, 证实气体运动论统一算法求解各流域气体流动问题的可行性. 尝试将统一算法进行HPF并行化程序设计, 基于对球体绕流及类``神舟'返回舱外形绕流问题进行HPF初步并行试算, 显示出统一算法具有很好的并行可扩展性, 可望建立起新型的能有效模拟各流域飞行器绕流HPF并行算法研究方向. 通过将气体运动论统一算法推广应用于微槽道流动计算研究, 已初步发展起可靠模拟二维短微槽道流动数值算法; 通过对Couette流、Poiseuille流、压力驱动的二维短槽道流数值模拟, 证实该算法对微槽道气体流动问题具有较强的模拟能力, 可望发展起基于Boltzmann模型方程能可靠模拟MEMS微流动问题气体运动论数值计算方法研究途径.      

七孔探针可压缩流场测量研究

介绍了七孔探针用于亚音速可压缩流的标定方法。作为一种可以同时获得流动速度大小、流动偏角、总压和静压的气动测量装置,七孔探针被广泛应用于各种流动测量,包括可压缩流动。但是它的校准过程周期很长,代价昂贵,影响了探针的推广。本文以数值计算为手段,对七孔探针进行亚音速可压缩流校准与测量的研究。结果表明,其校准拟合精度流动角为2%,内外区的总静压相对标准偏差都没有超过3%,高于相同状态下的实验校准精度。在实际应用中,本方法用于指导传统实验标定方法,可以节约大量的标定时间和成本,使七孔探针在亚音速可压缩流的测量变得简单可行。     

An interface-capturing method for incompressible two-phase flows. Validation and application to bubble dynamics

We report on the development and applications of an interface-capturing method aimed at computing three-dimensional incompressible two-phase flows involving high density and viscosity ratios, together with capillary effects. The numerical approach borrows some features to the Volume of Fluid method (since it is essentially based on the transport of the local volume fraction of the liquid) as well as to the Level Set technique (as no explicit reconstruction of the interface is carried out). The transport of the volume fraction is achieved by using a flux-limiting Zalesak scheme and the fronts are prevented from spreading in time by a specific strategy in which the velocity at nodes crossed by the interface is modified to keep the thickness of the transition region constant. As shown on several test cases, this algorithm allows the interface to deform properly while maintaining the numerical thickness of the transition region within three computational cells whatever the structure of the local flow field. The full set of governing equations is then used to investigate some fundamental aspects of bubble dynamics. More precisely we focus on the evolution of shape and rise velocity of a single bubble over a wide range of physical parameters and on head-on and side-by-side interactions between two rising bubbles.     

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.     

高阶谱元区域分解算法求解定常方腔驱动流

主要利用Jacobian-free的Newton-Krylov方法求解定常不可压缩Navier-Stokes方程,将基于高阶谱元法的区域分解Stokes算法的非定常时间推进步作为Newton迭代的预处理,回避了传统Newton方法Jacobian矩阵的显式装配,节省了程序内存,同时降低了Newton迭代线性系统的条件数,且没有非线性对流项的隐式求解,大大加快了收敛速度。对有分析解的Kovasznay流动的计算结果表明,本高阶谱元法在空间上有指数收敛的谱精度,且对定常解的Newton迭代是二次收敛的。本文模拟了二维方腔顶盖一致速度驱动流,同基准解符合得很好,表明本文方法是准确可靠的。本文还考虑了Re=800时方腔顶盖正弦速度驱动流,除得到已知的一个稳定对称解和一对稳定非对称解外,还获得了一对新的不稳定的非对称解。     

A Conceptual Study of Cavity Aeroacoustics Control Using Porous Media Inserts

In this study, an integrated flow simulation and aeroacoustics prediction methodology is applied to testing a sound control technique using porous inserts in an open cavity. Large eddy simulation (LES) combined with a three-dimensional Ffowcs Williams–Hawkings (FW–H) acoustic analogy is employed to predict the flow field, the acoustic sources and the sound radiation. The Darcy pressure – velocity law is applied to conceptually mimic the effect of porous media placed on the cavity floor and/or rear wall. Consequently, flow in the cavity could locally move in or out through these porous walls, depending on the local pressure differences. LES with “standard” subgrid-scale models for compressible flow is carried out to simulate the flow field covering the sound source and near fields, and the fully three-dimensional FW–H acoustic analogy is used to predict the sound field. The numerical results show that applying the conceptual porous media on cavity floor and/or rear wall could decrease the pressure fluctuations in the cavity and the sound pressure level in the far field. The amplitudes of the dominant oscillations (Rossiter modes) are suppressed and their frequencies are slightly modified. The dominant sound source is the transverse dipole term, which is significantly reduced due to the porous walls. As a result, the sound pressure in the far field is also suppressed. The preliminary study reveals that using porous-inserts is a promising technology for flow and sound radiation control.