Dr. Martin R. Pfaller is a postdoc in the group of Prof. Alison L. Marsden's Cardiovascular Biomechanics Computation Lab. He graduated with a PhD in Mechanical Engineering from the Technical University of Munich in 2019, working with Prof. Wolfgang A. Wall. His PhD research enhanced solid mechanics models of the heart by studying the interaction between the myocardium and the pericardium. He demonstrated how projection-based model order reduction could be used to speed up model personalization from patient data, such as cine MRI or blood pressure measurements. His current work focuses on physics-based and data-driven computational models for cardiovascular fluid and solid dynamics. Applications include congenital and acquired cardiovascular diseases, specifically the growth and remodeling of cardiac and cardiovascular tissue.
Honors & Awards
NIH Pathway to Independence Award (K99/R00), NHLBI (2022-2027)
Travel Award, Additional Ventures (2022)
Travel Award, Cardiovascular Institute (CVI) (2022)
Dissertation Award, Association of German Engineers (VDI) (2019)
Winner Science Slam, Technical University of Munich (2017)
Departmental Teaching Award, Technical University of Munich (2017)
Departmental Teaching Award, Technical University of Munich (2016)
Exchange Scholarship, German Academic Exchange Service (DAAD) (2013)
Exchange Scholarship, Prof. Dr.-Ing. Erich Müller-Stiftung (2013)
PhD, Technical University of Munich, Mechanical Engineering (2019)
MSc, Technical University of Munich, Mechanical Engineering (2013)
BSc, Technical University of Munich, Mechanical Engineering (2012)
Alison Marsden, Postdoctoral Faculty Sponsor
Automated generation of 0D and 1D reduced-order models of patient-specific blood flow.
International journal for numerical methods in biomedical engineering
Three-dimensional (3D) cardiovascular fluid dynamics simulations typically require hours to days of computing time on a high-performance computing cluster. One-dimensional (1D) and lumped-parameter zero-dimensional (0D) models show great promise for accurately predicting blood bulk flow and pressure waveforms with only a fraction of the cost. They can also accelerate uncertainty quantification, optimization, and design parameterization studies. Despite several prior studies generating 1D and 0D models and comparing them to 3D solutions, these were typically limited to either 1D or 0D and a singular category of vascular anatomies. This work proposes a fully automated and openly available framework to generate and simulate 1D and 0D models from 3D patient-specific geometries, automatically detecting vessel junctions and stenosis segments. Our only input is the 3D geometry; we do not use any prior knowledge from 3D simulations. All computational tools presented in this work are implemented in the open-source software platform SimVascular. We demonstrate the reduced-order approximation quality against rigid-wall 3D solutions in a comprehensive comparison with N=72 publicly available models from various anatomies, vessel types, and disease conditions. Relative average approximation errors of flows and pressures typically ranged from 1% to 10% for both 1D and 0D models, measured at the outlets of terminal vessel branches. In general, 0D model errors were only slightly higher than 1D model errors despite requiring only a third of the 1D runtime. Automatically generated ROMs can significantly speed up model development and shift the computational load from high-performance machines to personal computers. This article is protected by copyright. All rights reserved.
View details for DOI 10.1002/cnm.3639
View details for PubMedID 35875875
Model order reduction of flow based on a modular geometrical approximation of blood vessels.
Computer methods in applied mechanics and engineering
We are interested in a reduced order method for the efficient simulation of blood flow in arteries. The blood dynamics is modeled by means of the incompressible Navier-Stokes equations. Our algorithm is based on an approximated domain-decomposition of the target geometry into a number of subdomains obtained from the parametrized deformation of geometrical building blocks (e.g., straight tubes and model bifurcations). On each of these building blocks, we build a set of spectral functions by Proper Orthogonal Decomposition of a large number of snapshots of finite element solutions (offline phase). The global solution of the Navier-Stokes equations on a target geometry is then found by coupling linear combinations of these local basis functions by means of spectral Lagrange multipliers (online phase). Being that the number of reduced degrees of freedom is considerably smaller than their finite element counterpart, this approach allows us to significantly decrease the size of the linear system to be solved in each iteration of the Newton-Raphson algorithm. We achieve large speedups with respect to the full order simulation (in our numerical experiments, the gain is at least of one order of magnitude and grows inversely with respect to the reduced basis size), whilst still retaining satisfactory accuracy for most cardiovascular simulations.
View details for DOI 10.1016/j.cma.2021.113762
View details for PubMedID 34176992
On the Periodicity of Cardiovascular Fluid Dynamics Simulations.
Annals of biomedical engineering
Three-dimensional cardiovascular fluid dynamics simulations typically require computation of several cardiac cycles before they reach a periodic solution, rendering them computationally expensive. Furthermore, there is currently no standardized method to determine whether a simulation has yet reached that periodic state. In this work, we propose the use of an asymptotic error measurement to quantify the difference between simulation results and their ideal periodic state using open-loop lumped-parameter modeling. We further show that initial conditions are crucial in reducing computational time and develop an automated framework to generate appropriate initial conditions from a one-dimensional model of blood flow. We demonstrate the performance of our initialization method using six patient-specific models from the Vascular Model Repository. In our examples, our initialization protocol achieves periodic convergence within one or two cardiac cycles, leading to a significant reduction in computational cost compared to standard methods. All computational tools used in this work are implemented in the open-source software platform SimVascular. Automatically generated initial conditions have the potential to significantly reduce computation time in cardiovascular fluid dynamics simulations.
View details for DOI 10.1007/s10439-021-02796-x
View details for PubMedID 34169398
Using parametric model order reduction for inverse analysis of large nonlinear cardiac simulations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING
Predictive high-fidelity finite element simulations of human cardiac mechanics commonly require a large number of structural degrees of freedom. Additionally, these models are often coupled with lumped-parameter models of hemodynamics. High computational demands, however, slow down model calibration and therefore limit the use of cardiac simulations in clinical practice. As cardiac models rely on several patient-specific parameters, just one solution corresponding to one specific parameter set does not at all meet clinical demands. Moreover, while solving the nonlinear problem, 90% of the computation time is spent solving linear systems of equations. We propose to reduce the structural dimension of a monolithically coupled structure-Windkessel system by projection onto a lower-dimensional subspace. We obtain a good approximation of the displacement field as well as of key scalar cardiac outputs even with very few reduced degrees of freedom, while achieving considerable speedups. For subspace generation, we use proper orthogonal decomposition of displacement snapshots. Following a brief comparison of subspace interpolation methods, we demonstrate how projection-based model order reduction can be easily integrated into a gradient-based optimization. We demonstrate the performance of our method in a real-world multivariate inverse analysis scenario. Using the presented projection-based model order reduction approach can significantly speed up model personalization and could be used for many-query tasks in a clinical setting.
View details for DOI 10.1002/cnm.3320
View details for Web of Science ID 000515254400001
View details for PubMedID 32022424
Automatic mapping of atrial fiber orientations for patient-specific modeling of cardiac electromechanics using image registration
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING
2019; 35 (6): e3190
Knowledge of appropriate local fiber architecture is necessary to simulate patient-specific electromechanics in the human heart. However, it is not yet possible to reliably measure in vivo fiber directions especially in human atria. Thus, we present a method that defines the fiber architecture in arbitrarily shaped atria using image registration and reorientation methods based on atlas atria with fibers predefined from detailed histological observations. Thereby, it is possible to generate detailed fiber families in every new patient-specific geometry in an automated, time-efficient process. We demonstrate the good performance of the image registration and fiber definition on 10 differently shaped human atria. Additionally, we show that characteristics of the electrophysiological activation pattern that appear in the atlas atria also appear in the patients' atria. We arrive to analogous conclusions for coupled electro-mechano-hemodynamical computations.
View details for DOI 10.1002/cnm.3190
View details for Web of Science ID 000471315800008
View details for PubMedID 30829001
View details for PubMedCentralID PMC6619047
The importance of the pericardium for cardiac biomechanics: from physiology to computational modeling
BIOMECHANICS AND MODELING IN MECHANOBIOLOGY
2019; 18 (2): 503–29
The human heart is enclosed in the pericardial cavity. The pericardium consists of a layered thin sac and is separated from the myocardium by a thin film of fluid. It provides a fixture in space and frictionless sliding of the myocardium. The influence of the pericardium is essential for predictive mechanical simulations of the heart. However, there is no consensus on physiologically correct and computationally tractable pericardial boundary conditions. Here, we propose to model the pericardial influence as a parallel spring and dashpot acting in normal direction to the epicardium. Using a four-chamber geometry, we compare a model with pericardial boundary conditions to a model with fixated apex. The influence of pericardial stiffness is demonstrated in a parametric study. Comparing simulation results to measurements from cine magnetic resonance imaging reveals that adding pericardial boundary conditions yields a better approximation with respect to atrioventricular plane displacement, atrial filling, and overall spatial approximation error. We demonstrate that this simple model of pericardial-myocardial interaction can correctly predict the pumping mechanisms of the heart as previously assessed in clinical studies. Utilizing a pericardial model not only can provide much more realistic cardiac mechanics simulations but also allows new insights into pericardial-myocardial interaction which cannot be assessed in clinical measurements yet.
View details for DOI 10.1007/s10237-018-1098-4
View details for Web of Science ID 000461323600018
View details for PubMedID 30535650
An adaptive hybridizable discontinuous Galerkin approach for cardiac electrophysiology
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING
2018; 34 (5): e2959
Cardiac electrophysiology simulations are numerically challenging because of the propagation of a steep electrochemical wave front and thus require discretizations with small mesh sizes to obtain accurate results. In this work, we present an approach based on the hybridizable discontinuous Galerkin method (HDG), which allows an efficient implementation of high-order discretizations into a computational framework. In particular, using the advantage of the discontinuous function space, we present an efficient p-adaptive strategy for accurately tracking the wave front. The HDG allows to reduce the overall degrees of freedom in the final linear system to those only on the element interfaces. Additionally, we propose a rule for a suitable integration accuracy for the ionic current term depending on the polynomial order and the cell model to handle high-order polynomials. Our results show that for the same number of degrees of freedom, coarse high-order elements provide more accurate results than fine low-order elements. Introducing p-adaptivity further reduces computational costs while maintaining accuracy by restricting the use of high-order elements to resolve the wave front. For a patient-specific simulation of a cardiac cycle, p-adaptivity reduces the average number of degrees of freedom by 95% compared to the nonadaptive model. In addition to reducing computational costs, using coarse meshes with our p-adaptive high-order HDG method also simplifies practical aspects of mesh generation and postprocessing.
View details for DOI 10.1002/cnm.2959
View details for Web of Science ID 000431995600005
View details for PubMedID 29316340
Multiphysics Modeling of the Atrial Systole under Standard Ablation Strategies
CARDIOVASCULAR ENGINEERING AND TECHNOLOGY
2017; 8 (2): 205–18
The aim of this study was to develop a computational framework to compare the impact of standard ablation concepts on the mechanical performance of the atria, since different line combinations cannot be applied in practice to the same patient. For this purpuse, we coupled electro-mechano-hemodynamic mathematical models based on biophysical principles and simulate the contractile performance of the atria. We computed systolic pressures and volumes in two patient-specific atrial geometries (one of normal size and one hypertrophied) with various ablation concepts. We found that our computational model is able to detect the differences in the left atrial contractility and ejection fraction for various electrical activation sequences resulting from different ablation line combinations. We show that multiphysics modeling has the potential to quantify the hemodynamic performance of left atria for different ablation lines, which could be used as additional pre-operative clinical information for the choice of the ablation concept in the future.
View details for DOI 10.1007/s13239-017-0308-z
View details for Web of Science ID 000402142300008
View details for PubMedID 28512679
On the Role of Mechanics in Chronic Lung Disease
2013; 6 (12): 5639-5658
Progressive airflow obstruction is a classical hallmark of chronic lung disease, affecting more than one fourth of the adult population. As the disease progresses, the inner layer of the airway wall grows, folds inwards, and narrows the lumen. The critical failure conditions for airway folding have been studied intensely for idealized circular cross-sections. However, the role of airway branching during this process is unknown. Here, we show that the geometry of the bronchial tree plays a crucial role in chronic airway obstruction and that critical failure conditions vary significantly along a branching airway segment. We perform systematic parametric studies for varying airway cross-sections using a computational model for mucosal thickening based on the theory of finite growth. Our simulations indicate that smaller airways are at a higher risk of narrowing than larger airways and that regions away from a branch narrow more drastically than regions close to a branch. These results agree with clinical observations and could help explain the underlying mechanisms of progressive airway obstruction. Understanding growth-induced instabilities in constrained geometries has immediate biomedical applications beyond asthma and chronic bronchitis in the diagnostics and treatment of chronic gastritis, obstructive sleep apnea and breast cancer.
View details for DOI 10.3390/ma6125639
View details for Web of Science ID 000330297600014
View details for PubMedCentralID PMC5452755