An official website of the United States government
This study explores heat and turbulent modulation in three-dimensional multiphase Rayleigh–Bénard convection using direct numerical simulations. Two immiscible fluids with identical reference density undergo systematic variations in dispersed-phase volume fractions,$$0.0 \leq \varPhi \leq 0.5$$, and ratios of dynamic viscosity,$$\lambda _{\mu }$$, and thermal diffusivity,$$\lambda _{\alpha }$$, within the range$$[0.1\unicode{x2013}10]$$. The Rayleigh, Prandtl, Weber and Froude numbers are held constant at$$10^8$$,$$4$$,$$6000$$and$$1$$, respectively. Initially, when both fluids share the same properties, a 10 % Nusselt number increase is observed at the highest volume fractions. In this case, despite a reduction in turbulent kinetic energy, droplets enhance energy transfer to smaller scales, smaller than those of single-phase flow, promoting local mixing. By varying viscosity ratios, while maintaining a constant Rayleigh number based on the average mixture properties, the global heat transfer rises by approximately 25 % at$$\varPhi =0.2$$and$$\lambda _{\mu }=10$$. This is attributed to increased small-scale mixing and turbulence in the less viscous carrier phase. In addition, a dispersed phase with higher thermal diffusivity results in a 50 % reduction in the Nusselt number compared with the single-phase counterpart, owing to faster heat conduction and reduced droplet presence near walls. The study also addresses droplet-size distributions, confirming two distinct ranges dominated by coalescence and breakup with different scaling laws.
more »
« less
Turbulent Rayleigh–Bénard convection in non-colloidal suspensions
This study presents direct numerical simulations of turbulent Rayleigh–Bénard convection in non-colloidal suspensions, with special focus on the heat transfer modifications in the flow. Adopting a Rayleigh number of $10^8$ and Prandtl number of 7, parametric investigations of the particle volume fraction $$0\leq \varPhi \leq 40\,\%$$ and particle diameter $$1/20\leq d^*_p\leq 1/10$$ with respect to the cavity height, are carried out. The particles are neutrally buoyant, rigid spheres with physical properties that match the fluid phase. Up to $$\varPhi =25\,\%$$ , the Nusselt number increases weakly but steadily, mainly due to the increased thermal agitation that overcomes the decreased kinetic energy of the flow. Beyond $$\varPhi =30\,\%$$ , the Nusselt number exhibits a substantial drop, down to approximately 1/3 of the single-phase value. This decrease is attributed to the dense particle layering in the near-wall region, confirmed by the time-averaged local volume fraction. The dense particle layer reduces the convection in the near-wall region and negates the formation of any coherent structures within one particle diameter from the wall. Significant differences between $$\varPhi \leq 30\,\%$$ and 40 % are observed in all statistical quantities, including heat transfer and turbulent kinetic energy budgets, and two-point correlations. Special attention is also given to the role of particle rotation, which is shown to contribute to maintaining high heat transfer rates in moderate volume fractions. Furthermore, decreasing the particle size promotes the particle layering next to the wall, inducing a similar heat transfer reduction as in the highest particle volume fraction case.
more »
« less
- Award ID(s):
- 1854376
- PAR ID:
- 10341515
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 945
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)This study explores thermal convection in suspensions of neutrally buoyant, non-colloidal suspensions confined between horizontal plates. A constitutive diffusion equation is used to model the dynamics of the particles suspended in a viscous fluid and it is coupled with the flow equations. We employ a simple model that was proposed by Metzger, Rahli & Yin ( J. Fluid Mech. , vol. 724, 2013, pp. 527–552) for the effective thermal diffusivity of suspensions. This model considers the effect of shear-induced diffusion and gives the thermal diffusivity increasing linearly with the thermal Péclet number ( Pe ) and the particle volume fraction ( ϕ ). Both linear stability analysis and numerical simulation based on the mathematical models are performed for various bulk particle volume fractions $$({\phi _b})$$ ranging from 0 to 0.3. The critical Rayleigh number $$(R{a_c})$$ grows gradually by increasing $${\phi _b}$$ from the critical value $$(R{a_c} = 1708)$$ for a pure Newtonian fluid, while the critical wavenumber $$({k_c})$$ remains constant at 3.12. The transition from the conduction state of suspensions is subcritical, whereas it is supercritical for the convection in a pure Newtonian fluid $$({\phi _b} = 0)$$ . The heat transfer in moderately dense suspensions $$({\phi _b} = 0.2\text{--}0.3)$$ is significantly enhanced by convection rolls for small Rayleigh number ( Ra ) close to $$R{a_c}$$ . We also found a power-law increase of the Nusselt number ( Nu ) with Ra , namely, $$Nu\sim R{a^b}$$ for relatively large values of Ra where the scaling exponent b decreases with $${\phi _b}$$ . Finally, it turns out that the shear-induced migration of particles can modify the heat transfer.more » « less
-
Bounding walls or immersed surfaces are utilized in many industrial systems as the primary thermal source to heat a gas–solids mixture. Previous efforts to resolve the solids’ heat transfer near a boundary involve the extension of unbounded convection correlations into the near-wall region in conjunction with particle-scale theories for indirect conduction. Moreover, unbounded drag correlations are utilized in the near-wall region (without modification) to resolve the force exerted on a solid particle by the fluid. We rigorously test unbounded correlations and indirect conduction theory against outputs from direct numerical simulation of laminar flow past a hot plate and a static, cold particle. Here, local variables are utilized for consistency with unresolved computational fluid dynamics discrete element methods and lead to new unbounded correlations that are self-similar to those obtained with free-stream variables. The new drag correlation with local fluid velocity captures the drag force in both the unbounded system as well as the near-wall region while the classic, unbounded drag correlation with free-stream fluid velocity dramatically over-predicts the drag force in the near-wall region. Similarly, classic, unbounded convection correlations are found to under-predict the heat transfer occurring in the near-wall region. Inclusion of indirect conduction, in addition to unbounded convection, performs markedly better. To account for boundary effects, a new Nusselt correlation is developed for the heat transfer in excess of local, unbounded convection. The excess wall Nusselt number depends solely on the dimensionless particle–wall separation distance and asymptotically decays to zero for large particle–wall separation distances, seaming together the unbounded and near-wall regions.more » « less
-
null (Ed.)The movement of heat in a convecting system is typically described by the nondimensional Nusselt number, which involves an average over both space and time. In direct numerical simulations of turbulent flows, there is considerable variation in the contributions to the Nusselt number, both because of local spatial variations due to plumes and because of intermittency in time. We develop a statistical approach to more completely describe the structure of heat transfer, using an exit-distance extracted from Lagrangian tracer particles, which we call the Lagrangian heat structure. In a comparison between simulations of homogeneous turbulence driven by Boussinesq convection, the Lagrangian heat structure reveals significant non-Gaussian character, as well as a clear trend with Prandtl number and Rayleigh number. This has encouraging implications for simulations performed with the goal of understanding turbulent convection in natural settings such as Earth’s atmosphere and oceans, as well as planetary and stellar dynamos.more » « less
-
The central open question about Rayleigh–Bénard convection – buoyancy-driven flow in a fluid layer heated from below and cooled from above – is how vertical heat flux depends on the imposed temperature gradient in the strongly nonlinear regime where the flows are typically turbulent. The quantitative challenge is to determine how the Nusselt number $Nu$ depends on the Rayleigh number $Ra$ in the $$Ra\to \infty$$ limit for fluids of fixed finite Prandtl number $Pr$ in fixed spatial domains. Laboratory experiments, numerical simulations and analysis of Rayleigh's mathematical model have yet to rule out either of the proposed ‘classical’ $$Nu \sim Ra^{1/3}$$ or ‘ultimate’ $$Nu \sim Ra^{1/2}$$ asymptotic scaling theories. Among the many solutions of the equations of motion at high $Ra$ are steady convection rolls that are dynamically unstable but share features of the turbulent attractor. We have computed these steady solutions for $Ra$ up to $$10^{14}$$ with $Pr=1$ and various horizontal periods. By choosing the horizontal period of these rolls at each $Ra$ to maximize $Nu$ , we find that steady convection rolls achieve classical asymptotic scaling. Moreover, they transport more heat than turbulent convection in experiments or simulations at comparable parameters. If heat transport in turbulent convection continues to be dominated by heat transport in steady rolls as $$Ra\to \infty$$ , it cannot achieve the ultimate scaling.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2022.534
