A numerical study of turbulent line puffs via the renormalization group (RNG) k–ϵ model

The time evolution of a line puff, a turbulent non-buoyant element with significant momentum, is studied using the renormalization group (RNG) k–ϵ model. The numerical results show that the puff motion is characterized by a vortex pair flow; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless time of around 30. The added mass coefficient of the puff motion is found to be approximately unity. The predicted puff flow and mixing rate are substantially similar to those obtained from the standard k–ϵ model and are well supported by experimental data. The computed scalar field reveals significant secondary concentration peaks trailing behind in the wake of the puff. The present results suggest that the overall mixing rate of a puff is primarily determined by the large-scale motion and that streamline curvature probably plays a minor role. © 1998 John Wiley & Sons, Ltd.     

Calculation of turbulent fluid flow and heat transfer in ducts by a full Reynolds stress model

A computational method has been developed to predict the turbulent Reynolds stresses and turbulent heat fluxes in ducts by different turbulence models. The turbulent Reynolds stresses and other turbulent flow quantities are predicted with a full Reynolds stress model (RSM). The turbulent heat fluxes are modelled by a SED concept, the GGDH and the WET methods. Two wall functions are used, one for the velocity field and one for the temperature field. All the models are implemented for an arbitrary three‐dimensional channel. Fully developed condition is achieved by imposing cyclic boundary conditions in the main flow direction. The numerical approach is based on the finite volume technique with a non‐staggered grid arrangement. The pressure–velocity coupling is handled by using the SIMPLEC‐algorithm. The convective terms are treated by the van Leer scheme while the diffusive terms are handled by the central‐difference scheme. The hybrid scheme is used for solving the ε equation. The secondary flow generation using the RSM model is compared with a non‐linear k–ε model (non‐linear eddy viscosity model). The overall comparison between the models is presented in terms of the friction factor and Nusselt number. Copyright © 2003 John Wiley & Sons, Ltd.     

Turbulent mixing of accompanying and opposing flows

The solution of the problem of the edge of a turbulent jet, obtained in [1] starting from a model of the vortical motion of an ideal liquid, is generalized for the case of the turbulent mixing of two plane semibounded (accompanying or opposing) flows of an incompressible liquid. It is shown that the results of calculation are in qualitative agreement with experimental data and that they are close to them quantitatively. Some of the special characteristics of the method are discussed.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 92–101, November–December, 1973.     

Turbulent shear stress profiles in a bubbly channel flow assessed by particle tracking velocimetry

Particle tracking velocimetry (PTV) is applied to a bubbly two-phase turbulent flow in a horizontal channel at Re = 2 × 10⁴ to investigate the turbulent shear stress profile which had been altered by the presence of bubbles. Streamwise and vertical velocity components of liquid phase are obtained using a shallow focus imaging method under backlight photography. The size of bubbles injected through a porous plate in the channel ranged from 0.3 to 1.5 mm diameter, and the bubbles show a significant backward slip velocity relative to liquid flow. After bubbles and tracer particles are identified by binarizing the image, velocity of each phase and void fraction are profiled in a downstream region. The turbulent shear stress, which consists of three components in the bubbly two-phase flow, is computed by analysis of PTV data. The result shows that the fluctuation correlation between local void fraction and vertical liquid velocity provides a negative shear stress component which promotes frictional drag reduction in the bubbly two-phase layer. The paper also deals with the source of the negative shear stress considering bubble’s relative motion to liquid.     

BOUNDARY TREATMENT AND AN EFFICIENT PRESSURE ALGORITHM FOR INTERNAL TURBULENT FLOWS USING THE PDF METHOD

The generalized Langevin model, which is used to model the motion of stochastic particles in the velocity–composition joint probability density function (PDF) method for reacting turbulent flows, has been extended to incorporate solid wall effects. Anisotropy of Reynolds stresses in the near-wall region has been addressed. Numerical experiments have been performed to demonstrate that the forces in the near-wall region of a turbulent flow cause the stochastic particles approachi ng a solid wall to reverse their direction of motion normal to the wall and thereby, leave the near-wall layer. This new boundary treatment has subsequently been implemented in a full-scale problem to prove its validity. The test problem considered here is that of an isothermal, non-reacting turbulent flow in a two-dimensional channel with plug inflow and a fixed back-pressure. An efficient pressure correction method, developed in the spirit of the PISO algorithm, has been implemented. The pressure correction strategy is easy to implement and is completely consistent with the time- marching scheme used for the solution of the Lagrangian momentum equations. The results show remarkable agreement with both k–ϵ and algebraic Reynolds stress model calculations for the primary velocity. The secondary flow velocity and the turbulent moments are in better agreement with the algebraic Reynolds stress model predictions than the k– ϵ predictions. © 1997 by John Wiley & Sons, Ltd.     

Modelling a recirculating density-driven turbulent flow

A mathematical model of turbulent density-driven flows is presented and is solved numerically. A form of the k–? turbulence model is used to characterize the turbulent transport, and both this non-linear model and a sediment transport equation are coupled with the mean-flow fluid motion equations. A partitioned, Newton–Raphson-based solution scheme is used to effect a solution. The model is applied to the study of flow through a circular secondary sedimentation basin.     

Eulerian–Lagrangian 3-D simulations of unsteady two-phase gas–liquid flow in a rectangular column by considering bubble interactions

This work discusses the development of a three-dimensional Eulerian–Lagrangian CFD model for a gas–liquid flow in a rectangular column. The model resolves the time-dependent, three-dimensional motion of small gas bubbles in a liquid to simulate the dynamic characteristics of the oscillating bubble plume. Our model incorporates drag, gravity, buoyancy, lift, pressure gradient and virtual mass forces acting on a bubble rising in a liquid, and accounts for two-way momentum coupling between the phases. We use MUSIG model that provides a framework in which the population balance method together with the break up and coalescence models can be incorporated into three-dimensional CFD calculations. We use turbulent flow to describe liquid flow field. The standard κ–ε of turbulence is selected for calculating the properties of turbulent flow. The effect of aspect ratio of the column on the flow pattern, liquid velocity and gas hold-up profiles is discussed.     

Theory of the motion of rafts on a river and capsules in a pipeline. Prandtl's paradox

A method is proposed for calculating the motion of rafts floating on the surface of a river and for long capsules within delivery pipelines. The investigation gives a quantitative estimate of Prandtl's qualitative conclusions to the effect that floating streamlined bodies can travel faster than the flow of the river (for turbulent flow regime). The calculations showed, and an experiment confirmed that the motion of a raft or barge exceeds the flow velocity on the surface of a river by approximately 1.5 times. Similar conclusions are obtained for the motion of long capsules of neutral buoyancy in the case of turbulent flow of a liquid in a pipeline. In the case of a laminar flow regime of the liquid, effects of the opposite kind occur.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 88–93, November–December, 1980.     

Development of projection and artificial compressibility methodologies using the approximate factorization technique

Predictions for two-dimensional, steady, incompressible flows under both laminar and turbulent conditions are presented. The standard k-? turbulence model is used for the turbulent flows. 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. Mass conservation is enforced by either the pseudocompressibility method or the pressure correction method. Comparison of the two methods shows a superiority of the pressure correction method. Second- and fourth-order artifical dissipation terms are used in order to achieve good convergence and to handle the turbulence model equations efficiently. Several internal and external test cases are investigated, including attached and separated flows.     

Investigation of the flow in a circular cavity using stereo and tomographic particle image velocimetry

The turbulent flow over a circular cavity with an aspect ratio of D/H = 2 is investigated by multi-planar stereoscopic particle image velocimetry and with tomographic particle image velocimetry (PIV). The main aim of the study is the flow topology and the turbulent structure of the asymmetrical flow pattern that forms inside the cavity at these specific conditions. The flow field is measured in the vertical symmetry plane to describe the overall recirculation pattern in the cavity and the turbulent shear layer developing from the separation point. In this specific regime the shear layer fluctuations are recognized as those caused by instabilities together with the effect of the incoming boundary layer turbulence. Additional observations performed at several wall-parallel planes at different height inside the cavity allow to further evaluate the secondary flow circulation generated by this asymmetric regime. The observed flow pattern consists of a steady vortex, occupying the entire cavity volume and placed diagonally inside the cavity such to entrain the external flow from one side, capture it into a circulatory motion and eject it from the opposite side of the cavity. The spatial distribution of the turbulent fluctuations also reveals the same structure. The tomographic PIV measurement returns a visual inspection to the instantaneous three-dimensional structure of the turbulent fluctuations, which at the investigated height exhibit a low level of coherence with slightly elongated vortices in the recirculating flow inside the cavity.     

The size effect on the dispersion of a particle in a homogeneous isotropic turbulence

The mechanism of the response motion of a suspended particle to turbulent motion of its surrounding fluid is different according to size of turbulent eddies. The particle is dragged by the viscous force of large eddies, and meanwhile driven randomly by small eddies. Based on this understanding, the dispersion of a particle with finite size in a homogeneous isotropic turbulence is calculated in this study. Results show that there are two competing effects: when enhanced by the inertia of a particle, the long-term particle diffusivity is reduced by the finite size of the particle.     

Grid-averaged Lagrangian LES model for multiphase turbulent flow

A grid-averaged Lagrangian (GAL) model for dispersed particle motion in multiphase turbulent flow is presented to provide a large eddy simulation (LES) model for multiphase turbulent flow in which a quite large number of particles are involved. The GAL model is based on an averaging operation for a Lagrangian-type equation of motion of a particle over a computational grid volume and a procedure of reallocation of a dispersed particle cloud with its centroid movement to each grid. The model is therefore a mixed Eulerian–Lagrangian model which can effectively reduce computational time compared with existing Lagrangian-type models, without losing the advantage of Lagrangian-type models that they can properly describe the dynamical evolution of particles. Since the GAL model adopts the grid-volume averaging operation it can easily provide an effective SGS model for LES modeling of multiphase turbulent flow. The validity of the multiphase LES model developed, which is named the GAL-LES model, is confirmed through its application to a particle plume, in which the present model is found to simulate large-eddy motion usually observed in a jet and plume, and to give good agreements with experimental data.     

燃料消耗下充液航天器等效动力学建模与分析

在轨航天器贮腔内的液体可能表现出多种不同的运动模式, 主要包括液体相对于贮腔的整体性刚体运动、自由液面横向晃动、液体起旋后逐步发生明显的旋转晃动及液体自旋运动; 复合三自由度刚体摆晃动模型能够较为全面地描述这些液体运动模式, 同时为研究起旋阶段的液体晃动动力学问题提供了有效手段. 本文对非线性液体晃动刚体摆复合模型作进一步发展, 考虑模型等效参数随贮腔充液比的变化, 提出了变参数的刚体摆复合模型, 该模型适用于研究燃料消耗下非线性晃动类充液航天器大范围运动耦合动力学问题. 采用刚体摆复合模型对球形贮腔内的液体晃动进行等效后, 基于混合坐标意义下的拉格朗日方程推导了一类充液航天器轨道-姿态-晃动全耦合的动力学方程组, 并展开了充液航天器大角度三轴稳定姿态机动和零冲量轨道机动仿真以及航天器耦合动力学响应特性分析. 研究表明: 液体相对于贮腔的运动会造成航天器主刚体位置发生偏移, 当航天器在执行零冲量机动时, 燃料消耗会造成航天器的轨道平动速度无法收敛到零; 贮腔偏心布放时, 航天器在执行轨道机动过程中贮腔内液体易发生剧烈而且形式复杂的晃动行为, 进而可能造成航天器刚体运动的不稳定.      

Analysis and Modeling of the Turbulent Diffusion of Turbulent Kinetic Energy in Natural Convection

Buoyant flows often contain regions with unstable and stable thermal stratification from which counter gradient turbulent fluxes are resulting, e.g. fluxes of heat or of any turbulence quantity. Basing on investigations in meteorology an improvement in the standard gradient-diffusion model for turbulent diffusion of turbulent kinetic energy is discussed. The two closure terms of the turbulent diffusion, the velocity-fluctuation triple correlation and the velocity-pressure fluctuation correlation, are investigated based on Direct Numerical Simulation (DNS) data for an internally heated fluid layer and for Rayleigh–Bénard convection. As a result it is decided to extend the standard gradient-diffusion model for the turbulent energy diffusion by modeling its closure terms separately. Coupling of two models leads to an extended RANS model for the turbulent energy diffusion. The involved closure term, the turbulent diffusion of heat flux, is studied based on its transport equation. This results in a buoyancy-extended version of the Daly and Harlow model. The models for all closure terms and for the turbulent energy diffusion are validated with the help of DNS data for internally heated fluid layers with Prandtl number Pr = 7 and for Rayleigh–Bénard convection with Pr = 0.71. It is found that the buoyancy-extended diffusion model which involves also a transport equation for the variance of the vertical velocity fluctuation gives improved turbulent energy diffusion data for the combined case with local stable and unstable stratification and that it allows for the required counter gradient energy flux.     

Numerical prediction of turbulent thermocapillary convection in superposed fluid layers with a free interface

The present paper introduces a new numerical method for predicting the characteristics of thermocapillary turbulent convection in a differentially-heated rectangular cavity with two superposed and immiscible fluid layers. The unsteady Reynolds form of the Navier–Stokes equations and energy equation are solved by using the control volume approach on a staggered grid system using SIMPLE algorithm. The turbulence quantities are predicted by applying the standard k–ε turbulence model. The level set formulation is applied for predicting the topological changes of the interface separating the two fluid layers and to provide an accurate and robust modeling of the interfacial normal and tangential stresses. The computational results obtained showed good agreement when compared with the previous experimental, numerical and analytical benchmark data for different validation cases in both laminar and turbulent regimes. The present numerical method is then applied to predict the velocity and temperature distribution in two immiscible liquid layers with undeformable interface for a wide range of Marangoni numbers. The laminar-turbulent transition is demonstrated by obtaining the turbulence features at high interfacial temperature gradient which is characterized by high Marangoni number. The effect of increasing Marangoni number on the interface dynamics in turbulent regime is also investigated.     

EQUIVALENT DYNAMICS MODELING AND ANALYSIS OF LIQUID-FILLED SPACECRAFT WITH FUEL CONSUMPTION 1)

在轨航天器贮腔内的液体可能表现出多种不同的运动模式, 主要包括液体相对于贮腔的整体性刚体运动、自由液面横向晃动、液体起旋后逐步发生明显的旋转晃动及液体自旋运动; 复合三自由度刚体摆晃动模型能够较为全面地描述这些液体运动模式, 同时为研究起旋阶段的液体晃动动力学问题提供了有效手段. 本文对非线性液体晃动刚体摆复合模型作进一步发展, 考虑模型等效参数随贮腔充液比的变化, 提出了变参数的刚体摆复合模型, 该模型适用于研究燃料消耗下非线性晃动类充液航天器大范围运动耦合动力学问题. 采用刚体摆复合模型对球形贮腔内的液体晃动进行等效后, 基于混合坐标意义下的拉格朗日方程推导了一类充液航天器轨道-姿态-晃动全耦合的动力学方程组, 并展开了充液航天器大角度三轴稳定姿态机动和零冲量轨道机动仿真以及航天器耦合动力学响应特性分析. 研究表明: 液体相对于贮腔的运动会造成航天器主刚体位置发生偏移, 当航天器在执行零冲量机动时, 燃料消耗会造成航天器的轨道平动速度无法收敛到零; 贮腔偏心布放时, 航天器在执行轨道机动过程中贮腔内液体易发生剧烈而且形式复杂的晃动行为, 进而可能造成航天器刚体运动的不稳定.     

Diffusion of suspended particles in an isotropic turbulence field

A model for turbulent motion is proposed which makes it possible to evaluate the pulsation characteristics and the diffusion coefficients of the dispersed phase and also makes it possible to describe the effect of the suspended particles on the turbulence of the dispersing medium. Specific calculations are made for the situation when the undisturbed turbulent field is isotropic.The diffusion of an admisture having inertia in a turbulent stream has been studied previously on the assumption that the three-dimensional turbulence characteristics have practically no effect on the behavior of the suspended particles, so that the random motion of the latter is described by ordinary differential equations containing the natural independent variable the motion travel time [1–4]. In many cases this assumption is incorrect and the corresponding theory is obviously deficient. For example, a fundamental result of this theory, asserting that the turbulent diffusion coefficients of the particles and of the fluid moles are equal for a long diffusion time, is obviously incorrect if the relative motion of the particles is significant [5].     

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.     

K-ε-T model of dense liquid-solid two-phase turbulent flow and its application to the pipe flow

To predict the characteristics of dense liquid-solid two-phase flow, K-ε-T model is established, in which the turbulent flow of fluid phase is described with fluid turbulent kinetic energy K_f and its dissipation rate ε_f, and the particles random motion is described with particle turbulent energy K_p and its dissipation rate ε_p and pseudothermal temperature T_p. The governing equations are also derived. With K-ε-T model, numerical study of dense liquid-solid two-phase turbulent up-flow in a pipe is performed. The calculated results are in good agreement with experimental data of Alajbegovic et al. (1994), and some flow features are captured.     

Dynamics of large-scale eddies in turbulent flows

A model equation based on the equipartition of the turbulent dissipation is proposed for describing the dynamics of large-scale eddies in turbulent flows. The equation is reducible to the equation of motion of an inviscid fluid, so that the motion of the large-scale eddies can be described in terms of inviscid fluid dynamics. It is found that the large-scale eddies are always weakened by the background turbulence and their evolution is slowed down compared with the corresponding inviscid motion. In the case of turbulent mixing layer, its linear growth in downstream direction is accounted for by the exponential growth in time of the perturbation in an inviscid plane vortex sheet.