Numerical Simulation of Heavy Particle Dispersion—Scale Ratio and Flow Decay Considerations

1994 ◽  
Vol 116 (1) ◽  
pp. 154-163 ◽  
Author(s):  
Lian-Ping Wang ◽  
David E. Stock

Lagrangian statistical quantities related to the dispersion of heavy particles were studied numerically by following particle trajectories in a random flow generated by Fourier modes. An experimental fluid velocity correlation was incorporated into the flow. Numerical simulation was performed with the use of nonlinear drag. The simulation results for glass beads in a nondecaying turbulent air showed a difference between the horizontal dispersion coefficient and vertical dispersion coefficient. This difference was related to the differences of both the velocity scale and the time scale between the two direction. It was shown that for relatively small particle sizes the particle time scale ratio dominates the value of the diffusivity ratio. For large particles, the velocity scale ratio reaches a value of 1/2 and thus fully determines the diffusivity ratio. Qualitative explanation was provided to support the numerical findings. The dispersion data for heavy particles in grid-generated turbulences were successfully predicted by the simulation when flow decay was considered. As a result of the reduction in effective inertia and the increase in effective drift caused by the flow decay, the particle dispersion coefficient in decaying flow decreases with downstream location. The particle rms fluctuation velocity has a slower decay rate than the fluid rms velocity if the drift parameter is large. It was also found that the drift may substantially reduce the particle rms velocity.

1989 ◽  
Vol 111 (4) ◽  
pp. 492-494 ◽  
Author(s):  
Y. Zhuang ◽  
J. D. Wilson ◽  
E. P. Lozowski

For many purposes it is useful to be able to mimic the paths of heavy particles in a turbulent flow. This paper gives a simple procedure by which this may be achieved, provided particle spin is not important and under the restriction that the ratio of particle to fluid density exceeds about 1000. The procedure is related to the models of Faeth (1986) and Hunt and Nalpanis (1985). Simulation of the experiments of Snyder and Lumley (1971) yielded acceptable agreement with the observed rate of heavy particle dispersion.


1992 ◽  
Vol 114 (1) ◽  
pp. 100-106 ◽  
Author(s):  
Lian-Ping Wang ◽  
D. E. Stock

Numerical experiments can be used to study heavy particle dispersion by tracking particles through a numerically generated instantaneous turbulent flow field. In this manner, data can be generated to supplement physical experiments. To perform the numerical experiments efficiently and accurately, the time step used when tracking the particles through the fluid must be chosen correctly. After finding a suitable time step for one particular simulation, the time step must be reduced as the total integration time increases and as the free-fall velocity of the particle increases. Based on the numerical calculations, we suggest that the nonlinear drag be included in a numerical simulation if the ratio of the particle’s Stokes free-fall velocity to the fluid rms velocity is greater than two.


1997 ◽  
Vol 119 (1) ◽  
pp. 170-179 ◽  
Author(s):  
Renwei Mei ◽  
R. J. Adrian ◽  
T. J. Hanratty

The analysis of Reeks (1977) for particle dispersion in isotropic turbulence is extended so as to include a nonlinear drag law. The principal issue is the evaluation of the inertial time constants, βα−1, and the mean slip. Unlike what is found for the Stokesian drag, the time constants are functions of the slip velocity and are anisotropic. For settling velocity, VT, much larger than root-mean-square of the fluid velocity fluctuations, u0, the mean slip is given by VT. For VT→0, the mean slip is related to turbulent velocity fluctuation by assuming that fluctuations in βα are small compared to the mean value. An interpolation formula is used to evaluate βα and VT in regions intermediate between conditions of VT→0 and VT≫ u0. The limitations of the analysis are explored by carrying out a Monte-Carlo simulation for particle motion in a pseudo turbulence described by a Gaussian distribution and Kraichnan’s (1970) energy spectrum.


1980 ◽  
Vol 99 (3) ◽  
pp. 545-573 ◽  
Author(s):  
Z. Warhaft

The effect of homogeneous strain on passive scalar fluctuations, and the resultant evolution of the scalar field when the strain is removed, is experimentally studied by passing thermal fluctuations in decaying grid turbulence through a four-to-one axisymmetric contraction. Using amandoline(Warhaft & Lumley 1978a) to vary the scale size of the initial thermal fluctuations and hence the pre-contraction mechanical/thermal time-scale ratio,r, it is shown, for values ofrgreater than unity, that asris increased so is the post-contraction thermal decay rate, i.e. the contraction does not cause the thermal-fluctuation decay rate to equilibrate to a constant value. In these experiments the post-contraction thermal decay rate is always greater than the pre-contraction decay rate, i.e. the contraction accelerates the thermal-fluctuation decay. Moreover, the mechanical/thermal time-scale ratio in the post-contraction region is driven further from unity. In terms of scale size the uniform strain has the effect of increasing the thermal length scale by an amount equal in value to the contraction ratio if the pre-contraction thermal length scale is comparable to that of the pre-contraction velocity scale. However, if the pre-contraction thermal length scale is smaller than the pre-contraction velocity scale then the effect of the contraction on the thermal scale is less marked. The contraction induces significant negative cross-correlation ρuθbetween the longitudinal velocityuand thermal fluctuations θ even if the pre-contraction cross-correlation is close to zero. The magnitude of ρuθand hence the post-contraction heat flux is varied and the coherence structure is studied. It is shown that the thermal-fluctuation decay rate is insensitive to the magnitude of the heat flux, the latter of which decays rapidly compared to the relatively slow decay of turbulence energy in the post-contraction region. It is also shown that ρuθtends towards zero in this axisymmetric homogeneous flow at a faster rate than in isotropic turbulence. In accord with previous investigations, the return toward isotropy of the velocity field is very slow.


2006 ◽  
Vol 45 (9) ◽  
pp. 1332-1349 ◽  
Author(s):  
T. Bouvet ◽  
J. D. Wilson ◽  
A. Tuzet

Abstract This paper presents new observations of deposition of heavy particles (glass beads of gravitational settling velocity 8.7 cm s−1) within an undisturbed flow and within a flow disturbed by a porous windbreak fence. These data are then used to diagnose the capability of a Lagrangian stochastic (LS) particle trajectory model, which simulates heavy particle dispersion. The model is based on existing parameterizations and is coupled to a wind model based on a Reynolds stress turbulence closure that provides computed fields of wind statistics. The deposition rates, as simulated by the model, match the observation within E = 30% of accuracy, with E being the root-mean-square error normalized by the peak value on the deposition swath. These results suggest that the LS model handles properly the heterogeneities of the flow and that the heuristic adjustments made to account for the inertia of heavy particles are useful approximations. The model consequently proves to be a valuable tool to investigate the patterns of dispersion about an obstacle.


2021 ◽  
Vol 2021 (2) ◽  
Author(s):  
Arushi Bodas ◽  
Soubhik Kumar ◽  
Raman Sundrum

Abstract Non-analyticity in co-moving momenta within the non-Gaussian bispectrum is a distinctive sign of on-shell particle production during inflation, presenting a unique opportunity for the “direct detection” of particles with masses as large as the inflationary Hubble scale (H). However, the strength of such non-analyticity ordinarily drops exponentially by a Boltzmann-like factor as masses exceed H. In this paper, we study an exception provided by a dimension-5 derivative coupling of the inflaton to heavy-particle currents, applying it specifically to the case of two real scalars. The operator has a “chemical potential” form, which harnesses the large kinetic energy scale of the inflaton, $$ {\overset{\cdot }{\phi}}_0^{1/2}\approx 60H $$ ϕ ⋅ 0 1 / 2 ≈ 60 H , to act as an efficient source of scalar particle production. Derivative couplings of inflaton ensure radiative stability of the slow-roll potential, which in turn maintains (approximate) scale-invariance of the inflationary correlations. We show that a signal not suffering Boltzmann suppression can be obtained in the bispectrum with strength fNL ∼ $$ \mathcal{O} $$ O (0.01–10) for an extended range of scalar masses $$ \lesssim {\overset{\cdot }{\phi}}_0^{1/2} $$ ≲ ϕ ⋅ 0 1 / 2 , potentially as high as 1015 GeV, within the sensitivity of upcoming LSS and more futuristic 21-cm experiments. The mechanism does not invoke any particular fine-tuning of parameters or breakdown of perturbation-theoretic control. The leading contribution appears at tree-level, which makes the calculation analytically tractable and removes the loop-suppression as compared to earlier chemical potential studies of non-zero spins. The steady particle production allows us to infer the effective mass of the heavy particles and the chemical potential from the variation in bispectrum oscillations as a function of co-moving momenta. Our analysis sets the stage for generalization to heavy bosons with non-zero spin.


2014 ◽  
Vol 1008-1009 ◽  
pp. 850-860 ◽  
Author(s):  
Zhou Wei Zhang ◽  
Jia Xing Xue ◽  
Ya Hong Wang

A calculation method for counter-current type coil-wound heat exchanger is presented for heat exchange process. The numerical simulation method is applied to determine the basic physical parameters of wound bundles. By controlling the inlet fluid velocity varying in coil-wound heat exchanger to program and calculate the iterative process. The calculation data is analyzed by comparison of numerical result and the unit three dimensional pipe bundle model was built. Studies show that the introduction of numerical simulation can simplify the pipe winding process and accelerate the calculation and design of overall configuration in coil-wound heat exchanger. This method can be applied to the physical modeling and heat transfer calculation of pipe bundles in coil wound heat exchanger, program to calculate the complex heat transfer changing with velocity and other parameters, and optimize the overall design and calculation of spiral bundles.


1993 ◽  
Vol 115 (4) ◽  
pp. 751-759 ◽  
Author(s):  
Tsuneaki Ishima ◽  
Koichi Hishida ◽  
Masanobu Maeda

A particle dispersion has been experimentally investigated in a two-dimensional mixing layer with a large relative velocity between particle and gas-phase in order to clarify the effect of particle residence time on particle dispersion. Spherical glass particles 42, 72, and 135 μm in diameter were loaded directly into the origin of the shear layer. Particle number density and the velocities of both particle and gas phase were measured by a laser Doppler velocimeter with modified signal processing for two-phase flow. The results confirmed that the characteristic time scale of the coherent eddy apparently became equivalent to a shorter characteristic time scale due to a less residence time. The particle dispersion coefficients were well correlated to the extended Stokes number defined as the ratio of the particle relaxation time to the substantial eddy characteristic time scale which was evaluated by taking account of the particle residence time.


2016 ◽  
Vol 87 ◽  
pp. 80-89 ◽  
Author(s):  
Xiaoke Ku ◽  
Jianzhong Lin ◽  
Martin van Sint Annaland ◽  
Rob Hagmeijer

2011 ◽  
Vol 133 (8) ◽  
Author(s):  
Dong-Hyeog Yoon ◽  
Kyung-Soo Yang ◽  
Klaus Bremhorst

Characteristics of turbulent mass transfer around a rotating circular cylinder have been investigated by Direct Numerical Simulation. The concentration field was computed for three different cases of Schmidt number, Sc = 1, 10 and 100 at ReR* = 336. Our results confirm that the thickness of the Nernst diffusion layer decreases as Sc increases. Wall-limiting behavior within the diffusion layer was examined and compared with that of channel flow. Concentration fluctuation time scale was found to scale with r+2, while the time scale ratio nearly equals the Schmidt number throughout the diffusion layer. Scalar modeling closure constants based on gradient diffusion models were found to vary considerably within the diffusion layer. Results of an octant analysis show the significant role played by the ejection and sweep events just as is found for flat plate, channel, and pipe flow boundary layers. Turbulence budgets revealed a strong Sc dependence of turbulent scalar transport.


Sign in / Sign up

Export Citation Format

Share Document