Scroll Waves in Anisotropic Excitable Media with Application to the Heart. Sima Setayeshgar Department of Physics Indiana University
|
|
- Timothy Osborne
- 5 years ago
- Views:
Transcription
1 Scroll Waves in Anisotropic Excitable Media with Application to the Heart Sima Setayeshgar Department of Physics Indiana University KITP Cardiac Dynamics Mini-Program 1
2 Stripes, Spots and Scrolls KITP Cardiac Dynamics Mini-Program 2
3 Overview Common processes underlying natural patterns give rise to model equations capturing generic features of pattern-forming systems Theoretical approaches accessible near onset, but must resort to numerical tools and experiments far from onset Need for high-fidelity scientific computation to describe realistic physical systems as a bridge between theory and experiment KITP Cardiac Dynamics Mini-Program 3
4 KITP Cardiac Dynamics Mini-Program 4
5 The Heart as a Physical System Sudden cardiac failure is the leading cause of death in industrialized nations deaths/day in North America Growing experimental evidence that self-sustained patterns of electrical activity in cardiac tissue are related to fatal arrhythmias. Goal is to use analytical and numerical tools to study the dynamics of reentrant waves in the heart on physiologically realistic domains. And... The heart is an interesting arena for applying the ideas of pattern formation. KITP Cardiac Dynamics Mini-Program 5
6 Big Picture What are the mechanisms for transition from ventricular tachychardia to fibrillation? How can we control it? Tachycardia: Fibrillation: Courtesy of Sasha Panfilov, University of Utrecht Click for animation. Paradigm: Breakdown of single spiral to disordered state resulting from various mechanisms of spiral instability. KITP Cardiac Dynamics Mini-Program 6
7 Focus What is the role of the anisotropy inherent in the fiber architecture of the heart on scroll wave dynamics? Motivated by: A. T. Winfree, in Progress in Biophysics and Molecular Biology, D. Noble et al. eds., (1997). Numerical ``experiments'' In rectangular slab geometries: Panfilov, A. V. and Keener, J. P., Physica D 84, 545 (1995): Scroll wave breakup due to rotating anisotropy. Fenton, F. and Karma, A., Chaos 8, 20 (1998): Rotating anisotropy leading to ``twistons'', eventually destabilizing scroll filament. Analytical work Dynamics of scroll waves in isotropic excitable media, beginning with: Keener, J. P., Physica D 31, 269 (1988). Biktashev, V. N., Physica D 36, 167 (1989). KITP Cardiac Dynamics Mini-Program 7
8 Tissue Structure 3d conduction pathway with uniaxial anisotropy Propagation speeds: c = 0.5 m/s, c = 0.17 m/s From Textbook of Medical Physiology, by Guyton and Hall From Thomas, Am. J. Anatomy (1957). KITP Cardiac Dynamics Mini-Program 8
9 Credits Collaborator Andrew Bernoff, Mathematics Department, Harvey Mudd College Acknowledgements Alain Karma, Physics Department, Northeastern University Herb Keller, Applied Mathematics Department, Caltech Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 9
10 Rotating anisotropy D D D 2 D D D 2 from Streeter, et al., Circ. Res. 24, p. 339 (1969). Diffusion constants: D > D 1 D 2 Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 10
11 Coordinate System Natural coordinate system defined by fiber direction: x ỹ z } {{ } `new' = α }{{} S cos Θ(z) sin Θ(z) 0 sin Θ(z) cos Θ(z) }{{} R x y z } {{ } `old' S : rescaling, according to 2d anisotropy α (D 1 /D ) 1/2 R : rotation, according to fiber direction Θ(z) Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 11
12 Governing Equations Governing reaction-diffusion equation in new coordinates: u t = f( u) + D 2 u + D 2 u zz { [ + D 2 Θ 2 2 θ + 2 (α2 1)x 2 2 y + 2 2Θ [ θ Θ [ θ + (α 1)x y + (α 1)x y ( 1 )y α x 2 ) y ] u ( 1 α 1 ( ) 1 α 1 x ] y x z } u, ] u Only depends on fiber rotation rate, Θ (no explicit dependence on Θ(z)). For FitzHugh-Naguomo (FHN) kinetics: ( ) ( u u =, f u 3 + 3u v = v ɛ(u δ) ), D = ( D ), etc... Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 12
13 Peskin Fiber Distribution Profile Measured Derived from Streeter, et al., Circ. Res. 24, p. 339 (1969) from Peskin, et al., Comm. on Pure and Appl. Math 42, p. 79 (1989) Θ(z) = sin 1 (z/rl) r = cutoff parameter 2L = thickness of ventricular wall Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 13
14 Perturbation Analysis Consider the limit of `small rotating anisotropy' : Non-dimensional small parameter: ɛ 2 = D 2 ω 0 L 2 1 r 2 1 ( γ ) ( D 2 ω 0 ) 1/2 2L r γ = α + 1/α : transverse diffusion length, l : thickness of ventricular wall : cutoff parameter : `anisotropy' Seek a solution in the form of: where U 0 (r, θ ω 0 t) satisfies: u = U 0 (r, θ ω 0 t + Θ(z) + φ(z, t)) + ɛ 2 u 2, O(1) : U 0 t = f( U 0 ) + D 2 U 0 Scaling assumptions: u 2 O(1), φ z O(ɛ), φ t O(ɛ 2 ). Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 14
15 Validity of Perturbation Analysis? What is the value of the small parameter for the human ventricle? Parameter Value D 1.0 cm 2 s 1 D 0.1 cm 2 s 1 ω s 1 Θ 180 2L 1.0 cm r 1.5 ɛ Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 15
16 Phase Equation At O(ɛ 2 ), introducing Φ(z, t) ( c1 c 2 ) [φ(z, t) ( γ 2 1) Θ(z) ] : Φ t Φ 2 z Φ zz = A(γ, r) F (z; r), 1 < z < 1 Burgers' equation, with forcing given by fiber rotation: F (z; r) = 1 1/r2 1 (z/r) 2, ( A(γ, r) = Ã γ 2 ) r 2 1, Ã = ( ) 2 ( ) c1 4a1 c 2 c 1 1 (a i, c i ) given by inner products from the solvability condition, e.g., a 1 = Ỹθ, D 2 U 0θθ Seek asymptotic and numerical solutions, using constant frequency-shift ansatz: Φ(z, t) = z 1 k(z ) }{{} twist dz + λt + Φ 0 Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 16
17 Scroll Twist For a straight scroll: k(z, t) = ( ) ˆN z ˆN ẑ ˆN = u/ u normal to tip trajectory In new coordinates: k(z, t) = φ z (z, t) + Θ (z). In old coordinates: k( z, t) = Θ ( z) 2α ( φ z ( z, t) + Θ ( z) ) (α 2 1) cos [2 (ω 0 t φ( z, t) Θ( z))] + (1 + α 2 ). Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 17
18 Twist Solutions Diffusive Regime: Φ zz Φ 2 z Twist-dominated Regime: Φ zz Φ 2 z A > 0, A < 0: Maximum twist at boundary A > 0: Formation of large twist in boundary layer in bulk A < 0: Expulsion of large twist from bulk to boundaries Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 18
19 Relevance? Henzi, Lugosi, Winfree, Can. J. Phys. 68, 683 (1990): α = 1 Helical buckling (``sproing instability'') for twist > twist Tip Trajectory Scroll Period vs. Scroll Twist With rotating anisotropy, α 1, Θ z 0??? Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 19
20 Summary What has been done: Extension of asymptotics of scroll waves in isotropic media to include rotating anisotropy of cardiac tissue Phase dynamics (forced Burgers equation): nonconstant fiber rotation rate generates twist a Extensions: Coupling between twist and filament dynamics Extension to biodomain description of cardiac tissue a : Setayeshgar and Bernoff, Phys. Rev. Lett. (2002). Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 20
21 Summary (cont'd) Numerical sproing bifurcation diagram with rotating anisotropy Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 21
22 Filament motion Lab Frame: Fiber Frame: Moving Frame: Simplifying assumptions: Constant fiber rotation rate: Θ = constant φ(z, 0)/ z = 0, R c (z, 0)/ z = 0 in fiber frame No-flux or periodic boundary conditions at vertical boundaries Dynamics reduces to two dimensions: ``Helical'' buckling Motion of spiral center. Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 22
23 Additional coordinate transformation: Filament motion (cont'd) x x X c (t), t y y Y c (t), Phase equation: φ T = Θ 2 Z t + dx c dt x + dy c dt y, [ ( )] 1 c 3 (α) + r 1 α 2Y c 2 + α 2 2 X c Dynamics of the center: d dt ( Xc Y c ) = Θ 2 Z ( µ1 µ 2 ) } µ 3 µ 4 {{ } M ( Xc Y c ) where µ i = µ i (α): µ 1 (1) = µ 4 (1) = t 1 µ 2 (1) = µ 3 (1) = t 2. with: t 1 = Ỹx, D 2 x U 0 x = Ỹy, D 2 y U 0 y, etc. r 1 = Ỹ0, D 2 U 0 xx = Ỹ0, D 2 U 0 yy Notes: Symmetry: α 1/α At O(ɛ 2 ) with φ z = 0: Motion of center is decoupled from dynamics of phase. Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 23
24 Filament Tension Dynamics of spiral center governed by eigenvalues of M: R c (T ) = C + v + e λ + T + C v e λ T No anisotropy: α = 1 λ ± = t 1 ± it 2 Stability determined by filament tension, ( t 1 ). At O(ɛ 2 ), Θ Z determines only the scaling of filament dynamics. Weak anisotropy: α 1 δ, δ 1 λ ± t 1 ± t δ2 (B 2 A 2 ) Rotating anisotropy can lead to change in stability! (Necessary cond'n: B > A.) Dependence on reaction kinetics of α = 1: Filament tension (threshold to buckling) α 1: Possible destabilizing role of rotating anisotropy on buckling Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 24
25 Formally valid for: A 1 Twist-dominated Regime: Φ zz Φ 2 z Cole-Hopf transformation: k(z) = ψ z (z)/ψ(z), d 2 ψ/dz 2 + [ λ V (z)]ψ = 0, V (z) = A(γ, r)f (z; r). Ground state (smallest λ ) determined by potential in the vicinity of its minimum: Negative forcing: A < 0 1d harmonic oscillator equation, λ determined by behaviour at the origin: ( ) λ 0 = Ā/r 2, λ 1 = Ā 1/2 /r 2 Ā = Ã γ 2 /4 1 Positive forcing: A > 0 Airy equation, λ determined by behaviour at boundaries: λ 0 = Ā/(r 2 1), λ 1 = η ( 2Ā ) 2/3 /(r 2 1) 4/3 where η is the first zero of Ai (z). Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 25
26 Twist-dominated Regime: Φ zz Φ 2 z (cont'd) A > 0: Formation of large twist in boundary layer in bulk A < 0: Expulsion of large twist from bulk to boundaries Sima Setayeshgar, Indiana University, Bloomington KITP Cardiac Dynamics Mini-Program 26
Scroll Waves in Anisotropic Excitable Media with Application to the Heart. Sima Setayeshgar Department of Physics Indiana University
Scroll Waves in Anisotropic Excitable Media with Application to the Heart Sima Setayeshgar Department of Physics Indiana University KITP Cardiac Dynamics Mini-Program 1 Stripes, Spots and Scrolls KITP
More informationThe Physics of the Heart. Sima Setayeshgar
The Physics of the Heart Sima Setayeshgar Department of Physics Indiana University Indiana Unversity Physics REU Seminar, July 27 2005 1 Stripes, Spots and Scrolls Indiana Unversity Physics REU Seminar,
More informationThe Physics of the Heart. Sima Setayeshgar
The Physics of the Heart Sima Setayeshgar Department of Physics Indiana University Indiana Unversity Physics REU Seminar: August 1, 2007 1 Stripes, Spots and Scrolls Indiana Unversity Physics REU Seminar:
More informationIntroduction to Physiology V - Coupling and Propagation
Introduction to Physiology V - Coupling and Propagation J. P. Keener Mathematics Department Coupling and Propagation p./33 Spatially Extended Excitable Media Neurons and axons Coupling and Propagation
More informationApplication to Experiments
Chapter 12 Application to Experiments In the past three decades there have been a large number of papers reporting the experimental observation of chaos. In this chapter I describe two rather different
More informationROLE OF BIDOMAIN MODEL OF CARDIAC TISSUE IN THE DYNAMICS OF PHASE SINGULARITIES
ROLE OF BIDOMAIN MODEL OF CARDIAC TISSUE IN THE DYNAMICS OF PHASE SINGULARITIES Jianfeng Lv Submitted to the faculty of the University Graduate School in partial fulfillment of the requirement for the
More informationarxiv:patt-sol/ v1 3 Oct 1997
Propagation Failure in Excitable Media A. Hagberg Center for Nonlinear Studies and T-7, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545 E. Meron The Jacob Blaustein Institute
More informationPCMI Project: Resetting Reentrant Excitation Oscillations in Different Geometries
PCMI Project: Resetting Reentrant Excitation Oscillations in Different Geometries Elizabeth Doman mentor: John Milton summer 2005 PCMI Project:Resetting Reentrant ExcitationOscillations in Different Geometries
More informationBoundary-induced reentry in homogeneous excitable tissue
Boundary-induced reentry in homogeneous excitable tissue Fernando Siso-Nadal, 1 Niels F. Otani, 2 Robert F. Gilmour, Jr., 2 and Jeffrey J. Fox 1 1 Gene Network Sciences, Ithaca, New York 1485, USA 2 Department
More informationSpatiotemporal Chaos in Rayleigh-Bénard Convection
Spatiotemporal Chaos in Rayleigh-Bénard Convection Michael Cross California Institute of Technology Beijing Normal University June 2006 Janet Scheel, Keng-Hwee Chiam, Mark Paul Henry Greenside, Anand Jayaraman
More informationConduction velocity restitution in models of electrical excitation in the heart
Conduction velocity restitution in models of electrical excitation in the heart Who? From? When? School of Mathematics and Statistics University of Glasgow December 7, 2011 / Liverpool References and credits
More informationDrift velocity of rotating spiral waves in the weak deformation approximation
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 8 22 AUGUST 2003 Drift velocity of rotating spiral waves in the weak deformation approximation Hong Zhang a) Bambi Hu and Department of Physics, University
More informationCHEM 515: Chemical Kinetics and Dynamics
Alejandro J. Garza S01163018 Department of Chemistry, Rice University, Houston, TX email: ajg7@rice.edu, ext. 2657 Submitted December 12, 2011 Abstract Spontaneous antispiral wave formation was observed
More informationBifurcation analysis of a normal form for excitable media: Are stable dynamical alternans on a ring possible?
CHAOS 18, 013129 2008 Bifurcation analysis of a normal form for excitable media: Are stable dynamical alternans on a ring possible? Georg A. Gottwald a School of Mathematics and Statistics and Centre for
More informationLong-wave Instability in Anisotropic Double-Diffusion
Long-wave Instability in Anisotropic Double-Diffusion Jean-Luc Thiffeault Institute for Fusion Studies and Department of Physics University of Texas at Austin and Neil J. Balmforth Department of Theoretical
More informationLecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps
Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion Poincare-Bendixson theorem Hopf bifurcations Poincare maps Limit
More information2013 NSF-CMACS Workshop on Atrial Fibrillation
2013 NSF-CMACS Workshop on A Atrial Fibrillation Flavio H. Fenton School of Physics Georgia Institute of Technology, Atlanta, GA and Max Planck Institute for Dynamics and Self-organization, Goettingen,
More informationMultiscale Diffusion Modeling in Charged and Crowded Biological Environments
Multiscale Diffusion Modeling in Charged and Crowded Biological Environments Andrew Gillette Department of Mathematics University of Arizona joint work with Pete Kekenes-Huskey (U. Kentucky) and J. Andrew
More informationTensor Visualization. CSC 7443: Scientific Information Visualization
Tensor Visualization Tensor data A tensor is a multivariate quantity Scalar is a tensor of rank zero s = s(x,y,z) Vector is a tensor of rank one v = (v x,v y,v z ) For a symmetric tensor of rank 2, its
More informationSpiral waves on static and moving spherical domains
Journal of Computational and Applied Mathematics 182 2005) 472 486 www.elsevier.com/locate/cam Spiral waves on static and moving spherical domains Faridon Amdjadi, Jagnnathan Gomatam School of Computing
More informationWave propagation in an excitable medium with a negatively sloped restitution curve
CHAOS VOLUME 12, NUMBER 3 SEPTEMBER 2002 Wave propagation in an excitable medium with a negatively sloped restitution curve A. V. Panfilov a) Department of Theoretical Biology, Utrecht University, Padualaan
More informationProblems in diffusion and absorption: How fast can you hit a target with a random walk?
Problems in diffusion and absorption: How fast can you hit a target with a random walk? Andrew J. Bernoff Harvey Mudd College In collaboration with Alan Lindsay (Notre Dame) Thanks to Alan Lindsay, Michael
More informationMarch 9, :18 Int J. Bifurcation and Chaos/INSTRUCTION FILE Morfu2v2 EFFECT OF NOISE AND STRUCTURAL INHOMOGENEITIES IN REACTION DIFFUSION MEDIA
March 9, 2007 10:18 Int J. Bifurcation and Chaos/INSTRUCTION FILE Int J. Bifurcation and Chaos Submission Style EFFECT OF NOISE AND STRUCTURAL INHOMOGENEITIES IN REACTION DIFFUSION MEDIA S. Morfu Laboratoire
More information2011 NSF-CMACS Workshop on Atrial Fibrillation (5 th day )
2011 NSF-CMACS Workshop on Atrial Fibrillation (5 th day ) Flavio H. Fenton Department of Biomedical Sciences College of Veterinary Medicine, Cornell University, NY and Max Planck Institute for Dynamics
More informationScroll waves in spherical shell geometries
CHAOS VOLUME 11, NUMBER 4 DECEMBER 2001 Scroll waves in spherical shell geometries Francisco Chávez and Raymond Kapral Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Ontario
More informationModelling biological oscillations
Modelling biological oscillations Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Outline Outline of Topics Van der Pol equation Van
More informationMultiple Mechanisms of Spiral Wave Breakup in a Model of Cardiac Electrical Activity
Multiple Mechanisms of Spiral Wave Breakup in a Model of Cardiac Electrical Activity Flavio H. Fenton and Elizabeth M. Cherry Center for Arrhythmia Research at Hofstra University and The Heart Institute,
More informationPattern Formation and Spatiotemporal Chaos in Systems Far from Equilibrium
Pattern Formation and Spatiotemporal Chaos in Systems Far from Equilibrium Michael Cross California Institute of Technology Beijing Normal University May 2006 Michael Cross (Caltech, BNU) Pattern Formation
More informationQuasipatterns in surface wave experiments
Quasipatterns in surface wave experiments Alastair Rucklidge Department of Applied Mathematics University of Leeds, Leeds LS2 9JT, UK With support from EPSRC A.M. Rucklidge and W.J. Rucklidge, Convergence
More informationLecture Notes 8C120 Inleiding Meten en Modelleren. Cellular electrophysiology: modeling and simulation. Nico Kuijpers
Lecture Notes 8C2 Inleiding Meten en Modelleren Cellular electrophysiology: modeling and simulation Nico Kuijpers nico.kuijpers@bf.unimaas.nl February 9, 2 2 8C2 Inleiding Meten en Modelleren Extracellular
More informationMeandering Spiral Waves Induced by Time-Periodic Coupling Strength
Commun. Theor. Phys. 60 (2013) 545 550 Vol. 60, No. 5, November 15, 2013 Meandering Spiral Waves Induced by Time-Periodic Coupling Strength WANG Mao-Sheng ( á), SUN Run-Zhi (ê ì), HUANG Wan-Xia (á ), TU
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationParameters for Minimal Model of Cardiac Cell from Two Different Methods: Voltage-Clamp and MSE Method
Parameters for Minimal Model of Cardiac Cell from Two Different Methods: oltage-clamp and MSE Method Soheila Esmaeili 1, * and Bahareh beheshti 1 Department of Biomedical engineering, ran University of
More informationChimera states in networks of biological neurons and coupled damped pendulums
in neural models in networks of pendulum-like elements in networks of biological neurons and coupled damped pendulums J. Hizanidis 1, V. Kanas 2, A. Bezerianos 3, and T. Bountis 4 1 National Center for
More informationExpanding scroll rings and negative tension turbulence in a model of excitable media
PHYSICAL REVIEW E 70, 056201 (2004) Expanding scroll rings and negative tension turbulence in a model of excitable media S. Alonso, 1, * Ralf Kähler, 2 A. S. Mikhailov, 3 and F. Sagués 1 1 Departament
More informationAnswers to Problem Set Number MIT (Fall 2005).
Answers to Problem Set Number 5. 18.305 MIT Fall 2005). D. Margetis and R. Rosales MIT, Math. Dept., Cambridge, MA 02139). November 23, 2005. Course TA: Nikos Savva, MIT, Dept. of Mathematics, Cambridge,
More informationSuppression of Spiral Waves and Spatiotemporal Chaos Under Local Self-adaptive Coupling Interactions
Commun. Theor. Phys. (Beijing, China) 45 (6) pp. 121 126 c International Academic Publishers Vol. 45, No. 1, January 15, 6 Suppression of Spiral Waves and Spatiotemporal Chaos Under Local Self-adaptive
More informationI. Explore tip trajectories with rotating anisotropy;
Below is the summary of the recent work. Updates for V1.1: I. Added irregular behavior tables for Nz = 31. II. Added some comments of the irregular behaviors of scroll waves. III. Added the corrected phase
More informationBasic Theory of Dynamical Systems
1 Basic Theory of Dynamical Systems Page 1 1.1 Introduction and Basic Examples Dynamical systems is concerned with both quantitative and qualitative properties of evolution equations, which are often ordinary
More informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More informationSynchronization Transitions in Complex Networks
Synchronization Transitions in Complex Networks Y. Moreno 1,2,3 1 Institute for Biocomputation and Physics of Complex Systems (BIFI) University of Zaragoza, Zaragoza 50018, Spain 2 Department of Theoretical
More informationAPPM 2360 Project 3 Mathematical Investigation of Cardiac Dynamics
APPM 2360 Project 3 Mathematical Investigation of Cardiac Dynamics Due: Thursday, December 6, 2018 by 4:59 p.m. Submit as a PDF to Assignments on Canvas 1 Introduction Cardiac Arrhythmia, or irregular
More informationLYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR
LYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR Peter H. Baxendale Department of Mathematics University of Southern California Los Angeles, CA 90089-3 USA baxendal@math.usc.edu
More informationLecture 5: Travelling Waves
Computational Biology Group (CoBI), D-BSSE, ETHZ Lecture 5: Travelling Waves Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015 26. Oktober 2016 2 / 68 Contents 1 Introduction to Travelling Waves
More informationBifurcation analysis of a normal form for excitable media: Are stable dynamical alternans on a ring possible?
Bifurcation analysis of a normal form for excitable media: Are stable dynamical alternans on a ring possible? Georg A. Gottwald School of Mathematics & Statistics, University of Sydney, NSW 2006, Australia.
More informationComputational simulation of the heart Edmund Altman 2010/2011
Computational simulation of the heart Edmund Altman 2010/2011 The candidate confirms that the work submitted is their own and the appropriate credit has been given where reference has been made to the
More informationPHYSICAL REVIEW E VOLUME 62, NUMBER 6. Target waves in the complex Ginzburg-Landau equation
PHYSICAL REVIEW E VOLUME 62, NUMBER 6 DECEMBER 2000 Target waves in the complex Ginzburg-Landau equation Matthew Hendrey, Keeyeol Nam, Parvez Guzdar,* and Edward Ott University of Maryland, Institute for
More informationPattern formation in Nikolaevskiy s equation
Stephen Cox School of Mathematical Sciences, University of Nottingham Differential Equations and Applications Seminar 2007 with Paul Matthews, Nottingham Outline What is Nikolaevskiy s equation? Outline
More informationKey words. cable equation, Neumann boundary condition, voltage conservation, finite difference, phase field, cardiac model, monodomain
Abstract. We present a finite difference method for modeling the propagation of electrical waves in cardiac tissue using the cable equation with homogeneous Neumann boundary conditions. This method is
More informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More informationAnisotropic electron distribution functions and the transition between the Weibel and the whistler instabilities
Anisotropic electron distribution functions and the transition between the Weibel and the whistler instabilities F. Pegoraro, L. Palodhi, F. Califano 5 th INTERNATIONAL CONFERENCE ON THE FRONTIERS OF PLASMA
More informationarxiv:chao-dyn/ v1 12 Feb 1996
Spiral Waves in Chaotic Systems Andrei Goryachev and Raymond Kapral Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, ON M5S 1A1, Canada arxiv:chao-dyn/96014v1 12
More informationLecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit.
Lecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit. Steve Cowley UCLA. This lecture is meant to introduce the simplest ideas about kink modes. It would take many lectures to develop the
More informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
More informationA note on discontinuous rate functions for the gate variables in mathematical models of cardiac cells
Procedia Computer Science (2) (22) 6 945 95 Procedia Computer Science www.elsevier.com/locate/procedia International Conference on Computational Science ICCS 2 A note on discontinuous rate functions for
More informationHow do nerve cells behave? At
Excitable Media: The Belousov-Zhabotinsky Chemical Reaction and the Heart Harold M. Hastings Professor and Chairperson Department of Physics Figure 2. Simulated cardiac dynamics by Flavio Fenton and Elizabeth
More informationPredator-swarm interactions
Predator-swarm interactions Hypothesis: swarming behaviour is an evolutionary adaptation that confers certain benefits on the individuals or group as a whole [Parrish,Edelstein-Keshet 1999; Sumpter 2010,
More information6.3.4 Action potential
I ion C m C m dφ dt Figure 6.8: Electrical circuit model of the cell membrane. Normally, cells are net negative inside the cell which results in a non-zero resting membrane potential. The membrane potential
More information3-dimensional simulation of long QT syndrome: early afterdepolarizations and reentry Ray Huffaker, Boris Kogan
3-dimensional simulation of long QT syndrome: early afterdepolarizations and reentry Ray Huffaker, Boris Kogan Abstract Early afterdepolarizations (EADs), thought to be highly arrhythmogenic in the case
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationLaurette TUCKERMAN Numerical Methods for Differential Equations in Physics
Laurette TUCKERMAN laurette@pmmh.espci.fr Numerical Methods for Differential Equations in Physics Time stepping: Steady state solving: 0 = F(U) t U = LU + N(U) 0 = LU + N(U) Newton s method 0 = F(U u)
More informationLaurette TUCKERMAN Rayleigh-Bénard Convection and Lorenz Model
Laurette TUCKERMAN laurette@pmmh.espci.fr Rayleigh-Bénard Convection and Lorenz Model Rayleigh-Bénard Convection Rayleigh-Bénard Convection Boussinesq Approximation Calculation and subtraction of the basic
More informationWave Turbulence and Condensation in an Optical Experiment
Wave Turbulence and Condensation in an Optical Experiment S. Residori, U. Bortolozzo Institut Non Linéaire de Nice, CNRS, France S. Nazarenko, J. Laurie Mathematics Institute, University of Warwick, UK
More informationTheoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions
Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions Tom Elsden 1 Andrew Wright 1 1 Dept Maths & Stats, University of St Andrews DAMTP Seminar - 8th May 2017 Outline Introduction Coordinates
More informationMath 575-Lecture 26. KdV equation. Derivation of KdV
Math 575-Lecture 26 KdV equation We look at the KdV equations and the so-called integrable systems. The KdV equation can be written as u t + 3 2 uu x + 1 6 u xxx = 0. The constants 3/2 and 1/6 are not
More informationBasic mechanisms of arrhythmogenesis and antiarrhythmia
EHRA EDUCATIONAL REVIEW AND PREPARATORY COURSE ON INVASIVE CARDIAC ELECTROPHYSIOLOGY EUROPEAN HEART HOUSE, February 2011 Basic mechanisms of arrhythmogenesis and antiarrhythmia Antonio Zaza Università
More informationSpatio-temporal patterns in two-dimensional excitable media subject to Robin boundary conditions
Applied Mathematics and Computation 146 (2003) 55 72 www.elsevier.com/locate/amc Spatio-temporal patterns in two-dimensional excitable media subject to Robin boundary conditions J.I. Ramos E.T.S. Ingenieros
More informationClassifying Mechanisms of Spiral Wave Breakup Underlying Cardiac Fibrillation Using Quantitative Metrics
Rochester Institute of Technology RIT Scholar Works Theses Thesis/Dissertation Collections 8--04 Classifying Mechanisms of Spiral Wave Breakup Underlying Cardiac Fibrillation Using Quantitative Metrics
More informationThick Brane World. Seyen Kouwn Korea Astronomy and Space Science Institute Korea
Thick Brane World Seyen Kouwn Korea Astronomy and Space Science Institute Korea Introduction - Multidimensional theory 1 Why are the physically observed dimensions of our Universe = 3 + 1 (space + time)?
More informationFrustrated drift of an anchored scroll-wave filament and the geodesic principle
Frustrated drift of an anchored scroll-wave filament and the geodesic principle Marcel Wellner, 1,2 Christian Zemlin, 1 and Arkady M. Pertsov 1 1 Department of Pharmacology, SUNY Upstate Medical University,
More informationDynamical mechanism of atrial fibrillation: a topological approach Christopher D. Marcotte 1 and Roman O. Grigoriev 2
Dynamical mechanism of atrial fibrillation: a topological approach Christopher D. Marcotte and Roman O. Grigoriev 2 ) EPSRC Centre for Predictive Modelling in Healthcare, University of Exeter, EX44QJ,
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.
More informationOrigins of Spiral Wave Meander and Breakup in a Two-Dimensional Cardiac Tissue Model
Annals of Biomedical Engineering, Vol. 28, pp. 755 771, 2000 Printed in the USA. All rights reserved. 0090-6964/2000/28 7 /755/17/$15.00 Copyright 2000 Biomedical Engineering Society Origins of Spiral
More informationDynamical Systems and Chaos Part II: Biology Applications. Lecture 10: Coupled Systems. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part II: Biology Applications Lecture 10: Coupled Systems. Ilya Potapov Mathematics Department, TUT Room TD325 Foreword In order to model populations of physical/biological
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More informationCardiac Electrophysiology on the Moving Heart
Cardiac Electrophysiology on the Moving Heart (... or a short story on shape analysis) Felipe Arrate FOCUS PROGRAM ON GEOMETRY, MECHANICS AND DYNAMICS The Legacy of Jerry Marsden The Fields Institute July
More informationStochastic excitation of streaky boundary layers. Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden
Stochastic excitation of streaky boundary layers Jérôme Hœpffner Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden Boundary layer excited by free-stream turbulence Fully turbulent inflow
More informationElectromagnetic Waves
Electromagnetic Waves Our discussion on dynamic electromagnetic field is incomplete. I H E An AC current induces a magnetic field, which is also AC and thus induces an AC electric field. H dl Edl J ds
More informationLecture 3 : Bifurcation Analysis
Lecture 3 : Bifurcation Analysis D. Sumpter & S.C. Nicolis October - December 2008 D. Sumpter & S.C. Nicolis General settings 4 basic bifurcations (as long as there is only one unstable mode!) steady state
More informationCollective and Stochastic Effects in Arrays of Submicron Oscillators
DYNAMICS DAYS: Long Beach, 2005 1 Collective and Stochastic Effects in Arrays of Submicron Oscillators Ron Lifshitz (Tel Aviv), Jeff Rogers (HRL, Malibu), Oleg Kogan (Caltech), Yaron Bromberg (Tel Aviv),
More informationA Model of Evolutionary Dynamics with Quasiperiodic Forcing
paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 2-5, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth
More informationPeriodic Behavior in Cardiac Tissue: Dynamics of Spatially Discordant Calcium Alternans
University of Colorado, Boulder CU Scholar Applied Mathematics Graduate Theses & Dissertations Applied Mathematics Spring 1-1-213 Periodic Behavior in Cardiac Tissue: Dynamics of Spatially Discordant Calcium
More informationSolutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the
More informationInfluence of velocity slip conditions on MHD peristaltic flow of a Prandtl fluid in a non-uniform channel
Malaysian Journal of Mathematical Sciences 11): 35 47 16) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal Influence of velocity slip conditions on MHD peristaltic
More informationNON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM
CYBERNETICS AND PHYSICS, VOL. 5, NO. 3, 016, 91 95 NON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM Leonid I. Manevitch Polymer and Composite Materials Department N. N. Semenov Institute
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationBifurcation and Stability Analysis of a Prey-predator System with a Reserved Area
ISSN 746-733, England, UK World Journal of Modelling and Simulation Vol. 8 ( No. 4, pp. 85-9 Bifurcation and Stability Analysis of a Prey-predator System with a Reserved Area Debasis Mukherjee Department
More informationA Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term
A Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term Peter Sternberg In collaboration with Dmitry Golovaty (Akron) and Raghav Venkatraman (Indiana) Department of Mathematics
More informationPhys 7221 Homework # 8
Phys 71 Homework # 8 Gabriela González November 15, 6 Derivation 5-6: Torque free symmetric top In a torque free, symmetric top, with I x = I y = I, the angular velocity vector ω in body coordinates with
More informationDiagonalization of the Coupled-Mode System.
Diagonalization of the Coupled-Mode System. Marina Chugunova joint work with Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaborators: Mason A. Porter, California Institute
More informationQuestion: Total. Points:
MATH 308 May 23, 2011 Final Exam Name: ID: Question: 1 2 3 4 5 6 7 8 9 Total Points: 0 20 20 20 20 20 20 20 20 160 Score: There are 9 problems on 9 pages in this exam (not counting the cover sheet). Make
More informationMathematical Biology. Criterion for stable reentry in a ring of cardiac tissue. John W. Cain
J. Math. Biol. DOI 10.1007/s00285-007-0100-z Mathematical Biology Criterion for stable reentry in a ring of cardiac tissue John W. Cain Received: 11 September 2006 / Revised: 29 March 2007 Springer-Verlag
More informationSOLUTIONS, PROBLEM SET 11
SOLUTIONS, PROBLEM SET 11 1 In this problem we investigate the Lagrangian formulation of dynamics in a rotating frame. Consider a frame of reference which we will consider to be inertial. Suppose that
More informationStability conditions for the traveling pulse: Modifying the restitution hypothesis
CHAOS VOLUME 1, NUMBER 3 SEPTEMBER 00 Stability conditions for the traveling pulse: Modifying the restitution hypothesis Eric Cytrynbaum a) Institute of Theoretical Dynamics, University of California,
More informationTwo Fluid Dynamo and Edge-Resonant m=0 Tearing Instability in Reversed Field Pinch
1 Two Fluid Dynamo and Edge-Resonant m= Tearing Instability in Reversed Field Pinch V.V. Mirnov 1), C.C.Hegna 1), S.C. Prager 1), C.R.Sovinec 1), and H.Tian 1) 1) The University of Wisconsin-Madison, Madison,
More informationLearning Cycle Linear Hybrid Automata for Excitable Cells
Learning Cycle Linear Hybrid Automata for Excitable Cells Sayan Mitra Joint work with Radu Grosu, Pei Ye, Emilia Entcheva, I V Ramakrishnan, and Scott Smolka HSCC 2007 Pisa, Italy Excitable Cells Outline
More informationNonlinear Free Vibration of Nanobeams Subjected to Magnetic Field Based on Nonlocal Elasticity Theory
Nonlinear Free Vibration of Nanobeams Subjected to Magnetic Field Based on Nonlocal Elasticity Theory Tai-Ping Chang 1 and Quey-Jen Yeh 1 Department of Construction Engineering, National Kaohsiung First
More informationNumerical solutions of the small dispersion limit of KdV, Whitham and Painlevé equations
Numerical solutions of the small dispersion limit of KdV, Whitham and Painlevé equations Tamara Grava (SISSA) joint work with Christian Klein (MPI Leipzig) Integrable Systems in Applied Mathematics Colmenarejo,
More information