用二阶矩模型研究Richtmyer-Meshkov不稳定性诱导湍流混合

发展了考虑密度脉动和各向异性湍流的二阶矩模型,强调了涉及湍流能量产生项的关联。采用该模型对Poggi等的激波管实验进行了模拟。通过与实验结果的比较分析,验证了采用的模型封闭、模型常数、数值算法和程序实现是合适的。在此基础上,进一步探讨了冲击马赫数和Atwood数对混合的影响。     

EFFECT OF EMPIRICAL COEFFICIENTS ON SIMULATION IN TWO-SCALE SECOND-ORDER MOMENT PARTICLE-PHASE TURBULENCE MODEL

A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concept of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision. The proposed model is used to simulate gas-particle downer reactor flows. The computational results of both particle volume fraction and mean velocity are in agreement with the experimental results. After analyzing effects of empirical coefficient on prediction results, we can come to a conclusion that, inside the limit range of empirical coefficient, the predictions do not reveal a large sensitivity to the empirical coefficient in the downer reactor, but a relatively great change of the constants has important effect on the prediction.     

SIMULATION OF NOx FORMATION IN TURBULENT SWIRLINGCOMBUSTION USING A USM TURBULENCE—CHEMISTRY MODEL

A unified second-order moment (USM) turbulence-chemistry model for simulating Nox formation in turbulent combustion is proposed. All of correlations, including the correlation of the reaction-rate coefficient fluctuation with the concentration fluctuation, are closed by the transport equations in the same form. This model discards the approximation of series expansion of the exponential function or the approximation of using the product of several 1-D PDF‘s instead of a joint PDF. It is much simpler than other refined models, such as the PDF transport equation model and the condi-tional moment closure model. The proposed model is used to simulate methane-air swirling turbulent combustion and Nox formation. The prediction results are in good agreement with the experimental results.     

环形通道内湍流旋流流动的数值模拟

本文应用一种考虑湍流－旋流相互作用及湍流脉动各向异性的新的代数Ｒｅｙｎｏｌｄｓ应力模型,对环形通道内的湍流旋流流动进行了数值模拟,研究了改主为旋流流数,进口轴向速度及半径比等参数对环形通道内湍流流动的影响,以及对强化环形通道内传热的作用。     

Numerical simulation of turbulent impinging jet on a rotating disk

The calculations of quasi‐three‐dimensional momentum equations were carried out to study the influence of wall rotation on the characteristics of an impinging jet. The pressure coefficient, the mean velocity distributions and the components of Reynolds stress are calculated. The flow is assumed to be steady, incompressible and turbulent. The finite volume scheme is used to solve the continuity equation, momentum equations and k–ε model equations. The flow characteristics were studied by varying rotation speed ω for 0?ω?167.6 rad/s, the distance from nozzle to disk (H/d) was (3, 5, 8 and 10) and the Reynolds number Re base on V_J and d was 1.45 × 10⁴. The results showed that, the radial velocity and turbulence intensity increase by increasing the rotation speed and decrease in the impingement zone as nozzle to disk spacing increases. When the centrifugal force increases, the radial normal stresses and shear stresses increase. The location of maximum radial velocity decreases as the local velocity ratio (α) increases. The pressure coefficient depends on the centrifugal force and it decreases as the distance from nozzle to plate increases. In impingement zone and radial wall jet, the spread of flow increases as the angular velocity decreases The numerical results give good agreement with the experiment data of Minagawa and Obi (Int. J. of Heat and Fluid Flow 2004; 25 :759–766). Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical simulation of axisymmetric turbulent jets

The flow in axisymmetric turbulent jets is numerically simulated with the use of a semi-empirical second-order turbulence model including differential transport equations for the normal Reynolds stresses. Calculated results are demonstrated to agree with experimental data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 55–60, September–October, 2008.     

Direct numerical simulation of particle dispersion in a turbulent jet considering inter-particle collisions

To investigate the behaviour of inter-particle collision and its effects on particle dispersion, direct numerical simulation of a three-dimensional two-phase turbulent jet was conducted. The finite volume method and the fractional-step projection algorithm were used to solve the governing equations of the gas phase fluid and the Lagrangian method was applied to trace the particles. The deterministic hard-sphere model was used to describe the inter-particle collision. In order to allow an analysis of inter-particle collisions independent of the effect of particles on the flow, two-way coupling was neglected. The inter-particle collision occurs frequently in the local regions with higher particle concentration of the flow field. Under the influence of the local accumulation and the turbulent transport effects, the variation of the average inter-particle collision number with the Stokes number takes on a complex non-linear relationship. The particle distribution is more uniform as a result of inter-particle collisions, and the lateral and the spanwise dispersion of the particles considering inter-particle collision also increase. Furthermore, for the case of particles with the Rosin–Rammler distribution (the medial particle size is set d₅₀ = 36.7 μm), the collision number is significantly larger than that of the particles at the Stokes number of 10, and their effects on calculated results are also more significant.     

SECOND-ORDER MOMENT MODEL FOR DENSE TWO-PHASE TURBULENT FLOW OF BINGHAM FLUID WITH PARTICLES

The USM-θmodel of Bingham fluid for dense two-phase turbulent flow was developed, which combines the second-order moment model for two-phase turbulence with the particle kinetic theory for the inter-particle collision. In this model, phases interaction and the extra term of Bingham fluid yield stress are taken into account. An algorithm for USM-θmodel in dense two-phase flow was proposed, in which the influence of particle volume fraction is accounted for. This model was used to simulate turbulent flow of Bingham fluid single-phase and dense liquid-particle two-phase in pipe. It is shown USM-θmodel has better prediction result than the five-equation model, in which the particle-particle collision is modeled by the particle kinetic theory, while the turbulence of both phase is simulated by the two-equation turbulence model. The USM-θmodel was then used to simulate the dense two-phase turbulent up flow of Bingham fluid with particles. With the increasing of the yield stress, the velocities of Bingham and particle decrease near the pipe centre. Comparing the two-phase flow of Bingham-particle with that of liquid-particle, it is found the source term of yield stress has significant effect on flow.     

管内均相湍流燃烧加速的数值模拟

通过均相流体模型、湍流k 模型和EBU(EddyBreak Up) Arrhenius燃烧模型 ,选用Simple格式对管中戊烷和空气的燃烧实例进行了数值求解。其结果反映了燃烧导致的爆炸过程中管内流场各参数的变化规律 ,揭示了管内燃烧、流动、湍流之间的正反馈耦合关系 ,并与实验结果和相关数值结果基本相符。     

Interaction of a gas-droplet turbulent jet with a cocurrent high-velocity high-temperature gas flow

A mathematical model and a method for calculating a gas-droplet turbulent jet with allowance for velocity nonequilibrium and virtual mass of the condensed phase during turbulent fluctuations and also heat and mass transfer within the three-temperature scheme are developed. Methodical calculations are performed. The results of these calculations are in reasonable agreement with available experimental data. The structure of the gas-droplet jet in a cocurrent high-velocity high-temperature gas flow is studied by numerical methods. The ratio of intensities of heat and mass transfer between the phases and turbulent diffusion transfers of substances is found to be different at the initial, transitional, and basic segments of the jet. This difference is responsible for the nonmonotonic axial distribution of vapor density and the lines of the halved mass flow of the condensed phase. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 85–94, May–June, 2008.     

Schlieren measurement of turbulent structure in a diffusion flame

The regular and random mixing structures in a turbulent diffusion flame were investigated using the quantitative, dynamic crossed-beam schlieren method. Evidence was found close to the nozzle relating to the vortexlike structure of eddies surrounding the central fuel jet flow. The observations also make possible resolution of turbulent intensity, scales, convection, and spectra within the diffusion flame without the use of seeding or intrusion of measuring probes. It is found that length scales and other turbulence parameters in the diffusion flame progressively revert to values similar to those expected and observed in scalar passive mixing as the combustion reaction intensity reduces with axial distance from the nozzle system.     

Numerical simulation of turbulent flows in trapezoidal meandering compound open channels

Three‐dimensional (3D) numerical study is presented to investigate the turbulent flow in meandering compound open channels with trapezoidal cross‐sections. The flow simulation is carried out by solving the 3D Reynolds‐averaged continuity and Navier–Stokes equations with Reynolds stress equation model (RSM) for steady‐state flow. Finite volume method (FVM) is applied to numerically solve the governing equations of fluid flow. The velocity magnitude, tangential velocity, transverse velocity and Reynolds stresses are calculated for various flow conditions. Good agreement between the simulated and available laboratory measurements was obtained, indicating that the RSM can accurately predict the complicated flow phenomenon. Comparison of the calculated secondary currents of four cases (one being inbank flow and other three being overbank flow) with different water depths reveals that (i) the inbank flow exhibits different flow behaviors from that of the overbank flow does and (ii) the water depth has significant effects on the magnitude and direction of secondary currents. Copyright © 2010 John Wiley & Sons, Ltd.     

湍流加速火焰的三维数值模拟

火焰在设有障碍物的管内传播时会自身加速，并可能导致爆炸。本文基于湍流κ-ε模型和改进的EBU—Arrhenius反应模型，对该现象进行了三维空间的数值模拟。计算结果反映了障碍物、湍流和火焰之间相互作用的正反馈机理，描绘了火焰在管内加速传播的三维图像。     

Generation of internal waves by a turbulent jet in a stratified fluid

The process of generation of internal waves by an initially cylindrical, turbulent jet with a Gaussian profile of the average horizontal velocity component in a fluid with stable linear density stratification is investigated by direct numerical simulation. It is shown that on time intervals Nt < 30, where N is the buoyancy frequency, the vertical velocity pulsations collapse, which is accompanied by the generation of internal waves whose spatial period is close to the wavelength of the spiral mode of jet instability in a homogeneous fluid. The wave dynamics and kinematics can be satisfactorily described by the linear theory for a pulsed source and their parameters are in good agreement with the parameters of the “coherent” internal waves generated by a stratified wake in a laboratory experiment. At large times the wave generation ceases and the variations of the fluid density are localized in the neighborhood of the centers of large-scale vortices formed in the horizontal plane in the neighborhood of the jet.     

NUMERICAL SIMULATION OF THREE DIMENSIONAL TURBULENT FLOW IN SUDDENLY EXPANDED RECTANGULAR DUCT

ＮＵＭＥＲＩＣＡＬＳＩＭＵＬＡＴＩＯＮＯＦＴＨＲＥＥＤＩＭＥＮＳＩＯＮＡＬＴＵＲＢＵＬＥＮＴＦＬＯＷＩＮＳＵＤＤＥＮＬＹＥＸＰＡＮＤＥＤＲＥＣＴＡＮＧＵＬＡＲＤＵＣＴＮＵＭＥＲＩＣＡＬＳＩＭＵＬＡＴＩＯＮＯＦＴＨＲＥＥＤＩＭＥＮＳＩＯＮＡＬＴＵＲＢＵ...     

Large eddy simulation of turbulent channel flow using an algebraic model

In this paper an algebraic model from the constitutive equations of the subgrid stresses has been developed. This model has an additional term in comparison with the mixed model, which represents the backscatter of energy explicitly. The proposed model thus provides independent modelling of the different energy transfer mechanisms, thereby capturing the effect of subgrid scales more accurately. The model is also found to depict the flow anisotropy better than the linear and mixed models. The energy transfer capability of the model is analysed for the isotropic decay and the forced isotropic turbulence. The turbulent plane channel flow simulation is performed over three Reynolds numbers, Re_τ=180, 395 and 590, and the results are compared with that of the dynamic model, Smagorinsky model, and the DNS data. Both the algebraic and dynamic models are in good agreement with the DNS data for the mean flow quantities. However, the algebraic model is found to be more accurate for the turbulence intensities and the higher‐order statistics. The capability of the algebraic model to represent backscatter is also demonstrated. Copyright © 2005 John Wiley & Sons, Ltd.     

Large eddy simulation of turbulent premixed and stratified combustion using flame surface density model coupled with tabulation method

Large eddy simulations (LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density (FSD) model incorporated with two different wrinkling factor models, i.e., the Muppala and Charlette2 wrinkling factor models, is used to describe combustion/turbulence interaction, and the flamelet generated manifolds (FGM) method is employed to determine major scalars. This coupled sub-grid scale (SGS) combustion model is named as the FSD-FGM model. The FGM method can provide the detailed species in the flame which cannot be obtained from the origin FSD model. The LES results show that the FSD-FGM model has the ability of describing flame propagation, especially for stratified flames. The Charlette2 wrinkling factor model performs better than the Muppala wrinkling factor model in predicting the flame surface area change by the turbulence. The combustion characteristics are analyzed in detail by the flame index and probability distributions of the equivalence ratio and the orientation angle, which confirms that for the investigated stratified flame, the dominant combustion modes in the upstream and downstream regions are the premixed mode and the back-supported mode, respectively.     

Direct numerical simulation of turbulent pipe flow

Fully developed turbulent pipe flow at low Re-number is studied by means of direct numerical simulation (DNS). In contrast to many previous DNS's of turbulent flows in rectangular geometries, the present DNS code, developed for a cylindrical geometry, is based on the finite volume technique rather than being based on a spectral method. The statistical results are compared with experimental data obtained with two different experimental techniques. The agreement between numerical and experimental results is found to be good which indicates that the present DNS code is suited for this kind of numerical simulations.     

A prediction of the pressure field in a plane turbulent jet using an algebraic stress model

A numerical prediction is obtained for the mean pressure field in the similarity region of a plane turbulent jet. An algebraic stress model, which introduces non-isotropic relations for the Reynolds stress components, is used to close the mean momentum equation. The full two-dimensional form of the transport equations is retained and the resultant equation set solved elliptically. The numerical prediction simulates many of the characteristics of the pressure field measured by experimental studies. However, the overall level of the predicted field is lower than the experimental values. The level obtained for the mean pressure field depends strongly on the prediction for the transverse normal Reynolds stress component 〈u₂u₂〉. The pressure field is shown to represent a small negative contribution to the net strearnwise momentum balance.