Strongly Scalable Parallel Simulations of High-resolution Models in Computational Cardiology
|
|
- Karin Johnson
- 6 years ago
- Views:
Transcription
1 Strongly Scalable Parallel Simulations of High-resolution Models in Computational Cardiology Christoph Augustin, Gernot Plank in coop. G. Haase, M. Liebmann, O. Steinbach, G. Holzapfel, A. Neic, A. Prassl, T. Fastl, T. Eriksson, A. Crozier Medical University of Graz SFB Mathematical Optimization and Applications in Biomedical Science Modeling and Simulations in Biomechanics, September 15 th, / 23
2 Outline Electromechanical Modeling Parallel Strategies Configuration Numerical Examples Open Tasks and Perspectives 1 / 23
3 Outline Electromechanical Modeling Parallel Strategies Configuration Numerical Examples Open Tasks and Perspectives 2 / 23
4 Mechanics of Cardiovascular Tissues collagen fiber Find the displacement u such that Div FS(u, x) = 0 u(x) = u D (x) FS(u, x)n(x) = t N (x) matrix sheet-normal axis for x Ω, for x ΓD, for x ΓN. muscle fiber n0 S = Sp (u, x) + Sa (Vm, η, u, x) F = I + Grad u the deformation gradient Sp the passive 2nd Piola-Kirchhoff stress tensor1,2 sheet axis s0 f0 fiber axis Sa the active 2nd Piola-Kirchhoff stress tensor2,3 u D the prescribed displacement t N the prescribed traction n the normal vector Vm the transmembrane voltage and η state variables 1 Holzapfel and Ogden Philos. Trans. R. Soc. Lond. Ser. A, pp Eriksson et al Mathematics and Mechanics of Solids, pp Smith et al Acta Numer., pp C. Augustin Modeling and Simulations in Biomechan 3 / 23
5 Passive Stress Model Constitutive equation using the free-energy function Ψ S p = 2 Ψ(C), Ψ(C) = U(F)+Ψiso(C)+Ψaniso(C), Ψ (locally) convex C nearly incompressible: penalty with κ e.g., U(J) = κ 2 (J 1)2, J = det F isotropic components: ground matrix, elastin e.g., Ψ iso(c) = c 2 (J 2/3 tr(c) 3) anisotropic components 4,5 : fibers, sheets Ψ aniso(c) = a { exp[b(j 2/3 I f 1) 2 ] 1 } 2b invariant I f = Ff 0 Ff 0: stretch in fiber direction 4 Fung American Journal of Physiology, pp Eriksson, Gasser, Holzapfel, Ogden, Stiffness ds/d λ, kpa 600 Unloading Loading Stress S (F/A), kpa 4 / 23
6 Electical Activation in the Myocardium The Bidomain equations 6,7 describe the spread of cardiac electrical activity find (V m, φ e, η) such that (σ i + σ e)c 1 φ e = σ ic 1 V m, σ ic 1 V m = σ ic 1 φ e + β I m, V m I m = C m + I ion(v m, η, u), t η = f(vm, η) t Figure : Bidomain representation of cardiac tissue in 2D 6 Tung PhD thesis, 7 Vigmond et al Prog. Biophys. Mol. Biol., pp V m = φ i φ e the transmembrane voltage φ i, φ e intra- and extracellular potential η a vector of state variables σ i, σ e conductivity tensors I m(v m, η) transmembrane current flow I i, I e, I ion current densities β surface to volume ratio of cardiac cells Simplification: in our experiments we replace C 1 by the identity matrix 5 / 23
7 Active Stress Models Relaxed (+) (+) (-) Myosin (-) (+) (+) Contracted Actin + ATP + Ca 2+ Z disk Active stress is generated by electrical activation in the myocardium S a = S a(v m, η, u) I s f (f 0 f 0) with f 0 the myocyte fiber orientation s = 1 mathematical, s = 1 mechanical choice 2 Cell models to compute scalar-values stress term S a Weakly coupled electromechanics e.g. NPStress 9 : S a = ε(v m)(k Sa V m S a) Strongly Coupled Electromechanics e.g. Rice 10 : h(s a, S a, V m, η, λ, λ) = 0 8 Ambrosi and Pezzuto J. Elast., pp Nash and Panfilov Progress in Biophysics and Molecular Biology, pp Rice et al Biophysical Journal, pp / 23
8 Outline Electromechanical Modeling Parallel Strategies Configuration Numerical Examples Open Tasks and Perspectives 7 / 23
9 Motivation High-resolution geometries meshes with up to N = O(10 9 ) degrees of freedom direct solvers: solving time O(N 2 ), very high memory consumption iterative solvers: solving time O(N), lower memory consumption we require strongly scalable parallel algorithms using iterative solvers 8 / 23
10 Why do we need HR geometries? Easy to motivate in electrics: very small features influence wave propagation resolutions below 100 µm would be worthwhile i.e. cellular level, O(10 9 ) cells Harder to do so in mechanics: use same mesh as in electics no data mapping or mesh coarsening needed some phenomena of interest involve small-scale features infarcts and ischemic regions multiple tissue layers and thin structures as papillary muscles, heart strings, valves parallel framework is available improves computational time for smaller meshes as well 9 / 23
11 Global Problem FEM + Newton yields linearized system: K (u k ) u = F K(u k ) u k+1 = u k + u. p Decomposition: A i K i(u k i )Ai u i i=1 solve with algebraic multigrid method (AMG 11,12 ) Ω 10 / 23
12 Algebraic Multigrid Method FEM + Newton yields linearized system: Ω 1 Ω 3 Ω 2 Ω 4 K (u k ) u = F K(u k ) u k+1 = u k + u. p Decomposition: A i K i(u k i )Ai u i i=1 solve with algebraic multigrid method (AMG 11,12 ) Ω 5 Ω6 Ω 7 fine grid restriction base level prolongation presmoothing postsmoothing 11 Plank et al Biomedical Engineering, IEEE Transactions on, pp Neic et al IEEE Trans. Biomed. Eng., pp / 23
13 Finite Element Tearing and Interconnecting Ω 1 Ω 2 Ω 3 Ω 4 FETI: Finite element tearing and interconnecting Tearing: K 1,k u 1 K 1,k.... =. K p,k u p K p,k Ω 5 Ω 6 Ω 7 generally u i u j on the interface Γ C 13 Farhat and Roux Int. J. Numer. Methods Engrg., pp Klawonn and Rheinbach ZAMM Z. Angew. Math. Mech., pp Augustin, Holzapfel, and Steinbach Int. J. Numer. Meth. Engrg., pp / 23
14 Finite Element Tearing and Interconnecting Ω 5 Ω 1 Ω 2 Ω 3 Ω 4 Ω 6 Ω 7 FETI: Finite element tearing and interconnecting Interconnecting: with Lagrange multipliers λ and boolean matrices B i K 1,k B 1 u 1 K 1,k K p,k B p u =. p K p,k B 1 B p 0 λ 0 K i,k P a generalized inverse this yields p B ik i,k B i λ = P i=1 p B ik i,k f i. i=1 13 Farhat and Roux Int. J. Numer. Methods Engrg., pp Klawonn and Rheinbach ZAMM Z. Angew. Math. Mech., pp Augustin, Holzapfel, and Steinbach Int. J. Numer. Meth. Engrg., pp / 23
15 Software FETI: C. Pechstein, C. A. Cardiac electromechanics: CARP: G. Plank, E. Vigmond, C. A. PETSc: Krylov Solver Package with BoomerAMG (hypre) general AMG for various problems Parallel Toolbox (PT): M. Liebmann, A. Neic, G. Haase special AMG for electrics and mechanics Tarantula, ParaView, ParMetis, PaStiX, MUMPS, FEAP TIME++ S = S p + S a(v m, η, u, f) Start electrics Solver mechanics u 12 / 23
16 Outline Electromechanical Modeling Parallel Strategies Configuration Numerical Examples Open Tasks and Perspectives 13 / 23
17 Geometry - MR Images 16 from the Visible Heart Lab ( 14 / 23
18 Geometry - Smoothing 17 with K. Bredis and M. Holler (KFU Graz) 15 / 23
19 Myocardium Model Tetrahedral or hybrid meshes 18 Parameters 19, fiber directions 20 Decomposition with ParMETIS Computations on SuperMUC Leibniz Rechenzentrum Munich Nr 12 in Top500 List - June cores 18 Prassl et al IEEE Trans. Biomed. Engineering, pp Eriksson et al Mathematics and Mechanics of Solids, pp Bayer et al Ann. Biomed. Eng., pp / 23
20 Circulatory System - Pressure BC A to B: loading B to C: isovolumetric contraction C to D: ventricular ejection 16 D C D to E: isovolumetric relaxation E to B: refilling LV pressure (kpa) 3 0 E Challenges: A LV volume (ml) MRI images usually taken at point B unloading of the geometry needed in isovolumetric phases the cavity volumes have to stay constant pressure volume realtions in ejection and filling phase B 17 / 23
21 Fixed Location Boundary Conditions fix base (not physiological) contact problem 21 use bath from bid. model soft elastic material apply D-BC to bath heart covered by double layered membrane (pericard) space between layers is filled with fluid attached to diaphragm and pleura 21 Fritz et al Biomech Model Mechanobiol, 18 / 23
22 Fixed Location Boundary Conditions fix base (not physiological) contact problem 21 use bath from bid. model soft elastic material apply D-BC to bath heart covered by double layered membrane (pericard) space between layers is filled with fluid attached to diaphragm and pleura 21 Fritz et al Biomech Model Mechanobiol, 18 / 23
23 Fixed Location Boundary Conditions fix vessel in- and outlets heart covered by double layered membrane (pericard) space between layers is filled with fluid attached to diaphragm and pleura 21 Fritz et al Biomech Model Mechanobiol, 18 / 23
24 Outline Electromechanical Modeling Parallel Strategies Configuration Numerical Examples Open Tasks and Perspectives 19 / 23
25 Scaling for Electromechanics - AMG 0000 Assembling PC Setup Solver Monodomain Total nodes DOF tets 100 timesteps 5 NS each Number of processes [-] 20 / 23
26 GPU 21 / 23
27 Outline Electromechanical Modeling Parallel Strategies Configuration Numerical Examples Open Tasks and Perspectives 22 / 23
28 Open Tasks and Perspectives Ongoing work finish GPU assembling for mechanics TH-elements, dynamic model, adaptive timestepping,... block system solvers and preconditioning projections between fine and coarse mesh fitting parameters to experiments (Cardioproof project) Wish list contact problems, hemodynamics and FSI 23 / 23
29 Ambrosi, D and S Pezzuto (2012). Active stress vs. active strain in mechanobiology: constitutive issues. In: J. Elast , pp Augustin, CM, GA Holzapfel, and O Steinbach (2014). Classical and All-floating FETI Methods for the Simulation of Arterial Tissues. In: Int. J. Numer. Meth. Engrg. 99.4, pp Augustin, CM and G Plank (2013). Simulating the mechanics of myoardial tissue using strongly scalable parallel algorithms. In: Biomed Tech (Berl) i. Augustin, CM and O Steinbach (2013). FETI Methods for the Simulation of Biological Tissues. In: Domain Decomposition Methods in Science and Engineering XX, Lecture Notes in Computational Science and Engineering 91, pp Bayer, JD et al. (2012). A novel rule-based algorithm for assigning myocardial fiber orientation to computational heart models. In: Ann. Biomed. Eng. 40(10), pp Bols, J et al. (2011). A computational method to assess initial stresses and unloaded configuration of patient-specific blood vessels. In: 5th International conference on Advanced COmputational Methods in ENgineering (ACOMEN 2011). Université de Liège. Braess, Dietrich and Regina Sarazin (1997). An efficient smoother for the Stokes problem. In: Applied Numerical Mathematics 23.1, pp Eriksson, TSE et al. (2013). Influence of myocardial fiber/sheet orientations on left ventricular mechanical contraction. In: Mathematics and Mechanics of Solids 18, pp Farhat, C and FX Roux (1991). A method of finite element tearing and interconnecting and its parallel solution algorithm. In: Int. J. Numer. Methods Engrg. 32, pp / 23
30 Fritz, T et al. (2013). Simulation of the contraction of the ventricles in a human heart model including atria and pericardium : Finite element analysis of a frictionless contact problem. In: Biomech Model Mechanobiol. Fung, YC (1967). Elasticity of soft tissues in simple elongation. In: American Journal of Physiology 213, pp Gurev, V et al. (2011). Models of cardiac electromechanics based on individual hearts imaging data. In: Biomechanics and modeling in mechanobiology 10.3, pp Holzapfel, GA, TC Gasser, and RW Ogden (2000). A new constitutive framework for arterial wall mechanics and a comperative study of material models. In: J. Elasticity 61, pp Holzapfel, GA and RW Ogden (2009). Constitutive modelling of passive myocardium: a structurally based framework for material characterization. In: Philos. Trans. R. Soc. Lond. Ser. A , pp Kerckhoffs, RCP et al. (2007). Coupling of a 3D finite element model of cardiac ventricular mechanics to lumped systems models of the systemic and pulmonic circulation. In: Annals of biomedical engineering 35.1, pp Klawonn, A and O Rheinbach (2010). Highly scalable parallel domain decomposition methods with an application to biomechanics. In: ZAMM Z. Angew. Math. Mech. 90.1, pp Nash, MP and AV Panfilov (2004). Electromechanical model of excitable tissue to study reentrant cardiac arrhythmias. In: Progress in Biophysics and Molecular Biology 85, pp Neic, A et al. (2012). Accelerating cardiac bidomain simulations using graphics processing units. In: IEEE Trans. Biomed. Eng. 59.8, pp Pathmanathan, P and JP Whiteley (2009). A numerical method for cardiac mechanoelectric simulations. In: Ann. Biomed. Eng. 37.5, pp / 23
31 Plank, G, AJ Prassl, and CM Augustin (2013). Computational Challenges in Building Multi-Scale and Multi-Physics Models of Cardiac Electro-Mechanics. In: Biomed Tech (Berl). Plank, G et al. (2007). Algebraic multigrid preconditioner for the cardiac bidomain model. In: Biomedical Engineering, IEEE Transactions on 54.4, pp Prassl, AJ et al. (2009). Automatically Generated, Anatomically Accurate Meshes for Cardiac Electrophysiology Problems. In: IEEE Trans. Biomed. Engineering 56.5, pp Rice, JJ et al. (2008). Approximate Model of Cooperative Activation and Crossbridge Cycling in Cardiac Muscle Using Ordinary Differential Equations". In: Biophysical Journal 95.5, pp Rumpel, T and K Schweizerhof (2003). Volume-dependent pressure loading and its influence on the stability of structures. In: International Journal for Numerical Methods in Engineering 56.2, pp Smith, NP et al. (2004). Multiscale computational modelling of the heart. In: Acta Numer. 13, pp Tung, L (1978). A bi-domain model for describing ischemic myocardial D-C potentials. PhD thesis. MIT, Cambridge, MA. Vanka, SP (1986). Block-implicit multigrid solution of Navier-Stokes equations in primitive variables. In: Journal of Computational Physics 65.1, pp Vigmond, EJ et al. (2007). Solvers for the cardiac bidomain equations. In: Prog. Biophys. Mol. Biol , pp Westerhof, N, GIJS Elzinga, and P Sipkema (1971). An artificial arterial system for pumping hearts. In: Journal of applied physiology 31.5, pp / 23
32 Westerhof, N, J-W Lankhaar, and BE Westerhof (2009). The arterial windkessel. In: Medical & biological engineering & computing 47.2, pp / 23
FETI 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 informationBiomechanics. 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 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 informationConstitutive Modeling of Biological Soft Tissues
Constitutive Modeling of Biological Soft Tissues Attila P. Nagy 1, David J. Benson 1, Vikas Kaul 2, Mark Palmer 2 1 Livermore Software Technology Corporation, Livermore, CA 94551, USA 2 Medtronic plc,
More informationPreprint Alexander Heinlein, Axel Klawonn, and Oliver Rheinbach Parallel Two-Level Overlapping Schwarz Methods in Fluid-Structure Interaction
Fakultät für Mathematik und Informatik Preprint 2015-15 Alexander Heinlein, Axel Klawonn, and Oliver Rheinbach Parallel Two-Level Overlapping Schwarz Methods in Fluid-Structure Interaction ISSN 1433-9307
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 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 informationAdaptation of a Rabbit Myocardium Material Model for Use in a Canine Left Ventricle Simulation Study
Adaptation of a Rabbit Myocardium Material Model for Use in a Canine Left Ventricle Simulation Study Matthew G. Doyle e-mail: mdoyl009@uottawa.ca Stavros Tavoularis 1 e-mail: stavros.tavoularis@uottawa.ca
More informationA Parallel Newton-Krylov-FETI-DP Solver Based on FEAP
A Parallel Newton-Krylov-FETI-DP Solver Based on FEAP Large-Scale Applications and Scalability for Problems in the Mechanics of Soft Biological Tissues in Arterial Wall Structures Andreas Fischle Forschungsseminar
More informationTechnische Universität Graz
Technische Universität Graz Classical and all floating FETI methods for the simulation of arterial tissues C. M. Augustin, G. A. Holzapfel, O. Steinbach Berichte aus dem Institut für Numerische Mathematik
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 informationSpace time finite element methods in biomedical applications
Space time finite element methods in biomedical applications Olaf Steinbach Institut für Angewandte Mathematik, TU Graz http://www.applied.math.tugraz.at SFB Mathematical Optimization and Applications
More informationFull-field measurements and identification for biological soft tissues: application to arteries in vitro
Centre for Health Engineering CNRS UMR 5146 INSERM IFR 143 Prof. Stéphane Avril Full-field measurements and identification for biological soft tissues: application to arteries in vitro using single-gage
More informationElasticity Model for Blood Vessel
Elasticity Model for Blood Vessel Violaine Louvet To cite this version: Violaine Louvet. Elasticity Model for Blood Vessel. study of different models of the mechanical non-linear stress-strain behaviour
More informationAN ANISOTROPIC PSEUDO-ELASTIC MODEL FOR THE MULLINS EFFECT IN ARTERIAL TISSUE
XI International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XI E. Oñate, D.R.J. Owen, D. Peric and B. Suárez (Eds) AN ANISOTROPIC PSEUDO-ELASTIC MODEL FOR THE MULLINS
More informationA Simultaneous Augmented Lagrange Approach for UNCORRECTED PROOF. the Simulation of Soft Biological Tissue 3
1 AQ1 this figure will be printed in b/w A Simultaneous Augmented Lagrange Approach for 2 the Simulation of Soft Biological Tissue 3 Dirk Böse 1, Sarah Brinkhues 2, Raimund Erbel 1, Axel Klawonn 3,Oliver
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 informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationShort title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic
Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone
More informationA Recursive Trust-Region Method for Non-Convex Constrained Minimization
A Recursive Trust-Region Method for Non-Convex Constrained Minimization Christian Groß 1 and Rolf Krause 1 Institute for Numerical Simulation, University of Bonn. {gross,krause}@ins.uni-bonn.de 1 Introduction
More informationMathematical Model of Blood Flow in Carotid Bifurcation
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Mathematical Model of Blood Flow in Carotid Bifurcation E. Muraca *,1, V. Gramigna 1, and G. Fragomeni 1 1 Department of Experimental Medicine
More informationCOMPUTATION OF EIGENSTRESSES IN THREE-DIMENSIONAL PATIENT-SPECIFIC ARTERIAL WALLS
IS Computation - Biomechanics of eigenstresses and Computational in three-dimensional Modeling of Living patient-specific Tissue arterial walls XIII International Conference on Computational Plasticity.
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 informationConstitutive modelling of passive myocardium A structurally-based framework for material characterization
ME338A - Final project - Paper review - due in class Thu, March 12, 2009, 11am Constitutive modelling of passive myocardium A structurally-based framework for material characterization Gerhard A. Holzapfel
More informationThe All-floating BETI Method: Numerical Results
The All-floating BETI Method: Numerical Results Günther Of Institute of Computational Mathematics, Graz University of Technology, Steyrergasse 30, A-8010 Graz, Austria, of@tugraz.at Summary. The all-floating
More informationA monolithic FEM solver for fluid structure
A monolithic FEM solver for fluid structure interaction p. 1/1 A monolithic FEM solver for fluid structure interaction Stefan Turek, Jaroslav Hron jaroslav.hron@mathematik.uni-dortmund.de Department of
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 informationSIMULATION OF MECHANICAL TESTS OF COMPOSITE MATERIAL USING ANISOTROPIC HYPERELASTIC CONSTITUTIVE MODELS
Engineering MECHANICS, Vol. 18, 2011, No. 1, p. 23 32 23 SIMULATION OF MECHANICAL TESTS OF COMPOSITE MATERIAL USING ANISOTROPIC HYPERELASTIC CONSTITUTIVE MODELS Tomáš Lasota*, JiříBurša* This paper deals
More informationAn Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03024-7 An Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors Michel Lesoinne
More informationBETI for acoustic and electromagnetic scattering
BETI for acoustic and electromagnetic scattering O. Steinbach, M. Windisch Institut für Numerische Mathematik Technische Universität Graz Oberwolfach 18. Februar 2010 FWF-Project: Data-sparse Boundary
More informationarxiv: v1 [math.na] 19 Dec 2013
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH HYPERELASTIC MODELS I: A PARTITIONED APPROACH arxiv:1312.5561v1 [math.na] 19 Dec 2013 ULRICH LANGER AND HUIDONG YANG Abstract. In this
More informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
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 informationIMA Preprint Series # 2155
ON THE ACTIVE RESPONSE OF SOFT LIVING TISSUES By P. Nardinocchi and L. Teresi IMA Preprint Series # 2155 ( January 2007 ) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY OF MINNESOTA 400 Lind
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 informationDeformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera
Karel October Tůma 24, Simulation 2018 of a bovine eye 1/19 Deformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera Karel Tůma 1 joint work
More informationSpectral element agglomerate AMGe
Spectral element agglomerate AMGe T. Chartier 1, R. Falgout 2, V. E. Henson 2, J. E. Jones 4, T. A. Manteuffel 3, S. F. McCormick 3, J. W. Ruge 3, and P. S. Vassilevski 2 1 Department of Mathematics, Davidson
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 informationScalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems
Scalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems Pierre Jolivet, F. Hecht, F. Nataf, C. Prud homme Laboratoire Jacques-Louis Lions Laboratoire Jean Kuntzmann INRIA Rocquencourt
More informationEffects of Aging on the Mechanical Behavior of Human Arteries in Large Deformations
International Academic Institute for Science and Technology International Academic Journal of Science and Engineering Vol. 3, No. 5, 2016, pp. 57-67. ISSN 2454-3896 International Academic Journal of Science
More informationarxiv:q-bio/ v1 [q-bio.qm] 30 Nov 2003
Material forces in the context of biotissue remodelling arxiv:q-bio/0312002v1 [q-bio.qm] 30 Nov 2003 K. Garikipati, H. Narayanan, E. M. Arruda, K. Grosh, S. Calve University of Michigan, Ann Arbor, Michigan
More informationPALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations
2013 SIAM Conference On Computational Science and Engineering Boston, 27 th February 2013 PALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations U. Villa,
More informationModeling. EC-Coupling and Contraction
Bioeng 6460 Electrophysiology and Bioelectricity Modeling of EC-Coupling and Contraction Frank B. Sachse fs@cvrti.utah.edu Overview Recapitulation Group Work Excitation-Contraction Coupling Hill s Muscle
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 informationNumerical modelling of shear-thinning non-newtonian flows in compliant vessels
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2007; 00:1 [Version: 2002/09/18 v1.01] Numerical modelling of shear-thinning non-newtonian flows in compliant vessels M.
More informationUser-Materials in ANSYS
User-Materials in ANSYS Holzapfel-Model l lfor Soft Tissues Prof. Dr.-Ing A. Fritsch Possibilities of user programming ANSYS User Programmable Features (UPF) are capabilities you can use to write your
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 informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationAMG for a Peta-scale Navier Stokes Code
AMG for a Peta-scale Navier Stokes Code James Lottes Argonne National Laboratory October 18, 2007 The Challenge Develop an AMG iterative method to solve Poisson 2 u = f discretized on highly irregular
More informationAlgebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
More informationFast solvers for steady incompressible flow
ICFD 25 p.1/21 Fast solvers for steady incompressible flow Andy Wathen Oxford University wathen@comlab.ox.ac.uk http://web.comlab.ox.ac.uk/~wathen/ Joint work with: Howard Elman (University of Maryland,
More informationCompressible Navier-Stokes (Euler) Solver based on Deal.II Library
Compressible Navier-Stokes (Euler) Solver based on Deal.II Library Lei Qiao Northwestern Polytechnical University Xi an, China Texas A&M University College Station, Texas Fifth deal.ii Users and Developers
More informationHigher order optimization and adaptive numerical solution for optimal control of monodomain equations in cardiac electrophysiology
Higher order optimization and adaptive numerical solution for optimal control of monodomain equations in cardiac electrophysiology Chamakuri Nagaiah and Karl Kunisch Institute of Mathematics and Scientific
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationA NONLINEAR ELIMINATION PRECONDITIONED INEXACT NEWTON METHOD FOR HETEROGENEOUS HYPERELASTICITY
A NONLINEAR ELIMINATION PRECONDITIONED INEXACT NEWTON METHOD FOR HETEROGENEOUS HYPERELASTICITY SHIHUA GONG AND XIAO-CHUAN CAI Abstract. We propose and study a nonlinear elimination preconditioned inexact
More informationOPENING ANGLE OF HUMAN SAPHENOUS VEIN
Opening angle of human saphenous vein XIII International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIII E. Oñate, D.R.J. Owen, D. Peric and M. Chiumenti (Eds) OPENING
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 informationFETI domain decomposition method to solution of contact problems with large displacements
FETI domain decomposition method to solution of contact problems with large displacements Vít Vondrák 1, Zdeněk Dostál 1, Jiří Dobiáš 2, and Svatopluk Pták 2 1 Dept. of Appl. Math., Technical University
More informationNonlinear Elasticity, Anisotropy, Material Stability and Residual stresses in Soft Tissue
Nonlinear Elasticity, Anisotropy, Material Stability and Residual stresses in Soft Tissue R.W. Ogden Department of Mathematics, University of Glasgow Glasgow G 8QW, UK Email: rwo@maths.gla.ac.uk Home Page:
More informationNumerical Simulation of Magnetic Drug Targeting with Flow Structural Interaction in an Arterial Branching Region of Interest
Excerpt from the Proceedings of the COMSOL Conference 2010 Paris Numerical Simulation of Magnetic Drug Targeting with Flow Structural Interaction in an Arterial Branching Region of Interest Alexandru M.
More informationSimulating ventricular elastance with a heart-arterial interaction model
Simulating ventricular elastance with a heart-arterial interaction model Anita Gerstenmayer 1, Bernhard Hametner 2, Stephanie Parragh 1,2, Thomas Weber 3, Siegfried Wassertheurer 2 1 Department for Analysis
More informationConstitutive models. Constitutive model: determines P in terms of deformation
Constitutive models Constitutive model: determines P in terms of deformation Elastic material: P depends only on current F Hyperelastic material: work is independent of path strain energy density function
More informationin this web service Cambridge University Press
CONTINUUM MECHANICS This is a modern textbook for courses in continuum mechanics. It provides both the theoretical framework and the numerical methods required to model the behavior of continuous materials.
More informationApplication of Preconditioned Coupled FETI/BETI Solvers to 2D Magnetic Field Problems
Application of Preconditioned Coupled FETI/BETI Solvers to 2D Magnetic Field Problems U. Langer A. Pohoaţǎ O. Steinbach 27th September 2004 Institute of Computational Mathematics Johannes Kepler University
More informationAuxiliary space multigrid method for elliptic problems with highly varying coefficients
Auxiliary space multigrid method for elliptic problems with highly varying coefficients Johannes Kraus 1 and Maria Lymbery 2 1 Introduction The robust preconditioning of linear systems of algebraic equations
More informationPreprint ISSN
Fakultät für Mathematik und Informatik Preprint 2015-14 Allison H. Baker, Axel Klawonn, Tzanio Kolev, Martin Lanser, Oliver Rheinbach, and Ulrike Meier Yang Scalability of Classical Algebraic Multigrid
More informationComparison of the deflated preconditioned conjugate gradient method and algebraic multigrid for composite materials
Comput Mech (2012) 50:321 333 DOI 10.1007/s00466-011-0661-y ORIGINAL PAPER Comparison of the deflated preconditioned conjugate gradient method and algebraic multigrid for composite materials T. B. Jönsthövel
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationA FACE-BASED SMOOTHED FINITE ELEMENT METHOD FOR HYPERELASTIC MODELS AND TISSUE GROWTH
th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationModelling the effect of gap junctions on tissue-level cardiac electrophysiology
Modelling the effect of gap junctions on tissue-level cardiac electrophysiology Doug Bruce Pras Pathmanathan Jonathan P. Whiteley douglas.bruce@oriel.ox.ac.uk pras.pathmanathan@cs.ox.ac.uk jonathan.whiteley@cs.ox.ac.uk
More informationNonlocal computational methods applied to composites structures
Nonlocal computational methods applied to composites structures Norbert Germain, Frédéric Feyel, Jacques Besson To cite this version: Norbert Germain, Frédéric Feyel, Jacques Besson. Nonlocal computational
More informationThermo-elastic Response of Cutaneous and Subcutaneous Tissues to Noninvasive Radiofrequency Heating
Thermo-elastic Response of Cutaneous and Subcutaneous Tissues to Noninvasive Radiofrequency Heating Joel N. Jiménez Lozano, Paulino Vacas-Jacques, Walfre Franco. Excerpt from the Proceedings of the 2012
More informationNew Multigrid Solver Advances in TOPS
New Multigrid Solver Advances in TOPS R D Falgout 1, J Brannick 2, M Brezina 2, T Manteuffel 2 and S McCormick 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, P.O.
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationA User Friendly Toolbox for Parallel PDE-Solvers
A User Friendly Toolbox for Parallel PDE-Solvers Gundolf Haase Institut for Mathematics and Scientific Computing Karl-Franzens University of Graz Manfred Liebmann Mathematics in Sciences Max-Planck-Institute
More informationAlgebraic Multigrid Methods for the Oseen Problem
Algebraic Multigrid Methods for the Oseen Problem Markus Wabro Joint work with: Walter Zulehner, Linz www.numa.uni-linz.ac.at This work has been supported by the Austrian Science Foundation Fonds zur Förderung
More informationOn the Numerical Modelling of Orthotropic Large Strain Elastoplasticity
63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More informationSpline-Based Hyperelasticity for Transversely Isotropic Incompressible Materials
Paper 260 Civil-Comp Press, 2012 Proceedings of the Eleventh International Conference on Computational Structures Technology, B.H.V. Topping, (Editor), Civil-Comp Press, Stirlingshire, Scotland Spline-Based
More informationDaniel Ch. Haspinger 1 Sae-Il Murtada 1,2 Justyna A. Niestrawska 1 Gerhard A. Holzapfel 1,3 1 INTRODUCTION
Received: 23 April 2018 Revised: 9 August 2018 Accepted: 5 September 2018 DOI: 10.1002/zamm.201800113 SPECIAL ISSUE PAPER Numerical analyses of the interrelation between extracellular smooth muscle orientation
More informationPropagation of dissection in a residually-stressed artery model
Biomech Model Mechanobiol (2017) 16:139 149 DOI 10.1007/s10237-016-0806-1 ORIGINAL PAPER Propagation of dissection in a residually-stressed artery model Lei Wang 1 Steven M. Roper 1 Nicholas A. Hill 1
More informationA new strain energy function for the hyperelastic modelling of ligaments and tendons
A new strain energy function for the hyperelastic modelling of ligaments and tendons University of Manchester BMC-BAMC 2015 Anterior cruciate ligament reconstruction surgery Ligament and tendon hierarchical
More informationANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow
ANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow Michel Kern Inria Paris Rocquencourt Maison de la Simulation C2S@Exa Days, Inria Paris Centre, Novembre 2016 M. Kern
More informationLoss of mass and performance in skeletal muscle tissue: a continuum model
Communications in Applied and Industrial Mathematics ISSN 2038-0909 Research Article Loss of mass and performance in skeletal muscle tissue: a continuum model Giulia Giantesio 1, Alfredo Marzocchi 1, Alessandro
More informationMultigrid and Iterative Strategies for Optimal Control Problems
Multigrid and Iterative Strategies for Optimal Control Problems John Pearson 1, Stefan Takacs 1 1 Mathematical Institute, 24 29 St. Giles, Oxford, OX1 3LB e-mail: john.pearson@worc.ox.ac.uk, takacs@maths.ox.ac.uk
More informationarxiv: v1 [math.na] 3 May 2018
; :1 25 Fluid-structure interaction simulations of venous valves: a monolithic ALE method for large structural displacements arxiv:1805.01321v1 [math.na] 3 May 2018 S. Calandrini 1 *, E. Aulisa 1 1 Department
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More informationPeristaltic Pump. Introduction. Model Definition
Peristaltic Pump Introduction In a peristaltic pump, rotating rollers are squeezing a flexible tube. As the pushed down rollers move along the tube, the fluid in the tube follows the motion. The main advantage
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
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 informationForce Traction Microscopy: an inverse problem with pointwise observations
Force Traction Microscopy: an inverse problem with pointwise observations Guido Vitale Dipartimento di Matematica Politecnico di Torino, Italia Cambridge, December 2011 Force Traction Microscopy: a biophysical
More informationCellular Electrophysiology and Biophysics
BIOEN 6003 Cellular Electrophysiology and Biophysics Modeling of Force Development in Myocytes II Frank B. Sachse, University of Utah Overview Experimental Studies Sliding Filament Theory Group work Excitation-Contraction
More information14 - hyperelastic materials. me338 - syllabus hyperelastic materials hyperelastic materials. inflation of a spherical rubber balloon
me338 - syllabus holzapfel nonlinear solid mechanics [2000], chapter 6, pages 205-305 1 2 isotropic hyperelastic materials inflation of a spherical rubber balloon cauchy stress invariants in terms of principal
More informationSchwarz-type methods and their application in geomechanics
Schwarz-type methods and their application in geomechanics R. Blaheta, O. Jakl, K. Krečmer, J. Starý Institute of Geonics AS CR, Ostrava, Czech Republic E-mail: stary@ugn.cas.cz PDEMAMIP, September 7-11,
More informationAn Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84
An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84 Introduction Almost all numerical methods for solving PDEs will at some point be reduced to solving A
More informationInternational Journal of Scientific & Engineering Research, Volume 4, Issue 7, July-2013 ISSN
648 Local sensitivity analysis of cardiovascular system parameters R. Gul and S. Bernhard, Fachbereich Mathematik, FU Berlin, Germany. Dept. of Electrical Engineering and Information Technology, Pforzheim
More informationAdaptive Analysis of Bifurcation Points of Shell Structures
First published in: Adaptive Analysis of Bifurcation Points of Shell Structures E. Ewert and K. Schweizerhof Institut für Mechanik, Universität Karlsruhe (TH), Kaiserstraße 12, D-76131 Karlsruhe, Germany
More information