skip to main content


Title: Hyperuniformity and phase enrichment in vortex and rotor assemblies
Abstract

Ensembles of particles rotating in a two-dimensional fluid can exhibit chaotic dynamics yet develop signatures of hidden order. Such rotors are found in the natural world spanning vastly disparate length scales — from the rotor proteins in cellular membranes to models of atmospheric dynamics. Here we show that an initially random distribution of either driven rotors in a viscous membrane, or ideal vortices with minute perturbations, spontaneously self assemble into a distinct arrangement. Despite arising from drastically different physics, these systems share a Hamiltonian structure that sets geometrical conservation laws resulting in prominent structural states. We find that the rotationally invariant interactions isotropically suppress long-wavelength fluctuations — a hallmark of a disordered hyperuniform material. With increasing area fraction, the system orders into a hexagonal lattice. In mixtures of two co-rotating populations, the stronger population will gain order from the other and both will become phase enriched. Finally, we show that classical 2D point vortex systems arise as exact limits of the experimentally accessible microscopic membrane rotors, yielding a new system through which to study topological defects.

 
more » « less
Award ID(s):
2004469
NSF-PAR ID:
10362701
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Nature Communications
Volume:
13
Issue:
1
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    This article is devoted to stationary solutions of Euler’s equation on a rotating sphere, and to their relevance to the dynamics of stratospheric flows in the atmosphere of the outer planets of our solar system and in polar regions of the Earth. For the Euler equation, under appropriate conditions, rigidity results are established, ensuring that the solutions are either zonal or rotated zonal solutions. A natural analogue of Arnold’s stability criterion is proved. In both cases, the lowest mode Rossby–Haurwitz stationary solutions (more precisely, those whose stream functions belong to the sum of the first two eigenspaces of the Laplace-Beltrami operator) appear as limiting cases. We study the stability properties of these critical stationary solutions. Results on the local and global bifurcation of non-zonal stationary solutions from classical Rossby–Haurwitz waves are also obtained. Finally, we show that stationary solutions of the Euler equation on a rotating sphere are building blocks for travelling-wave solutions of the 3D system that describes the leading order dynamics of stratospheric planetary flows, capturing the characteristic decrease of density and increase of temperature with height in this region of the atmosphere.

     
    more » « less
  2. View Video Presentation: https://doi.org/10.2514/6.2022-3218.vid The ability to accurately and rapidly assess unsteady interactional aerodynamics is a shortcoming and bottleneck in the design of various next-generation aerospace systems: from electric vertical takeoff and landing (eVTOL) aircraft to airborne wind energy (AWE) and wind farms. In this study, we present a meshless CFD framework based on the reformulated vortex particle method (rVPM) for the analysis of complex interactional aerodynamics. The rVPM is a large eddy simulation (LES) solving the Navier-Stokes equations in their vorticity form. It uses a meshless Lagrangian scheme, which not only avoids the hurdles of mesh generation, but it also conserves the vortical structure of wakes over long distances with minimal numerical dissipation, while being 100x faster than conventional mesh-based LES. Wings and rotating blades are introduced in the computational domain through actuator line and actuator surface models. Simulations are coupled with an aeroacoustics solver to predict tonal and broadband noise radiated by rotors. The framework, called FLOWUnsteady, is hereby released as an open-source code and extensively validated. Validation studies published in previous work by the authors are summarized, showcasing rotors across operating conditions with a rotor in hover, propellers, a wind turbine, and two side-by-side rotors in hover. Validation of rotor-wing interactions is presented simulating a tailplane with tip-mounted propellers and a blown wing with propellers mounted mid-span. The capabilities of the framework are showcased through the simulation of a tiltwing eVTOL vehicle and an AWE wind-harvesting aircraft, featuring rotors with variable RPM, variable pitch, tilting of wings and rotors, non-trivial flight paths, and complex aerodynamic interactions. 
    more » « less
  3. null (Ed.)
    In this study, we investigate and develop scaling laws as a function of external non-dimensional control parameters for heat and momentum transport for non-rotating, slowly rotating and rapidly rotating turbulent convection systems, with the end goal of forging connections and bridging the various gaps between these regimes. Two perspectives are considered, one where turbulent convection is viewed from the standpoint of an applied temperature drop across the domain and the other with a viewpoint in terms of an applied heat flux. While a straightforward transformation exist between the two perspectives indicating equivalence, it is found the former provides a clear set of connections that bridge between the three regimes. Our generic convection scalings, based upon an Inertial-Archimedean balance, produce the classic diffusion-free scalings for the non-rotating limit (NRL) and the slowly rotating limit (SRL). This is characterized by a free-falling fluid parcel on the global scale possessing a thermal anomaly on par with the temperature drop across the domain. In the rapidly rotating limit (RRL), the generic convection scalings are based on a Coriolis-Inertial-Archimedean (CIA) balance, along with a local fluctuating-mean advective temperature balance. This produces a scenario in which anisotropic fluid parcels attain a thermal wind velocity and where the thermal anomalies are greatly attenuated compared to the total temperature drop. We find that turbulent scalings may be deduced simply by consideration of the generic non-dimensional transport parameters --- local Reynolds $Re_\ell = U \ell /\nu$; local P\'eclet $Pe_\ell = U \ell /\kappa$; and Nusselt number $Nu = U \vartheta/(\kappa \Delta T/H)$ --- through the selection of physically relevant estimates for length $\ell$, velocity $U$ and temperature scales $\vartheta$ in each regime. Emergent from the scaling analyses is a unified continuum based on a single external control parameter, the convective Rossby number\JMA{,} $\RoC = \sqrt{g \alpha \Delta T / 4 \Omega^2 H}$, that strikingly appears in each regime by consideration of the local, convection-scale Rossby number $\Rol=U/(2\Omega \ell)$. Thus we show that $\RoC$ scales with the local Rossby number $\Rol$ in both the slowly rotating and the rapidly rotating regimes, explaining the ubiquity of $\RoC$ in rotating convection studies. We show in non-, slowly, and rapidly rotating systems that the convective heat transport, parameterized via $Pe_\ell$, scales with the total heat transport parameterized via the Nusselt number $Nu$. Within the rapidly-rotating limit, momentum transport arguments generate a scaling for the system-scale Rossby number, $Ro_H$, that, recast in terms of the total heat flux through the system, is shown to be synonymous with the classical flux-based `CIA' scaling, $Ro_{CIA}$. These, in turn, are then shown to asymptote to $Ro_H \sim Ro_{CIA} \sim \RoC^2$, demonstrating that these momentum transport scalings are identical in the limit of rapidly rotating turbulent heat transfer. 
    more » « less
  4. Abstract

    We analyze the directional dependence of the gravitational wave (GW) emission from 15 3D neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Using spin weighted spherical harmonics, we develop a new analytic technique to quantify the evolution of the distribution of GW emission over all angles. We construct a physics-informed toy model that can be used to approximate GW distributions for general ellipsoid-like systems, and use it to provide closed form expressions for the distribution of GWs for different CCSN phases. Using these toy models, we approximate the protoneutron star (PNS) dynamics during multiple CCSN stages and obtain similar GW distributions to simulation outputs. When considering all viewing angles, we apply this new technique to quantify the evolution of preferred directions of GW emission. For nonrotating cases, this dominant viewing angle drifts isotropically throughout the supernova, set by the dynamical timescale of the PNS. For rotating cases, during core bounce and the following tens of milliseconds, the strongest GW signal is observed along the equator. During the accretion phase, comparable—if not stronger—GW amplitudes are generated along the axis of rotation, which can be enhanced by the lowT/∣W∣ instability. We show two dominant factors influencing the directionality of GW emission are the degree of initial rotation and explosion morphology. Lastly, looking forward, we note the sensitive interplay between GW detector site and supernova orientation, along with its effect on detecting individual polarization modes.

     
    more » « less
  5. Dipolar spin ensembles with random spin positions attract much attention currently because they help to understand decoherence as it occurs in solid state quantum bits in contact with spin baths. Also, these ensembles are systems which may show many-body localization, at least in the sense of very slow spin dynamics. We present measurements of the autocorrelations of spins on diamond surfaces in a doubly-rotating frame which eliminates local disorder. Strikingly, the time scales in the longitudinal and the transversal channel differ by more than one order of magnitude which is a factor much greater than one would have expected from simulations of spins on lattices. A previously developed dynamic mean-field theory for spins (spinDMFT) fails to explain this phenomenon. Thus, we improve it by extending it to clusters (CspinDMFT). This theory does capture the striking mismatch up to two orders of magnitude for random ensembles. Without positional disorder, however, the mismatch is only moderate with a factor below 4. The pivotal role of positional disorder suggests that the strong mismatch is linked to precursors of many-body localization. 
    more » « less