Eur. Phys. J. E (2019) 42: 126
https://doi.org/10.1140/epje/i2019-11882-y

Regular Article

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

¹ Fluid Dynamics Laboratory, Department of Applied Physics, Eindhoven University of Technology, P.O. Box 513, 5600 MB, Eindhoven, The Netherlands
² Centre of Analysis, Scientific Computing, and Applications W&I, Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, 5600 MB, Eindhoven, The Netherlands
³ Istituto per le Applicazioni del Calcolo, Consiglio Nazionale delle Ricerche, Via dei Taurini 19, 00185, Rome, Italy

^* e-mail: k.m.j.alards@tue.nl
^** e-mail: F.Toschi@tue.nl

Received: 28 March 2018
Accepted: 5 August 2019
Published online: 12 September 2019

Abstract

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, , and the thermal expansion coefficient of particles relative to that of the fluid, . In particular, the residence time is constant for small thermal response times, , and increasing with for larger thermal response times, . 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 for larger response times, thus explaining the trends in the DNS data well.