Spectral analysis of energy transfer in turbulent flows laden with heated particles

2017 ◽  
Vol 813 ◽  
pp. 1156-1175 ◽  
Author(s):  
H. Pouransari ◽  
H. Kolla ◽  
J. H. Chen ◽  
A. Mani

In this study we consider particle-laden turbulent flows with significant heat transfer between the two phases due to sustained heating of the particle phase. The sustained heat source can be due to particle heating via an external radiation source as in the particle-based solar receivers or an exothermic reaction in the particles. Our objective is to investigate the effects of fluid heating by a dispersed phase on the turbulence evolution. An important feature in such settings is the preferential clustering phenomenon which is responsible for non-uniform distribution of particles in the fluid medium. Particularly, when the ratio of particle inertial relaxation time to the turbulence time scale, namely the Stokes number, is of order unity, particle clustering is maximized, leading to thin regions of heat source similar to the flames in turbulent combustion. However, unlike turbulent combustion, a particle-laden system involves a wide range of clustering scales that is mainly controlled by particle Stokes number. To study these effects, we considered a decaying homogeneous isotropic turbulence laden with heated particles over a wide range of Stokes numbers. Using a low-Mach-number formulation for the fluid energy equation and a Lagrangian framework for particle tracking, we performed numerical simulations of this coupled system. We then applied a high-fidelity framework to perform spectral analysis of kinetic energy in a variable-density fluid. Our results indicate that particle heating can considerably influence the turbulence cascade. We show that the pressure-dilatation term introduces turbulent kinetic energy at a range of scales consistent with the scales observed in particle clusters. This energy is then transferred to high wavenumbers via the energy transfer term. For low and moderate levels of particle heating intensity, quantified by a parameter $\unicode[STIX]{x1D6FC}$ defined as the ratio of eddy time to mean temperature increase time, turbulence modification occurs primarily in the dilatational modes of the velocity field. However, as the heating intensity is increased, the energy transfer term converts energy from dilatational modes to divergence-free modes. Interestingly, as the heating intensity is increased, the net modification of turbulence by heating is observed dominantly in divergence-free modes; the portion of turbulence modification in dilatational modes diminishes with higher heating. Moreover, we show that modification of divergence-free modes is more pronounced at intermediate Stokes numbers corresponding to the maximum particle clustering. We also present the influence of heating intensity on the energy transfer term itself. This term crosses over from negative to positive values beyond a threshold wavenumber, showing the cascade of energy from large scales to small scales. The threshold is shown to shift to higher wavenumbers with increasing heating, indicating a growth in the total energy transfer from large scales to small scales. The fundamental energy transfer analysis presented in this paper provides insightful guidelines for subgrid-scale modelling and large-eddy simulation of heated particle-laden turbulence.

Author(s):  
Amilcare Porporato ◽  
Milad Hooshyar ◽  
Andrew D. Bragg ◽  
Gabriel Katul

Turbulent flows are out-of-equilibrium because the energy supply at large scales and its dissipation by viscosity at small scales create a net transfer of energy among all scales. This energy cascade is modelled by approximating the spectral energy balance with a nonlinear Fokker–Planck equation consistent with accepted phenomenological theories of turbulence. The steady-state contributions of the drift and diffusion in the corresponding Langevin equation, combined with the killing term associated with the dissipation, induce a stochastic energy transfer across wavenumbers. The fluctuation theorem is shown to describe the scale-wise statistics of forward and backward energy transfer and their connection to irreversibility and entropy production. The ensuing turbulence entropy is used to formulate an extended turbulence thermodynamics.


2020 ◽  
Vol 6 (4) ◽  
pp. eaaw6256
Author(s):  
Rory T. Cerbus ◽  
Chien-chia Liu ◽  
Gustavo Gioia ◽  
Pinaki Chakraborty

Turbulent flows are not only everywhere, but every turbulent flow is the same at small scales. The extraordinary simplification engendered by this “small-scale universality” is a hallmark of turbulence theory. However, on the basis of the restrictive assumptions invoked by A. N. Kolmogorov to demonstrate this universality, it is widely thought that only idealized turbulent flows conform to this framework. Using experiments and simulations that span a wide range of Reynolds number, we show that small-scale universality governs the spectral structure of a class of flows with no apparent ties to the idealized flows: transitional pipe flows. Our results not only extend the universality of Kolmogorov’s framework beyond expectation but also establish an unexpected link between transitional pipe flows and Kolmogorovian turbulence.


2009 ◽  
Author(s):  
Balaji Gopalan ◽  
Edwin Malkiel ◽  
Joseph Katz

The addition of dispersants, water and oil soluble surfactants that lower the interfacial tension of the crude oil, along with oceanic turbulence can breakdown oil spills into droplets. Knowledge of the dispersion rate of these droplets by oceanic turbulence is essential for the development of better models to assess the environmental impact of spills. The objective of this research is to study, experimentally, the dispersion of oil droplets in turbulent flows. The measurements are performed in a specialized laboratory facility that enables generation of carefully controlled, isotropic, homogeneous turbulence at a wide range of fully characterized intensities and length scales. The oil dispersion is visualized using high-speed inline digital holographic cinematography. Holographic data has been analyzed and Lagrangian statistics of droplet velocity, dispersion and acceleration has been calculated. As the relative size of the droplet diameter to the Kolmogorov length scale and its Stokes number increases, the acceleration autocorrelation shifts from dropping to zero faster than the fluid particles to slower.


Materials ◽  
2021 ◽  
Vol 14 (7) ◽  
pp. 1740
Author(s):  
Konrad Jakubowski ◽  
Manfred Heuberger ◽  
Rudolf Hufenus

The increasing interest in luminescent waveguides, applied as light concentrators, sensing elements, or decorative illuminating systems, is fostering efforts to further expand their functionality. Yarns and textiles based on a combination of distinct melt-spun polymer optical fibers (POFs), doped with individual luminescent dyes, can be beneficial for such applications since they enable easy tuning of the color of emitted light. Based on the energy transfer occurring between differently dyed filaments within a yarn or textile, the collective emission properties of such assemblies are adjustable over a wide range. The presented study demonstrates this effect using multicolor, meltspun, and photoluminescent POFs to measure their superimposed photoluminescent emission spectra. By varying the concentration of luminophores in yarn and fabric composition, the overall color of the resulting photoluminescent textiles can be tailored by the recapturing of light escaping from individual POFs. The ensuing color space is a mean to address the needs of specific applications, such as decorative elements and textile illumination by UV down-conversion.


2020 ◽  
Vol 234 (7-9) ◽  
pp. 1233-1250 ◽  
Author(s):  
Arrke J. Eskola ◽  
Mark A. Blitz ◽  
Michael J. Pilling ◽  
Paul W. Seakins ◽  
Robin J. Shannon

AbstractThe rate coefficient for the unimolecular decomposition of CH3OCH2, k1, has been measured in time-resolved experiments by monitoring the HCHO product. CH3OCH2 was rapidly and cleanly generated by 248 nm excimer photolysis of oxalyl chloride, (ClCO)2, in an excess of CH3OCH3, and an excimer pumped dye laser tuned to 353.16 nm was used to probe HCHO via laser induced fluorescence. k1(T,p) was measured over the ranges: 573–673 K and 0.1–4.3 × 1018 molecule cm−3 with a helium bath gas. In addition, some experiments were carried out with nitrogen as the bath gas. Ab initio calculations on CH3OCH2 decomposition were carried out and a transition-state for decomposition to CH3 and H2CO was identified. This information was used in a master equation rate calculation, using the MESMER code, where the zero-point-energy corrected barrier to reaction, ΔE0,1, and the energy transfer parameters, ⟨ΔEdown⟩ × Tn, were the adjusted parameters to best fit the experimental data, with helium as the buffer gas. The data were combined with earlier measurements by Loucks and Laidler (Can J. Chem.1967, 45, 2767), with dimethyl ether as the third body, reinterpreted using current literature for the rate coefficient for recombination of CH3OCH2. This analysis returned ΔE0,1 = (112.3 ± 0.6) kJ mol−1, and leads to $k_{1}^{\infty}(T)=2.9\times{10^{12}}$ (T/300)2.5 exp(−106.8 kJ mol−1/RT). Using this model, limited experiments with nitrogen as the bath gas allowed N2 energy transfer parameters to be identified and then further MESMER simulations were carried out, where N2 was the buffer gas, to generate k1(T,p) over a wide range of conditions: 300–1000 K and N2 = 1012–1025 molecule cm−3. The resulting k1(T,p) has been parameterized using a Troe-expression, so that they can be readily be incorporated into combustion models. In addition, k1(T,p) has been parametrized using PLOG for the buffer gases, He, CH3OCH3 and N2.


2001 ◽  
Author(s):  
X. Ai ◽  
B. Q. Li

Abstract Turbulent magnetically flows occur in a wide range of material processing systems involving electrically conducting melts. This paper presents a parallel higher order scheme for the direct numerical simulation of turbulent magnetically driven flows in induction channels. The numerical method is based on the higher order finite difference algorithm, which enjoys the spectral accuracy while minimizing the computational intensity. This, coupled with the parallel computing strategy, provides a very useful means to simulate turbulent flows. The higher order finite difference formulation of magnetically driven flow problems is described in this paper. The details of the parallel algorithm and its implementation for the simulations on parallel machines are discussed. The accuracy and numerical performance of the higher order finite difference scheme are assessed in comparison with the spectral method. The examples of turbulent magnetically driven flows in induction channels and pressure gradient driven flows in regular channels are given, and the computed results are compared with experimental measurements wherever possible.


2021 ◽  
Author(s):  
Stefano Berti ◽  
Guillaume Lapeyre

<p>Oceanic motions at scales larger than few tens of km are quasi-horizontal due to seawater stratification and Earth’s rotation and are characterized by quasi-two-dimensional turbulence. At scales around 300 km (in the mesoscale range), coherent vortices contain most of the kinetic energy in the ocean. At scales around 10 km (in the submesoscale range) the flow is populated by smaller eddies and filamentary structures associated with intense gradients (e.g. of temperature), which play an important role in both physical and biogeochemical budgets. Such small scales are found mainly in the weakly stratified mixed layer, lying on top of the more stratified thermocline. Submesoscale dynamics should strongly depend on the seasonal cycle and the associated mixed-layer instabilities. The latter are particularly relevant in winter and are responsible for the generation of energetic small scales that are not trapped at the surface, as those arising from mesoscale-driven processes, but extend down to the thermocline. The knowledge of the transport properties of oceanic flows at depth, which is essential to understand the coupling between surface and interior dynamics, however, is still limited.</p><p>By means of numerical simulations, we explore Lagrangian pair dispersion in turbulent flows from a quasi-geostrophic model consisting in two coupled fluid layers (representing the mixed layer and the thermocline) with different stratification. Such a model has been previously shown to give rise to both meso and submesoscale instabilities and subsequent turbulent dynamics that compare well with observations of wintertime submesoscale flows. We focus on the identification of different dispersion regimes and on the possibility to relate the characteristics of the spreading process at the surface and at depth, which is relevant to assess the possibility of inferring the dynamical features of deeper flows from the experimentally more accessible (e.g. by satellite altimetry) surface ones.</p><p>Using different statistical indicators, we find a clear transition of dispersion regime with depth, which is generic and can be related to the statistical features of the turbulent flows. The spreading process is local (namely, governed by eddies of the same size as the particle separation distance) at the surface. In the absence of a mixed layer it rapidly changes to nonlocal (meaning essentially driven by the largest eddies) at small depths, while in the opposite case this only occurs at larger depths, below the mixed layer. We then identify the origin of such behavior in the existence of fine-scale energetic structures due to mixed-layer instabilities. We further discuss the effect of vertical shear and address the properties of the relative motion of subsurface particles with respect to surface ones. In the absence of a mixed layer, the properties of the spreading process are found to rapidly decorrelate from those at the surface, but the relation between the surface and subsurface dispersion appears to be largely controlled by vertical shear. In the presence of mixed-layer instabilities, instead, the statistical properties of dispersion at the surface are found to be a good proxy for those in the whole mixed layer.</p>


2015 ◽  
Vol 725-726 ◽  
pp. 1255-1260
Author(s):  
Tamara Daciuk ◽  
Vera Ulyasheva

Numerical experiment has been successfully used during recent 10-15 years to solve a wide range of thermal and hydrogasodynamic tasks. Application of mathematical modeling used to design the ventilation systems for production premises characterized by heat emission may be considered to be an effective method to obtain reasonable solutions. Results of calculation performed with numerical solution of ventilation tasks depend on turbulence model selection. Currently a large number of different turbulence models used to calculate turbulent flows are known. Testing and definition of applicability limits for semiempirical models of turbulence should be considered to be a preliminary stage of calculation. This article presents results of test calculations pertaining to thermal air process modeling in premises characterized by presence of heat emission sources performed with employment of different models of turbulence. Besides, analysis of calculation results and comparison with field measurements data are presented.


1991 ◽  
Vol 113 (1) ◽  
pp. 34-41 ◽  
Author(s):  
G. J. Yoo ◽  
R. M. C. So ◽  
B. C. Hwang

Internal rotating boundary-layer flows are strongly influenced by large circumferential strain and the turbulence field is anisotropic. This is especially true in the entry region of a rotating pipe where the flow is three dimensional, the centrifugal force due to fluid rotation is less important, and the circumferential strain created by surface rotation has a significant effect on the turbulence field near the wall. Consequently, viscous effects cannot be neglected in the near-wall region. Several low-Reynolds-number turbulence closures are proposed for the calculation of developing rotating pipe flows. Some are two-equation closures with and without algebraic stress correction, while others are full Reynolds-stress closures. It is found that two-equation closures with and without algebraic stress correction are totally inadequate for this three-dimensional flow, while Reynolds-stress closures give results that are in good agreement with measurements over a wide range of rotation numbers.


Sign in / Sign up

Export Citation Format

Share Document