Hostname: page-component-586b7cd67f-dsjbd Total loading time: 0 Render date: 2024-12-08T19:27:48.193Z Has data issue: false hasContentIssue false

Autonomous propulsion of nanorods trapped in an acoustic field – CORRIGENDUM

Published online by Cambridge University Press:  07 February 2022

Jesse F. Collis
Affiliation:
ARC Centre of Excellence in Exciton Science, School of Mathematics and Statistics, The University of Melbourne, VIC3010, Australia
Debadi Chakraborty
Affiliation:
ARC Centre of Excellence in Exciton Science, School of Mathematics and Statistics, The University of Melbourne, VIC3010, Australia
John E. Sader*
Affiliation:
ARC Centre of Excellence in Exciton Science, School of Mathematics and Statistics, The University of Melbourne, VIC3010, Australia
*
Email address for correspondence: jsader@unimelb.edu.au

Abstract

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

In Collis, Chakraborty & Sader (Reference Collis, Chakraborty and Sader2017), the effect of Lagrangian motion of the particle on its propulsion velocity, $U_{prop}$, was ignored and is now included. Qualitative features of the propulsion remain unchanged, while some quantitative changes arise. All symbols are as defined in Collis et al. (Reference Collis, Chakraborty and Sader2017), unless otherwise specified.

(i) We rederive the key results of § 2, which replaces the derivation in Appendix A. Consider a general rigid-body particle motion where the time-dependent position vector of a material point on the particle surface (in an inertial frame), $\boldsymbol {r}_p(t) \in S_p (t)$, is

(1.1)\begin{equation} \boldsymbol{r}_p(t) = \boldsymbol{X}(t) + \boldsymbol{\mathsf{A}}(t)\boldsymbol{\cdot}\boldsymbol{r}_0, \end{equation}

where $\boldsymbol {X}(t)$ is the displacement of the particle's geometric centre, $\boldsymbol{\mathsf{A}}(t)$ is the rotation tensor and $\boldsymbol {r}_0 \in S_0$ is the initial position of the same material point. The velocity of this material point is therefore

(1.2)\begin{equation} \boldsymbol{U}_p(t) = \boldsymbol{U}(t) + \boldsymbol{{\varOmega}}(t)\times(\boldsymbol{\mathsf{A}}(t)\boldsymbol{\cdot} \boldsymbol{r}_0), \end{equation}

where $\boldsymbol {U}$ and $\boldsymbol {{\varOmega }}$ are the particle's linear and angular velocities, respectively.

Axisymmetry of the particle ensures rotation can only occur about the $y$-axis and small-amplitude motion in the far field constrains the angular displacement to be

(1.3)\begin{equation} \theta(t) = \epsilon \, {\rm Re}[- \text{i} \bar{{\varOmega}}_p \text{e}^{-\text{i}t}] + o(\epsilon), \end{equation}

and hence

(1.4)\begin{equation} \boldsymbol{{\varOmega}}(t) = {\rm Re}[-\bar{{\varOmega}}_p \text{e}^{-\text{i}t}] \, \boldsymbol{j} + o(\epsilon), \end{equation}

where $\bar {{\varOmega }}_p$ is defined in (2.6). Note that the $O(\epsilon )$ term in (1.4), i.e. for steady and twice-frequency rotations, is zero. Similarly,

(1.5a,b)\begin{equation} \boldsymbol{X}(t) = \epsilon \, {\rm Re}[\text{i} \bar{W}_p \text{e}^{-\text{i}t}]\boldsymbol{k} + o(\epsilon) , \quad \boldsymbol{U}(t) = {\rm Re}[\bar{W}_p \text{e}^{-\text{i}t}]\boldsymbol{k} + \epsilon\boldsymbol{U}^{(1)}_p + o(\epsilon), \end{equation}

with $\bar {W}_p$ specified in (2.6). Substituting (1.3) and (1.5a) into (1.1), and Taylor expanding in $\epsilon$, gives

(1.6)\begin{equation} \boldsymbol{r}_p(t) = \boldsymbol{r}_0 + \epsilon \, {\rm Re}[\bar{\boldsymbol{r}}_1\text{e}^{-\text{i}t}] + o(\epsilon), \end{equation}

where $\bar {\boldsymbol {r}}_1 = \text {i}(\bar {W}_p \boldsymbol {k} + \bar {{\varOmega }}_p \boldsymbol {r}_0\times \boldsymbol {j})$. Taylor expanding the fluid velocity field at the particle surface yields

(1.7)\begin{equation} \boldsymbol{u}|_{\boldsymbol{r}_p} = \boldsymbol{u}|_{\boldsymbol{r}_0} + \epsilon \, {\rm Re}[\bar{\boldsymbol{r}}_1\text{e}^{-\text{i}t}]\boldsymbol{\cdot}\boldsymbol{\nabla}\boldsymbol{u}|_{\boldsymbol{r}_0} + o(\epsilon). \end{equation}

Substituting (1.7) into the asymptotic expansion for the fluid velocity field, (2.4a), gives the boundary conditions for the leading-order flow, (2.6), and the first-order flow,

(1.8)\begin{equation} \bar{\boldsymbol{u}}^{(1)}|_{\boldsymbol{r}_0} = U_{prop}\boldsymbol{i} - \tfrac{1}{4} (\bar{\boldsymbol{r}}_1\boldsymbol{\cdot}\boldsymbol{\nabla} \bar{\boldsymbol{u}}^{(0)*}|_{\boldsymbol{r}_0} + \bar{\boldsymbol{r}}_{1}^{*} \boldsymbol{\cdot}\boldsymbol{\nabla} \bar{\boldsymbol{u}}^{(0)}|_{\boldsymbol{r}_0}), \end{equation}

where $\bar {\boldsymbol {u}}^{(1)}$ and $U_{prop}$ are the steady components of $\boldsymbol {u}^{(1)}$ and $\boldsymbol {i} \boldsymbol {\cdot } \boldsymbol {U}^{(1)}_p$, respectively.

To determine the propulsion velocity, $U_{prop}$, we first calculate the force acting on the moving particle, i.e. $\boldsymbol {F} \equiv \int _{S_p} \boldsymbol {n}\boldsymbol {\cdot } \boldsymbol {\sigma } \,{\rm d}S$. Expanding the normal vector, $\boldsymbol {n}$, gives

(1.9)\begin{equation} \boldsymbol{n} = \boldsymbol{n}_0 + \epsilon \, {\rm Re}[\bar{\boldsymbol{n}}_1 \text{e}^{-\text{i}t}] + o(\epsilon), \end{equation}

where $\bar {\boldsymbol {n}}_1 = \text {i} \bar {{\varOmega }}_p \boldsymbol {n}_0 \times \boldsymbol {j}$ and $\boldsymbol {n}_0$ is the normal vector to the surface, $S_0$, defined in the text immediately following (1.1). Substituting the asymptotic expansion for the fluid variables, (2.4a,b), into the resulting expression for the force, and then Taylor expanding the surface integral, gives $\boldsymbol {F} = \boldsymbol {F}^{(0)} + \epsilon \boldsymbol {F}^{(1)} + o(\epsilon )$. The Fourier component of $\boldsymbol {F}^{(0)}$ is non-zero in the $z$-direction only, which is given by the left-hand side of (2.7a). Moreover, the steady component of $\boldsymbol {F}^{(1)}$ is

(1.10)\begin{align} \bar{\boldsymbol{F}}^{(1)} = \int_{S_0} \boldsymbol{n}_0\boldsymbol{\cdot}\bar{\boldsymbol{\sigma}}^{(1)} + \tfrac{1}{4}[\bar{\boldsymbol{n}}_1\boldsymbol{\cdot}\bar{\boldsymbol{\sigma}}^{(0)*} +\bar{\boldsymbol{n}}_1^*\boldsymbol{\cdot}\bar{\boldsymbol{\sigma}}^{(0)} +\boldsymbol{n}_0\boldsymbol{\cdot}(\bar{\boldsymbol{r}}_1\boldsymbol{\cdot} \boldsymbol{\nabla}\bar{\boldsymbol{\sigma}}^{(0)*} +\bar{\boldsymbol{r}}_1^*\boldsymbol{\cdot}\boldsymbol{\nabla} \bar{\boldsymbol{\sigma}}^{(0)})]\,{\rm d}S, \end{align}

which is non-zero in the $x$-direction only.

To evaluate $\boldsymbol {i}\boldsymbol {\cdot }\int _{S_0} \boldsymbol {n}_0 \boldsymbol {\cdot }\bar {\boldsymbol {\sigma }}^{(1)}\,{\rm d}S$, we use the generalised Lorentz reciprocal theorem on the fluid domain, $V_0$, which is specified as the exterior of the surface defined by the particle's initial position, i.e. $S_0$; not the true (time-dependent) fluid domain. The two flows used in the reciprocal theorem are (a) the first-order steady Stokes flow with body force as in (2.9b), but with boundary condition (1.8), and (b) an auxiliary (body-force-free) Stokes flow, $(\boldsymbol {u}', \boldsymbol {\sigma }')$, specified by unitary translation of the particle in the $x$-direction. This gives

(1.11)\begin{align} \int_{S_0} \boldsymbol{n}_0\boldsymbol{\cdot}(\boldsymbol{i} \boldsymbol{\cdot}\bar{\boldsymbol{\sigma}}^{(1)} - \bar{\boldsymbol{u}}^{(1)} \boldsymbol{\cdot} \boldsymbol{\sigma}')\,{\rm d}S ={-}\frac{\beta}{4}\int_{V_0} \boldsymbol{u}'\boldsymbol{\cdot}(\bar{\boldsymbol{u}}^{(0)} \boldsymbol{\cdot}\boldsymbol{\nabla}\bar{\boldsymbol{u}}^{(0)*} + \bar{\boldsymbol{u}}^{(0)*}\boldsymbol{\cdot} \boldsymbol{\nabla}\bar{\boldsymbol{u}}^{(0)})\,{\rm d}V. \end{align}

Conservation of the particle's linear momentum (to first order in $\epsilon$) requires $\bar {\boldsymbol {F}}^{(1)} = 0$. Substituting (1.8) into (1.11) and its result into (1.10) gives the corrected expression for (2.10):

(1.12)\begin{align} U_{{prop}} &= \frac{\beta}{4F_p} \int_{V_0} \boldsymbol{u}'\boldsymbol{\cdot}(\boldsymbol{\bar{u}}^{(0)} \boldsymbol{\cdot}\boldsymbol{\nabla}\bar{\boldsymbol{u}}^{(0)*} +\bar{\boldsymbol{u}}^{(0)*} \boldsymbol{\cdot}\boldsymbol{\nabla}\bar{\boldsymbol{u}}^{(0)})\,\mathrm{d}V \nonumber\\ & \quad + \frac{1}{4F_p}\int_{S_0} (\boldsymbol{n}_{0}\boldsymbol{\cdot}\boldsymbol{\sigma}') \boldsymbol{\cdot}(\boldsymbol{\bar{r}}_{1}\boldsymbol{\cdot}\boldsymbol{\nabla} \bar{\boldsymbol{u}}^{(0)*} + \bar{\boldsymbol{r}}^{*}_{1}\boldsymbol{\cdot}\boldsymbol{\nabla} \bar{\boldsymbol{u}}^{(0)}) \nonumber\\ &\quad - \boldsymbol{i}\boldsymbol{\cdot}[ \boldsymbol{n}_{0} \boldsymbol{\cdot}(\boldsymbol{\bar{r}}_{1}\boldsymbol{\cdot}\boldsymbol{\nabla} \bar{\boldsymbol{\sigma}}^{(0)*} + \bar{\boldsymbol{r}}^{*}_{1}\boldsymbol{\cdot} \boldsymbol{\nabla}\bar{\boldsymbol{\sigma}}^{(0)}) + \bar{\boldsymbol{n}}_{1} \boldsymbol{\cdot}\bar{\boldsymbol{\sigma}}^{(0)*} + \bar{\boldsymbol{n}}^{*}_{1} \boldsymbol{\cdot} \bar{\boldsymbol{\sigma}}^{(0)} ] \, \mathrm{d}S; \end{align}

the surface integral in (1.12) must also be added to (2.11).

(ii) Corrected asymptotic forms for the propulsion velocity of the dumbbell, (3.6) and (3.7), are, respectively,

(1.13) \begin{equation} U_{prop} = \frac{(\kappa - 1)(\gamma_1 - 1)^2}{A}\left\{ \begin{array}{@{}ll} \displaystyle -\dfrac{1+\kappa^4}{243} \beta_1^{3}, & \beta_1 \ll 1, \\ \displaystyle \dfrac{45\sqrt{2}}{(1+2\gamma_1)^3 \kappa(1+\kappa)} \beta_1^{{-}1/2} , & \beta_1 \gg 1, \end{array} \right. \end{equation}

and

(1.14) \begin{equation} U_{prop} = \frac{(\gamma_2 - \gamma_1)}{A} \left\{ \begin{array}{@{}ll} \displaystyle -\dfrac{\gamma_1+\gamma_2-2}{486} \beta_1^{3}, & \beta_1 \ll 1,\\ \displaystyle \dfrac{3(4\gamma_1\gamma_2 -\gamma_1 - \gamma_2 - 2)}{(1 + 2\gamma_1)^2 (1 + 2\gamma_2)^2} , & \beta_1 \gg 1. \end{array} \right. \end{equation}

(iii) Corrected figures 2 and 4 are similar to Collis et al. (Reference Collis, Chakraborty and Sader2017), with quantitative changes to the zero propulsion point. Corrected table 1 shows closer agreement with experiments, which indicate that the concave and light end lead the motion; similar changes to figure 5 (not shown).

Figure 2. Corrected: axis scales are identical to Collis et al. (Reference Collis, Chakraborty and Sader2017).

Figure 4. Corrected: axis scales are identical to Collis et al. (Reference Collis, Chakraborty and Sader2017).

Table 1. Corrected.

(iv) Changes to the caption of figure 3: (c) Non-zero contribution from the body force is cancelled by that of the surface integrals in (2.10), giving $U_{{prop}} = 0$. (f) Combining these two flows produces a jet that propels the sphere in the negative $x$-direction, with $U_{prop} = -0.039$.

Declaration of interests

The authors report no conflict of interest.

References

REFERENCES

Collis, J.F., Chakraborty, D. & Sader, J.E. 2017 Autonomous propulsion of nanorods trapped in an acoustic field. J. Fluid Mech. 825, 2948.CrossRefGoogle Scholar
Figure 0

Figure 2. Corrected: axis scales are identical to Collis et al. (2017).

Figure 1

Figure 4. Corrected: axis scales are identical to Collis et al. (2017).

Figure 2

Table 1. Corrected.