Hostname: page-component-586b7cd67f-vdxz6 Total loading time: 0 Render date: 2024-11-28T09:01:05.337Z Has data issue: false hasContentIssue false

Vortex dynamics on a Möbius strip

Published online by Cambridge University Press:  22 July 2021

Jacques Vanneste*
Affiliation:
School of Mathematics and Maxwell Institute for Mathematical Sciences, University of Edinburgh, King's Buildings, EdinburghEH9 3FD, UK
*
Email address for correspondence: [email protected]

Abstract

We consider the dynamics of a two-dimensional incompressible perfect fluid on a Möbius strip embedded in $\mathbb {R}^{3}$. The vorticity–stream function formulation of the Euler equations is derived from an exterior-calculus form of the momentum equation. The non-orientability of the Möbius strip and the distinction between forms and pseudo-forms this introduces lead to unusual properties: a boundary condition is provided by the conservation of circulation along the single boundary of the strip, and there is no integral conservation for the vorticity density or for any odd function thereof. A finite-difference numerical implementation is used to illustrate the Möbius-strip realisation of familiar phenomena: translation of vortices along boundaries, shear instability and decaying turbulence.

Type
JFM Papers
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2021. Published by Cambridge University Press

1. Introduction

In this paper, we examine the dynamics of a two-dimensional incompressible inviscid fluid confined to a Möbius strip. The main motivation is a curiosity about the impact of the non-orientability of the Möbius strip on the phenomenology of vortex dynamics. A secondary, more frivolous motivation is the wish to produce interesting fluid-dynamical pictures. The system is physically realisable: soap films, for example, can be used as proxies for two-dimensional fluids (e.g. Couder & Basdevant Reference Couder and Basdevant1986; Couder, Chomaz & Rabaud Reference Couder, Chomaz and Rabaud1989); with a suitable wire frame they take the topology of a Möbius strip (e.g. Courant Reference Courant1940; Goldstein et al. Reference Goldstein, Moffatt, Pesci and Ricca2010). The experiments in this paper are, however, exclusively numerical.

The Möbius strip is a topological object, but fluid dynamics requires a geometry, which we need to choose. The simplest geometry is that of a flat strip, imposing antisymmetry conditions across a ‘seam’ to achieve the necessary twist. However, the flat Möbius strip cannot be embedded in $\mathbb {R}^{3}$, a necessary condition for physical realisability. We therefore have chosen another widely used model, the ruled surface constructed by rotating the midpoint of a segment along a horizontal circle while rotating the segment in the vertical at half the rate. This surface, which we denote by $\mathcal {M}$, is parameterised as

(1.1a)$$\begin{gather} x=(1+\zeta \cos \theta) \cos(2 \theta), \end{gather}$$
(1.1b)$$\begin{gather}y=(1+\zeta \cos \theta) \sin(2 \theta), \end{gather}$$
(1.1c)$$\begin{gather}z=\zeta \sin \theta, \end{gather}$$

with $-a \le \zeta \le a$ and $0 \le \theta \le {\rm \pi}$ (or, rather, $\theta$ interpreted as taken modulo ${\rm \pi}$). The property $\boldsymbol {x}(\zeta ,\theta +{\rm \pi} )=\boldsymbol {x}(-\zeta ,\theta )$ reflects the twist of the Möbius-strip topology. The parameterisation and shape of $\mathcal {M}$ are illustrated in figure 1. The mean curvature of $\mathcal {M}$ does not vanish, so it is not a minimal surface as formed by a soap film. Its Gaussian curvature does not vanish either, so it is not developable and cannot be made of an inextensible material. Closed-form expressions are available for Möbius strips in $\mathbb {R}^{3}$ that are either minimal (Odehnal Reference Odehnal2016) or developable (see Schwarz Reference Schwarz1990) and could replace (1.1) to gain realism. The equilibrium shape taken by an inextensible material twisted as a Möbius strip can also be obtained, but this requires solving an energy-minimisation problem numerically (Starostin & van der Heijden Reference Starostin and van der Heijden2007).

Figure 1. Möbius strip in $\mathbb {R}^{3}$ and its parameter space for $a=0.25$. The Möbius strip is coloured according to the value of $|g|^{1/2}$, with $|g|$ the metric determinant in (2.3), which ranges from approximately $1.5$ (dark blue) to $2.5$ (light yellow) and whose variations are associated with curvature effects.

The non-orientability of $\mathcal {M}$ is manifested in the fact that the basis $(\partial _{\zeta } \boldsymbol {x},\partial _\theta \boldsymbol {x})$ of its tangent space has opposite orientations for $\theta =0$ and $\theta ={\rm \pi}$. Equivalently, if a normal vector $\boldsymbol {N}$ is added to $\partial _{\zeta } \boldsymbol {x}$ and $\partial _\theta \boldsymbol {x}$ to make up a right-handed basis of $\mathbb {R}^{3}$, it satisfies $\boldsymbol {N}(\theta ={\rm \pi} )=-\boldsymbol {N}(\theta =0)$. Non-orientability is a global property, not discernible at the level of a local patch of the strip; incompressible fluid dynamics is also global, through the pressure field. Nonetheless, for narrow strips such as the one in figure 1, and with the scale of fluid structures set by the strip's width, we can expect the non-orientability to have a limited impact on the dynamics, which differs from that in a flat cylinder (channel) primarily because of the strip's curvature. The main effect of non-orientability arises from the fact that the sense of rotation of vortices is reversed when they are transported once along the entire length of the strip, that is, when $\theta$ increases by ${\rm \pi}$. We will illustrate this by considering the dynamics of a single vortex propagating along the strip's boundary. A mathematical consequence is that vorticity cannot be meaningfully integrated over the Möbius strip; as a result, the familiar connection between circulation and the vorticity integral is lost, and the conservation of integrals of functions of the vorticity (Casimirs) is restricted to even functions.

The paper is structured as follows. In § 2, we derive the vorticity–stream function formulation of the Euler equations for an incompressible perfect fluid on $\mathcal {M}$ in terms of the coordinates $(\zeta ,\theta )$. We start with a convenient formulation of the Euler equations that highlights the dual role of the vector and 1-form fields associated with velocity, and treats vorticity as a 2-form. We carefully distinguish between forms (such as the vorticity 2-form) and pseudo-forms (such as the area 2-form, the stream function and a pseudo-scalar version of the vorticity which we term ‘vorticity density’). The distinction is crucial for a non-orientable surface such as $\mathcal {M}$. We derive the boundary conditions for the vorticity–stream function formulation; these include a condition involving the circulation along the (single) boundary of $\mathcal {M}$. In § 3 we describe a numerical implementation of the equations which we employ to examine the propagation of a vortex along the strip's boundary, shear instability and decaying turbulence. Section 4 summarises the main features of vortex dynamics on a non-orientable surface. Two appendices provide mathematical details.

2. Euler equations on the Möbius strip

A coordinate expression for the Euler equations governing the dynamics of an incompressible perfect fluid on an arbitrary surface was derived by Zermelo (Reference Zermelo1902) (see Zermelo (Reference Zermelo2013) for a commented English translation). Modern formulations with a focus on point vortices and closed surfaces can be found in, for example, Kimura (Reference Kimura1999), Boatto & Koiller (Reference Boatto and Koiller2015), Dritschel & Boatto (Reference Dritschel and Boatto2015), Ragazzo & de Barros Viglioni (Reference Ragazzo and de Barros Viglioni2017) or Gustafsson (Reference Gustafsson2019). While Zermelo's expressions can be readily particularised to the Möbius strip, it is instructive to start from a general, coordinate-free version of the Euler equations. A convenient form is

(2.1a)$$\begin{gather} (\partial_t + \mathcal{L}_{\boldsymbol{u}}) \nu ={-} {\rm d} \left( p - \tfrac{1}{2} | \boldsymbol{u}|^{2} \right), \end{gather}$$
(2.1b)$$\begin{gather}{\mathrm{div}} \boldsymbol{u} = 0. \end{gather}$$

Here $\boldsymbol {u}$ is the velocity (vector) field, $\nu$ is the momentum 1-form, obtained from $\boldsymbol {u}$ and the metric $g$ as $\nu = g(\boldsymbol {u},\cdot )$, $\mathcal {L}_{\boldsymbol {u}}$ is the Lie derivative and $\textrm {d}$ the exterior derivative (see e.g. Schutz Reference Schutz1980; Arnold & Khesin Reference Arnold and Khesin1998; Besse & Frisch Reference Besse and Frisch2017; Gilbert & Vanneste Reference Gilbert and Vanneste2018).

The metric on the Möbius strip $\mathcal {M}$ is found by pulling back the Euclidean metric ${\textrm {d}x} \otimes {\textrm {d}x} + {\textrm {d}y} \otimes {\textrm {d}y} + \textrm {d} z \otimes \textrm {d} z$ of $\mathbb {R}^{3}$, that is, by using (1.1) to express ${\textrm {d}x}$, ${\textrm {d}y}$ and $\textrm {d} z$ in terms of $\textrm {d} \zeta$ and $\textrm {d} \theta$ to find

(2.2)\begin{equation} g = \textrm{d} \zeta \otimes \textrm{d} \zeta + \left(4 + 8 \cos \theta \zeta + (3+2 \cos(2 \theta)) \zeta^{2}\right) \textrm{d} \theta\otimes \textrm{d} \theta. \end{equation}

This is a diagonal tensor, indicating that the coordinates $(\zeta ,\theta )$ are orthogonal (i.e. the basis vectors $\partial _\zeta$ and $\partial _\theta$ are orthogonal), with determinant

(2.3)\begin{equation} |g| = 4 + 8 \cos \theta \zeta + \left(3+2 \cos(2 \theta)\right) \zeta^{2}. \end{equation}

It follows from (2.2) that the vectors $\partial _\zeta$ and $\partial _\theta$ have norms $|\partial _\zeta | =1$ and $|\partial _\theta |=|g|^{1/2}$. Writing the velocity field in components as

(2.4)\begin{equation} \boldsymbol{u} = u^{\zeta} \partial_\zeta + u^{\theta} \partial_\theta, \end{equation}

the corresponding momentum 1-form is found from $\nu =g(\boldsymbol {u},\cdot )$ as

(2.5)\begin{equation} \nu = \nu_\zeta \,\textrm{d} \zeta + \nu_\theta \,\textrm{d} \theta = u^{\zeta} \,\textrm{d} \zeta + |g| u^{\theta} \,\textrm{d} \theta. \end{equation}

Introducing (2.5) into (2.1a) with $\mathcal {L}_{\boldsymbol {u}} = u^{\zeta } \partial _\zeta + u^{\theta } \partial _\theta$ and using the commutation $\mathcal {L}_{\boldsymbol {u}} \,\textrm {d} = \textrm {d} \mathcal {L}_{\boldsymbol {u}}$ readily gives the momentum equation in terms of $u^{\zeta }$ and $u^{\theta }$. This is worked out in Appendix A.

The divergence in (2.1b) is defined by the relation ${\mathrm {div}} \boldsymbol {u} \,\mu = \textrm {d} (\boldsymbol {u} \lrcorner \mu ) = \mathcal {L}_{\boldsymbol {u}} \mu$, where $\mu$ is the volume form associated with the metric $g$ (here area form since $\mathcal {M}$ is two-dimensional) and $\lrcorner$ denotes the interior product or contraction. The area form is given by

(2.6)\begin{equation} \mu ={\pm} |g|^{1/2} \,\textrm{d} \zeta \wedge \textrm{d} \theta. \end{equation}

This is a pseudo-2-form, which changes sign when the orientation of the basis $(\partial_\zeta,\partial_\theta)$ changes, as indicated by the notation $\pm$. Since $\mathcal {M}$ is non-orientable, a change in the orientation of the basis is unavoidable; it occurs across $\theta = 0 \mod {\rm \pi}$ in the parameterisation (1.1). It can be verified that (2.6) is the projection of the Euclidean volume ${\textrm {d}x} \wedge {\textrm {d}y} \wedge \textrm {d} z$ on the normal $\boldsymbol {N}=(N^{1},N^{2},N^{3})$ to the Möbius strip, that is, $\mu = {\boldsymbol {N}} \lrcorner \, {\textrm {d}x} \wedge {\textrm {d}y} \wedge \textrm {d} z = N^{1} {\textrm {d}y} \wedge \textrm {d} z + N^{2} \,\textrm {d} z \wedge \mathrm {d} x + N^{3} \,\mathrm {d} x \wedge {\textrm {d}y}$. With (2.6) we compute

(2.7)\begin{equation} \boldsymbol{u} \lrcorner \mu ={\pm} |g|^{1/2} (u^{\zeta} \,\textrm{d} \theta - u^{\theta} \,\textrm{d} \zeta), \end{equation}

and apply ${\rm d}$ to obtain the incompressibility condition in the explicit form

(2.8)\begin{equation} {\mathrm{div}} \boldsymbol{u} = {|g|^{{-}1/2}} \left(\partial_\zeta ( |g|^{1/2} u^{\zeta}) + \partial_\theta (|g|^{1/2} u^{\theta}) \right) = 0. \end{equation}

The momentum and incompressibility equations (A3) and (2.8) can be conveniently replaced by the vorticity equation. This is best derived directly from (2.1a): applying ${\rm d}$ and using that ${\rm d}$ and $\mathcal {L}_{\boldsymbol {u}}$ commute and that ${\rm d}^{2}=0$ gives

(2.9)\begin{equation} (\partial_t + \mathcal{L}_{\boldsymbol{u}}) \,\textrm{d} \nu = 0, \end{equation}

showing that the vorticity 2-form $\textrm {d} \nu$ is transported by the flow. A vorticity density $\omega$ can be defined by

(2.10)\begin{equation} \omega \mu = \textrm{d} \nu. \end{equation}

This is a pseudo-scalar, which, like $\mu$, changes sign with the orientation of the basis on $\mathcal {M}$. Applying ${\rm d}$ to (2.4) and using (2.6) gives the explicit expression

(2.11)\begin{equation} \omega ={\pm} |g|^{{-}1/2} \left(\partial_\zeta (|g| u^{\theta}) - \partial_\theta u^{\zeta} \right). \end{equation}

Incompressibility ensures that $(\partial _t + \mathcal {L}_{\boldsymbol {u}}) \mu = 0$, hence $\omega$ is also transported by the flow:

(2.12)\begin{equation} (\partial_t + \mathcal{L}_{\boldsymbol{u}}) \omega =0. \end{equation}

The velocity $\boldsymbol {u}$ can be related to the vorticity by means of a stream function $\psi$, with

(2.13a,b)\begin{equation} u^{\zeta} ={-} |g|^{{-}1/2} \partial_\theta \psi \quad \text{and}\quad u^{\theta} = |g|^{{-}1/2} \partial_\zeta \psi, \end{equation}

so as to satisfy the incompressibility condition (2.8). Equivalently, this can be written as

(2.14)\begin{equation} \textrm{d} \psi ={-} \boldsymbol{u} \lrcorner \mu . \end{equation}

With (2.13a,b), (2.12) becomes

(2.15)\begin{equation} \partial_t \omega + |g|^{{-}1/2} \partial(\psi,\omega)=0, \end{equation}

where $\partial (\psi ,\omega ) = \partial _\zeta \psi \partial _\theta \omega - \partial _\theta \psi \partial _\zeta \omega$ is the Jacobian operator. Substituting (2.13a,b) into (2.11) then gives the vorticity–stream function relation

(2.16)\begin{equation} \omega = {\rm \Delta} \psi = {|g|^{{-}1/2}} \left(\partial_\zeta \left(|g|^{1/2} \partial_\zeta \psi\right) + \partial_\theta \left(|g|^{{-}1/2} \partial_\theta \psi \right) \right), \end{equation}

where we have identified the Laplacian ${\rm \Delta}$. A completely geometric derivation of (2.13a,b) and (2.16), which avoids coordinates but requires the introduction of the Hodge * operator, is given in Appendix B.

Equations (2.15) and (2.16), together with the form (2.3) of the metric determinant $|g|$ constitute the vorticity–stream function formulation of the Euler equations on the Möbius strip. In the familiar way, the vorticity density $\omega$ is the dynamical variable from which the stream function $\psi$ is derived by inverting the Poisson equation (2.16). This requires boundary conditions. Because $\omega$ and $\psi$ are pseudo-scalars, they satisfy

(2.17a,b)\begin{equation} \omega(\zeta,\theta+{\rm \pi},t) ={-}\omega(-\zeta,\theta,t) \quad \text{and}\quad \psi(\zeta,\theta+{\rm \pi},t) ={-}\psi(-\zeta,\theta,t). \end{equation}

In particular,

(2.18)\begin{equation} \psi(\zeta,{\rm \pi},t)={-}\psi(-\zeta,0,t), \end{equation}

which provides a first boundary condition for (2.16). The no-normal-flow condition

(2.19)\begin{equation} u^{\zeta}(\zeta={\pm} a,\theta,t)=0 \end{equation}

and (2.13a,b) imply that $\psi$ depends only on $t$ along the boundary of the strip. Because of the pseudo-scalar nature of $\psi$, this constancy translates into

(2.20)\begin{equation} \psi (\zeta ={\pm} a, \theta,t)={\pm} C(t). \end{equation}

This provides the second boundary condition for the Poisson equation (2.16). To determine $C(t)$, we need to revert to the momentum formulation (2.1). It implies conservation of the circulation along the (single) boundary, given by

(2.21)\begin{align} \varGamma &= \int_{\partial \mathcal{M}} \nu = \int_0^{\rm \pi} \left(\nu_\theta(a,\theta,t)+\nu_\theta({-}a,\theta,t)\right) \textrm{d} \theta \nonumber\\ &= \int_0^{\rm \pi} \left({|g(a,\theta)|^{1/2}} \partial_\zeta \psi (a,\theta,t)+{|g({-}a,\theta)|^{1/2}} \partial_\zeta \psi ({-}a,\theta,t)\right) \textrm{d} \theta. \end{align}

The equation

(2.22)\begin{equation} \frac{\mathrm{d} \varGamma}{\mathrm{d} t}=0 \end{equation}

completes the formulation. The description of our numerical implementation in § 3 makes it clear how (2.22) determines $C(t)$ for use in (2.20).

The existence of infinitely many conserved integrals is a key feature of two-dimensional inviscid incompressible fluids. Here the non-orientability of the Möbius strip leads to important differences compared with the familiar, orientable case. This is because, while pseudo-2-forms can be integrated over non-orientable surfaces, proper 2-forms cannot – the result of such an integration would not be coordinate independent (e.g. Frankel Reference Frankel2004). Thus, while the area of the Möbius strip $\int _\mathcal {M} \mu$ is well defined, the integral $\int _\mathcal {M} {\rm d} \nu$ of the vorticity 2-form is not (and cannot be related to the circulation $\varGamma$). As a result, conserved Casimirs functions

(2.23)\begin{equation} \mathcal{C}_f = \int_{\mathcal{M}} f(\omega) \mu = \int_{{-}a}^{a} \int_0^{\rm \pi} f(\omega) |g|^{1/2} \,\textrm{d} \zeta \,\textrm{d} \theta, \end{equation}

can only be defined for arbitrary even functions $f$. The evenness of $f$ is required to ensure that the change of sign of $\omega$ which accompanies a change of basis orientation has no effect on $f(\omega )$, so that $f(\omega ) \mu$ is a pseudo-2-form. The energy

(2.24)\begin{equation} \mathcal{E} ={-} \tfrac{1}{2} \int_M \psi \omega \mu ={-} \tfrac{1}{2} \int_{{-}a}^{a} \int_0^{\rm \pi} \psi \omega |g|^{1/2} \,\textrm{d} \zeta\,\textrm{d} \theta \end{equation}

is also well defined and conserved.

3. Numerical simulations

3.1. Numerical implementation

We solve the vorticity–stream function equations (2.15) and (2.16) numerically. To control the generation of fine scales that inevitably accompanies non-trivial flows, we modify (2.15), adding a small-scale dissipation to obtain

(3.1)\begin{equation} \partial_t \omega + |g|^{{-}1/2} \partial(\psi,\omega) = \varepsilon {\rm \Delta} \omega, \end{equation}

with $\varepsilon >0$ small. This requires an additional boundary condition, taken as $\omega =0$ on the strip boundary $\partial \mathcal {M}$, corresponding to a no-stress condition. With dissipation, the circulation $\varGamma$ along $\partial \mathcal {M}$ changes. We derive the form of $\mathrm {d} \varGamma /\mathrm {d} t$ in Appendix B where it is given as (B6). Note that the dissipation model chosen is not standard Newtonian friction: the complete Navier–Stokes vorticity equation involves geometric terms omitted from the right-hand side of (3.1) (see Appendix B, Gilbert, Riedinger & Thuburn (Reference Gilbert, Riedinger and Thuburn2014) and Gilbert & Vanneste (Reference Gilbert and Vanneste2021)). In the simulations of § 3.2, $\varepsilon$ is taken small enough that the results are representative of the limit $\varepsilon \to 0$ except for scales close to the grid size.

To solve (2.16) and (3.1) numerically, we discretise $\omega$ and $\psi$ on a regular grid in $(\zeta ,\theta )$ coordinates. The Laplacian ${\rm \Delta}$ is discretised using a straightforward five-point stencil, accounting for the condition (2.18) at the seam $\theta =0$ and assuming homogeneous conditions for $\zeta =\pm a$ corresponding to the strip's boundary $\partial \mathcal {M}$. The resulting matrix is sparse and has an interesting structure displayed in figure 2, with the standard band structure modified by antidiagonal blocks in the top-right and bottom-left corners. These result from the twist of condition (2.18). The boundary condition (2.20) is implemented as follows. We first compute a numerical approximation to the (time-independent) harmonic function $\psi _{har}$ solving

(3.2)\begin{equation} {\rm \Delta} \psi_{har} = 0 \quad \textrm{with}\ \psi_{har}(\zeta={\pm} a,\theta)={\pm} 1, \end{equation}

then, at each time step, solve the Poisson equation (2.16) with homogeneous boundary condition

(3.3)\begin{equation} {\rm \Delta} \psi_{hom} = \omega \quad \textrm{with}\ \psi_{hom}(\zeta={\pm} a,\theta,t)=0. \end{equation}

We compute the circulations $\varGamma _{har}$ and $\varGamma _{hom}(t)$ associated with each solution by direct evaluation of (2.21), and we obtain the total circulation $\varGamma (t)$ from the differential equation (B6). The solution of the Poisson equation satisfying the boundary condition (2.20) is then calculated as

(3.4)\begin{equation} \psi = \psi_{hom} + C(t) \psi_{har}, \end{equation}

with $C(t)$ given by

(3.5)\begin{equation} C(t) = \left( \varGamma(t) - \varGamma_{hom}(t) \right)/\varGamma_{har} \end{equation}

to ensure that the circulation associated with $\psi$ is $\varGamma (t)$, and (2.20) holds.

Figure 2. Visualisation of the non-zero entries in the matrix discretising the Laplacian operator on the Möbius strip. The unusual antidiagonal bands in the top-right and bottom-left corners result from the twist of the Möbius strip.

We use a splitting scheme to advance (3.1) in time, with a three-step Adams–Bashforth scheme for the nonlinear advection relying on the Arakawa (Reference Arakawa1966) discretisation of the Jacobian, and a backward Euler scheme for the viscous diffusion. For the simulations reported below, the half-width of the Möbius strip is set to $a=0.25$; the spatial discretisation uses a $200 \times 1000$ grid in $(\zeta ,\theta )$; the time step is $0.2$ and the dissipation parameter is $\epsilon = 4 \times 10^{-7}$. With this small value, energy is approximately conserved over the duration of the simulations presented next.

3.2. Simulations

We start by simulating the dynamics of a single vortex patch initialised near the boundary of the Möbius strip. The vortex is expected to travel along the boundary, forming a dipole with its notional image as is familiar in the planar case. The Möbius strip case is intriguing in that the strip has a single boundary, and the direction of rotation of the fluid can be reversed by transporting the vortex once along the strip. The motion of the vortex is in fact straightforward: the vortex travels continuously along the entire boundary. When the along-strip coordinate $\theta$ of the vortex has increased by ${\rm \pi}$, so that the vortex returns to the same segment $\theta ={const.}$, the sign of $\omega$ is reversed but the vortex initially near the boundary at $\zeta =a$ is then near $\zeta =-a$ and the propagation direction is unchanged. This is illustrated in the snapshots shown in figure 3 and (better) in the accompanying movie 1. These visualise the vorticity 2-form ${\rm d} \nu$ which, unlike $\omega$, is a true geometric object, independent of the choice of coordinate basis. The field ${\rm d} \nu$ can be represented as the (true) vector field $\omega \boldsymbol {N}$, normal to $\mathcal {M}$ and defined intrinsically as dual to ${\rm d}\nu$ via the $\mathbb {R}^{3}$ volume-form ${\textrm {d}x} \wedge {\textrm {d}y} \wedge \textrm {d} z$. The visualisation in figure 3 and movie 1 shows $\omega \boldsymbol {N}$ where $|\omega |$ is maximum, that is, at the centre of the vortex, as well as the scalar field $|\omega |$ on the strip. The propagation of a vortex along the edge of the Möbius strip makes the single-sidedness of the strip manifest and can be viewed as a fluid-dynamical equivalent of Escher's famous print Möbius Strip II showing ants crawling along the strip.

Figure 3. Propagation of a vortex along the boundary of the Möbius strip: the vortex travels anticlockwise along the entire boundary, with a vorticity vector field $\omega \boldsymbol {N}$ (shown by the arrow in the centre of the vortex) reversing direction as the along-strip coordinate $\theta$ increases by ${\rm \pi}$ (see also supplementary movie 1 available at https://doi.org/10.1017/jfm.2021.581).

We next examine a version of shear instability by simulating the evolution of the initial vorticity density

(3.6)\begin{equation} \omega=\zeta (a^{2}-\zeta^{2}). \end{equation}

This is not a stationary solution of the Euler equations because the associated stream function depends on $\theta$ through the metric terms that appear in the Poisson equation (2.16). Its evolution is nonetheless very slow compared with that which occurs when a small $\theta$-dependent perturbation is added – we use a perturbation of the form $\omega ' \propto (a^{2}-\zeta ^{2}) \sin (8 \theta )$. The rapid evolution in the presence of a perturbation can be interpreted as resulting from shear instability. Note the relative complexity of (3.6) compared with the quadratic vorticity profile that is often employed to demonstrate shear instability in planar flows (a quadratic vorticity has a gradient vanishing at a single point, the inflection point of the velocity profile). This is imposed by the topological constraint that 2-forms, here $\textrm {d} \nu = \omega \,{\textrm {d}x} \wedge {\textrm {d}y}$, must vanish somewhere on non-orientable surfaces. Figure 4 shows a series of snapshots of the pseudo-scalar vorticity density $\omega$ at equal time intervals. For distributed vorticity distributions, this field is easier to visualise than the vector field $\omega \boldsymbol {N}$ shown in figure 3, but it has the drawback of an abrupt sign change across the seam $\theta = 0$, visible near the top right in the representation of the Möbius strip in figure 4. We emphasise that this discontinuity is artificial and that the physical fields $\nu$ and $\textrm {d} \nu$ (or equivalently $\omega \boldsymbol {N}$) and the direction of rotation of the fluid are continuous. The snapshots of figure 4 and the accompanying movie 2 show that shear instability on the Möbius strip evolves very much as in planar channel, with the amplification of the initial wavy perturbation of the constant-vorticity lines leading to the formation of rows of counter-rotating vortices.

Figure 4. Shear instability on a Möbius strip: the (pseudo-scalar) vorticity density $\omega$ is shown at six successive time intervals, starting with a perturbation of the field (3.6) (see also movie 2).

Our final simulation, with results displayed in figure 5 and movie 3, illustrates the dynamics of decaying turbulence on the Möbius strip. We build a small-scale, random initial vorticity density $\omega _0(\zeta ,\theta )$, enforcing its pseudo-scalar property by first constructing a random function $f(\zeta ,\theta )$ that is $4a$-periodic in $\zeta$ and $2{\rm \pi}$-periodic in $\theta$ using a standard Fourier series method. The initial vorticity density is then taken as

(3.7)\begin{equation} \omega_0(\zeta,\theta) = f(\zeta,\theta)-f(2a-\zeta,\theta)+f(2a+\zeta,\theta+{\rm \pi})-f(-\zeta,\theta+{\rm \pi}). \end{equation}

The periodicity of $f$ ensures that $\omega _0$ vanishes at the boundary and satisfies condition (2.17a,b) at the seam $\theta ={\rm \pi}$, as required for a pseudo-scalar. The choice of peak wavenumbers $8{\rm \pi} /a$ and $80{\rm \pi}$ for $f$ leads to the reasonably isotropic, small-scale $\omega _0$ shown in the first panel of figure 5. In addition to an initial vorticity density, we also impose a non-zero (weak) initial circulation which leads to an overall drift of the vortices best seen in movie 3. The evolution shown in figure 5 and movie 3 is characteristic of an inverse energy cascade, with vortex mergers that are local and hence unaffected by the topology of the strip. Over time scales much larger than those explored in our simulation, the scale of the vorticity field may become such that the evolution is more strongly influenced by the topology of $\mathcal {M}$.

Figure 5. Decaying turbulence on a Möbius strip: the (pseudo-scalar) vorticity density $\omega$ is shown at six successive equal time intervals (see also movie 3).

4. Summary

This paper examines how the Euler equations governing the dynamics of two-dimensional inviscid fluids can be formulated on a non-orientable surface, using a Möbius strip embedded in $\mathbb {R}^{3}$ as a concrete example. The velocity vector field $\boldsymbol {u}$, the corresponding momentum 1-form $\nu$ and the pressure field $p$ are all geometrically intrinsic objects, which are independent of the choice of coordinate basis regardless of whether the surface is orientable or not. Pseudo-fields, which depend on coordinates through the orientation of the basis, only appear when the Euler equations are simplified to their vorticity–stream function formulation. While the vorticity $\textrm {d} \nu$ is a proper 2-form, the more convenient $\omega$ is a pseudo-scalar because its definition via $\textrm {d} \nu = \omega \mu$ involves the area pseudo-2-form $\mu$. Similarly, the stream function $\psi$ is defined via the relation $\textrm {d} \psi = - \boldsymbol {u} \lrcorner \mu$ and is a pseudo-scalar. With our choice of coordinates $(\zeta ,\theta )$, the difference between pseudo-scalars and scalars is simply that the former satisfy the twist condition (2.18) across $\theta = 0 \mod {\rm \pi}$ instead of continuity. In the vorticity–stream function formulation of the Euler equations, the condition that $\psi$ depends on time only along each connected component of the boundary must be supplemented by additional conditions tracing back to the original momentum formulation. On the Möbius strip, the boundary is simply connected and the single function of time $C(t)$ that needs to be determined is found straightforwardly from conservation of the circulation along the single boundary. The main dynamical effect of non-orientability is that the sign of $\omega$ can be changed by transport: a vortex patch rotating clockwise transported once along the strip returns rotating anticlockwise, as demonstrated in the numerical experiment of figure 3. This is why odd functions of the vorticity density are not conserved in an integral sense. Loosely speaking, the number of integral conservation laws on the Möbius strip is half of what it is on an orientable surface. It would be interesting to examine the consequences of this for the long-time dynamics of turbulence, for example as described by statistical mechanics. It would also be interesting to study purely non-dissipative dynamics on a Möbius strip by considering point vortices or contour dynamics.

Supplementary movies

Supplementary movies are available at https://doi.org/10.1017/jfm.2021.581.

Acknowledgements

J.V. thanks S.G. Llewellyn Smith for pointing out useful references.

Declaration of interests

The authors report no conflict of interest.

Appendix A. Momentum equation

We derive the momentum equation in coordinates $(\zeta ,\theta )$. The derivation illustrates the convenience of working directly with forms rather than their components, as advocated by Frankel (Reference Frankel2004). Starting from the coordinate-free momentum equation (2.1a) we compute

(A1)\begin{align} (\partial_t + \mathcal{L}_{\boldsymbol{u}}) \nu &= (\partial_t + \mathcal{L}_{\boldsymbol{u}}) (u^{\zeta} \,\textrm{d} \zeta + |g| u^{\theta} \,\textrm{d} \theta) \nonumber\\ &= (\partial_t + \mathcal{L}_{\boldsymbol{u}}) (u^{\zeta}) \,\textrm{d} \zeta + u^{\zeta} \,\textrm{d} (\mathcal{L}_{\boldsymbol{u}} \zeta) + (\partial_t + \mathcal{L}_{\boldsymbol{u}}) (|g| u^{\theta}) \,\textrm{d} \theta + |g| u^{\theta} \,\textrm{d}(\mathcal{L}_{\boldsymbol{u}} \theta) \nonumber\\ &= (\partial_t + \mathcal{L}_{\boldsymbol{u}}) (u^{\zeta}) \,\textrm{d} \zeta + u^{\zeta} \,\textrm{d} u^{\zeta} + (\partial_t + \mathcal{L}_{\boldsymbol{u}}) (|g| u^{\theta}) \,\textrm{d} \theta + |g| u^{\theta} \,\textrm{d} u^{\theta} \end{align}

and note that

(A2)\begin{align} \textrm{d} \left(p - \tfrac{1}{2} |\boldsymbol{u}|^{2}\right) &= \textrm{d} \left(p - \tfrac{1}{2} (u^{\zeta})^{2} - \tfrac{1}{2} |g| (u^{\theta})^{2}\right) \nonumber\\ &= \textrm{d} p - u^{\zeta} \,\textrm{d} u^{\zeta} - |g| u^{\theta} \,\textrm{d} u^{\theta} - \tfrac{1}{2} (u^{\theta})^{2} (\partial_\zeta |g| \,\textrm{d} \zeta + \partial_\theta |g| \,\textrm{d} \theta) \end{align}

to find the two components

(A3a)$$\begin{gather} (\partial_t + \mathcal{L}_{\boldsymbol{u}}) u^{\zeta} - \tfrac{1}{2} (u^{\theta})^{2} \partial_\zeta |g| ={-} \partial_\zeta p, \end{gather}$$
(A3b)$$\begin{gather}(\partial_t + \mathcal{L}_{\boldsymbol{u}})(|g| u^{\theta}) - \tfrac{1}{2} (u^{\theta})^{2} \partial_\theta |g| ={-} \partial_\theta p, \end{gather}$$

where $\mathcal {L}_{\boldsymbol {u}} = u^{\zeta } \partial _\zeta + u^{\theta } \partial _\theta$ since it acts on scalars.

Appendix B. Geometric derivation

The vorticity–stream function formulation can be derived without resorting to coordinates. Following Kimura (Reference Kimura1999) and Ragazzo & de Barros Viglioni (Reference Ragazzo and de Barros Viglioni2017), this is best carried out using the Hodge * operator which, in two dimensions, is determined by the following:

(B1ac)\begin{equation} *1 = \mu,\quad * \alpha = \alpha^{\sharp} \lrcorner \mu \quad \text{and}\quad *\mu=1, \end{equation}

where $\alpha$ is a 1-form and $\alpha ^{\sharp }$ its dual vector, with $\alpha =g(\alpha ^{\sharp },\cdot )$. In coordinates, we have

(B2)\begin{equation} * (\alpha_\zeta \,\textrm{d} \zeta + \alpha_\theta \,\textrm{d} \theta) = |g|^{1/2} \alpha_\zeta \,\textrm{d} \theta - |g|^{{-}1/2} \alpha_\theta \,\textrm{d} \zeta. \end{equation}

Noting that $** \alpha = - \alpha$ (e.g. Frankel Reference Frankel2004), we can rewrite (2.14) and (2.16) as

(B3a,b)\begin{equation} \nu = *\textrm{d} \psi \quad \text{and}\quad \omega = * \textrm{d} \nu = * \textrm{d} * \textrm{d} \psi, \end{equation}

which identifies the Laplacian acting on (pseudo-)scalars as ${\rm \Delta} = * \textrm {d} * \textrm {d}$.

This form of the Laplacian is also useful to compute the change in the circulation along the boundary that is introduced by dissipation. For simplicity, we take the our dissipative model to be

(B4)\begin{equation} (\partial_t + \mathcal{L}_u) \nu = \varepsilon {\rm \Delta} \nu, \end{equation}

where ${\rm \Delta} = \textrm {d}*\textrm {d}* + *\textrm {d}*\textrm {d}$ is the Laplace–de Rham operator. We emphasise that this is not the standard (Navier–Stokes) viscous dissipation which, instead, involves a different Laplacian defined via the divergence of the stress tensor (Gilbert et al. Reference Gilbert, Riedinger and Thuburn2014; Gilbert & Vanneste Reference Gilbert and Vanneste2021). Integrating over the (closed) boundary $\partial \mathcal {M}$ and noting that the outer exact differential in ${\rm \Delta}$ integrates to zero gives

(B5)\begin{equation} \frac{\mathrm{d} \varGamma}{\mathrm{d} t} = \varepsilon \int_{\partial \mathcal{M}} *\textrm{d} * \textrm{d} \nu = \varepsilon \int_{\partial \mathcal{M}} *\textrm{d} \omega, \end{equation}

leading to the coordinate expression

(B6)\begin{equation} \frac{\mathrm{d} \varGamma}{\mathrm{d} t} = \varepsilon \int_0^{\rm \pi} \left( |g|^{1/2}({a},\theta) \partial_\zeta \omega({a},\theta,t) + |g|^{1/2}(-{a},\theta) \partial_\zeta \omega(-{a},\theta,t) \right) \textrm{d} \theta \end{equation}

on using (B2) and that $\omega =0$ on $\partial \mathcal {M}$.

References

REFERENCES

Arakawa, A. 1966 Computational design for long-term numerical integration of the equations of fluid motion: two-dimensional incompressible flow. Part I. J. Comput. Phys. 1, 119143.10.1016/0021-9991(66)90015-5CrossRefGoogle Scholar
Arnold, V.I. & Khesin, B.A. 1998 Topological Methods in Hydrodynamics, Applied Mathematical Sciences, vol. 125. Springer.10.1007/b97593CrossRefGoogle Scholar
Besse, N. & Frisch, U. 2017 Geometric formulation of the Cauchy invariants for incompressible euler flow in flat and curved spaces. J. Fluid Mech. 825, 412478.10.1017/jfm.2017.402CrossRefGoogle Scholar
Boatto, S. & Koiller, J. 2015 Vortices on closed surfaces. In Geometry, Mechanics, and Dynamics (ed. G. Patrick, D.Chang, D. Holm & T. Ratiu), Fields Institute Communications, vol. 73, pp. 185–238. Springer.10.1007/978-1-4939-2441-7_10CrossRefGoogle Scholar
Couder, Y. & Basdevant, C. 1986 Experimental and numerical study of vortex couples in two-dimensional flows. J. Fluid Mech. 173, 225251.10.1017/S0022112086001155CrossRefGoogle Scholar
Couder, Y., Chomaz, J.M. & Rabaud, M. 1989 On the hydrodynamics of soap films. Physica D 37 (1), 384405.10.1016/0167-2789(89)90144-9CrossRefGoogle Scholar
Courant, R. 1940 Soap film experiments with minimal surfaces. Am. Math. Mon. 47 (3), 167174.10.1080/00029890.1940.11990957CrossRefGoogle Scholar
Dritschel, D.G. & Boatto, S. 2015 The motion of point vortices on closed surfaces. Proc. R. Soc. Lond. A 471, 20140890.Google Scholar
Frankel, T. 2004 The Geometry of Physics, 2nd edn. Cambridge University Press.Google Scholar
Gilbert, A.D., Riedinger, X. & Thuburn, J. 2014 On the form of the viscous term for two dimensional Navier–Stokes flows. Q. J. Mech. Appl. Maths 67, 205228.10.1093/qjmam/hbu004CrossRefGoogle Scholar
Gilbert, A.D. & Vanneste, J. 2018 Geometric generalised Lagrangian-mean theories. J. Fluid Mech. 839, 95134.10.1017/jfm.2017.913CrossRefGoogle Scholar
Gilbert, A.D. & Vanneste, J. 2021 A geometric look at momentum flux and stress in fluid mechanics. arXiv:1911.06613.Google Scholar
Goldstein, R.E., Moffatt, H.K., Pesci, A.I. & Ricca, R.L. 2010 Soap-film Möbius strip changes topology with a twist singularity. Proc. Natl Acad. Sci. USA 107 (51), 2197921984.10.1073/pnas.1015997107CrossRefGoogle Scholar
Gustafsson, B. 2019 Vortex motion and geometric function theory: the role of connections. Phil. Trans. R. Soc. Lond. A 377, 20180341.Google Scholar
Kimura, Y. 1999 Vortex motion on surfaces with constant curvature. Proc. R. Soc. Lond. A 455, 245259.10.1098/rspa.1999.0311CrossRefGoogle Scholar
Odehnal, B. 2016 A rational minimal Möbius strip. In Proceedings of the 17th International Conference on Geometry and Graphics, paper 70.Google Scholar
Ragazzo, C.G. & de Barros Viglioni, H.H. 2017 Hydrodynamics vortex on surfaces. J. Nonlinear Sci. 27, 16091640.10.1007/s00332-017-9380-7CrossRefGoogle Scholar
Schutz, B. 1980 Geometrical Methods of Mathematical Physics. Cambridge University Press.10.1017/CBO9781139171540CrossRefGoogle Scholar
Schwarz, G.E. 1990 The dark side of the Moebius strip. Am. Math. Mon. 97, 890897.10.1080/00029890.1990.11995680CrossRefGoogle Scholar
Starostin, E.L. & van der Heijden, G.H.M. 2007 The shape of a Möbius strip. Nat. Mater. 6, 563567.10.1038/nmat1929CrossRefGoogle ScholarPubMed
Zermelo, E. 1902 Hydrodynamische untersuchungen über die Wirbelbewegungen in einer kugelfläche. Z. Math. Phys. 47, 201237.Google Scholar
Zermelo, E. 2013 Hydrodynamical investigations of vortex motions in the surface of a sphere. In Ernst Zermelo – Collected works (ed. H.-D. Ebbinghaus & A. Kanamori), vol. II, pp. 300–391. Springer.Google Scholar
Figure 0

Figure 1. Möbius strip in $\mathbb {R}^{3}$ and its parameter space for $a=0.25$. The Möbius strip is coloured according to the value of $|g|^{1/2}$, with $|g|$ the metric determinant in (2.3), which ranges from approximately $1.5$ (dark blue) to $2.5$ (light yellow) and whose variations are associated with curvature effects.

Figure 1

Figure 2. Visualisation of the non-zero entries in the matrix discretising the Laplacian operator on the Möbius strip. The unusual antidiagonal bands in the top-right and bottom-left corners result from the twist of the Möbius strip.

Figure 2

Figure 3. Propagation of a vortex along the boundary of the Möbius strip: the vortex travels anticlockwise along the entire boundary, with a vorticity vector field $\omega \boldsymbol {N}$ (shown by the arrow in the centre of the vortex) reversing direction as the along-strip coordinate $\theta$ increases by ${\rm \pi}$ (see also supplementary movie 1 available at https://doi.org/10.1017/jfm.2021.581).

Figure 3

Figure 4. Shear instability on a Möbius strip: the (pseudo-scalar) vorticity density $\omega$ is shown at six successive time intervals, starting with a perturbation of the field (3.6) (see also movie 2).

Figure 4

Figure 5. Decaying turbulence on a Möbius strip: the (pseudo-scalar) vorticity density $\omega$ is shown at six successive equal time intervals (see also movie 3).

Vanneste supplementary movie 1

See pdf file for movie caption

Download Vanneste supplementary movie 1(Video)
Video 1.1 MB

Vanneste supplementary movie 2

See pdf file for movie caption

Download Vanneste supplementary movie 2(Video)
Video 6.6 MB

Vanneste supplementary movie 3

See pdf file for movie caption

Download Vanneste supplementary movie 3(Video)
Video 8.6 MB
Supplementary material: PDF

Vanneste supplementary material

Captions for movies 1-3

Download Vanneste supplementary material(PDF)
PDF 11.1 KB