Decaying turbulence in a passive mean scalar gradient

We consider the generation of passive scalar fluctuations by decaying isotropic turbulence in the presence of a uniform mean scalar gradient. At high Reynolds numbers, two distinct similarity states may be established depending on the form of the energy spectrum at low wavenumber magnitude (k). In the first similarity state characterized by a low wavenumber magnitude energy spectrum proportional tok ², the mean-square scalar fluctuation grows liket ^4/5, while in the second similarity state characterized by a spectrum proportional tok ⁴, the mean-square scalar fluctuation grows approximately liket ^4/7. These two high Reynolds number asymptotic similarity states have been subsequently confirmed by large-eddy numerical simulations. As a consequence of the decreasing flow Reynolds number as the turbulence decays, these similarity states do not continue indefinitely. At very long times, a final period of decay of the turbulence occurs, and in this final period, the mean-square scalar fluctuation in the first state continues to grow liket ^1/2, while that in the second state ultimately decays liket ^–1/2.     

A RANS simulation toward the effect of turbulence and cavitation on spray propagation and combustion characteristics

A multidimensional computational fluid dynamic code was developed and integrated with probability density function combustion model to give the detailed account of multiphase fluid flow. The vapor phase within injector domain is treated with Reynolds-averaged Navier–Stokes technique. A new parameter is proposed which is an index of plane-cut spray propagation and takes into account two parameters of spray penetration length and cone angle at the same time. It was found that spray propagation factor (SPI) tends to increase at lower r/d ratios, although the spray penetration tends to decrease. The results of SPI obtained by empirical correlation of Hay and Jones were compared with the simulation computation as a function of respective r/d ratio. Based on the results of this study, the spray distribution on plane area has proportional correlation with heat release amount, NO _x emission mass fraction, and soot concentration reduction. Higher cavitation is attributed to the sharp edge of nozzle entrance, yielding better liquid jet disintegration and smaller spray droplet that reduces soot mass fraction of late combustion process. In order to have better insight of cavitation phenomenon, turbulence magnitude in nozzle and combustion chamber was acquired and depicted along with spray velocity.     

High Reynolds number gravity currents along V-shaped valleys

The motion of saline gravity currents propagating horizontally in a tank of rectangular upper cross section and lower V-shaped valley is investigated both by lock-exchange experiments and a box model. The experiments were performed for equal depths of heavy and light fluid on both sides of the lock gate. The density ratio of the heavy fluid to the light fluid was in the range 1.04–1.13 and the lock height to length aspect ratios ranged from 0.5 to 1.6. We show that a box model with the Froude number of the head defined using the distance from the top of the current to the bottom of the valley predicts the position of the head in close agreement with the experiments. The presence of the valley results in three major differences in the gravity current compared to that flowing along a flat bottom. These are (a) the front of the current is approximately parabolic with radius of curvature proportional to the initial depth of the current, (b) for sufficiently large time t, the velocity of the current in the V-shaped valley varies as t^−1/5 compared to t^−1/3 in the flat bottom case, and (c) the width of the current in the V-shaped valley decreases with time t according to t^−2/5. Based on the box model, we predict that the steeper the flanks of the valley the faster the flow.     

Comparative study of Euler/Euler and Euler/Lagrange approaches simulating evaporation in a turbulent gas–liquid flow

This study deals with the Reynolds‐averaged Navier–Stokes simulation of evaporation in a turbulent gas–liquid flow in a three‐dimensional duct, focussing on the results obtained by a four‐equation turbulence model within the framework of the Euler/Euler approach for multiphase flow calculations: in addition to the two‐equation k?ε model describing the turbulence of the continuous (C) phase, the computational model employs transport equations for the turbulence kinetic energy of the disperse (D) phase and for the velocity covariance q=〈{u}^D{u}^C〉^D. In the present study, the evaporation model according to Abramzon and Sirignano (Int. J. Heat Mass Transfer 1989; 32 :1605–1618) has been extended by introducing an additional transport equation for a newly defined quantity ā, defined as the phase‐interface surface fraction. This allows the change in the drop diameter to be quantified in terms of a probability density function. The source term in the equation describing the dynamics of the volumetric fraction of the dispersed phase α^D is related to the evaporation time scale τ_Γ. The performance of the new model is evaluated by performing a comparative analysis of the results obtained by simulating a polydispersed spray in a three‐dimensional duct configuration with the results of the Euler/Lagrange calculations performed in parallel. Prior to these calculations, some selected (solid) particle‐laden flow configurations were computationally examined with respect to the validation of the background, four‐equation, eddy‐viscosity‐based turbulence model. Copyright © 2008 John Wiley & Sons, Ltd.     

On the Stress Tensor Jump at the Interface between a Polyatomic Gas and a Condensed Phase

For the components of the hydrodynamic stress tensor the boundary condition at the interface between a polyatomic gas and a condensed phase is obtained. The boundary-value problem is solved within the framework of the previously proposed kinetic model by the method of semi-spatial moments, taking into account the rotational degrees of freedom of the gas molecules. The gas-kinetic coefficient C _p entering into the boundary condition for the stress tensor components depends on the accommodation coefficient of the tangential momentum, q, the accommodation coefficients of the translational, _t , and rotational, _r , energy components, and the Prandtl number. This coefficient is calculated for several polyatomic gases.     

Shock waves and turbulence in a hypersonic inlet

A numerical investigation for an axisymmetric hypersonic turbulent inlet flow field of a perfect gas is presented for a three-shock configuration consisting of a biconic and a cowl. An upwind parabolized Navier-Stokes solver based on Roe's scheme is used to compute an oncoming flow Mach numberM =8, temperatureT =216 K, and pressureP =5.5293×10³ N/m². In order to assess the flow quantities, the interaction between shock and turbulence, and the inlet efficiency, three different flow calculations — laminar, turbulent with incompressible and compressible two-equationk- turbulence models — have been performed in this work.Computational results show that turbulence is markedly enhanced across an oblique shock with step-like increases in turbulence kinetic energy and dissipation rate. This enhancement is at the expense of the mean kinetic energy of the flow. Therefore, the velocity behind the shock is smaller in turbulent flow and hence the shock becomes stronger. The entropy increase through a shock is caused not only by the amplification of random molecular motion, but also by the enhancement of the chaotic turbulent flow motion. However, only the compressiblek- turbulence model can properly predict a decrease in turbulence length scale across a shock. Our numerical simulation reveals that the incompressiblek- turbulence model exaggerates the interaction between shock and turbulence with turbulence kinetic energy and dissipation rate remaining high and almost undissipated far beyond the shock region. It is shown that proper modeling of turbulence is essential for a realistic prediction of hypersonic inlet flowfield. The performed study shows that the viscous effect is not restricted in the boundary layer but extends into the main flow behind a shock wave. The loss of the available energy in the inlet performance therefore needs to be determined from the shock-turbulence interaction. The present study predicts that the inlet efficiency becomes relatively lower when turbulence is taken into account.     

One dimensional lattice gases with rapidly decreasing interaction

We consider the pressure and the correlation functions of a one dimensional lattice gas in which the mutual interaction decreases as r exp-n ^t, (r, t>0), when the interparticle distance n. We prove that such a system cannot show phase transitions of order k1 in the sense that the pressure and the correlation functions are infinitely differentiable with respect to any relevant parameter (such as the temperature or the chemical potential).     

Transient positively and negatively buoyant turbulent round jets

 The transient character of the jet issuing from an upward nozzle centered at the bottom of a vertical cylindrical tank into bulk liquid of a different density was measured using flow visualization and PIV for varying densimetric Froude numbers by varying the jet Reynolds numbers and the ratios of fluid densities. Positively buoyant jets penetrate to the free surface, driven by both momentum and buoyancy in the upward direction. The lighter jet fluid stratifies in a layer above the bulk liquid. Upon starting, a negatively buoyant jet has three stages. First the jet penetrates to its maximum height in the tank. Then the jet penetration decreases due to the downward backflow of heavier fluid surrounding the jet, which reduces the jet’s upward momentum. Finally the jet penetration height fluctuates around a mean value about 70% the maximum height of penetration. For small negative Froude numbers, the flow is fountain-like. The downward flow turns radially outward as it reaches the bottom of the tank and eventually an annular recirculation zone forms at the bottom of the tank with vortical motion opposite the vorticity of the jet. For large negative Froude numbers, the spreading of the jet extends far enough so the annular downward flow is along the walls of the tank resulting in a large annular recirculation zone. The penetration depth, h, and time, t, scale with buoyancy flux, F, and the jet momentum flux, M, as hM ^-3/4∣F∣^1/2 and t∣F∣/M to collapse the transient jet penetration height data onto a single curve over a wide range of Froude numbers for either positively or negatively buoyant jets. Received: 8 June 1998/Accepted: 3 February 1999     

Mixing with first-order decay in variable-velocity porous media flow

The transport and mixing of solutes undergoing first-order decay is central to many problems in groundwater hydrology. Mixing in porous media flow occurs due to advective dilution, hydrodynamic dispersion, and molecular diffusion. Mixing is stronger in regions of higher velocity, and weaker in slower-moving regions. Two-dimensional numerical experiments show that concentration profiles normal to the flow direction are displaced toward regions of slower flow in a flow field with a velocity gradient. Variable-velocity flow fields occur in subsurface flow around permeability heterogeneities, between recirculation cells, and in flow driven by natural convection. We examine two-dimensional solute concentration fields rather than breakthrough curves since for many complicated flow patterns, the breakthrough curve cannot discern important details of the concentration field. For the case of a species undergoing first-order decay, the effect of parent accumulation in regions of low velocity is enhanced for the daughter species because (i) the rate of daughter production is proportional to the local concentration of parent, and (ii) mixing is proportional to the local velocity. The resulting displacement of concentration profiles toward low-velocity regions may have important consequences for subsurface radionuclide transport and also for flows in chemically reactive systems and strongly coupled systems.Nomenclature d daughter component - d, ^(k) molecular diffusivity m² s^-1 - D dispersion coefficient m² s^-1 - dispersion tensor m²s^-1 - g acceleration of gravity vector m s^-2 - F Darcy flux vector kg m² s^-1 - identity matrix - k permeability m² - Kd distribution coefficient m³ kg^-1 - M mass accumulation term kg m^-3 - MW molecular weight kg mol^-1 - n outward unit normal vector - NK number of mass components (species) - p parent component - P total pressure Pa - q source term kg m^-3 s^-1 - R retardation factor - t time s, years - t _1/2 half-life s, years - u Darcy velocity in Y-direction m s^-1 - u Darcy velocity vector m s^-1 - v pore velocity m s^-1 - V volume m³ - X mass fraction - Y Y-coordinate - Z Z-coordinate (positive upward) - intrinsic dispersivity m - x distance interval m - surface area m² - decay constant s^-1 - dynamic viscosity kg m^-1 s^-1 - Ø porosity - fluid density kg m^-3 - tortuosity factor - d dispersion - L longitudinal - 0 reference value - T transverse - k mass component index     

Diffusion of moisture in deformable porous media. Anisotropic fibrous materials

This paper presents a study on the deformation of anisotropic fibrous porous media subjected to moistening by water in the liquid phase. The deformation of the medium is studied by applying the concept of effective stress. Given the structure of the medium, the displacement of the solid matrix is not taken into account with respect to the displacement of the liquid phase. The transport equations are derived from the model proposed by Narasimhan. The transport coefficients and the relation between the variation in apparent density and effective stress are obtained by test measurements. A numerical model has been established and applied for studying drip moistening of mineral wool samples capable or incapable of deformation.Nomenclature D mass diffusion coefficient [L²t^–1] - e void fraction - g gravity acceleration [Lt^–2] - J mass transfer density [ML^–2t^–1] - K hydraulic conductivity [Lt^–1] - K _s hydraulic conductivity of the solid phase [Lt^–1] - K ^* hydraulic conductivity of the deformable porous medium [Lt^–1] - P pressure of moistening liquid [ML^–1 t^–2] - S degree of saturation - t time [t] - V speed [Lt^–1] - X horizontal coordinate [L] - Z vertical coordinate measured from the bottom of porous medium [L] - z z-coordinate [L] Greek Letters porosity - ₁ total hydric potential [L] - _g gas density [ML^–3] - ₁ liquid density [ML^–3] - ₀ apparent density [ML^–3] - _s density of the solid phase [ML^–3] - density of the moist porous medium [ML^–3] - external load [ML^–1t^–2] - effective stress [ML^–1t^–2] - bishop's parameter - matrix potential or capillary suction [L] Indices g gas - 1 moistening liquid - p direction perpendicular to fiber planes - s solid matrix - t direction parallel to fiber planes - v pore Exponent * movement of solid particles taken into account     

The effect of leaching on mixing in porous media with a single heterogeneity

This paper presents an exploratory study of the effect of leaching on mixing in a porous medium containing a single heterogeneity to investigate the effect of the heterogeneity and time-dependent pore structure on dispersion. A percolation-convection simulation (PCS) model is used along with laboratory model experiments to study the mixing. The results show that mixing changes when the pores of the models are leached and that there is a change in regime influence during leaching. The simulation represents the mixing through a first leach for homogeneous media and for heterogeneous media with significant changes in permeability. If the pore structure is changing with time, prediction of mixing must include effects of heterogeneity and regime influence. Although the experimental results are representative of idealized laboratory sized systems they provide insight into the effects of leaching in heterogeneous media. Further the simulation may be useful on a field scale.Nomenclature b molecular weight, gm/mol - C concentration, mol/cm³ - C ₀ initial concentration, mol/cm³ - d _rms root-mean-squared distance, cm - d ₅₀ 50% grain size, cm - D channel depth, cm - f _n fraction of input tracer in effluent at time t _n - K ₁ permeability of flow field outside of heterogeneity, cm² - k ₂ permeability of heterogeneity, cm² - k _S reaction rate constant, cm/min - K _L microscopic dispersion coefficient in the longitudinal direction, cm²/sec - K _O overall dispersion coefficient in the longitudinal direction, cm²/sec - K _T microscopic dispersion coefficient in the transverse direction, cm²/sec - L length of channel, cm - n exponent for velocity - P pressure, N/M² - Pe Peclet number, Lv/K _O - P _ext local pressure outside heterogeneity, N/M² - P _int local pressure inside heterogeneity, N/M² - Q volumetric flow rate, cm³/sec - R channel half width, cm - t time, sec - W _c channel width, cm - W _c0 initial channel width, cm - v interstitial fluid velocity, cm/sec - v _k macroscopic velocity in transverse direction, cm/sec - v _y macroscopic velocity in longitudinal direction, cm/sec - v fluid velocity entering the medium, cm/sec - x _i transverse location of parcel at time t _i, cm - y _i longitudinal location of parcel at time t _i, cm - x microscopic movement in transverse direction, cm - y microscopic movement in longitudinal direction, cm Greek Letters t time increment, sec - ₀ overall dispersivity, cm - ² longitudinal variance of the distribution, cm² - porosity - _B bulk density, gm/cm³ - fractional grade of leachable material Currently with Center for Naval Analysis.     

Numerical simulation of the three-dimensional micro-scale dispersion of air-pollutants in regions with complex topography

 A computational method for predicting the dispersion of gaseous emissions in regions with complex topography is presented. Different meteorological conditions, source data, and gas compositions and specific weights are taken into consideration. The three-dimensional governing equations of fluid- and thermodynamics are numerically solved by means of the finite volume method. The kɛ turbulence model is utilized to account for the turbulent nature of the flow. The numerical results obtained by the proposed method show satisfactory agreement with both results obtained by other numerical methods and experimental and/or measured data. The presented method is applied to four test cases, including steady and unsteady problems, in order to illustrate its usefulness. Received on 1 March 2001     

不同重力场下颗粒冲击过程的离散元分析

开展不同重力场下颗粒材料冲击动力学研究有助于加深对颗粒运动机制的理解和深空探测任务的实施。本文采用离散元模拟对颗粒介质受球形冲击物的冲击过程进行了数值模拟,并与地球重力下冲击的试验结果进行对比验证。在此基础上,进一步研究了重力加速度对冲击物动力学的影响规律。计算结果表明,在所有重力加速度下,冲击物的穿透深度d与冲击速度v₀的关系可以用Poncelet模型表达;d与冲击物下落的总高度H表现为d～Hⁿ的幂律关系,当H<10 m时,d与H的幂率标度为0.322,而H>10 m时,d与H的幂率标度下降到0.211。此外,穿透深度小于冲击物半径时,重力加速度对冲击物减速过程无影响。在所有的重力加速度下,当冲击速度大于5 m/s时,冲击物的持续碰撞时间t_c是恒定的,且与重力的-1/2次方呈正比。     

Velocity and turbulence measurements in combustion systems

A laser-Doppler velocimeter is used in the measurement of high-temperature gas flows. A two-stage fluidization particle generator provides magnesium oxide particles to serve as optical scattering centers. The one-dimensional dual-beam system is frequency shifted to permit measurements of velocities up to 300 meters per second and turbulence intensities greater than 100 percent.Exiting flows from can-type gas turbine combustors and burners with pre-mixed oxy-acetylene flames are described in terms of the velocity, turbulence intensity, and temperature profiles.The results indicate the influence of the combustion process on turbulence.List of Symbols A exit area of combustor or burner (m²) - A/F mass air-fuel ratio - D exit diameter of combustor or burner (m) - M mass flow rate of gases (kg/s) - N _D number of Doppler bursts used in each velocity measurement - Q volumetric flow rate at T _r (m³/s) - R exit radius of combustor or burner (m) - R _1/2 distance from centerline to radius where the velocity is one-half of the local centerline velocity (m) - Re exit Reynolds number based on cold flow, QD/A - r distance from centerline of flow (m) - T temperature (°C) - T _CL centerline temperature (°C) - T _r inlet (cold) air temperature of combustor or burner (°C) - T.I. turbulence intensity, - mean velocity (m/s) - U _i instantaneous velocity individually realized by LDV (m/s) - mean velocity at centerline of flow (m/s) - mean square velocity fluctuation (m²/s² - x distance along centerline downstream of exit (m) - absolute viscosity at T _r (kg/(ms)) - density at T _r (kg/m³)     

Modeling Turbulent Droplet Motion in Vaporizing Sprays

The present article is concerned with the influence of turbulent gas-velocity fluctuations on both droplet dispersion and droplet-gas slip velocity in the context of spray simulation. The role of turbulence in generating slip and thus enhancing interphase heat and mass transfer has so far received little attention and is investigated in this work. A model for turbulent gas-velocity fluctuations along droplet trajectories is presented and is first tuned to reproduce elementary dispersion phenomena. It is then shown to give good results for more general dispersion problems as well as for slip velocities. As a fundamental source of information and for the purpose of model validation and comparison, direct numerical simulation (DNS) of droplet motion in homogeneous isotropic steady turbulence (HIST) is used. Dispersion of “injected” droplets (i.e. droplets under the influence of drift due to high injection velocity) as well as slip velocities for linear and nonlinear droplet drag are studied, and reasonable agreement is found with the model. The distributions of the slip velocity are found to be very similar for linear and highly nonlinear drag law. The present model is also used to investigate the influence of turbulence on droplet penetration. Comparison is made with an eddy-interaction model (the KIVA-2 model), which reveals various weaknesses of this model, in particular the underprediction of average slip velocity. The influence of slip due to turbulence on vaporization is shown for a fuel spray injected into a premix gas-turbine combustor. The classical eddy-interaction model is seen to underestimate the rate of vaporization due to the underprediction of slip. This revised version was published online in July 2006 with corrections to the Cover Date.     

Analysis of unsteady gas–liquid flows in a rectangular tank: Comparison of Euler–Eulerian and Euler–Lagrangian simulations

We study the dynamics of gas–liquid flows experimentally and computationally in a rectangular bubble column where the gas source is introduced at the corner. The flow in this reactor is complex and inherently unsteady in nature. The two-dimensional liquid phase velocity field is calculated by an Eulerian approach solving the unsteady Reynolds Averaged Navier Stokes equations. The conservation equations are closed using a two parameter turbulence model. The two-way coupling was accounted for by adding source terms in the conservation equations of the continuous phase to take into account the interaction with the dispersed phase. Bubble tracking is achieved through a Lagrangian approach. Here the equations of motion are solved taking into account the drag, pressure, buoyancy and gravity forces. The time-averaged flows along with the variables which characterize turbulence are analyzed for a wide range of gas flow-rates using Euler–Lagrangian simulations. These simulation predictions are validated with Euler–Eulerian simulations where the gas-phase distribution is captured as a void fraction and PIV experiments. The motion of bubbles induces turbulence in the flow. The applicability of two parameter models for turbulence like the standard k–ε model on time-averaged flow properties is addressed. From the results of the time averaged velocity field, turbulence intensity, turbulent viscosity and gas hold-up profiles, it is concluded that the Euler–Lagrangian model is applicable at lower gas flow-rates. The Euler–Eulerian approach was found to be valid at lower as well as higher gas flow-rates.     

Separation in a centrifuge-model

Summary The gas centrifuges which were used for isotope separation ^1,2) have cylindrical rotors. In consequence of the equality of the relaxation times for diffusion and gas motion, and taking into account the pressure-dependence of the diffusion constant, the time-dependent equation for a cylindrical centrifuge is rather complicated ³). Therefore a simple centrifuge-model where gas motion does not need to be considered is of interest. Though the theoretical maximum peripheral velocity for this model is up to 2-times higher than for a cylindrical rotor ¹), this interest is only mathematical.     

Penetration height correlations for non-aerated and aerated transverse liquid jets in supersonic cross flow

Experimental results on the mixing of non-aerated and aerated transverse liquid jet in supersonic cross flow (M = 1.5) are presented in this paper. The goal of this study is to investigate the effect of the gas/liquid mass ratio on the penetration and atomization of an aerated liquid jet in high speed cross flow and to develop correlations for the penetration heights. High speed imaging system was used in this study for the visualization of the injection of aerated liquid jet. The results show the effect of jet/cross flow momentum flux ratio, the gas/liquid mass ratio and the Ohnesorge number on the penetration of aerated liquid jet in supersonic cross-flow. New correlations of the spray penetration height for the non-aerated liquid jet (GLR = 0) and the net gain in spray penetration height for the aerated liquid jet (GLR > 0) are presented.     

Free Convection Boundary Layer Flow from a Heated Upward Facing Horizontal Flat Plate Embedded in a Porous Medium Filled by a Nanofluid with Convective Boundary Condition

The steady laminar incompressible free convective flow of a nanofluid over a permeable upward facing horizontal plate located in porous medium taking into account the thermal convective boundary condition is studied numerically. The nanofluid model used involves the effect of Brownian motion and the thermophoresis. Using similarity transformations the continuity, the momentum, the energy, and the nanoparticle volume fraction equations are transformed into a set of coupled similarity equations, before being solved numerically, by an implicit finite difference numerical method. Our analysis reveals that for a true similarity solution, the convective heat transfer coefficient related with the hot fluid and the mass transfer velocity must be proportional to x ^−2/3, where x is the horizontal distance along the plate from the origin. Effects of the various parameters on the dimensionless longitudinal velocity, the temperature, the nanoparticle volume fraction, as well as on the rate of heat transfer and the rate of nanoparticle volume fraction have been presented graphically and discussed. It is found that Lewis number, the Brownian motion, and the convective heat transfer parameters increase the heat transfer rate whilst the thermophoresis decreases the heat transfer rate. It is also found that Lewis number and the convective heat transfer parameter enhance the nanoparticle volume fraction rate whilst the thermophoresis parameter decreases nanoparticle volume fraction rate. A very good agreement is found between numerical results of the present article for special case and published results. This close agreement supports the validity of our analysis and the accuracy of the numerical computations.     

Lyapunov Functional Method for 1D Radiative and Reactive Viscous Gas Dynamics

We consider the compressible Navier–Stokes system for 1D-flows of a viscous heat-conducting gas, with the pressure law and a one-order kinetics to include radiative effects and reactive processes. The mass force and the ignition phenomenon are also taken into account. For large data and under general assumptions on the heat conductivity, we establish global-in-time bounds and exponential stabilization for solutions in L^q and H¹ norms. To this end, we construct new global Lyapunov functionals and show that they describe the dynamics of solutions for any t≧0. A short proof of the corresponding global existence is also included for completeness.