Papers › Residence time of inertial particles in 3D thermal convection: implications for magma...
Residence time of inertial particles in 3D thermal convection: implications for magma reservoirs
Vojtěch Patočka, Nicola Tosi, Enrico Calzavarini
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
The dynamic behavior of crystals in convecting fluids determines how magma bodies solidify. In particular, it is often important to estimate how long crystals stay in suspension in the host liquid before being deposited at its bottom (or top, for light particles). We perform a systematic 3D numerical study of particle-laden Rayleigh-Benard convection, and derive a robust model for the particle residence time. For Rayleigh numbers higher than 10^7, inertial particles' trajectories exhibit a monotonic transition from fluid tracer-like to free-fall dynamics, the control parameter being the ratio between particle Stokes velocity and the fluid velocity. The average settling rate is proportional to the particle Stokes velocity in both the end-member regimes, but the distribution of the residence times differs markedly from one to the other. For lower Rayleigh numbers (<10^7), an interaction between large-scale circulation and particle motion emerges, increasing the settling rates on average. Nevertheless, the mean residence time does not exceed the terminal time, i.e. the settling time from a quiescent fluid, by a factor larger than four. An exception are simulations with only a slightly super-critical Rayleigh number (~ 10^4), for which stationary convection develops and some particles become trapped indefinitely. 2D simulations of the same problem overestimate the flow-particle interaction - and hence the residence time - for both high and low Rayleigh numbers, which stresses the importance of using 3D geometries for simulating particle-laden flows. We outline how our model can be used to explain depth changes of crystal size distribution in sedimentary layers of magmatic intrusions that are thought to have formed via settling of a crystal cargo, and discuss how the micro-structural observations of solidified intrusions can be used to infer the past convective velocity of magma.
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections