Heat transfer mechanisms in bubbly Rayleigh-Bénard convection

Analysis of the heat transfer fluctuations in the Rayleigh-Bénard convection of concentrated emulsions with finite-size droplets

Employing numerical simulations, we provide an accurate insight into the heat transfer mechanism in the Rayleigh-Bénard convection of concentrated emulsions with finite-size droplets. We focus on the unsteady dynamics characterizing the thermal convection of these complex fluids close to the transition from conductive to convective states, where the heat transfer phenomenon, expressed in terms of the Nusselt number Nu, is characterized by pronounced fluctuations triggered by collective droplet motion [F. Pelusi et al., Soft Matter, 2021, 17(13), 3709-3721]. By systematically increasing the droplet concentration, we show how these fluctuations emerge along with the segregation of "extreme events" in the boundary layers, causing intermittent bursts in the heat flux fluctuations. Furthermore, we quantify the extension S and the duration of the coherent droplet motion accompanying these extreme events via a suitable statistical analysis involving the droplet displacements. We show how the increase in droplet concentration results in a power-law behaviour of the probability distribution function of S and and how this outcome is robust at changing the analysis protocol. Our work offers a comprehensive picture, linking macroscopic heat transfer fluctuations with the statistics of droplets at the mesoscale.

Rayleigh-Bénard convection of a model emulsion: anomalous heat-flux fluctuations and finite-size droplet effects

We present mesoscale numerical simulations of Rayleigh-Bénard (RB) convection in a two-dimensional model emulsion. The systems under study are constituted of finite-size droplets, whose concentration Φ₀ is systematically varied from small (Newtonian emulsions) to large values (non-Newtonian emulsions). We focus on the characterisation of the heat transfer properties close to the transition from conductive to convective states, where it is well known that a homogeneous Newtonian system exhibits a steady flow and a time-independent heat flux. In marked contrast, emulsions exhibit non-steady dynamics with fluctuations in the heat flux. In this paper, we aim at the characterisation of such non-steady dynamics via detailed studies on the time-averaged heat flux and its fluctuations. To quantitatively understand the time-averaged heat flux, we propose a side-by-side comparison between the emulsion system and a single-phase (SP) system, whose viscosity is suitably constructed from the shear rheology of the emulsion. We show that such local closure works well only when a suitable degree of coarse-graining (at the droplet scale) is introduced in the local viscosity. To delve deeper into the fluctuations in the heat flux, we furthermore propose a side-by-side comparison between a Newtonian emulsion (i.e., with a small droplet concentration) and a non-Newtonian emulsion (i.e., with a large droplet concentration), at fixed time-averaged heat flux. This comparison elucidates that finite-size droplets and the non-Newtonian rheology cooperate to trigger enhanced heat-flux fluctuations at the droplet scales. These enhanced fluctuations are rooted in the emergence of space correlations among distant droplets, which we highlight via direct measurements of the droplets displacement and the characterisation of the associated correlation function. The observed findings offer insights on heat transfer properties for confined systems possessing finite-size constituents.

Linear stability analysis of bubble-induced convection in a horizontal liquid layer

We investigate with a linear analysis the stability of a horizontal liquid layer subjected to injection of gas bubbles through a bottom wall. The injection is assumed uniform in space and constant in time. Injected bubbles ascend in the liquid layer due to the Archimedean buoyancy force and are ejected from the top free surface of the liquid layer. Modeling this two-phase flow system as two interpenetrating liquid and gas continua, we show that homogeneous upward gas flows become unstable at large gas fluxes. We determine the critical conditions of this homogeneous-heterogeneous regime transition and show that the critical modes are made of stationary convection rolls, either multi- or whole-layered depending on liquid viscosity, the radius of bubbles, and the thickness of liquid layer. By examining the energy transfer from base to perturbation flows, we indicate that liquid convective motion is driven by the buoyancy on heterogeneously distributed bubbles. We also reveal that the lift forces on bubbles have significant stabilizing effects by homogenizing bubble distribution close to the bottom wall.

Statistical properties of thermally expandable particles in soft-turbulence Rayleigh-Bénard convection

The dynamics of inertial particles in Rayleigh-Bénard convection, where both particles and fluid exhibit thermal expansion, is studied using direct numerical simulations (DNS) in the soft-turbulence regime. We consider the effect of particles with a thermal expansion coefficient larger than that of the fluid, causing particles to become lighter than the fluid near the hot bottom plate and heavier than the fluid near the cold top plate. Because of the opposite directions of the net Archimedes' force on particles and fluid, particles deposited at the plate now experience a relative force towards the bulk. The characteristic time for this motion towards the bulk to happen, quantified as the time particles spend inside the thermal boundary layers (BLs) at the plates, is shown to depend on the thermal response time, [Formula: see text], and the thermal expansion coefficient of particles relative to that of the fluid, [Formula: see text]. In particular, the residence time is constant for small thermal response times, [Formula: see text], and increasing with [Formula: see text] for larger thermal response times, [Formula: see text]. Also, the thermal BL residence time is increasing with decreasing K. A one-dimensional (1D) model is developed, where particles experience thermal inertia and their motion is purely dependent on the buoyancy force. Although the values do not match one-to-one, this highly simplified 1D model does predict a regime of a constant thermal BL residence time for smaller thermal response times and a regime of increasing residence time with [Formula: see text] for larger response times, thus explaining the trends in the DNS data well.

Clusterlike instabilities in bubble-plume-driven flows

Continuous release of gas bubbles in large numbers from a localized source in a liquid column, popularly known as "bubble plumes", is very relevant in nature and industries. The bubble plumes morphologically consist of a long continuous stem supporting a dispersed head. Through our direct numerical simulations using two-way coupled Euler-Lagrangian framework, we show that a bubble plume rising in a quiescent liquid column develops clusterlike instabilities for the Grashof numbers, Gr>145. For levels Gr<100, the stem is continuous with a small plume head, whereas at high buoyancy (Gr>350), the plume stem shows intermittently passing puffing instabilities in the form of bubble clusters. The clusters are a group of bubbles localized in space with high concentration that travel upward with speed C_{ph}=0.45U_{C} and are separated by a distance of at least 5L_{0}, where U_{C} is the characteristic velocity and L_{0} is the characteristic length based on the injection conditions. The bubble rise Reynolds numbers in the steady state for both the plume head and the stem shows Re∝Gr^{0.45±0.03}, and the proportionality constant is ten times higher in the plume stem than in the plume head. In the plume core, the spatial acceleration due to the bubble motion generates the turbulent production, whereas, at the plume edge, the small-scale fluctuations generate the mean vorticity. At high Gr, the clusters evolve due to the lift forces acting on the bubbles as a result of increase in the mean vorticity. While rising, bubbles entrain the liquid from the surroundings, and we found that the entrainment rate is not as strong as compared to the classical thermal plumes.

Effects of particle settling on Rayleigh-Bénard convection

The effect of particles falling under gravity in a weakly turbulent Rayleigh-Bénard gas flow is studied numerically. The particle Stokes number is varied between 0.01 and 1 and their temperature is held fixed at the temperature of the cold plate, of the hot plate, or the mean between these values. Mechanical, thermal, and combined mechanical and thermal couplings between the particles and the fluid are studied separately. It is shown that the mechanical coupling plays a greater and greater role in the increase of the Nusselt number with increasing particle size. A rather unexpected result is an unusual kind of reverse one-way coupling, in the sense that the fluid is found to be strongly influenced by the particles, while the particles themselves appear to be little affected by the fluid, despite the relative smallness of the Stokes numbers. It is shown that this result derives from the very strong constraint on the fluid behavior imposed by the vanishing of the mean fluid vertical velocity over the cross sections of the cell demanded by continuity.

Heat transport in bubbling turbulent convection

Boiling is an extremely effective way to promote heat transfer from a hot surface to a liquid due to numerous mechanisms, many of which are not understood in quantitative detail. An important component of the overall process is that the buoyancy of the bubble compounds with that of the liquid to give rise to a much-enhanced natural convection. In this article, we focus specifically on this enhancement and present a numerical study of the resulting two-phase Rayleigh-Bénard convection process in a cylindrical cell with a diameter equal to its height. We make no attempt to model other aspects of the boiling process such as bubble nucleation and detachment. The cell base and top are held at temperatures above and below the boiling point of the liquid, respectively. By keeping this difference constant, we study the effect of the liquid superheat in a Rayleigh number range that, in the absence of boiling, would be between 2 × 10(6) and 5 × 10(9). We find a considerable enhancement of the heat transfer and study its dependence on the number of bubbles, the degree of superheat of the hot cell bottom, and the Rayleigh number. The increased buoyancy provided by the bubbles leads to more energetic hot plumes detaching from the cell bottom, and the strength of the circulation in the cell is significantly increased. Our results are in general agreement with recent experiments on boiling Rayleigh-Bénard convection.

New perspectives in turbulent Rayleigh-Bénard convection

Recent experimental, numerical and theoretical advances in turbulent Rayleigh-Bénard convection are presented. Particular emphasis is given to the physics and structure of the thermal and velocity boundary layers which play a key role for the better understanding of the turbulent transport of heat and momentum in convection at high and very high Rayleigh numbers. We also discuss important extensions of Rayleigh-Bénard convection such as non-Oberbeck-Boussinesq effects and convection with phase changes.

Convection in multiphase fluid flows using lattice Boltzmann methods

We present high-resolution numerical simulations of convection in multiphase flows (boiling) using a novel algorithm based on a lattice Boltzmann method. We first study the thermodynamical and kinematic properties of the algorithm. Then, we perform a series of 3D numerical simulations changing the mean properties in the phase diagram and compare convection with and without phase coexistence at Rayleigh number Ra∼10(7). We show that in the presence of nucleating bubbles non-Oberbeck-Boussinesq effects develop, the mean temperature profile becomes asymmetric, and heat-transfer and heat-transfer fluctuations are enhanced, at all Ra studied. We also show that small-scale properties of velocity and temperature fields are strongly affected by the presence of the buoyant bubble leading to high non-gaussian profiles in the bulk.

Rayleigh-Bénard convection with phase changes in a Galerkin model

The transition to turbulence in Rayleigh-Bénard convection with phase changes and the resulting convective patterns are studied in a three-dimensional Galerkin model. Our study is focused on the conditionally unstable regime of moist convection in which the stratification is stable for unsaturated air parcels and unstable for saturated parcels. We perform a comprehensive statistical analysis of the transition to convection that samples the dependence of attractors (or fixed points) in the phase space of the model on the dimensionless parameters. Conditionally unstable convection can be initiated either from a fully unsaturated linearly stable equilibrium or a fully saturated linearly unstable equilibrium. Highly localized moist convection can be found in a steady state, in an oscillating recharge-discharge regime, or turbulent in dependence of the aspect ratio and the degree of stable stratification of the unsaturated air. Our phase-space analysis predicts parameter ranges for which self-sustained convective regimes in the case of subcritical conditional instability can be observed. The observed regime transitions for moist convection bear some similarities to transitions to turbulence in simple shear flows.

Effect of vapor bubbles on velocity fluctuations and dissipation rates in bubbly Rayleigh-Bénard convection

Numerical results for kinetic and thermal energy dissipation rates in bubbly Rayleigh-Bénard convection are reported. Bubbles have a twofold effect on the flow: on the one hand, they absorb or release heat to the surrounding liquid phase, thus tending to decrease the temperature differences responsible for the convective motion; but on the other hand, the absorbed heat causes the bubbles to grow, thus increasing their buoyancy and enhancing turbulence (or, more properly, pseudoturbulence) by generating velocity fluctuations. This enhancement depends on the ratio of the sensible heat to the latent heat of the phase change, given by the Jakob number, which determines the dynamics of the bubble growth.

Volumetric formulation for a class of kinetic models with energy conservation

We analyze a volumetric formulation of lattice Boltzmann for compressible thermal fluid flows. The velocity set is chosen with the desired accuracy, based on the Gauss-Hermite quadrature procedure, and tested against controlled problems in bounded and unbounded fluids. The method allows the simulation of thermohydrodyamical problems without the need to preserve the exact space-filling nature of the velocity set, but still ensuring the exact conservation laws for density, momentum, and energy. Issues related to boundary condition problems and improvements based on grid refinement are also investigated.