Jacobian-free variational method for computing connecting orbits in nonlinear dynamical systems.

One approach for describing spatiotemporal chaos is to study the unstable invariant sets embedded in the chaotic attractor of the system. While equilibria, periodic orbits, and invariant tori can be computed using existing methods, the numerical identification of heteroclinic and homoclinic connections between them remains challenging. We propose a robust matrix-free variational method for computing connecting orbits between equilibrium solutions. Instead of a common shooting-based approach, we view the identification of a connecting orbit as a minimization problem in the space of smooth curves in the state space that connect the two equilibria. In this approach, the deviation of a connecting curve from an integral curve of the vector field is penalized by a non-negative cost function. Minimization of the cost function deforms a trial curve until, at a global minimum, a connecting orbit is obtained. The method has no limitation on the dimension of the unstable manifold at the origin equilibrium and does not suffer from exponential error amplification associated with time-marching a chaotic system. Owing to adjoint-based minimization techniques, no Jacobian matrices need to be constructed. Therefore, the memory requirement scales linearly with the size of the problem, allowing the method to be applied to high-dimensional dynamical systems. The robustness of the method is demonstrated for the one-dimensional Kuramoto-Sivashinsky equation.

Predicting chaotic statistics with unstable invariant tori.

It has recently been speculated that long-time average quantities of hyperchaotic dissipative systems may be approximated by weighted sums over unstable invariant tori embedded in the attractor, analogous to equivalent sums over periodic orbits, which are inspired by the rigorous periodic orbit theory and which have shown much promise in fluid dynamics. Using a new numerical method for converging unstable invariant two-tori in a chaotic partial differential equation (PDE), and exploiting symmetry breaking of relative periodic orbits to detect those tori, we identify many quasiperiodic, unstable, invariant two-torus solutions of a modified Kuramoto-Sivashinsky equation. The set of tori covers significant parts of the chaotic attractor and weighted averages of the properties of the tori-with weights computed based on their respective stability eigenvalues-approximate average quantities for the chaotic dynamics. These results are a step toward exploiting higher-dimensional invariant sets to describe general hyperchaotic systems, including dissipative spatiotemporally chaotic PDEs.

Invariant tori in dissipative hyperchaos.

One approach to understand the chaotic dynamics of nonlinear dissipative systems is the study of non-chaotic yet dynamically unstable invariant solutions embedded in the system's chaotic attractor. The significance of zero-dimensional unstable fixed points and one-dimensional unstable periodic orbits capturing time-periodic dynamics is widely accepted for high-dimensional chaotic systems, including fluid turbulence, while higher-dimensional invariant tori representing quasiperiodic dynamics have rarely been considered. We demonstrate that unstable 2-tori are generically embedded in the hyperchaotic attractor of a dissipative system of ordinary differential equations; tori can be numerically identified via bifurcations of unstable periodic orbits and their parameteric continuation and characterization of stability properties are feasible. As higher-dimensional tori are expected to be structurally unstable, 2-tori together with periodic orbits and equilibria form a complete set of relevant invariant solutions on which to base a dynamical description of chaos.

Nondestructive Prediction of the Buckling Load of Imperfect Shells.

From soda cans to space rockets, thin-walled cylindrical shells are abundant, offering exceptional load carrying capacity at relatively low weight. However, the actual load at which any shell buckles and collapses is very sensitive to imperceptible defects and cannot be predicted, which challenges the of such structures. Consequently, probabilistic descriptions in terms of empirical design rules are used and designing reliable structures requires the use of conservative strength estimates. We introduce a nonlinear description where finite-amplitude perturbations trigger buckling. Drawing from the analogy between imperfect shells which buckle and imperfect pipe flow which becomes turbulent, we experimentally show that lateral probing of cylindrical shells reveals their strength nondestructively. A new ridge-tracking method is applied to commercial cylinders with a hole showing that when the location where buckling nucleates is known we can accurately predict the buckling load of each individual shell, within ±5%. Our study provides a new promising framework to understand shell buckling, and more generally, imperfection-sensitive instabilities.

Stability Landscape of Shell Buckling.

We measure the response of cylindrical shells to poking and identify a stability landscape, which fully characterizes the stability of perfect shells and imperfect ones in the case where a single defect dominates. We show that the landscape of stability is independent of the loading protocol and the poker geometry. Our results suggest that the complex stability of shells reduces to a low dimensional description. Tracking ridges and valleys of this landscape defines a natural phase-space coordinates for describing the stability of shells.

Increasing lifetimes and the growing saddles of shear flow turbulence.

In linearly stable shear flows, turbulence spontaneously decays with a characteristic lifetime that varies with Reynolds number. The lifetime sharply increases with Reynolds number so that a possible divergence marking the transition to sustained turbulence at a critical point has been discussed. We present a mechanism by which the lifetimes increase: in the system's state space, turbulent motion is supported by a chaotic saddle. Inside this saddle a locally attracting periodic orbit is created and undergoes a traditional bifurcation sequence generating chaos. The formed new "turbulent bubble" is initially an attractor supporting persistent chaotic dynamics. Soon after its creation, it collides with its own boundary, by which it becomes leaky and dynamically connected with the surrounding structures. The complexity of the chaotic saddle that supports transient turbulence hence increases by incorporating the remnant of a new bubble. As a a result, the time it takes for a trajectory to leave the saddle and decay to the laminar state is increased. We demonstrate this phenomenon in plane Couette flow and show that characteristic lifetimes vary nonsmoothly and nonmonotonically with Reynolds number.

Finite lifetime of turbulence in shear flows.

Generally, the motion of fluids is smooth and laminar at low speeds but becomes highly disordered and turbulent as the velocity increases. The transition from laminar to turbulent flow can involve a sequence of instabilities in which the system realizes progressively more complicated states, or it can occur suddenly. Once the transition has taken place, it is generally assumed that, under steady conditions, the turbulent state will persist indefinitely. The flow of a fluid down a straight pipe provides a ubiquitous example of a shear flow undergoing a sudden transition from laminar to turbulent motion. Extensive calculations and experimental studies have shown that, at relatively low flow rates, turbulence in pipes is transient, and is characterized by an exponential distribution of lifetimes. They also suggest that for Reynolds numbers exceeding a critical value the lifetime diverges (that is, becomes infinitely large), marking a change from transient to persistent turbulence. Here we present experimental data and numerical calculations covering more than two decades of lifetimes, showing that the lifetime does not in fact diverge but rather increases exponentially with the Reynolds number. This implies that turbulence in pipes is only a transient event (contrary to the commonly accepted view), and that the turbulent and laminar states remain dynamically connected, suggesting avenues for turbulence control.

Constructing periodic orbits of high-dimensional chaotic systems by an adjoint-based variational method.

Chaotic dynamics in systems ranging from low-dimensional nonlinear differential equations to high-dimensional spatiotemporal systems including fluid turbulence is supported by nonchaotic, exactly recurring time-periodic solutions of the governing equations. These unstable periodic orbits capture key features of the turbulent dynamics and sufficiently large sets of orbits promise a framework to predict the statistics of the chaotic flow. Computing periodic orbits for high-dimensional spatiotemporally chaotic systems remains challenging as known methods either show poor convergence properties because they are based on time-marching of a chaotic system causing exponential error amplification, or they require constructing Jacobian matrices which is prohibitively expensive. We propose a new matrix-free method that is unaffected by exponential error amplification, is globally convergent, and can be applied to high-dimensional systems. The adjoint-based variational method constructs an initial value problem in the space of closed loops such that periodic orbits are attracting fixed points for the loop dynamics. We introduce the method for general autonomous systems. An implementation for the one-dimensional Kuramoto-Sivashinsky equation demonstrates the robust convergence of periodic orbits underlying spatiotemporal chaos. Convergence does not require accurate initial guesses and is independent of the period of the respective orbit.

Algorithm for a microfluidic assembly line.

Microfluidic technology has revolutionized the control of flows at small scales giving rise to new possibilities for assembling complex structures on the microscale. We analyze different possible algorithms for assembling arbitrary structures, and demonstrate that a sequential assembly algorithm can manufacture arbitrary 3D structures from identical constituents. We illustrate the algorithm by showing that a modified Hele-Shaw cell with 7 controlled flow rates can be designed to construct the entire English alphabet from particles that irreversibly stick to each other.

Superspreading events suggest aerosol transmission of SARS-CoV-2 by accumulation in enclosed spaces.

Viral transmission pathways have profound implications for public safety; it is thus imperative to establish a complete understanding of viable infectious avenues. Mounting evidence suggests SARS-CoV-2 can be transmitted via the air; however, this has not yet been demonstrated. Here we quantitatively analyze virion accumulation by accounting for aerosolized virion emission and destabilization. Reported superspreading events analyzed within this framework point towards aerosol mediated transmission of SARS-CoV-2. Virion exposure calculated for these events is found to trace out a single value, suggesting a universal minimum infective dose (MID) via aerosol that is comparable to the MIDs measured for other respiratory viruses; thus, the consistent infectious exposure levels and their commensurability to known aerosol-MIDs establishes the plausibility of aerosol transmission of SARS-CoV-2. Using filtration at a rate exceeding the destabilization rate of aerosolized SARS-CoV-2 can reduce exposure below this infective dose.

Snakes and ladders: localized solutions of plane Couette flow.

We demonstrate the existence of a large number of exact solutions of plane Couette flow, which share the topology of known periodic solutions but are localized in one spatial dimension. Solutions of different size are organized in a snakes-and-ladders structure strikingly similar to that observed for simpler pattern-forming partial differential equations. These new solutions are a step towards extending the dynamical systems view of transitional turbulence to spatially extended flows.

Buckling instabilities and spatio-temporal dynamics of active elastic filaments.

Biological filaments driven by molecular motors tend to experience tangential propulsive forces also known as active follower forces. When such a filament encounters an obstacle, it deforms, which reorients its follower forces and alters its entire motion. If the filament pushes a cargo, the friction on the cargo can be enough to deform the filament, thus affecting the transport properties of the cargo. Motivated by cytoskeletal filament motility assays, we study the dynamic buckling instabilities of a two-dimensional slender elastic filament driven through a dissipative medium by tangential propulsive forces in the presence of obstacles or cargo. We observe two distinct instabilities. When the filament's head is pinned or experiences significant translational but little rotational drag from its cargo, it buckles into a steadily rotating coiled state. When it is clamped or experiences both significant translational and rotational drag from its cargo, it buckles into a periodically beating, overall translating state. Using minimal analytically tractable models, linear stability theory and fully nonlinear computations, we study the onset of each buckling instability, characterize each buckled state, and map out the phase diagram of the system. Finally, we use particle-based Brownian dynamics simulations to show our main results are robust to moderate noise and steric repulsion. Overall, our results provide a unified framework to understand the dynamics of tangentially propelled filaments and filament-cargo assemblies.

Exact invariant solution reveals the origin of self-organized oblique turbulent-laminar stripes.

Wall-bounded shear flows transitioning to turbulence may self-organize into alternating turbulent and laminar regions forming a stripe pattern with non-trivial oblique orientation. Different experiments and flow simulations identify oblique stripe patterns as the preferred solution of the well-known Navier-Stokes equations, but the origin of stripes and their oblique orientation remains unexplained. In concluding his lectures, Feynman highlights the unexplained stripe pattern hidden in the solution space of the Navier-Stokes equations as an example demonstrating the need for improved theoretical tools to analyze the fluid flow equations. Here we exploit dynamical systems methods and demonstrate the existence of an exact equilibrium solution of the fully nonlinear 3D Navier-Stokes equations that resembles oblique stripe patterns in plane Couette flow. The stripe equilibrium emerges from the well-studied Nagata equilibrium and exists only for a limited range of pattern angles. This suggests a mechanism selecting the non-trivial oblique orientation angle of turbulent-laminar stripes.

Origin of localized snakes-and-ladders solutions of plane Couette flow.

Spatially localized invariant solutions of plane Couette flow are organized in a snakes-and-ladders structure strikingly similar to that observed for simpler pattern-forming partial differential equations [Schneider, Gibson, and Burke, Phys. Rev. Lett. 104, 104501 (2010)PRLTAO0031-900710.1103/PhysRevLett.104.104501]. We demonstrate the mechanism by which these snaking solutions originate from well-known periodic states of the Taylor-Couette system. They are formed by a localized slug of wavy-vortex flow that emerges from a background of Taylor vortices via a modulational sideband instability. This mechanism suggests a close connection between pattern-formation theory and Navier-Stokes flow.

Rational design of a high-throughput droplet sorter.

The high-throughput selection of individual droplets is an essential function in droplet-based microfluidics. Fluorescence-activated droplet sorting is achieved using electric fields triggered at rates up to 30 kHz, providing the ultra-high throughput relevant in applications where large libraries of compounds or cells must be analyzed. To achieve such sorting frequencies, electrodes have to create an electric field distribution that generates maximal actuating forces on the droplet while limiting the induced droplet deformation and avoid disintegration. We propose a metric characterizing the performance of an electrode design relative to the theoretical optimum and analyze existing devices using full 3D simulations of the electric fields. By combining parameter optimization with numerical simulation we derive rational design guidelines and propose optimized electrode configurations. When tested experimentally, the optimized design show significantly better performance than the standard designs.

Lifetime statistics in transitional pipe flow.

Several experimental and numerical studies have shown that turbulent motions in circular pipe flow near transitional Reynolds numbers may not persist forever, but may decay. We study the properties of these decaying states within direct numerical simulations for Reynolds numbers up to 2200 and in pipes with lengths equal to 5, 9, and 15 times the diameter. We show that the choice of the ensemble of initial conditions affects the short time parts of lifetime distributions, but does not change the characteristic decay rate for long times. Comparing lifetimes for pipes of different length we notice a linear increase in the characteristic lifetime with length, which reproduces the experimental results when extrapolated to 30 diameters, the length of an equilibrium turbulent puff at these Reynolds numbers.

Laminar-turbulent boundary in plane Couette flow.

We apply the iterated edge-state tracking algorithm to study the boundary between laminar and turbulent dynamics in plane Couette flow at Re=400. Perturbations that are not strong enough to become fully turbulent or weak enough to relaminarize tend toward a hyperbolic coherent structure in state space, termed the edge state, which seems to be unique up to obvious continuous shift symmetries. The results reported here show that in cases where a fixed point has only one unstable direction, such as for the lower-branch solution in plane Couette flow, the iterated edge tracking algorithm converges to this state. They also show that the choice of initial state is not critical and that essentially arbitrary initial conditions can be used to find the edge state.

High-throughput multiplexed fluorescence-activated droplet sorting.

Fluorescence-activated droplet sorting (FADS) is one of the most important features provided by droplet-based microfluidics. However, to date, it does not allow to compete with the high-throughput multiplexed sorting capabilities offered by flow cytometery. Here, we demonstrate the use of a dielectrophoretic-based FADS, allowing to sort up to five different droplet populations simultaneously. Our system provides means to select droplets of different phenotypes in a single experimental run to separate initially heterogeneous populations. Our experimental results are rationalized with the help of a numerical model of the actuation of droplets in electric fields providing guidelines for the prediction of sorting designs for upscaled or downscaled microsystems.

Statistical analysis of coherent structures in transitional pipe flow.

Numerical and experimental studies of transitional pipe flow have shown the prevalence of coherent flow structures that are dominated by downstream vortices. They attract special attention because they contribute predominantly to the increase of the Reynolds stresses in turbulent flow. In the present study we introduce a convenient detector for these coherent states, calculate the fraction of time the structures appear in the flow, and present a Markov model for the transition between the structures. The fraction of states that show vortical structures exceeds 24% for a Reynolds number of about Re=2200 , and it decreases to about 20% for Re=2500 . The Markov model for the transition between these states is in good agreement with the observed fraction of states, and in reasonable agreement with the prediction for their persistence. It provides insight into dominant qualitative changes of the flow when increasing the Reynolds number.

Taking the Perfect Selfie: Investigating the Impact of Perspective on the Perception of Higher Cognitive Variables.

Taking selfies is now becoming a standard human habit. However, as a social phenomenon, research is still in the fledgling stage and the scientific framework is sparse. Selfies allow us to share social information with others in a compact format. Furthermore, we are able to control important photographic and compositional aspects, such as perspective, which have a strong impact on the assessment of a face (e.g., demonstrated by the height-weight illusion, effects of gaze direction, faceism-index). In Study 1, we focused on the impact of perspective (left/right hemiface, above/below vs. frontal presentation) on higher cognitive variables and let 172 participants rate the perceived attractiveness, helpfulness, sympathy, dominance, distinctiveness, and intelligence, plus important information on health issues (e.g., body weight), on the basis of 14 3D faces. We could show that lateral snapshots yielded higher ratings for attractiveness compared to the classical frontal view. However, this effect was more pronounced for left hemifaces and especially female faces. Compared to the frontal condition, 30° right hemifaces were rated as more helpful, but only for female faces while faces viewed from above were perceived as significant less helpful. Direct comparison between left vs. right hemifaces revealed no effect. Relating to sympathy, we only found a significant effect for 30° right male hemifaces, but only in comparison to the frontal condition. Furthermore, female 30° right hemifaces were perceived as more intelligent. Relating to body weight, we replicated the so-called "height-weight illusion." Other variables remained unaffected. In Study 2, we investigated the impact of a typical selfie-style condition by presenting the respective faces from a lateral (left/right) and tilted (lower/higher) vantage point. Most importantly, depending on what persons wish to express with a selfie, a systematic change of perspective can strongly optimize their message; e.g., increasing their attractiveness by shooting from above left, and in contrast, decreasing their expressed helpfulness by shooting from below. We could further extent past findings relating to the height-weight illusion and showed that an additional rotation of the camera positively affected the perception of body weight (lower body weight). We discuss potential explanations for perspective-related effects, especially gender-related ones.