Inertia induces strong orientation fluctuations of non-spherical atmospheric particles

Abstract

The orientation of non-spherical particles in the atmosphere, such as volcanic ash and ice crystals, influences their residence times, and the radiative properties of the atmosphere. Here, we demonstrate experimentally that the orientation of heavy submillimeter spheroids settling in still air exhibits decaying oscillations, whereas it relaxes monotonically in liquids. Theoretical analysis shows that these oscillations are due to particle inertia, caused by the large particle-fluid mass-density ratio. This effect must be accounted for to model solid particles in the atmosphere.

Figures (4)

• Figure 1: Experimental setup. (a) Optical table with top cameras (TX and TY), and bottom cameras (BX and BY) named after the shown coordinate system ($Z$ is the direction of gravity $\hbox{\boldmath$g$}$), the settling chamber (SC), the synchronized pulsed LED unit (LED), controlled with a waveform generator (WG), and the oscilloscope (OSC). (b) Snapshots of a settling prolate spheroid recorded at 2932 frames per second. The tilt-angle -- the angle between the particle symmetry axis and gravity -- is shown in $5.1 {\rm ms}$ intervals in units of degrees. See Supplemental Material (SM) sm for details.
• Figure 2: Comparison between experiments and theory. (a) Time evolution of tilt-angle $\varphi$ from experiment (blue) and its model prediction (red) for spheroids with aspect ratios $\lambda=0.2$ and $5$ from group I (Table \ref{['tab:sd']}). (b) Steady-state settling speed $v^\ast_g$, frequency, and decay rate of the tilt-angle oscillations as functions of the aspect ratio $\lambda$. Markers show averages obtained for all experiments with error bars indicating 95% confidence bounds for groups I (red/$\circ$), II (green/$\Box$), and III (blue/$\diamond$) (Table \ref{['tab:sd']}). Solid lines show large-time asymptotes from a linear-stability analysis of the model described in Appendix A. Shaded regions indicate how much the theoretical predictions change as the measured settling speed varies along the particle trajectory. See Appendices A and B for details. Dashed lines show results of a linear-stability analysis of a harmonic-oscillator approximation, Eq. (\ref{['eq:linear']}) in Appendix D
• Figure 3: (a) Bifurcation diagram. The dashed horizontal line distinguishes decay without oscillations ($\Delta >0$) from decay with oscillations ($\Delta < 0$). Particles from group I in Table \ref{['tab:sd']} are shown as red/$\circ$, group II as green/$\Box$, group III as blue/$\diamond$. Particles from Ref. Cabrera_2022 as $+$, and fibers from Ref. Newsom_Bruce_1994 as a gray region. We approximated cylindrical fibers as slender prolate spheroids and estimated $\Delta$ by setting $C_T = C_F = 1$, since the corresponding ${\rm Re}_{\rm p}$ are very small. (b) Standard deviation of tilt-angle fluctuations for spheroids settling in weakly turbulent air with dissipation rate $\varepsilon$, as a function of ${\rm Re}_{\rm p}$. Shown are simulation results of the model described in Appendix E in symbols, for spheroids with $\lambda = 0.2$ but different volumes, $V_{\rm p}$. Solid lines correspond to the overdamped approximation, neglecting particle inertia (Eq. (S6) in SM sm). The harmonic-oscillator bifurcation ($\Delta = 0$) is shown as a dashed vertical line.
• Figure 4: Empirical coefficients $C_F$ and $C_T$ in Eq. (\ref{['eq:para']}) for the parameters in Fig. \ref{['fig:Fig2']}, obtained as described in Appendix A. (a) Force coefficient $C_F$ as a function of aspect ratio $\lambda$. The groups refer to Table \ref{['tab:sd']}. (b) Torque coefficient $C_T$.