## Abstract

Recent experiments demonstrate the importance of substrate curvature for actively forced fluid dynamics. Yet, the covariant formulation and analysis of continuum models for nonequilibrium flows on curved surfaces still poses theoretical challenges. Here, we introduce and study a generalized covariant Navier-Stokes model for fluid flows driven by active stresses in nonplanar geometries. The analytical tractability of the theory is demonstrated through exact stationary solutions for the case of a spherical bubble geometry. Direct numerical simulations reveal a curvature-induced transition from a burst phase to an anomalous turbulent phase that differs distinctly from externally forced classical 2D Kolmogorov turbulence. This new type of active turbulence is characterized by the self-assembly of finite-size vortices into linked chains of antiferromagnetic order, which percolate through the entire fluid domain, forming an active dynamic network. The coherent motion of the vortex chain network provides an efficient mechanism for upward energy transfer from smaller to larger scales, presenting an alternative to the conventional energy cascade in classical 2D turbulence.

Original language | English (US) |
---|---|

Article number | 164503 |

Journal | Physical review letters |

Volume | 120 |

Issue number | 16 |

DOIs | |

State | Published - Apr 20 2018 |

Externally published | Yes |

## ASJC Scopus subject areas

- Physics and Astronomy(all)

## Access to Document

## Other files and links

## Fingerprint

Dive into the research topics of 'Anomalous Chained Turbulence in Actively Driven Flows on Spheres'. Together they form a unique fingerprint.## Cite this

- APA
- Author
- BIBTEX
- Harvard
- Standard
- RIS
- Vancouver

*Physical review letters*,

*120*(16), [164503]. https://doi.org/10.1103/PhysRevLett.120.164503

**Anomalous Chained Turbulence in Actively Driven Flows on Spheres**. In: Physical review letters. 2018 ; Vol. 120, No. 16.

}

*Physical review letters*, vol. 120, no. 16, 164503. https://doi.org/10.1103/PhysRevLett.120.164503

**Anomalous Chained Turbulence in Actively Driven Flows on Spheres.** / Mickelin, Oscar; Słomka, Jonasz; Burns, Keaton J. et al.

Research output: Contribution to journal › Article › peer-review

TY - JOUR

T1 - Anomalous Chained Turbulence in Actively Driven Flows on Spheres

AU - Mickelin, Oscar

AU - Słomka, Jonasz

AU - Burns, Keaton J.

AU - Lecoanet, Daniel

AU - Vasil, Geoffrey M.

AU - Faria, Luiz M.

AU - Dunkel, Jörn

N1 - Funding Information: Mickelin Oscar 1 Słomka Jonasz 1 Burns Keaton J. 2 Lecoanet Daniel 3 Vasil Geoffrey M. 4 Faria Luiz M. 1 Dunkel Jörn 1 Department of Mathematics, 1 Massachusetts Institute of Technology , Cambridge, Massachusetts 02139-4307, USA Department of Physics, 2 Massachusetts Institute of Technology , Cambridge, Massachusetts 02139-4307, USA Princeton Center for Theoretical Science, 3 Princeton University , Princeton, New Jersey 08544, USA School of Mathematics and Statistics, 4 University of Sydney , Sydney, New South Wales 2006, Australia 20 April 2018 20 April 2018 120 16 164503 16 October 2017 © 2018 American Physical Society 2018 American Physical Society Recent experiments demonstrate the importance of substrate curvature for actively forced fluid dynamics. Yet, the covariant formulation and analysis of continuum models for nonequilibrium flows on curved surfaces still poses theoretical challenges. Here, we introduce and study a generalized covariant Navier-Stokes model for fluid flows driven by active stresses in nonplanar geometries. The analytical tractability of the theory is demonstrated through exact stationary solutions for the case of a spherical bubble geometry. Direct numerical simulations reveal a curvature-induced transition from a burst phase to an anomalous turbulent phase that differs distinctly from externally forced classical 2D Kolmogorov turbulence. This new type of active turbulence is characterized by the self-assembly of finite-size vortices into linked chains of antiferromagnetic order, which percolate through the entire fluid domain, forming an active dynamic network. The coherent motion of the vortex chain network provides an efficient mechanism for upward energy transfer from smaller to larger scales, presenting an alternative to the conventional energy cascade in classical 2D turbulence. Princeton Center for Theoretical Science Australian Research Council 10.13039/501100000923 DE140101960 Alfred P. Sloan Foundation 10.13039/100000879 National Science Foundation 10.13039/100000001 CBET-1510768 James S. McDonnell Foundation 10.13039/100000913 Substrate geometry profoundly affects dynamics and energy transport in complex fluids flowing far from equilibrium [1–3] . Examples range from magnetohydrodynamic turbulence on stellar surfaces [1] to the rich microscale dynamics of topological defects in active nematic vesicles [2,3] . Studying the interplay between spatial curvature and actively driven fluid flows is also essential for understanding microbial locomotion [4] , biofilm formation [5] and bioremediation [6] in soils [7] , tissues [8] , and water [9–11] . Over the past two decades, important breakthroughs have been made in characterizing active-stress driven matter flows in planar Euclidean geometries both theoretically [12–15] and experimentally [16–18] . More recently, theoretical work has begun to focus on incorporating curvature effects into active matter models [19–25] . Despite some promising progress, the hydrodynamic description of pattern-forming nonequilibrium liquids in non-Euclidean spaces continues to pose conceptual challenges, attributable to the difficulty of formulating exactly solvable continuum models and devising efficient spectral methods in curved geometries. Aiming to help improve upon these two issues, we introduce and investigate here the covariant extension of a generalized Navier-Stokes (GNS) model [26–29] describing incompressible active fluid flow on an arbitrarily curved surface. Focusing on a spherical “bubble” geometry, we derive exact stationary solutions and numerically explore the effects of curvature on the steady-state flow dynamics, using the open-source spectral code Dedalus [30] . The numerically obtained phase diagrams, energy spectra, and flux curves predict an anomalous turbulent phase when the spectral bandwidth of the active stresses becomes sufficiently narrow. This novel type of 2D turbulence supports an unexpected upward energy transfer mechanism, mediated by the large-scale collective dynamics of self-organized vortex chains, akin to actively moving antiferromagnetic spin chains. At high curvature, the anomalous turbulence transforms into a quasistationary burst phase, whereas for broadband spectral forcing the flow dynamics transitions to classical 2D Kolmogorov turbulence, accumulating energy in a few large-scale vortices. We next motivate and define the covariant GNS model for an arbitrary 2D surface; analytical and numerical results for the sphere case will be discussed subsequently. Recent experiments have investigated the collective dynamics of swimming bacteria [16] and algae [31] in thin quasi-2D soap films held by a coplanar wire frame. Generalizing to non-Euclidean geometries [2,3] , which can be realized with soap bubbles or curved wire frames [32] , we consider here a freestanding nonplanar 2D film in which the fluid flow is driven by active stresses, as in suspensions of swimming bacteria [33,34] or ATP-driven microtubule networks [35,36] . On a curved manifold, the fluid velocity field components v a satisfy incompressibility and Cauchy momentum conservation [37,38] , (1) ∇ a v a = 0 , (1a) ∂ t v a + v b ∇ b v a = ∇ a σ + ∇ b T a b , (1b) where ∇ b v a denotes the covariant derivative of v a , a , b = 1 , 2 and σ is the (surface) tension. The stress tensor T a b includes passive and active contributions from the solvent fluid viscosity and the stresses exerted by the microswimmers on the fluid. Below, we study the covariant version of the linear active-stress model [26–29] T a b = f ( ∇ 2 ) ( ∇ a u b + ∇ b u a ) , f ( ∇ 2 ) = Γ 0 - Γ 2 ∇ 2 + Γ 4 ∇ 2 ∇ 2 , (1c) where ∇ 2 = ∇ a ∇ a is the tensor Laplacian. In qualitative agreement with experimental observations for active suspensions [16,18,36,39] , the polynomial ansatz for f in Eq. (1c) generates vortices of characteristic size Λ and growth time τ , provided that Γ 2 < 0 , which introduces a bandwidth κ of linearly unstable modes [28] . General mathematical stability considerations demand Γ 0 , Γ 4 > 0 . The phenomenological model defined in Eq. (1) is minimal in the sense that it assumes the active stresses create to leading order a linear instability, while neglecting energy transfer within the active component. As verified in Ref. [29] , the linear active-stress model, Eq. (1c) , suffices to quantitatively reproduce the experimentally measured velocity distributions and flow correlations in 3D bacterial [18] and ATP-driven microtubule [36] suspensions. More generally, closely related GNS models have also been studied in the context of soft-mode turbulence and seismic waves [26,27] . Numerical solutions of the GNS Eqs. (1) show significant phenomenological similarities with magnetohydrodynamic (MHD) flows driven by electromagnetic stresses [40] , suggesting that the results below may also apply to astrophysical systems. Exact stationary solutions of Eqs. (1) for a sphere of radius R can be constructed from the vorticity-stream function formulation (Supplemental Material [41] ) (2) Δ ψ = - ω , (2a) ∂ t ω + { ω , ψ } = f ( Δ + 4 K ) ( Δ + 2 K ) ω , (2b) where ψ and ω are the stream function and vorticity. The advection term in spherical coordinates ( θ , ϕ ) reads { ω , ψ } = ( ∂ θ ω ∂ ϕ ψ - ∂ ϕ ω ∂ θ ψ ) / ( R 2 sin θ ) . K = R - 2 is the Gaussian curvature and Δ the standard spherical Laplacian. Since the spherical harmonics Y ℓ m diagonalize the Laplacian, Δ Y ℓ m = - R - 2 ℓ ( ℓ + 1 ) Y ℓ m for integers ℓ , m such that ℓ ≥ 0 and | m | ≤ ℓ , an arbitrary superposition ψ = ∑ | m | ≤ ℓ ψ m ℓ Y ℓ m (3) solves the system Eq. (2) exactly, provided that the eigenvalue ℓ is an integer root of f [ - ℓ ( ℓ + 1 ) + 4 ] = 0 (Supplemental Material [41] ). As usual, the velocity field is tangent to the level sets of the stream function. Two particular exact solutions are shown in Fig. 1 . The first example, Fig. 1(a) , is reminiscent of the square lattice solutions found earlier in the flat 2D case [28] . The second example in Fig. 1(b) illustrates a flow field with fivefold symmetry, obtained by applying the superposition procedure of Ref. [46] . Although these exact solutions are not stable, they provide some useful intuition about the instantaneous flow patterns expected in dynamical simulations (Fig. 2 ), similar to exact coherent structures [47] in conventional turbulence [48] . 1 10.1103/PhysRevLett.120.164503.f1 FIG. 1. Stationary solutions of Eqs. (2) are superpositions of the form Eq, (3) with f [ - ℓ ( ℓ + 1 ) + 4 ] = 0 . (a) An exact stationary solution with ℓ = 6 which is also approximately realized as a transient state in the time-dependent burst solution of Fig. 2 (movie 1). (b) Complex symmetric solutions can be constructed by choosing the expansion coefficients ψ m ℓ accordingly [46] . In both panels, the stream functions are normalized by their maxima; see Supplemental Material [41] for coefficients ψ m ℓ . 2 10.1103/PhysRevLett.120.164503.f2 FIG. 2. Phase diagrams [(a),(b)] and representative still images [(c)–(e)] from simulations showing quasistationary burst dynamics ( B phase), anomalous vortex-network turbulence ( A phase), and classical 2D turbulence ( T phase). (a), (b) The A and T phase are approximately separated by the condition κ Λ = 1 (vertical dashed line) and differ by the average number of vortices (a), the branch geometry of the tension field (b), and the energy spectra (Fig. 3 ). The B phase arises for narrow-band energy injection κ R ≲ 1 when only a single ℓ mode is active (region right below the dash-dotted line); decreasing κ further gives a passive fluid (white region). (c)–(e) Top: Instantaneous vorticity fields normalized by their maxima. Bottom: Surface tension fields normalized by the maximum deviation from the mean. (c) Quasistationary preburst state from movie 1 resembling the exact solution in Fig. 1(a) ; see movies 2 and 3 for additional examples labeled by * in panel (a). (d) For subcritical curvature and intermediate energy injection bandwidths, R - 1 < κ < Λ - 1 , the flows develop a percolating vortex-chain network structure (movie 4), with accumulation of tension and vorticity along the edges. (e) For broadband energy injection κ Λ > 1 , smaller eddies merge to create larger vortices, as typical of classical 2D turbulence (movie 5). Parameters: (a) α ω = 0.5 ; (c) R / Λ = 2 , τ = 4.9 s , κ Λ = 0.29 ; (d) R / Λ = 10 , τ = 14.9 s , κ Λ = 0.5 ; (e) R / Λ = 10 , τ = 11.7 s , κ Λ = 2.0 . Panels (a),(b) show steady-state time averages over [ 50 τ , 100 τ ] . Solid curves in (c)–(e) indicate stream lines of the velocity fields. To find and analyze time-dependent solutions of Eqs. (1) , we performed numerical simulations using Dedalus [30] , an open-source framework for solving differential equations with spectral methods. The Eqs. (1) were solved directly as a coupled partial differential-algebraic system for the scalar tension σ and vector velocity v a . To spatially discretize the system, we used spin-weighted spherical harmonics, which are a parameterized family of basis functions that correctly capture the analytical behavior of spin-weighted functions on the sphere (Supplemental Material [41] ). Under this spectral expansion, the system Eq. (1) is reduced to a set of coupled ordinary differential-algebraic equations for the time evolution of the expansion coefficients. We solve these equations using mixed implicit-explicit time stepping, in which the linear terms of the evolution equations are integrated implicitly, the linear constraints are enforced implicitly, and the nonlinear terms are integrated explicitly. This allows us to simultaneously evolve the velocity field while enforcing the incompressibility constraint, and with a time step that is limited by the advective Courant-Friedrichs-Lewy time condition rather than the diffusive time at any scale. The parameters ( Γ 0 , Γ 2 , Γ 4 ) in Eqs. (1) define a characteristic time scale τ , a characteristic vortex diameter Λ , and a characteristic spectral bandwidth κ , which can be directly inferred from experimental data [29] ; explicit expressions are derived in the Supplemental Material [41] . Given a sphere of radius R , fixing ( τ , Λ , κ ) uniquely determines the parameters ( Γ 0 , Γ 2 , Γ 4 ). To explore the interplay between curvature and activity, we run 351 simulations, using R / Λ ∈ [ 2 , 10 ] and κ · Λ ∈ [ 0.1 , 2.0 ] . Typical vortex diameters for bacterial and microtubule suspensions are Λ ∼ 50 – 100 μ m with τ of the order of seconds [16,18,36,39] . Time steps were in the range [ 5 × 10 - 4 τ , 5 × 10 - 3 τ ] with a total simulation time 100 τ , allowing the system to fully develop its dynamics after an initial relaxation phase during which active stresses inject energy until the viscous dissipation and activity balance on average. In the remainder, it will be convenient to regard Λ as reference length and compare the flow topologies across the ( κ , R ) parameter plane. Our simulations reveal three qualitatively distinct flow regimes (Fig. 2 ): a quasistationary burst phase for κ R ≲ 1 [domain B in Fig. 2(a) ; movies 1–3], an anomalous turbulence for R - 1 < κ < Λ - 1 [domain A in Fig. 2(a) ; movie 4], and normal 2D turbulence for κ Λ > 1 [domain T in Fig. 2(a) ; movie 5]. Representative vorticity and tension fields from the corresponding steady-state dynamics are shown in Figs. 2(c)–2(e) . In the B phase, the energy injection bandwidth κ is close to the wave number spacing set by the sphere curvature R - 1 , leaving only a single active wave number ℓ . Decreasing κ further completely suppresses active modes resulting in globally damped fluid motion [white domain in Fig. 2(a) ]. The B phase is characterized by the formation of intermittent quasistationary flow patterns that lie in the vicinity of the exact stationary solutions, Eq. (3) , cf. Fig. 1(a) and Fig. 2(c) . Once formed, the amplitude of these flow patterns grows exponentially (Fig. S3) until nonlinear advection becomes dominant and eventually causes energy to be released through a rapid burst. Afterwards, the dynamics becomes quasilinear again with the flow settling into a new quasistationary pattern. These burst cycles are continuously repeated (movies 1–3). The two turbulent phases A and T in Fig. 2(a) can be distinguished through topological, geometric, and spectral measures. We demonstrate this by determining the topology of the vorticity fields, the geometry of the high-tension domains, and the energy spectra for each simulation after flows had reached the chaotic steady state. To study the vortex topology, we fix a threshold α ω ∈ [ 0 , 1 ] and identify regions in which the vorticity is larger (or smaller) than α ω times the maximum (or minimum) vorticity (Supplemental Material [41] ). This thresholding divides the sphere into patches of high absolute vorticity (Fig. S1). The number of connected domains, given by the zeroth Betti number, counts the vortices in the system. For a fixed pair ( κ , R ) , we denote the vortex number at time t by N ω ( κ , R ; t ) . Although more sophisticated methods for vortex detection exist [49] , the thresholding criterion proved to be sufficient for our analysis (Fig. S2). To normalize vortex numbers across the parameter space, we fix a reference value κ * = 0.3 / Λ . With this, we can define a normalized Betti number as Betti ω ( κ , R ) = ⟨ N ω ( κ , R ; t ) - N ω ( κ * , R ; t ) ⟩ ⟨ N ω ( κ * , R ; t ) ⟩ , (4) where the time average ⟨ ⋯ ⟩ is taken after the initial relaxation period. Intuitively, large values of Betti ω indicate many vortices of comparable circulation, whereas small values suggest the presence of a few dominant eddies. The variation of Betti ω in the ( κ , R ) -parameter plane is color coded in Fig. 2(a) . In the anomalous turbulent A phase, vortices of diameter ≈ Λ eventually cover the surface of the sphere, with stronger vortices forming chains of anti-ferromagnetic order [Fig. 2(d) top; movie 4]. By contrast, in the T phase characterized by broadband energy injection κ > Λ - 1 , smaller eddies merge to create a small number of larger vortices, as typical of classical 2D turbulence [50] [Fig. 2(e) top; movie 5]. Interestingly, the A phase shares phenomenological similarities with the low-entropy states found in quasi-2D superfluid models [51] , while the vortex condensation in the T phase corresponds approximately to the negative “temperature” regime in Onsager’s statistical hydrodynamics [52] . Moreover, the upper region of Fig. 2(a) , which corresponds to the small-curvature limit R / Λ ≫ 1 , suggests that the two phases extend to planar geometries, provided boundary effects remain negligible. To obtain a more detailed geometric characterization of the turbulent A and T phases, we next consider the corresponding tension fields. Analogously to the case of vorticity above, we focus on regions where the local tension σ ( t , x ) is larger than the instantaneous global mean value. For each connected component of the identified high-tension regions, we denote by A its total area and by ∂ A its total boundary area in pixels. The ratio ∂ A / A is a measure of chainlike structures in the tension fields, a large value signaling a highly branched structure, whereas smaller values indicate less branching. Denoting the instantaneous sum of the ratios ∂ A / A over all connected high-tension domains by A σ ( κ , R ; t ) , a normalized branching index can then be defined by (Supplemental Material [41] ) Branch σ ( κ , R ) = ⟨ A σ ( κ , R ; t ) - A σ ( κ * , R ; t ) ⟩ ⟨ A σ ( κ * , R ; t ) ⟩ , (5) where the time average is again taken after the initial relaxation. As evident from the phase diagram in Fig. 2(b) and the corresponding tension fields in Figs. 2(d) and 2(e) , and movies 4 and 5, the geometric characterization confirms the existence of an anomalous turbulent phase, in which vortices combine to form percolating dynamic networks with high tension being localized along the edges [Fig. 2(d) bottom; movie 4]. To compare the energy transport in the anomalous turbulent phase with classical 2D turbulence, we analyze the energy spectra and fluxes for the A and T phases. Expanding in spherical harmonics, ψ = ∑ m , ℓ ψ m ℓ Y ℓ m , the energy of mode ℓ is E ( ℓ ) = ∑ | m | ≤ ℓ ℓ ( ℓ + 1 ) | ψ m ℓ | 2 . The corresponding mean energy flux across ℓ in the statistically stationary state is obtained as (Supplemental Material [41] ) Π ( ℓ ) = - 2 ∑ ℓ ′ ≥ ℓ f [ 4 - ℓ ′ ( ℓ ′ + 1 ) ] [ 2 - ℓ ′ ( ℓ ′ + 1 ) ] ⟨ E ℓ ′ ⟩ , (6) where f is the polynomial defined in Eq. (1c) . Figure 3 shows the numerically obtained energy spectra E ( ℓ ) and fluxes Π ( ℓ ) for four active bandwidths κ . In all four cases, the kinetic energy produced in the injection range ( ℓ ∼ π R / Λ ) propagates to both large ( ℓ < π R / Λ ) and small ( ℓ > π R / Λ ) scales, as indicated by negative and positive values of Π ( ℓ ) , respectively. Energy transfer to large scales is a prominent feature of classical 2D turbulence [50,53,54] and our results show that it also occurs in active turbulence. However, the transfer mechanisms can be dramatically different, as already implied by the preceding analysis of the vorticity and tension fields. For broadband spectral forcing κ Λ ≫ 1 , the classical 2D turbulence picture of vortex mergers and energy condensation at large scales prevails [Fig. 2(e) ; movie 3]. For κ Λ ≲ 2 the spectrum follows a k - 1 scaling, indicating the formation of a dilute-vortex system [55] . For even larger values of κ , additional large-scale dissipation is needed to bound the upward energy transfer, in which case the spectrum is expected to approach the Kolmogorov k - 5 / 3 scaling [50] . By contrast, for narrow-band driving κ Λ ≲ 1 , the upward energy transfer is realized through the coherent motion of high-tension vortex chains. Interestingly, only this anomalous type of inverse energy cascade appears to persist in 3D active bulk fluids [29] , where it is sustained by spontaneous chiral symmetry breaking [56] . 3 10.1103/PhysRevLett.120.164503.f3 FIG. 3. Time-averaged energy spectra and fluxes indicate two qualitatively different types of upward energy transport. (a) For narrow-band energy injection κ Λ < 1 , the energy spectrum exhibits a peak corresponding to the dominant vortex size Λ (red curve). For broadband injection κ Λ ∼ 2 , the spectra decay monotonically (blue and black curves). (b) In all four examples, the fluxes confirm inverse energy transport, albeit with different origins. For broadband energy injection (blue and black curves), the upward energy flux to larger scales is due to vortex mergers [Fig. 2(e) ; movie 5]. By contrast, for narrow-band injection (red curve), a relatively stronger upward energy flux arises from the collective motion of vortex chains [Fig. 2(d) ; movie 4]. The shaded regions indicate the energy injection ranges with colors matching those of the corresponding curves, respectively. Parameters: R / Λ = 10 for a unit sphere, τ = 11.7 s , time step 5 × 10 - 4 τ , total simulation time 500 τ . Spectra and fluxes were determined after relaxation by averaging over [ 150 τ , 500 τ ] . For κ Λ ≫ 1 , energy steadily accumulates at larger scales and the absence of a large-scale dissipative mechanism leads to a divergent total enstrophy and kinetic energy on the sphere. In summary, we have presented analytical and numerical solutions for generalized Navier-Stokes equations describing actively driven nonequilibrium flows on a sphere. Our calculations predict that spectrally localized active stresses can induce a novel turbulent phase, in which finite-size vortices self-organize into chain complexes of antiferromagnetic order that percolate through the surface [57] . The collective motion of these chain networks enables a significant upward energy transport and may thus provide a basis for efficient fluid mixing in quasi-2D active and magnetohydrodynamic flows. Future generalizations to rotating spheres could thus promise insights into pattern formation in planetary and stellar atmospheres [58] . D. L. is supported by a Princeton Center for Theoretical Science Fellowship, and a Lyman Spitzer Jr. Fellowship. G. M. V. acknowledges support from The Australian Research Council Project No. DE140101960. This work was supported by an Alfred P. Sloan Research Fellowship (J. D.), an Edmund F. Kelly Research Award (J. D.), NSF Grant No. CBET-1510768 (J. S. and J. D.) and a James S. McDonnell Foundation Complex Systems Scholar Award (J. D.). O. M. and J. S. contributed equally to this work. O. M. performed the numerical simulations and J. S. the analytical calculations. K. J. B, D. L., and G. M. V. developed the numerical simulations. [1] 1 S. Cranmer and A. van Ballegooijen , Astrophys. J. Suppl. Ser. 156 , 265 ( 2005 ). APJSA2 1538-4365 10.1086/426507 [2] 2 F. C. Keber , E. Loiseau , T. Sanchez , S. J. DeCamp , L. Giomi , M. J. Bowick , M. C. Marchetti , Z. Dogic , and A. R. Bausch , Science 345 , 1135 ( 2014 ). SCIEAS 0036-8075 10.1126/science.1254784 [3] 3 R. Zhang , Y. Zhou , M. Rahimi , and J. J. de Pablo , Nat. Commun. 7 , 13483 ( 2016 ). NCAOBW 2041-1723 10.1038/ncomms13483 [4] 4 O. Sipos , K. Nagy , R. Di Leonardo , and P. Galajda , Phys. Rev. Lett. 114 , 258104 ( 2015 ). PRLTAO 0031-9007 10.1103/PhysRevLett.114.258104 [5] 5 Y.-W. Chang , A. A. Fragkopoulos , S. M. Marquez , H. D. Kim , T. E. Angelini , and A. Fernandez-Nieves , New J. Phys. 17 , 033017 ( 2015 ). NJOPFM 1367-2630 10.1088/1367-2630/17/3/033017 [6] 6 K. Das and A. K. Mukherjee , Bioresour. Technol. 98 , 1339 ( 2007 ). BIRTEB 0960-8524 10.1016/j.biortech.2006.05.032 [7] 7 H. C. Bold , Bulletin of the Torrey Botanical Club 76 , 101 ( 1949 ). BTBCAL 0040-9618 10.2307/2482218 [8] 8 J. W. Costerton , P. S. Stewart , and E. P. Greenberg , Science 284 , 1318 ( 1999 ). SCIEAS 0036-8075 10.1126/science.284.5418.1318 [9] 9 E. Rosenberg , R. Legmann , A. Kushmaro , R. Taube , E. Adler , and E. Z. Ron , Biodegradation 3 , 337 ( 1992 ). BIODEG 0923-9820 10.1007/BF00129092 [10] 10 M. Rosenberg , FEMS Microbiol. Lett. 262 , 129 ( 2006 ). FMLED7 0378-1097 10.1111/j.1574-6968.2006.00291.x [11] 11 E. Rosenberg , E. F. DeLong , S. Lory , E. Stackebrandt , and F. Thompson , in The Prokaryotes ( Springer , Berlin, Heidelberg, 2013 ), pp. 201–214 . [12] 12 M. C. Marchetti , J. Joanny , S. Ramaswamy , T. Liverpool , J. Prost , M. Rao , and R. A. Simha , Rev. Mod. Phys. 85 , 1143 ( 2013 ). RMPHAT 0034-6861 10.1103/RevModPhys.85.1143 [13] 13 C. W. Wolgemuth , Biophys. J. 95 , 1564 ( 2008 ). BIOJAU 0006-3495 10.1529/biophysj.107.118257 [14] 14 D. Saintillan and M. J. Shelley , Phys. Fluids 20 , 123304 ( 2008 ). PHFLE6 1070-6631 10.1063/1.3041776 [15] 15 T. Brotto , J.-B. Caussin , E. Lauga , and D. Bartolo , Phys. Rev. Lett. 110 , 038101 ( 2013 ). PRLTAO 0031-9007 10.1103/PhysRevLett.110.038101 [16] 16 A. Sokolov , I. S. Aranson , J. O. Kessler , and R. E. Goldstein , Phys. Rev. Lett. 98 , 158102 ( 2007 ). PRLTAO 0031-9007 10.1103/PhysRevLett.98.158102 [17] 17 M. Wu , J. W. Roberts , S. Kim , D. L. Koch , and M. P. DeLisa , Appl. Environ. Microbiol. 72 , 4987 ( 2006 ). AEMIDF 0099-2240 10.1128/AEM.00158-06 [18] 18 J. Dunkel , S. Heidenreich , K. Drescher , H. H. Wensink , M. Bär , and R. E. Goldstein , Phys. Rev. Lett. 110 , 228102 ( 2013 ). PRLTAO 0031-9007 10.1103/PhysRevLett.110.228102 [19] 19 L. M. Janssen , A. Kaiser , and H. Löwen , Sci. Rep. 7 , 5667 ( 2017 ). SRCEC3 2045-2322 10.1038/s41598-017-05569-6 [20] 20 G. Salbreux and F. Jülicher , Phys. Rev. E 96 , 032404 ( 2017 ). PRESCM 2470-0045 10.1103/PhysRevE.96.032404 [21] 21 S. Henkes , M. C. Marchetti , and R. Sknepnek , arXiv:1705.05166 . [22] 22 D. Khoromskaia and G. P. Alexander , New J. Phys. 19 , 103043 ( 2017 ). NJOPFM 1367-2630 10.1088/1367-2630/aa89aa [23] 23 F. Alaimo , C. Köhler , and A. Voigt , Sci. Rep. 7 , 5211 ( 2017 ). SRCEC3 2045-2322 10.1038/s41598-017-05612-6 [24] 24 Y. Fily , A. Baskaran , and M. F. Hagan , Soft Matter 10 , 5609 ( 2014 ). SMOABF 1744-683X 10.1039/C4SM00975D [25] 25 Y. Fily , A. Baskaran , and M. F. Hagan , Phys. Rev. E 91 , 012125 ( 2015 ). PRESCM 1539-3755 10.1103/PhysRevE.91.012125 [26] 26 I. A. Beresnev and V. N. Nikolaevskiy , Physica (Amsterdam) 66D , 1 ( 1993 ). PDNPDT 0167-2789 10.1016/0167-2789(93)90217-O [27] 27 M. I. Tribelsky and K. Tsuboi , Phys. Rev. Lett. 76 , 1631 ( 1996 ). PRLTAO 0031-9007 10.1103/PhysRevLett.76.1631 [28] 28 J. Słomka and J. Dunkel , Phys. Rev. Fluids 2 , 043102 ( 2017 ). PLFHBR 2469-990X 10.1103/PhysRevFluids.2.043102 [29] 29 J. Słomka and J. Dunkel , Proc. Natl. Acad. Sci. U.S.A. 114 , 2119 ( 2017 ). PNASA6 0027-8424 10.1073/pnas.1614721114 [30] 30 K. J. Burns , G. M. Vasil , J. S. Oishi , D. Lecoanet , B. P. Brown , and E. Quataert (to be published). [31] 31 J. S. Guasto , K. A. Johnson , and J. P. Gollub , Phys. Rev. Lett. 105 , 168102 ( 2010 ). PRLTAO 0031-9007 10.1103/PhysRevLett.105.168102 [32] 32 R. E. Goldstein , J. McTavish , H. K. Moffatt , and A. I. Pesci , Proc. Natl. Acad. Sci. U.S.A. 111 , 8339 ( 2014 ). PNASA6 0027-8424 10.1073/pnas.1406385111 [33] 33 E. Lauga and T. R. Powers , Rep. Prog. Phys. 72 , 096601 ( 2009 ). RPPHAG 0034-4885 10.1088/0034-4885/72/9/096601 [34] 34 K. Drescher , J. Dunkel , L. H. Cisneros , S. Ganguly , and R. E. Goldstein , Proc. Natl. Acad. Sci. U.S.A. 108 , 10940 ( 2011 ). PNASA6 0027-8424 10.1073/pnas.1019079108 [35] 35 S. Ramaswamy , Annu. Rev. Condens. Matter Phys. 1 , 323 ( 2010 ). ARCMCX 1947-5454 10.1146/annurev-conmatphys-070909-104101 [36] 36 T. Sanchez , D. T. N. Chen , S. J. DeCamp , M. Heymann , and Z. Dogic , Nature (London) 491 , 431 ( 2012 ). NATUAS 0028-0836 10.1038/nature11591 [37] 37 L. E. Scriven , Chem. Eng. Sci. 12 , 98 ( 1960 ). CESCAC 0009-2509 10.1016/0009-2509(60)87003-0 [38] 38 R. Aris , Vectors, Tensors and the Basic Equations of Fluid Mechanics ( Dover Publications, Inc. , New York, 1989 ). [39] 39 C. Dombrowski , L. Cisneros , S. Chatkaew , R. E. Goldstein , and J. O. Kessler , Phys. Rev. Lett. 93 , 098103 ( 2004 ). PRLTAO 0031-9007 10.1103/PhysRevLett.93.098103 [40] 40 G. M. Vasil and M. G. P. Cassell (to be published). [41] 41 See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.120.164503 which contains Refs. [42–45], for details of derivations and numerical implementation. [42] 42 G. Schwarz , Hodge Decomposition—A Method for Solving Boundary Value Problems ( Springer-Verlag , Berlin, 1995 ). [43] 43 E. Delay , Manuscr. Math. 123 , 147 ( 2007 ). MSMHB2 1432-1785 10.1007/s00229-007-0089-z [44] 44 A. L. Besse , Einstein Manifolds ( Springer Science & Business Media , New York, 2007 ). [45] 45 V. I. Arnold and B. A. Khesin , Topological Methods in Hydrodynamics ( Springer Science & Business Media , New York, 1999 ), Vol. 125 . [46] 46 W. Prandl , P. Schiebel , and K. Wulf , Acta Crystallogr. Sect. A 52 , 171 ( 1996 ). ACACEQ 0108-7673 10.1107/S0108767395012578 [47] 47 F. Waleffe , J. Fluid Mech. 435 , 93 ( 2001 ). JFLSA7 0022-1120 10.1017/S0022112001004189 [48] 48 F. Waleffe , in Turbulence and Interactions , edited by M. Deville , T. H. Lê , and P. Sagaut , Notes on Numerical Fluid Mechanics and Multidisciplinary Design , Vol. 105 ( Springer , Berlin, Heidelberg, 2009 ). [49] 49 M. Jiang , R. Machiraju , and D. Thompson , in The Visualization Handbook , edited by C. D. Hansen and C. R. Johnson ( Academic Press , New York, 2011 ). [50] 50 G. Boffetta and R. E. Ecke , Annu. Rev. Fluid Mech. 44 , 427 ( 2012 ). ARVFA3 0066-4189 10.1146/annurev-fluid-120710-101240 [51] 51 T. Simula , M. J. Davis , and K. Helmerson , Phys. Rev. Lett. 113 , 165302 ( 2014 ). PRLTAO 0031-9007 10.1103/PhysRevLett.113.165302 [52] 52 L. Onsager , Nuovo Cimento Suppl. 6 , 279 ( 1949 ). NUCUAF 0550-3868 10.1007/BF02780991 [53] 53 R. H. Kraichnan and D. Montgomery , Rep. Prog. Phys. 43 , 547 ( 1980 ). RPPHAG 0034-4885 10.1088/0034-4885/43/5/001 [54] 54 S. Sukoriansky , B. Galperin , and N. Dikovskaya , Phys. Rev. Lett. 89 , 124501 ( 2002 ). PRLTAO 0031-9007 10.1103/PhysRevLett.89.124501 [55] 55 T. Kusumura , H. Takeuchi , and M. Tsubota , J. Low Temp. Phys. 171 , 563 ( 2013 ). JLTPAC 0022-2291 10.1007/s10909-012-0827-9 [56] 56 L. Biferale , S. Musacchio , and F. Toschi , Phys. Rev. Lett. 108 , 164501 ( 2012 ). PRLTAO 0031-9007 10.1103/PhysRevLett.108.164501 [57] 57 A. Doostmohammadi , T. N. Shendruk , K. Thijssen , and J. M. Yeomans , Nat. Commun. 8 , 15326 ( 2017 ). NCAOBW 2041-1723 10.1038/ncomms15326 [58] 58 L. N. Fletcher , Science 319 , 79 ( 2008 ). SCIEAS 0036-8075 10.1126/science.1149514 Publisher Copyright: © 2018 American Physical Society.

PY - 2018/4/20

Y1 - 2018/4/20

N2 - Recent experiments demonstrate the importance of substrate curvature for actively forced fluid dynamics. Yet, the covariant formulation and analysis of continuum models for nonequilibrium flows on curved surfaces still poses theoretical challenges. Here, we introduce and study a generalized covariant Navier-Stokes model for fluid flows driven by active stresses in nonplanar geometries. The analytical tractability of the theory is demonstrated through exact stationary solutions for the case of a spherical bubble geometry. Direct numerical simulations reveal a curvature-induced transition from a burst phase to an anomalous turbulent phase that differs distinctly from externally forced classical 2D Kolmogorov turbulence. This new type of active turbulence is characterized by the self-assembly of finite-size vortices into linked chains of antiferromagnetic order, which percolate through the entire fluid domain, forming an active dynamic network. The coherent motion of the vortex chain network provides an efficient mechanism for upward energy transfer from smaller to larger scales, presenting an alternative to the conventional energy cascade in classical 2D turbulence.

AB - Recent experiments demonstrate the importance of substrate curvature for actively forced fluid dynamics. Yet, the covariant formulation and analysis of continuum models for nonequilibrium flows on curved surfaces still poses theoretical challenges. Here, we introduce and study a generalized covariant Navier-Stokes model for fluid flows driven by active stresses in nonplanar geometries. The analytical tractability of the theory is demonstrated through exact stationary solutions for the case of a spherical bubble geometry. Direct numerical simulations reveal a curvature-induced transition from a burst phase to an anomalous turbulent phase that differs distinctly from externally forced classical 2D Kolmogorov turbulence. This new type of active turbulence is characterized by the self-assembly of finite-size vortices into linked chains of antiferromagnetic order, which percolate through the entire fluid domain, forming an active dynamic network. The coherent motion of the vortex chain network provides an efficient mechanism for upward energy transfer from smaller to larger scales, presenting an alternative to the conventional energy cascade in classical 2D turbulence.

UR - http://www.scopus.com/inward/record.url?scp=85045892660&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85045892660&partnerID=8YFLogxK

U2 - 10.1103/PhysRevLett.120.164503

DO - 10.1103/PhysRevLett.120.164503

M3 - Article

C2 - 29756929

AN - SCOPUS:85045892660

VL - 120

JO - Physical Review Letters

JF - Physical Review Letters

SN - 0031-9007

IS - 16

M1 - 164503

ER -