Application of a Volume Averaged k-ε Model to Particle-Laden Turbulent Channel Flow

2009 ◽  
Vol 131 (10) ◽  
Author(s):  
J. D. Schwarzkopf ◽  
C. T. Crowe ◽  
P. Dutta
Keyword(s):  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
Carrier Phase ◽  
Critical Role ◽  
Turbulent Channel Flow ◽  
Mean Velocity ◽  
Closed Set ◽  
Production Term ◽  
Velocity Gradients ◽  
Additional Production

A closed set of volume averaged equations for modeling turbulence in the carrier phase of particle-laden flows is presented. The equations incorporate a recently developed dissipation transport equation that contains an additional production of dissipation term due to particle surfaces. 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. To assess the effects of this additional production term, a numerical model was developed and applied to particles falling in a channel of downward turbulent air flow. Boundary conditions were 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 coefficients associated with the production of dissipation due to mean velocity gradients and particle surfaces were varied to assess the effects of the dispersed phase on the carrier phase turbulent kinetic energy across the channel. The results show that the model predicts a decrease in turbulent kinetic energy near the wall with increased particle loading, and that the dissipation coefficients play a critical role in predicting the turbulent kinetic energy in particle-laden turbulent flows.

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.

Interaction between a spatially growing turbulent boundary layer and embedded streamwise vortices

1996 ◽  
Vol 326 ◽  
pp. 151-179 ◽  
Author(s):  
Junhui Liu ◽  
Ugo Piomelli ◽  
Philippe R. Spalart
Keyword(s):  
Boundary Layer ◽  
Kinetic Energy ◽  
Turbulent Boundary Layer ◽  
Turbulent Kinetic Energy ◽  
Vortex Core ◽  
Eddy Viscosity ◽  
Shear Stresses ◽  
Streamwise Vortices ◽  
Mean Velocity ◽  
Production Term

The interaction between a zero-pressure-gradient turbulent boundary layer and a pair of strong, common-flow-down, streamwise vortices with a sizeable velocity deficit is studied by large-eddy simulation. The subgrid-scale stresses are modelled by a localized dynamic eddy-viscosity model. The results agree well with experimental data. The vortices drastically distort the boundary layer, and produce large spanwise variations of the skin friction. The Reynolds stresses are highly three-dimensional. High levels of kinetic energy are found both in the upwash region and in the vortex core. The two secondary shear stresses are significant in the vortex region, with magnitudes comparable to the primary one. Turbulent transport from the immediate upwash region is partly responsible for the high levels of turbulent kinetic energy in the vortex core; its effect on the primary stress 〈u′v′〉 is less significant. The mean velocity gradients play an important role in the generation of 〈u′v′〉 in all regions, while they are negligible in the generation of turbulent kinetic energy in the vortex core. The pressure-strain correlations are generally of opposite sign to the production terms except in the vortex core, where they have the same sign as the production term in the budget of 〈u′v′〉. The results highlight the limitations of the eddy-viscosity assumption (in a Reynolds-averaged context) for flows of this type, as well as the excessive diffusion predicted by typical turbulence models.

NUMERICAL RESEARCH ON THE EFFECT OF FIBERS ON THE PROPERTY OF TURBULENT FIBER SUSPENSION IN A CHANNEL

2009 ◽  
Vol 23 (03) ◽  
pp. 509-512 ◽  
Author(s):  
SUHUA SHEN ◽  
JIANZHONG LIN
Keyword(s):  
Reynolds Number ◽  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
Channel Flow ◽  
Turbulent Channel Flow ◽  
Fiber Suspension ◽  
Additional Stress ◽  
Mean Velocity ◽  
Navier Stokes Equation ◽  
Fiber Concentration

To explore the rheological property in turbulent channel flow of fiber suspensions, the equation of probability distribution function for mean fiber orientation and the Reynolds averaged Navier-Stokes equation with the term of additional stress resulted from fibers were solved with numerical methods to get the distributions of the mean velocity and turbulent kinetic energy. The simulation results show that the effect of fibers on turbulent channel flow is equivalent to an additional viscosity. The turbulent velocity profiles of fiber suspension become gradually sharper by increasing the fiber concentration and/or decreasing the Reynolds number. The turbulent kinetic energy will increase with increasing Reynolds number and fiber concentration.

LDV Measurements of Periodic Fully Developed Main and Secondary Flows in a Channel With Rib-Disturbed Walls

1993 ◽  
Vol 115 (1) ◽  
pp. 109-114 ◽  
Author(s):  
T.-M. Liou ◽  
Y.-Y. Wu ◽  
Y. Chang
Keyword(s):  
Kinetic Energy ◽  
Secondary Flow ◽  
Turbulent Kinetic Energy ◽  
Secondary Flows ◽  
Laser Doppler Velocimeter ◽  
Parameter Distribution ◽  
Reynolds Stresses ◽  
Mean Velocity ◽  
Production Term ◽  
Basic Assumptions

Laser-Doppler velocimeter measurements of mean velocities, turbulence intensities, and Reynolds stresses are presented for periodic fully developed flows in a channel with square rib-disturbed walls on two opposite sides. Quantities such as the vorticity thickness and turbulent kinetic energy are used to characterize the flow. The investigated flow was periodic in space. The Reynolds number based on the channel hydraulic diameter was 3.3×104. The ratios of pitch to rib-height and rib-height to chamber-height were 10 and 0.133, respectively. Regions where maximum and minimum Reynolds stress and turbulent kinetic energy occurred were identified from the results. The growth rate of the shear layers of the present study was compared with that of a backward-facing step. The measured turbulence anisotropy and structure parameter distribution were used to examine the basic assumptions embedded in the k–ε and k–ε–A models. For a given axial station, the peak axial mean-velocity was found not to occur at the center point. The secondary flow was determined to be Prandtl’s secondary flow of the second kind according to the measured streamwise mean vorticity and its production term.

Rotating Effect on Fluid Flow in a Smooth Duct With a 180-Deg Sharp Turn

2000 ◽  
Author(s):  
Chung-Chu Chen ◽  
Tong-Miin Liou
Keyword(s):  
Heat Transfer ◽  
Kinetic Energy ◽  
Secondary Flow ◽  
Turbulent Kinetic Energy ◽  
Turbine Blades ◽  
Outer Part ◽  
Intensity Level ◽  
Mean Velocity ◽  
High Heat Transfer ◽  
Sharp Turn

Laser-Doppler velocimetry (LDV) measurements are presented of turbulent flow in a two-pass square-sectioned duct simulating the coolant passages employed in gas turbine blades under rotating and non-rotating conditions. For all cases studied, the Reynolds number characterized by duct hydraulic diameter (Dh) and bulk mean velocity (Ub) was fixed at 1 × 104. The rotating case had a range of rotation number (Ro = ΩDh/Ub) from 0 to 0.2. It is found that both the skewness of streamwise mean velocity and magnitude of secondary-flow velocity increase linearly, and the magnitude of turbulence intensity level increases non-linearly with increasing Ro. As Ro is increased, the curvature induced symmetric Dean vortices in the turn for Ro = 0 is gradually dominated by a single vortex most of which impinges directly on the outer part of leading wall. The high turbulent kinetic energy is closely related to the dominant vortex prevailing inside the 180-deg sharp turn. For the first time, the measured flow characteristics account for the reported spanwise heat transfer distributions in the rotating channels, especially the high heat transfer enhancement on the leading wall in the turn. For both rotating and non-rotating cases, the direction and strength of the secondary flow with respect to the wall are the most important fluid dynamic factors affecting local heat transfer distributions inside a 180-deg sharp turn. The role of the turbulent kinetic energy in affecting the overall enhancement of heat transfer is well addressed.

scholarly journals Discussion on improved method of turbulence model for supercritical water flow and heat transfer

2020 ◽  
Vol 24 (5 Part A) ◽  
pp. 2729-2741
Author(s):  
Zhenchuan Wang ◽  
Guoli Qi ◽  
Meijun Li
Keyword(s):  
Heat Transfer ◽  
Reynolds Number ◽  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
Turbulence Model ◽  
Buoyancy Effect ◽  
Improved Method ◽  
Production Term ◽  
Flow And Heat Transfer ◽  
Effect Model

The turbulence model fails in supercritical fluid-flow and heat transfer simulation, owing to the drastic change of thermal properties. The inappropriate buoyancy effect model and the improper turbulent Prandtl number model are several of these factors lead to the original low-Reynolds number turbulence model unable to predict the wall temperature for vertically heated tubes under the deteriorate heat transfer conditions. This paper proposed a simplified improved method to modify the turbulence model, using the generalized gradient diffusion hypothesis approximation model for the production term of the turbulent kinetic energy due to the buoyancy effect, using a turbulence Prandtl number model for the turbulent thermal diffusivity instead of the constant number. A better agreement was accomplished by the improved turbulence model compared with the experimental data. The main reason for the over-predicted wall temperature by the original turbulence model is the misuse of the buoyancy effect model. In the improved model, the production term of the turbulent kinetic energy is much higher than the results calculated by the original turbulence model, especially in the boundary-layer. A more accurate model for the production term of the turbulent kinetic energy is the main direction of further modification for the low Reynolds number turbulence model.

Coherent Velocity Structures in the Mixed Layer: Characteristics, Energetics, and Turbulent Kinetic Energy Budget

2021 ◽  
Author(s):  
Ewa Jarosz ◽  
Hemantha W. Wijesekera ◽  
David W. Wang
Keyword(s):  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
Mixed Layer ◽  
Vertical Velocity ◽  
Mixed Layer Depth ◽  
Vertical Transport ◽  
Rate Of Change ◽  
Production Term ◽  
Forcing Mechanisms ◽  
Surface Buoyancy

AbstractVelocity, hydrographic, and microstructure observations collected under moderate to high winds, large surface waves, and significant ocean-surface heat losses were utilized to examine coherent velocity structures (CVS) and turbulent kinetic energy (TKE) budget in the mixed layer on the outer shelf in the northern Gulf of Mexico in February 2017. The CVS exhibited larger downward velocities in downweling regions and weaker upward velocities in broader upwelling regions, elevated vertical velocity variance, vertical velocity maxima in the upper part of the mixed layer, and phasing of crosswind velocities relative to vertical velocities near the base of the mixed layer. Temporal scales ranged from 10 min to 40 min and estimated lateral scales ranged from 90 m to 430 m, which were 1.5 – 6 times larger than the mixed layer depth. Nondimensional parameters, Langmuir and Hoenikker numbers, indicated that plausible forcing mechanisms were surface-wave driven Langmuir vortex and destabilizing surface buoyancy flux. The rate of change of TKE, shear production, Stokes production, buoyancy production, vertical transport of TKE, and dissipation in the TKE budget were evaluated. The shear and Stokes productions, dissipation, and vertical transport of TKE were the dominant terms. The buoyancy production term was important at the sea surface, but it decreased rapidly in the interior. A large imbalance term was found under the unstable, high wind, and high-sea state conditions. The cause of this imbalance cannot be determined with certainty through analyses of the available observations; however, our speculation is that the pressure transport is significant under these conditions.

Rotating Effect on Fluid Flow in Two Smooth Ducts Connected by a 180-Degree Bend

2003 ◽  
Vol 125 (1) ◽  
pp. 138-148 ◽  
Author(s):  
Tong-Miin Liou ◽  
Chung-Chu Chen ◽  
Meng-Yu Chen
Keyword(s):  
Heat Transfer ◽  
Kinetic Energy ◽  
Secondary Flow ◽  
Turbulent Kinetic Energy ◽  
Separation Bubble ◽  
Turbine Blades ◽  
Outer Part ◽  
Intensity Level ◽  
Mean Velocity ◽  
Sharp Turn

Laser Doppler velocimetry (LDV) measurements are presented of turbulent flow in a two-pass square-sectioned smooth duct simulating the coolant passages employed in gas turbine blades under rotating and nonrotating conditions. For all cases studied, the Reynolds number characterized by duct hydraulic diameter and bulk mean velocity was fixed at 1×104. The rotation number Ro was varied from 0 to 0.2. It is found that as Ro is increased, both the skewness (SK) of streamwise mean velocity and magnitude of secondary-flow velocity increase linearly, SK=2.3 Ro and U2+V2¯/Uh=2.3 Ro+0.4, and the magnitude of turbulence intensity level increases exponentially. As Ro is increased, the curvature induced symmetric Dean vortices in the turn for Ro=0 is gradually dominated by a single vortex most of which impinges directly on the outer part of leading wall. The high turbulent kinetic energy is closely related to the dominant vortex prevailing inside the 180-deg sharp turn. The size of separation bubble immediately after the turn is found to diminish to null as Ro is increased from 0 to 0.2. A simple correlation is developed between the bubble size and Ro. A critical range of Ro responsible for the switch of faster moving flow from near the outer wall to the inner wall is identified. For both rotating and nonrotating cases, the direction and strength of the secondary flow with respect to the wall are the most important fluid dynamic factors affecting local the heat transfer distributions inside a 180-deg sharp turn. The role of the turbulent kinetic energy in the overall enhancement of heat transfer is well addressed.

The evolution of a uniformly sheared thermally stratified turbulent flow

1997 ◽  
Vol 334 ◽  
pp. 61-86 ◽  
Author(s):  
PAUL PICCIRILLO ◽  
CHARLES W. VAN ATTA
Keyword(s):  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
Richardson Number ◽  
Initial Conditions ◽  
Velocity Shear ◽  
Spectral Energy ◽  
Small Scale ◽  
Mean Flow ◽  
Mean Velocity ◽  
Small Scales

Experiments were carried out in a new type of stratified flow facility to study the evolution of turbulence in a mean flow possessing both uniform stable stratification and uniform mean shear.The new facility is a thermally stratified wind tunnel consisting of ten independent supply layers, each with its own blower and heaters, and is capable of producing arbitrary temperature and velocity profiles in the test section. In the experiments, four different sized turbulence-generating grids were used to study the effect of different initial conditions. All three components of the velocity were measured, along with the temperature. Root-mean-square quantities and correlations were measured, along with their corresponding power and cross-spectra.As the gradient Richardson number Ri = N2/(dU/dz)2 was increased, the downstream spatial evolution of the turbulent kinetic energy changed from increasing, to stationary, to decreasing. The stationary value of the Richardson number, Ricr, was found to be an increasing function of the dimensionless shear parameter Sq2/∈ (where S = dU/dz is the mean velocity shear, q2 is the turbulent kinetic energy, and ∈ is the viscous dissipation).The turbulence was found to be highly anisotropic, both at the small scales and at the large scales, and anisotropy was found to increase with increasing Ri. The evolution of the velocity power spectra for Ri [les ] Ricr, in which the energy of the large scales increases while the energy in the small scales decreases, suggests that the small-scale anisotropy is caused, or at least amplified, by buoyancy forces which reduce the amount of spectral energy transfer from large to small scales. For the largest values of Ri, countergradient buoyancy flux occurred for the small scales of the turbulence, an effect noted earlier in the numerical results of Holt et al. (1992), Gerz et al. (1989), and Gerz & Schumann (1991).

Simulation of a propelled wake with moderate excess momentum in a stratified fluid

2011 ◽  
Vol 692 ◽  
pp. 28-52 ◽  
Author(s):  
Matthew B. de Stadler ◽  
Sutanu Sarkar
Keyword(s):  
Kinetic Energy ◽  
Velocity Profile ◽  
Turbulent Kinetic Energy ◽  
Vertical Direction ◽  
Qualitative Change ◽  
Stratified Fluid ◽  
Mean Velocity ◽  
Moderate Amount ◽  
Excess Momentum ◽  
The Mean

AbstractDirect numerical simulation is used to simulate the turbulent wake behind an accelerating axisymmetric self-propelled body in a stratified fluid. Acceleration is modelled by adding a velocity profile corresponding to net thrust to a self-propelled velocity profile resulting in a wake with excess momentum. The effect of a small to moderate amount of excess momentum on the initially momentumless self-propelled wake is investigated to evaluate if the addition of excess momentum leads to a large qualitative change in wake dynamics. Both the amount and shape of excess momentum are varied. Increasing the amount of excess momentum and/or decreasing the radial extent of excess momentum was found to increase the defect velocity, mean kinetic energy, shear in the velocity gradient and the wake width. The increased shear in the mean profile resulted in increased production of turbulent kinetic energy leading to an increase in turbulent kinetic energy and its dissipation. Slightly larger vorticity structures were observed in the late wake with excess momentum although the differences between vorticity structures in the self-propelled and 40 % excess momentum cases was significantly smaller than suggested by previous experiments. Buoyancy was found to preserve the doubly inflected velocity profile in the vertical direction, and similarity for the mean velocity and turbulent kinetic energy was found to occur in both horizontal and vertical directions. While quantitative differences were observed between cases with and without excess momentum, qualitatively similar evolution was found to occur.

Sign in / Sign up

Export Citation Format

Share Document