1. Introduction
The dynamics of objects in fluid flow can be incredibly intricate. Even restricting consideration to a simple shear flow in the Stokes regime, which arises as a key component in the local linearisation of general laminar flows, the corresponding equations of motion can be difficult to analyse explicitly. However, there are certain classes of object for which analytic solutions are known, beginning with the celebrated work of Jeffery (Reference Jeffery1922). Jeffery's seminal study considered the orientation of ellipsoids in a shear flow, deriving an exact dynamical system of three governing equations for a general ellipsoid. Moreover, for a spheroid (also referred to as an ellipsoid of revolution), the governing equations simplify to what are now known as Jeffery's equations, for which exact solutions can be found explicitly. These solutions were shown to be in agreement with experiments by Taylor (Reference Taylor1923) one year later.
Given the impact of the contribution by Jeffery (Reference Jeffery1922), the solutions to the derived system of equations are known as Jeffery's orbits. Furthermore, these governing equations apply far beyond the stated geometric restrictions of Jeffery's original study. In particular, the work of Bretherton (Reference Bretherton1962) and Brenner (Reference Brenner1964) showed that Jeffery's equations also hold for the motion of general axisymmetric objects, with the equations of angular motion reducing simply to those that describe the evolution of a spheroid in flow. That is, the rotational dynamics of general axisymmetric objects behave as effective spheroids in shear flow, with a corresponding effective aspect ratio that depends on the particular object.
In all the work discussed above, the bodies are passive. However, many bodies of interest are active, such as living microswimmers or externally driven particles. The consideration of active bodies contributes additional propulsion and rotation to the dynamics, resulting in behaviours that are significantly more intricate than their passive counterparts (Wittkowski & Löwen Reference Wittkowski and Löwen2012; Junot et al. Reference Junot, Figueroa-Morales, Darnige, Lindner, Soto, Auradou and Clément2019; Fung, Bearon & Hwang Reference Fung, Bearon and Hwang2022). Active matter has been the subject of much recent research, including studies of bacterial locomotion across scales (Hyon et al. Reference Hyon, Marcos, Powers, Stocker and Fu2012; Constantino et al. Reference Constantino, Jabbarzadeh, Fu and Bansil2016; Lauga Reference Lauga2016; Aranson Reference Aranson2022), the swimming behaviours of spermatozoon (Gaffney, Ishimoto & Walker Reference Gaffney, Ishimoto and Walker2021; Gong et al. Reference Gong, Rode, Gompper, Kaupp, Elgeti, Friedrich and Alvarez2021) and in active suspensions more generally (Saintillan & Shelley Reference Saintillan and Shelley2013, Reference Saintillan and Shelley2015, Reference Saintillan and Shelley2008). Driven by advances in micromanipulation and microfabrication, there has also been the theoretical and experimental study of externally driven artificial active matter, such as the rigid, magnetic helices explored by Ghosh & Fischer (Reference Ghosh and Fischer2009) and Zhang et al. (Reference Zhang, Abbott, Dong, Kratochvil, Bell and Nelson2009) (see also Zhou et al. (Reference Zhou, Mayorga-Martinez, Pané, Zhang and Pumera2021) for a recent review). In many active swimmer systems, including both biological and artificial active matter, motion often proceeds over separated time scales. Rapid movements (often oscillatory) dominate small-scale dynamics over shorter time scales, and these give way to larger scale, emergent trajectories over longer time scales, as may be readily inferred by inspecting the length and time scales of observed swimmer trajectories. Some examples include the trajectories of Crithidia deanei in figure 1(b) of Gadelha, Wickstead & Gull (Reference Gadelha, Wickstead and Gull2007); mouse sperm in figure 13(a) of Woolley (Reference Woolley2003); a magnetically driven bacterium-like artificial swimmer in figure 4(d) of Huang et al. (Reference Huang, Sakar, Petruska, Pane and Nelson2016); and bacteria in figure S3 of Turner et al. (Reference Turner, Ping, Neubauer and Berg2016). In particular, while bacteria flagellum rotation rates can vary extensively with viscosity, number of flagella and other factors, some bacteria exhibit more than $10^3$ flagellum rotations per second, as reviewed by Berg (Reference Berg2000). One example is smooth swimming monotrichous Vibrio alginolyticus mutants lacking lateral flagella (Magariyama et al. Reference Magariyama, Sugiyama, Muramoto, Kawagishi, Imae and Kudot1995), where the longer length- and time-scale dynamics of smooth swimming emerges from rapid rotation. Furthermore, theoretical studies also highlight the emergence of larger-scale trajectories from rapid small-scale movements, for example, the direct numerical simulations of Hyon et al. (Reference Hyon, Marcos, Powers, Stocker and Fu2012) and Park, Kim & Lim (Reference Park, Kim and Lim2019) for bacteria.
Recently, examples of such multiscale dynamics have been explored using the asymptotic method of multiple scales (Hinch Reference Hinch1991; Bender & Orszag Reference Bender and Orszag1999), yielding systematically simplified dynamical systems that govern the dominant behaviours of objects and swimmers in flow (Gaffney et al. Reference Gaffney, Dalwadi, Moreau, Ishimoto and Walker2022; Ma, Pujara & Thiffeault Reference Ma, Pujara and Thiffeault2022; Walker et al. Reference Walker, Ishimoto, Moreau, Gaffney and Dalwadi2022a). Motivated by the recent successes of the method of multiple scales in different swimmer contexts, which have been predominantly limited to motion in two dimensions (with the notable exception of the three-dimensional (3-D) model in the appendix of Ma et al. Reference Ma, Pujara and Thiffeault2022), the primary aim of this study is to exploit separated time scales in a 3-D model of an active particle to derive effective governing equations. Here, drawing inspiration from a general class of rapidly spinning bodies that include bacteria-like swimmers, we investigate the dynamics of objects whose self-generated spinning motion is fast compared with the flow and self-generated translation time scales in the swimming problem. In particular, we focus on objects whose geometry can be reasonably approximated as constant in time, noting that this assumption is often a feature of externally driven rigid objects and simple models of shape-changing swimmers (Lauga & Powers Reference Lauga and Powers2009; Lauga Reference Lauga2020). This approach is also consistent with the popular squirmer model, in which a self-generated slip velocity drives the motion of a rigid object (Magar & Pedley Reference Magar and Pedley2005; Ishimoto & Gaffney Reference Ishimoto and Gaffney2013; Ito, Omori & Ishikawa Reference Ito, Omori and Ishikawa2019).
Seeking generality, we abstract away from any particular application or interpretation to develop an appropriate general framework in Part 1, before using this to study spheroidal objects. In Reference Dalwadi, Moreau, Gaffney, Ishimoto and WalkerPart 2 (Dalwadi et al. Reference Dalwadi, Moreau, Gaffney, Ishimoto and Walker2023) we apply the framework we develop in Part 1 to study more general-shaped objects (chiral particles with helicoidal symmetry). Hence, in Part 1 we are essentially investigating active counterparts of the passive spheroids in shear flow analysed by Jeffery (Reference Jeffery1922) to understand the effect of self-driven translation and fast rotation on their dynamics. In order to retain tractability, we study particles with minimally complex active characteristics. Specifically, we assume that our objects generate a prescribed constant linear velocity and fast angular velocity, each fixed in a reference frame that evolves with the object orientation. Even with this minimal additional complexity, the governing equations we derive are a strongly nonlinear coupled system that requires a non-standard and technically intricate multiple scales analysis. In particular, we must solve a nonlinear coupled 3-D system at leading order and, using the method of multiple scales for systems (see, for example, pp. 127–128 of Dalwadi (Reference Dalwadi2014) or p. 22 of Dalwadi et al. Reference Dalwadi, Chapman, Oliver and Waters2018), a non-self-adjoint coupled 3-D system at next order.
Nevertheless, we are able to make significant analytic progress to derive emergent equations of motion that systematically account for the rapid spinning. We show that there is a natural and fundamental correspondence between these emergent equations and the classic equations of Jeffery (Reference Jeffery1922). Specifically, in terms of quantities that we derive analytically, we show that the emergent equations of motion are equivalent to Jeffery's equations of motion in appropriately transformed variables, with effective coefficients that describe the effective hydrodynamic spheroidal shape in terms of the original shape parameters and the rapid rotation.
To summarise, in this first part of our two-part study, we develop the multiscale theoretical framework to consider a rapidly spinning self-propelled particle in Stokes flow, and apply it to the specific problem of a spheroidal particle in 3-D shear flow. We start by formulating the equations of motion for such a swimmer in § 2. To introduce the multiple scales notation for the full system, we first analyse a more intuitive limit of the full system in § 3, limiting the axis of rapid spinning to the symmetry axis of the spheroid. This also allows us to illustrate the fundamental steps involved in a standard multiple scales analysis before using an appropriate modification for systems to analyse the full problem in § 4, which includes general off-axis rapid spinning. We then summarise and discuss the key physical implications and conclusions of our technical analysis in § 5 in a non-technical manner, for the benefit of those that wish to skip the details of the analysis. Finally, we discuss our results and their wider implications in § 6.
2. Equations of motion
We consider the motion of a rigid, self-propelled spheroid in Stokes flow with imposed far-field shear flow, as illustrated in figure 1. Henceforth, we non-dimensionalise time with the inverse shear rate and space with the equatorial radius of the spheroid, so all subsequent variables and parameters are dimensionless with respect to these scalings. The spheroid is therefore defined by a dimensionless equatorial radius of $1$ and distance from centre to pole of $r$ along the spheroid symmetry axis. In a quiescent fluid the self-generated swimming velocity of the spheroid is $\boldsymbol {V}$, and its self-generated angular velocity is $\boldsymbol {\varOmega }$. We assume that both $\boldsymbol {V}$ and $\boldsymbol {\varOmega }$ are constant in a swimmer-fixed reference frame, and that their generating forces and torques are unchanged by the imposition of the external flow. The swimmer-fixed frame will depend on the orientation of the object relative to the laboratory frame, and this orientation depends on $\boldsymbol {\varOmega }$. We define these two frames below.
We define the spheroidal axis of symmetry via a swimmer-fixed axis $\hat {\boldsymbol {e}}_{1}$. Then, without loss of generality, we choose $\hat {\boldsymbol {e}}_{2}$ such that $\boldsymbol {\varOmega }$ is in a plane spanned by $\hat {\boldsymbol {e}}_{1}$ and $\hat {\boldsymbol {e}}_{2}$, where $\boldsymbol {\varOmega }$ makes an angle of $\alpha \in [-{\rm \pi} /2,{\rm \pi} /2]$ with $\hat {\boldsymbol {e}}_{1}$. We therefore write $\boldsymbol {\varOmega } = \varOmega _{\parallel }\hat {\boldsymbol {e}}_{1} + \varOmega _{\perp }\hat {\boldsymbol {e}}_{2}$, with $\varOmega _{\parallel }$ and $\varOmega _{\perp }$ being the components of angular velocity that are parallel and perpendicular to the symmetry axis, respectively, and hence, $\tan \alpha = \varOmega _{\perp }/\varOmega _{\parallel }$. Then $\hat {\boldsymbol {e}}_{3} = \hat {\boldsymbol {e}}_{1} \times \hat {\boldsymbol {e}}_{2}$. (If $\boldsymbol {\varOmega } \parallel \hat {\boldsymbol {e}}_{1}$ then $\hat {\boldsymbol {e}}_{2}$ and $\hat {\boldsymbol {e}}_{3}$ can be chosen arbitrarily to complete the orthonormal triad.) In this swimmer basis, we write $\boldsymbol {V} = V_1\hat {\boldsymbol {e}}_{1} + V_2\hat {\boldsymbol {e}}_{2} + V_3\hat {\boldsymbol {e}}_{3}$, while the centre of the spheroid is given by $\boldsymbol {X}(t) = X(t) \boldsymbol {e}_{1} + Y(t) \boldsymbol {e}_{2} + Z(t) \boldsymbol {e}_{3}$ with respect to the orthonormal basis $\{\boldsymbol {e}_{1},\boldsymbol {e}_{2},\boldsymbol {e}_{3}\}$ of the laboratory frame. We illustrate these vectors in figure 1.
Finally, we impose a far-field shear flow with velocity field $\boldsymbol {u}(x,y,z) = y \boldsymbol {e}_{3}$, where $x,y,z$ are the spatial coordinates in the laboratory frame. We are specifically interested in the resultant motion of the swimmer in the presence of this shear flow.
The dynamics for the orientation of the swimmer frame are given in terms of the Euler angles $(\theta, \psi, \phi )$, which we define formally in Appendix A. Intuitively, one can consider $\theta \in [0, {\rm \pi}]$ to be the pitch, $\psi \mod {2 {\rm \pi}}$ to be the roll and $\phi \mod {2 {\rm \pi}}$ to be the yaw. By combining the classic equations of Jeffery (Reference Jeffery1922) with the appropriate rotation axis in terms of Euler angles, the orientation dynamics is governed by
where the functions $f_i$ encode the hydrodynamic interaction effects of the flow on the swimmer. For a spheroidal swimmer, Jeffery (Reference Jeffery1922) showed that these functions are
where $B = (r^2 - 1)/(r^2 + 1)$ is a constant defined in terms of the spheroidal aspect ratio $r$, and commonly referred to as the Bretherton parameter (Bretherton Reference Bretherton1962). For brevity, we suppress the explicit dependence of the $f_i$ on $B$ unless specifically important. For the chiral swimmers we consider in Reference Dalwadi, Moreau, Gaffney, Ishimoto and WalkerPart 2, the functions (2.2) must be modified appropriately to account for the additional generality in swimmer shape.
Given the symmetries of a spheroidal swimmer, the governing equation for the position of the swimmer in the laboratory frame, $\boldsymbol {X} = X \boldsymbol {e}_{1} + Y \boldsymbol {e}_{2} + Z \boldsymbol {e}_{3}$, is
We emphasise that the simple form of (2.3) conceals the full coupling of the system; $\boldsymbol {V} = V_1\hat {\boldsymbol {e}}_{1} + V_2\hat {\boldsymbol {e}}_{2} + V_3\hat {\boldsymbol {e}}_{3}$ is only constant in the swimmer frame, so the system (2.3) is coupled to the angular dynamics through the dependence of the swimmer basis on the Euler angles, which evolve via (2.1). We also note that the intrinsic symmetries of the spheroid significantly simplify the form of (2.3) so that it does not include any contributions involving the rate of strain. We derive and account for these contributions in Reference Dalwadi, Moreau, Gaffney, Ishimoto and WalkerPart 2, where they are important for the more complex swimmer shapes we investigate therein. As one might expect, those equations for more complex shapes reduce to (2.1)–(2.3) in the appropriate limits.
To gain some intuition for the full system (2.1)–(2.3), it is helpful to briefly consider two limiting cases. The first is the absence of rotational propulsion, i.e. $\boldsymbol {\varOmega }=\boldsymbol {0}$. In this limit, (2.1)–(2.3) reduce to the Jeffery's equations with translation, so the spheroid follows a classic Jeffery's orbit. The other limit is a rotating swimmer (i.e. $\boldsymbol {\varOmega }\neq \boldsymbol {0}$) in quiescent flow (i.e. $f_i = 0$). Then, invoking the helix theorem derived (non-constructively) in Shapere & Wilczek (Reference Shapere and Wilczek1989) for a general swimmer moving periodically in a quiescent flow, the body will exactly follow a helical path in the laboratory frame. The analytic results we derive in this paper allow us to derive this result constructively in Appendix B for the rotating swimmer we consider. We provide explicit expressions for the radius and pitch of the helical trajectory (dependent on $\boldsymbol {\varOmega }$ and $\boldsymbol {V}$) that may be of interest for modelling purposes. The goal of this study is to describe the coupled dynamics arising from the complex interaction between the effects of rotation and flow.
Specifically, we systematically investigate how the well-known Jeffery's orbits (which occur when $\boldsymbol {\varOmega }=\boldsymbol {0}$) should be modified to account for the emergent dynamics of a rapidly rotating active swimmer, which would trace out a helical trajectory in the absence of external flow (which occurs when $f_i = 0$). We are interested in understanding this interaction in the limit where the shear rate of the externally imposed flow is much smaller than the rotation rate of the swimmer, corresponding to $|\boldsymbol {\varOmega }| = |\varOmega _{\parallel }\hat {\boldsymbol {e}}_{1} + \varOmega _{\perp }\hat {\boldsymbol {e}}_{2}| \gg 1$. Additionally, for the translational dynamics, we are interested in the limit where the time scale of translation over the swimmer body radius is much larger than the rotation time scale of the swimmer, corresponding to $|\boldsymbol {\varOmega }| \gg |\boldsymbol {V}|$. These regimes are relevant to many biological and artificial swimmers, specifically those that achieve propulsion through rapid spinning, as noted in § 1 for specific examples.
To investigate this interaction, we examine the dynamics that emerge from the nonlinear, autonomous dynamical system (2.1)–(2.3). We treat $\varOmega _{\parallel } > 0$ without loss of generality, allowing $\varOmega _{\perp }$ to take any real value. The rapid nature of the spinning implies that either $\varOmega _{\parallel }$ or $|\varOmega _{\perp }|$ (typically both) must be large. In the main analysis of this paper (§ 4), we consider rapid generic spinning, i.e. $\varOmega _{\parallel }$, $|\varOmega _{\perp }| \gg 1$ with all other parameters of $\mathit {O}(1)$. Specifically, we treat $\varOmega _{\perp } = \mathit {O}(\varOmega _{\parallel })$ with $\varOmega _{\parallel } \gg 1$ (which will give the same information as taking $\left \lvert \boldsymbol {\varOmega } \right \rvert \gg 1$ with $\alpha = \mathit {O}(1)$), which is a distinguished asymptotic limit of the system, i.e. a general case from which regular asymptotic sublimits can be distilled.
However, since the general analysis is fairly technical, requiring a modification of standard multiple scales analysis, we first consider the sublimit of $\varOmega _{\parallel } \gg 1$ and $\varOmega _{\perp } = \mathit {O}(1)$ in § 3 (equivalent to $\left \lvert \alpha \right \rvert \ll 1$). In contrast to the general case, the multiple scales analysis required for this sublimit is standard, and will serve to demonstrate both our notation and the standard methodology. Additionally, this case of rapid axial rotation aligns with a regime common in bacterial swimming, as many bacteria achieve locomotion through the rapid axial rotation of helical flagellar filaments. Readers familiar with the standard method of multiple scales may wish to skip § 3 and proceed directly to § 4, which captures the full range of behaviours of rapidly spinning swimmers via a more general technical analysis. Additionally, readers more interested in the physical implications of our analysis may wish to skip directly to § 5, where we summarise and discuss our results in terms of their physical relevance.
3. Emergent behaviour in the bacterial sublimit
In this section we consider the sublimit of fast bacterial spinning, corresponding to the system (2.1)–(2.3) for $\varOmega _{\parallel } \gg 1$ and $\varOmega _{\perp } = \mathit {O}(1)$ (equivalent to $\left \lvert \alpha \right \rvert \ll 1$), with all other parameters of $\mathit {O}(1)$ (justified on physical grounds). We animate the full dynamics of this system in supplementary movie 1 available at https://doi.org/10.1017/jfm.2023.923. We analyse this system using the standard method of multiple scales, exploiting large $\varOmega _{\parallel }$. As noted above, this sublimit is mathematically straightforward but introduces notation that will be helpful for the more technical general analysis in § 4.
3.1. Angular dynamics
The translational dynamics (2.3) decouple from the rest of the system. Hence, we first investigate the angular dynamics given by the system (2.1)–(2.2) in the $\varOmega _{\parallel } \gg 1$ and $\varOmega _{\perp } = \mathit {O}(1)$ limit. To this end, we introduce the ‘fast’ time scale $T$ via
and refer to the original time scale $t$ as the ‘slow’ time scale. Proceeding via the standard method of multiple scales (Hinch Reference Hinch1991; Bender & Orszag Reference Bender and Orszag1999), we treat each dependent variable as a function of both the fast- and slow-time scales, converting ordinary differential equations (ODEs) into partial differential equations. The additional degrees of freedom this introduces will be removed later by imposing appropriate periodicity constraints in the fast-time scale.
With the new time scale (3.1), the time derivative becomes
which transforms the system (2.1) into
We expand each dependent variable as an asymptotic series in inverse powers of $\varOmega _{\parallel }$ as follows:
Using the asymptotic expansions (3.4) in the transformed governing equations (3.3), we obtain the leading-order (i.e. $\mathit {O}(\varOmega _{\parallel })$) system
The system (3.5a–c) is straightforward to directly integrate. This will not be the case for the more general problem in § 4, when we have to solve a non-trivial nonlinear problem at leading order. Directly integrating (3.5a–c), we obtain the solutions
Since $\psi$ and $\phi$ are angles interpreted modulo $2 {\rm \pi}$, the leading-order solutions (3.6a–c) are $2{\rm \pi}$-periodic in the fast time $T$. Importantly, $\bar {\vartheta }$, $\bar {\varPsi }$ and $\bar {\varphi }$ are as-of-yet undetermined functions of the slow time, directly related to $\theta _0$, $\psi _0$ and $\phi _0$, respectively. Our remaining goal is to determine the governing equations for these slow-time functions. To do this, we must proceed to the next asymptotic order and determine the solvability conditions.
After posing the asymptotic expansions (3.4) and moving the slow-time derivatives to the right-hand side, the $\mathit {O}(1)$ terms in (3.3) are
now accompanied by fast-time periodicity constraints over any $2{\rm \pi}$ period in $T$. These remove the additional degrees of freedom introduced by initially taking $T$ and $t$ to be independent. For the system (3.7), it is straightforward to determine the requisite solvability conditions. This is because the linear operators acting on $\theta _1$, $\psi _1$ and $\phi _1$ (on the left-hand sides) decouple from one another, even though the full problem (2.1) is fully coupled. This will not be the case for the more general problem in § 4, when we have to derive and solve a non-trivial adjoint problem.
Nevertheless, for the sublimit of bacterial spinning, we consider in this section, this decoupling means that we can straightforwardly obtain our desired solvability conditions by integrating (3.7) with respect to $T$ from $0$ to $2 {\rm \pi}$. The integrated fast-time derivatives on the left-hand side vanish due to the imposed periodicity, and the terms multiplied by $\varOmega _{\perp }$ on the right-hand sides also vanish due to their specific trigonometric forms. Therefore, this integration yields the system
which provide the governing equations we seek for the slow-time functions in the leading-order solutions (3.6a–c). Importantly, we see that the slow-time system (3.8a–c) is equivalent to the full system (2.1) without the rotation terms (i.e. $\varOmega _{\parallel } = \varOmega _{\perp } = 0$) using the equivalence $(\theta, \psi, \phi ) \leftrightarrow (\bar {\vartheta }, \bar {\varPsi }, \bar {\varphi })$. One consequence of this is that $\bar {\varPsi }$ decouples from $\bar {\vartheta }$ and $\bar {\varphi }$, since $\psi$ decouples from $\phi$ and $\theta$ in the absence of rotation. This property is maintained in the more general analysis we present in § 4. Moreover, we note that the Bretherton parameter $B$ remains unchanged in the emergent behaviour in this limit. This property is not maintained in § 4.
To dwell on the general point for clarity, in the sublimit of fast bacterial spinning (with $\varOmega _{\parallel } \gg 1$ and $\varOmega _{\perp } = \mathit {O}(1)$) the emergent behaviour is the same as it would be without any fast spinning. While $\psi$ does vary rapidly due to this fast spinning, the leading-order emergent angular dynamics (up to and including time scales of $t = \mathit {O}(1)$) are not otherwise affected by the fast spinning. Importantly, this emergent indifference to fast spinning is not the case for the more general rotation we consider in § 4.
Figure 2 illustrates these observations through numerical solution of the full and the systematically averaged governing equations, displaying the evolution of the swimmer's orientation in several different cases. In this figure we illustrate the changing orientation of the object by plotting a trajectory on the surface of the unit sphere, which $\theta$ and $\phi$ naturally parameterise via standard spherical coordinates. With the rapid evolution of $\psi$ not illustrated, the particle's orientation traces out a closed orbit on the sphere, corresponding to a Jeffery's orbit that depends on the Bretherton parameter $B$ and the initial conditions. The orbits of passive particles, with $\varOmega _{\perp }=\varOmega _{\parallel }=0$, are illustrated in (a), with varying shape parameters and initial conditions ($\theta (0)={\rm \pi} /18$, $\psi (0)=-{\rm \pi} /4$ in all columns, and $\phi (0)={\rm \pi} /2$ in all but the final column, where $\phi (0) = 0$). We see from (a,b) that the emergent angular dynamics is not affected by the fast spinning as long as $\varOmega _{\perp }$ is not significant, and so active spinning particles generate the same orbits as passive particles if the off-axis spinning is not significant. However, we see from (c) that this result breaks down once $\varOmega _{\perp }$ starts to become a little larger; for an active particle with a significant component of off-axis spinning, the evolution of $\phi$ and $\psi$ differs from that of a passive particle.
In more detail, in figure 2(b) the swimmer spins with $\varOmega _{\parallel }=15$ and $\varOmega _{\perp }=0.5$ and its orientation over time is represented as a thin dark blue line. The trajectory of the leading-order dynamics is shown as a thick green line, and can be seen to be in excellent agreement with the average evolution of the swimmer's orientation; moreover, it is indeed essentially identical to the non-spinning Jeffery orbit from figure 2(a), in line with the conclusions of our asymptotic analysis.
The analysis in this section starts to lose validity when the magnitude of the off-axis rotation $\varOmega _{\perp }$ increases, since this voids the requirement $\varOmega _{\perp } = \mathit {O}(1)$. This can be seen in figure 2(c), where the dynamics predicted by (3.8a–c) start to diverge from the actual emergent behaviour of the swimmer, here with $\varOmega _{\parallel }=15$ and $\varOmega _{\perp }=3$. As before, the blue line shows the full dynamics of the swimmer, and the green line shows the emergent dynamics predicted by the analysis of this section. The latter fails to well approximate the average evolution of the orientation since $\varOmega _{\perp }$ is now significant. In § 4 we generalise our analysis to capture the emergent behaviour outside the bacterial sublimit considered in this section. The predictions of this later analysis are shown as red curves in figure 2(c), which exhibit excellent agreement with the averaged behaviour of the full system.
We have now calculated the emergent behaviour of the angular dynamics in the sublimit of bacterial spinning. We close this problem by considering the emergent behaviour of the translational dynamics.
3.2. Translational dynamics
The governing equations for translation are given in (2.3). Using the transformation (3.2), the translational governing equations (2.3) become
with the $\hat {\boldsymbol {e}}_{i}$ denoting the swimmer frame basis vectors, given in terms of the laboratory frame basis vectors $\boldsymbol {e}_{i}$ in (A1).
Following (3.4) and expanding each dependent variable as an asymptotic series in inverse powers of $\varOmega _{\parallel }$, we obtain the $\mathit {O}(\varOmega _{\parallel })$ (leading-order) system
Hence, the vector position $\boldsymbol {X}$ is independent of the fast time $T$ at leading order. That is, $\boldsymbol {X}_0 = \boldsymbol {X}_0(t)$. Our remaining goal is to derive the dependence of $\boldsymbol {X}_0$ on the slow time via the calculation of solvability conditions at higher order.
At next order (i.e. ${O}(1)$), the translational governing equations (3.9) are
along with fast-time periodicity constraints over any $2 {\rm \pi}$ period in $T$. It is straightforward to obtain the appropriate solvability conditions by integrating (3.11) with respect to the fast time $T$ from $0$ to $2 {\rm \pi}$. In the same way as in § 3.1, the fast-time derivative on the left-hand side vanishes due to the imposition of fast-time periodicity. Substituting the leading-order angular dynamics solutions (3.6a–c) into the explicit forms of $\hat {\boldsymbol {e}}_{i}$ given in (A1), we find that
using the notation $\left \langle \boldsymbol {\cdot } \right \rangle$, defined as the average of its argument over one fast-time oscillation
Additionally, we note that $\tilde {\boldsymbol {e}}_{1}(\bar {\vartheta },\bar {\varphi })$ can be considered equivalent to the (hatted) basis vector $\hat {\boldsymbol {e}}_{1}$ in (A1), but with argument $(\theta,\phi )$ replaced by $(\bar {\vartheta },\bar {\varphi })$.
Hence, using the results (3.12a–c), taking the fast-time average of (3.11) yields
In particular, comparing the slow-time equation (3.14) to the full translational dynamics of (2.3) highlights that the structure of the translational dynamics is broadly unchanged by the rapid spinning, with basis vectors being replaced by their averages. However, this means that the self-induced velocity components in the original $\hat {\boldsymbol {e}}_{2}(\theta,\psi,\phi )$ and $\hat {\boldsymbol {e}}_{3}(\theta,\psi,\phi )$ directions are cancelled out by the fast spinning of the swimmer, so that only the velocity component in the original $\hat {\boldsymbol {e}}_{1}(\theta,\phi )$ direction contributes to the average velocity alongside the shear flow. Hence, this analysis shows that one can effectively neglect the off-axis components of the self-induced propulsive velocity of such a swimmer, at least at leading order and for $t = \mathit {O}(1)$. Interpreted physically, this result establishes that rapid axial spinning, in the absence of significant off-axis rotation, is sufficient to justify modelling the translation of such a swimmer using only the axial component of propulsive linear velocity, since the other self-velocity contributions cancel out at leading order. We show in § 4.5 that this intuitive result requires modification when off-axis rotation is significant.
We show some examples of the swimming behaviours of actively spinning particles in figure 3 and supplementary movie 1. Each panel displays the evolution of a swimmer's position as a black line. A superposed ribbon illustrates the evolution of the intrinsic orientation angle $\psi$; its shade of grey indicates the value of $\psi \in [-{\rm \pi}, {\rm \pi}]$ according to the colourbar on the right-hand side of the figure.
Figure 3(a) showcases the relatively simple dynamics of a particle that is not actively spinning ($\varOmega _{\parallel }=\varOmega _{\perp }=0$) and is only propulsing along its symmetry axis ($V_1 = 1, V_2 = V_3 = 0$), while figures 3(b)–3(f) illustrate the behaviours of spinning particles with various propulsive velocities, having fixed $V_1 = 1$ in each panel and taking $\varOmega _{\parallel }=15$ and $\varOmega _{\perp }=0$ in all but figure 3(a). In each panel, the trajectories predicted by the asymptotic analysis are shown as thick red lines, showcasing excellent agreement with the full dynamics in all but the final case, as expected, wherein the components of propulsive velocity become comparable in magnitude to the fast rate of active rotation.
The near-indistinguishable nature of the red and black lines in panel (b), where $V_2 = V_3 = 0$, highlights that the rapid spinning illustrated by the twisting ribbon need not significantly modify the emergent translational dynamics when there is no off-axis propulsion. In panels (c–e) we allow $(V_2,V_3)\neq\ (0,0)$, which causes far more complex fast-scale dynamics, as can be seen by the cork-screwing of the black line. However, as predicted by our asymptotic analysis, the combination of active spinning and non-trivial self-propulsion does not materially modify the predicted average emergent dynamics (red line). A slight exception to this is shown in figure 3(f), which corresponds to a case where the propulsive velocity is no longer asymptotically smaller than the rate of the rapid rotation. This starts to break the requirement of a separation of scales, giving rise to a discrepancy between the predictions of the asymptotic analysis and the full dynamics.
4. Deriving the general emergent behaviour
We now consider the more general problem, the distinguished asymptotic limit where $\varOmega _{\parallel }$, $|\varOmega _{\perp }| \gg 1$ with $\varOmega _{\perp } = \mathit {O}(\varOmega _{\parallel })$ ($\alpha = \mathit {O}(1)$), and all other parameters are of $\mathit {O}(1)$. We animate the full dynamics of this system in supplementary movie 2. As before, we treat $\varOmega _{\parallel } > 0$ without loss of generality, but allow $\varOmega _{\perp }$ to take any real value. In this distinguished limit the subcases of $|\varOmega _{\perp }| \ll \varOmega _{\parallel }$ and $|\varOmega _{\perp } | \gg \varOmega _{\parallel }$ can be obtained through taking appropriate limits of the results we derive. Given the above, it will be helpful to introduce the notation $\omega = \mathit {O}(1)$ such that
and to formally consider the single asymptotic limit $\varOmega _{\parallel } \gg 1$. We emphasise that $\omega$ can take any real value, and that the angle between the symmetry and rotational axes $\alpha$ is related to $\omega$ through $\tan \alpha = \omega$. Finally, we treat all other parameters as $\mathit {O}(1)$, justified on physical grounds.
We therefore analyse the system (2.1)–(2.3) using a modified method of multiple scales for systems in the limit of large $\varOmega _{\parallel }$, treating $\varOmega _{\perp } = \mathit {O}(\varOmega _{\parallel })$. Our goal is to derive effective equations for the emergent dynamics, in terms of the system parameters. To this end, we introduce the fast-time scale $T$ as
where we use
for notational convenience, and we refer to the original time scale $t$ as the slow-time scale. Under the transformations (4.2), the time derivative becomes
We note that the fast-time notation $T$ in (4.2) is slightly more general than $T$ from the previous section (defined in (3.1)), but that there is a regular limit in which (4.2) reduces to (3.1). Specifically, the bacterial sublimit in the previous section can be obtained formally from the analysis in this section by considering the limit $\omega \to 0$, and hence, $\lambda \to 1$. We proceed as before, treating each dependent parameter as a function of both new time scales. Again, we later remove the additional degrees of freedom this introduces by imposing appropriate periodicity constraints in the fast-time scale.
4.1. Leading-order angular dynamics
For our analysis of the general problem, we start by considering the angular dynamics. The time-derivative transformation (4.4) converts the system (2.1) into
We expand each dependent variable as an asymptotic series in inverse powers of $\varOmega _{\parallel }$, following (3.4). Then the leading-order (i.e. $\mathit {O}(\varOmega _{\parallel })$) version of (4.5) is
The leading-order system (4.6) is forced solely by the rapid spinning of the swimmer, and constitutes a nonlinear 3-D dynamical system for $\boldsymbol {x}_0 := (\theta _0, \psi _0, \phi _0)^{{\rm T}}$, which is notation we will use later.
A key difference in our analysis of the general problem in this section is that the nonlinear, coupled leading-order system (4.6) is not straightforward to solve. This is in direct contrast to the linear, decoupled leading-order system (3.5a–c) of the previous section. One consequence of this is that the important slow-time variables are not immediately obvious.
Of great help in our task is the fact that we can solve the leading-order system (4.6) in closed form. This is abetted by two key observations. Firstly, (4.6c) decouples from the other two subequations, and hence, we treat it last. Secondly, (4.6a) and (4.6b) can be combined and integrated to admit a first integral, which constitutes one of the three slow-time functions we expect to arise. To see this, we first divide (4.6b) by (4.6a) and multiply through by $\omega \cos \psi _0$ to yield
Since the left-hand side of (4.7) can be written as $\omega \partial (\sin \psi _0)/\partial \theta _0$, (4.7) can be rewritten as a first-order linear differential equation for $\sin \psi _0$ as a function of $\theta _0$, i.e. implicitly making the change of variables $(T,t) \mapsto (\theta _0,t)$. Multiplying through by the appropriate integrating factor of $\sin \theta _0$ and grouping everything on the same side, we may rewrite (4.7) as
Hence, the required first integral is obtained by directly integrating (4.8) to yield
thus generating the function of integration $H(t)$, emphasising that $H$ is independent of the fast time $T$. As will become clear once we complete the full analysis, it will turn out to be particularly convenient to redefine $H(t) = \lambda \cos \bar {\vartheta }(t)$, so we formally write
where $\bar {\vartheta } = \bar {\vartheta }(t)$ is the first of three slow-time functions (each marked with an overbar) we obtain from our leading-order analysis in this section. We note that our definition of $\bar {\vartheta }$ in (4.10) in the limit $\omega \to 0$ coincides with its former definition in (3.6a). The goal of the next-order analysis in § 4.2 is to derive the governing equations satisfied by $\bar {\vartheta }$ and the two other appropriate slow-time functions we derive in the remainder of this section.
Substituting the fast-time conserved quantity (4.10) into (4.6a), we obtain the following nonlinear first-order differential equation:
Since (4.11) is an autonomous first-order differential equation for $\theta _0$ as a function of $T$, it admits a solution in the implicit form
where $\bar {\varPsi } = \bar {\varPsi }(t)$ is the second of the three slow-time functions we obtain from our leading-order analysis. In fact, we can evaluate (4.12) explicitly using the substitution $\lambda \cos s = \cos \bar {\vartheta } + u \sin \bar {\vartheta }$ to convert (4.12) into
omitting the transformed limit of the integral for notational convenience. It is then straightforward to integrate the left-hand side of (4.13) and rearrange to deduce that
We note that the $\mp$ in (4.13) (as well as the choice of $\cos (T +\bar {\varPsi })$ in (4.14a) instead of $\sin (T +\bar {\varPsi })$) can be absorbed into an initial phase shift in $\bar {\varPsi }$. Therefore, we have taken the negative root in moving from (4.13) to (4.14a) for later convenience, without loss of generality. Rearranging (4.14a), we deduce that
It will also be helpful later to note that the combination of (4.10) and (4.14a) yields the relationship
and that differentiating (4.14a) with respect to $T$ and invoking (4.6a) yields the relationship
Given the relationships (4.14), we see that the trigonometric quantities $\cos \theta _0$, $\sin \theta _0 \sin \psi _0$ and $\sin \theta _0 \cos \psi _0$ are fast oscillations in time whose mean, amplitude and phase are slowly varying, coupled quantities.
The final result required from the leading-order system is the solution of (4.6c), and the subsequent appearance of the third (and final) slow-time function. Rewriting the right-hand side of (4.6c) as $\sin \theta _0 \sin \psi _0/\sin ^2 \theta _0$, we can use the expressions (4.14b)–(4.14c) to rewrite (4.6c) as
Then, observing that
integrating (4.15) directly with respect to the variable $\omega \sin (T + \bar {\varPsi }) / (\omega \cos \bar {\vartheta } \cos (T + \bar {\varPsi }) + \sin \bar {\vartheta })$ yields the solution
where $\bar {\varphi } = \bar {\varphi }(t)$ is the final slow-time function we obtain from our leading-order analysis. For later convenience, it is also helpful to expand out the left-hand side of (4.17) and rearrange to obtain the relationship
We have now fully solved the nonlinear, coupled leading-order problem (4.6) to obtain the fast-time solutions (4.14) and (4.18), in terms of the three slow-time functions of integration $\bar {\vartheta }(t)$, $\bar {\varPsi }(t)$ and $\bar {\varphi }(t)$. Importantly for the method of multiple scales, the leading-order solutions given in (4.14) and (4.18) are $2{\rm \pi}$-periodic in the fast time $T$ (appropriately interpreting $\psi _0$ and $\phi _0$ modulo $2 {\rm \pi}$). (Despite starting with a nonlinear oscillator system, we are able to use multiple scales instead of the more general method of Kuzmak (Reference Kuzmak1959) because the period of the fast-time oscillation is independent of the slow time.) The slow-time functions $\bar {\vartheta }(t)$, $\bar {\varPsi }(t)$ and $\bar {\varphi }(t)$ are currently undetermined, and the goal of our next-order analysis in § 4.2 is to derive the governing equations for these functions. For generic initial conditions of the full system, one can also derive the appropriate initial conditions for these slow-time functions. We present this derivation in Appendix C, noting that some care is required in order to ensure that the correct branches are taken.
Over the slow time, one can think of $\bar {\vartheta }$ as controlling some emergent amplitude or mean of oscillation, $\bar {\varPsi }$ as controlling some emergent phase of oscillation and $\bar {\varphi }$ as an emergent drift in yawing. Additionally, in a sense to be made precise later, we can interpret each slow-time function as being associated with an underlying variable. Specifically, and in the same way as their equivalent variables in § 3, we can associate $\bar {\vartheta }$ with $\theta$, $\bar {\varPsi }$ with $\psi$ and $\bar {\varphi }$ with $\phi$.
4.2. Solvability conditions from the next-order system
Our remaining goal is to determine the governing equations satisfied by the slow-time functions $\bar {\vartheta }(t)$, $\bar {\varPsi }(t)$ and $\bar {\varphi }(t)$. We note that the procedure to do this is non-standard for our problem, as it involves a non-trivial system of equations.
Specifically, when using the method of multiple scales for standard problems involving first-order ODEs, procedures to obtain governing equations for the slow-time functions often involve solving the next-order system of the general form
where $\boldsymbol {x}_1$ is periodic in $T$ (say $2 {\rm \pi}$-periodic, without loss of generality), $L$ is a differential-algebraic linear operator (generally dependent on $\boldsymbol {x}_0(T,t)$) and $\boldsymbol {F}$ is a function of both $\boldsymbol {x}_0$ and its slow-time derivative. Then one proceeds by visually identifying secular terms in the solution and forcing their coefficients to vanish, generating the appropriate slow-time governing equations (Hinch Reference Hinch1991; Bender & Orszag Reference Bender and Orszag1999). However, this procedure can be challenging for more complicated systems where full solutions of the next-order system are prohibitively difficult to obtain, or where secular terms are difficult to identify.
In such systems it can often be simpler to obtain the requisite governing equations by deriving solvability conditions for the next-order system instead, through the imposition of periodicity on the fast scale. This can be simpler because it circumvents the need to solve an inhomogeneous problem. The appropriate way to derive these solvability conditions is via the method of multiple scales for systems (see, for example, pp. 127–128 of Dalwadi (Reference Dalwadi2014) or p. 22 of Dalwadi et al. Reference Dalwadi, Chapman, Oliver and Waters2018), which follows from the Fredholm alternative theorem applied to a periodic system. To briefly summarise here: a next-order linear system such as (4.19) yields appropriate slow-time evolution equations via its solvability condition
where $\boldsymbol {P}$ is the fast-time periodic solution of the homogeneous adjoint problem
with fast-time periodicity conditions, where $L^*$ is the differential-algebraic linear adjoint operator. The solvability conditions (4.20) will generate the appropriate slow-time governing equations, reducing our task to explicitly solving the homogeneous 3-D system (4.21) instead of the inhomogeneous 3-D system (4.19). Generally, each linearly independent solution of the adjoint problem (4.21) will contribute one solvability condition.
For the specific problem we consider, the next-order system (4.19) constitutes a 3-D dynamical system for $\boldsymbol {x}_1(T,t) := (\theta _1, \psi _1, \phi _1)^{{\rm T}}$, and is generated from the $\mathit {O}(1)$ terms in (4.5) after posing the asymptotic expansions (3.4). Specifically, the linear operator $L$ and the right-hand side $\boldsymbol {F}$ are
where $\partial _{T} \equiv \partial /\partial T$. The homogeneous adjoint operator $L^*$ is the transpose of the matrix operator taking the adjoint of each element, and is therefore defined as
We emphasise that $L \neq L^*$, i.e. this problem is not self-adjoint.
To generate the slow-time equations we are seeking for the three slow-time functions $\bar {\vartheta }(t)$, $\bar {\varPsi }(t)$ and $\bar {\varphi }(t)$, we must derive the appropriate solvability conditions (4.20) by solving the homogeneous adjoint problem (4.21) with $L^*$ defined in (4.24). Given that the coefficients in (4.24) depend on the fast time $T$, this is a non-trivial task. Nevertheless, we are able to solve (4.21) in Appendix D, where we deduce that the general solution is
for arbitrary constants $C_1$, $C_2$ and $C_3$.
Therefore, we obtain our required solvability conditions by substituting the adjoint solutions (4.25) into the general solvability condition (4.20), with $\boldsymbol {F}$ defined in (4.23), and setting the resulting coefficients of $C_i$ ($i = 1, 2, 3$) to zero. This procedure yields the following three solvability conditions:
Here we use the subscript $t$ to denote partial differentiation with respect to $t$, and we use the notation $\left \langle \boldsymbol {\cdot } \right \rangle$ to denote the average of its argument over one fast-time oscillation, as defined in (3.13).
Our remaining task is to evaluate the averages in (4.26) in terms of the slow-term functions $\bar {\vartheta }$, $\bar {\varPsi }$ and $\bar {\varphi }$, and this will lead us to the emergent slow-time equations for the system. The left-hand sides of (4.26) are currently written in terms of the original variables $\theta _0(T,t)$, $\psi _0(T,t)$ and $\phi _0(T,t)$. To proceed, we need to write them in terms of the slow-time functions $\bar {\vartheta }(t)$, $\bar {\varPsi }(t)$ and $\bar {\varphi }(t)$ and then take appropriate averages. We carry out this task in Appendix E, which allows us to rewrite (4.26) as
Equations (4.27) represent the appropriate solvability conditions for any $f_i$ acting on a rapidly spinning object in Stokes flow. In the remainder of Part 1, we evaluate these solvability conditions for the $f_i$ defined in (2.2), which represent the hydrodynamic interaction of shear flow with a spheroidal object. In Reference Dalwadi, Moreau, Gaffney, Ishimoto and WalkerPart 2 we consider the interaction of shear flow with general chiral objects with helicoidal symmetry. This requires evaluating a different set of functions on the right-hand side of (4.27), resulting in fundamentally different emergent equations to Part 1.
4.3. Evaluating the right-hand side of the solvability conditions
Our final task is to evaluate the right-hand side of the three solvability conditions in (4.27), using the $f_i$ defined in (2.2). While (4.14) provide helpful expressions for the required trigonometric functions of $\theta _0$ and $\psi _0$, currently our only expression for $\phi _0$ is (4.18), given in terms of $\tan \phi _0$. Given that the $f_i$ defined in (2.2) involve $\sin 2 \phi _0$ and $\cos 2 \phi _0$, it is helpful to derive expressions for these quantities. To do this, we first calculate expressions involving $\sin \phi _0$ and $\cos \phi _0$ from $\tan \phi _0$ defined in (4.18), i.e.
using the shorthand notation
for algebraic convenience. The factors of $\lambda \sin \theta _0$ in (4.28) arise from imposing the constraint $\sin ^2 \phi _0 + \cos ^2 \phi _0 = 1$ and using (4.14b) to evaluate $\sin ^2 \theta _0$. While it is not essential, since the sign differences would cancel with one another in our subsequent analysis, for definiteness, we choose the positive square root, essentially exploiting the fact that $\sin \theta _0 > 0$ (since $\theta _0 \in [0, {\rm \pi}]$), and using our choice of $\bar {\varphi }(0)$ up to an additive multiple of ${\rm \pi}$ (see Appendix C). Then, from the expressions (4.28), we may use appropriate double-angle formulae to deduce that
where we have introduced the functions
To evaluate the fast-time averages on the right-hand side of (4.27), it is helpful to exploit the parity of various quantities. In particular, we note that $\mathcal {C}$ and $\mathcal {S}$ are even and odd, respectively, around $\sigma = {\rm \pi}$, which follows immediately from their definitions (4.30c)–(4.30d). Since $\sin ^2 \theta _0$ is also even, as follows from (4.14b), this means that $\cos 2 \phi _0$ and $\sin 2 \phi _0$ defined in (4.30a)–(4.30b) are naturally decomposed into their odd and even parts. Additionally, since the expressions in (4.14) tell us that $\cos \theta _0$ and $\sin \theta _0 \sin \psi _0$ are even, and that $\sin \theta _0 \cos \psi _0$ is odd (each around $\sigma = {\rm \pi}$), several terms on the right-hand side of (4.27) will vanish immediately under fast-time averaging. Exploiting this parity, we find that
To evaluate the required average in (4.31a) and (4.31b) (which is the same), we use the relationships (4.14), (4.30c)–(4.30d), and simplify the resulting expressions to deduce that
Similarly, the required average in (4.31c) can be evaluated using (4.14b) and (4.30c)–(4.30d), to obtain
Then using the results (4.32) in (4.31), we can write
Equation (4.33) concludes our evaluation of the right-hand sides of the solvability conditions. The next step is to piece this all together to obtain the emergent angular dynamics.
4.4. The emergent angular dynamics
We now have the requisite information to obtain the governing equations for $\bar {\vartheta }$, $\bar {\varPsi }$ and $\bar {\varphi }$, which are the main results we have been seeking through this analysis. To this end, we substitute the right-hand side results (4.33) into the solvability conditions (4.27), then rearrange to obtain the following nonlinear system:
Here $\bar {\vartheta }$, $\bar {\varPsi }$ and $\bar {\varphi }$ feed into the full leading-order asymptotic solution through (4.14) and (4.18).
Remarkably, (4.34) can be rewritten in terms of the functions $f_i$, which are defined in (2.2), to yield
where we define $\hat {B}$ to be the effective Bretherton parameter, given explicitly by
Writing the system (4.34) in the form (4.35) leads to several key deductions. Firstly, by comparison with the full system (2.1)–(2.2), we see that we can formally identify each slow-time function with an underlying variable: $\bar {\vartheta }$ with $\theta$, $\bar {\varPsi }$ with $\psi$, and $\bar {\varphi }$ with $\phi$. Moreover, the compact relationship described by (4.35) retroactively motivates our notational choice to define the effective fast-time conserved quantity as $\lambda \cos \bar {\vartheta }(t)$ in (4.10) instead of $H(t)$ as in (4.9).
Secondly, the effective Bretherton parameter is given in (4.35b), and plotted with respect to $\omega$ in figure 4(a), from which we note that the sign of the effective Bretherton parameter is different to that of the actual Bretherton parameter if $|\omega | > \sqrt {2}$. Moreover, the effective Bretherton parameter vanishes (i.e. the effective particle behaves as a sphere) if $|\omega | = \sqrt {2}$. More generally, since (4.35) is exactly equivalent to Jeffery's equations, we have found that the rapid spinning of spheroidal swimmers generates Jeffery's orbits with an effective Bretherton parameter that we have calculated analytically. Therefore, the effect of rapid spinning is to modify the effective shape of the spheroidal swimmer, as quantified through the relationship (4.35b).
Finally, it is of interest to briefly consider the limiting cases of $\omega \to 0$ and $|\omega | \to \infty$. We have already considered the limit of $\omega \to 0$ (where the rotation axis coincides with the spheroidal symmetry axis) explicitly in § 3. It is reassuring to note that in this limit, our generalised effective Bretherton parameter $\hat {B} \to B$ from (4.35b), and hence, the system (4.35) reduces to
in agreement. That is, this type of spinning leaves the system essentially unchanged, recapturing the explicit results of § 3.
In the limit of $|\omega | \to \infty$ (where the rotation axis is perpendicular to the spheroidal symmetry axis), our generalised effective Bretherton parameter $\hat {B} \to -B/2$ from (4.35b), and the system (4.35) reduces to
That is, the sign of the effective Bretherton parameter changes (so an active oblate spheroid behaves as a passive prolate spheroid, and vice versa) and its magnitude halves (making it behave more like a sphere, in some sense).
We illustrate our asymptotic results in figure 5, showing the rotational dynamics of the swimmer orientation $(\theta,\phi )$ via the same representation on the unit sphere as in figure 2. With figure 5(a) showing the Jeffery's orbits of passive particles for reference, it is clear that rapid spinning with non-zero $\omega$ can significantly impact on the emergent dynamics, with $\omega$ and the Bretherton parameter $B$ varying between panels. Within each column, the initial configuration is fixed and shown as a black dot for reference. The average trajectory predicted by the asymptotic analysis is shown as a red curve in each plot, and the full numerical solution is shown as a blue curve. In each case, the full asymptotic solution obtained from combining the leading-order fast-time solutions (4.14) and (4.18) with the slow-time equations (4.35) successfully captures the leading-order behaviour of the full oscillatory dynamics, even when it is markedly complex, though we omit this comparison in figure 5 to allow the details of the trajectories to be seen cleanly.
As predicted by (4.35), the emergent dynamics of the spinning objects are simply Jeffery's orbits in transformed variables with modified Bretherton parameters. That is, a rapidly spinning object behaves as though it is a differently shaped spin-free object. In particular, with $\left \lvert \hat {B} \right \rvert < \left \lvert B \right \rvert$, the emergent Jeffery's orbits are somewhat simpler, in that they more closely resemble the straightforward behaviour exhibited by a sphere. This is particularly evident as $\left \lvert \omega \right \rvert \rightarrow \sqrt {2}$, with the dynamics approaching those of a sphere for all particles irrespective of their shape, as shown in figure 5(d).
Our final task is to use our results for the emergent angular dynamics to calculate the emergent translational dynamics.
4.5. The emergent translational dynamics
In this final subsection for the general problem, we calculate the emergent translational dynamics, for which the full governing equations are given in (2.3). Our aim is to derive emergent (effective) governing equations for the translation, again using the method of multiple scales. Using the time-derivative transformation (4.4), the governing equations (2.3) become
We proceed as before, taking an asymptotic expansion of the dependent variables in inverse powers of $\varOmega _{\parallel }$, as in (3.4). At leading order (i.e. ${O}(\varOmega _{\parallel })$), we obtain
Therefore, in the same way as for the bacterial sublimit in § 3, $\boldsymbol {X}_0 = \boldsymbol {X}_0(t)$ and the vector position $\boldsymbol {X}$ is independent of the fast time $T$ at leading order.
At next order (i.e. ${O}(1)$), substituting the leading-order fast-time angular dynamics solutions (4.14) and (4.18) in (4.38), we obtain the system
with a fast-time periodicity constraint on $\boldsymbol {X}_1$. The solvability conditions are straightforward to obtain by integrating (4.40) with respect to the fast time $T$ from $0$ to $2 {\rm \pi}$. As before, the fast-time derivative vanishes due to periodicity. In Appendix F we calculate the fast-time averages of the unit vectors as
where $\tilde {\boldsymbol {e}}_{1}(\bar {\vartheta },\bar {\varphi })$ can be considered equivalent to the (hatted) basis vector $\hat {\boldsymbol {e}}_{1}$ in (A1), but with argument $(\theta,\phi )$ replaced by $(\bar {\vartheta },\bar {\varphi })$.
Hence, substituting (4.41a–c) into the fast-time averaged version of (4.40), the emergent governing equation for translation is
Of note, and generalising the case explored in § 3, we now retain a contribution from an off-axis direction of swimming, namely $V_2$. That is, this off-axis swimming component does not average to zero in the general case. The evolution of the effective swimming velocity $\hat {V}$ (normalised by $|V| = \sqrt {V_1^2+V_2^2}$) with respect to $\omega$ is shown in figure 4(b) for different values of off-axis velocity. Perhaps surprisingly, the coefficient of $V_2$ in (4.42) is odd in $\omega$, suggesting that the behaviour of a swimmer with $\omega <0$ can be markedly different to that of one with $\omega >0$, recalling that the sign of $\omega$ is determined solely by the relative directions of spinning via (4.1). Thus, supposing that $V_1=\left \lvert \omega \right \rvert V _2$, a swimmer with both $\varOmega _{\perp }>0$ and $\varOmega _{\parallel }>0$ propels itself at an effective speed of $2\left \lvert V_1 \right \rvert /\lambda$, while swapping either one of the directions of spinning gives an effective swimming speed of zero. The role of $V_2$ is further highlighted in the limit $\left \lvert \omega \right \rvert \to \infty$, in which the contribution of $V_1$ vanishes. Hence, in this limit, which might be considered the ‘opposite’ to the bacterial case, the translational dynamics are only modulated by $V_2$ and the shape-independent advection by the shear.
The conditional significance of $V_2$ is illustrated in figures 6, 7 and supplementary movie 2, which explore the emergent translational dynamics more generally. Figures 6(a) and 6(b) highlight how $V_3$ has no effect on the leading-order translation, with the average trajectories in each panel being qualitatively indistinct, while we see an overall trend that increasing $\omega$ decreases the effective swimming speed when $V_2=0$, as predicted by our analytic result (4.42). In figures 6(c) and 6(d), however, we see that increasing $V_2$ from zero does modify the swimming trajectory, with these modifications increasing in significance as $\left \lvert \omega \right \rvert$ increases, in line with the predictions of (4.42). We emphasise that $\hat {V} \to V_2$ as $\left \lvert \omega \right \rvert \to \infty$, in line with physical intuition, despite the apparent importance of the sign of $\omega$. This is because $\tilde {\boldsymbol {e}}_{1}(\bar {\vartheta },\bar {\varphi })$ depends on $\omega$ and, in fact, tends to exact opposite directions in the limits $\omega \to \pm \infty$ (when the initial conditions derived in Appendix C are taken into account). Figure 7 further emphasises the significant differences that can be observed on the translational dynamics upon switching the sign of $\omega$.
5. Results and conclusions
In this section we summarise the main results and conclusions that can be drawn from our detailed asymptotic analysis.
5.1. Dynamics in the bacterial limit
In the ‘bacterial’ limit, fast rotation occurs only about the spheroidal axis of symmetry. In § 3.1, which corresponds to the limit $\omega \to 0$ in the general case considered in § 4, we showed that rapid rotation in the bacterial limit does not materially impact on the emergent angular dynamics. In particular, the leading-order dynamics of all but the spin angle $\psi$ are exactly as in the rotation-free case, so that rapid rotation leaves the direction of the swimmer symmetry axis unchanged. Furthermore, if we consider the leading-order average dynamics, then the spin angle evolves precisely as if there were no rapid rotation. Hence, the average orientation of the swimmer is not altered by rapid spinning about the symmetry axis, and hence, one can safely neglect rapid rotation without consequence when considering leading-order angular dynamics in the bacterial limit.
However, one must be careful with the effect of rapid spinning in this limit when concerned with the translational dynamics of the swimmer. In § 3.2 we showed that the leading-order dynamics of self-propulsion are dominated by the component of velocity along the axis of symmetry. Physically, this is because the other components of linear velocity cancel out over a period of the fast axial rotation, so one can safely neglect off-axis components of translation at leading order. Moreover, the effective direction of self-propulsion follows the average orientation of the symmetry axis, which evolves as if there were no rapid spinning.
Thus, overall, our analysis demonstrates that a spheroid that spins rapidly about its symmetry axis can be reliably modelled as a non-spinning version of the same spheroid in a shear flow. Moreover, we have shown that any relatively slow, off-axis rotation can be neglected, as can any components of translational velocity that are not aligned with the axis of symmetry. In the context of the idealised bacterial swimmers that motivated this limit, our analysis supports the use of simple, effective models of bacterial swimmers that do not include any details of off-axis rotation or off-axis translation, which might plausibly arise due to the finite length of a propulsive flagellum or fascicle, for instance.
5.2. Dynamics for general rapid spinning
For general rapid spinning, the axis of rotation does not coincide with the symmetry axis, and the straightforward simplifications of the bacterial limit no longer hold. Nevertheless, from our analysis in § 4 for $\omega = \mathit {O}(1)$, we are still able to characterise and quantify the emergent behaviours that arise from the multiscale dynamics, generalising those seen in the bacterial case. For a more intuitive interpretation of our results, it is helpful to work with $\alpha$, the angle between the rotational and symmetry axes, instead of $\omega$. These quantities are related through (the principal branch of) $\tan \alpha = \omega$.
A key result of our general analysis is that rapid spinning of spheroidal objects always leads to effective Jeffery's orbits, though now in transformed variables that can be associated with the original angular variables. In addition, we have analytically calculated the parameter that characterises these orbits: the effective Bretherton parameter $\hat {B}$, defined in (4.35b), which depends nonlinearly on $\alpha$. This effective parameter $\hat {B}$ can be sufficiently different from the original Bretherton parameter $B = (r^2 - 1)/(r^2 + 1)$ (where $r$ is the length-to-diameter ratio of the spheroid) so as to completely alter the character of the orbit. Associating the effective Bretherton parameter $\hat {B}$ with an effective aspect ratio $\hat {r}$ via $\hat {B} = (\hat {r}^2 - 1)/(\hat {r}^2 + 1)$, and using our key result for $\hat {B}$ in (4.35b), we may deduce an explicit expression for the effective aspect ratio:
We illustrate the effective aspect ratio $\hat {r}$ in figure 8 as a function of $\alpha$ and $r$, noting that $\hat {r}$ is even in $\alpha$ (as one might expect intuitively). Our results show that rapid rotation modifies the effective hydrodynamic shape of the spheroid through its effective aspect ratio and the natural parameterisation of its orientation. In supplementary movie 3 we show how the average rotational dynamics of a rapidly spinning spheroid are well represented by the rotational dynamics of a passive spheroid with effective Bretherton parameter (4.35b) (which corresponds to an effective aspect ratio (5.1)).
We see that $\hat {r} = r$ for $\alpha = 0$, in agreement with the results of § 3 summarised in § 5.1. As $\alpha$ increases, the effective aspect ratio $\hat {r}$ changes monotonically, with the effective shape of prolate spheroids ($r > 1$) becoming less prolate and eventually oblate, and vice versa for oblate spheroids ($r < 1$). Thus, the distortion of the effective shape is most pronounced when $\alpha = {\rm \pi}/2$, i.e. when the rotation and symmetry axes are perpendicular. In this scenario, we see from (5.1) that $\hat {r}^2 = (3 + r^2)/(1 + 3 r^2)$. In particular, this tells us that an extremely prolate spheroid rapidly rotating along an axis perpendicular to its symmetry axis will behave as an oblate spheroid with effective aspect ratio $1/\sqrt {3}$, no matter the original aspect ratio of the spheroid. Similarly, an extremely oblate spheroid will behave as a prolate spheroid with effective ratio $\sqrt {3}$. Hence, the effective hydrodynamic shape is not simply the envelope of the rotating shape.
A particularly striking result of our analysis is the scenario in which $\tan \alpha = \sqrt {2}$, which leads to an effective aspect ratio of $\hat {r} = 1$ for all values of the original aspect ratio $r$. That is, if the rotational axis is at an angle $\alpha = \arctan \sqrt {2} \approx {54.7}^{\circ }$ from the symmetry axis of a spheroid, the spinning spheroid will behave hydrodynamically as an effective sphere in the shear flow. This universal reduction can be seen in figure 8. It can also be seen in figure 5(d), with the orbits of differently shaped spheroids collapsing onto the same dynamics for $\omega = \sqrt {2}$ (i.e. $\alpha = \arctan \sqrt {2}$). This identification of an active spheroidal particle behaving like a sphere in 3-D shear flow complements recent work in which a passive non-spheroidal particle (a symmetric boomerang with angle $\arccos (2 \sqrt {3} - 3) \approx 62.3^{\circ }$) was shown to behave like a sphere in two-dimensional (2-D) shear flow (Roggeveen & Stone Reference Roggeveen and Stone2022).
Now considering the translational dynamics, we again see that effective propulsion occurs along the average symmetry axis. However, in this general case, from (4.42) the effective speed of propulsion is $\hat {V} = V_1 \cos \alpha + V_2 \sin \alpha$. This gives an intuitive explanation for the effective speed of propulsion; it is exactly the magnitude of the projection of $\boldsymbol {V}$ along the rotational axis $\boldsymbol {\varOmega }$. That is, $\hat {V} = (\boldsymbol {V} \boldsymbol {\cdot } \boldsymbol {\varOmega })/|\boldsymbol {\varOmega }|$. This means that the effective speed is now influenced by both the on-axis speed $V_1$ and the component of propulsion in the direction of off-axis spinning $V_2$, with their relative contributions modulated by their relative contributions along the rotational axis. In addition, as $\alpha \rightarrow {\rm \pi}/2$, the off-axis component of propulsion $V_2$ dominates the on-axis propulsion $V_1$ (though the maximal speed of self-propulsion is always bounded above by $\sqrt {V_1^2 + V_2^2}$). Hence, in contrast to the bacterial limit, off-axis propulsion cannot simply be neglected for spheroidal objects undergoing off-axis spinning; here, the effects of propulsion parallel to the direction of off-axis spinning can combine with the off-axis spinning itself to contribute to the effective propulsion.
Indeed, the relative signs of the translational and angular components of propulsion can significantly alter the emergent behaviour. For example, if $V_1\approx -V_2$, then spinning with $\alpha = {\rm \pi}/4$ can lead to propulsion having essentially no leading-order contribution. Increasing $\alpha$ beyond this value highlights that the direction of effective propulsion can strongly depend on $\alpha$ and, therefore, on the relative direction of the rotational axis in comparison to the symmetry axis. In particular, it is possible to be in a scenario where the reversal of precisely one of the components of rotation gives rise to both a reversal and a change of magnitude of the effective propulsive velocity.
6. Discussion
In this first of two studies, we have explored the behaviours of rapidly spinning spheroidal objects in 3-D Stokes shear flow. Using the method of multiple scales for systems, we have derived effective governing equations that, when written in terms of appropriately transformed variables, are simply the classic Jeffery's equations (Bretherton Reference Bretherton1962; Jeffery Reference Jeffery1922), with appropriately modified Bretherton parameter $\hat {B}$. This result is in a similar spirit to Bretherton (Reference Bretherton1962) and Brenner (Reference Brenner1964), who showed that Jeffery's equations hold for the motion of general axisymmetric objects, with the Bretherton parameter being interpreted as a measure of an effective aspect ratio. Indeed, our work shows that Jeffery's equations also hold for active axisymmetric (not just spheroidal) objects. That is, the rotational dynamics of rapidly spinning axisymmetric objects also behave as effective passive spheroids in shear flow, with a corresponding effective aspect ratio that depends on the particular object, and is not just the envelope of the spinning object. Through our analysis, we analytically calculate how the rapid spinning modifies this effective parameter in the effective equations, and hence, the effective aspect ratio (and, therefore, shape) of the spheroid. This identification of modified, effective parameterisations of the original 3-D governing equations is similar to those found in other recent applications of multi-time-scale analysis to planar swimming problems (Walker et al. Reference Walker, Ishimoto, Moreau, Gaffney and Dalwadi2022a,Reference Walker, Ishimoto, Gaffney, Moreau and Dalwadib), suggesting that such reductions may be commonplace in similarly posed problems of Stokes flow.
In two opposite limits of the rapid spinning problem, where spinning is dominated by either the on- or off-axis component, we obtain substantial simplifications to our general results. The case where spinning about the symmetry axis is dominant can be considered as a simplified model of a rotating bacterial swimmer. (Acknowledging that much of the detail of the complex swimming problem has been neglected in this idealised study of rigid bodies.) In this bacterial limit, our results show that one can effectively ignore any off-axis propulsion and all components of rapid spinning, so that the swimmer can be effectively modelled as a non-spinning object with axial propulsion. In the other limit, where the dominant spinning is perpendicular to the symmetry axis, while the angular dynamics now exhibit strong (leading-order) oscillations, the long-term angular dynamics is still governed by effective Jeffery's equations in terms of appropriate transformed variables. Moreover, the effective propulsion is still along the ‘average’ symmetry axis, as in the bacterial case, though the effective propulsive velocity is now dominated by the off-axis component of instantaneous propulsion, perhaps counterintuitively.
In Reference Dalwadi, Moreau, Gaffney, Ishimoto and WalkerPart 2 of this two-part study we will consider swimmers of a more general shape. Specifically, we will study helicoidal swimmers with chirality, generalising our results beyond the spheroidal swimmers considered in Part 1. As one might expect, this additional generality complicates the analysis we have presented here. Nevertheless, we will see that it remains analytically tractable via the techniques employed in this work.
Supplementary material
Supplementary movies are available at https://doi.org/10.1017/jfm.2023.923. Minimal computer code for exploring the dynamics, as well as the scripts used to generate the figures in this study are available at https://github.com/Clementmoreau/spinningswimmers.
Funding
M.P.D. is supported by the UK Engineering and Physical Sciences Research Council (grant no. EP/W032317/1). C.M. is a JSPS Postdoctoral Fellow (P22023) and acknowledges support by the JSPS-KAKENHI Grant-in Aid for JSPS Fellows (grant no. 22F22023). K.I. acknowledges JSPS-KAKENHI for Young Researchers (grant no. 18K13456), JSPS-KAKENHI for Transformative Research Areas (grant no. 21H05309), JST, PRESTO (grant no. JPMJPR1921) and JST, FOREST (grant no. JPMJFR212N). B.J.W. is supported by the Royal Commission for the Exhibition of 1851.
Declaration of interests
The authors report no conflict of interest.
Appendix A. Defining the Euler angles
Here, we derive the Euler angles used in the main text. We relate the bases of the laboratory frame and the swimmer frame by an $xyx$-Euler angle transformation. In particular, with $c_\theta$, $s_\theta$ denoting $\cos \theta$, $\sin \theta$, and similarly for other angles, we have
as illustrated in figure 9. Here, $\boldsymbol {A}$ is the matrix for a rotation of angle $\phi$ about $\boldsymbol {e}_{1}$, $\boldsymbol {B}$ is the matrix for a rotation of angle $\theta$ about the $y$ axis of the resulting reference frame, with associated basis vector
and $\boldsymbol {C}$ is the matrix for a rotation of angle $\psi$ about the body fixed axis, $\hat {\boldsymbol {e}}_{1}$. Intermediate basis vectors are defined in figure 9.
The inverse of this mapping is generated by applying the rotations in reverse order with negated angles. Thus, by applying the transformation $(\phi,\theta,\psi )\mapsto (-\psi,-\theta,-\phi )$ to the matrix in (A1), we have
Furthermore, the Euler angle transformation also gives the relation between the angular velocity of the swimmer frame in the presence of flow, denoted $\boldsymbol {\varOmega }^f$, and the time derivatives of the Euler angles via
which simplifies to
Appendix B. Helical trajectories in the absence of shear flow
In this appendix we discuss the helical trajectory of a spinning microswimmer in the absence of any background flow. We rederive the well-known helix theorem of Shapere & Wilczek (Reference Shapere and Wilczek1989), demonstrating that a swimmer with constant swimming velocity $\boldsymbol {V}$ and angular velocity $\boldsymbol {\varOmega }$ follows a helical trajectory. In particular, we provide a relation between the swimming velocity and the properties of the helix, such as its radius and pitch.
In the absence of any background flow, there is a single time scale present in the angular dynamics, so that they are in fact captured exactly by (4.6), interpreted as ODEs in time $T$ and with $\bar {\vartheta }$, $\bar {\varPsi }$ and $\bar {\varphi }$ being constants. For clarity, we denote the solutions to these helical angular dynamics by ${\theta }_h$, ${\psi }_h$ and ${\phi }_h$. Since there is no background flow, without loss of generality we may set the initial orientation of the swimmer to be ${\theta }_h(0)={\rm \pi} /2$, ${\phi }_h(0)=\arctan (-1/\omega )$ and ${\psi }_h(0)=0$ for later convenience, and we may take $\hat {\boldsymbol {e}}_{1}$ to be parallel to the swimming velocity, so that $\boldsymbol {V} = V\hat {\boldsymbol {e}}_{1}$, where $V$ is the swimming speed.
From the analysis of § 4.1, the now constant $\bar {\vartheta }$ obeys (4.10). Given our initial conditions, we therefore have $\cos \bar {\vartheta }=0$, so that $\bar {\vartheta }=\pm {\rm \pi}/2$. Substituting this into (4.14b) gives the time evolution of ${\theta }_h$ as
and we obtain the second constant of integration as $\bar {\varPsi }=\pm {\rm \pi}/2$ from the initial conditions. The choice of sign in $\bar {\vartheta }$ and $\bar {\varPsi }$ is determined by the initial values of the derivative in (4.6), giving $\bar {\vartheta } = -{\rm \pi} /2$, yielding
The dynamics of ${\psi }_h$ are then obtained from (4.14c) and (4.14d) as
The third and final constant of integration $\bar {\varphi }$ is obtained using (4.17), which here reads
The initial condition for ${\phi }_h$ then gives $\bar {\varphi }=- {\rm \pi}/2$, yielding the dynamics of ${\phi }_h$ as
With reference to the Euler angle parameterisation given explicitly in (A3), the translational dynamics $\mathrm {d}\boldsymbol {X}/\mathrm {d}t=V\hat {\boldsymbol {e}}_{1}$ can be written in the lab frame as
By substituting (B2), (B3) and (B5) into these translational evolution equations gives, after integration in time,
which describe a helical trajectory. The radius $\rho$ and pitch $b$ of the helix are then given by
Appendix C. Initial conditions for the slow-time functions
In this appendix we provide the appropriate initial conditions for the slow-time functions, i.e. $\bar {\vartheta }(0)$, $\bar {\varPsi }(0)$ and $\bar {\varphi }(0)$, given generic initial conditions for the original variables, i.e. $\theta (0)$, $\psi (0)$ and $\phi (0)$.
We first use (4.10) to determine that
We then combine (4.14a,c,d) to deduce that
Since $\cos$ and $\tan$ are not monotonic, (C1) and (C2) do not uniquely prescribe $\bar {\vartheta }(0)$ and $\bar {\varPsi }(0)$. Therefore, we must prescribe additional consistency conditions using results from § 4. To this end, we generate the following two additional consistency conditions from (4.14a,c,d):
Then, we perform the following procedure to determine $\bar {\vartheta }(0)$ and $\bar {\varPsi }(0)$. We take the principle values of $\arccos$ and $\arctan$ to calculate $\bar {\vartheta }(0) \in [0, {\rm \pi}]$ from (C1), and to generate a potential value of $\bar {\varPsi }(0) \in [-{\rm \pi} /2, {\rm \pi}/2]$ from (C2). We then check whether this potential value of $\bar {\varPsi }(0)$ is consistent with (C3). If it is not, we are on the wrong branch of $\bar {\varPsi }$, and so we instead use $\bar {\varPsi }(0) \mapsto \bar {\varPsi }(0) + {\rm \pi}$, which means $\bar {\varPsi }(0) \in [{\rm \pi} /2, 3{\rm \pi} /2]$.
Finally, we must determine