Source author record

Michael D. Graham

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

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

16 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.

preprint2016arXiv

Low-dimensional representations of exact coherent states of the Navier-Stokes equations from the resolvent model of wall turbulence

We report that many exact invariant solutions of the Navier-Stokes equations for both pipe and channel flows are well represented by just few modes of the model of McKeon & Sharma J. Fl. Mech. 658, 356 (2010). This model provides modes that act as a basis to decompose the velocity field, ordered by their amplitude of response to forcing arising from the interaction between scales. The model was originally derived from the Navier-Stokes equations to represent turbulent flows and has been used to explain coherent structure and to predict turbulent statistics. This establishes a surprising new link between the two distinct approaches to understanding turbulence.

preprint2015arXiv

Exact coherent states and connections to turbulent dynamics in minimal channel flow

Several new families of nonlinear three-dimensional travelling wave solutions to the Navier-Stokes equation, also known as exact coherent states, are computed for Newtonian plane Poiseuille flow. The symmetries and streak/vortex structures are reported and their possible connections to critical layer dynamics examined. While some of the solutions clearly display fluctuations that are localized around the critical layer (the surface on which the streamwise velocity matches the wave speed of the solution), for others this connection is not as clear. Dynamical trajectories along unstable directions of the solutions are computed. Over certain ranges of Reynolds number, two solution families are shown to lie on the basin boundary between laminar and turbulent flow. Direct comparison of nonlinear travelling wave solutions to turbulent flow in the same channel is presented. The state-space dynamics of the turbulent flow are organized around one of the newly-identified travelling wave families, and in particular the lower branch solutions of this family are closely approached during transient excursions away from the dominant behaviour. These observations provide a firm dynamical-systems foundation for prior observations that minimal channel turbulence displays time intervals of "active" turbulence punctuated by brief periods of "hibernation" (see e.g. Xi, L. and Graham, M. D., Phys. Rev. Lett., 104, 218301 (2010)). The hibernating intervals are approaches to lower branch nonlinear travelling waves. Representing these solutions on a Prandtl-von Karman plot illustrates how their bulk flow properties are related to those of Newtonian turbulence as well as the universal asymptotic state called maximum drag reduction (MDR) found in viscoelastic turbulent flow.

preprint2012arXiv

Accelerated boundary integral method for multiphase flow in non-periodic geometries

An accelerated boundary integral method for Stokes flow of a suspension of deformable particles is presented for an arbitrary domain and implemented for the important case of a planar slit geometry. The computational complexity of the algorithm scales as O(N) or $O(N\log N$), where $N$ is proportional to the product of number of particles and the number of elements employed to discretize the particle. This technique is enabled by the use of an alternative boundary integral formulation in which the velocity field is expressed in terms of a single layer integral alone, even in problems with non-matched viscosities. The density of the single layer integral is obtained from a Fredholm integral equation of the second kind involving the double layer integral. Acceleration in this implementation is provided by the use of General Geometry Ewald-like method (GGEM) for computing the velocity and stress fields driven by a set of point forces in the geometry of interest. For the particular case of the slit geometry, a Fourier-Chebyshev spectral discretization of GGEM is developed. Efficient implementations employing the GGEM methodology are presented for the resulting single and the double layer integrals. The implementation is validated with test problems on the velocity of rigid particles and drops between parallel walls in pressure driven flow, the Taylor deformation parameter of capsules in simple shear flow and the particle trajectory in pair collisions of capsules in shear flow. The computational complexity of the algorithm is verified with results from several large scale multiparticle simulations.

preprint2012arXiv

Mechanism of margination in confined flows of blood and other multicomponent suspensions

Flowing blood displays a phenomenon called margination, in which leukocytes and platelets are preferentially found near blood vessel walls, while erythrocytes are depleted from these regions. Here margination is investigated using direct hydrodynamic simulations of a binary suspension of stiff (s) and floppy (f) capsules, as well as a stochastic model that incorporates the key particle transport mechanisms in suspensions -- wall-induced hydrodynamic migration and shear-induced pair collisions. The stochastic model allows the relative importance of these two mechanisms to be directly evaluated and thereby indicates that margination, at least in the dilute case, is largely due to the differential dynamics of homogeneous (e.g. s-s) and heterogeneous (s-f) collisions

preprint2011arXiv

Segregation by membrane rigidity in flowing binary suspensions of elastic capsules

Spatial segregation in the wall normal direction is investigated in suspensions containing a binary mixture of Neo-Hookean capsules subjected to pressure driven flow in a planar slit. The two components of the binary mixture have unequal membrane rigidities. The problem is studied numerically using an accelerated implementation of the boundary integral method. The effect of a variety of parameters was investigated, including the capillary number, rigidity ratio between the two species, volume fraction, confinement ratio, and the number fraction of the more floppy particle $X_f$ in the mixture. It was observed that in suspensions of pure species, the mean wall normal positions of the stiff and the floppy particles are comparable. In mixtures, however, the stiff particles were found to be increasingly displaced towards the walls with increasing $X_f$, while the floppy particles were found to increasingly accumulate near the centerline with decreasing $X_f$. The origins of this segregation is traced to the effect of the number fraction $X_f$ on the localization of the stiff and the floppy particles in the near wall region -- the probability of escape of a stiff particle from the near wall region to the interior is greatly reduced with increasing $X_f$, while the exact opposite trend is observed for a floppy particle with decreasing $X_f$. Simple model studies on heterogeneous pair collisions involving a stiff and a floppy particle mechanistically explain this observation. The key result in these studies is that the stiff particle experiences much larger cross-stream displacement in heterogeneous collisions than the floppy particle. A unified mechanism incorporating the wall-induced migration of deformable particles and the particle fluxes associated with heterogeneous and homogeneous pair collisions is presented.

preprint2010arXiv

Active and hibernating turbulence in minimal channel flow of Newtonian and polymeric fluids

Turbulent channel flow of drag-reducing polymer solutions is simulated in minimal flow geometries. Even in the Newtonian limit, we find intervals of "hibernating" turbulence that display many features of the universal maximum drag reduction (MDR) asymptote observed in polymer solutions: weak streamwise vortices, nearly nonexistent streamwise variations and a mean velocity gradient that quantitatively matches experiments. As viscoelasticity increases, the frequency of these intervals also increases, while the intervals themselves are unchanged, leading to flows that increasingly resemble MDR.

preprint2008arXiv

Dynamics of confined suspensions of swimming particles

Low Reynolds number direct simulations of large populations of hydrodynamically interacting swimming particles confined between planar walls are performed. The results of simulations are compared with a theory that describes dilute suspensions of swimmers. The theory yields scalings with concentration for diffusivities and velocity fluctuations as well as a prediction of the fluid velocity spatial autocorrelation function. Even for uncorrelated swimmers, the theory predicts anticorrelations between nearby fluid elements that correspond to vortex-like swirling motions in the fluid with length scale set by the size of a swimmer and the slit height. Very similar results arise from the full simulations indicating either that correlated motion of the swimmers is not significant at the concentrations considered or that the fluid phase autocorrelation is not a sensitive measure of the correlated motion. This result is in stark contrast with results from unconfined systems, for which the fluid autocorrelation captures large-scale collective fluid structures. The additional length scale (screening length) introduced by the confinement seems to prevent these large-scale structures from forming.