Researcher profile

Michael D. Graham

Michael D. Graham contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 published item(s)

preprint2023arXiv

Dynamics of a data-driven low-dimensional model of turbulent minimal Couette flow

Because the Navier-Stokes equations are dissipative, the long-time dynamics of a flow in state space are expected to collapse onto a manifold whose dimension may be much lower than the dimension required for a resolved simulation. On this manifold, the state of the system can be exactly described in a coordinate system parameterizing the manifold. Describing the system in this low-dimensional coordinate system allows for much faster simulations and analysis. We show, for turbulent Couette flow, that this description of the dynamics is possible using a data-driven manifold dynamics modeling method. This approach consists of an autoencoder to find a low-dimensional manifold coordinate system and a set of ordinary differential equations defined by a neural network. Specifically, we apply this method to minimal flow unit turbulent plane Couette flow at $\textit{Re}=400$, where a fully resolved solutions requires $\mathcal{O}(10^5)$ degrees of freedom. Using only data from this simulation we build models with fewer than $20$ degrees of freedom that quantitatively capture key characteristics of the flow, including the streak breakdown and regeneration cycle. At short-times, the models track the true trajectory for multiple Lyapunov times, and, at long-times, the models capture the Reynolds stress and the energy balance. For comparison, we show that the models outperform POD-Galerkin models with $\sim$2000 degrees of freedom. Finally, we compute unstable periodic orbits from the models. Many of these closely resemble previously computed orbits for the full system; additionally, we find nine orbits that correspond to previously unknown solutions in the full system.

preprint2022arXiv

Data-driven discovery of intrinsic dynamics

Dynamical models underpin our ability to understand and predict the behavior of natural systems. Whether dynamical models are developed from first-principles derivations or from observational data, they are predicated on our choice of state variables. The choice of state variables is driven by convenience and intuition, and in the data-driven case the observed variables are often chosen to be the state variables. The dimensionality of these variables (and consequently the dynamical models) can be arbitrarily large, obscuring the underlying behavior of the system. In truth, these variables are often highly redundant and the system is driven by a much smaller set of latent intrinsic variables. In this study, we combine the mathematical theory of manifolds with the representational capacity of neural networks to develop a method that learns a system's intrinsic state variables directly from time series data, and also learns predictive models for their dynamics. What distinguishes our method is its ability to reduce data to the intrinsic dimensionality of the nonlinear manifold they live on. This ability is enabled by the concepts of charts and atlases from the theory of manifolds, whereby a manifold is represented by a collection of patches that are sewn together -- a necessary representation to attain intrinsic dimensionality. We demonstrate this approach on several high-dimensional systems with low-dimensional behavior. The resulting framework provides the ability to develop dynamical models of the lowest possible dimension, capturing the essence of a system.

preprint2022arXiv

Wrinkling and multiplicity in the dynamics of deformable sheets in uniaxial extensional flow

The processing of thin-structured materials in a fluidic environment, from nearly inextensible but flexible graphene sheets to highly extensible polymer films, arises in many applications. So far, little is known about the dynamics of such thin sheets freely suspended in fluid. In this work, we study the dynamics of freely suspended soft sheets in uniaxial extensional flow. Elastic sheets are modeled with a continuum model that accounts for in-plane deformation and out-of-plane bending, and the fluid motion is computed using the method of regularized Stokeslets. We explore two types of sheets: "stiff" sheets that strongly resist bending deformations and always stay flat, and "flexible" sheets with both in-plane and out-of plane deformability that can wrinkle. For stiff sheets, we observe a coil-stretch-like transition, similar to what has been observed for long-chain linear polymers under extension as well as elastic sheets under planar extension: in a certain range of capillary number (flow strength relative to in-plane deformability), the sheets exhibit either a compact or a highly stretched conformation, depending on deformation history. For flexible sheets, sheets with sufficiently small bending stiffness wrinkle to form various conformations. Here, the compact-stretched bistability still occurs, but is strongly modified by the wrinkling instability: a highly-stretched planar state can become unstable and wrinkle, after which it may dramatically shrink in length due to hydrodynamic screening associated with wrinkling. Therefore, wrinkling renders a shift in the bistability regime. In addition, we can predict and understand the nonlinear long-term dynamics for some parameter regimes with linear stability analysis of the flat steady states.

preprint2021arXiv

Data-Driven Reduced-Order Modeling of Spatiotemporal Chaos with Neural Ordinary Differential Equations

Dissipative partial differential equations that exhibit chaotic dynamics tend to evolve to attractors that exist on finite-dimensional manifolds. We present a data-driven reduced order modeling method that capitalizes on this fact by finding the coordinates of this manifold and finding an ordinary differential equation (ODE) describing the dynamics in this coordinate system. The manifold coordinates are discovered using an undercomplete autoencoder -- a neural network (NN) that reduces then expands dimension. Then the ODE, in these coordinates, is approximated by a NN using the neural ODE framework. Both of these methods only require snapshots of data to learn a model, and the data can be widely and/or unevenly spaced. We apply this framework to the Kuramoto-Sivashinsky for different domain sizes that exhibit chaotic dynamics. With this system, we find that dimension reduction improves performance relative to predictions in the ambient space, where artifacts arise. Then, with the low-dimensional model, we vary the training data spacing and find excellent short- and long-time statistical recreation of the true dynamics for widely spaced data (spacing of ~0.7 Lyapunov times). We end by comparing performance with various degrees of dimension reduction, and find a "sweet spot" in terms of performance vs. dimension.

preprint2020arXiv

Coil-stretch-like transition of elastic sheets in extensional flows

The conformation of a long linear polymer dissolved in fluid and exposed to an extensional flow is well-known to exhibit a "coil-stretch" transition, which for sufficiently long chains can lead to bistability. The present work reports computations indicating that an analogous "compact-stretched" transition arises in the dynamics of a thin elastic sheet. Sheets of nominally circular, square or rectangular shape are simulated in planar and biaxial flows using a finite element method for the sheet conformations and a regularized Stokeslet method for the fluid flow. If a neo-Hookean constitutive model is used for the sheet elasticity, the sheets will stretch without bound once a critical extension rate, as characterized nondimensionally by a capillary number, is exceeded. Nonlinear elasticity, represented with the Yeoh model, arrests the stretching, leading to a highly-stretched steady state once the critical capillary number is exceeded. For all shapes and in both planar and biaxial extension, a parameter regime exists in which both weakly stretched (compact) and strongly stretched states can be found, depending on initial conditions. I.e. this parameter regime displays bistability. As in the long-chain polymer case, the bistable behavior arises from the hydrodynamic interaction between distant elements of the sheet, and vanishes if these interactions are artificially screened by use of a Brinkman model for the fluid motion. While the sheets can transiently display wrinkled shapes, all final shapes in planar and biaxial extension are planar.

preprint2020arXiv

Deep learning to discover and predict dynamics on an inertial manifold

A data-driven framework is developed to represent chaotic dynamics on an inertial manifold (IM), and applied to solutions of the Kuramoto-Sivashinsky equation. A hybrid method combining linear and nonlinear (neural-network) dimension reduction transforms between coordinates in the full state space and on the IM. Additional neural networks predict time-evolution on the IM. The formalism accounts for translation invariance and energy conservation, and substantially outperforms linear dimension reduction, reproducing very well key dynamic and statistical features of the attractor.

preprint2020arXiv

Discovering multiscale and self-similar structure with data-driven wavelets

Many materials, processes, and structures in science and engineering have important features at multiple scales of time and/or space; examples include biological tissues, active matter, oceans, networks, and images. Explicitly extracting, describing, and defining such features are difficult tasks, at least in part because each system has a unique set of features. Here, we introduce an analysis method that, given a set of observations, discovers an energetic hierarchy of structures localized in scale and space. We call the resulting basis vectors a "data-driven wavelet decomposition". We show that this decomposition reflects the inherent structure of the dataset it acts on, whether it has no structure, structure dominated by a single scale, or structure on a hierarchy of scales. In particular, when applied to turbulence---a high-dimensional, nonlinear, multiscale process---the method reveals self-similar structure over a wide range of spatial scales, providing direct, model-free evidence for a century-old phenomenological picture of turbulence. This approach is a starting point for the characterization of localized hierarchical structures in multiscale systems, which we may think of as the building blocks of these systems.

preprint2020arXiv

Flow-induced segregation and dynamics of red blood cells in sickle cell disease

Blood flow in sickle cell disease (SCD) can substantially differ from normal blood flow due to significant alterations in the physical properties of the red blood cells (RBCs). Chronic complications, such as inflammation of the endothelial cells lining blood vessel walls, are associated with SCD, for reasons that are unclear. Here, detailed boundary integral simulations are performed to investigate an idealized model flow flow in SCD, a binary suspension of flexible biconcave discoidal fluid-filled capsules and stiff curved prolate capsules that represent healthy and sickle RBCs, respectively, subjected to pressure-driven flow in a planar slit. The key observation is that, unlike healthy RBCs that concentrate around the center of the channel and form an RBC-depleted layer (i.e. cell-free layer) next to the walls, sickle cells are largely drained from the bulk of the suspension and aggregate inside the cell-free layer, displaying strong margination. These cells are found to undergo a rigid-body-like rolling orbit near the walls. Furthermore, the additional shear stress on the walls induced by the presence of marginated cells is computed. Compared to the small fluctuations in wall shear stress for a suspension of healthy RBCs, large local peaks in wall shear stress are observed for the binary suspensions, due to the proximity of the marginated stiff cells to the walls. As endothelial cells are known to mechanotransduce physical forces such as aberrations in shear stress and convert them to physiological processes such as activation of inflammatory signals, these results may aid in understanding mechanisms for endothelial dysfunction associated with SCD.

preprint2020arXiv

Self-sustained elastoinertial Tollmien-Schlichting waves

Direct simulations of two-dimensional plane channel flow of a viscoelastic fluid at Reynolds number Re = 3000 reveal the existence of a family of attractors whose structure closely resembles the linear Tollmien-Schlichting (TS) mode, and in particular exhibits strongly localized stress fluctuations at the critical layer position of the TS mode. At the parameter values chosen, this solution branch is not connected to the nonlinear TS solution branch found for Newtonian flow, and thus represents a solution family that is nonlinearly self-sustained by viscoelasticity. The ratio between stress and velocity fluctuations is in quantitative agreement for the attractor and the linear TS mode, and increases strongly with Weissenberg number, Wi. For the latter, there is a transition in the scaling of this ratio as Wi increases, and the Wi at which the nonlinear solution family comes into existence is just above this transition. Finally, evidence indicates that this branch is connected through an unstable solution branch to two-dimensional elastoinertial turbulence (EIT). These results suggest that, in the parameter range considered here, the bypass transition leading to EIT is mediated by nonlinear amplification and self-sustenance of perturbations that excite the Tollmien-Schlichting mode.