https://www.nature.com/articles/s41467-024-47672-z Skip to main content Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript. Advertisement Advertisement Nature Communications * View all journals * Search * Log in * Explore content * About the journal * Publish with us * Sign up for alerts * RSS feed 1. nature 2. nature communications 3. articles 4. article Levitation and dynamics of bodies in supersaturated fluids Download PDF Download PDF * Article * Open access * Published: 09 May 2024 Levitation and dynamics of bodies in supersaturated fluids * Saverio E. Spagnolie ORCID: orcid.org/0000-0002-8484-1573^1,2, * Samuel Christianson^1 & * Carsen Grote^1 Nature Communications volume 15, Article number: 3910 (2024) Cite this article * 1029 Accesses * 453 Altmetric * Metrics details Subjects * Chemical engineering * Fluid dynamics * Volcanology Abstract A body immersed in a supersaturated fluid like carbonated water can accumulate a dynamic field of bubbles upon its surface. If the body is mobile, the attached bubbles can lift it upward against gravity, but a fluid-air interface can clean the surface of these lifting agents and the body may plummet. The process then begins anew, and continues for as long as the concentration of gas in the fluid supports it. In this work, experiments using fixed and free immersed bodies reveal fundamental features of force development and gas escape. A continuum model which incorporates the dynamics of a surface buoyancy field is used to predict the ranges of body mass and size, and fluid properties, for which the system is most dynamic, and those for which body excursions are suppressed. Simulations are then used to probe systems which are dominated by a small number of large bubbles. Body rotations at the surface are critical for driving periodic vertical motions of large bodies, which in turn can produce body wobbling, rolling, and damped surface 'bouncing' dynamics. Similar content being viewed by others [41586_2020] Floating under a levitating liquid Article 02 September 2020 [41598_2020] The dynamics of gas-bubble formation at saturated conditions in porous media flow Article Open access 06 August 2020 [41598_2022] Coalescence of two growing bubbles in a Hele-Shaw cell Article Open access 24 January 2022 Introduction A fluid containing dissolved gas may become supersaturated upon a rapid change in temperature or pressure, leading to bubble formation and eventual gas escape. This phenomenon is most commonly observed when opening a can of sparkling water, or other carbonated beverages. When the fluid is under sufficient pressure, gas accumulation is inhibited, preventing bubble formation^1. Upon a rapid reduction of pressure, bubbles form on or near any containing or immersed surfaces, then detach and depart towards the fluid-air interface, leading to the eventual escape of the gas to the environment^2,3,4. In everyday settings, rather than forming directly upon container walls, bubbles commonly form inside small cellulose fibers left behind during cleaning^5,6. These cavities ('lumen') are remarkable sites for consistent and rapid bubble growth and release since a pocket of gas remains behind after each pinch-off event to seed the next growth. The coalescence of diffusively growing bubbles itself presents a classical modeling challenge^7,8, as is the process by which coalescing bubbles depart from a wall^9,10,11,12. See in particular the reviews by Liger-Belair^13 and Lohse^14. Supersaturated fluids also appear in geophysical settings. Explosive fragmentation of particulate matter in magma can lead to volcanic eruptions. A Strombolian eruption is caused by bubbles which coalesce into rising 'gas slugs', transporting gas and entraining magma flow towards the surface^15,16,17. Such liquids also appear in industrial processes like deacidification and fractionation of oils^18,19, and in biological settings (e.g., blood and tissues during decompression) ^20,21. Large-scale flows associated with bubble motion depend on the geometry of the container, and can result in peculiar downward bubble motion and even bubble waves and cascades in fluids from stout beers to magma^22,23,24. When a free body is introduced into a supersaturated environment it presents new sites at which bubbles can nucleate and grow - a field of such bubbles on a body surface can result in surprising dynamics. Using nothing more than carbonated water and raisins (Fig. 1), periodic vertical body motions can be observed for nearly two hours, though the time between excursions slows considerably after the first 20 min (see Movie S1). This phenomenon has earned numerous playful names by a variety of fizziologists^25, from dancing raisins to divers^26,27 and fizz balls^28. Similar body oscillations can also be generated by reactions in chemical gardens^29,30. Recently, this phenomenon has been investigated using peanuts and beer, based on a common practice in Argentinian drinking establishments^31. Different contact angles of individual attached bubbles show the importance of surface roughness on bubble formation^31. Such effects add to a growing list of fundamental scientific findings which have emerged from the kitchen^32,33. Fig. 1: Dancing raisins. figure 1 Raisins in carbonated water present numerous folds conducive to bubble nucleation and growth; these bubbles may then lift the body upward against gravity, only to release it upon arrival at the free surface (see Movie S1). Full size image In this paper, we explore these oscillating dynamics using experiments and simulations and match numerous features to theoretical predictions. Experiments are used to measure the force development on a spherical body fixed in carbonated water, and its oscillatory dynamics when free. Using a discrete bubble model, and a continuum model which incorporates the dynamics of a surface buoyancy field, ranges of body and fluid properties are provided for which the system is most dynamic, and those for which body excursions are suppressed. Simulations are used to explore the dependence on system parameters in a more controlled setting. Body rotations are found to be critical for the onset of periodic vertical motions of large bodies, which in turn can produce body wobbling, rolling, and multi-period surface return in a damped bouncing dynamics. Results Experiments Mass loss of a supersaturated fluid We first measured the change in gas concentration in a supersaturated fluid upon depressurization by examining the fluid's mass loss over time. An empty glass vessel with a square cross-section of side length L = 8.9 cm was filled with a just-opened can of Klarbrunn-brand carbonated water which was stored at room temperature (21.6 degC). The fluid was poured into the vessel and then left alone for two hours at the same room temperature. Time, denoted by t, is measured in minutes, and t = 0 corresponds to the time just after depressurization (the moment when the can was opened). The mean mass loss across five experiments at each time is shown in Fig. 2a as a thin red curve. In each experiment the initial volume of fluid, V[f], was approximately 355 cm^3 (12 oz); the total mass loss, M[e](t), over the course of two hours was roughly 0.3% of the initial mass. Fig. 2: Gas loss from an empty container of carbonated water. figure 2 a The mass loss in grams of one can of carbonated water during the first 80 min after depressurization (thin red line) at time t averaged over four independent experiments, and the mass loss after subtracting a constant linear loss due to evaporation (dark blue line). Error bars show the standard deviation across experiments at each time. b The CO[2] gas supersaturation ratio, minus a stationary minimum value below which bubbles do not form on the container, S(t) - S[mc], inferred from the data in (a) (dark blue line). The dashed line is a best-fit curve using (1). Full size image Rapid mass loss over the first 20 min was due to the growth and detachment of bubbles on the container surface, and diffusive transport of gas from the free surface, as discussed below. The remaining 100 min revealed linear behavior due to water evaporation^ 34. At long times we observed that dM[e]/dt [?] k[e] for relatively large t, where k[e] = 4 x 10^-3 g/min (the same value was found using only tap water). The solid blue curve in Fig. 2a shows M[e](t)-k[e]t, expected to be the mass loss of CO[2], and the standard deviation at each time across experiments is shown using error bars. Denoting the volume-averaged gas concentration at time t by \(\bar{c} (t)\), with units of g/cm^3, the mass of CO[2] in the vessel is given by \(\bar{c}(t){V}_{f}\) (an adjustment due to evaporative volume change is negligible). The mass loss from the empty vessel is written as \({M}_{e}(t)={V}_{f}(\bar{c}(0)-\bar{c}(t))+{k}_{e}t\), where \(\ bar{c}(0)\) is the initial concentration of CO[2]. The supersaturation ratio, \(S(t)=(\bar{c}(t)/{c}_{s}-1)\), measures the gas concentration level compared to a critical value c[s], an equilibrium concentration corresponding to a partial pressure of gaseous CO[2] at 1 atm^13,35. At room temperature in water this value is c[s] = 1.48 g/L^36,37. The gas concentration decreases until S(t) reaches a value below which bubbles no longer form on or near the container walls; in this case, most likely in the lumen of adhered fibrous matter^5,6. We denote this minimum value associated with the container as S[mc] and discuss its genesis more fully in SSII B. The difference S(t) - S[mc] may then be inferred from the mass loss data, and is plotted in Fig. 2b. The dynamics of S are modeled by a Riccati equation, \(\dot{S}=-(S-{S}_{mc})/{T}_{r}-q{(S-{S}_{mc})}^{2}\) with the dot denoting a time derivative, for reasons to be discussed, which results in a predicted evolution of the form $$S(t)={S}_{mc}+\frac{({S}_{0}-{S}_{mc})\exp (-t/{T}_{r})}{1+\chi (1- \exp (-t/{T}_{r}))},$$ (1) where S[0] = S(0) is the supersaturation ratio at t = 0, and kh = qT[r ](S[0] - S[mc]). This function is plotted atop the experimental measurements in Fig. 2b as a dashed curve using fitted parameters S [0] - S[mc] = 1.66, kh = 2.5, and T[r] = 36.2 min. With S[0] - S[mc] [?] S[0], the observations here are in line with carbonation levels reported in other sparkling beverages^4,38, particularly given the substantial gas loss while pouring^4,39. Surface buoyancy growth with a fixed body Next, a body was held fixed in the fluid, and we measured the force development due to bubble growth on its surface. A sphere of radius 1 cm composed of Polylactic acid (PLA) was printed with an Ultimaker 3+ 3D printer using a 0.15 mm nozzle. A schematic of the experimental setup is shown in Fig. 3a. The fluid was depressurized and poured gently into the container, and placed upon the digital scale. The test sphere was then lowered into the carbonated water at different insertion times, denoted by t[0], and affixed to the table below. As bubbles grew on the surface (see Fig. 3e) they imposed an upward vertical force on the stationary sphere, and the opposing downward force was registered by the scale. The sphere was briefly removed from the fluid every 4 min, then reinserted while still wet. Videos of bubble growth, coalescence, and arrival at a fluctuating steady state at insertion times t[0] = 1 min and t[0] = 10 min are included as Movies S2-S3. Fig. 3: Measuring the surface buoyancy on a fixed body. figure 3 a Experimental setup for measuring the surface buoyancy. A spherical 3D-printed body is held fixed in place in carbonated water for 4 min, then removed from and reinserted into the fluid. Bubbles form, grow, merge and detach, and the force is registered by the scale. b The surface buoyancy force, B(t; t[0]) during a single experiment. The force increases at a rate l(t[0]) which depends on the insertion time t[0], tapering off to a value B[s](t[0]) (measured later, at \({t}_ {0}+4\min\)). The rate l(t[0]) is computed using the lines shown. c The growth rate, l(t[0]). The mean and standard deviation across five experiments are shown as symbols and error bars. The solid curve is based on a model which assumes l [?] S^3/2 and using (1). d The stabilized surface buoyancy force, B[s], measured at \({t}_{0}+4\min \), along with a linear fit. e Bubble configuration on the spherical surface at three different moments after an early insertion time of t [0] = 1 min. Full size image The scale-registered weight of the system, denoted by F(t), was decomposed as F(t) = F(0) + M[e](t)g + B(t; t[0]). F(0) includes a constant buoyancy force due to the volume of the displaced fluid; M[e ](t) is the mass change of a body-free fluid previously described; and B(t; t[0]) is the surface buoyancy force, the contribution to the buoyancy by the bubbles since the insertion time, t[0]. The weight loss due to gas escape from an empty vessel is small, but it is of comparable magnitude to the force on the body for roughly the first 10 min, and must be accounted for here. F(t) - F(0) was measured directly, which then provided an indirect measurement of the added buoyancy force, B(t; t[0]). Figure 3b shows the surface buoyancy force as a function of time since insertion, t - t[0], for insertion times up to 60 min. Upon insertion the force increased, at a rate which diminished as the insertion time became larger and the fluid became calmer. As discussed in SSII B, the radius a(t) of an individual bubble is expected to grow roughly as \(a(t)\, \approx \, {(2DS({t}_{0})(t-{t}_ {0}))}^{1/2}\), where D is the diffusion constant for CO[2] at room temperature. The associated volumetric growth suggests a total force on the body proportional to \({(t-{t}_{0})}^{3/2}\); however, once heterogeneous bubble coalescence begins the growth appears closer to linear in time (seen for insertion times less than 30 min in Fig. 3 b). This growth rate, denoted by l(t[0]), is shown in Fig. 3c. The solid curve corresponds to a prediction \(\lambda ({t}_{0})=kS{({t}_ {0})}^{3/2}\), where k is a proportionality constant, revealing an excellent match to the experimentally measured growth rates. Equivalently, defining l[0] to be the linear growth rate at t = 0, we may write $$\lambda ({t}_{0})={\lambda }_{0}{\left(\frac{S({t}_{0})}{{S}_{0}}\ right)}^{3/2}={\lambda }_{0}{\left(\frac{(S({t}_{0})-{S}_{mc})+{S}_ {mc}}{({S}_{0}-{S}_{mc})+{S}_{mc}}\right)}^{3/2}.$$ (2) Using S(t) - S[mc] from the previous experiment with no immersed body (see Fig. 2b), the model best fits the data on a logarithmic scale using l[0] = 58.4 dyn/s and S[mc] = 0.020. By comparing with the previous experiment, we find S[0] [?] 1.68. The surface buoyancy force eventually stabilized to a value B[s](t [0]) measured four minutes after reinsertion time t[0]. This value is shown in Fig. 3d along with a linear fit, B[s](t[0]) [?] 712 dyn- (6.9 dyn/min)t. The saturated value B[s](t[0]) diminished far more slowly than did the growth rate l(t[0]) - roughly, at later insertion times, bubbles merely take longer to accumulate and grow until reaching a critical size, at which time they pinch off and depart alone. Whether the slow decay was due to surface wetting or other phenomena remains unclear. A simple method for estimating B[s] is included in the Supplementary Information. After the buoyancy force stabilized, large fluctuations were observed. They were most apparent at smaller insertion times due to large bubble surface sliding and detachment events, which carry numerous other bubbles away at the same time. This surface cleaning effect has recently seen more specific attention^12. Detachment is expected beyond a critical 'Fritz radius' where buoyancy overwhelms capillary forces^31,40,41. Bubbles have other opportunities to depart during merging events via self-propelled detachment^9. The same force development measurement was performed using a skewer of 8 Sunmaid raisins, as reported in the Supplementary Information. The initial maximum force per raisin due to bubbles was B[s] [?] 100 dyn, and the initial growth rate per raisin was l[0] [?] 20 dyn/s. Levitation and dynamics of a freely moving sphere In a final set of experiments, the printed sphere was placed into the fluid, free to move (Fig. 4a). The body, with radius A = 1 cm and mass m = 4.25 g, was inserted into the fluid just after depressurization and pouring, and its motion was recorded with a Nikon D7000 DSLR camera with a 24 fps framerate for two hours (see Movies S4-S6). The body's vertical position at each time, denoted by z = A Z(t) (with the fluid surface at z = 0), was recovered using an image tracking code written in Matlab. The experiment was performed five times. Fig. 4: Long-time dynamics of a free body. figure 4 a A freely moving 3D-printed sphere of radius 1 cm and mass 4.25 g in carbonated water, rotates and returns to the free surface. b The locally averaged position over a window of +- 5 min for a single experiment is shown as a dark curve; the standard deviation is shown with error bars. c The instantaneous vertical position from (b). The body tends to visit the surface numerous times in succession, clearing off more and more of the surface, before plummeting. d At later times, body rising events are often cut short by premature detachment of large individual bubbles. e The dancing frequency, f [?] 1/D(t), where D(t) measures the averaged excursion time (over a +- 5 min window) from surface-departure to surface return. The solid curve shows the mean over three experiments, error bars show the standard deviation across experiments, and the dashed lines are from (11). Inset: the same, on a log-linear scale. Full size image A slow decay of the highly oscillatory dynamics is observed upon locally averaging the vertical position. Figure 4b shows the mean vertical position at time t over a window [t - 5 min, t + 5 min] during one experimental run. The free surface is located at z = 0, and the body is just touching the free surface when Z = -1. For the first 5 min the body spends most of its time near the free surface-bubbles cleaned from the body are rapidly replaced. Eventually, the body begins to spend more time in the bulk fluid, and then on the container floor, where the body resides for longer and longer 'charging' times before rising again. The instantaneous vertical position starting at t = 24.7 min is shown in Fig. 4c, revealing a type of damped 'bouncing' dynamics from the free surface. Multiple visits to the fluid-air interface in succession are commonly observed; each visit clears different portions of the body surface, resulting eventually in a longer excursion into the bulk. Figure 4d shows common late-stage behavior: long periods on the container floor, punctuated by increasingly rare rising events. Rising often ends prematurely due to the detachment and departure of single large bubbles. The late-stage surface forces tend to be dominated by such individuals, and the body's fate can be dictated by their singular activity. An excursion is defined as having occurred if the body's center of mass drops below one diameter beneath the free surface before returning upward and crossing the same vertical threshold. The times at which the jth excursion begins and ends are denoted respectively by t[j] (Z(t[j]) = -2 and \({Z}^{{\prime} }({t}_{j}) \, < \, 0\)), and p[j] (Z(p[j]) = -2 and \({Z}^{{\prime} }({p}_{j}) \, > \, 0\)). The jth excursion time is given by D[j] [?] p[j] - t[j], for j >= 1. The body undergoes roughly 300-700 excursions during the two-hour experiment. A smoothed excursion time, D(t), is defined as another moving average using +- 5 min on either side of time t. The dancing frequency, f [?] 1/D, is shown in Fig. 4e, and on a log-linear scale in the inset, including error bars representing one standard deviation at that time across three independent experiments. A clear monotonic decrease in the frequency is observed with nearly exponential decay over time. The dashed curve corresponds to a theoretical prediction discussed below. For bodies with a dense surface bubble coverage, body rotation at the surface plays a critical role. Early experiments, not shown here, constrained the motion to translation along and rotation about the vertical axis only. This dramatically inhibited surface bubble removal and associated dancing. Only a rare event of a large bubble leaving the surface could result in an excursion. Rather, if the body is free to rotate, then upon reaching the surface and losing a portion of the surface-bound bubbles, the body becomes unstable to rotation due to the remaining bubbles on the underbelly of the sphere. Once the body starts to rotate, the body is cleared of a larger number of buoyancy-conferring bubbles, and an excursion becomes far more likely (see Movie S4). The previously noted bouncing at the free surface is an additional consequence. Fluctuations from the active fluid, or other immersed bodies (Movie S5), which nudge the rotational instability to develop appear to be important as well. For bodies like raisins, whose trajectories are influenced more by individual bubble growth, lift, and removal, vertical dancing may be observed without need for such rotations. Mathematical model Gas escape There are a number of mechanisms by which gas may escape from the system: formation of bubbles on the container walls and on the immersed body, which eventually exits into the surrounding environment, and by diffusive transport through the fluid-air interface^13. For the first mechanism, assuming simply that the frequency of surface bubble growth and bubble volume upon exit each to be linear in the supersaturation ratio, the gas loss from this process follows \(\dot{S}=-q{(S-{S}_{mc})}^{2}\), where q is a constant which encapsulates the number of bubble growth sites and their geometrical features, and S[mc] is the minimum value below which bubbles no longer form on or near the container walls (discussed below). This quadratic law is supported by the data in Fig. 2, and it suggests that this mechanism of gas loss is dominant for the first one or two minutes of the experiment. Gas loss due to growth on the immersed body and direct delivery to the surface may be neglected (see Supplementary Information). The steady stream of bubbles rising from (near) the container walls drives a large-scale circulation flow on the scale of the container's lateral measure akin to intrinsic convection in particle sedimentation^39,42,43. Such a flow is expected to continually replenish the region near the fluid-air interface with the well-mixed concentration from the bulk fluid. Balancing advection and diffusion near the free surface, the concentration is predicted to decay from its volume-averaged value inside the fluid to approximately zero outside the fluid across a small boundary layer of size \(\delta=\ sqrt{DH/(2{U}_{b})}\), where H is the fluid depth and U[b] is the velocity of bubble rise (see Supplementary Information). This motivates a model for the supersaturation ratio, \(\dot{S}=-(S-{S}_ {mc})/{T}_{r}\), where T[r] = V[f]d/(DS[f]) is a relaxation timescale, with S[f] the free surface area and V[f] the fluid volume. Using the dimensions of the experiment, L = 8.9 cm and H = 4.5 cm, the diffusion constant for CO[2] in water D = 1.85 x 10^-5 cm^2/s^44, and an observed bubble rise velocity of U[b] [?] 1 cm/s, this gives d = 65 mm. Then with S[f] = 79 cm^2, and fluid volume V[f] = 355 cm^3, the predicted relaxation time T[r] is 26 min, within range of the best-fit value used in Fig. 2. Combining the mechanisms above, \(\dot {S}=-(S-{S}_{mc})/{T}_{r}-q{(S-{S}_{mc})}^{2}\), produces the expression in (1). Discrete and continuum buoyancy growth models We consider two models of bubble/buoyancy growth, a discrete model and a continuum model, as each can be more appropriate depending on the body size and surface properties. In the discrete model, each of N bubbles are assumed to grow independently according to the bubble growth law attributed to Scriven, which builds upon the Rayleigh-Plesset equation: for an isolated bubble of radius a(t), we have $$\dot{a}=\frac{D}{a}\left(S-\frac{2\sigma }{pa}\right),$$ (3) where D is a diffusion constant (D = 1.85 x 10^-5 cm^2/s for CO[2] in water^44), S(t) is the supersaturation ratio, s is the surface tension of water (s [?] 70 mN/m = 70 dyn/cm at room temperature), and p is the pressure near the bubble (p [?] 1 atm = 10^6 dyn/cm^2)^1,45,46, 47. The force conferred to the body is the buoyancy experienced by the bubble. In the discrete model, a bubble is removed from the surface when its position on the body exits the fluid, or when it reaches a critical size for pinch-off, a[p], at which time it is immediately rebirthed (once in the fluid) with size a[0], a characteristic scale of surface roughness. Of particular importance at longer times, there is a supersaturation ratio below which bubbles are overwhelmed by pressure and surface tension and cease to form on the body, S = 2s/pa[0]. The initial bubble size, a[0], depends upon the surface and the nature of bubble nucleation there. The length scale of surface roughness on the body may be associated with such initial bubble sizes; the lumen diameter of the fibers left behind during cleaning on the container walls is another. This value for the bubble formation near the surface was denoted S[mc], and inferred from the experiments in SSII to be S[mc] [?] 0.020. This suggests that the surface roughness (or remnant material) scale on the container wall is 2s/(pS[mc]) [?] 30 mm, which matches the diameter of cellulose fibers^5. On the body, meanwhile, using the resolution of the 3D printer of 0.15 mm to estimate a[0][?]0.015 cm, a minimal supersaturation ratio for bubble growth on the body, denoted by S[mb], may be roughly 0.010. That S[mb] < S[mc] indicates that bubbles should continue to form on the body after they cease to form on/near the container walls. Indeed, after a carbonated fluid was left alone for hours until it appeared 'flat' (bubble-free), inserting at that time a 3D-printed body or a raisin resulted in bubble growth on the body and the onset of body rising events. A detailed determination of S[mb] generally involves not only the surface roughness but also its chemistry, as the contact angle has been shown to be important for bubble nucleation^31. The second model considered is a continuum model, more appropriate when the body is covered in a large number of bubbles which are continually growing, merging, and detaching. With the local maximum of the added traction given by \({B}_{s}/(4\pi {A}^{2})\hat{{{{{{{{\ bf{z}}}}}}}}}\), the instantaneous portion of this local contribution is defined as \(b{B}_{s}/(4\pi {A}^{2})\hat{{{{{{{{\bf{z}}}}}}}}}\), where b(x, t) [?] [0, 1] and x is a point on the body surface. The evolution of this buoyancy fraction is then modeled as $$\frac{d}{dt}b({{{{{{{\bf{x}}}}}}}},t)=\frac{\lambda (t)}{{B}_{s}}\ left(1-b({{{{{{{\bf{x}}}}}}}},t)\right),$$ (4) which holds pointwise at every position x on the body surface. Spatial variations in b(x, t) are possible, and generically produce a body torque. Holding the supersaturation ratio fixed and neglecting surface tension and pressure, the bubble growth law suggests a bubble radius growth \(a(t) \sim {({a}_{0}^{2}+2DS({t}_{0})(t-{t}_{0}))}^{1/ 2}\). Since the contributed buoyancy is proportional to a(t)^3, this suggests a growth rate l(t) which scales as S^3/2 for appreciable times, consistent with the measured data. This leads to the model for the growth rate given in (2). In both models, the instantaneous lifting force and torque in the lab frame are written as \({B}_{s}{F}_{B}[b]\hat{{{{{{{{\bf{z}}}}}}}}}\) and AB[s]L[B][b], respectively. In the continuum model, \({F}_{B}[b]= \langle b\rangle :={(4\pi )}^{-1}{\int}_{{S}_{0}}b({{{{{{{\bf {X}}}}}}}},t)\,d{S}_{0}\), where dS[0] is the surface area element in the body frame, and \({{{{{{{{\bf{L}}}}}}}}}_{B}[b]=\left({{{{{{{\bf {Q}}}}}}}}\langle b{{{{{{{\bf{X}}}}}}}}\rangle \right)\times \hat {{{{{{{{\bf{z}}}}}}}}}\). Or, defining the center of surface buoyancy in the body frame as \({{{{{{{\bf{M}}}}}}}}=\langle b{{{{{{{\bf {X}}}}}}}}\rangle={(4\pi )}^{-1}{\int}_{{S}_{0}}{{{{{{{\bf{X}}}}}}}}b ({{{{{{{\bf{X}}}}}}}},t)\,d{S}_{0}\) and m = QM the same in the lab frame, we may write \({{{{{{{{\bf{L}}}}}}}}}_{B}[b]=({{{{{{{\bf {Q}}}}}}}}{{{{{{{\bf{M}}}}}}}})\times \hat{{{{{{{{\bf{z}}}}}}}}}\). Additional details are given in the Supplementary Information. Equations of motion and dimensionless groups Taking \(T:=\sqrt{A/g}\) to be a characteristic time and \(A/T=\sqrt {Ag}\) to be a characteristic speed, we write the position of the body centroid as \({{{{{{{\bf{r}}}}}}}}(t)=A\cdot Z(t)\hat{{{{{{{{\bf {z}}}}}}}}}\), the vertical velocity as (A/T)W, and the body rotation rate as T^-1O. Force and torque balance, with the dimensionless time s [?] t/T, are expressed as $$\frac{dW}{ds}=-1+\frac{\alpha (s)}{{{{{{{{\mathcal{M}}}}}}}}}+\beta {F}_{B}[b]-\frac{9{C}_{T}}{2{{{{{{{\mathcal{M}}}}}}}}{{{{{{{\rm {Re}}}}}}}}}W,$$ (5) $$\frac{d}{ds}\left({I}_{R}{{{{{{{\mathbf{\Omega }}}}}}}}\right)=\ beta {{{{{{{{\bf{L}}}}}}}}}_{B}[b]-\frac{6{C}_{R}}{{{{{{{{\mathcal {M}}}}}}}}{{{{{{{\rm{Re}}}}}}}}}{{{{{{{\mathbf{\Omega }}}}}}}},$$ (6) where we have introduced the following dimensionless numbers, $${{{{{{{\mathcal{M}}}}}}}}=\frac{m}{\rho V},\,\,\,\beta=\frac{{B}_ {s}}{mg},\,\,\,\Lambda=\frac{{(A/g)}^{1/2}{\lambda }_{0}}{{B}_{s}},\, \,\,\tau=\frac{{T}_{r}}{{(A/g)}^{1/2}}.$$ (7) Here, V is the body volume, V[S](s) is the submerged volume at time s , a(s) [?] V[S](s)/V [?] [0, 1], and mA^2I[R] is the moment of inertia, where I[R] = 2/5 for a rigid sphere. Hydrodynamic drag and torque coefficients, respectively, are given by \({C}_{T}=1+.0183\,{{{{{{{\ rm{Re}}}}}}}}\,| W|\), and \({C}_{R}=1+0.0044\sqrt{{{{{{{{\rm {Re}}}}}}}}\,| {{{{{{{\boldsymbol{\Omega }}}}}}}}| }\), with \ ({{{{{{{\rm{Re}}}}}}}}=\rho {A}^{3/2}{g}^{1/2}/\mu\) the Reynolds number (The Reynolds number is defined as \({{{{{{{\rm{Re}}}}}}}}=\ rho UA/\mu\), with r the fluid density, U a characteristic speed, A a characteristic length scale, and m the fluid viscosity, and gives a measure of the importance of inertia relative to viscous dissipation.). Bubbles generally affect the drag on the body but we neglect this detail here. The surface buoyancy fraction b evolves at each point on the body as $$\frac{d}{ds}b({{{{{{{\bf{x}}}}}}}},s)=\Lambda g(s)\left(1-b ({{{{{{{\bf{x}}}}}}}},s)\right),$$ (8) where \(g(s)={\left(S(Ts)/{S}_{0}\right)}^{3/2}\), and using Eqn. (1 ), $$g(s)={\left(\frac{{S}_{mc}}{{S}_{0}}+\frac{(1-{S}_{mc}/{S}_{0})\exp (-s/\tau )}{1+\chi (1-\exp (-s/\tau ))}\right)}^{3/2}.$$ (9) Any points x on the body surface which are outside of the fluid are given the value b = 0; only once those points reenter the fluid does the growth there begin again according to (8). The system is closed by tracking the body's position and orientation: \(\dot{Z}=W\) and \ (\dot{{{{{{{{\bf{Q}}}}}}}}}={{{{{{{\bf{Q}}}}}}}}\,\hat{\Omega }\) where dots denote derivatives upon the dimensionless time s, and \(\ hat{\Omega }{{{{{{{\bf{q}}}}}}}}:={{{{{{{\mathbf{\Omega }}}}}}}}\ times {{{{{{{\bf{q}}}}}}}}\) for any vector q. Initial conditions are generally taken to be Z(0) = - 1, W(0) = 0, O(0) = 0, Q(0) = I, and b (x, 0) = 0. The system is thus characterized by a mass ratio, \({{{{{{{\mathcal {M}}}}}}}}\), the relative lifting force b, which we term the fizzy lifting number ("Fizzy lifting drinks! They fill you with bubbles, and the bubbles are full of a special kind of gas, and this gas is so terrifically lifting that it lifts you right off the ground just like a balloon, and up you go until you're bumping against the ceiling!" -Charlie and the Chocolate Factory, by Roald Dahl.), initial bubble growth rate, L, and relaxation time, t, the minimal supersaturation ratio for bubble growth along the container, S[mc], the initial supersaturation ratio, S[0], and the Reynolds number. For the 3D-printed body used above, using A = 1 cm, m = 4.25 g, B[s] = 700 dyn, l[0] = 58.4 dyn/s, and T[r] = 36.2 min, and for water r = 1 g/cm^3 and m = 0.01 g/(cm s), we find \(({{{{{{{\mathcal {M}}}}}}}},\, \beta,\, \Lambda,\, \tau,\, {{{{{{{\rm{Re}}}}}}}})= (1.015,\, 0.17,\, 2.6\cdot 1{0}^{-3}, 6.8\cdot 1{0}^{4},\, 3.1\cdot 1 {0}^{3})\). For a raisin, using a prolate spheroidal body with semi-major axis length A = 0.6 cm and semi-minor axis lengths 0.4 cm, mass m = 0.45 g, B[s] = 100 dyn, and l[0] = 20 dyn/s, we find \ (({{{{{{{\mathcal{M}}}}}}}},\, \beta,\, \Lambda,\, {{{{{{{\rm {Re}}}}}}}})=(1.12,\, 0.23,\, 4.9\cdot 1{0}^{-3},\, 1.5\cdot 1{0}^ {3})\). These values of b and L make raisins particularly strong dancers. The body is positively buoyant and floats without bubbles if \ ({{{{{{{\mathcal{M}}}}}}}} < 1\). If \({{{{{{{\mathcal{M}}}}}}}} > 1 \), the body can only be lifted upward against gravity if the fizzy lifting number b is sufficiently large; namely \(\beta+{{{{{{{{\ mathcal{M}}}}}}}}}^{-1}-1\) must be positive. This gives a range of masses for which oscillatory dynamics are expected to reside: \ ({{{{{{{\mathcal{M}}}}}}}}\in \left(1,1/(1-\beta )\right)\) if b < 1, and \({{{{{{{\mathcal{M}}}}}}}}\in (1,\infty )\) if b >= 1. The 3D-printed bodies are predicted to dance with density ratios r[s]/r [?] (1, 1.20); raisins are expected to dance in a similar range, r[s]/r [?] (1, 1.29). With B[s] generally scaling with the surface area and m scaling with the body volume, b is generically larger for very small bodies, and bubble-induced lifting is expected for such bodies to be more immediate. Dancing frequency We approximate the solution to (8) by neglecting the early period of rapid gas escape and assuming S[mc] [?] S[0]. Taking \(g(s)\approx {\ left(\exp (-s/\tau )/(1+\chi )\right)}^{3/2}\), with b(x, s[0]) = b [0], we find $$b({{{{{{{\bf{x}}}}}}}},s)=1-(1-{b}_{0})\exp \left(\frac{2\Lambda \ tau }{3{(1+\chi )}^{3/2}}\left[{e}^{-3s/2\tau }-{e}^{-3{s}_{0}/2\tau }\right]\right).$$ (10) The dimensionless 'charging time' before the total buoyancy overcomes gravity is given by s[charge] [?] s[c] - s[0] such that \(\beta \langle b\rangle ({s}_{c})+1/{{{{{{{\mathcal{M}}}}}}}}-1=0\), which from (10) is given by $${s}_{charge}=\frac{2\tau }{3}\log \left(\frac{1}{{e}^{-3{s}_{0}/2\ tau }-\frac{3{(1+\chi )}^{3/2}}{2\Lambda \tau }{L}_{\beta }}\right)- {s}_{0},$$ (11) where \({L}_{\beta }=\log \left(\frac{\beta {{{{{{{\mathcal{M}}}}}}}} (1-{b}_{0})}{1-(1-\beta ){{{{{{{\mathcal{M}}}}}}}}}\right)\). For a sufficiently small container the charging time serves as a proxy for the excursion time, with transit from one surface to another playing only a small role. Two curves corresponding to f = 1/D [?] 1/t[charge] (with t[charge] [?] Ts[charge]) are included in Fig. 3e, one with b[0] = 0 and one which best fits the data with b[0] = 0.073. The frequency is sensitive to b[0], which points yet again to the importance of rotations, and how many bubbles are removed upon each surface visit. For insertion times s[0] > s[fun], where \({s}_{fun}=(2\tau /3)\log (2\Lambda \tau /(3(1+\chi ){L}_{\beta }))\), the charging time in (11 ) is infinite. Using the values from the experiments the associated dimensionless time is t[fun] [?] Ts[fun] = 170 min. A separate approximation starting from (8) is more appropriate at long times, assuming S[mc] > S[mb] (i.e., bubbles continue to form on the body after they cease to form along the container walls). In this case, as s[0]/t - [?], we have \(g({s}_{0}) \sim {({S}_{mc}/{S}_{0})}^{3/2}=:{g} _{\infty }\), leading to a constant charging time of \({s}_{charge}= {(\Lambda {g}_{\infty })}^{-1}{L}_{\beta }\). Using the experimental parameters this gives a final dimensional charging time of between 14 min and 2.5 min, for initial coverages b[0] = 0 and b[0] = 0.073, respectively. Once S <= S[mc], bubbles no longer form along the container and the primary mechanism driving gas escape is removed. The body can continue to form bubbles and perform its low-frequency dance, even in an otherwise quiescent fluid. Since we neglect the days-long timescale of pure diffusive transport, and convective diffusion affected by the body motion, this low-frequency dancing is predicted to carry on indefinitely. Dancing with a mean frequency of 1.5 min^-1 was indeed observed in one experimental run for the last 4 h of a 5 h run. Simulations To examine the dynamics in a more controlled environment we solve (5- 8) numerically. Figure 5 shows the dynamics of a body with bubbles at the vertices of a regular icosahedron (N = 12), using the discrete bubble model, with \(({{{{{{{\mathcal{M}}}}}}}},\, \beta,\, \Lambda, \, {{{{{{{\rm{Re}}}}}}}})=(1.02,\, 0.17,\, 0.003,\, 3100)\). The rebirth and pinch-off radii used are a[0]/A = 0.03 and a[p]/A = 0.24, respectively. The bubbles first grow to sufficient stature to lift the body to the surface, where the topmost bubbles are released. The body begins to descend slowly before undergoing a rotation, returning soon after to the surface. Each visit to the surface includes bubble removal, a body rotation, additional bubble removal, and then a plummeting to deeper waters. Fig. 5: Simulations using the discrete bubble model. figure 5 A body with \(({{{{{{{\mathcal{M}}}}}}}},\beta,\Lambda,{{{{{{{\rm {Re}}}}}}}})=(1.02,0.17,0.003,3100)\), and bubbles at the vertices of an inscribed icosahedron (N = 12), undergoes "bouncing" dynamics near the surface. After the bubbles nearest to the surface are removed, the center of surface buoyancy rests below the center of mass, resulting in a torque and eventual body rotation. Only after a few returns to the surface to clear off more bubbles does the body begin a large excursion back towards the container floor. See Movie S7. Full size image Figure 6a shows the dimensionless vertical dancing frequency, Tf, using the discrete bubble model as a function of the number of bubbles, all placed at the vertices of a regular polyhedron, all with \(({{{{{{{\mathcal{M}}}}}}}},\, \beta,\, \Lambda )=(1.015,\, 0.17,\, 0.0016)\). To isolate the role of bubble position we adjust the maximum bubble size \({a}_{p}/A={(3{B}_{s}/(4\pi \rho gN))}^{1/3}\) so that each body achieves the same maximal surface buoyancy (b) if all bubbles are at their maximal (pinch-off) size. The frequency increases with N, as bodies with a greater number of bubbles can begin their descent while still maintaining partial surface coverage, and less bubble growth is needed before the body becomes positively buoyant again. Fig. 6: Numerical simulations using the discrete and continuum models. figure 6 a The (dimensionless) dancing frequency, which increases with the number of bubbles, N, for fixed maximal surface lifting force (b = 0.017). b The mean rotation rate increases with the growth rate, L, in the continuum model (solid line). For the discrete model with N = 12 bubbles (circles), the mean rotation rate increases with L until premature bubble detachment becomes important. c The mean rotation rate as a function of b, for a selection of growth rates L, using the discrete model with N = 12 bubbles. Rotations are damped at large b and L due to the body emerging from the fluid. d Same as (c) but using the continuum model. Full size image To further explore body rotations, Fig. 6b shows the mean body rotation rate as a function of the growth rate, L, using both the discrete bubble model (shown as symbols) and the continuum model (as a solid curve). With increasing L the body spends more time at the surface, and experiences opportunities to rotate more frequently, if not more rapidly. The growth is approximately logarithmic in L, and for small L there is close agreement between the discrete and continuum models. The discontinuity in the discrete model is due to the onset of premature bubble detachment at large bubble growth rates (see Movie S7). For sufficiently large growth rates (L > 0.008) the bubbles grow from their initial size a[0] to the pinch-off size a[p] on a shorter timescale than the body's excursion time. Consequently, most, if not all, bubbles are removed near the time that the body reaches the surface, nearly eliminating the torques on the body through pinch-off alone and thus dampening rotations. Figure 6c shows the mean rotation rate instead as a function of b for a selection of growth rates, L, using the discrete bubble model, while Fig. 6d shows the same using the continuum model. Generally, larger values of b are associated with larger torques, and thus faster body spinning. The discrete model shows again the importance of premature bubble pinch-off and departure at large bubble growth rates. For L = 0.016, if b is small, it is common for most of the bubbles to pinch off before the body traverses the full length of the container. For large b, the body emerges completely out of the fluid in a dramatic jump, and all bubbles are removed leaving none to generate a torque. Wobbling and rolling For bodies which are large relative to the maximum bubble size, body rotations are commonly observed, as are a related dynamics: wobbling. As a coarse approximation we consider bubbles to have been removed from one half of the spherical surface, which then regrow with rate L. Writing the center of surface buoyancy in the body frame as \ ({{{{{{{\bf{M}}}}}}}}={{{{{{{\bf{M}}}}}}}}(0)\exp (-\Lambda s)\), with \({{{{{{{\bf{M}}}}}}}}(0)={(4\pi )}^{-1}\langle b{{{{{{{\bf {X}}}}}}}}\rangle=(1/4)\hat{{{{{{{{\bf{z}}}}}}}}}\), then with \ ({{{{{{{{\bf{q}}}}}}}}}_{1}=\cos (\theta )\hat{{{{{{{{\bf{x}}}}}}}}}+ \sin (\theta )\hat{{{{{{{{\bf{z}}}}}}}}}\) and \({{{{{{{{\bf {q}}}}}}}}}_{2}=\hat{{{{{{{{\bf{y}}}}}}}}}\), the (dimensionless) torque in the lab frame is \({{{{{{{{\bf{L}}}}}}}}}_{B}=-(1/4)\exp (- \Lambda s)\sin (\theta )\hat{{{{{{{{\bf{z}}}}}}}}}\). Writing \ ({{{{{{{\boldsymbol{\Omega }}}}}}}}=\dot{\theta }\hat{{{{{{{{\bf {y}}}}}}}}}\), and neglecting the nonlinear part of the hydrodynamic moment, we arrive at \(\ddot{\theta }=-{a}_{1}(s)\sin (\theta )-{a}_ {2}\dot{\theta }\), where \({a}_{1}(s)=5\beta \exp (-\Lambda s)/8\) and \({a}_{2}=15/({{{{{{{\mathcal{M}}}}}}}}\,{{{{{{{\rm{Re}}}}}}}})\) , the equation for a damped nonlinear oscillator with diminishing torque. Wobbling diminishes either by viscous damping or by bubble (re) growth. A characteristic initial (s = 0) wobbling frequency from the above is \({f}_{wobble}\approx \sqrt{{a}_{1}}/(2\pi )\, \approx \, {\ left(5\beta /8\right)}^{1/2}/(2\pi )\). For the 3D-printed body, using b = 0.17, the initial frequency is roughly f[wobble]/T [?] 1.6 Hz. A raisin, meanwhile, due to its larger value of b, oscillates with a higher frequency of just over 2.4 Hz. These values are consistent with the experimental observations (see Movies S1, S4, and S5). This effect is similar to that seen in Quincke rotor dynamics, where surface charging is driven by electrohydrodynamics^48,49,50,51, 52,53,54. A transient rolling mode was also observed. Bubbles on the surface of a rolling body begin to grow upon reentry and are larger just before they exit, producing a sustained rolling torque. For this we consider a cylindrical body and a two-dimensional cross-section. If the body is fixed at a vertical position \(Z=-\cos ({\theta }^{*})\), with th^ * [?] [-p, p], then the steady state lifting distribution assuming \(\ dot{b}=\Lambda\) (and using \(\dot{b}=\Omega {b}_{\theta }\)) is b(th) = LO^-1(th - th^*), where \(\Omega=\dot{\theta }=| {{{{{{{\boldsymbol {\Omega }}}}}}}}|\). The resulting dimensionless torque is \({{{{{{{\ bf{m}}}}}}}}=2\Lambda {\Omega }^{-1}{(\pi -{\theta }^{*})}^{2}\hat {{{{{{{{\bf{y}}}}}}}}}\). Balancing with a viscous drag \(-\eta {({{{{{{{\mathcal{M}}}}}}}}{{{{{{{\rm{Re}}}}}}}})}^{-1}\Omega \hat {{{{{{{{\bf{y}}}}}}}}}\) yields a steady rotation rate \(\Omega={(2 {(\pi -{\theta }^{*})}^{2}\Lambda {{{{{{{\mathcal{M}}}}}}}}{{{{{{{\rm {Re}}}}}}}}/\eta )}^{1/2}\), where e <= 15, since part of the body sits outside of the water. A more detailed study of this rotational drag, like that performed by Hunt et al.^55, is needed. Discussion Supersaturated fluids present an accessible playground for exploring the dynamics of bodies and their relationship to a complex fluid environment. In the framework proposed by Spagnolie & Underhill^56 this would appear as either a Type I or Type II system - the body is much larger than the 'obstacles' (be they bubbles or gas molecules), and the fluid exhibits a natural relaxation time. Among the unexpected findings in this system, we have observed a critical dependence of the dynamics on body rotations for large body-to-bubble size ratios, and multi-period oscillatory dynamics when a single surface interaction is insufficient to clean the body surface of its lifting agents. Another intriguing feature is that bubbles can continue to form on the body long after they cease to form at the container walls when S[mc] > S[mb] (when the surface roughness or fibrous material scale is smaller on the container walls than on the body). Raisins inserted into a fluid which was left out for hours and appeared motionless were indeed observed to dance, albeit at a leisurely pace. Relatedly, Pereira et al.^31 identified that a smaller contact angle between a bubble and a surface decreases the energy needed to form bubbles there. Preliminary results also suggest that the body's presence can affect fluid degassing, even potentially slowing it by disturbing the large-scale convective flow responsible for gas escape, reminiscent of how moving boundaries affect flows which promote heat transport^57,58,59,60. Additional constraints on the dynamics are expected in general. At greater depths, bubbles are less likely to grow due to the increase in hydrostatic pressure. Should a body plummet sufficiently far below the surface, bubble-assisted levitation may vanish and the plummeting will continue unresisted. Surfactants may also adjust the range of bodies possible to levitate in this manner, since their presence can affect the nature of bubble pinch-off and coalescence, in competition with the pressure and surface roughness scale via S[mb]^61,62. The shape of the container and temperature can also affect the rate of CO [2] loss^63. Another uncharacterized but potentially important feature is wetting. A return to the surface releases bubbles from the body surface, but interaction with the air above may also help to nucleate other bubbles by drawing additional gas into small cavities. Initially, dry bodies danced for far longer than initially wet bodies. The interaction with the free surface appears to damp rotations as well, and as we have seen, any inhibition of the rotation of large bodies tends to inhibit vertical dancing. A number of directions lie ahead based on additional observations not described above. Preliminary studies suggest a substantial encouragement of excursions for more elongated and asymmetric body shapes. Also, the behavior of multiple bodies in the system can result in stable rafts of bubble-sharing bodies at the surface. Since rotation is critical for triggering an excursion from the surface, and sharing bubbles inhibits rotation, the system transitions from exhibiting oscillatory to overdamped behavior with increased particle volume fraction. Additional bodies, however, also increase the fluctuations in the system, which can encourage the rolling and plummeting of others. This instability to rotation and cooperative effects are reminiscent of iceberg capsize dynamics^64,65,66. Exploration of the optimal number of dancing partners is under current investigation. A more detailed study of the fluid flow itself will be highly informative regarding the nature of degassing and the effect of the immersed body. Theoretical advances are also needed. Models of growing arrays of bubbles date back to classical works by Lifschitz and Slyozov^67 and Wagner^68. Transient coarsening kinetics depend on numerous simultaneous mass transfer mechanisms, resulting in overlapping scaling behaviors in time^69,70,71,72. Collective formation and dissolution of bubbles on a regular patterned grid have recently provided some insight on these coupled effects^73,74,75. The additional presence of a flowing environment presents a substantial new challenge. Future work exploring body shape, multi-body dynamics, and fluid-structure interactions is likely to prove... fruitful. Data availability Datasets for force growth on fixed bodies, and vertical positions during vertical oscillations, are available at https://doi.org/ 10.6084/m9.figshare.25302136. References 1. Scriven, L. E. On the dynamics of phase growth. Chem. Eng. Sci. 10, 1 (1959). Article CAS Google Scholar 2. Barker, G. S., Jefferson, B. & Judd, S. J. The control of bubble size in carbonated beverages. Chem. Eng. Sci. 57, 565 (2002). Article CAS Google Scholar 3. Zhang, Y. & Xu, Z. "Fizzics" of bubble growth in beer and champagne. Elements 4, 47 (2008). Article ADS CAS Google Scholar 4. Liger-Belair, G., Sternenberg, F., Brunner, S., Robillard, B. & Cilindre, C. Bubble dynamics in various commercial sparkling bottled waters. J. Food Eng. 163, 60 (2015). Article CAS Google Scholar 5. Liger-Belair, G., Voisin, C. & Jeandet, P. Modeling nonclassical heterogeneous bubble nucleation from cellulose fibers: application to bubbling in carbonated beverages. J. Phys. Chem. B 109, 14573 (2005). Article CAS PubMed Google Scholar 6. Uzel, S., Chappell, M. A. & Payne, S. J. Modeling the cycles of growth and detachment of bubbles in carbonated beverages. J. Phys. Chem. B 110, 7579 (2006). Article CAS PubMed Google Scholar 7. Soto, A. M., Maddalena, T., Fraters, A., Van Der Meer, D. & Lohse, D. Coalescence of diffusively growing gas bubbles. J. Fluid Mech. 846, 143 (2018). Article ADS MathSciNet CAS Google Scholar 8. Fardin, M. A., Hautefeuille, M. & Sharma, V. Spreading, pinching, and coalescence: the Ohnesorge units. Soft Matter 18, 3291-3303 (2022). Article ADS CAS PubMed Google Scholar 9. Lv, P. et al. Self-propelled detachment upon coalescence of surface bubbles. Phys. Rev. Lett. 127, 235501 (2021). Article ADS CAS PubMed Google Scholar 10. Iwata, R. et al. How coalescing bubbles depart from a wall. Langmuir 38, 4371 (2022). Article CAS PubMed Google Scholar 11. Zhao, P., Hu, Z., Cheng, P., Huang, R. & Gong, S. Coalescence-induced bubble departure: effects of dynamic contact angles. Langmuir 38, 10558 (2022). Article CAS PubMed Google Scholar 12. Hooshanginejad, A. et al. Effect of angle in removing proteins or bacteria on a tilted surface using air bubbles. Phys. Rev. Fluids 8, 043602 (2023). Article ADS Google Scholar 13. Liger-Belair, G. The physics behind the fizz in champagne and sparkling wines. Euro. Phys. J. Spec. Top. 201, 1 (2012). Article ADS CAS Google Scholar 14. Lohse, D. Bubble puzzles: from fundamentals to applications. Phys. Rev. Fluids 3, 110504 (2018). Article ADS Google Scholar 15. Sparks, R. S. J. The dynamics of bubble formation and growth in magmas: a review and analysis. J. Volcanol. Geotherm. Res. 3, 1 (1978). Article ADS CAS Google Scholar 16. Gilbert, J. S. and Sparks, R. S. J. The Physics of Explosive Volcanic Eruptions (Geological Society of London, 1998). 17. Lensky, N. G., Navon, O. & Lyakhovsky, V. Bubble growth during decompression of magma: experimental and theoretical investigation. J. Volcanol. Geotherm. Res. 129, 7 (2004). Article ADS CAS Google Scholar 18. Goncalves, M., Vasconcelos, A. M. P., Gomes de Azevedo, E. J. S., Chaves das Neves, H. J. & Nunes da Ponte, M. On the application of supercritical fluid extraction to the deacidification of olive oils. J. Am. Oil Chem. Soc. 68, 474 (1991). Article Google Scholar 19. Markom, M., Singh, H. & Hasan, M. Supercritical co2 fractionation of crude palm oil. J. Supercrit. Fluids 20, 45 (2001). Article CAS Google Scholar 20. Harvey, E. N. et al. Bubble formation in animals. I. Physical factors. J. Cell. Comp. Physiol. 24, 1 (1944). Article CAS Google Scholar 21. Papadopoulou, V., Tang, M.-X., Balestra, C., Eckersley, R. J. & Karapantsios, T. D. Circulatory bubble dynamics: from physical to biological aspects. Adv. Colloid Interface Sci. 206, 239 (2014). Article CAS PubMed Google Scholar 22. Manga, M. Waves of bubbles in basaltic magmas and lavas. J. Geophys. Res. Solid Earth 101, 17457 (1996). Article Google Scholar 23. Robinson, M., Fowler, A. C., Alexander, A. J. & O'Brien, S. B. G. Waves in Guinness. Phys. Fluids 20, 067101 (2008). Article ADS Google Scholar 24. Watamura, T. et al. Bubble cascade in Guinness beer is caused by gravity current instability. Sci. Rep. 9, 5718 (2019). Article ADS PubMed PubMed Central Google Scholar 25. Planinsic, G. Fizziology. Phys. Educ. 39, 65 (2004). Article ADS Google Scholar 26. Cordry, S. M. Finicky clay divers. Phys. Teach. 36, 82 (1998). Article ADS Google Scholar 27. Derr, B. J., Derr, H. R. & Lewis, T. Gas me up, or, A baking powder diver. J. Chem. Educ. 77, 171 (2000). Article CAS Google Scholar 28. Moinester, M., Gerland, L., Liger-Belair, G. & Ocherashvili, A. Fizz-ball fizzics. Phys. Teacher 50, 284 (2012). Article ADS Google Scholar 29. Maselko, J. et al. Emergence of complex behavior in chemical cells: the system AlCl3-NaOH. Langmuir 30, 5726 (2014). Article CAS PubMed Google Scholar 30. Wang, Q. & Steinbock, O. Bobbing chemical garden tubes: oscillatory self-motion from buoyancy and catalytic gas production. Soft Matter 19, 2138 (2023). Article ADS CAS PubMed Google Scholar 31. Pereira, L. et al. The physics of dancing peanuts in beer. Roy. Soc. Open Sci. 10, 230376 (2023). Article ADS CAS Google Scholar 32. Zenit, R. & Rodriguez-Rodriguez, J. The fluid mechanics of bubbly drinks. Phys. Today 71, 44 (2018). Article Google Scholar 33. Mathijssen, A. J. T. M., Lisicki, M., Prakash, V. N. & Mossige, E. J. L. Culinary fluid mechanics and other currents in food science. Rev. Mod. Phys. 95, 025004 (2023). Article ADS MathSciNet CAS Google Scholar 34. Hisatake, K., Tanaka, S. & Aizawa, Y. Evaporation rate of water in a vessel. J. Appl. Phys. 73, 7395 (1993). Article ADS CAS Google Scholar 35. Lubetkin, S. & Blackwell, M. The nucleation of bubbles in supersaturated solutions. J. Colloid Inter. Sci. 126, 610 (1988). Article ADS CAS Google Scholar 36. Carroll, J. J., Slupsky, J. D. & Mather, A. E. The solubility of carbon dioxide in water at low pressure. J. Phys. Chem. Ref. Data 20, 1201 (1991). Article ADS CAS Google Scholar 37. Sander, R. Compilation of henry's law constants (version 4.0) for water as solvent. Atmos. Chem. Phys. 15, 4399 (2015). Article ADS CAS Google Scholar 38. Ashurst, P. R. Chemistry and Technology of Soft Drinks and Fruit Juices (John Wiley & Sons, 2016). 39. Liger-Belair, G. Effervescence in champagne and sparkling wines: from grape harvest to bubble rise. Euro. Phys. J. Spec. Top. 226, 3 (2017). Article ADS CAS Google Scholar 40. Fritz, W. Berechnung des maximalvolumes von dampfblasen. Physik. Z. 36, 379 (1935). Google Scholar 41. Oguz, H. N. & Prosperetti, A. Dynamics of bubble growth and detachment from a needle. J. Fluid Mech. 257, 111 (1993). Article ADS CAS Google Scholar 42. Peysson, Y. & Guazzelli, E. An experimental investigation of the intrinsic convection in a sedimenting suspension. Phys. Fluids 10 , 44 (1998). Article ADS CAS Google Scholar 43. Guazzelli, E. and Morris, J. F. A Physical Introduction to Suspension Dynamics, Vol. 45 (Cambridge Univ. Press, 2011). 44. Liger-Belair, G., Prost, E., Parmentier, M., Jeandet, P. & Nuzillard, J.-M. Diffusion coefficient of CO2 molecules as determined by 13C NMR in various carbonated beverages. J. Agric. Food Chem. 51, 7560 (2003). Article CAS PubMed Google Scholar 45. Rayleigh, L. Viii. on the pressure developed in a liquid during the collapse of a spherical cavity. Lond. Edinb. Dublin Philos. Mag. J. Sci. 34, 94 (1917). Article Google Scholar 46. Plesset, M. S. The dynamics of cavitation bubbles. J. Appl. Mech. 16, 277 (1949). Article ADS Google Scholar 47. Epstein, P. S. & Plesset, M. S. On the stability of gas bubbles in liquid-gas solutions. J. Chem. Phys. 18, 1505 (1950). Article ADS CAS Google Scholar 48. Quincke, G. Ueber Rotationen im constanten electrischen Felde. Ann. Phys. Chem. 59, 417 (1896). Article ADS Google Scholar 49. Jones, T. B. Quincke rotation of spheres. IEEE Trans. Ind. Appl. 4, 845 (1984). Article Google Scholar 50. Turcu, I. Electric field induced rotation of spheres. J. Phys. A 20, 3301 (1987). Article ADS Google Scholar 51. Lemaire, E. & Lobry, L. Chaotic behavior in electro-rotation. Phys. A 314, 663 (2002). Article Google Scholar 52. Peters, F., Lobry, L. & Lemaire, E. Experimental observation of Lorenz chaos in the Quincke rotor dynamics. Chaos 15, 013102 (2005). 53. Vlahovska, P. M. Electrohydrodynamics of drops and vesicles. Annu. Rev. Fluid Mech. 51, 305 (2019). Article ADS MathSciNet Google Scholar 54. Zhang, Z., Yuan, H., Dou, Y., De La Cruz, M. O. & Bishop, K. J. Quincke oscillations of colloids at planar electrodes. Phys. Rev. Lett. 126, 258001 (2021). Article ADS CAS PubMed Google Scholar 55. Hunt, R. et al. Drag on a partially immersed sphere at the capillary scale. Phys. Rev. Fluids 8, 084003 (2023). Article ADS Google Scholar 56. Spagnolie, S. E. & Underhill, P. T. Swimming in complex fluids. Annu. Rev. Condens. Mat. Phys. 14, 381 (2023). Article ADS Google Scholar 57. Zhong, J.-Q. & Zhang, J. Thermal convection with a freely moving top boundary. Phys. Fluids 17, 115105 (2005). 58. Liu, B. et al. Self-induced cyclic reorganization of free bodies through thermal convection. Phys. Rev. Lett. 100, 244501 (2008). Article ADS PubMed Google Scholar 59. Lappa, M. Thermal Convection: Patterns, Evolution and Stability (John Wiley & Sons, 2009). 60. Ahlers, G., Grossmann, S. & Lohse, D. Heat transfer and large scale dynamics in turbulent Rayleigh-Benard convection. Rev. Mod. Phys. 81, 503 (2009). Article ADS CAS Google Scholar 61. King, L. & Sadhal, S. S. Effect of surfactants on the growth and departure of bubbles from solid surfaces. Heat Mass Transf. 50, 373 (2014). Article ADS CAS Google Scholar 62. Raza, M. Q., Kumar, N. & Raj, R. Surfactants for bubble removal against buoyancy. Sci. Rep. 6, 19113 (2016). Article ADS CAS PubMed PubMed Central Google Scholar 63. Liger-Belair, G., Bourget, M., Pron, H., Polidori, G. & Cilindre, C. Monitoring gaseous CO2 and ethanol above champagne glasses: flute versus coupe, and the role of temperature. PLoS ONE 7, e30628 (2012). Article ADS CAS PubMed PubMed Central Google Scholar 64. Macayeal, D. R., Scambos, T. A., Hulbe, C. L. & Fahnestock, M. A. Catastrophic ice-shelf break-up by an ice-shelf-fragment-capsize mechanism. J. Glaciol. 49, 22 (2003). Article ADS Google Scholar 65. Burton, J. C. et al. Laboratory investigations of iceberg capsize dynamics, energy dissipation and tsunamigenesis. J. Geophys. Res. Earth Surf. 117, F01007 (2012). 66. Bonnet, P. et al. Modelling capsizing icebergs in the open ocean. Geophys. J. Int. 223, 1265 (2020). Article ADS CAS Google Scholar 67. Lifshitz, I. M. & Slyozov, V. V. The kinetics of precipitation from supersaturated solid solutions. J. Phys. Chem. Solids 19, 35 (1961). Article ADS Google Scholar 68. Wagner, C. Theorie der alterung von niederschlagen durchumlosen, Zeitschrift Zeitschr. Elektrochemie 655, 81 (1961). Google Scholar 69. Slezov, V. V., Sagalovich, V. V. & Tanatarov, L. V. Theory of diffusive decomposition of supersaturated solid solution under the condition of simultaneous operation of several mass-transfer mechanisms. Phys. Chem. Solids 39, 705 (1978). Article ADS CAS Google Scholar 70. Ratke, L. Simultaneous coarsening of dispersions by growth and coagulation. Colloid Interface Sci. 119, 391 (1987). Article ADS CAS Google Scholar 71. Alexandrov, D. V. On the theory of Ostwald ripening in the presence of different mass transfer mechanisms. Phys. Chem. Solids 91, 48 (2016). Article ADS CAS Google Scholar 72. Li, Y. et al. Revisiting transient coarsening kinetics: a new framework in the Lifshitz-Slyozov-Wagner space. Acta Mater. 237, 118196 (2022). Article CAS Google Scholar 73. Peng, S., Mega, T. L. & Zhang, X. Collective effects in microbubble growth by solvent exchange. Langmuir 32, 11265 (2016). Article CAS PubMed Google Scholar 74. Zhu, X., Verzicco, R., Zhang, X. & Lohse, D. Diffusive interaction of multiple surface nanobubbles: shrinkage, growth, and coarsening. Soft Matter 14, 2006 (2018). Article ADS CAS PubMed Google Scholar 75. Michelin, S., Guerin, E. & Lauga, E. Collective dissolution of microbubbles. Phys. Rev. Fluids 3, 043601 (2018). Article ADS Google Scholar Download references Acknowledgements Support for this research was provided by the Office of the Vice Chancellor for Research and Graduate Education with funding from the Wisconsin Alumni Research Foundation, and by donations to the AMEP program (Applied Math, Engineering, and Physics), at the University of Wisconsin-Madison. S.E.S. gratefully acknowledges conversations with Hongyi Huang, Thomas G. J. Chandler, and Jean-Luc Thiffeault, and preliminary work with Carina Spagnolie. Author information Authors and Affiliations 1. Department of Mathematics, University of Wisconsin-Madison, Madison, WI, 53706, USA Saverio E. Spagnolie, Samuel Christianson & Carsen Grote 2. Department of Chemical & Biological Engineering, University of Wisconsin-Madison, Madison, WI, 53706, USA Saverio E. Spagnolie Authors 1. Saverio E. Spagnolie View author publications You can also search for this author in PubMed Google Scholar 2. Samuel Christianson View author publications You can also search for this author in PubMed Google Scholar 3. Carsen Grote View author publications You can also search for this author in PubMed Google Scholar Contributions S.E.S. conceived and coordinated the research. S.E.S. and S.C. performed physical experiments, S.E.S. and C.G. performed numerical experiments. The manuscript was prepared by S.E.S. Corresponding author Correspondence to Saverio E. Spagnolie. Ethics declarations Competing interests The authors declare no competing interests. Peer review Peer review information Nature Communications thanks Thodoris Karapantsios, Gerard Liger-Belair, and the other, anonymous, reviewer for their contribution to the peer review of this work. A peer review file is available. Additional information Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary information Supplementary Information Peer Review File Description of Supplementary Files Movie S1 Movie S2 Movie S3 Movie S4 Movie S5 Movie S6 Movie S7 Rights and permissions Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions About this article Check for updates. Verify currency and authenticity via CrossMark Cite this article Spagnolie, S.E., Christianson, S. & Grote, C. Levitation and dynamics of bodies in supersaturated fluids. Nat Commun 15, 3910 (2024). https://doi.org/10.1038/s41467-024-47672-z Download citation * Received: 04 October 2023 * Accepted: 08 April 2024 * Published: 09 May 2024 * DOI: https://doi.org/10.1038/s41467-024-47672-z Share this article Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative Comments By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Download PDF Advertisement Advertisement Explore content * Research articles * Reviews & Analysis * News & Comment * Videos * Collections * Subjects * Follow us on Facebook * Follow us on Twitter * Sign up for alerts * RSS feed About the journal * Aims & Scope * Editors * Journal Information * Open Access Fees and Funding * Calls for Papers * Editorial Values Statement * Journal Metrics * Editors' Highlights * Contact * Editorial policies * Top Articles Publish with us * For authors * For Reviewers * Language editing services * Submit manuscript Search Search articles by subject, keyword or author [ ] Show results from [All journals] Search Advanced search Quick links * Explore articles by subject * Find a job * Guide to authors * Editorial policies Nature Communications (Nat Commun) ISSN 2041-1723 (online) nature.com sitemap About Nature Portfolio * About us * Press releases * Press office * Contact us Discover content * Journals A-Z * Articles by subject * protocols.io * Nature Index Publishing policies * Nature portfolio policies * Open access Author & Researcher services * Reprints & permissions * Research data * Language editing * Scientific editing * Nature Masterclasses * Research Solutions Libraries & institutions * Librarian service & tools * Librarian portal * Open research * Recommend to library Advertising & partnerships * Advertising * Partnerships & Services * Media kits * Branded content Professional development * Nature Careers * Nature Conferences Regional websites * Nature Africa * Nature China * Nature India * Nature Italy * Nature Japan * Nature Middle East * Privacy Policy * Use of cookies * Your privacy choices/Manage cookies * Legal notice * Accessibility statement * Terms & Conditions * Your US state privacy rights Springer Nature (c) 2024 Springer Nature Limited Close Nature Briefing Sign up for the Nature Briefing newsletter -- what matters in science, free to your inbox daily. Email address [ ] Sign up [ ] I agree my information will be processed in accordance with the Nature and Springer Nature Limited Privacy Policy. Close Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing * *