Droplet–turbulence interactions in low-Mach-number homogeneous shear two-phase flows

1998 ◽  
Vol 367 ◽  
pp. 163-203 ◽  
Author(s):  
FARZAD MASHAYEK
Keyword(s):  
Kinetic Energy ◽  
Turbulent Flows ◽  
Carrier Phase ◽  
Mass Fraction ◽  
Turbulence Kinetic Energy ◽  
Droplet Temperature ◽  
Low Mach Number ◽  
Homogeneous Shear ◽  
Initial Droplet ◽  
The Mean

Several important issues pertaining to dispersion and polydispersity of droplets in turbulent flows are investigated via direct numerical simulation (DNS). The carrier phase is considered in the Eulerian context, the dispersed phase is tracked in the Lagrangian frame and the interactions between the phases are taken into account in a realistic two-way (coupled) formulation. The resulting scheme is applied for extensive DNS of low-Mach-number, homogeneous shear turbulent flows laden with droplets. Several cases with one- and two-way couplings are considered for both non-evaporating and evaporating droplets. The effects of the mass loading ratio, the droplet time constant, and thermodynamic parameters, such as the droplet specific heat, the droplet latent heat of evaporation, and the boiling temperature, on the turbulence and the droplets are investigated. The effects of the initial droplet temperature and the initial vapour mass fraction in the carrier phase are also studied. The gravity effects are not considered as the numerical methodology is only applicable in the absence of gravity. The evolution of the turbulence kinetic energy and the mean internal energy of both phases is studied by analysing various terms in their transport equations. The results for the non-evaporating droplets show that the presence of the droplets decreases the turbulence kinetic energy of the carrier phase while increasing the level of anisotropy of the flow. The droplet streamwise velocity variance is larger than that of the fluid, and the ratio of the two increases with the increase of the droplet time constant. Evaporation increases both the turbulence kinetic energy and the mean internal energy of the carrier phase by mass transfer. In general, evaporation is controlled by the vapour mass fraction gradient around the droplet when the initial temperature difference between the phases is negligible. In cases with small initial droplet temperature, on the other hand, the convective heat transfer is more important in the evaporation process. At long times, the evaporation rate approaches asymptotic values depending on the values of various parameters. It is shown that the evaporation rate is larger for droplets residing in high-strain-rate regions of the flow, mainly due to larger droplet Reynolds numbers in these regions. For both the evaporating and the non-evaporating droplets, the root mean square (r.m.s.) of the temperature fluctuations of both phases becomes independent of the initial droplet temperature at long times. Some issues relevant to modelling of turbulent flows laden with droplets are also discussed.

scholarly journals Turbulence Adjustment and Scaling in an Offshore Convective Internal Boundary Layer: A CASPER Case Study

2020 ◽  
Vol 77 (5) ◽  
pp. 1661-1681
Author(s):  
Qingfang Jiang ◽  
Qing Wang ◽  
Shouping Wang ◽  
Saša Gaberšek
Keyword(s):  
Boundary Layer ◽  
Surface Layer ◽  
Kinetic Energy ◽  
Wind Shear ◽  
Vertical Wind Shear ◽  
Internal Boundary Layer ◽  
Turbulence Kinetic Energy ◽  
Internal Boundary ◽  
Field Observations ◽  
The Mean

Abstract The characteristics of a convective internal boundary layer (CIBL) documented offshore during the East Coast phase of the Coupled Air–Sea Processes and Electromagnetic Ducting Research (CASPER-EAST) field campaign has been examined using field observations, a coupled mesoscale model (i.e., Navy’s COAMPS) simulation, and a couple of surface-layer-resolving large-eddy simulations (LESs). The Lagrangian modeling approach has been adopted with the LES domain being advected from a cool and rough land surface to a warmer and smoother sea surface by the mean offshore winds in the CIBL. The surface fluxes from the LES control run are in reasonable agreement with field observations, and the general CIBL characteristics are consistent with previous studies. According to the LESs, in the nearshore adjustment zone (i.e., fetch < 8 km), the low-level winds and surface friction velocity increase rapidly, and the mean wind profile and vertical velocity skewness in the surface layer deviate substantially from the Monin–Obukhov similarity theory (MOST) scaling. Farther offshore, the nondimensional vertical wind shear and scalar gradients and higher-order moments are consistent with the MOST scaling. An elevated turbulent layer is present immediately below the CIBL top, associated with the vertical wind shear across the CIBL top inversion. Episodic shear instability events occur with a time scale between 10 and 30 min, leading to the formation of elevated maxima in turbulence kinetic energy and momentum fluxes. During these events, the turbulence kinetic energy production exceeds the dissipation, suggesting that the CIBL remains in nonequilibrium.

An explicit algebraic Reynolds-stress and scalar-flux model for stably stratified flows

2013 ◽  
Vol 723 ◽  
pp. 91-125 ◽  
Author(s):  
W. M. J. Lazeroms ◽  
G. Brethouwer ◽  
S. Wallin ◽  
A. V. Johansson
Keyword(s):  
Heat Flux ◽  
Kinetic Energy ◽  
Temperature Gradient ◽  
Reynolds Stress ◽  
Turbulent Heat Flux ◽  
Turbulent Heat ◽  
Homogeneous Shear ◽  
Self Consistent ◽  
The Mean ◽  
Homogeneous Shear Flow

AbstractThis work describes the derivation of an algebraic model for the Reynolds stresses and turbulent heat flux in stably stratified turbulent flows, which are mutually coupled for this type of flow. For general two-dimensional mean flows, we present a correct way of expressing the Reynolds-stress anisotropy and the (normalized) turbulent heat flux as tensorial combinations of the mean strain rate, the mean rotation rate, the mean temperature gradient and gravity. A system of linear equations is derived for the coefficients in these expansions, which can easily be solved with computer algebra software for a specific choice of the model constants. The general model is simplified in the case of parallel mean shear flows where the temperature gradient is aligned with gravity. For this case, fully explicit and coupled expressions for the Reynolds-stress tensor and heat-flux vector are given. A self-consistent derivation of this model would, however, require finding a root of a polynomial equation of sixth-order, for which no simple analytical expression exists. Therefore, the nonlinear part of the algebraic equations is modelled through an approximation that is close to the consistent formulation. By using the framework of a$K\text{{\ndash}} \omega $model (where$K$is turbulent kinetic energy and$\omega $an inverse time scale) and, where needed, near-wall corrections, the model is applied to homogeneous shear flow and turbulent channel flow, both with stable stratification. For the case of homogeneous shear flow, the model predicts a critical Richardson number of 0.25 above which the turbulent kinetic energy decays to zero. The channel-flow results agree well with DNS data. Furthermore, the model is shown to be robust and approximately self-consistent. It also fulfils the requirements of realizability.

A k-ε Model for Particle-Laden Turbulent Flows

2009 ◽  
Author(s):  
J. D. Schwarzkopf ◽  
C. T. Crowe ◽  
P. Dutta
Keyword(s):  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
Turbulent Flows ◽  
Carrier Phase ◽  
Critical Role ◽  
New Production ◽  
Mean Velocity ◽  
Production Term ◽  
Dissipation Term ◽  
Velocity Gradients

A dissipation transport equation for the carrier phase of particle-laden turbulent flows was recently developed. This equation shows a new production of dissipation term due to the presence of particles that is related to the velocity difference between the particle and the surrounding fluid. In the development, it was assumed that each coefficient was the sum of the coefficient for single phase flow and a coefficient quantifying the contribution of the particulate phase. The coefficient for the new production term (due to the presence of particles) was found from homogeneous turbulence generation by particles and the coefficient for the dissipation of dissipation term was analyzed using DNS. A numerical model was developed and applied to particles falling in a channel of downward turbulent air flow. Boundary conditions were also developed to ensure that the production of turbulent kinetic energy due to mean velocity gradients and particle surfaces balanced with the turbulent dissipation near the wall. The turbulent kinetic energy is compared with experimental data. The results show attenuation of turbulent kinetic energy with increased particle loading; however the model does under predict the turbulent kinetic energy near the center of the channel. To understand the effect of this additional production of dissipation term (due to particles), the coefficients associated with the production of dissipation due to mean velocity gradients and particle surfaces are varied to assess the effects of the dispersed phase on the carrier phase turbulent kinetic energy across the channel. The results show that this additional term plays a significant role in predicting the turbulent kinetic energy and a reason for under predicting the turbulent kinetic energy near the center of the channel is discussed. It is concluded that the dissipation coefficients play a critical role in predicting the turbulent kinetic energy in particle-laden turbulent flows.

scholarly journals Energy-consistent entrainment relations for jets and plumes

2015 ◽  
Vol 782 ◽  
pp. 333-355 ◽  
Author(s):  
Maarten van Reeuwijk ◽  
John Craske
Keyword(s):  
Kinetic Energy ◽  
Richardson Number ◽  
Turbulence Kinetic Energy ◽  
Mean Flow ◽  
Unified Framework ◽  
First Order ◽  
Entrainment Coefficient ◽  
The Mean ◽  
Flow Turbulence ◽  
Pure Plume

We discuss energetic restrictions on the entrainment coefficient${\it\alpha}$for axisymmetric jets and plumes. The resulting entrainment relation includes contributions from the mean flow, turbulence and pressure, fundamentally linking${\it\alpha}$to the production of turbulence kinetic energy, the plume Richardson number$\mathit{Ri}$and the profile coefficients associated with the shape of the buoyancy and velocity profiles. This entrainment relation generalises the work by Kaminskiet al. (J. Fluid Mech., vol. 526, 2005, pp. 361–376) and Fox (J. Geophys. Res., vol. 75, 1970, pp. 6818–6835). The energetic viewpoint provides a unified framework with which to analyse the classical entrainment models implied by the plume theories of Mortonet al.(Proc. R. Soc. Lond.A, vol. 234, 1955, pp. 1–23) and Priestley & Ball (Q. J. R. Meteorol. Soc., vol. 81, 1954, pp. 144–157). Data for pure jets and plumes in unstratified environments indicate that to first order the physics is captured by the Priestley and Ball entrainment model, implying that (1) the profile coefficient associated with the production of turbulence kinetic energy has approximately the same value for pure plumes and jets, (2) the value of${\it\alpha}$for a pure plume is roughly a factor of$5/3$larger than for a jet and (3) the enhanced entrainment coefficient in plumes is primarily associated with the behaviour of the mean flow and not with buoyancy-enhanced turbulence. Theoretical suggestions are made on how entrainment can be systematically studied by creating constant-$\mathit{Ri}$flows in a numerical simulation or laboratory experiment.

Full-coverage film cooling. Part 2. Prediction of the recovery-region hydrodynamics

1980 ◽  
Vol 101 (1) ◽  
pp. 159-178 ◽  
Author(s):  
S. Yavuzkurt ◽  
R. J. Moffat ◽  
W. M. Kays
Keyword(s):  
Boundary Layer ◽  
Kinetic Energy ◽  
Film Cooling ◽  
Turbulence Kinetic Energy ◽  
Turbulent Shear ◽  
Mixing Length ◽  
Two Dimensional ◽  
Mean Velocity ◽  
Full Coverage ◽  
The Mean

Hydrodynamic data are reported in the companion paper (Yavuzkurt, Moffat & Kays 1980) for a full-coverage film-cooling situation, both for the blown and the recovery regions. Values of the mean velocity, the turbulent shear stress, and the turbulence kinetic energy were measured at various locations, both within the blown region and in the recovery region. The present paper is concerned with an analysis of the recovery region only. Examination of the data suggested that the recovery-region hydrodynamics could be modelled by considering that a new boundary layer began to grow immediately after the cessation of blowing. Distributions of the Prandtl mixing length were calculated from the data using the measured values of mean velocity and turbulent shear stresses. The mixing-length distributions were consistent with the notion of a dual boundary-layer structure in the recovery region. The measured distributions of mixing length were described by using a piecewise continuous but heuristic fit, consistent with the concept of two quasi-independent layers suggested by the general appearance of the data. This distribution of mixing length, together with a set of otherwise normal constants for a two-dimensional boundary layer, successfully predicted all of the observed features of the flow. The program used in these predictions contains a one-equation model of turbulence, using turbulence kinetic energy with an algebraic mixing length. The program is a two-dimensional, finite-difference program capable of predicting the mean velocity and turbulence kinetic energy profiles based upon initial values, boundary conditions, and a closure condition.

scholarly journals Biases in Thorpe-Scale Estimates of Turbulence Dissipation. Part II: Energetics Arguments and Turbulence Simulations

2015 ◽  
Vol 45 (10) ◽  
pp. 2522-2543 ◽  
Author(s):  
Alberto Scotti
Keyword(s):  
Kinetic Energy ◽  
Potential Energy ◽  
Turbulent Flows ◽  
Companion Paper ◽  
Mean Flow ◽  
Available Potential Energy ◽  
Ozmidov Scale ◽  
Thorpe Scale ◽  
The Mean ◽  
Shear Driven Flows

AbstractThis paper uses the energetics framework developed by Scotti and White to provide a critical assessment of the widely used Thorpe-scale method, which is used to estimate dissipation and mixing rates in stratified turbulent flows from density measurements along vertical profiles. This study shows that the relevant displacement scale in general is not the rms value of the Thorpe displacement. Rather, the displacement field must be Reynolds decomposed to separate the mean from the turbulent component, and it is the turbulent component that ought to be used to diagnose mixing and dissipation. In general, the energetics of mixing in an overall stably stratified flow involves potentially complex exchanges among the available potential energy and kinetic energy associated with the mean and turbulent components of the flow. The author considers two limiting cases: shear-driven mixing, where mixing comes at the expense of the mean kinetic energy of the flow, and convective-driven mixing, which taps the available potential energy of the mean flow to drive mixing. In shear-driven flows, the rms of the Thorpe displacement, known as the Thorpe scale is shown to be equivalent to the turbulent component of the displacement. In this case, the Thorpe scale approximates the Ozmidov scale, or, which is the same, the Thorpe scale is the appropriate scale to diagnose mixing and dissipation. However, when mixing is driven by the available potential energy of the mean flow (convective-driven mixing), this study shows that the Thorpe scale is (much) larger than the Ozmidov scale. Using the rms of the Thorpe displacement overestimates dissipation and mixing, since the amount of turbulent available potential energy (measured by the turbulent displacement) is only a fraction of the total available potential energy (measured by the Thorpe scale). Corrective measures are discussed that can be used to diagnose mixing from knowledge of the Thorpe displacement. In a companion paper, Mater et al. analyze field data and show that the Thorpe scale can indeed be much larger than the Ozmidov scale.

Simulation of Wind Effects Embracing a Complex Shape Super High-Rise Steel TV Tower Using CFD

2011 ◽  
Vol 201-203 ◽  
pp. 2807-2813
Author(s):  
Ya Ya Tan ◽  
Chun Xiang Li
Keyword(s):  
Kinetic Energy ◽  
Complex Shape ◽  
Measurement Data ◽  
Turbulence Kinetic Energy ◽  
High Rise ◽  
Distribution Rule ◽  
Pressure Region ◽  
Mean Pressure ◽  
The Mean ◽  
Inflow Conditions

Full scale Ningxia super high-rise steel TV tower (NXTVT) model is established in the numerical simulation to investigate the flow field around it. And then, the UDF program of mean wind section plane, turbulence kinetic energy and turbulence dissipation rate were compiled. The flow fields under these two typical inflow conditions, 30 degree is chosen as the field measurement data. The results show that both the RNG model and the Realizable model can simulate the distribution rule and similarity. Both can give reasonable turbulence kinetic energy distribution in the region of windward sharp edge, where flow separation occurs. Subsequently, the mean pressure field with certain precision can be obtained in the highly turbulent negative pressure region.

scholarly journals Multiscale Shear Forcing of Turbulence in the Nocturnal Boundary Layer: A Statistical Analysis

2020 ◽  
Author(s):  
Vyacheslav Boyko ◽  
Nikki Vercauteren
Keyword(s):  
Boundary Layer ◽  
Kinetic Energy ◽  
Nocturnal Boundary Layer ◽  
Turbulence Kinetic Energy ◽  
Eddy Correlation ◽  
Model Parameters ◽  
Intermittent Turbulence ◽  
Velocity Signal ◽  
The Mean ◽  
Vertical Gradients

AbstractThe lower nocturnal boundary layer is governed by intermittent turbulence which is thought to be triggered by sporadic activity of so-called sub-mesoscale motions in a complex way. We analyze intermittent turbulence based on an assumed relation between the vertical gradients of the sub-mean scales and turbulence kinetic energy. We analyze high-resolution nocturnal eddy-correlation data from 30-m tower collected during the Fluxes over Snow Surfaces II field program. The non-turbulent velocity signal is decomposed using a discrete wavelet transform into three ranges of scales interpreted as the mean, jet and sub-mesoscales. The vertical gradients of the sub-mean scales are estimated using finite differences. The turbulence kinetic energy is modelled as a discrete-time autoregressive process with exogenous variables, where the latter ones are the vertical gradients of the sub-mean scales. The parameters of the discrete model evolve in time depending on the locally-dominant turbulence-production scales. The three regimes with averaged model parameters are estimated using a subspace-clustering algorithm which illustrates a weak bimodal distribution in the energy phase space of turbulence and sub-mesoscale motions for the very stable boundary layer. One mode indicates turbulence modulated by sub-mesoscale motions. Furthermore, intermittent turbulence appears if the sub-mesoscale intensity exceeds $$10 \%$$ 10 % of the mean kinetic energy in strong stratification.

Sign in / Sign up

Export Citation Format

Share Document