Dynamic coupling between carrier and dispersed phases in Rayleigh–Bénard convection laden with inertial isothermal particles

2021 ◽  
Vol 930 ◽  
Author(s):  
Wenwu Yang ◽  
Yi-Zhao Zhang ◽  
Bo-Fu Wang ◽  
Yuhong Dong ◽  
Quan Zhou
Keyword(s):  
Rayleigh Number ◽  
Kinetic Energy ◽  
Stokes Number ◽  
Point Particle ◽  
Inertial Particles ◽  
Dynamic Coupling ◽  
Exact Relation ◽  
Turbulent Momentum ◽  
Two Phases ◽  
Rayleigh Benard

We investigate the dynamic couplings between particles and fluid in turbulent Rayleigh–Bénard (RB) convection laden with isothermal inertial particles. Direct numerical simulations combined with the Lagrangian point-particle mode were carried out in the range of Rayleigh number $1\times 10^6 \le {Ra}\le 1 \times 10^8$ at Prandtl number ${Pr}=0.678$ for three Stokes numbers ${St_f}=1 \times 10^{-3}$ , $8 \times 10^{-3}$ and $2.5 \times 10^{-2}$ . It is found that the global heat transfer and the strength of turbulent momentum transfer are altered a small amount for the small Stokes number and large Stokes number as the coupling between the two phases is weak, whereas they are enhanced a large amount for the medium Stokes number due to strong coupling of the two phases. We then derived the exact relation of kinetic energy dissipation in the particle-laden RB convection to study the budget balance of induced and dissipated kinetic energy. The strength of the dynamic coupling can be clearly revealed from the percentage of particle-induced kinetic energy over the total induced kinetic energy. We further derived the power law relation of the averaged particles settling rate versus the Rayleigh number, i.e. $S_p/(d_p/H)^2{\sim} Ra^{1/2}$ , which is in remarkable agreement with our simulation. We found that the settling and preferential concentration of particles are strongly correlated with the coupling mechanisms.

scholarly journals Concentration and velocity statistics of inertial particles in upward and downward pipe flow

2017 ◽  
Vol 822 ◽  
pp. 640-663 ◽  
Author(s):  
J. L. G. Oliveira ◽  
C. W. M. van der Geld ◽  
J. G. M. Kuerten
Keyword(s):  
Liquid Velocity ◽  
Fluid Velocity ◽  
Three Dimensional ◽  
Stokes Number ◽  
Terminal Velocity ◽  
Point Particle ◽  
Inertial Particles ◽  
Carrier Fluid ◽  
Velocity Statistics ◽  
The Difference

Three-dimensional particle tracking velocimetry is applied to particle-laden turbulent pipe flows at a Reynolds number of 10 300, based on the bulk velocity and the pipe diameter, for developed fluid flow and not fully developed flow of inertial particles, which favours assessment of the radial migration of the inertial particles. Inertial particles with Stokes number ranging from 0.35 to 1.11, based on the particle relaxation time and the radial-dependent Kolmogorov time scale, and a ratio of the root-mean-square fluid velocity to the terminal velocity of order 1 have been used. Core peaking of the concentration of inertial particles in up-flow and wall peaking in down-flow have been found. The difference in mean particle and Eulerian mean liquid velocity is found to decrease to approximately zero near the wall in both flow directions. Although the carrier fluid has all of the characteristics of the corresponding turbulent single-phase flow, the Reynolds stress of the inertial particles is different near the wall in up-flow. These findings are explained from the preferential location of the inertial particles with the aid of direct numerical simulations with the point-particle approach.

Settling of heated particles in homogeneous turbulence

2016 ◽  
Vol 792 ◽  
pp. 869-893 ◽  
Author(s):  
Ari Frankel ◽  
H. Pouransari ◽  
F. Coletti ◽  
A. Mani
Keyword(s):  
Thermal Radiation ◽  
Settling Velocity ◽  
Stokes Number ◽  
Point Particle ◽  
Homogeneous Turbulence ◽  
Turbulence Structure ◽  
Inertial Particles ◽  
Order Unity ◽  
Buoyant Plumes ◽  
Heated Particle

We study the case of inertial particles heated by thermal radiation while settling by gravity through a turbulent transparent gas. We consider dilute and optically thin regimes in which each particle receives the same heat flux. Numerical simulations of forced homogeneous turbulence are performed taking into account the two-way coupling of both momentum and temperature between the dispersed and continuous phases. Particles much smaller than the smallest flow scales are considered and the point-particle approximation is adopted. The particle Stokes number (based on the Kolmogorov time scale) is of order unity, while the nominal settling velocity is up to an order of magnitude larger than the Kolmogorov velocity, marking a critical difference with previous two-way coupled simulations. It is found that non-heated particles enhance turbulence when their settling velocity is sufficiently high compared to the Kolmogorov velocity. Energy spectra show that the non-heated particle settling impacts both the very small and very large flow scales, while the intermediate scales are weakly affected. When heated, particles shed plumes of buoyant gas, further modifying the turbulence structure. At the considered radiation intensities, clustering is strong but the classic mechanism of preferential concentration is modified, while preferential sweeping is eliminated or even reversed. Particle heating also causes a significant reduction of the mean settling velocity, which is caused by rising buoyant plumes in the vicinity of particle clusters. The turbulent kinetic energy is affected non-monotonically as the radiation intensity is increased due to the competing effects of the downward gravitational force and the upward buoyancy force. The thermal radiation influences all scales of the turbulence. The effects of settling and buoyancy on the turbulence anisotropy are also discussed.

Turbulent Rotating Rayleigh–Benard Convection: Spatiotemporal and Statistical Study

2008 ◽  
Vol 131 (2) ◽  
Author(s):  
A. Husain ◽  
M. F. Baig ◽  
H. Varshney
Keyword(s):  
Rayleigh Number ◽  
Kinetic Energy ◽  
Rotation Rate ◽  
Large Scale ◽  
Coriolis Forces ◽  
Large Scale Structures ◽  
Vertical Heat ◽  
Scale Structures ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

The present study involves a 3D numerical investigation of rotating Rayleigh–Benard convection in a large aspect-ratio (8:8:1) rectangular enclosure. The rectangular cavity is rotated about a vertical axis passing through the center of the cavity. The governing equations of mass, momentum, and energy for a frame rotating with the enclosure, subject to generalized Boussinesq approximation applied to the body and centrifugal force terms, have been solved on a collocated grid using a semi-implicit finite difference technique. The simulations have been carried out for liquid metal flows having a fixed Prandtl number Pr=0.01 and fixed Rayleigh number Ra=107 while rotational Rayleigh number Raw and Taylor number Ta are varied through nondimensional rotation rate (Ω) ranging from 0 to 104. Generation of large-scale structures is observed at low-rotation (Ω=10) rates though at higher-rotation rates (Ω=104) the increase in magnitude of Coriolis forces leads to redistribution of buoyancy-induced vertical kinetic energy to horizontal kinetic energy. This brings about inhibition of vertical fluid transport, thereby leading to reduced vertical heat transfer. The magnitude of rms velocities remains unaffected with an increase in Coriolis forces from Ω=0 to 104. An increase in rotational buoyancy (Raw), at constant rotation rate (Ω=104), on variation in Raw/Ta from 10−3 to 10−2 results in enhanced breakup of large-scale structures with a consequent decrease in rms velocities but with negligible reduction in vertical heat transport.

scholarly journals Two-way coupled stochastic model for dispersion of inertial particles in turbulence

2012 ◽  
Vol 700 ◽  
pp. 29-62 ◽  
Author(s):  
Madhusudan G. Pai ◽  
Shankar Subramaniam
Keyword(s):  
Kinetic Energy ◽  
Stochastic Model ◽  
Turbulent Kinetic Energy ◽  
Time Scales ◽  
Time Scale ◽  
Stokes Number ◽  
Particle Dispersion ◽  
Point Particle ◽  
Two Phase Flows ◽  
Two Phase

AbstractTurbulent two-phase flows are characterized by the presence of multiple time and length scales. Of particular interest in flows with non-negligible interphase momentum coupling are the time scales associated with interphase turbulent kinetic energy transfer (TKE) and inertial particle dispersion. Point-particle direct numerical simulations (DNS) of homogeneous turbulent flows laden with sub-Kolmogorov size particles report that the time scale associated with the interphase TKE transfer behaves differently with Stokes number than the time scale associated with particle dispersion. Here, the Stokes number is defined as the ratio of the particle momentum response time scale to the Kolmogorov time scale of turbulence. In this study, we propose a two-way coupled stochastic model (CSM), which is a system of two coupled Langevin equations for the fluctuating velocities in each phase. The basis for the model is the Eulerian–Eulerian probability density function formalism for two-phase flows that was established in Pai & Subramaniam (J. Fluid Mech., vol. 628, 2009, pp. 181–228). This new model possesses the unique capability ofsimultaneouslycapturing the disparate dependence of the time scales associated with interphase TKE transfer and particle dispersion on Stokes number. This is ascertained by comparing predicted trends of statistics of turbulent kinetic energy and particle dispersion in both phases from CSM, for varying Stokes number and mass loading, with point-particle DNS datasets of homogeneous particle-laden flows.

Modulation of turbulence by dispersed solid particles in a spatially developing flat-plate boundary layer

2016 ◽  
Vol 802 ◽  
pp. 359-394 ◽  
Author(s):  
Dong Li ◽  
Kun Luo ◽  
Jianren Fan
Keyword(s):  
Boundary Layer ◽  
Flat Plate ◽  
Fluid Velocity ◽  
Plate Boundary ◽  
Stokes Number ◽  
Point Particle ◽  
Solid Particles ◽  
Inertial Particles ◽  
Mass Loading ◽  
Turbulence Modulation

Direct numerical simulations of particle-laden spatially developing turbulent boundary layers over a flat plate have been performed to investigate the effect of inertial particles on turbulence modulation, using the Eulerian–Lagrangian point-particle approach with two-way coupling. The particles are smaller than the Kolmogorov length scale of the dilute flow, and inter-particle collisions are not considered. The simulation results show that the addition of small solid particles increases the mean streamwise fluid velocity, which in turn leads to a reduction in the boundary layer integral parameters and an increase in the skin-friction drag. These effects become more pronounced as the particle Stokes number and mass loading increase. The streamwise turbulence intensity is slightly enhanced in the close vicinity of the wall but damped in the outer layer. In contrast, the Reynolds stress and the turbulence intensities in the wall-normal and spanwise directions are substantially attenuated across the entire boundary layer, and the levels of attenuation increase monotonically with both particle Stokes number and mass loading. The exchange of kinetic energy between particles and fluid indicates that particle–fluid interactions cause extra energy dissipation, which plays a crucial role in turbulence modulation.

Interpretation of shadowgraph patterns in Rayleigh-Bénard convection

1988 ◽  
Vol 190 ◽  
pp. 451-469 ◽  
Author(s):  
D. R. Jenkins
Keyword(s):  
Rayleigh Number ◽  
Linear Theory ◽  
Vertical Component ◽  
Fluid Flows ◽  
Bénard Convection ◽  
Benard Convection ◽  
Higher Order Terms ◽  
The Relationship ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

The relationship between observations of cellular Rayleigh-Bénard convection using shadowgraphs and theoretical expressions for convection planforms is considered. We determine the shadowgraphs that ought to be observed if the convection is as given by theoretical expressions for roll, square or hexagonal planforms and compare them with actual experiments. Expressions for the planforms derived from linear theory, valid for low supercritical Rayleigh number, produce unambiguous shadowgraphs consisting of cells bounded by bright lines, which correspond to surfaces through which no fluid flows and on which the vertical component of velocity is directed downwards. Dark spots at the centre of cells, indicating regions of hot, rising fluid, are not accounted for by linear theory, but can be produced by adding higher-order terms, predominantly due to the temperature dependence of a material property of the fluid, such as its viscosity.

Snowflakes in the atmospheric surface layer: observation of particle–turbulence dynamics

2017 ◽  
Vol 814 ◽  
pp. 592-613 ◽  
Author(s):  
Andras Nemes ◽  
Teja Dasari ◽  
Jiarong Hong ◽  
Michele Guala ◽  
Filippo Coletti
Keyword(s):  
Large Scale ◽  
Isotropic Turbulence ◽  
Volume Fraction ◽  
Response Times ◽  
Field Measurements ◽  
Stokes Number ◽  
Inertial Particles ◽  
Natural Setting ◽  
Scale Particle ◽  
Particle Imaging

We report on optical field measurements of snow settling in atmospheric turbulence at $Re_{\unicode[STIX]{x1D706}}=940$. It is found that the snowflakes exhibit hallmark features of inertial particles in turbulence. The snow motion is analysed in both Eulerian and Lagrangian frameworks by large-scale particle imaging, while sonic anemometry is used to characterize the flow field. Additionally, the snowflake size and morphology are assessed by digital in-line holography. The low volume fraction and mass loading imply a one-way interaction with the turbulent air. Acceleration probability density functions show wide exponential tails consistent with laboratory and numerical studies of homogeneous isotropic turbulence. Invoking the assumption that the particle acceleration has a stronger dependence on the Stokes number than on the specific features of the turbulence (e.g. precise Reynolds number and large-scale anisotropy), we make inferences on the snowflakes’ aerodynamic response time. In particular, we observe that their acceleration distribution is consistent with that of particles of Stokes number in the range $St=0.1{-}0.4$ based on the Kolmogorov time scale. The still-air terminal velocities estimated for the resulting range of aerodynamic response times are significantly smaller than the measured snow particle fall speed. This is interpreted as a manifestation of settling enhancement by turbulence, which is observed here for the first time in a natural setting.

Plume formation and resonant bifurcations in porous-media convection

1994 ◽  
Vol 272 ◽  
pp. 67-90 ◽  
Author(s):  
Michael D. Graham ◽  
Paul H. Steen
Keyword(s):  
Boundary Layer ◽  
Porous Media ◽  
Rayleigh Number ◽  
Scaling Laws ◽  
Thermal Plumes ◽  
Scaling Behaviour ◽  
Parametric Resonances ◽  
Classical Boundary ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

The classical boundary-layer scaling laws proposed by Howard for Rayleigh–Bénard convection at high Rayleigh number extend to the analogous case of convection in saturated porous media. We computationally study two-dimensional porous-media convection near the onset of this scaling behaviour. The main result of the paper is the observation and study of instabilities that lead to deviations from the scaling relations.At Rayleigh numbers below the scaling regime, boundary-layer fluctuations born at a Hopf bifurcation strengthen and eventually develop into thermal plumes. The appearance of plumes corresponds to the onset of the boundary-layer scaling behaviour of the oscillation frequency and mean Nusselt number, in agreement with the classical theory. As the Rayleigh number increases further, the flow undergoes instabilities that lead to ‘bubbles’ in parameter space of quasi-periodic flow, and eventually to weakly chaotic flow. The instabilities disturb the plume formation process, effectively leading to a phase modulation of the process and to deviations from the scaling laws. We argue that these instabilities correspond to parametric resonances between the timescale for plume formation and the characteristic convection timescale of the flow.

scholarly journals A RANS K-? SIMULATION OF 2D TURBULENT NATURAL CONVECTION IN AN ENCLOSURE WITH HEATING SOURCES

2019 ◽  
Vol 20 (1) ◽  
pp. 229-244
Author(s):  
Mehdi Ahmadi ◽  
Seyed Ali Agha Mirjalily ◽  
Seyed Amir Abbas Oloomi
Keyword(s):  
Heat Transfer ◽  
Natural Convection ◽  
Rayleigh Number ◽  
Kinetic Energy ◽  
Aspect Ratio ◽  
Heat Transfer Rate ◽  
Transfer Rate ◽  
Turbulent Natural Convection ◽  
Cold Source ◽  
Thermal Sources

ABSTRACT: This study is conducted to investigate turbulent natural convection flow in an enclosure with thermal sources using the low-Reynolds number (LRN) k-? model. This enclosure has a cold source with temperature Tc and a hot source with temperature Th as thermal sources, other walls of the enclosure are adiabatic. The aim of this study is to predict the effect of change in Rayleigh number, repositioning of cold and hot sources, and thermal sources aspect ratio on the flow field, temperature, and rate of heat transfer. To achieve this aim, the equations of continuity, momentum, energy, turbulent kinetic energy, and kinetic energy dissipation are employed in the case of 2D turbulence with constant thermo-physical properties except the density in the buoyancy term (Boussinesq approximation). To numerically solve these equations, the finite volume method and SIMPLE algorithm are used. According to the modeling results, the most optimal temperature distribution in the enclosure is seen when the hot source is below the cold source. With decreasing distance between hot and cold sources, heat transfer rate increases. The maximal heat transfer rate is derived via study of the heating sources aspect ratio. In constant positions of cold and hot sources on a wall, the heat transfer rate increases with increasing Rayleigh number (Ra=109-1011). ABSTAK: Kajian ini dijalankan bagi mengkaji perubahan semula jadi aliran perolakan dalam tempat tertutup dengan sumber haba menggunakan model k-? nombor Reynolds-rendah (LRN). Bekas tertutup ini mempunyai dua sumber haba iaitu sumber sejuk dengan suhu Tc dan sumber panas dengan suhu Th, manakala dinding lain bekas ini adalah adiabatik. Tujuan kajian ini adalah bagi mengesan perubahan nombor Rayleigh, mengubah sumber sejuk dan panas dan nisbah sumber haba kepada kawasan aliran, suhu dan halaju perubahan haba. Bagi mencapai tujuan tersebut, persamaan sambungan, momentum, tenaga, tenaga kinetik perolakan, dan pengurangan tenaga kinetik telah dilaksanakan dalam kes perolakan 2D dengan sifat fizikal-haba berterusan (malar) kecuali isipadu terma keapungan (anggaran Boussinesq). Bagi menyelesaikan persamaan ini secara berangka, kaedah isipadu terhad dan algorithma MUDAH telah digunakan. Berdasarkan keputusan model, suhu distribusi optimal dalam bekas tertutup dilihat apabila sumber panas adalah kurang daripada sumber sejuk. Dengan pengurangan jarak antara sumber panas dan sejuk, kadar pertukaran haba meningkat. Kadar pertukaran haba maksima telah diperoleh melalui kajian nisbah  aspek sumber pemanasan. Kadar pertukaran haba bertambah dengan bertambahnya nombor Rayleigh  (Ra=109-1011), pada posisi tetap sumber sejuk dan panas pada dinding bekas.

Sign in / Sign up

Export Citation Format

Share Document