RESUMO
A continuum description is essential for understanding a variety of collective phenomena in active matter. However, building quantitative continuum models of active matter from first principles can be extremely challenging due to both the gaps in our knowledge and the complicated structure of nonlinear interactions. Here, we use a physically informed data-driven approach to construct a complete mathematical model of an active nematic from experimental data describing kinesin-driven microtubule bundles confined to an oil-water interface. We find that the structure of the model is similar to the Leslie-Ericksen and Beris-Edwards models, but there are appreciable and important differences. Rather unexpectedly, elastic effects are found to play no role in the experiments considered, with the dynamics controlled entirely by the balance between active stresses and friction stresses.
RESUMO
Despite a long and rich history of scientific investigation, fluid turbulence remains one of the most challenging problems in science and engineering. One of the key outstanding questions concerns the role of coherent structures that describe frequently observed patterns embedded in turbulence. It has been suggested, but not proved, that coherent structures correspond to unstable, recurrent solutions of the governing equation of fluid dynamics. Here, we present experimental and numerical evidence that three-dimensional turbulent flow tracks, episodically but repeatedly, the spatial and temporal structure of multiple such solutions. Our results provide compelling evidence that coherent structures, grounded in the governing equations, can be harnessed to predict how turbulent flows evolve.
RESUMO
Machine learning offers an intriguing alternative to first-principle analysis for discovering new physics from experimental data. However, to date, purely data-driven methods have only proven successful in uncovering physical laws describing simple, low-dimensional systems with low levels of noise. Here we demonstrate that combining a data-driven methodology with some general physical principles enables discovery of a quantitatively accurate model of a non-equilibrium spatially extended system from high-dimensional data that is both noisy and incomplete. We illustrate this using an experimental weakly turbulent fluid flow where only the velocity field is accessible. We also show that this hybrid approach allows reconstruction of the inaccessible variables - the pressure and forcing field driving the flow.
RESUMO
In laboratory studies and numerical simulations, we observe clear signatures of unstable time-periodic solutions in a moderately turbulent quasi-two-dimensional flow. We validate the dynamical relevance of such solutions by demonstrating that turbulent flows in both experiment and numerics transiently display time-periodic dynamics when they shadow unstable periodic orbits (UPOs). We show that UPOs we computed are also statistically significant, with turbulent flows spending a sizable fraction of the total time near these solutions. As a result, the average rates of energy input and dissipation for the turbulent flow and frequently visited UPOs differ only by a few percent.
RESUMO
Sparse regression has recently emerged as an attractive approach for discovering models of spatiotemporally complex dynamics directly from data. In many instances, such models are in the form of nonlinear partial differential equations (PDEs); hence sparse regression typically requires the evaluation of various partial derivatives. However, accurate evaluation of derivatives, especially of high order, is infeasible when the data are noisy, which has a dramatic negative effect on the result of regression. We present an alternative and rather general approach that addresses this difficulty by using a weak formulation of the problem. For instance, it allows accurate reconstruction of PDEs involving high-order derivatives, such as the Kuramoto-Sivashinsky equation, from data with a considerable amount of noise. The flexibility of our approach also allows reconstruction of PDE models that involve latent variables which cannot be measured directly with acceptable accuracy. This is illustrated by reconstructing a model for a weakly turbulent flow in a thin fluid layer, where neither the forcing nor the pressure field is known.
RESUMO
This paper investigates how models of spatiotemporal dynamics in the form of nonlinear partial differential equations can be identified directly from noisy data using a combination of sparse regression and weak formulation. Using the 4th-order Kuramoto-Sivashinsky equation for illustration, we show how this approach can be optimized in the limits of low and high noise, achieving accuracy that is orders of magnitude better than what existing techniques allow. In particular, we derive the scaling relation between the accuracy of the model, the parameters of the weak formulation, and the properties of the data, such as its spatial and temporal resolution and the level of noise.
RESUMO
In spatially extended systems, it is common to find latent variables that are hard, or even impossible, to measure with acceptable precision but are crucially important for the proper description of the dynamics. This substantially complicates construction of an accurate model for such systems using data-driven approaches. The present paper illustrates how physical constraints can be employed to overcome this limitation using the example of a weakly turbulent quasi-two-dimensional Kolmogorov flow driven by a steady Lorenz force with an unknown spatial profile. Specifically, the terms involving latent variables in the partial differential equations governing the dynamics can be eliminated at the expense of raising the order of that equation. We show that local polynomial interpolation combined with sparse regression can handle data on spatiotemporal grids that are representative of typical experimental measurement techniques such as particle image velocimetry. However, we also find that the reconstructed model is sensitive to measurement noise and trace this sensitivity to the presence of high-order spatial and/or temporal derivatives.
RESUMO
Recent studies suggest that unstable recurrent solutions of the Navier-Stokes equation provide new insights into dynamics of turbulent flows. In this study, we compute an extensive network of dynamical connections between such solutions in a weakly turbulent quasi-two-dimensional Kolmogorov flow that lies in the inversion-symmetric subspace. In particular, we find numerous isolated heteroclinic connections between different types of solutions-equilibria, periodic, and quasiperiodic orbits-as well as continua of connections forming higher-dimensional connecting manifolds. We also compute a homoclinic connection of a periodic orbit and provide strong evidence that the associated homoclinic tangle forms the chaotic repeller that underpins transient turbulence in the symmetric subspace.
RESUMO
The motion of and interaction between phase singularities that lie at the centers of spiral waves capture many qualitative and, in some cases, quantitative features of complex dynamics in excitable systems. Being able to accurately reconstruct their position is thus quite important, even if the data are noisy and sparse, as in electrophysiology studies of cardiac arrhythmias, for instance. A recently proposed global topological approach [Marcotte and Grigoriev, Chaos 27, 093936 (2017)] promises to meaningfully improve the quality of the reconstruction compared with traditional, local approaches. Indeed, we found that this approach is capable of handling noise levels exceeding the range of the signal with minimal loss of accuracy. Moreover, it also works successfully with data sampled on sparse grids with spacing comparable to the mean separation between the phase singularities for complex patterns featuring multiple interacting spiral waves.
Assuntos
Arritmias Cardíacas/fisiopatologia , Fibrilação Atrial/patologia , Fibrilação Atrial/terapia , Ablação por Cateter , Coração/fisiopatologia , Simulação por Computador , Humanos , Modelos Cardiovasculares , Nó Sinoatrial/fisiologiaRESUMO
Recent studies suggest that unstable, nonchaotic solutions of the Navier-Stokes equation may provide deep insights into fluid turbulence. In this article, we present a combined experimental and numerical study exploring the dynamical role of unstable equilibrium solutions and their invariant manifolds in a weakly turbulent, electromagnetically driven, shallow fluid layer. Identifying instants when turbulent evolution slows down, we compute 31 unstable equilibria of a realistic two-dimensional model of the flow. We establish the dynamical relevance of these unstable equilibria by showing that they are closely visited by the turbulent flow. We also establish the dynamical relevance of unstable manifolds by verifying that they are shadowed by turbulent trajectories departing from the neighborhoods of unstable equilibria over large distances in state space.
RESUMO
The mechanisms underlying cardiac fibrillation have been investigated for over a century, but we are still finding surprising results that change our view of this phenomenon. The present study focuses on the transition from normal rhythm to spiral wave chaos associated with a gradual increase in the pacing rate. While some of our findings are consistent with existing experimental, numerical, and theoretical studies of this problem, one result appears to contradict the accepted picture. Specifically we show that, in a two-dimensional model of paced homogeneous atrial tissue, transition from discordant alternans to conduction block, wave breakup, reentry, and spiral wave chaos is associated with the transient growth of finite amplitude disturbances rather than a conventional instability. It is mathematically very similar to subcritical, or bypass, transition from laminar fluid flow to turbulence, which allows many of the tools developed in the context of fluid turbulence to be used for improving our understanding of cardiac arrhythmias.
Assuntos
Função Atrial/fisiologia , Estimulação Cardíaca Artificial , Modelos Cardiovasculares , Fatores de TempoRESUMO
While spiral wave breakup has been implicated in the emergence of atrial fibrillation, its role in maintaining this complex type of cardiac arrhythmia is less clear. We used the Karma model of cardiac excitation to investigate the dynamical mechanisms that sustain atrial fibrillation once it has been established. The results of our numerical study show that spatiotemporally chaotic dynamics in this regime can be described as a dynamical equilibrium between topologically distinct types of transitions that increase or decrease the number of wavelets, in general agreement with the multiple wavelets' hypothesis. Surprisingly, we found that the process of continuous excitation waves breaking up into discontinuous pieces plays no role whatsoever in maintaining spatiotemporal complexity. Instead, this complexity is maintained as a dynamical balance between wave coalescence-a unique, previously unidentified, topological process that increases the number of wavelets-and wave collapse-a different topological process that decreases their number.
Assuntos
Fibrilação Atrial/fisiopatologia , Modelos Cardiovasculares , Sistema de Condução Cardíaco/fisiopatologia , Análise Numérica Assistida por Computador , Fatores de TempoRESUMO
The existence and dynamical role of particular unstable solutions (exact coherent structures) of the Navier-Stokes equation is revealed in laboratory studies of weak turbulence in a thin, electromagnetically driven fluid layer. We find that the dynamics exhibit clear signatures of numerous unstable equilibrium solutions, which are computed using a combination of flow measurements from the experiment and fully resolved numerical simulations. We demonstrate the dynamical importance of these solutions by showing that turbulent flows visit their state space neighborhoods repeatedly. Furthermore, we find that the unstable manifold associated with one such unstable equilibrium predicts the evolution of turbulent flow in both experiment and simulation for a considerable period of time.
RESUMO
Channel flow of an incompressible fluid at Reynolds numbers above 2400 possesses a number of different spatially localized solutions that approach laminar flow far upstream and downstream. We use one such relative time-periodic solution, which corresponds to a spatially localized version of a Tollmien-Schlichting wave, to illustrate how the upstream and downstream asymptotics can be computed analytically. In particular, we show that for these spanwise uniform states the asymptotics predict the exponential localization that has been observed for numerically computed solutions of several canonical shear flows but never properly understood theoretically.
RESUMO
This paper introduces a numerical method for computing the spectrum of adjoint (left) eigenfunctions of spiral wave solutions to reaction-diffusion systems in arbitrary geometries. The method is illustrated by computing over a hundred eigenfunctions associated with an unstable time-periodic single-spiral solution of the Karma model on a square domain. We show that all leading adjoint eigenfunctions are exponentially localized in the vicinity of the spiral tip, although the marginal modes (response functions) demonstrate the strongest localization. We also discuss the implications of the localization for the dynamics and control of unstable spiral waves. In particular, the interaction with no-flux boundaries leads to a drift of spiral waves which can be understood with the help of the response functions.
RESUMO
This paper investigates the properties of unstable single-spiral wave solutions arising in the Karma model of two-dimensional cardiac tissue. In particular, we discuss how such solutions can be computed numerically on domains of arbitrary shape and study how their stability, rotational frequency, and spatial drift depend on the size of the domain as well as the position of the spiral core with respect to the boundaries. We also discuss how the breaking of local Euclidean symmetry due to finite size effects as well as the spatial discretization of the model is reflected in the structure and dynamics of spiral waves. This analysis allows identification of a self-sustaining process responsible for maintaining the state of spiral chaos featuring multiple interacting spirals.
Assuntos
Modelos Cardiovasculares , Miocárdio , Dinâmica não Linear , Animais , HumanosRESUMO
Unstable nonchaotic solutions embedded in the chaotic attractor can provide significant new insight into chaotic dynamics of both low- and high-dimensional systems. In particular, in turbulent fluid flows, such unstable solutions are referred to as exact coherent structures (ECS) and play an important role in both initiating and sustaining turbulence. The nature of ECS and their role in organizing spatiotemporally chaotic dynamics, however, is reasonably well understood only for systems on relatively small spatial domains lacking continuous Euclidean symmetries. Construction of ECS on large domains and in the presence of continuous translational and/or rotational symmetries remains a challenge. This is especially true for models of excitable media which display spiral turbulence and for which the standard approach to computing ECS completely breaks down. This paper uses the Karma model of cardiac tissue to illustrate a potential approach that could allow computing a new class of ECS on large domains of arbitrary shape by decomposing them into a patchwork of solutions on smaller domains, or tiles, which retain Euclidean symmetries locally.
Assuntos
Potenciais de Ação/fisiologia , Relógios Biológicos/fisiologia , Sistema de Condução Cardíaco/fisiologia , Modelos Cardiovasculares , Dinâmica não Linear , Animais , Simulação por Computador , HumanosRESUMO
Alternans-an arrhythmic response of cardiac tissue to periodic pacing-often serves as a precursor to a more dangerous, and potentially lethal, state of fibrillation. Suppression of alternans using feedback control may be a plausible method to prevent fibrillation. Several approaches based on impulsive control have been proposed previously, where feedback is applied for a brief instance of time during each pacing interval. This paper presents a continuous-time approach, where feedback current is applied at all times, which is capable of suppressing alternans in fibers of significantly greater length (up to at least 4 cm), compared with impulsive control (less than 1 cm), and for a wide range of pacing cycle lengths.
Assuntos
Fibrilação Atrial/fisiopatologia , Modelos Cardiovasculares , Modelos Neurológicos , Ramos Subendocárdicos/fisiopatologia , Fibrilação Atrial/terapia , Estimulação Cardíaca Artificial , HumanosRESUMO
In this article, we compare quantitatively the efficiency of three different protocols commonly used in commercial defibrillators. These are based on monophasic and both symmetric and asymmetric biphasic shocks. A numerical one-dimensional model of cardiac tissue using the bidomain formulation is used in order to test the different protocols. In particular, we performed a total of 4.8 × 10(6) simulations by varying shock waveform, shock energy, initial conditions, and heterogeneity in internal electrical conductivity. Whenever the shock successfully removed the reentrant dynamics in the tissue, we classified the mechanism. The analysis of the numerical data shows that biphasic shocks are significantly more efficient (by about 25%) than the corresponding monophasic ones. We determine that the increase in efficiency of the biphasic shocks can be explained by the higher proportion of newly excited tissue through the mechanism of direct activation.
Assuntos
Arritmias Cardíacas/fisiopatologia , Cardioversão Elétrica , Modelos Cardiovasculares , Miocárdio , Arritmias Cardíacas/terapia , HumanosRESUMO
This paper describes a systematic approach to suppressing cardiac alternans in simulated Purkinje fibers using localized current injections. We investigate the controllability and observability of the periodically paced Noble model for different locations of the recording and control electrodes. In particular, we show that the loss of controllability causes the failure of the control approach introduced by Echebarria and Karma [Chaos 12, 923 (2002)] for longer fiber lengths. Furthermore, we explain how the optimal locations for the recording and control electrodes and the timing of the feedback current can be selected, accounting for both linear and nonlinear effects, effectively doubling the length of fibers that can be controlled with previous methods.