BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

George Haller

George Haller appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.DSphysics.flu-dynnlin.CDeess.SYmath-phmath.MPphysics.ao-phSystems and Controlphysics.geo-phMachine Learningmath.OCnlin.PSRoboticsastro-ph.EPComputational Engineering, Finance, and ScienceGraphicsmath.APnlin.SIphysics.plasm-ph

Catalog footprint

What is connected

41works
19topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2026arXiv

Data-Driven Modeling of Global Bifurcations and Chaos in a Mechanical System under Delayed and Quantized Control

We illustrate how the recent theory of Spectral Submanifolds (SSM) can capture global bifurcations and complex dynamics in mechanical systems even under delay and spatial discretization. Specifically, we build a parameter-dependent SSM-reduced model that predicts global heteroclinic and local bifurcations in a Furuta pendulum under control with delay, and verify these predictions numerically. Under additional spatial discretization of the digital controller, we also obtain an SSM-reduced model that correctly reproduces a numerically and experimentally observed microchaotic attractor in the system.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Nonlinear Spectral Modeling and Control of Soft-Robotic Muscles from Data

Artificial muscles are essential for compliant musculoskeletal robotics but complicate control due to nonlinear multiphysics dynamics. Hydraulically amplified electrostatic (HASEL) actuators, a class of soft artificial muscles, offer high performance but exhibit memory effects and hysteresis. Here we present a data-driven reduction and control strategy grounded in spectral submanifold (SSM) theory. In the adiabatic regime, where inputs vary slowly relative to intrinsic transients, trajectories rapidly converge to a low-dimensional slow manifold. We learn an explicit input-to-output map on this manifold from forced-response trajectories alone, avoiding decay experiments that can trigger hysteresis. We deploy the SSM-based model for real-time control of an antagonistic HASEL-clutch joint. This approach yields a substantial reduction in tracking error compared to feedback-only and feedforward-only baselines under identical settings. This record-and-control workflow enables rapid characterization and high-performance control of soft muscles and muscle-driven joints without detailed physics-based modeling.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Model reduction for nonlinearizable dynamics via delay-embedded spectral submanifolds

Delay embedding is a commonly employed technique in a wide range of data-driven model reduction methods for dynamical systems, including the Dynamic mode decomposition (DMD), the Hankel alternative view of the Koopman decomposition (HAVOK), nearest-neighbor predictions and the reduction to spectral submanifolds (SSMs). In developing these applications, multiple authors have observed that delay embedding appears to separate the data into modes, whose orientations depend only on the spectrum of the sampled system. In this work, we make this observation precise by proving that the eigenvectors of the delay-embedded linearized system at a fixed point are determined solely by the corresponding eigenvalues, even for multi-dimensional observables. This implies that the tangent space of a delay-embedded invariant manifold can be predicted a priori using an estimate of the eigenvalues. We apply our results to three datasets to identify multimodal SSMs and analyse their nonlinear modal interactions. While SSMs are the focus of our study, these results generalize to any delay-embedded invariant manifold tangent to a set of eigenvectors at a fixed point. Therefore, we expect this theory to be applicable to a number of data-driven model reduction methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Data-Driven Modeling and Prediction of Non-Linearizable Dynamics via Spectral Submanifolds

We develop a methodology to construct low-dimensional predictive models from data sets representing essentially nonlinear (or non-linearizable) dynamical systems with a hyperbolic linear part that are subject to external forcing with finitely many frequencies. Our data-driven, sparse, nonlinear models are obtained as extended normal forms of the reduced dynamics on low-dimensional, attracting spectral submanifolds (SSMs) of the dynamical system. We illustrate the power of data-driven SSM reduction on high-dimensional numerical data sets and experimental measurements involving beam oscillations, vortex shedding and sloshing in a water tank. We find that SSM reduction trained on unforced data also predicts nonlinear response accurately under additional external forcing.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dynamics-based machine learning of transitions in Couette flow

We derive low-dimensional, data-driven models for transitions among exact coherent states (ECSs) in one of the most studied canonical shear flows, the plane Couette flow. These one- or two-dimensional nonlinear models represent the leading-order reduced dynamics on attracting spectral submanifolds (SSMs), which we construct using the recently developed SSMLearn algorithm from a small number of simulated transitions. We find that the energy input and output rates provide efficient parametrizations for the most important SSMs. By restricting the dynamics to these SSMs, we obtain reduced-order models that also reliably predict nearby, off-SSM transitions that were not used in their training.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fast data-driven model reduction for nonlinear dynamical systems

We present a fast method for nonlinear data-driven model reduction of dynamical systems onto their slowest nonresonant spectral submanifolds (SSMs). We use observed data to locate a low-dimensional, attracting slow SSM and compute a maximally sparse approximation to the reduced dynamics on it. The recently released SSMLearn algorithm uses implicit optimization to fit a spectral submanifold to data and reduce the dynamics to the normal form. Here, we present two simplified algorithms, which reformulate manifold fitting and normal form computation as explicit problems under certain assumptions. We show on both numerical and experimental datasets that these algorithms yield accurate and sparse rigorous models for essentially nonlinear (or non-linearizable) dynamics. The new algorithms are significantly simplified and provide a speedup of several orders of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nonlinear analysis of forced mechanical systems with internal resonance using spectral submanifolds, Part I: Periodic response and forced response curve

We show how spectral submanifold theory can be used to construct reduced-order models for harmonically excited mechanical systems with internal resonances. Efficient calculations of periodic and quasi-periodic responses with the reduced-order models are discussed in this paper and its companion, Part II, respectively. The dimension of a reduced-order model is determined by the number of modes involved in the internal resonance, independently of the dimension of the full system. The periodic responses of the full system are obtained as equilibria of the reduced-order model on spectral submanifolds. The forced response curve of periodic orbits then becomes a manifold of equilibria, which can be easily extracted using parameter continuation. To demonstrate the effectiveness and efficiency of the reduction, we compute the forced response curves of several high-dimensional nonlinear mechanical systems, including the finite-element models of a von Kármán beam and a plate.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nonlinear analysis of forced mechanical systems with internal resonance using spectral submanifolds, Part II: Bifurcation and quasi-periodic response

In Part I of this paper, we have used spectral submanifold (SSM) theory to construct reduced-order models for harmonically excited mechanical systems with internal resonances. In that setting, extracting forced response curves formed by periodic orbits of the full system was reduced to locating the solution branches of equilibria of the corresponding reduced-order model. Here we use bifurcations of the equilibria of the reduced-order model to predict bifurcations of the periodic response of the full system. Specifically, we identify Hopf bifurcations of equilibria and limit cycles in reduced models on SSMs to predict the existence of two-dimensional and three-dimensional quasi-periodic attractors and repellers in periodically forced mechanical systems of arbitrary dimension. We illustrate the accuracy and efficiency of these computations on finite-element models of beams and plates.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Objective Momentum Barriers in Wall Turbulence

We use the recent frame-indifferent theory of diffusive momentum transport to identify internal barriers in wall-bounded turbulence. Formed by the invariant manifolds of the Laplacian of the velocity field, the barriers block the viscous part of the instantaneous momentum flux in the flow. We employ the level sets of single-trajectory Lagrangian diagnostic tools, the trajectory rotation average and trajectory stretching exponent, to approximate both vortical and internal wall-parallel momentum transport barrier (MTB) interfaces. These interfaces provide frame-indifferent alternatives to classic velocity-gradient-based vortices and high-shear boundaries between uniform momentum zones (UMZs). Indeed, we find that these elliptic manifold approximations and MTBs outperform standard vortices and UMZ interfaces in blocking diffusive momentum transport, suggesting our momentum barriers are physical features that may be the cause of coherence signatures in statistical and non-objective diagnostics. We also introduce normalized trajectory metrics that provide unprecedented visualizations of objective coherent structures by avoiding strong turbulence biases.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quasi-Objective Coherent Structure Diagnostics from Single Trajectories

We derive measures of local material stretching and rotation that are computable from individual trajectories without reliance on other trajectories or on an underlying velocity field. Both measures are quasi-objective: they approximate objective (i.e., observer-independent) coherence diagnostics in frames satisfying a certain condition. This condition requires the trajectory accelerations to dominate the angular acceleration induced by the spatial mean vorticity. We illustrate in examples how quasi-objective coherence diagnostics highlight elliptic and hyperbolic Lagrangian coherent structures even from very sparse trajectory data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Objective Deformation Component of a Velocity Field

For an arbitrary velocity field $\mathbf{v}$ defined on a finite, fixed spatial domain, we find the closest rigid-body velocity field $\mathbf{v}_{RB}$ to $\mathbf{v}$ in the $L^2$ norm. The resulting deformation velocity component, $\mathbf{v}_{d}=\mathbf{v-\mathbf{v}}_{RB}$, turns out to be frame-indifferent and physically observable. Specifically, if $\mathbf{Q}_{\text{RB}}(t)$ is the rotation tensor describing the motion of the closest rigid body frame, then $\mathbf{v}$ is seen as $\mathbf{Q}_{\text{RB}}^{T}\mathbf{v}_{d}$ by an observer in that frame. As a consequence, the momentum, energy, vorticity, enstrophy, and helicity of the flow all become frame-indifferent when computed from the deformation velocity component $\mathbf{v}_{d}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Using Spectral Submanifolds for Nonlinear Periodic Control

Very high dimensional nonlinear systems arise in many engineering problems due to semi-discretization of the governing partial differential equations, e.g. through finite element methods. The complexity of these systems present computational challenges for direct application to automatic control. While model reduction has seen ubiquitous applications in control, the use of nonlinear model reduction methods in this setting remains difficult. The problem lies in preserving the structure of the nonlinear dynamics in the reduced order model for high-fidelity control. In this work, we leverage recent advances in Spectral Submanifold (SSM) theory to enable model reduction under well-defined assumptions for the purpose of efficiently synthesizing feedback controllers.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Data-driven Nonlinear Model Reduction to Spectral Submanifolds in Mechanical Systems

While data-driven model reduction techniques are well-established for linearizable mechanical systems, general approaches to reducing non-linearizable systems with multiple coexisting steady states have been unavailable. In this paper, we review such a data-driven nonlinear model reduction methodology based on spectral submanifolds (SSMs). As input, this approach takes observations of unforced nonlinear oscillations to construct normal forms of the dynamics reduced to very low dimensional invariant manifolds. These normal forms capture amplitude-dependent properties and are accurate enough to provide predictions for non-linearizable system response under the additions of external forcing. We illustrate these results on examples from structural vibrations, featuring both synthetic and experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Explicit Unsteady Navier-Stokes Solutions and their Analysis via Local Vortex Criteria

We construct a class of spatially polynomial velocity fields that are exact solutions of the planar unsteady Navier-Stokes equation. These solutions can be used as simple benchmarks for testing numerical methods or verifying the feasibility of flow-feature identification principles. We use examples from the constructed solution family to illustrate deficiencies of streamlines-based feature detection and of the Okubo-Weiss criterion, which is the common two-dimensional version of the broadly used Q-, Delta-, Lambda-2- and Lambda-Ci-criteria for vortex-detection. Our planar polynomial solutions also extend directly to explicit, three-dimensional unsteady Navier-Stokes solutions with a symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How do Conservative Backbone Curves Perturb into Forced Responses? A Melnikov Function Analysis

Weakly damped mechanical systems under small periodic forcing tend to exhibit periodic response in a close vicinity of certain periodic orbits of their conservative limit. Specifically, amplitude frequency plots for the conservative limit have often been noted, both numerically and experimentally, to serve as backbone curves for the near resonance peaks of the forced response. In other cases, such a relationship between the unforced and forced response was not observed. Here we provide a systematic mathematical analysis that predicts which members of conservative periodic orbit families will serve as backbone curves for the forced-damped response. We also obtain mathematical conditions under which approximate numerical and experimental approaches, such as energy balance and force appropriation, are justifiable. Finally, we derive analytic criteria for the birth of isolated response branches (isolas) whose identification is otherwise challenging from numerical continuation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Universal Upper Estimate for Prediction Errors under Moderate Model Uncertainty

We derive universal upper estimates for model-prediction error under moderate but otherwise unknown model uncertainty. Our estimates give upper bounds on the leading order trajectory-uncertainty arising along model trajectories, solely as functions of the invariants of the known Cauchy-Green strain tensor of the model. Our bounds turn out to be optimal, which means that they cannot be improved for general systems. The quantity relating the leading-order trajectory-uncertainty to the model uncertainty is the Model Sensitivity, which we find to be a useful tool for a quick global assessment of the impact of modeling uncertainties in various domains of the phase space. Examining the expectation that Finite-Time Lyapunov Exponents capture sensitivity to modeling errors, we show that this does not generally follow. However, we find that certain important features of the FTLE persist in the MS field.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Barriers to the Transport of Diffusive Scalars in Compressible Flows

Our recent work identifies material surfaces in incompressible flows that extremize the transport of an arbitrary, weakly diffusive scalar field relative to neighboring surfaces. Such barriers and enhancers of transport can be located directly from the deterministic component of the velocity field without diffusive or stochastic simulations. Here we extend these results to compressible flows and to diffusive concentration fields affected by sources or sinks, as well as by spontaneous decay. We construct diffusive transport extremizers with and without constraining them on a specific initial concentration distribution. For two-dimensional flows, we obtain explicit differential equations and a diagnostic scalar field that identify the most observable extremizers with pointwise uniform transport density. We illustrate our results by uncovering diffusion barriers and enhancers in analytic, numerical, and observational velocity fields.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Explicit Third-Order Model Reduction Formulas for General Nonlinear Mechanical Systems

For general nonlinear mechanical systems, we derive closed-form, reduced-order models up to cubic order based on rigorous invariant manifold results. For conservative systems, the reduction is based on Lyapunov Subcenter Manifold (LSM) theory, whereas for damped-forced systems, we use Spectral Submanifold (SSM) theory. To evaluate our explicit formulas for the reduced model, no coordinate changes are required beyond an initial linear one. The reduced-order models we derive are simple and depend only on physical and modal parameters, allowing us to extract fundamental characteristics, such as backbone curves and forced-response curves, of multi-degree-of-freedom mechanical systems. To numerically verify the accuracy of the reduced models, we test the reduction formulas on several mechanical systems, including a higher-dimensional nonlinear Timoshenko beam.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Search and rescue at sea aided by hidden flow structures

Every year hundreds of people die at sea because of vessel and airplane accidents. A key challenge in reducing the number of these fatalities is to make Search and Rescue (SAR) algorithms more efficient. Here we address this challenge by uncovering hidden TRansient Attracting Profiles (TRAPs) in ocean-surface velocity data. Computable from a single velocity-field snapshot, TRAPs act as short-term attractors for all floating objects. In three different ocean field experiments, we show that TRAPs computed from measured as well as modelled velocities attract deployed drifters and manikins emulating people fallen in the water. TRAPs, which remain hidden to prior flow diagnostics, thus provide critical information for hazard responses, such as SAR and oil spill containment, and hence have the potential to save lives and limit environmental disasters.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Vortex boundaries as barriers to diffusive vorticity transport in two-dimensional flows

We put forward the idea of defining vortex boundaries in planar flows as closed material barriers to the diffusive transport of vorticity. Such diffusive vortex boundaries minimize the leakage of vorticity from the fluid mass they enclose when compared to other nearby material curves. Building on recent results on passive diffusion barriers, we develop an algorithm for the automated identification of such structures from general, two-dimensional unsteady flow data. As examples, we identify vortex boundaries as vorticity diffusion barriers in two flows: an explicitly known laminar flow and a numerically generated turbulent Navier--Stokes flow.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Spectral Clustering Approach to Lagrangian Vortex Detection

One of the ubiquitous features of real-life turbulent flows is the existence and persistence of coherent vortices. Here we show that such coherent vortices can be extracted as clusters of Lagrangian trajectories. We carry out the clustering on a weighted graph, with the weights measuring pairwise distances of fluid trajectories in the extended phase space of positions and time. We then extract coherent vortices from the graph using tools from spectral graph theory. Our method locates all coherent vortices in the flow simultaneously, thereby showing high potential for automated vortex tracking. We illustrate the performance of this technique by identifying coherent Lagrangian vortices in several two- and three-dimensional flows.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An Autonomous Dynamical System Captures all LCSs in Three-Dimensional Unsteady Flows

Lagrangian coherent structures (LCSs) are material surfaces that shape finite-time tracer patterns in flows with arbitrary time dependence. Depending on their deformation properties, elliptic and hyperbolic LCSs have been identified from different variational principles, solving different equations. Here we observe that, in three dimensions, initial positions of all variational LCSs are invariant manifolds of the same autonomous dynamical system, generated by the intermediate eigenvector field, $ξ_{2}(x_{0})$, of the Cauchy-Green strain tensor. This $ξ_{2}$-system allows for the detection of LCSs in any unsteady flow by classic methods, such as Poincaré maps, developed for autonomous dynamical systems. As examples, we consider both steady and time-aperiodic flows, and use their dual $ξ_{2}$-system to uncover both hyperbolic and elliptic LCSs from a single computation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Defining Coherent Vortices Objectively from the Vorticity

Rotationally coherent Lagrangian vortices are formed by tubes of deforming fluid elements that complete equal bulk material rotation relative to the mean rotation of the deforming fluid volume. We show that initial positions of such tubes coincide with tubular level surfaces of the Lagrangian-Averaged Vorticity Deviation (LAVD), the trajectory integral of the normed difference of the vorticity from its spatial mean. LAVD-based vortices are objective, i.e., remain unchanged under time-dependent rotations and translations of the coordinate frame. In the limit of vanishing Rossby numbers in geostrophic flows, cyclonic LAVD vortex centers are precisely the observed attractors for light particles. A similar result holds for heavy particles in anticyclonic LAVD vortices. We also establish a relationship between rotationally coherent Lagrangian vortices and their instantaneous Eulerian counterparts. The latter are formed by tubular surfaces of equal material rotation rate, objectively measured by the Instantaneous Vorticity Deviation (IVD). We illustrate the use of the LAVD and the IVD to detect rotationally coherent Lagrangian and Eulerian vortices objectively in several two- and three-dimensional flows.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Exact Model Reduction by a Slow-Fast Decomposition of Nonlinear Mechanical Systems

We derive conditions under which a general nonlinear mechanical system can be exactly reduced to a lower-dimensional model that involves only the most flexible degrees of freedom. This Slow-Fast Decomposition (SFD) enslaves exponentially fast the stiff degrees of freedom to the flexible ones as all oscillations converge to the reduced model defined on a slow manifold. We obtain an expression for the domain boundary beyond which the reduced model ceases to be relevant due to a generic loss of stability of the slow manifold. We also find that near equilibria, the SFD gives a mathematical justification for two modal-reduction methods used in structural dynamics: static condensation and modal derivatives. These formal reduction procedures, however, are also found to return incorrect results when the SFD conditions do not hold. We illustrate all these results on mechanical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Level set formulation of two-dimensional Lagrangian vortex detection methods

We propose here the use of the variational level set methodology to capture Lagrangian vortex boundaries in 2D unsteady velocity fields. This method reformulates earlier approaches that seek material vortex boundaries as extremum solutions of variational problems. We demonstrate the performance of this technique for two different variational formulations built upon different notions of coherence. The first formulation uses an energy functional that penalizes the deviation of a closed material line from piecewise uniform stretching [Haller and Beron-Vera, J. Fluid Mech. 731, R4 (2013)]. The second energy function is derived for a graph-based approach to vortex boundary detection [Hadjighasem et al., Phys. Rev. E 93, 063107 (2016)]. Our level-set formulation captures an a priori unknown number of vortices simultaneously at relatively low computational cost. We illustrate the approach by identifying vortices from different coherence principles in several examples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonlinear normal modes and spectral submanifolds: Existence, uniqueness and use in model reduction

We propose a unified approach to nonlinear modal analysis in dissipative oscillatory systems. This approach eliminates conflicting definitions, covers both autonomous and time-dependent systems, and provides exact mathematical existence, uniqueness and robustness results. In this setting, a nonlinear normal mode (NNM) is a set filled with small-amplitude recurrent motions: a fixed point, a periodic orbit or the closure of a quasiperiodic orbit. In contrast, a spectral submanifold (SSM) is an invariant manifold asymptotic to a NNM, serving as the smoothest nonlinear continuation of a spectral subspace of the linearized system along the NNM. The existence and uniqueness of SSMs turns out to depend on a spectral quotient computed from the real part of the spectrum of the linearized system. This quotient may well be large even for small dissipation, thus the inclusion of damping is essential for firm conclusions about NNMs, SSMs and the reduced-order models they yield.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamic Rotation and Stretch Tensors from a Dynamic Polar Decomposition

The local rigid-body component of continuum deformation is typically characterized by the rotation tensor, obtained from the polar decomposition of the deformation gradient. Beyond its well-known merits, the polar rotation tensor also has a lesser known dynamical inconsistency: it does not satisfy the fundamental superposition principle of rigid-body rotations over adjacent time intervals. As a consequence, the polar rotation diverts from the observed mean material rotation of fibers in fluids, and introduces a purely kinematic memory effect into computed material rotation. Here we derive a generalized polar decomposition for linear processes that yields a unique, dynamically consistent rotation component, the dynamic rotation tensor, for the deformation gradient. The left dynamic stretch tensor is objective, and shares the principal strain values and axes with its classic polar counterpart. Unlike its classic polar counterpart, however, the dynamic stretch tensor evolves in time without spin. The dynamic rotation tensor further decomposes into a spatially constant mean rotation tensor and a dynamically consistent relative rotation tensor that is objective for planar deformations. We also obtain simple expressions for dynamic analogues of Cauchy's mean rotation angle that characterize a deforming body objectively.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Explicit form of spatially linear Navier-Stokes velocity fields

We show that a smooth linear unsteady velocity field $u(x,t)=A(t)x+f(t)$ solves the incompressible Navier--Stokes equation if and only if the matrix $A(t)$ has zero trace, and $\dot{A}(t)+A^{2}(t)$ is symmetric. In two dimensions, these constraints imply that $A(t)$ is the sum of an arbitrary time-dependent traceless symmetric matrix and an arbitrary constant skew-symmetric matrix. One can, therefore, verify by inspection if an unsteady spatially linear vector field is a Navier--Stokes solution. In three dimensions, we obtain a simple ordinary differential equation that $A(t)$ must solve. Our formulas enable the construction of simple yet unsteady and dynamically consistent flows for testing numerical schemes and verifying coherent structure criteria.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Geodesic Transport Barriers in Jupiter's Atmosphere: A Video-Based Analysis

Jupiter's zonal jets and Great Red Spot are well known from still images. Yet the planet's atmosphere is highly unsteady, which suggests that the actual material transport barriers delineating its main features should be time-dependent. Rare video footages of Jupiter's clouds provide an opportunity to verify this expectation from optically reconstructed velocity fields. Available videos, however, provide short-time and temporally aperiodic velocity fields that defy classical dynamical systems analyses focused on asymptotic features. To this end, we use here the recent theory of geodesic transport barriers to uncover finite-time mixing barriers in the wind field extracted from a video captured by NASA's Cassini space mission. More broadly, the approach described here provides a systematic and frame-invariant way to extract dynamic coherent structures from time-resolved remote observations of unsteady continua.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Objective Eulerian Coherent Structures

We define objective Eulerian Coherent Structures (OECSs) in two-dimensional, non-autonomous dynamical systems as instantaneously most influential material curves. Specifically, OECSs are stationary curves of the averaged instantaneous material stretching-rate or material shearing-rate functionals. From these objective (frame-invariant) variational principles, we obtain explicit differential equations for hyperbolic, elliptic and parabolic OECSs. As illustration, we compute OECSs in an unsteady ocean velocity data set. In comparison to structures suggested by other common Eulerian diagnostic tools, we find OECSs to be the correct short-term cores of observed trajectory deformation patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Polar rotation angle identifies elliptic islands in unsteady dynamical systems

We propose rotation inferred from the polar decomposition of the flow gradient as a diagnostic for elliptic (or vortex-type) invariant regions in non-autonomous dynamical systems. We consider here two- and three-dimensional systems, in which polar rotation can be characterized by a single angle. For this polar rotation angle (PRA), we derive explicit formulas using the singular values and vectors of the flow gradient. We find that closed level sets of the PRA reveal elliptic islands in great detail, and singular level sets of the PRA uncover centers of such islands. Both features turn out to be objective (frame-invariant) for two-dimensional systems. We illustrate the diagnostic power of PRA for elliptic structures on several examples.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Asymptotic dynamics of inertial particles with memory

Recent experimental and numerical observations have shown the significance of the Basset--Boussinesq memory term on the dynamics of small spherical rigid particles (or inertial particles) suspended in an ambient fluid flow. These observations suggest an algebraic decay to an asymptotic state, as opposed to the exponential convergence in the absence of the memory term. Here, we prove that the observed algebraic decay is a universal property of the Maxey--Riley equation. Specifically, the particle velocity decays algebraically in time to a limit that is $\mathcal O(ε)$-close to the fluid velocity, where $0<ε\ll 1$ is proportional to the square of the ratio of the particle radius to the fluid characteristic length-scale. These results follows from a sharp analytic upper bound that we derive for the particle velocity. For completeness, we also present a first proof of existence and uniqueness of global solutions to the Maxey--Riley equation, a nonlinear system of fractional-order differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Automated detection of coherent Lagrangian vortices in two-dimensional unsteady flows

Coherent boundaries of Lagrangian vortices in fluid flows have recently been identified as closed orbits of line fields associated with the Cauchy-Green strain tensor. Here we develop a fully automated procedure for the detection of such closed orbits in large-scale velocity data sets. We illustrate the power of our method on ocean surface velocities derived from satellite altimetry.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Do Finite-Size Lyapunov Exponents Detect Coherent Structures?

Ridges of the Finite-Size Lyapunov Exponent (FSLE) field have been used as indicators of hyperbolic Lagrangian Coherent Structures (LCSs). A rigorous mathematical link between the FSLE and LCSs, however, has been missing. Here we prove that an FSLE ridge satisfying certain conditions does signal a nearby ridge of some Finite-Time Lyapunov Exponent (FTLE) field, which in turn indicates a hyperbolic LCS under further conditions. Other FSLE ridges violating our conditions, however, are seen to be false positives for LCSs. We also find further limitations of the FSLE in Lagrangian coherence detection, including ill-posedness, artificial jump-discontinuities, and sensitivity with respect to the computational time step.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

How coherent are the vortices of two-dimensional turbulence?

We use recent developments in the theory of finite-time dynamical systems to objectively locate the material boundaries of coherent vortices in two-dimensional Navier--Stokes turbulence. We show that these boundaries are optimal in the sense that any closed curve in their exterior will lose coherence under material advection. Through a detailed comparison, we find that other available Eulerian and Lagrangian techniques significantly underestimate the size of each coherent vortex.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Lagrangian Coherent Structures from Video Streams of Jupiter

Jupiter's fast rotation - one rotation over 10 hours - creates strong jet streams, smearing its clouds into linear bands of dark and light zonal belts that circle the planet on lines of almost constant latitude. Such a high degree of axisymmetry is absent in our own atmosphere. Moreover, Jupiter has the largest and longest-living known atmospheric vortex, the Great Red Spot (GRS). Such vortices abound in nature, but GRS's size, long-term persistence, and temporal longitudinal oscillations make it unique. Here, we uncover, for the first time, unsteady material structures that form the cores of zonal jets and the boundary of the GRS in Jupiter's atmosphere. We perform our analysis on a velocity field extracted from a video footage acquired by the NASA Cassini spacecraft.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Maxey-Riley Equation: Existence, Uniqueness and Regularity of Solutions

The Maxey--Riley equation describes the motion of an inertial (i.e., finite-size) spherical particle in an ambient fluid flow. The equation is a second-order, implicit integro-differential equation with a singular kernel, and with a forcing term that blows up at the initial time. Despite the widespread use of the equation in applications, the basic properties of its solutions have remained unexplored. Here we fill this gap by proving local existence and uniqueness of weak solutions. For certain initial velocities between the particle and the fluid, the results extend to strong solutions. We also prove continuous differentiability of the weak and strong solutions with respect to their initial conditions. This justifies the search for coherent structures in inertial flows using the Cauchy--Green strain tensor.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Attracting and repelling Lagrangian coherent structures from a single computation

Hyperbolic Lagrangian Coherent Structures (LCSs) are locally most repelling or most attracting material surfaces in a finite-time dynamical system. To identify both types of hyperbolic LCSs at the same time instance, the standard practice has been to compute repelling LCSs from future data and attracting LCSs from past data. This approach tacitly assumes that coherent structures in the flow are fundamentally recurrent, and hence gives inconsistent results for temporally aperiodic systems. Here we resolve this inconsistency by showing how both repelling and attracting LCSs are computable at the same time instance from a single forward or a single backward run. These LCSs are obtained as surfaces normal to the weakest and strongest eigenvectors of the Cauchy-Green strain tensor.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Detecting invariant manifolds, attractors, and generalized KAM tori in aperiodically forced mechanical systems

We show how the recently developed theory of geodesic transport barriers for fluid flows can be used to uncover key invariant manifolds in externally forced, one-degree-of-freedom mechanical systems. Specifically, invariant sets in such systems turn out to be shadowed by least-stretching geodesics of the Cauchy-Green strain tensor computed from the flow map of the forced mechanical system. This approach enables the finite-time visualization of generalized stable and unstable manifolds, attractors and generalized KAM curves under arbitrary forcing, when Poincare maps are not available. We illustrate these results by detailed visualizations of the key finite-time invariant sets of conservatively and dissipatively forced Duffing oscillators.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hyperbolic and Elliptic Transport Barriers in Three-Dimensional Unsteady Flows

We develop a general theory of transport barriers for three-dimensional unsteady flows with arbitrary time-dependence. The barriers are obtained as two-dimensional Lagrangian Coherent Structures (LCSs) that create locally maximal deformation. Along hyperbolic LCSs, this deformation is induced by locally maximal normal repulsion or attraction. Along shear LCSs, the deformation is created by locally maximal tangential shear. Hyperbolic LCSs, therefore, play the role of generalized stable and unstable manifolds, while closed shear LCSs (elliptic LCSs) act as generalized KAM tori or KAM-type cylinders. All these barriers can be computed from our theory as explicitly parametrized surfaces. We illustrate our results by visualizing two-dimensional hyperbolic and elliptic barriers in steady and unsteady versions of the ABC flow.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Shearless transport barriers in unsteady two-dimensional flows and maps

We develop a variational principle that extends the notion of a shearless transport barrier from steady to general unsteady two-dimensional flows and maps defined over a finite time interval. This principle reveals that hyperbolic Lagrangian Coherent Structures (LCSs) and parabolic LCSs (or jet cores) are the two main types of shearless barriers in unsteady flows. Based on the boundary conditions they satisfy, parabolic barriers are found to be more observable and robust than hyperbolic barriers, confirming widespread numerical observations. Both types of barriers are special null-geodesics of an appropriate Lorentzian metric derived from the Cauchy--Green strain tensor. Using this fact, we devise an algorithm for the automated computation of parabolic barriers. We illustrate our detection method on steady and unsteady non-twist maps and on the aperiodically forced Bickley jet.

0 citations0 reviewsNo code linkedOpen access