Cardiac Electrophysiology on the Moving Heart
|
|
- Theresa Cain
- 5 years ago
- Views:
Transcription
1 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 12, 2012
2 2 Outline 1 Introduction 2 Model: Images + Electrophysiology 3 Electrophysiology 4 Image Model: Tracking shapes Tracking Fibers Shape matching Tracking the fibers 5 MeshFree Model
3 3 Cardiovascular diseases Heart diseases are today responsible for the 28% of deaths in western countries (Cancer is 30%). gure 2.1: Anatomical structure of the human heart. The illustration shows the heart in anterior-oblique view after having removed the anterior ventricular wall and the right rial appendage. Heart failure is responsible of 400,000 deaths per year in Europe fatigue insufficient ventricular contraction Sudden Death affects 1 in 10, 000 per year in developed countries electrical disorder ventricular fibrillation
4 3 Cardiovascular diseases Heart diseases are today responsible for the 28% of deaths in western countries (Cancer is 30%). gure 2.1: Anatomical structure of the human heart. The illustration shows the heart in anterior-oblique view after having removed the anterior ventricular wall and the right rial appendage. Individualized model for cardiac electrophysiology Diagnosis. Efficacy of defibrillation in infarcted hearts. Ablation location for correcting arrhythmias.
5 4 Outline 1 Introduction 2 Model: Images + Electrophysiology 3 Electrophysiology 4 Image Model: Tracking shapes Tracking Fibers Shape matching Tracking the fibers 5 MeshFree Model
6 5 Model: Images + Electrophysiology, but... and... Where is the mechanic model? The Heart Mechanics involves the detailed knowledge of elastic constants at different levels inside the myocardium. No-negligible importance of internal fibers.
7 5 Model: Images + Electrophysiology, but... and... Where is the mechanic model? The Heart Mechanics involves the detailed knowledge of elastic constants at different levels inside the myocardium. No-negligible importance of internal fibers. Can we infer more from a study mainly based on medical images? Is it possible to infer good/bad behaviour from motion itself?. Idea: Given a geometry, the wave is going to move like this, contracting the muscle like this, but is actually doing this... The solution requires less modeling.
8 6 Outline 1 Introduction 2 Model: Images + Electrophysiology 3 Electrophysiology 4 Image Model: Tracking shapes Tracking Fibers Shape matching Tracking the fibers 5 MeshFree Model
9 7 Maxwell s equations, but... Figure 2.3: Anatomical and fibrous structure of the human heart shown in a posterioranterior oblique view (from [64]). Note the piece of epicardium cut off the left ventricle, demonstrating the counterclockwise change of the myocardial fiber structure with increasing depth. 1 Electrical field strengths are not too high biological tissue is assumed to behave linear with regard to its electrical properties. 2 Cardiac electrical activity is reflected by low frequency components only ( several khz) derivatives with respect to time can be neglected ( quasi-static approximation of Maxwell s equations)
10 Maxwell: E = B t J = DE and... B t = 0 E = u Figure 2.3: Anatomical and fibrous structure of the human heart shown in a posterioranterior oblique view (from [64]). Note the piece of epicardium cut off the left ventricle, demonstrating the counterclockwise change of the myocardial fiber structure with increasing depth. Monodomain Model: 1) du dt = c 1f(u, w) + (D u) 2) dw dt = ɛ(u γw) 3)f(u, w)= c 1u(u α)(u 1) + c 2uw FitzHugh-Nagumo Model I ion=f(u, w) dw dt =ɛ(u γw) 7
11 8 Outline 1 Introduction 2 Model: Images + Electrophysiology 3 Electrophysiology 4 Image Model: Tracking shapes Tracking Fibers Shape matching Tracking the fibers 5 MeshFree Model
12 9 Outline 1 Introduction 2 Model: Images + Electrophysiology 3 Electrophysiology 4 Image Model: Tracking shapes Tracking Fibers Shape matching Tracking the fibers 5 MeshFree Model
13 10 The Method of coordinates: On Growth and Form (1917) In a very large part of morphology, our essential task lies in the comparison of related forms rather than in the precise definition of each; and the deformation of a complicated figure may be a phenomenon easy of comprehension, though the figure itself may have to be left unanalyzed and undefined....this method is the Method of Coordinates, on which is based the Theory of Transformations.
14 10 The Method of coordinates: On Growth and Form (1917) In a very large part of morphology, our essential task lies in the comparison of related forms rather than in the precise definition of each; and the deformation of a complicated figure may be a phenomenon easy of comprehension, though the figure itself may have to be left unanalyzed and undefined....this method is the Method of Coordinates, on which is based the Theory of Transformations.
15 10 The Method of coordinates: On Growth and Form (1917) In a very large part of morphology, our essential task lies in the comparison of related forms rather than in the precise definition of each; and the deformation of a complicated figure may be a phenomenon easy of comprehension, though the figure itself may have to be left unanalyzed and undefined....this method is the Method of Coordinates, on which is based the Theory of Transformations.
16 11 Our Method of coordinates Template Target E(ϕ T emp, T arg) = d V(id, ϕ) 2 }{{} Regulariz. + U(ϕ T emp, T arg) }{{} Dissimilarity
17 11 Our Method of coordinates Template Target or E(ϕ T emp, T arg) = d V(id, ϕ) 2 }{{} Regulariz. Ẽ(v) = 1 0 v 2 V dt + λu(ϕ(1)) + U(ϕ T emp, T arg) }{{} Dissimilarity
18 11 Our Method of coordinates Template Target... or ( µ(t) w) := ( ρ 0 [Dϕ(t, x)] 1 w ) x d 1 ( ) E(µ(t)) = µ(t) ( µ(t) K(ϕ(t, x), ϕ(t, y))ei) x e i dt + λu(ϕ(1)) ( dµ(t) dt dϕ(t, y) dt i=1 = 0 d i=1 ) d w = ( ) µ(t) K i (ϕ(t, x), ϕ(t, y)) e i x i=1 ( ( ) ) µ(t) µ(t) D 2K i (ϕ(t, z), ϕ(t, x)). w e i z x y
19 11 Our Method of coordinates Template Target For points... E(x, α) = 1 N α i(t) T α j(t) γ(x i(t), x j(t)) 2 i,j=1 dx s(t) dt dα s(t) dt = N γ (x k (t), x s(t)) α k (t) k=1 = N k=1 { } α s(t) T α k (t) 2γ (x k (t), x s(t))
20 12 Gradient Descent on the initial momentum Energy on the initial momentum: E(ρ 0) = ( ρ 0 Kρ 0) + λu(ϕ(1)) Variation on the initial momentum: ρ 0 ρ 0 + δρ 0 ( ) δu δe(ρ 0) = 2 ( δρ 0 Kρ 0) + λ δϕ (ϕ(1)) δϕ(1) (1) 1) From EPDiff... Linearized model: { } { } δϕ δϕ t = J δµ ϕ,µ δµ δϕ(0) = 0, δµ(0) = δρ 0 2) Adjoint system: ξ ϕ := δϕ, ξ µ := δµ 3) So (1) becomes { } ξϕ t = Jϕ,µ ξ µ { ξϕ ξ µ E(ρ 0) = 2ρ 0 + λk 1 ξ µ(0) } δe(ρ 0) = ( δρ 0 2Kρ 0 + λξ µ(0)) ξ ϕ(1) = δu (ϕ(1)), ξµ(1) = 0 δφ
21 13... Gradient Descent on the initial momentum
22 14 Outline 1 Introduction 2 Model: Images + Electrophysiology 3 Electrophysiology 4 Image Model: Tracking shapes Tracking Fibers Shape matching Tracking the fibers 5 MeshFree Model
23 15 Tracking the in-vivo fibers: Piecewise Affine + Diffeomorphic 1 Affine Registration: Translations + Rotations + Scale. Dissimilarity measure: Normalized correlation. Pipeline: Performs an affine registration, including scaling. Follows a PPD procedure to affine register the tensors. It makes easier (faster, more accurate registration) for lddmm to perform the nonlinear registration once the affine deformed template is much closer to the target. It allows the user to input the obtained deformed tensors into DtiStudio. 2 Nonlinear Registration: Diffeomorphic matching.
24 Introduction Model Electrophysiology Tracking MeshFree Tracking the in-vivo fibers: Piecewise Affine + Diffeomorphic 3 Results: Tensor Translation + Rotation + Nonlinear Deformation. 1 1 Vadakkumpadan, Trayanova, N., et al., Image-based models of cardiac structure in health and disease, Wiley Interdisciplinary Reviews: Systems Biology and Medicine, Vol. 2, (No Scale, Closest Neighbor Interpolation, Reorienting Vectors not complete tensors) 15
25 16 Outline 1 Introduction 2 Model: Images + Electrophysiology 3 Electrophysiology 4 Image Model: Tracking shapes Tracking Fibers Shape matching Tracking the fibers 5 MeshFree Model
26 17 Mesh or... Meshless Method? AND... Finite element methods have been extensively used for the spatial discretization of the myocardium. Complicated meshing procedures and element-based interpolation functions often result in algorithms which are either easy to implement, but numerically inaccurate, or accurate but labor-intensive The meshfree platform is more adaptive to different cardiac geometries and thus beneficial to individualized analysis.
27 18 Ω ϕ 3 ϕ(ω)
28 18 Ω ϕ 3 ϕ(ω)
29 18 Ω ϕ ϕ(ω) Particle Methods Complicated volume meshing procedures are excluded. No re-meshing is needed for improving spatial accuracy when deformation occurs.
30 19 Moving Least Squares (MLS) Approximation 1. {x 1(t),..., x N (t)} nodes (particles) in Ω R 3 2. p T (x) = [p 1(x),..., p m(x)] polynomial basis. Governing equation Given: locations x i (t), for i = 1,..., N, values u(x i, t), for i = 1,..., N Solve the associated ODE system { } d u = Φ(u, w) dt w and obtain new locations x i (t + t), for i = 1,..., N, new values u(x i, t + t), for i = 1,..., N Moving Least Squares Approximate the solution by: m u(x) = p k (x)a k (x) k=1 Minimizing the functional J = N i=1 ] 2 w(x x i ) [p T (x i )a(x) u i (w(x x i ) weighing function with compact support) Solve the MLS problem A(x)a(x) = B(x)u
31 MLS Approximation: Monodomain Model 1. Monodomain model u - membrane potential w - recovery variable. u t = c 1f(u, w) + (D u) w = ɛ(u γw) t f(u, w) = c 1 u(u α)(u 1) c 2 uw 2. Weak formulation φ - regular test function φ u Ω t = φ c 1 f(u, w) φ T (D u) Ω Ω φ v Ω t = φ ɛ(u γw) Ω 3. Meshfree approximation Φ = [φ 1 (x),..., φ N (x)] - shape function u Φu w Φw [ ] [ ] u Φ T Φ Ω t = Φ T Φ f(u, w) Ω [ ] Φ T D Φ u Ω [ ] [ ] w Φ T Φ Ω t = Φ T Φ ɛ(u γw) Ω 4. ODE system u t = f(u, w) + M 1 K u w = ɛ(u γw) t f(u, w) = c 1 u (u α) (u 1) c 2 u w where M = Φ T Φ K = Φ T D Φ Ω Ω 20
32 20 MLS Approximation: Monodomain Model 1. Monodomain model u - membrane potential w - recovery variable. u t = c 1f(u, w) + (D u) w = ɛ(u γw) t f(u, w) = c 1 u(u α)(u 1) c 2 uw 2. Weak formulation φ - regular test function φ u Ω t = φ c 1 f(u, w) φ T (D u) Ω Ω φ v Ω t = φ ɛ(u γw) Ω W R O N G!!!
33 20 MLS Approximation: Monodomain Model 1. Monodomain model u - membrane potential w - recovery variable. u t = c 1f(u, w) + (D u) w = ɛ(u γw) t f(u, w) = c 1 u(u α)(u 1) c 2 uw 3. Meshfree Approximation M u [ ] t + Φ T [(JΦ)ẋ] T u = Mf(u, w) Ku Ω [ M w t + ] Φ T [(JΦ)ẋ] T w = Mɛ(u γw) Ω 2. Weak formulation φ - regular test function φ u Ω t = φ c 1 f(u, w) φ T (D u) Ω Ω φ v Ω t = φ ɛ(u γw) Ω
34 21...Some preliminary results... fixed heart
35 22 Thanks!
Biomechanics. Soft Tissue Biomechanics
Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic
More information3D constitutive modeling of the passive heart wall
3D constitutive modeling of the passive heart wall Thomas S.E. Eriksson Institute of Biomechanics - Graz University of Technology Department of Biophysics - Medical University of Graz SFB status 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, July 27 2005 1 Stripes, Spots and Scrolls Indiana Unversity Physics REU Seminar,
More informationThe Significance of Reduced Conductance in the Sino-Atrial Node for Cardiac Rythmicity
The Significance of Reduced Conductance in the Sino-Atrial Node for Cardiac Rythmicity Will Nesse May 14, 2005 Sino-Atrial Node (SA node): The Pacemaker of the Heart 1 The SA node and the Atria: How the
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 informationSimulating Cardiac Electromechanics using Abaqus UEL
1 Simulating Cardiac Electromechanics using Abaqus UEL Introduction From a finite elements point of view, modeling the complex beating of heart tissues involves solving strongly coupled electromechanical
More informationMedical Image Analysis
Medical Image Analysis Instructor: Moo K. Chung mchung@stat.wisc.edu Lecture 3. Deformation-based Morphometry (DBM) January 30, 2007 Deformation based Morphometry (DBM) It uses deformation fields obtained
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 informationScroll 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 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 informationSegmentation by discrete watersheds Part 1: Watershed cuts
Segmentation by discrete watersheds Part 1: Watershed cuts Jean Cousty Four-Day Course on Mathematical Morphology in image analysis Bangalore 19- October 010 J. Serra, J. Cousty, B.S. Daya Sagar : Course
More informationArmin Toffler, The whole is more than the sum of the parts. Aristotle, BC. Prigogine I. and Stengers I. Order out of Chaos (1984)
1 OUTLINE 1. Introduction 2. Complexity 3. Waves in microstructured materials internal variables examples 4. Biophysics / biomechanics 5. Complexity around snow 6. Final remarks 2 One of the most highly
More informationWARHEAD FRAGMENTATION MODELING WITH PERIDYNAMICS
23 RD INTERNATIONAL SYMPOSIUM ON BALLISTICS TARRAGONA, SPAIN 16-20 APRIL 2007 WARHEAD FRAGMENTATION MODELING WITH PERIDYNAMICS Paul N. Demmie 1, Dale S. Preece 1 and Stewart A. Silling 1 1 Sandia National
More informationOblique derivative problems for elliptic and parabolic equations, Lecture II
of the for elliptic and parabolic equations, Lecture II Iowa State University July 22, 2011 of the 1 2 of the of the As a preliminary step in our further we now look at a special situation for elliptic.
More informationScroll 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 informationHervé DELINGETTE. INRIA Sophia-Antipolis
Hervé DELINGETTE INRIA Sophia-Antipolis French National Research Institute in Computer Science and Automatic Control Asclepios : 3D Segmentation, Simulation Platform, Soft Tissue Modeling http://www.inria.fr/asclepios
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 informationLectures on. Constitutive Modelling of Arteries. Ray Ogden
Lectures on Constitutive Modelling of Arteries Ray Ogden University of Aberdeen Xi an Jiaotong University April 2011 Overview of the Ingredients of Continuum Mechanics needed in Soft Tissue Biomechanics
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
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 informationMATH3203 Lecture 1 Mathematical Modelling and ODEs
MATH3203 Lecture 1 Mathematical Modelling and ODEs Dion Weatherley Earth Systems Science Computational Centre, University of Queensland February 27, 2006 Abstract Contents 1 Mathematical Modelling 2 1.1
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationCalculating Accurate Waveforms for LIGO and LISA Data Analysis
Calculating Accurate Waveforms for LIGO and LISA Data Analysis Lee Lindblom Theoretical Astrophysics, Caltech HEPL-KIPAC Seminar, Stanford 17 November 2009 Results from the Caltech/Cornell Numerical Relativity
More informationTract-Specific Analysis for DTI of Brain White Matter
Tract-Specific Analysis for DTI of Brain White Matter Paul Yushkevich, Hui Zhang, James Gee Penn Image Computing & Science Lab Department of Radiology University of Pennsylvania IPAM Summer School July
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 informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More informationNumerical simulations of cardiac arrhythmias and defibrillation
Numerical simulations of cardiac arrhythmias and defibrillation M. Hanslien, J. Sundnes and G.T. Lines Simula Research Laboratory and Department of Informatics, University of Oslo Abstract This paper is
More informationNeuroimage Processing
Neuroimage Processing Instructor: Moo K. Chung mkchung@wisc.edu Lecture 10-11. Deformation-based morphometry (DBM) Tensor-based morphometry (TBM) November 13, 2009 Image Registration Process of transforming
More informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
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 informationComputer Modeling in Bioengineering
Computer Modeling in Bioengineering Theoretical Background, Examples and Software Milos Kojic Harvard School of Public Health, USA University of Kragujevac, Serbia University of Texas Health Science Center
More informationMagnetoHemoDynamics in MRI devices
MagnetoHemoDynamics in MRI devices Conference on Mathematics of Medical Imaging June -4, Jean-Frédéric Gerbeau INRIA Paris-Rocquencourt & Laboratoire J-L. Lions France Joint work with A. Drochon (UTC),
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationSegmented medical images based simulations of Cardiac electrical activity and electrocardiogram: a model comparison
Segmented medical images based simulations of Cardiac electrical activity and electrocardiogram: a model comparison Charles Pierre, Olivier Rousseau, Yves Bourgault To cite this version: Charles Pierre,
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 informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationFast Geodesic Regression for Population-Based Image Analysis
Fast Geodesic Regression for Population-Based Image Analysis Yi Hong 1, Polina Golland 2, and Miaomiao Zhang 2 1 Computer Science Department, University of Georgia 2 Computer Science and Artificial Intelligence
More informationCIND Pre-Processing Pipeline For Diffusion Tensor Imaging. Overview
CIND Pre-Processing Pipeline For Diffusion Tensor Imaging Overview The preprocessing pipeline of the Center for Imaging of Neurodegenerative Diseases (CIND) prepares diffusion weighted images (DWI) and
More informationRitz Method. Introductory Course on Multiphysics Modelling
Ritz Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences Warsaw Poland
More informationExperiences with Model Reduction and Interpolation
Experiences with Model Reduction and Interpolation Paul Constantine Stanford University, Sandia National Laboratories Qiqi Wang (MIT) David Gleich (Purdue) Emory University March 7, 2012 Joe Ruthruff (SNL)
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
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 informationMathematical Models. Chapter Modelling the Body as a Volume Conductor
Chapter 2 Mathematical Models 2.1 Modelling the Body as a Volume Conductor As described in the previous chapter, the human body consists of billions of cells, which may be connected by various coupling
More informationCS 468. Differential Geometry for Computer Science. Lecture 17 Surface Deformation
CS 468 Differential Geometry for Computer Science Lecture 17 Surface Deformation Outline Fundamental theorem of surface geometry. Some terminology: embeddings, isometries, deformations. Curvature flows
More informationFETI Methods for the Simulation of Biological Tissues
SpezialForschungsBereich F 32 Karl Franzens Universita t Graz Technische Universita t Graz Medizinische Universita t Graz FETI Methods for the Simulation of Biological Tissues Ch. Augustin O. Steinbach
More informationarxiv: v1 [math.na] 1 May 2013
arxiv:3050089v [mathna] May 03 Approximation Properties of a Gradient Recovery Operator Using a Biorthogonal System Bishnu P Lamichhane and Adam McNeilly May, 03 Abstract A gradient recovery operator based
More information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More informationIsogeometric Analysis of the electrophysiology in the human heart: numerical simulation of the bidomain equations on the atria
MOX-Report No. 22/218 Isogeometric Analysis of the electrophysiology in the human heart: numerical simulation of the bidomain equations on the atria Pegolotti, L.; Dede', L.; Quarteroni, A. MOX, Dipartimento
More informationFluid Extensions for Optical Flow and Diffusion-based Image Registration
Fluid Extensions for Optical Flow and Diffusion-based Image Registration Jens-Peer Kuska 1, Patrick Scheibe 2, Ulf-Dietrich Braumann 2 1 Interdisciplinary Centre for Bioinformatics, University Leipzig,
More informationConservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body
Continuum Mechanics We ll stick with the Lagrangian viewpoint for now Let s look at a deformable object World space: points x in the object as we see it Object space (or rest pose): points p in some reference
More informationLecture 18: Bistable Fronts PHYS 221A, Spring 2017
Lecture 18: Bistable Fronts PHYS 221A, Spring 2017 Lectures: P. H. Diamond Notes: Xiang Fan June 15, 2017 1 Introduction In the previous lectures, we learned about Turing Patterns. Turing Instability is
More informationCilia-driven flow fields for transport and selective capture of bacteria
Cilia-driven flow fields for transport and selective capture of bacteria 0.1mm Janna Nawroth, Ph.D. Wyss Institute for Biologically Inspired Engineering, Harvard University Active Matter Kavli Institute
More informationFrom Diffusion Data to Bundle Analysis
From Diffusion Data to Bundle Analysis Gabriel Girard gabriel.girard@epfl.ch Computational Brain Connectivity Mapping Juan-les-Pins, France 20 November 2017 Gabriel Girard gabriel.girard@epfl.ch CoBCoM2017
More informationLearning Spectral Graph Segmentation
Learning Spectral Graph Segmentation AISTATS 2005 Timothée Cour Jianbo Shi Nicolas Gogin Computer and Information Science Department University of Pennsylvania Computer Science Ecole Polytechnique Graph-based
More informationNumerical Solutions of Geometric Partial Differential Equations. Adam Oberman McGill University
Numerical Solutions of Geometric Partial Differential Equations Adam Oberman McGill University Sample Equations and Schemes Fully Nonlinear Pucci Equation! "#+ "#* "#) "#( "#$ "#' "#& "#% "#! "!!!"#$ "
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 informationFinite volume method for CFD
Finite volume method for CFD Indo-German Winter Academy-2007 Ankit Khandelwal B-tech III year, Civil Engineering IIT Roorkee Course #2 (Numerical methods and simulation of engineering Problems) Mentor:
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationHearts ablaze! Radio Frequency Ablation as Treatment for Cardiac Arrhythmia
Hearts ablaze! Radio Frequency Ablation as Treatment for Cardiac Arrhythmia Kelly Gillen, Jason Greenberg, Sam Pevzner, Becky Rounds, Quinn Schneller Bee 453: Computer-Aided Engineering Professor Ashim
More informationQuestion 9: PDEs Given the function f(x, y), consider the problem: = f(x, y) 2 y2 for 0 < x < 1 and 0 < x < 1. x 2 u. u(x, 0) = u(x, 1) = 0 for 0 x 1
Question 9: PDEs Given the function f(x, y), consider the problem: 2 u x 2 u = f(x, y) 2 y2 for 0 < x < 1 and 0 < x < 1 u(x, 0) = u(x, 1) = 0 for 0 x 1 u(0, y) = u(1, y) = 0 for 0 y 1. a. Discuss how you
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More information2 A: The Shallow Water Equations
2 A: The Shallow Water Equations 2.1 Surface motions on shallow water Consider two-dimensional (x-z) motions on a nonrotating, shallow body of water, of uniform density, as shown in Fig. 1 below. The ow
More informationThe momentum map representation of images
The momentum map representation of images M. Bruveris 1, F. Gay-Balmaz, D. D. Holm 1, and T. S. Ratiu 3 Revised version August 1, 1 arxiv:91.99v [nlin.cd] 8 Apr 15 Abstract This paper discusses the mathematical
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationNetwork ABSTRACT. 1. Introduction. Yajuan Zhang. complex networks. representing. The method is. Atrial Fibrillation. constructs complex.
Engineering, 23, 5, 22-26 http://dx.doi.org/.4236/eng.23.5b5 Published Online October 23 (http://www.scirp.org/ournal/eng) Identification of Atrial Fibrillation Using Complex Network Similarity Yauan Zhang,2,
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More informationStudy of rotation of ellipsoidal particles in combined simple shear flow and magnetic fields
Study of rotation of ellipsoidal particles in combined simple shear flow and magnetic fields Jie Zhang 1, Cheng Wang 1 1. Department of Mechanical and Aerospace Engineering, Missouri University of Science
More informationMorphometry. John Ashburner. Wellcome Trust Centre for Neuroimaging, 12 Queen Square, London, UK.
Morphometry John Ashburner Wellcome Trust Centre for Neuroimaging, 12 Queen Square, London, UK. Morphometry is defined as: Measurement of the form of organisms or of their parts. The American Heritage
More informationThe Generalized Interpolation Material Point Method
Compaction of a foam microstructure The Generalized Interpolation Material Point Method Tungsten Particle Impacting sandstone The Material Point Method (MPM) 1. Lagrangian material points carry all state
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 informationNonlinear Equations for Finite-Amplitude Wave Propagation in Fiber-Reinforced Hyperelastic Media
Nonlinear Equations for Finite-Amplitude Wave Propagation in Fiber-Reinforced Hyperelastic Media Alexei F. Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada
More informationInterpolation in h-version finite element spaces
Interpolation in h-version finite element spaces Thomas Apel Institut für Mathematik und Bauinformatik Fakultät für Bauingenieur- und Vermessungswesen Universität der Bundeswehr München Chemnitzer Seminar
More informationCardiac Electrophysiology using LS-DYNA
Cardiac Electrophysiology using LS-DYNA Pierre L Eplattenier 1, Sarah Bateau-Meyer 1, Dave Benson 1, Vikas Kaul 2, Carl Schu 2 Mark Palmer 2, Darrell Swenson 2, Joshua Blauer 2 1 Livermore Software Technology
More information(TRAVELLING) 1D WAVES. 1. Transversal & Longitudinal Waves
(TRAVELLING) 1D WAVES 1. Transversal & Longitudinal Waves Objectives After studying this chapter you should be able to: Derive 1D wave equation for transversal and longitudinal Relate propagation speed
More informationIntroduction and preliminaries Wouter Zomervrucht, Februari 26, 2014
Introduction and preliminaries Wouter Zomervrucht, Februari 26, 204. Introduction Theorem. Serre duality). Let k be a field, X a smooth projective scheme over k of relative dimension n, and F a locally
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationNUMERICAL SIMULATION OF CARDIAC MUSCLES IN A RAT BIVENTRICULAR MODEL
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11 15 June 2018, Glasgow, UK NUMERICAL SIMULATION OF CARDIAC MUSCLES IN A RAT
More informationAn inverse problem for a system of steady-state reaction-diffusion equations on a porous domain using a collage-based approach
Journal of Physics: Conference Series PAPER OPEN ACCESS An inverse problem for a system of steady-state reaction-diffusion equations on a porous domain using a collage-based approach To cite this article:
More informationSchild s Ladder for the parallel transport of deformations in time series of images
Schild s Ladder for the parallel transport of deformations in time series of images Lorenzi Marco 1,2, Nicholas Ayache 1, Xavier Pennec 1, and the Alzheimer s Disease Neuroimaging Initiative 1 Project
More informationFast Needle Insertion to Minimize Tissue Deformation and Damage
29 IEEE International Conference on Robotics and Automation Kobe International Conference Center Kobe, Japan, May 2-7, 29 Fast Needle Insertion to Minimize Tissue Deformation and Damage Mohsen Mahvash
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 informationMECh300H Introduction to Finite Element Methods. Finite Element Analysis (F.E.A.) of 1-D Problems
MECh300H Introduction to Finite Element Methods Finite Element Analysis (F.E.A.) of -D Problems Historical Background Hrenikoff, 94 frame work method Courant, 943 piecewise polynomial interpolation Turner,
More informationLinear Algebra Review (Course Notes for Math 308H - Spring 2016)
Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationDynamics of beams. Modes and waves. D. Clouteau. September 16, Department of Mechanical and Civil Engineering Ecole Centrale Paris, France
Dynamics of and waves Department of Mechanical and Civil Engineering Ecole Centrale Paris, France September 16, 2008 Outline 1 2 Pre-stressed Bernoulli-Euler in moving frames Kinematical assumption Resultant
More informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More informationDetermination of pressure data in aortic valves
Determination of pressure data in aortic valves Helena Švihlová a, Jaroslav Hron a, Josef Málek a, K.R.Rajagopal b and Keshava Rajagopal c a, Mathematical Institute, Charles University, Czech Republic
More informationMultivariate Statistical Analysis of Deformation Momenta Relating Anatomical Shape to Neuropsychological Measures
Multivariate Statistical Analysis of Deformation Momenta Relating Anatomical Shape to Neuropsychological Measures Nikhil Singh, Tom Fletcher, Sam Preston, Linh Ha, J. Stephen Marron, Michael Wiener, and
More informationNetworks and Heterogeneous Media, Volume 1, Number 1, March 2006, pp
Manuscript submitted to AIMS Journals Volume X, Number X, XX 2X Website: http://aimsciences.org pp. X XX Networks and Heterogeneous Media, Volume 1, Number 1, March 26, pp. 185-218 ANALYSIS OF A CLASS
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 informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
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 informationPassive and Active Deformation Processes of 3D Fibre-Reinforced Caricatures of Cardiovascular Tissues
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Passive and Active Deformation Processes of 3D Fibre-Reinforced Caricatures of Cardiovascular Tissues Antonio DiCarlo 1, Paola Nardinocchi
More informationApplication of diffusion MRI to cancer, heart and brain connectome imaging
Colloquium @ Department of Physics, NTU Application of diffusion MRI to cancer, heart and brain connectome imaging March 11, 2014 Wen-Yih Isaac Tseng MD, PhD Advanced Biomedical MRI Lab Center for Optoelectronic
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationREDUCED ORDER METHODS
REDUCED ORDER METHODS Elisa SCHENONE Legato Team Computational Mechanics Elisa SCHENONE RUES Seminar April, 6th 5 team at and of Numerical simulation of biological flows Respiration modelling Blood flow
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture 1-20 September, 2017
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture 1-20 September, 2017 Institute of Structural Engineering Method of Finite Elements II 1 Course
More information