A Model of the Blood Flow in the Retina of the Eye
|
|
- Eugenia Wade
- 5 years ago
- Views:
Transcription
1 A Model of the Blood Flow in the Retina of the Eye Edmond Rusjan with Andrea Dziubek, Giovanna Guidoboni, Alon Harris, Anil Hirani, William Thistleton SUNY Poly, IUPUI, University of Illinois at Urbana-Champaign Vanderbilt, April 14, 2017
2 Abstract: Understanding the blood flow in the retina of the eye may provide insights into several eye pathologies and ultimately lead to better treatments. We model the flow as a hierarchical Darcy flow on a curved surface and solve the model numerically using discrete exterior calculus and finite element exterior calculus. Results support the hypothesis that changes in the shape of the retina cause significant changes in the ocular blood flow, which may play a role in the dynamics of open angle glaucoma.
3 Outline Modeling the Eye: Goals and Motivation A Multiscale Model of Blood Flow in the Retina Hierarchical Darcy Flow Vascular Tree Architecture Computational Method Hierarchical Discretization Spatial Discretization Curved Surfaces Discrete Exterior Calculus (DEC) and Finite Element Exterior Calculus (FEEC) Results Ocular Curvature May Play a Role in Glaucoma Summary and Future Plans Modeling the Eye: Goals and Motivation 3
4 The Eye The eye is the organ of sight: the sensory part is the retina. Modeling the Eye: Goals and Motivation 4
5 Goals Advance the theory: most accurate models become theory Help reason: simplified models help us think more clearly Simulate experiments: experiments may be hard, or even impossible Competing goals: tension between simplicity and accuracy Solution: a collection of models at various levels of refinement Modeling the Eye: Goals and Motivation 5
6 Motivation Preserve the vision Window into the body Part of the brain Opportunity to study flow on curved surfaces Modeling the Eye: Goals and Motivation 6
7 Outline Modeling the Eye: Goals and Motivation A Multiscale Model of Blood Flow in the Retina Hierarchical Darcy Flow Vascular Tree Architecture Computational Method Hierarchical Discretization Spatial Discretization Curved Surfaces Discrete Exterior Calculus (DEC) and Finite Element Exterior Calculus (FEEC) Results Ocular Curvature May Play a Role in Glaucoma Summary and Future Plans A Multiscale Model of Blood Flow in the Retina 7
8 Fundus Images A healthy right eye (left) and a healthy left eye (right) A Multiscale Model of Blood Flow in the Retina 8
9 Regions of Particular Clinical Interest in the Retina optic nerve head (ONH) inferior quadrant (I) superior quadrant (S) nasal quadrant (N) temporal quadrant (T) fovea i At the level of capillaries and small arterioles and venules, we don t see the details of the vascular tree architecture The tissue appears to be a porous medium Think of the vascular tree as a hierarchical structure Parameterize these hierarchical levels by a continuous parameter ϑ A Multiscale Model of Blood Flow in the Retina 9
10 Hierarchical Darcy Flow Consider averaged (smeared) quantities at position x in some spatial domain Ω for each hierarchy level ϑ Velocity v(x, ϑ) Hierarchical velocity ω(x, ϑ) Pressure p(x, ϑ) Model each ϑ level by Darcy flow equations and couple vertically with hierarchical flow (n b v) + ϑ (n bω) = 0 n b v = K p n b ω = α p ϑ x Ω, ϑ (0,1) x Ω, ϑ (0,1) x Ω, ϑ (0,1) A Multiscale Model of Blood Flow in the Retina 10
11 Parameters and Boundary Conditions Parameters Tissue porosity n b (x, ϑ) Tissue perfusion n b ω Tissue permeability tensor K(x, ϑ) Hierarchical permeability α(x, ϑ) Hydraulic conductivities G v, G a Venous and arterial pressures p v(x) and p a(x) Boundary conditions n b ω(x, 0) = G v (p(x, 0) p v (x)) n b ω(x, 1) = G a (p a (x) p(x, 1)) n b v(x, ϑ) n = 0 x Ω x Ω x Ω, ϑ [0,1] Couple to feeding arteries and draining veins G v = α v δ v (x) and G a = α a δ a (x) α v and α a are venous and arterial conductances A Multiscale Model of Blood Flow in the Retina 11
12 Vascular Tree Architecture δ v (x) and δ a (x) are delta functions which differ from zero only where the feeding arteries and draining veins are located, respectively Need to extract feeding arteries and draining veins from fundus images Alternatively, model the vascular tree Young (1808): Compare parent and daughters vessels Assume two identical daughters r Radius ratio: p r d = = 1.26 A Area ratio: p A d = Murray (1926): Generalize to asymmetric trees Similar to Pythagora, except cubes: r 3 p = r 3 d Derives from optimality (variational principle) Forgotten for many years, then rediscovered A Multiscale Model of Blood Flow in the Retina 12
13 Optimal Transport = Murray s Law P f power needed to maintain the blood flow (viscous losses) P m power needed to metabolically maintain the blood and the vessel P t = P f + P m total power Minimize P t! Assume cylindrical vessels: f = cp f is volumetric flow rate p is pressure difference c is conductance coefficient Assume Poiseuille flow: c = πr4 8µl µ is blood viscosity l is vessel length Min P f = pf = af 2 r 4 = r? a = 8µl π Min P m = mv = mπr 2 l = br 2 = r 0? m is a metabolic coefficient b = mπl Min P t = dpt dr = 0 = 4af 2 r 5 + 2br = 0 = f = kr 3! k = ( b 2a )1/2 A Multiscale Model of Blood Flow in the Retina 13
14 Recent Experimental Data and Models Takahashi (2009, 2013): f = kr m Theory (Murray): m = 3 Data (Takahashi 2009): m 2.85 Data (Takahashi 2013): significant deviations in m between various organs In addition to radius, interested in length and angle Theory: open problem! Data: rather scarce - open problem! Model (Takahashi 2009): L = 7.4r 1.15 Open question: what optimality? A Multiscale Model of Blood Flow in the Retina 14
15 Outline Modeling the Eye: Goals and Motivation A Multiscale Model of Blood Flow in the Retina Hierarchical Darcy Flow Vascular Tree Architecture Computational Method Hierarchical Discretization Spatial Discretization Curved Surfaces Discrete Exterior Calculus (DEC) and Finite Element Exterior Calculus (FEEC) Results Ocular Curvature May Play a Role in Glaucoma Summary and Future Plans Computational Method 15
16 Hierarchical Discretization Standard FE procedure Rewrite in weak form Introduce piecewise linear basis functions {ϕ i } with i = 0,.., n on [0, 1] corresponding to nodes {θ i } and we let ϕ = ϕ i. We assume permeabilities and pressure are piecewise linear interpolations of the hierarchical variable n n K(x, θ) = K(x, θ k )ϕ k (θ), α(x, θ) = α(x, θ k )ϕ k (θ), n j=0 Ω k=0 p(x, θ) = k=0 n p(x, θ k )ϕ k (θ) k=0 ( ) K ij (x) p j (x), q(x) dx n j=0 Ω α ij (x) p j (x) q(x) dx Computational Method = f i (x) q(x) dx 16 Ω
17 Hierarchical Discretization Cont. where K ij (x) = α ij (x) = n K(x, θ k ) ϕ k (θ) ϕ i (θ) ϕ j (θ) dθ [0,1] n α(x, θ k ) ϕ k (θ) ϕ i(θ) ϕ j (θ) dθ (1) [0,1] θ θ k=0 k=0 f i (x) = α v p v δ i0 + α a p a δ in The hierarchical discretization leads to a tridiagonal system for the pressures p i (x, θ i ). Computational Method 17
18 Example: 3 layer model For three levels i = 0, 1, 2 and with isotropic and space independent permeability K ij = K ij I per level we have K ij = K ij This yields K 00 + α 00 K 01 + α 01 K 10 + α 10 K 11 + α 11 K 12 + α 12 p 0 p 1 = K 21 + α 21 K 22 + α 22 p 2 Permeability matrix K ij, conductivity matrix α ij, and right hand side vector f i as given in equation (1). f 0 f 1 f 2. (2) Computational Method 18
19 The K and α matrices in the 3 layer model The scalar K ij is the ij entry of the matrix 3K K 1 K 0 + K 1 0 K 0 + K 1 K 0 + 6K 1 + K 2 K 1 + K 2, 24 0 K 1 + K 2 K 1 + 3K 2 The scalar α ij is the ij entry of the matrix α 0 + α 1 + G v (α 0 + α 1 ) 0 (α 0 + α 1 ) α 0 + 2α 1 + α 2 (α 1 + α 2 ), 0 (α 1 + α 2 ) α 1 + α 2 + G a The scalar f i is the i component of the vector [ Gv p v, 0, G a p a ] T. Here α 0, α 1, α 2 and K 0, K 1, K 2 are the conductivities and permeabilities at the three hierarchical levels and are given in equation (1). Computational Method 19
20 Exterior Algebra Motivation: generalize the cross product from R 3 to R n V an n dimensional real vector space with elements u, v,... u v exterior (wedge) product of u and v anti-symmetric: u v = v u linear in each factor: (au + bv) w = au w + bv w 0 V = R 1 V = V 2 V = V V, 2-vectors continue recursively, demanding associativity p V p-vectors, dim( p V ) = ( ) n p dim( n V ) = 1 V inner product space: natural isomorphism between V and its dual V flat: : V V, u (v) = (u, v) sharp: : V V, = 1 inner product on V induces inner product on p V (the notion of length induces notions of area, volume, p-volume) Computational Method 20
21 Hodge Motivation: u v is a 2-vector, so not quite u v yet Solution: so map u v to u v - call this map Hodge Then: u v = (u v) V oriented inner product space: : p V (n p) V natural isomorphism Let {e 1, e 2,..., e n } be a basis for V Let σ = e 1 e 2 e n be a chosen orientation Hodge : u v = ( u, v)σ, v (n p) V Intuition: complementary n p vector to the given p vector orthogonal consistent with the orientation Side result: u v = (u, v) = (u v) = (v u) Computational Method 21
22 Hodge Example V = R 2 : V = R 3 : 1 = e 1 e 2 e 1 = e 2 e 2 = e 1 (e 1 e 2 ) = 1 1 = e 1 e 2 e 3 e 1 = e 2 e 3... =... (e 1 e 2 ) = e 3... =... (e 1 e 2 e 3 ) = 1 Computational Method 22
23 Exterior Calculus Differential forms - generalizations of differentials Let M be a manifold 0-forms are functions 1-forms are differentials k-forms at a point P are k vectors (elements in k TP M) Exterior Derivative d df = f dx i x i d(fdx 1 dx k ) = df dx 1 dx k Adjoint derivative: δ = d, with (δα, β) = (α, dβ) Laplacian: = (d + δ) 2 Interior Product: (i X α)(v 2,..., v k ) = α(x, v 2,..., v k ) Lie Derivative (using Cartan s Magic Formula): L X α = di X α + i X dα Computational Method 23
24 Curved Surfaces Numerically, a major issue is: retina is a curved surface Flat: ( x i ) = g ij dx j, (v ) i = g ij v j Sharp: (dx i ) = g ij x j Hodge : ( α) i1,i 2,...,i n k = 1 k! αj1,...,j k det g ɛ j1,...,j k,i 1,...,i n k Computational Method 24
25 Algebraic Topology Approach to Space Discretization Simplices Simplicial Complexes Chains Cochains Boundary operator Discrete exterior derivative D = T Discrete version of Stokes theorem Computational Method 25
26 M DEC : DEC Hodge * 1 Marsden, Hirani, Desbrun, et al.: σ σ α = 1 σ σ α [M DEC ] ij = σk i σ k i δ ij Advantage: M DEC diagonal Disadvantage: M DEC not positive definite, because circumcenter may be outside the simplex Example: standard 2-simplex; vertices: (0, 0), (1, 0), (0, 1) 1 M1 DEC = Computational Method 26
27 M W hit : FEEC Hodge * Arnold, Falk, Winther. Whitney 0-forms are barycentric coordinates: Whitney 1-forms: η k = λ i dλ j λ j dλ i [M W hit ] ij = η i η j da Advantage: M W hit is positive definite, because barycenter is always inside the simplex Disadvantage: M W hit not diagonal, (M W hit ) 1 not sparse Example: standard 2-simplex; vertices: (0, 0), (1, 0), (0, 1) M W hit 1 = Computational Method 27
28 PyDEC - A Python Implementation of DEC Authors: Hirani and Bell Efficient operator implementation in terms of sparse matrices Data structures Simplicial Complex Regular Cube Complex Operators Discrete Exterior Derivative Diagonal Sparse Matrix Discrete Hodge First Order Whitney Hodge Provides a nice playground for numerical experiments Computational Method 28
29 Data Figure: Retina Fundus Image webvision.med.utah.edu/book/part-i-foundations/simple-anatomy-of-the-retina/ Computational Method 29
30 Spatial Discretization Figure: Triangularization by gmsh. Computational Method 30
31 Example: 3 layer model geometry Figure: Three layer discretization The θ = 1 layer represents arterioles, the θ = 1 2 capillaries and θ = 0 represents venules. layer represents Computational Method 31
32 Outline Modeling the Eye: Goals and Motivation A Multiscale Model of Blood Flow in the Retina Hierarchical Darcy Flow Vascular Tree Architecture Computational Method Hierarchical Discretization Spatial Discretization Curved Surfaces Discrete Exterior Calculus (DEC) and Finite Element Exterior Calculus (FEEC) Results Ocular Curvature May Play a Role in Glaucoma Summary and Future Plans Results 32
33 Capillary Velocity Field Figure: Capillary Blood Velocity Results 33
34 Application: Myopia Figure: Capillary Blood Velocity in Myopia The myopic eye is extended along the optical axis. Results 34
35 Eye Deformations Figure: Eye Shapes Results 35
36 Velocity Field Deformations Figure: Velocity Difference Caution: to compare the velocity field difference in a meaningful way, all velocity fields are first pulled back to the sphere. Results 36
37 Outline Modeling the Eye: Goals and Motivation A Multiscale Model of Blood Flow in the Retina Hierarchical Darcy Flow Vascular Tree Architecture Computational Method Hierarchical Discretization Spatial Discretization Curved Surfaces Discrete Exterior Calculus (DEC) and Finite Element Exterior Calculus (FEEC) Results Ocular Curvature May Play a Role in Glaucoma Summary and Future Plans Summary and Future Plans 37
38 Summary The hierarchical Darcy equation models blood flow in the retina as a hierarchical porous medium. Exterior Calculus generalizes vector calculus to manifold and is convenient for formulating and solving models on surfaces. The discrete exterior derivative d is unique and determined by the discrete Stokes theorem. The discrete Hodge * operator is an open question and a topic of current research. DEC and FEEC are related; they differ in the Hodge * discretization. DEC and FEEC can be mixed with traditional FEs. Curvature may play an important role in the development of glaucoma. Summary and Future Plans 38
39 Future Plans Generalize the model to space dependent parameters. Couple the model with elasticity. Use the model to make predictions for other eye pathologies. The eye is the window into the body. What can we learn about conditions in the rest of the body? Personalized medicine: couple the model with image analysis front end. Summary and Future Plans 39
40 Acknowledgments NSF DMS and NSF CCF NIH 1R21EY A1 Slocum Dickson Foundation Research to Prevent Blindness (RPB, NY, USA) Indiana University Collaborative Research Grant of the Office of the Vice President for Research Chair Gutenberg funds SUNY STEM Summer Research SUNY Poly School of Arts and Sciences 40
41 References Abraham, Marsden, Ratiu, Manifolds, Tensor Analysis and Applications. Flanders, Differential Forms with Applications to the Physical Sciences. Frankel, The Geometry of Physics. T. Sherman, On Connecting Large Vessels to Small, The Meaning of Murray s Law, J. Gen. Phisiol, 78, , (1981). T. Takahashi et. al., A mathematical model for the distribution of hemodynamic parameters in the human retinal microvascular network, J. Biorrheol, 23, 77-86, (2009). Desbrun, A. N. Hirani, J. E. Marsden: Discrete Exterior Calculus for Variational Problems in Computer Graphics and Vision, Proceedings of the 42nd IEEE Conference on Decision and Control, Honolulu 5 (2003), Arnold D., Falk R., Winther R. : Finite Element Exterior Calculus: from Hodge Theory to Numerical Stability, Bull. Amer. Math. Soc. 47,2 (2010),
42 References Cont. Anil N. Hirani and Kaushik Kalyanaraman: Numerical Experiments for Darcy Flow on a Surface Using Mixed Exterior Calculus Methods, E-print arxiv: v1 [cs.na] on arxiv.org (2011). A. N. Hirani, K. B. Nakshatrala, J. H. Chaudhry: Numerical method for Darcy flow derived using Discrete Exterior Calculus, E-print arxiv: [math.na] on arxiv.org (2008) Nathan Bell and Anil N. Hirani: PyDEC: Software and Algorithms for Discretization of Exterior Calculus, ACM Transactions on Mathematical Software, Vol 39 1 (2012), A. Dziubek, E. Rusjan, W. Thistleton: Challenges in Blood Flow Simulation: Numerical Methods and Image Processing Tools, Proceedings of the ASME/FDA, 1st Annual Frontiers in Medical Devices, Washington DC (2013). A. Dziubek, G. Guidoboni, A. Harris, A.H. Hirani, E. Rusjan, W. Thistleton, Effect of ocular shape and vascular geometry on retinal hemodynamics: a computational model, Biomechanics and Modeling in Mechanobiology (2015),
43 Thank You! 43
New Model Stability Criteria for Mixed Finite Elements
New Model Stability Criteria for Mixed Finite Elements Andrew Gillette Department of Mathematics Institute of Computational Engineering and Sciences University of Texas at Austin http://www.math.utexas.edu/users/agillette
More informationFrom point cloud data to the continuum model of geometry
From point cloud data to the continuum model of geometry J. Harrison University of California, Berkeley July 22, 2007 Continuum model of geometry By the continuum model of geometry we refer to smooth manifolds
More informationA Generalization for Stable Mixed Finite Elements
A Generalization for Stable Mixed Finite Elements Andrew Gillette joint work with Chandrajit Bajaj Department of Mathematics Institute of Computational Engineering and Sciences University of Texas at Austin,
More informationDifferential Geometry: Discrete Exterior Calculus
Differential Geometry: Discrete Exterior Calculus [Discrete Exterior Calculus (Thesis). Hirani, 23] [Discrete Differential Forms for Computational Modeling. Desbrun et al., 25] [Build Your Own DEC at Home.
More informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More informationCompatible Discretization Schemes
Compatible Discretization Schemes Marc Gerritsma 1 1 Collaborators Artur Palha, Guido Oud, Jasper Kreeft & Mick Bouman Marc Gerritsma Burgers Course CFD I, 25-29 January 2010 1 / 81 What if... Why mimetic
More informationIn collaboration with Mathieu Desbrun, Anil N. Hirani, and Jerrold E. Marsden. Abstract
73 Chapter 3 Discrete Exterior Calculus In collaboration with Mathieu Desbrun, Anil N. Hirani, and Jerrold E. Marsden. Abstract We present a theory and applications of discrete exterior calculus on simplicial
More informationDISCRETE EXTERIOR CALCULUS
DISCRETE EXTERIOR CALCULUS MATHIEU DESBRUN, ANIL N. HIRANI, MELVIN LEOK, AND JERROLD E. MARSDEN Abstract. We present a theory and applications of discrete exterior calculus on simplicial complexes of arbitrary
More informationTowards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity
Towards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity Melvin Leok Mathematics, University of Michigan, Ann Arbor. Joint work with Mathieu Desbrun, Anil Hirani, and Jerrold
More informationA case study in Interaction cohomology Oliver Knill, 3/18/2016
A case study in Interaction cohomology Oliver Knill, 3/18/2016 Simplicial Cohomology Simplicial cohomology is defined by an exterior derivative df (x) = F (dx) on valuation forms F (x) on subgraphs x of
More informationA posteriori error estimates in FEEC for the de Rham complex
A posteriori error estimates in FEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094
More informationElectromagnetic Theory with Discrete Exterior Calculus
Progress In Electromagnetics Research, Vol. 159, 59 78, 2017 Electromagnetic Theory with Discrete Exterior Calculus Shu C. Chen 1 and Weng C. Chew 2, * Abstract A self-contained electromagnetic theory
More informationNotes on Discrete Exterior Calculus
Notes on Discrete Exterior Calculus Andrew Gillette August 27, 2009 These notes are aimed at mathematics graduate students familiar with introductory algebraic and differential topology. The notes were
More informationEffect of ocular shape and vascular geometry on retinal hemodynamics: a computational model
Biomech Model Mechanobiol DOI 10.1007/s10237-015-0731-8 ORIGINAL PAPER Effect of ocular shape and vascular geometry on retinal hemodynamics: a computational model Andrea Dziubek 1 Giovanna Guidoboni 2,3,4
More informationMath 225B: Differential Geometry, Final
Math 225B: Differential Geometry, Final Ian Coley March 5, 204 Problem Spring 20,. Show that if X is a smooth vector field on a (smooth) manifold of dimension n and if X p is nonzero for some point of
More informationLecture 5 - Lie Algebra Cohomology II
Lecture 5 - Lie Algebra Cohomology II January 28, 2013 1 Motivation: Left-invariant modules over a group Given a vector bundle F ξ G over G where G has a representation on F, a left G- action on ξ is a
More informationCS 468 (Spring 2013) Discrete Differential Geometry
CS 468 (Spring 2013) Discrete Differential Geometry Lecture 13 13 May 2013 Tensors and Exterior Calculus Lecturer: Adrian Butscher Scribe: Cassidy Saenz 1 Vectors, Dual Vectors and Tensors 1.1 Inner Product
More informationarxiv: v3 [math.na] 27 May 2009
GEOMETRIC COMPUTATIONAL ELECTRODYNAMICS WITH VARIATIONAL INTEGRATORS AND DISCRETE DIFFERENTIAL FORMS arxiv:0707.4470v3 [math.na] 7 May 009 ARI STERN, YIYING TONG, MATHIEU DESBRUN, AND JERROLD E. MARSDEN
More informationSection 2. Basic formulas and identities in Riemannian geometry
Section 2. Basic formulas and identities in Riemannian geometry Weimin Sheng and 1. Bianchi identities The first and second Bianchi identities are R ijkl + R iklj + R iljk = 0 R ijkl,m + R ijlm,k + R ijmk,l
More informationVectors. January 13, 2013
Vectors January 13, 2013 The simplest tensors are scalars, which are the measurable quantities of a theory, left invariant by symmetry transformations. By far the most common non-scalars are the vectors,
More informationHemodynamics II. Aslı AYKAÇ, PhD. NEU Faculty of Medicine Department of Biophysics
Hemodynamics II Aslı AYKAÇ, PhD. NEU Faculty of Medicine Department of Biophysics Laplace s Law Relates the pressure difference across a closed elastic membrane on liquid film to the tension in the membrane
More informationMATH 4030 Differential Geometry Lecture Notes Part 4 last revised on December 4, Elementary tensor calculus
MATH 4030 Differential Geometry Lecture Notes Part 4 last revised on December 4, 205 Elementary tensor calculus We will study in this section some basic multilinear algebra and operations on tensors. Let
More informationDr. Allen Back. Dec. 3, 2014
Dr. Allen Back Dec. 3, 2014 forms are sums of wedge products of the basis 1-forms dx, dy, and dz. They are kinds of tensors generalizing ordinary scalar functions and vector fields. They have a skew-symmetry
More informationDifferential Geometry II Lecture 1: Introduction and Motivation
Differential Geometry II Lecture 1: Introduction and Motivation Robert Haslhofer 1 Content of this lecture This course is on Riemannian geometry and the geometry of submanifol. The goal of this first lecture
More informationSpace-Time Finite-Element Exterior Calculus and Variational Discretizations of Gauge Field Theories
21st International Symposium on Mathematical Theory of Networks and Systems July 7-11, 2014. Space-Time Finite-Element Exterior Calculus and Variational Discretizations of Gauge Field Theories Joe Salamon
More informationMATH Linear Algebra
MATH 304 - Linear Algebra In the previous note we learned an important algorithm to produce orthogonal sequences of vectors called the Gramm-Schmidt orthogonalization process. Gramm-Schmidt orthogonalization
More informationNOTES ON DIFFERENTIAL FORMS. PART 3: TENSORS
NOTES ON DIFFERENTIAL FORMS. PART 3: TENSORS 1. What is a tensor? Let V be a finite-dimensional vector space. 1 It could be R n, it could be the tangent space to a manifold at a point, or it could just
More informationDiscrete Exterior Calculus
Discrete Exterior Calculus and the Averaged Template Matching Equations Anil N. Hirani Jet Propulsion Laboratory California Institute of Technology hirani at caltech dot edu http://www.cs.caltech.edu/
More informationInvariant differential operators on the sphere
Third CoE Talk at the University of Tokyo p. 1/21 Invariant differential operators on the sphere Michael Eastwood Australian National University Third CoE Talk at the University of Tokyo p. 2/21 The round
More informationVariational Integrators for Maxwell s Equations with Sources
PIERS ONLINE, VOL. 4, NO. 7, 2008 711 Variational Integrators for Maxwell s Equations with Sources A. Stern 1, Y. Tong 1, 2, M. Desbrun 1, and J. E. Marsden 1 1 California Institute of Technology, USA
More informationThe Hodge Star Operator
The Hodge Star Operator Rich Schwartz April 22, 2015 1 Basic Definitions We ll start out by defining the Hodge star operator as a map from k (R n ) to n k (R n ). Here k (R n ) denotes the vector space
More informationSurvey on exterior algebra and differential forms
Survey on exterior algebra and differential forms Daniel Grieser 16. Mai 2013 Inhaltsverzeichnis 1 Exterior algebra for a vector space 1 1.1 Alternating forms, wedge and interior product.....................
More informationJustin Solomon MIT, Spring 2017
Fisher et al. Design of Tangent Vector Fields (SIGGRAPH 2007) Justin Solomon MIT, Spring 2017 Original version from Stanford CS 468, spring 2013 (Butscher & Solomon) Divergence Theorem Green s Theorem
More informationCALCULUS ON MANIFOLDS. 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M =
CALCULUS ON MANIFOLDS 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M = a M T am, called the tangent bundle, is itself a smooth manifold, dim T M = 2n. Example 1.
More informationComplex manifolds, Kahler metrics, differential and harmonic forms
Complex manifolds, Kahler metrics, differential and harmonic forms Cattani June 16, 2010 1 Lecture 1 Definition 1.1 (Complex Manifold). A complex manifold is a manifold with coordinates holomorphic on
More informationThe symplectic structure on moduli space (in memory of Andreas Floer)
The symplectic structure on moduli space (in memory of Andreas Floer) Alan Weinstein Department of Mathematics University of California Berkeley, CA 94720 USA (alanw@math.berkeley.edu) 1 Introduction The
More informationAposteriorierrorestimatesinFEEC for the de Rham complex
AposteriorierrorestimatesinFEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094 and a
More informationMathematical Model. M. Umar Qureshi, Mitchel J. Colebank, and Mette S. Olufsen
Mathematical Model M. Umar Qureshi, Mitchel J. Colebank, and Mette S. Olufsen Department of Mathematics, North Carolina State University, Raleigh, North Carolina 27695 Friday 7 th September, 2018 The 1D
More informationArterial Macrocirculatory Hemodynamics
Arterial Macrocirculatory Hemodynamics 莊漢聲助理教授 Prof. Han Sheng Chuang 9/20/2012 1 Arterial Macrocirculatory Hemodynamics Terminology: Hemodynamics, meaning literally "blood movement" is the study of blood
More informationIs a dual mesh really necessary?
Is a dual mesh really necessary? Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at ICIAM, Vancouver, 2011. Joint work with Ari Stern, UCSD. 19 July 2011
More informationVector fields and differential forms. William G. Faris
Vector fields and differential forms William G. Faris September 25, 2008 ii Contents 1 Forms 1 1.1 The dual space............................ 1 1.2 Differential 1-forms.......................... 2 1.3
More informationDISCRETE DIFFERENTIAL GEOMETRY: AN APPLIED INTRODUCTION Keenan Crane CMU /858B Fall 2017
DISCRETE DIFFERENTIAL GEOMETRY: AN APPLIED INTRODUCTION Keenan Crane CMU 15-458/858B Fall 2017 LECTURE 4: DIFFERENTIAL FORMS IN R n DISCRETE DIFFERENTIAL GEOMETRY: AN APPLIED INTRODUCTION Keenan Crane
More informationFOUNDATIONS OF COMPUTATIONAL GEOMETRIC MECHANICS
FOUNDATIONS OF COMPUTATIONAL GEOMETRIC MECHANICS MELVIN LEOK Received the SIAM Student Paper Prize and the Leslie Fox Prize in Numerical Analysis (second prize). Abstract. Geometric mechanics involves
More informationAlgebraic Multigrid for k-form Laplacians
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2; :1 6 [Version: 22/9/18 v1.2] Algebraic Multigrid for k-form Laplacians Nathan Bell and Luke N. Olson Siebel Center for Computer
More informationDifferential Geometry and Lie Groups with Applications to Medical Imaging, Computer Vision and Geometric Modeling CIS610, Spring 2008
Differential Geometry and Lie Groups with Applications to Medical Imaging, Computer Vision and Geometric Modeling CIS610, Spring 2008 Jean Gallier Department of Computer and Information Science University
More informationChapter 1. Preliminaries. The purpose of this chapter is to provide some basic background information. Linear Space. Hilbert Space.
Chapter 1 Preliminaries The purpose of this chapter is to provide some basic background information. Linear Space Hilbert Space Basic Principles 1 2 Preliminaries Linear Space The notion of linear space
More informationDiscrete Differential Geometry. Discrete Laplace-Beltrami operator and discrete conformal mappings.
Discrete differential geometry. Discrete Laplace-Beltrami operator and discrete conformal mappings. Technische Universität Berlin Geometric Methods in Classical and Quantum Lattice Systems, Caputh, September
More informationPublished in Computer Methods in Applied Mechanics and Engineering 198 (2009)
Published in Computer Methods in Applied Mechanics and Engineering 198 (2009) 1660 1672 GEOMETRIC DECOMPOSITIONS AND LOCAL BASES FOR SPACES OF FINITE ELEMENT DIFFERENTIAL FORMS DOUGLAS N. ARNOLD, RICHARD
More information4 Riemannian geometry
Classnotes for Introduction to Differential Geometry. Matthias Kawski. April 18, 2003 83 4 Riemannian geometry 4.1 Introduction A natural first step towards a general concept of curvature is to develop
More informationThe Chain Collocation Method: A Spectrally Accurate Calculus of Forms
The Chain Collocation Method: A Spectrally Accurate Calculus of Forms Dzhelil Rufat, Gemma Mason, Patrick Mullen, Mathieu Desbrun Applied Geometry Lab, Caltech Abstract Preserving in the discrete realm
More informationGEOMETRICAL TOOLS FOR PDES ON A SURFACE WITH ROTATING SHALLOW WATER EQUATIONS ON A SPHERE
GEOETRICAL TOOLS FOR PDES ON A SURFACE WITH ROTATING SHALLOW WATER EQUATIONS ON A SPHERE BIN CHENG Abstract. This is an excerpt from my paper with A. ahalov [1]. PDE theories with Riemannian geometry are
More informationHandout to Wu Characteristic Harvard math table talk Oliver Knill, 3/8/2016
Handout to Wu Characteristic Harvard math table talk Oliver Knill, 3/8/2016 Definitions Let G = (V, E) be a finite simple graph with vertex set V and edge set E. The f-vector v(g) = (v 0 (G), v 1 (G),...,
More informationMONDAY - TALK 4 ALGEBRAIC STRUCTURE ON COHOMOLOGY
MONDAY - TALK 4 ALGEBRAIC STRUCTURE ON COHOMOLOGY Contents 1. Cohomology 1 2. The ring structure and cup product 2 2.1. Idea and example 2 3. Tensor product of Chain complexes 2 4. Kunneth formula and
More informationPolitecnico di Milano School of Industrial and Information Engineering Master of Science in Mathematical Engineering
Politecnico di Milano School of Industrial and Information Engineering Master of Science in Mathematical Engineering PACS Project Coupled problems on one-dimensional networks with application to microcirculation
More informationMath 868 Final Exam. Part 1. Complete 5 of the following 7 sentences to make a precise definition (5 points each). Y (φ t ) Y lim
SOLUTIONS Dec 13, 218 Math 868 Final Exam In this exam, all manifolds, maps, vector fields, etc. are smooth. Part 1. Complete 5 of the following 7 sentences to make a precise definition (5 points each).
More informationINTRODUCTION TO ALGEBRAIC GEOMETRY
INTRODUCTION TO ALGEBRAIC GEOMETRY WEI-PING LI 1 Preliminary of Calculus on Manifolds 11 Tangent Vectors What are tangent vectors we encounter in Calculus? (1) Given a parametrised curve α(t) = ( x(t),
More informationSolution: We can cut the 2-simplex in two, perform the identification and then stitch it back up. The best way to see this is with the picture:
Samuel Lee Algebraic Topology Homework #6 May 11, 2016 Problem 1: ( 2.1: #1). What familiar space is the quotient -complex of a 2-simplex [v 0, v 1, v 2 ] obtained by identifying the edges [v 0, v 1 ]
More informationSymplectic and Poisson Manifolds
Symplectic and Poisson Manifolds Harry Smith In this survey we look at the basic definitions relating to symplectic manifolds and Poisson manifolds and consider different examples of these. We go on to
More informationMath 114: Course Summary
Math 114: Course Summary Rich Schwartz September 25, 2009 General Information: Math 114 is a course in real analysis. It is the second half of the undergraduate series in real analysis, M113-4. In M113,
More informationBrief Review of Vector Algebra
APPENDIX Brief Review of Vector Algebra A.0 Introduction Vector algebra is used extensively in computational mechanics. The student must thus understand the concepts associated with this subject. The current
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationFFTs in Graphics and Vision. The Laplace Operator
FFTs in Graphics and Vision The Laplace Operator 1 Outline Math Stuff Symmetric/Hermitian Matrices Lagrange Multipliers Diagonalizing Symmetric Matrices The Laplacian Operator 2 Linear Operators Definition:
More informationAn abstract Hodge Dirac operator and its stable discretization
An abstract Hodge Dirac operator and its stable discretization Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at ICCA, Tartu,. Joint work with Ari Stern,
More informationNONLOCAL EXTERIOR CALCULUS ON RIEMANNIAN MANIFOLDS
The Pennsylvania State University The Graduate School Department of Mathematics NONLOCAL EXTERIOR CALCULUS ON RIEMANNIAN MANIFOLDS A Dissertation in Mathematics by Thinh Duc Le c 2013 Thinh Duc Le Submitted
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems Professor Dr E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8, 2014 1 / 56 We study some basic
More informationCHARACTERISTIC CLASSES
1 CHARACTERISTIC CLASSES Andrew Ranicki Index theory seminar 14th February, 2011 2 The Index Theorem identifies Introduction analytic index = topological index for a differential operator on a compact
More informationAnalysis II: The Implicit and Inverse Function Theorems
Analysis II: The Implicit and Inverse Function Theorems Jesse Ratzkin November 17, 2009 Let f : R n R m be C 1. When is the zero set Z = {x R n : f(x) = 0} the graph of another function? When is Z nicely
More information22 : Hilbert Space Embeddings of Distributions
10-708: Probabilistic Graphical Models 10-708, Spring 2014 22 : Hilbert Space Embeddings of Distributions Lecturer: Eric P. Xing Scribes: Sujay Kumar Jauhar and Zhiguang Huo 1 Introduction and Motivation
More informationA simplicial gauge theory on spacetime
A simplicial gauge theory on spacetime Tore G. Halvorsen, NTNU, Trondheim Norway. Joint work with Snorre H. Christiansen, UiO. 1 Abstract In this talk I will introduce the variational form of the SU(N)
More informationINSTANTON MODULI AND COMPACTIFICATION MATTHEW MAHOWALD
INSTANTON MODULI AND COMPACTIFICATION MATTHEW MAHOWALD () Instanton (definition) (2) ADHM construction (3) Compactification. Instantons.. Notation. Throughout this talk, we will use the following notation:
More informationarxiv: v2 [math.ds] 26 Mar 2010
STRUCTURE-PRESERVING DISCRETIZATION OF INCOPRESSIBLE FLUIDS arxiv:0912.3989v2 [math.ds] 26 ar 2010 D. PAVLOV, P. ULLEN, Y. TONG, E. KANSO, J.E. ARSDEN,. DESBRUN CALTECH SU USC Abstract. The geometric nature
More informationDifferential Geometry qualifying exam 562 January 2019 Show all your work for full credit
Differential Geometry qualifying exam 562 January 2019 Show all your work for full credit 1. (a) Show that the set M R 3 defined by the equation (1 z 2 )(x 2 + y 2 ) = 1 is a smooth submanifold of R 3.
More informationLet V, W be two finite-dimensional vector spaces over R. We are going to define a new vector space V W with two properties:
5 Tensor products We have so far encountered vector fields and the derivatives of smooth functions as analytical objects on manifolds. These are examples of a general class of objects called tensors which
More informationDiscriminative Direction for Kernel Classifiers
Discriminative Direction for Kernel Classifiers Polina Golland Artificial Intelligence Lab Massachusetts Institute of Technology Cambridge, MA 02139 polina@ai.mit.edu Abstract In many scientific and engineering
More informationHamiltonian Toric Manifolds
Hamiltonian Toric Manifolds JWR (following Guillemin) August 26, 2001 1 Notation Throughout T is a torus, T C is its complexification, V = L(T ) is its Lie algebra, and Λ V is the kernel of the exponential
More informationBarometer Fluid rises until pressure at A, due its weight, equals atmospheric pressure at B. Unit: mm Hg (millimeters that mercury rises)
FLUID MECHANICS The study of the properties of fluids resulting from the action forces. Fluid a liquid, gas, or plasma We will only consider incompressible fluids i.e. liquids Pressure P F A (normal force)
More informationVortex Equations on Riemannian Surfaces
Vortex Equations on Riemannian Surfaces Amanda Hood, Khoa Nguyen, Joseph Shao Advisor: Chris Woodward Vortex Equations on Riemannian Surfaces p.1/36 Introduction Yang Mills equations: originated from electromagnetism
More informationarxiv: v1 [math.gr] 8 Nov 2008
SUBSPACES OF 7 7 SKEW-SYMMETRIC MATRICES RELATED TO THE GROUP G 2 arxiv:0811.1298v1 [math.gr] 8 Nov 2008 ROD GOW Abstract. Let K be a field of characteristic different from 2 and let C be an octonion algebra
More informationSome Elliptic and Subelliptic Complexes from Geometry
Geometric and Nonlinear Partial Differential Equations, Xi An Jiaotong University p. 1/14 Some Elliptic and Subelliptic Complexes from Geometry Michael Eastwood [ based on joint work with Robert Bryant,
More informationOn the singular elements of a semisimple Lie algebra and the generalized Amitsur-Levitski Theorem
On the singular elements of a semisimple Lie algebra and the generalized Amitsur-Levitski Theorem Bertram Kostant, MIT Conference on Representations of Reductive Groups Salt Lake City, Utah July 10, 2013
More informationA classification of equivariant gerbe connections
A classification of equivariant gerbe connections Byungdo Park (KIAS) joint work with Corbett Redden (LIU Post) Topology in Australia and South Korea IBS Center for Geometry and Physics 24.04.2018 Outline
More informationHigh Order Differential Form-Based Elements for the Computation of Electromagnetic Field
1472 IEEE TRANSACTIONS ON MAGNETICS, VOL 36, NO 4, JULY 2000 High Order Differential Form-Based Elements for the Computation of Electromagnetic Field Z Ren, Senior Member, IEEE, and N Ida, Senior Member,
More informationMath 797W Homework 4
Math 797W Homework 4 Paul Hacking December 5, 2016 We work over an algebraically closed field k. (1) Let F be a sheaf of abelian groups on a topological space X, and p X a point. Recall the definition
More information14 Higher order forms; divergence theorem
Tel Aviv University, 2013/14 Analysis-III,IV 221 14 Higher order forms; divergence theorem 14a Forms of order three................ 221 14b Divergence theorem in three dimensions.... 225 14c Order four,
More informationLECTURE 8: THE SECTIONAL AND RICCI CURVATURES
LECTURE 8: THE SECTIONAL AND RICCI CURVATURES 1. The Sectional Curvature We start with some simple linear algebra. As usual we denote by ( V ) the set of 4-tensors that is anti-symmetric with respect to
More informationAnalysis-3 lecture schemes
Analysis-3 lecture schemes (with Homeworks) 1 Csörgő István November, 2015 1 A jegyzet az ELTE Informatikai Kar 2015. évi Jegyzetpályázatának támogatásával készült Contents 1. Lesson 1 4 1.1. The Space
More informationarxiv: v2 [math.na] 8 Sep 2015
OMPLEXES OF DISRETE DISTRIBUTIONAL DIFFERENTIAL FORMS AND THEIR HOMOLOGY THEORY MARTIN WERNER LIHT arxiv:1410.6354v2 [math.na] 8 Sep 2015 Abstract. omplexes of discrete distributional differential forms
More informationDIFFERENTIAL FORMS AND COHOMOLOGY
DIFFERENIAL FORMS AND COHOMOLOGY ONY PERKINS Goals 1. Differential forms We want to be able to integrate (holomorphic functions) on manifolds. Obtain a version of Stokes heorem - a generalization of the
More informationarxiv: v1 [math.dg] 26 Apr 2014
WHITNEY FORMS FOR SPHERICAL TRIANGLES I: THE EULER, CAGNOLI, AND TUYNMAN AREA FORMULAS, BARYCENTRIC COORDINATES, AND CONSTRUCTION WITH THE EXTERIOR CALCULUS arxiv:1404.659v1 [math.dg] 6 Apr 014 DAVID W.
More informationGravitation: Tensor Calculus
An Introduction to General Relativity Center for Relativistic Astrophysics School of Physics Georgia Institute of Technology Notes based on textbook: Spacetime and Geometry by S.M. Carroll Spring 2013
More informationFinite element exterior calculus: A new approach to the stability of finite elements
Finite element exterior calculus: A new approach to the stability of finite elements Douglas N. Arnold joint with R. Falk and R. Winther Institute for Mathematics and its Applications University of Minnesota
More informationDifferential Geometry MTG 6257 Spring 2018 Problem Set 4 Due-date: Wednesday, 4/25/18
Differential Geometry MTG 6257 Spring 2018 Problem Set 4 Due-date: Wednesday, 4/25/18 Required problems (to be handed in): 2bc, 3, 5c, 5d(i). In doing any of these problems, you may assume the results
More informationA P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS
A P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS SHANGYOU ZHANG DEDICATED TO PROFESSOR PETER MONK ON THE OCCASION OF HIS 6TH BIRTHDAY Abstract. On triangular grids, the continuous
More informationComputation of Maxwell s equations on Manifold using DEC
arxiv:0908.4448v6 [math.na] 9 Dec 009 Computation of Maxwell s equations on Manifold using DEC Zheng Xie Yujie Ma. Center of Mathematical Sciences, Zhejiang University (3007),China. Key Laboratory of Mathematics
More information7 Curvature of a connection
[under construction] 7 Curvature of a connection 7.1 Theorema Egregium Consider the derivation equations for a hypersurface in R n+1. We are mostly interested in the case n = 2, but shall start from the
More informationMinimal Element Interpolation in Functions of High-Dimension
Minimal Element Interpolation in Functions of High-Dimension J. D. Jakeman, A. Narayan, D. Xiu Department of Mathematics Purdue University West Lafayette, Indiana Random Model as a function We can consider
More informationA POROUS MODEL OF TUMOR IN HYPERTHERMIA THERAPY WITH CONVECTION OF BLOOD FLOW
ISTP-16, 005, PRAGUE 16 TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA A POROUS MODEL OF TUMOR IN HYPERTHERMIA THERAPY WITH CONVECTION OF BLOOD FLOW Ping Yuan Department of Mechanical Engineering Lee-Ming
More informationGorenstein rings through face rings of manifolds.
Gorenstein rings through face rings of manifolds. Isabella Novik Department of Mathematics, Box 354350 University of Washington, Seattle, WA 98195-4350, USA, novik@math.washington.edu Ed Swartz Department
More informationNormed & Inner Product Vector Spaces
Normed & Inner Product Vector Spaces ECE 174 Introduction to Linear & Nonlinear Optimization Ken Kreutz-Delgado ECE Department, UC San Diego Ken Kreutz-Delgado (UC San Diego) ECE 174 Fall 2016 1 / 27 Normed
More informationCOMPUTING THE POISSON COHOMOLOGY OF A B-POISSON MANIFOLD
COMPUTING THE POISSON COHOMOLOGY OF A B-POISSON MANIFOLD MELINDA LANIUS 1. introduction Because Poisson cohomology is quite challenging to compute, there are only very select cases where the answer is
More information