Hostname: page-component-78c5997874-ndw9j Total loading time: 0 Render date: 2024-11-08T05:35:37.044Z Has data issue: false hasContentIssue false

Equilibrated crater: fragmentation and mixing

Published online by Cambridge University Press:  13 May 2022

E. Villermaux*
Affiliation:
Aix Marseille Université, CNRS, Centrale Marseille, IRPHE UMR 7342, 13384 Marseille, France Institut Universitaire de France, 75005 Paris, France
*
Email address for correspondence: [email protected]

Abstract

The composition of the Earth's mantle, built by impacts of metallic core planetesimals, depends on the time offered to metal–silicate equilibration (concentration equalization by mixing) as the impactor fragments settle down. Smaller fragments equilibrate while bigger ones remain segregated and accumulate in the – therefore iron rich – Earth's core. To understand the primary instability setting the interpenetration depth and wavelength between the impactor and the mantle, Lherm et al. (J. Fluid Mech., vol. 937, 2022, A20) analyse the intermingled phenomena occurring when a drop impacts a substrate of a pool of another liquid, digging a crater. The density difference between the phases triggers a Rayleigh–Taylor instability when the drop is denser than the pool. The overall planet chemical composition relies on a tradeoff between fragmentation and mixing; this study provides elements to decipher their relative importance.

Type
Focus on Fluids
Copyright
© The Author(s), 2022. Published by Cambridge University Press

1. Introduction

We have known since Newton that in order to move a massive object from rest, or slow it down to immobility, a force must be applied to it (Newton Reference Newton1687, Law 1). In fluids, forces are measured by pressure gradients, and Lherm et al. (Reference Lherm, Deguen, Alboussière and Landeau2022) recall that when the density and pressure gradients are in opposite directions at an interface between two fluids of different densities, then the interface is unstable. This is the Rayleigh (Reference Rayleigh1883)–Taylor (Reference Taylor1950) instability paradigm, a ubiquitous phenomenon with an immense range of applications when it comes to understanding how different phases interpenetrate, mix or fragment in non-Galilean (i.e. accelerated or decelerated) frames.

According to this principle, when a deformable object (a liquid drop) impacts a deformable substrate (a pool of another liquid), a Rayleigh–Taylor instability is likely to occur if the impacting drop is denser than the pool, a consequence of the drop/pool interface deceleration consecutive to the impact. This is actually what happens, as Lherm et al. (Reference Lherm, Deguen, Alboussière and Landeau2022) have nicely shown (see figure 1b) by documenting and analysing step by step the many concomitant, interconnected processes at play in this phenomenon, which is simple in its principle, but turns out to be complicated (as the author's analysis is) when all the imbricated details are handled quantitatively. Elements of the phenomenology had been identified long ago (Thomson & Newall Reference Thomson and Newall1885; Arecchi et al. Reference Arecchi, Buah-Bassuah, Francini and Residori1996) in what the authors call the ‘static’ limit, when destabilization is due to permanent body forces (gravity g, see figure 1a).

Figure 1. (a) A drop of heavy liquid deposited at the surface of a lighter pool forms a vortex ring which soon fragments into secondary rings; it would not do so if it were lighter than the pool, as first noticed by Thomson & Newall (Reference Thomson and Newall1885). (b) Patterns of the cavity (crater) at maximal extension following the impact of a drop of density $\rho _1$ with radius $R_0$ on a pool with density $\rho _2$ at velocity $u$, defining the Froude number $Fr=u^{2}/gR_0$. Lherm et al. (Reference Lherm, Deguen, Alboussière and Landeau2022) clearly demonstrate here that, when $\rho _1/\rho _2<1$, the cavity is smooth while, for $\rho _1/\rho _2>1$, corrugations have developed at its surface, which are more pronounced for larger $\rho _1/\rho _2$. (c) Collapse of a hydrogen–oxygen mixture ignited in a bubble immersed in a liquid, and Rayleigh–Taylor corrugations of its interface as it rebounds on the compressed burnt gases (adapted from Duplat & Villermaux Reference Duplat and Villermaux2015).

The broad picture is that the momentum carried by the drop is absorbed by digging a cavity – a crater in the author's terminology – underneath the pool's surface, moving the fluid around in a close-to-radial fashion so that the drop crunches at the inner surface of a radially expanding, nearly spherical bubble. The expansion decelerates continuously such that if $R(t)$ is the radius of the cavity, $\dot {R}>0$ but $\ddot {R}<0$. This is one out of the many examples of a Rayleigh–Taylor instability in a spherical geometry, other notorious examples including the explosion of shell cases and the implosion–explosion of supernovae (Villermaux Reference Villermaux2020, see § 5.1 for a review), or collapsing nearly empty bubbles which, at rebound when their interior has been sufficiently compressed, develop surface corrugations, a phenomenon known to compromise the efficiency of some inertial confinement fusion processes (Smalyuk et al. Reference Smalyuk2014; Duplat & Villermaux Reference Duplat and Villermaux2015, and figure 1c). Note that the instability in question here operates at the interface between the cavity and the pool; it should not be confused with the one affecting the border of the liquid sheet ejected from a droplet impacting a solid surface (Allen Reference Allen1975), although both instabilities are of the same nature, and both occur on a decelerating, radially expanding substrate.

2. Overview

The authors interest is, however, in geophysics. Planets like the Earth were built by high-energy impacts of planetesimals with metallic cores of their own. The composition of the Earth's mantle depends critically on the time offered to metal–silicate chemical equilibration (meaning concentration equalization through mixing) as the impactor material settles towards the Earth's core. Efficient equilibration is achieved with smaller fragments, allowing fast metal–silicate mass transfer, while the fate of bigger ones is to accumulate in the – therefore iron rich – core (Deguen, Landeau & Olson Reference Deguen, Landeau and Olson2014; Landeau, Deguen & Olson Reference Landeau, Deguen and Olson2014). The overall planets chemical composition thus relies on the initial fragmentation state of the crunched impactors since this initial division state will determine in all cases the subsequent molecular mixing kinetics of the fragments in the molten silicate magma. This is the reason why documenting primary instabilities, such as the typical wavelength of the contacting Rayleigh–Taylor instability, which sets the length scale at which the impactor and the mantle interpenetrate, and the net depth of interpenetration, is paramount.

Through delicate but precise enough experiments in the spirit of those of the early contributors to the field (Lewis Reference Lewis1950; Emmons, Chang & Watson Reference Emmons, Chang and Watson1960; Niederhaus & Jacobs Reference Niederhaus and Jacobs2003), the authors embark on this endeavour. Admittedly, this is not an easy task since, in this problem, the instability is caused by a time varying deceleration ($\ddot {R}<0$ decays in time but not fast enough to be considered as impulsional) on a expanding substrate ($\dot {R}>0$) and applies to a fluid layer (the crunched drop pancake covering the cavity) which itself gets thinner in time. The impacting drop inertia triggers the pool motion, which is counterbalanced by buoyancy because gravity tends to level out the cavity. With the help of fudge factors introduced to compensate for the lack of energy conservation in the formulation, which is basically conservative in a close-to-spherical but not exactly spherical geometry, the authors nevertheless make sense of the cavity dynamics, and derive from it a fair caricature of the crunched drop layer destabilization; its wavelength is controlled by inertia, buoyancy and the viscosity of the pool, a result which was not obvious from the start.

The only ingredient the authors neglect is capillarity. Considering the application they have in mind, namely the formation of a crater by the impact of an asteroid falling on Earth (for which the Weber number is astronomical), this neglect is legitimate; but in view of their much more gentle experimental conditions, the non-importance of surface tension is less obvious. There is a special limit of the Rayleigh–Taylor instability for thin layers (as opposed to an interface separating two semi-infinite phases) when surface tension operates (Keller & Kolodner Reference Keller and Kolodner1954; Bremond & Villermaux Reference Bremond and Villermaux2005). Here, a dense decelerated layer is sandwiched between two lighter phases (the pool and air). We know that the proximity of the two interfaces on each side of the layer not only affects the growth rate, but also the instability mode selection, both being a function of the layer thickness. The authors disregard this effect because the preferred wavelength in their experiments is of the order of, or smaller than, and in any case independent of it. But this effect does exist, and has been shown to fully operate in systems very close to theirs, in a true spherical geometry (Vledouts et al. Reference Vledouts, Quinard, Vandenberghe and Villermaux2016).

3. Future

The study of Lherm et al. (Reference Lherm, Deguen, Alboussière and Landeau2022) culminates with an estimate of the pool material mass entrained through the instability development. It can be ten times as large as the impacting drop mass; this provides a measure of the dilution level the metal/silicate mixture will have once molecular mixing (equilibration) is achieved. This entrainment index is the equivalent of the ‘mixing layer thickness’ in shear flows; it gives an idea of the extent of the inter-dispersion region between the phases, but says nothing about their fine-scale interpenetration at the molecular scale, which is mixing per se. The phases could well remain segregated from each other, however, chemically equilibrated implies molecularly mixed. The authors are well aware of the distinction between stirring and mixing (Villermaux Reference Villermaux2019) and know, in some prototypical situations, how to articulate one with the other (Lherm & Deguen Reference Lherm and Deguen2018). This important aspect remains to be quantified in the present configuration, a goal which should be attainable given the quality of their diagnostics.

Also, knowing the instability length scale in the interpenetrated region between the phases is useful because it sets the typical stirring scale for describing the mixing kinetics further, but we also know that fragmentation phenomena give rise to broad distributions of fragments sizes (Lhuissier et al. Reference Lhuissier, Sun, Prosperetti and Lohse2013, although capillary aggregation certainly plays a role there). Identifying the proper instability at play, and the emergence of its preferred length scale is one good thing, but describing the fragmentation process in its diversity and the nature of the induced, competing mixing of the fragments assembly is another level, certainly on the author's agenda for the future.

Declaration of interests

The author reports no conflict of interest.

References

REFERENCES

Allen, R.F. 1975 The role of surface tension in splashing. J. Colloid Interface Sci. 51 (2), 350351.CrossRefGoogle Scholar
Arecchi, F.T., Buah-Bassuah, P.K., Francini, F. & Residori, S. 1996 Fragmentation of a drop as it falls in a lighter miscible fluid. Phys. Rev. E 54 (1), 424429.CrossRefGoogle Scholar
Bremond, N. & Villermaux, E. 2005 Bursting thin liquid films. J. Fluid Mech. 524, 121130.CrossRefGoogle Scholar
Deguen, R., Landeau, M. & Olson, P. 2014 Turbulent metal–silicate mixing, fragmentation, and equilibration in magma oceans. Earth Planet. Sci. Lett. 391, 274287.CrossRefGoogle Scholar
Duplat, J. & Villermaux, E. 2015 Luminescence from collapsing centimeter bubbles expanded by chemical reaction. Phys. Rev. Lett. 115, 094501.CrossRefGoogle ScholarPubMed
Emmons, H.W., Chang, C.T. & Watson, B.C. 1960 Taylor instability of finite surfaces waves. J. Fluid Mech. 7 (2), 177193.CrossRefGoogle Scholar
Keller, J.B. & Kolodner, I. 1954 Instability of liquid surfaces and the formation of drops. J. Appl. Phys. 25 (7), 918921.CrossRefGoogle Scholar
Landeau, M., Deguen, R. & Olson, P. 2014 Experiments on the fragmentation of a buoyant liquid volume in another liquid. J. Fluid Mech. 749, 478518.CrossRefGoogle Scholar
Lewis, D.J. 1950 The instability of liquid surfaces when accelerated in a direction perpendicular to their planes. II. Proc. R. Soc. Lond. A 202, 8196.Google Scholar
Lherm, V. & Deguen, R. 2018 Small-scale metal/silicate equilibration during core formation: the influence of stretching enhanced diffusion on mixing. J. Geophys. Res. 123 (12), 1049610516.CrossRefGoogle Scholar
Lherm, V., Deguen, R., Alboussière, T. & Landeau, M. 2022 Rayleigh–Taylor instability in impact cratering experiments. J. Fluid Mech. 937, A20.CrossRefGoogle Scholar
Lhuissier, H., Sun, C., Prosperetti, A. & Lohse, D. 2013 Drop fragmentation at impact onto a bath of an immiscible liquid. Phys. Rev. Lett. 110, 264503.CrossRefGoogle ScholarPubMed
Newton, I. 1687 Principia Mathematica.Google Scholar
Niederhaus, C.E. & Jacobs, J.W. 2003 Experimental study of the Richtmyer-Meshkov instability of incompressible fluids. J. Fluid Mech. 485, 243277.CrossRefGoogle Scholar
Rayleigh, Lord 1883 Investigation on the character of the equilibrium of an incompressible heavy fluid of variable density. Proc. R. Soc. Lond. A 14, 170177.Google Scholar
Smalyuk, V.A., et al. 2014 Measurements of an ablator-gas atomic mix in indirectly driven implosions at the national ignition facility. Phys. Rev. Lett. 112, 025002.CrossRefGoogle ScholarPubMed
Taylor, G.I. 1950 The instability of liquid surfaces when accelerated in a direction perpendicular to their planes. I. Proc. R. Soc. Lond. A 201 (1065), 192196.Google Scholar
Thomson, J.J. & Newall, H.F. 1885 On the formation of vortex rings by drops falling into liquids, and some allied phenomena. Proc. R. Soc. Lond. 39, 417436.Google Scholar
Villermaux, E. 2019 Mixing versus stirring. Annu. Rev. Fluid Mech. 51, 245273.CrossRefGoogle Scholar
Villermaux, E. 2020 Fragmentation versus cohesion. J. Fluid Mech. 898, P1.CrossRefGoogle Scholar
Vledouts, A., Quinard, J., Vandenberghe, N. & Villermaux, E. 2016 Explosive fragmentation of liquid shells. J. Fluid Mech. 788, 246273.CrossRefGoogle Scholar
Figure 0

Figure 1. (a) A drop of heavy liquid deposited at the surface of a lighter pool forms a vortex ring which soon fragments into secondary rings; it would not do so if it were lighter than the pool, as first noticed by Thomson & Newall (1885). (b) Patterns of the cavity (crater) at maximal extension following the impact of a drop of density $\rho _1$ with radius $R_0$ on a pool with density $\rho _2$ at velocity $u$, defining the Froude number $Fr=u^{2}/gR_0$. Lherm et al. (2022) clearly demonstrate here that, when $\rho _1/\rho _2<1$, the cavity is smooth while, for $\rho _1/\rho _2>1$, corrugations have developed at its surface, which are more pronounced for larger $\rho _1/\rho _2$. (c) Collapse of a hydrogen–oxygen mixture ignited in a bubble immersed in a liquid, and Rayleigh–Taylor corrugations of its interface as it rebounds on the compressed burnt gases (adapted from Duplat & Villermaux 2015).