Dynamics of freely rising spheres: the effect of moment of inertia

The goal of this study is to elucidate the effect the particle moment of inertia (MOI) has on the dynamics of spherical particles rising in a quiescent and turbulent fluid. To this end, we performed experiments with varying density ratios $\Gamma$, the ratio of the particle density and fluid density, ranging from $0.37$ up to $0.97$. At each $\Gamma$ the MOI was varied by shifting mass between the shell and the center of the particle to vary $I^*$ (the particle MOI normalised by the MOI of particle with the same weight and a uniform mass distribution). Helical paths are observed for low, and `3D chaotic' trajectories at higher values of $\Gamma$. The present data suggests no influence of $I^*$ on the critical value for this transition $0.42<\Gamma_{\textrm{crit}}<0.52$. For the `3D chaotic' rise mode we identify trends of decreasing particle drag coefficient ($C_d$) and amplitude of oscillation with increasing $I^*$. Due to limited data it remains unclear if a similar dependence exists in the helical regime as well. Path oscillations remain finite for all cases studied and no `rectilinear' mode is encountered, which may be the consequence of allowing for a longer transient distance in the present compared to earlier work. Rotational dynamics did not vary significantly between quiescent and turbulent surroundings, indicating that these are predominantly wake driven.


Introduction
It is widely known that freely rising spheres can exhibit a host of different and complex path oscillations. Numerous studies have been devoted to this topic, which is of interest e.g. as a paradigmatic case for fluid-structure interactions. Canonically, the independent parameters considered are the density ratio Γ = ρ p /ρ f and the particle Reynolds number Re = v z D/ν (or a related quantity such as the Galileo number Ga = |1 − Γ |gD 3 /ν). Here, ρ p and ρ f denote the particle and fluid densities, respectively, · indicates a time or ensemble average and v i is the velocity component of the particle velocity v in direction i (which in the definition of the Galileo number is replaced by the buoyancy velocity V b = |1 − Γ |gD). Further, D is the sphere diameter, ν the kinematic viscosity of the fluid, and g is the acceleration due to gravity. Both parameters, Γ and Re, are related to the vertical momentum balance where F f is the fluid forcing on the body, m f is the particle mass, and e z is a unit vector pointing opposite to the direction of gravity. The Reynolds number dependence enters implicitly in (1.1) via the fluid forcing F f on the sphere. Once Re 200 (Jenny et al. 2003), vortex shedding sets in in the particle wake, which results in an approximately periodic forcing and a complex dynamical coupling between particle motion and the surrounding flow field (Bearman 1984;Parkinson 1989;Williamson & Govardhan 2004;Govardhan & Williamson 2005). The most comprehensive investigation of the Γ -Re parameter space reported to date is by Horowitz & Williamson (2010). Based on their experiments, these authors conclude that a critical density ratio Γ crit exists, which governs the onset of path oscillations. The value of Γ crit was shown to exhibit a Re dependence and path oscillations did not occur for Γ crit > 0.36 for 260 < Re < 1550 and Γ crit > 0.6 at Re > 1550. Horowitz & Williamson (2010) also state that the presence of path oscillations is associated with a high drag regime, for which the values of the drag coefficient C d significantly exceed values reported for a fixed sphere at similar Re. However, there remain fundamental and largely unexplained discrepancies in the literature on the topic. This is most evident in the spread of reported C d values (see figure 3 (b) and the corresponding discussion), but also manifests in differences in the reported rise modes. Whereas Horowitz & Williamson (2010) reported only planar ('zigzaging') trajectories, other studies find helical or spiralling motions (Preukschat 1962;Shafrir 1965;Auguste & Magnaudet 2018;Will et al. 2021;Will & Krug 2021) for comparable parameter values. Also the 'rectilinear mode' described in Horowitz & Williamson (2010) for Γ > Γ crit , in which particles rise straight without path oscillations, is not observed consistently elsewhere; e.g. Preukschat (1962) report a reduction in oscillating amplitude with increasing Γ within this regime, but did not encounter perfectly 'vertical' (i.e. non-oscillating) trajectories.
There certainly are a host of possible explanations for these differences and the origin of some of them may well be linked to the precise experimental conditions. The latter include the precision of the particle fabrication, residual disturbances in the flow, and the size of the tank among potentially many more. However, recent findings (Namkoong et al. 2008;Mathai et al. 2017Mathai et al. , 2018Will & Krug 2021) also suggest a more systematic cause as they point to an additional relevance of the rotational dynamics in setting the overall particle dynamics. Rotations of the sphere are governed by where, ω is the angular velocity of the sphere and T f the torque induced by the fluid. The additional independent parameter introduced by (1.2) is the dimensionless moment of inertia (MOI), I * = I p /I Γ with I Γ = π/60ρ p D 5 the moment of inertia of a particle with a uniform material density of ρ p . Note that the definition of I * is chosen such that its value is entirely determined by the mass distribution within the particle and independent of Γ (i.e. fluid properties). For a homogeneous sphere I * = 1 and I * < 1 (I * > 1) if the mass is accumulated towards (away from) the centre. Note, however, that the dynamically relevant parameter implied by (1.2) is given by the product I * Γ . While there is no explicit coupling between (1.1) and (1.2), the two degrees of freedom can interact via the flow field, e.g. through a Magnus force F m ∼ ω × v. A potential relevance of the MOI as an additional parameter was already mentioned by Ryskin & G a = 6 0 0 0 5 0 0 0 4 0 0 0 3 0 0 0 0.37 < Γ < 0.42 0.52 < Γ < 0.57 0.66 < Γ < 0.70 0.79 < Γ < 0.85 0.90 < Γ < 0.97 3000 < Ga < 3500 4000 < Ga < 4500 4500 < Ga < 5000 5000 < Ga < 5500 Leal (1984). For the case of cylinders, its importance has been established via systematic studies in two-dimensional simulations (Namkoong et al. 2008;Mathai et al. 2017). More recently, Mathai et al. (2018) also uncovered a regime transition induced by a variation in the MoI for spheres rising in a turbulent flow. These authors also reported differences when the particles were rising in still fluid, but these observations remained qualitative and limited to two different values of I * at a single density ratio (Γ ). Moreover, recent experiments by Will & Krug (2021) confirmed the general relevance of rotational dynamics for rising or settling spheres. By introducing a center of mass offset, these authors selectively varied the rotational dynamics while keeping Γ and Ga constant. This led to a resonant behaviour between particle rotation and wake shedding with significant impact on parameters such as oscillation amplitude and C d . Further and remarkably, also horizontal path oscillations ceased almost fully once the offset got large enough to mostly suppress rotational motion.
On this basis, it is the goal of this study to systematically explore the effect variations in I * have on the rise behaviour of light spheres. Little can be gleaned from existing data sets (most of them based on particles with a non-uniform mass distribution) to answer this question, as this parameter is generally not reported. Therefore, we designed and manufactured particles to perform new experiments exploring the parameter space systematically. Details on this can be found in § 2. Afterwards, we present and discuss the results of the experiments in § 3. Additionally, we investigate the effects of background disturbances in the fluid and the effect of the time between experiments (waiting time) (see § 5) and finally conclude in § 6.

Particle design and manufacture
We aim to vary the MOI while keeping Ga and Γ nominally constant. This is achieved by shifting the weight between the outer shell of the particle and a metal ball of varying size at its center as required, see figure 1 (a). The shells are designed using 3D CAD software and 3D printed on a RapidShape 30L printer with a horizontal resolution of 21 µm and a layer thickness of 25 µm. The print is performed in two halves, which are then glued together and sanded to smoothen the surface. In a last step, a pattern is painted on the particle to enable the rotation tracking, resulting in the final particles shown in figure 1 (b). The mass of the paint and glue contributed less than 0.5% of the total particle mass and hence does not induce a significant centre of mass offset nor variation of I * . Measured values of the final particle weight and diameter are used to update the CAD model in order to obtain an more accurate value of the MOI. The final particles were smooth to the touch and we expect that the residual surface roughness will not affect the outcome of the experiments significantly. This is based on the fact that in the present range of Re, the flow over the roughness elements is laminar (Achenbach 1972) and the skin friction will thus not depend on the roughness height (Moody 1944).
The sphericity achieved with this method is high with diameter measurements (taken with a caliper) differing by less than 1% of the diameter at different cross sections. An overview of the parameter space covered in this study is shown in figure 1 (c,d ). Throughout this work the marker colour will designate the different Γ regimes and the different marker types indicate the value of Ga. Line color is used to indicate ranges of I * Γ . For each Γ -regime, I * is varied as much as physically possible, the resulting ranges are shown in figure 1 (d ) in terms of I * Γ .

Experimental setup and methods
All experiments were performed in the approximately 3 m high, water filled, test section of the Twente Water Tunnel facility. The setup is schematically shown in figure 2 (a). The lab is temperature controlled at 20 • C, thus we assume constant fluid properties ρ f = 998 kg m −3 and ν = 1.0035 × 10 −6 m 2 s. The particles were released using a specifically built release mechanism located approximately 1.8 m below the measurement region. The release mechanism, depicted in figure 2 (b), consists of a pipe with a cutout from which the particle can be released. Particles can be inserted into this pipe without draining the tank using a basket. Once pushed to a cutout in the pipe in the center of the tank, this basket can be turned, thereby releasing the trapped air through a grate while keeping the particle inside. After the bubbles have risen to the top of the tank, the water is left to settle for at least 8 minutes before the particle is released by gently tilting the basket further. Doing so did not cause significant rotation of the spheres upon leaving the release mechanism. In §5, we validate the dependence of the behaviour on the waiting time, since this was previously found to be critical to the rise behaviour (Horowitz & Williamson 2010).
Once the particle has entered the measurement section, the particle is recorded by two pairs of perpendicularly placed high-speed cameras (PHOTRON Fastcam AX200 with 1024×1024 pixels at 256 grey levels, fitted with ZEISS Milvus 100mm lenses). The cameras are placed more than 3 meters away from the centre of the tunnel in order to get a near isometric view of the pattern on the particles. Stacking two camera pairs (see figure 2 (a)) allows tracking of the particles over a distance exceeding 1 meter in the vertical direction. Grey background panels are used to contrast with the white and black pattern on the surface of the particles. The frame rate of the cameras was adjusted depending on the rise-velocity in order to keep the inter-frame translation between 2-6 pixels.
Based on the calibration of the camera position and taking into account parallax effects, the 3D particle position in space was reconstructed from the two orthogonal views. The origin of the coordinate system is located at the base of the measurement domain in the centre of the tunnel. The directions are defined as depicted in figure 2 (a), with x and y spanning the horizontal plane and z pointing upward, i.e. opposite to the direction of Schematic detailing the release mechanism and procedure. Left: particle is inserted in the basket and pushed to the centre of the tank. Centre: basket is rotated in the pipe to expose the grate and to allow the air to escape the basket. Right: after waiting for at least 8 minutes, the particle is released by rotating the basket slightly further. (c) An example of a particle trajectory (red) and its center-line (blue). Also shown are projections of the particle path onto the sides and bottom of the domain. Note that horizontal axes are re-scaled with respect to the vertical to highlight the path oscillations. The properties of this particle are: Γ = 0402, I * = 1.140, D = 16.2 mm. (d ) An example of a particle trajectory for a secondary particle with properties: Γ = 0.666, I * = 0.827, D = 19.3 mm.
gravity. The obtained position data is smoothed by convolution with a Gaussian kernel to obtain the trajectories shown as red curves in figures 2 (c, d ). We further obtain the velocity and acceleration of the particle by convolution with the derivatives of a Gaussian kernel (Mordant et al. 2004). The window sizes and standard deviations of the kernels were varied based on the particle size and were determined to minimize the noise, while leaving the underlying signal intact (Mathai et al. 2016). Additionally, we track the orientation of the spheres by matching the particle images to rendered projections of the patterns at different orientations (Will et al. 2021). With the sequence of orientations known, the angular velocity and acceleration can be derived as described in Will & Krug (2021).
From the raw data for the trajectories we then determined the dominant frequency of oscillation by considering the horizontal velocity signals. This result is then used to determine the centre-line of the trajectory indicated by the blue curves in figures 2 (c, d ) (Will et al. 2021). The amplitude of the path oscillations (â) is based on the distance between the trajectory and its centre-line, i.e. the trajectory data are drift-corrected.
Finally, in §4 we perform experiments with a mean downward flow present and active grid generated turbulence in the same facility. In these measurement we only use the top Will & Krug 2021Will et al. 2021 rising settling In all panels, the symbols correspond to ranges in particle Galileo number and the colours indicate ranges in density ratio Γ two cameras and the measurement region is 0.8m downstream of the active grid. In these experiments we release the particle using the same mechanism but with a downward flow in the channel, balancing out the particle rise velocity. The mean flow velocity in the channel is measured using a magnetic flow meter and is kept constant during the experiment. The active turbulence grid at the top of the channel is turned on and when the particle is in the measurement domain the recording is started. The particle can stay in this region for a long time producing recordings in access of a duration of 30 seconds. Typical measurement times easily exceeded 30s, resulting in very good statistics for this configuration. The setup is identical to that used by Mathai et al. (2018) with Re λ ≈ 300 and Ξ = D/η ≈ 100, with η the Kolmogorov length scale of the turbulent flow.

Particle drag coefficient
We start by considering the drag coefficient as a function of the dimensionless moment of inertia I * Γ in figure 3 (a). Most saliently, these results cluster into a high-drag regime (C d ≈ 0.7) and a low-drag regime with C d ≈ 0.45. Albeit not as pronounced as the difference between these regimes, there is further a distinct trend of decreasing C d with increasing I * Γ within the low-drag regime.
In figure 3 (a) we also included relevant data from Will et al. (2021); Will & Krug (2021) for both rising and settling particles. These data points are largely in line with the low-drag mode in the present data set. Only for the lowest values of Γ considered here (0.37 < Γ < 0.42) do we encounter the high-drag regime. Given the limited data points, no conclusions on the dependence of C d on I * can be drawn in this case. Finally, we also note that the transition between the two drag regimes appears independent of I * in the present data.
To establish how the observed trends -and in particular the dependence of C d on I * Γ -relate to literature data, we compare our data to published values of C d (Preukschat 1962;MacCready & Jex 1964;Shafrir 1965;Kuwabara et al. 1983;Jenny et al. 2004;Karamanev et al. 1996;Stringham et al. 1969;Veldhuis & Biesheuvel 2007;Veldhuis et al. 2009;Allen 1900;Liebster 1927;Lunnon 1928;Boillat & Graf 1981;Will et al. 2021;Will & Krug 2021) in figure 3 (b). The existence of two different drag states as a function of Γ has been observed and documented before by Horowitz & Williamson (2010). They attributed the regimes to a transition from a "vibrating mode", with large path oscillations at low Γ , to a rectilinear mode with almost no path oscillations. The respective C d values for these two regimes largely match our results (see also inset in figure 3 (a)). However, the threshold density ratio in the present data (0.42 Γ 0.52) is significantly lower compared to the critical value of Γ = 0.61 determined in Horowitz & Williamson (2010) (for unknown I * ). Regarding the dependence of C d on I * , it is noteworthy that the C d values of the rectilinear mode -and at the same time also those for a stationary sphere-are best matched at high I * Γ . The increase in C d with decreasing I * Γ then leads to a deviation from these reference data. Some part, but certainly not all of the spread in the drag data reported in the literature, might therefore indeed be attributed to differences in the moment of inertia. However, the effect appears less strong compared to the Γ -dependence, which remains the dominant parameter in governing the particle drag.
Finally, we also document the particle Reynolds number defined as Re = v z n D/ν. These results are shown in figure 3 (c) as a function of Ga. The Re-range for all our experiments are in the Newtonian drag regime, Re 1000 (Clift & Gauvin 1971), for which the drag is pressure dominated (C d ∼ D 2 ) and mostly independent of Re. The Reynolds and Galileo numbers are linked to the particle drag coefficient as shown in (3.1) and indicated by diagonal dashed lines in the figure with bi-logarithmic scales. The variation in Re is small between most particles and there appears to be no systematic relation between C d and Re. This indicates that the variation in C d observed in figure 3 (a) in the low drag mode is indeed an I * effect and not related to Re. In the following, we will investigate in more detail how varying Γ and the MOI changes properties of the particle trajectories.

Particle trajectories
In figure 4 a number of representative trajectories for the five ranges of Γ and for three values of I * are shown. To provide a sense of the variability, the main panel shows the most commonly observed trajectory, while the inset shows a second trajectory for the same particle that is (visually) most different from the typical one. Surprisingly, the horizontal trajectories corresponding to the high drag-regime (0.37 < Γ < 0.42) are almost circular, indicating that the path is helical for these cases. This result is unlike the planar 'zigzag' observed by Horowitz & Williamson (2008, 2010 for spheres with a low density ratio at similar Ga. However, helical trajectories are not uncommon and have been observed for spheres at low density ratios (Preukschat 1962;Karamanev & Nikolov 1992;Karamanev et al. 1996;Veldhuis et al. 2009) as well as for bubbles (Ellingsen & Risso 2001;Mougin & Magnaudet 2001, 2006. In the work by Veldhuis et al. (2009), these helical trajectories were connected to a different drag scaling, which could be attributed to lift-induced drag resulting from the shedding of additional vorticity in the wake. This non standard drag behaviour for spiralling trajectories was also noted in the numerical . Representative particle trajectories in the horizontal (x − y) plane. The trajectories are colour coded using the normalized vertical velocity vz/V b , symbols mark the start (circle) and end (cross) of the recorded path. Γ increases from left to right (columns (a-e) and vertically in each column the value of I * increases from top to bottom. For each case two trajectories are shown: the main window shows the most frequently occurring trajectory and the inset shows the most aberrant run for that same particle. work by Auguste & Magnaudet (2018). They found a similar spiralling regime at low Γ , but their data set is limited to Ga 700. Karamanev & Nikolov (1992) state that all particles with Re > 130 and Γ 0.3 rose in a spiralling trajectory. These spirals featured a constant angle pitch angle between v and e z of about ±29 • . For the spiralling regime here, we find this angle to be around 18.5 • independent of I * . In this context it is also important to note that trajectories became more circular for spheres with centre of mass offset in resonance with their natural frequency (Will & Krug 2021). However, this can be ruled out as a factor here since the required offset of 0.03D certainly exceeds our fabrication tolerance. Furthermore, it is extremely unlikely to randomly hit resonance for all particles across a range of I * values. We are therefore convinced that the present trajectories reflect the genuine particle behaviour at the present values of Ga, Γ and I * .
The dynamics in the low drag-regime for Γ > 0.52 are distinctly different from the spiralling motion at low density ratios. While there remains periodicity in the sequence of direction changes, the turning angles appear random. The resulting behaviour is characteristic of the '3D chaotic' regime (Auguste & Magnaudet 2018), which applies to all particles with Γ > 0.52 here. These results are at odds with the findings of Horowitz & Williamson (2010), who found a vertical rise regime for all particles larger than Γ ≈ 0.61 at comparable Ga.
Finally, there is no clear trend visible in the shape of the horizontal trajectories for varying MOIs. We will proceed to investigate the horizontal motion more quantitatively in order to elucidate such effects.

Fluctuations of the horizontal velocity
To uncover the time-varying dynamics of the rising spheres, it is useful to consider the probability density functions (PDFs) of the horizontal velocity. The distributions of the normalised horizontal velocity ( v 2 x + v 2 y /V b ) are shown for various values of I * Γ at different Γ in figure 5 (a-e). As a most obvious trend, we note a change in skewness from negative for Γ 0.42 (figure 5 (a)) to positive skew at Γ 0.52 ( figure 5 (b-e)). This is indicative of the transition from spiralling motion, for which the horizontal velocity is generally high, to the '3D chaotic' state, for which strong horizontal translation occurs more intermittently.
The effect of varying I * Γ (indicated by the line colour) is negligible at 0.37 < Γ < 0.42, as evidenced by figure 5 (a). Given the limited (by physical constraints) range of I * at this density ratio, it remains unclear to what extent this indicates a lesser importance of the rotational dynamics in the spiralling regime (see also § 3.5 on this).
For Γ 0.52 (i.e. in the '3D chaotic' state), however, a clear dependence of the PDFs of v 2 x + v 2 y /V b on I * Γ emerges. We observe that for low I * Γ , the distribution is in general broader and extends further towards high velocities, whereas at high I * Γ the distributions is narrow and the peak shift towards lower velocities. This trend is most evident in figure 5 (c, e), showing that at very low I * the distributions become rather flat, reminiscent of a fluttering behaviour with more pronounced horizontal motion.
It should be noted that the trends discussed here for v 2 x + v 2 y /V b also manifest in the statistics of the vertical velocity v z (not shown), albeit with inverse effects regarding the skewness and the MOI dependence. While there is a distinct negative correlation between instantaneous vertical velocity and instantaneous horizontal velocity, which is also obvious from the colour coding in figure 4 (b-e)), the magnitude of the velocity, ||v|| is not constant but fluctuates quasi-periodically for all cases.

Oscillation frequency and amplitude
In all experiments performed here, we observed significant oscillations in the trajectory of the rising spheres. In this section, we will characterize these in terms of the frequency and amplitude of the horizontal oscillations. In figure 6 (a), we start by plotting the frequency f as a function of the dimensionless moment of inertia I * Γ . It is evident from this figure that there is no significant dependence of f on I * Γ as the data points with the same colour (i.e. constant Γ ) are at a near constant frequency. Note that the error bars on the frequency, indicating the spread in the data, are very small for all density ratios besides 0.52 < Γ < 0.57. This indicates that even though some of the motion appears quite random, there exists a strong dominant frequency associated with the Will & Krug 2021Will et al. 2021 rising settling 3000 < Ga < 3500 4000 < Ga < 4500 4500 < Ga < 5000 5000 < Ga < 5500 0.37 < Γ < 0.42 0.52 < Γ < 0.57 0.66 < Γ < 0.70 0.79 < Γ < 0.85 0.90 < Γ < 0.97 (a) (b) (c) vortex shedding. The outlier at 0.52 < Γ < 0.57 is most likely related to this case falling within the transitional regime between helical paths and 3D chaotic patterns. In figure 6 (b), we show the frequency in dimensionless form in terms of the Strouhal number, defined as: This normalization separates the data into two regimes, akin to those encountered for C d and for the trajectories. For Γ > 0.66 we find that the Strouhal number takes a constant value of approximately Str = 0.04-0.05, while for 0.37 < Γ < 0.42 we find Str ≈ 0.09. Only the data for 0.52 < Γ < 0.57 does not completely fall in line with this decomposition and lies at a slightly higher value of Str ≈ 0.06, consistent with the transitional behavior of this case mentioned above. Furthermore, also the Strouhal number appears to be rather insensitive to changes in I * . The value of Str ≈ 0.09 at low Γ matches the results by Horowitz & Williamson (2010) (figure 31) closely. They find for identical Ga and Γ also Str ≈ 0.09, however, the sphere is zigzagging in their case instead of spiralling as observed here. The value of Str ≈ 0.04 − 0.05 for Γ > 0.66 is in line with the results by Preukschat (1962) (see fig. 19 and page 20), who for 0.582 Γ 0.875 also find Str as low as 0.05.
Finally, we report the normalized amplitude of the path oscillationsâ/D as a function of the I * Γ (see figure 6 (c)). Despite the scatter in these data, there appears to be a consistent trend ofâ/D decreasing from approximately 0.8 to 0.3 with increasing I * Γ over the full range accessible here. Surprisingly, the values of the amplitude are congruent for both the helical regime as well as for the more chaotic trajectories. Our results imply that decreasing I * Γ by either changing the internal structure of the particle or the density ratio, results in larger amplitude path oscillations. In fact, it appears that a similar trend with consistent amplitudes exists in the data of Horowitz & Williamson (2010) for varying Γ in their 'zigzag' regime. However, in their experiments I * is not monitored nor controlled, which renders a quantitative comparison of both datasets impossible.

Rotational dynamics and translational coupling
The mechanism by which the particle moment of inertia I * affects the particle kinematics and dynamics is solely through the rotational equation of motion 1.2, effectively scaling the particle rotation to the fluid torques. Rotational dynamics in turn affect the B N 1 0 T 3000 < Ga < 3500 4000 < Ga < 4500 4500 < Ga < 5000 5000 < Ga < 5500 0.37 < Γ < 0.42 0.52 < Γ < 0.57 0.66 < Γ < 0.70 0.79 < Γ < 0.85 0.90 < Γ < 0.97 T is parallel to the instantaneous direction of motion, N is parallel to the curvature of the path and B is defined parallel to T × N . Definition of the direction of rotation ω using azimuth (φ) and elevation (θ) with respect to the TNB coordinate system. (c) Normalized histograms of the alignment of ω in terms of φ and θ in the TNB coordinate system as a function of particle density ratio (Γ ). These figures contain the data for all I * -values within this Γ range. In these figures θ = 90 is aligned with T , φ = −90, θ = 0 is aligned with N , and φ = 0, θ = 0 is aligned with B.
flow field around the body, thereby inducing a coupling with the lateral motion, e.g. via Magnus lift type forcing. Particle rotation can also affect vortex detachment and consequently the flow structure in the wake of the particle, an effect that is believed to be at the heart of the regime transition observed by Mathai et al. (2018). Therefore, the most direct parameter in investigating the impact of the MOI is the rotation rate of the body. This quantity is explored in figure 7 (a), where we plot the mean dimensionless rotation rate ω * = ||ω|| D/V b versus I * Γ . We find that the dimensionless rotation rate (ω * ), similar to the drag and amplitude of the oscillation of the trajectory, shows a slight dependence on I * Γ . The particles with lower rotational inertia indeed rotate more vigorously compared to their higher I * Γ counterparts. We further note that the normalization of ||ω|| with D and V b collapses the results across all density ratios convincingly. It is important to note that the rotation rate is not affected significantly by the change from spiralling to the 3D chaotic regime.
In order to explore the importance of rotational dynamics further, we examine the alignment of the rotation vector ω with respect to the particle acceleration along the curvature of the path. Doing so allows us to establish the relevance of Magnus lift forcing on the particle dynamics. To this end, we consider ω in the Frenet-Serret (TNB) coordinate system (Zimmermann et al. 2011). As is shown in figure 7 (b), T points in the direction of the instantaneous velocity v of the sphere, N is aligned with the curvature of the path (the acceleration of the sphere that is non-parallel to the direction of motion; a ⊥T ), and B = T × N . Thus, when ω B the induced Magnus force F m a ⊥T , making this coordinate system very useful to study the effect of rotation on the horizontal motion. The alignment of ω within the TNB coordinate system is given in terms of two angles: the azimuth φ and the elevation θ as indicated in figure 7 (b). The directions of N and B are also indicated in figure 7 (d ), where we show normalized histograms of the alignment 0.37 < Γ < 0.42 of ω in the φ − θ plane for different values of Γ . Most striking about these results is the enormous difference in rotational alignments between the low-Γ and high-Γ regimes. For the low-Γ regime (0.37 < Γ < 0.42), we find that, depending on the direction of the spiralling motion (clockwise or counter-clockwise), the alignment of ω is either θ ≈ −40 • or 40 • , respectively. As illustrated in figure 7 (b) for the counter-clockwise spiralling case (but also true for clockwise rotation), this alignment implies that the Magnus lift force (F m ∼ ω × v) predominantly acts downward. This mechanism therefore leads to lift induced drag, in a manner similar to that suggested by Mougin & Magnaudet (2006) and Veldhuis et al. (2009). Furthermore, with |φ| = 90 • the Magnus force acts perpendicular to a ⊥v , and is hence not responsible for the lateral acceleration for the lowest Γ . Thus, the spiralling motion encountered for these particles must have a different origin, e.g. a rotation of the vortex shedding position relative to the direction of motion (Karamanev et al. 1996). On the other hand, for Γ > 0.52 the distribution of orientations of ω is centred around B. The histograms of all cases with Γ > 0.52 are remarkably similar underlining that these belong to the same dynamical regime. The data clusters around θ ≈ 0, i.e. ω pointing normal to v. This is consistent with vortex shedding being a main driver of particle rotation, as this induces a torque perpendicular to T . Further, θ ≈ 0 on average implies that there is no net contribution of the Magnus force in the vertical direction.
So far, we have only considered how the alignment statistics of ω depend on Γ . Potential variations with I * Γ are masked by plotting data with different MOI together in figure 7 (d ). To elaborate on the influence of the MOI, we present additional orientation statistics of ω for 0.90 < Γ < 0.97 at five different I * Γ values separately in figure 8 (a). At the lowest MOI (I * Γ ≈ 0.48), the alignment of ω and B is very strong. This behaviour is associated with a strong coupling of particle rotation and lateral acceleration resulting in a "fluttering" type of behaviour (Mathai et al. 2018). The alignment with B progressively weakens for higher values of I * and the distribution is approximately flat in terms of φ at I * Γ ≈ 1.10. We note that the trend observed in figure 8 (a) is representative also for the lower Γ cases in the 3D chaotic regime. However, the alignment with B at the respectively lowest I * values is not equally as pronounced as for 0.90 < Γ < 0.97. It hence appears that the 'flutter' type motion occurs at different values of I * Γ at different Γ .
A related analysis, that allows to more readily compare different combinations of I * and Γ , is to evaluate the autocorrelation C ω of an arbitrary horizontal component of ω as was done previously in Mathai et al. (2018). Such results are shown for 3 different ranges of the density ratio in figure 8 (b-d ), where the correlation coefficient is plotted against the dimensionless time lags t/τ vs , with the vortex-shedding timescale τ vs ∼ D/ v z n . The influence of the MOI seen in figure 8 (a) is reflected also in the corresponding results in figure 8 (d ). In this case, it manifests itself by a gradual transition from a periodic behaviour of C ω (low I * Γ , blue line) to a slow decorrelation without significant oscillation at large I * Γ (red lines). Note, however, that the drop to 0 at the largest t/τ vs shown is an artefact of the finite observation time in all cases. A similar influence of the MOI as seen at Γ ≈ 0.9 is also found at the lowest Γ in the 3D chaotic regime in figure 8 (c). Here, the correlation timescale outgrows the vortex shedding time scale τ vs at larger I * Γ , while at the same time the correlation also drops in magnitude. For the lowest values of Γ ( figure 8 (b)), for which spiralling trajectories are observed, the periodicity in C ω is most pronounced of all the cases shown. Typical periods are about 0.5 t/τ vs with no obvious dependence on I * Γ . Non-periodic rotations are observed most prominently at high Γ and high MOI, which is consistent with the mostly random orientation of the ω vector in the TNB coordinate system (see figure 8 (a) for these cases. Such behaviour was previously only reported in turbulent flow (Mathai et al. 2018) and we will study this aspect in more detail in the next section.

On the effect of turbulence
Previously, Mathai et al. (2018) investigated the effect of MOI on spheres rising in a downward turbulent flow. They showed that depending on their MOI the trajectories of spheres with Ga ≈ 6000 and Γ ≈ 0.88 differed significantly: For I * Γ = 0.6 a 'zig-zag' (or 'flutter') motion was observed without significant horizontal drift, whereas for I * Γ = 1, the behaviour changed to a tumbling motion featuring a strong mean drift.
To test if such a transition can also be observed for the present set of particles, we conducted measurements in turbulent flow. The turbulence properties Re λ ≈ 300 and Ξ = D/η ≈ 100, with η the Kolmogorov length scale, were chosen to match those of Mathai et al. (2018). Further, we restricted measurements to two particles, both with Γ ≈ 0.9 and I * Γ = 0.52 and I * Γ = 1.18, respectively, again matching the literature conditions closely.
Clearly and consistently with Mathai et al. (2018), the two particles behave differently in the turbulent flow as evidenced by the results for C ω in figure 9 (a). However, discounting the differences caused by limited observation times in the still fluid, there appears to be little difference for the same particle in either quiescent or turbulent surroundings. This holds for C ω , but also for the fact that the tumbling behaviour reported for the turbulence measurements in Mathai et al. (2018) was not observed here. Consequently, we also see no difference in the mean drift rate of the two particles (29.07 mm/s and 26.67 mm/s for low and high MOI, respectively). It therefore appears that whether or  not a sphere tumbles in turbulence might depend quite sensitively on the flow or particle parameters (e.g. values of D/η differ by about 10% between the two studies). Other MOI related effects, such as the decorrelation of the rotation rate and arguably the overall rise pattern, appear largely independent of the flow state.
Naturally, small differences between the turbulent and the quiescent case arise for specific parameters and some of these are documented in figures 9 (b-e). The drag coefficient C d (figure 9 (b)) is approximately 20% higher in turbulence for both particles. Similarly, also the typical oscillation amplitude is higher in turbulent surroundings (figure 9 (c)). Interestingly, the mean rotation rate (see figure 9 (d )) is barely affected by the ambient flow, which may indicate that these dynamics remain wake driven even in turbulence. For all these quantities the trend that higher values are observed at lower I * Γ is also preserved in turbulence. The only exception to this might be the Strouhal number in figure 9 (e), for which no MOI dependence is visible for the turbulence results. However, with a standard deviation of 0.1 the spread in Str is much larger for these data compared to those from the quiescent experiments. The difference at a given I * Γ is therefore only of the order of the typical fluctuations of Str in turbulence.

The effect of disturbances in fluid
Apart from the effect of artificially introduced (turbulent) velocity fluctuations, another point of interest is the effect of residual background disturbances in the nominally quiescent fluid. This point was stressed in the work by Horowitz & Williamson (2010). They found that the emergence of the rectilinear regime was very sensitive to the flow conditions. Any remaining disturbances in the tank caused by the insertion of the release mechanism, or caused by previous experiments, was sufficient to result in path oscillations. Since this rectilinear regime was not observed here in the parameter range specified by Horowitz & Williamson (2010), it is therefore imperative to check if this might be related to flow disturbances.
For this purpose, we performed a series of experiments varying the waiting time between the evacuation of the bubbles from the release mechanism (figure 2 (b)) and the particle release. All measurements were performed with a single particle with properties (Γ = 0.907, Ga = 4800 and I * = 0.822) well within the rectilinear regime of Horowitz & Williamson (2010). The waiting times considered were 2,4,8,16,32,64,128 and > 360 minutes with 5 repeats performed for each of these cases. None of these experiments displayed a rectilinear rise behaviour, and the influence of the settling time on the outcomes of the experiments was generally very limited. This is exemplified in figure 10 (a, b), presenting two representative trajectories, taken after waiting for 8 minutes and for more than 6 hours, respectively. Visually, there is no discernible difference between the two. For a more quantitative analysis, we present in figure 10 (c, d, e) results for the oscillation frequency f and amplitudeâ, as well as for the particle drag coefficient as a function of the waiting time. Both f andâ appear insensitive to the waiting time with differences within the range of the error bars, which primarily reflect the period-to-period variations that occur in these experiments. There is a slight initial decrease for C d with longer waiting times. However, around a waiting time of 10 minutes the curve levels off and C d stays constant for even longer waiting times. In summary, we therefore conclude that with waiting times > 8 minutes our results presented in § 3 are not affected by residual disturbances in any significant manner. An interesting observation in our experiments was that close to particle release, there was very little rotation of the particle and it appeared to rise more vertically compared to higher up in the measurement section of the water tank. The puzzling differences between the present results and those of Horowitz & Williamson (2010) might therefore be related to a relatively long transient period, in which the spheres build up rotational kinetic energy. This is in line with a transient over 60-80D in the present study, whereas the particles only travelled for ≈ 10D in Horowitz & Williamson (2010) before entering the measurement section.

Conclusion
The systematic study performed here showed mixed results regarding the relevance of MOI variations for free rising spheres. On the one side, we identified a dependence of the drag coefficient on the MOI in the 3D chaotic regime. Namely, we found that C d increases with decreasing I * Γ . This is accompanied by a trend of stronger rotations at lower I * Γ and therefore in line with findings of Will & Krug (2021), who suggested a correlation between rotation rate and drag. Furthermore, our results indicate an increasing amplitude of horizontal path oscillations at lower I * Γ , which confirms similar qualitative observations in Mathai et al. (2018). Finally, at lower values of the MOI there is a stronger alignment between the rotation vector and the acceleration normal to the path. This indicates a more dominant rotational-translational coupling in this case and leads to fluttering-type trajectories.
On the other hand, the effect of MOI variations was limited to altering parameters within the 3D chaotic rise mode. It did not lead to regime changes in the overall behaviour, not even in case of turbulent flow, for which a transition from flutter to tumble has been reported previously (Mathai et al. 2018). The variation, e.g in C d , induced by variations of I * Γ is hence limited and is not adequate to explain the scatter in the literature data. It appears that the rotational dynamics induced by centre of mass offsets, which can lead to drastic changes in the drag (Will & Krug 2021), play a more important role in this context.
Apart from the effect of the MOI, our results are also useful in revisiting the broader studied Γ -dependence, which together with Re is the dominant parameter in governing the rise mode. Our analysis revealed the existence of a spiralling mode for Γ 0.42 for which the drag is significantly elevated (C d ≈ 0.69) similar to that found by Auguste & Magnaudet (2018) at lower Ga. Since neither the lateral amplitude nor the rotation rate excel similarly in this regime, it appears likely that the higher drag originates from a change in the wake structure. 3D chaotic motion was found for Γ > 0.52 with C d ≈ 0.47 (and the I * Γ dependence outlined above). This drag value is close to that of a fixed sphere at comparable Re and matches also the results of Preukschat (1962) in the same range of density ratios (Γ ) and the cases in Horowitz & Williamson (2010) for Γ 0.61. Horowitz & Williamson (2010) also identified a density ratio triggered regime transition consistent with our observations, albeit with a different critical density ratio Γ crit ≈ 0.61 in the same Ga regime. Other relevant differences between the present results and those of Horowitz & Williamson (2010) relate to the existence of a rectilinear regime at larger Γ , which -potentially due to a longer initial transient (see § 5)-is not observed in the present study. Based on varying the waiting time between experiments, we could rule out residual disturbances as a cause for this discrepancy. In this regard our results appear consistent with Preukschat (1962), who report a gradually decreasing but finite amplitude as Γ → 1, in good agreement with the I * Γ dependence established here. Finally, even though the spiralling trajectories found at Γ 0.42 are in line with Karamanev & Nikolov (1992) and Auguste & Magnaudet (2018) at lower Ga, they are at odds with the 2D zigzag reported in Horowitz & Williamson (2010)

Declaration of Interests
The authors report no conflict of interest.