Parametric model order reduction and its application to inverse analysis of large nonlinear coupled cardiac problems

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 a novel approach to reduce only the structural dimension of the 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. To incorporate changes in the parameter set into our reduced order model, we provide a comparison of subspace interpolation methods. We further show how projection-based model order reduction can be easily integrated into a gradient-based optimization and demonstrate its performance 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.

[1]  Felipe Alonso Atienza,et al.  Inverse Problem of Electrocardiography: Estimating the Location of Cardiac Ischemia in a 3D Realistic Geometry , 2015, FIMH.

[2]  A. Quarteroni,et al.  A reduced computational and geometrical framework for inverse problems in hemodynamics , 2013, International journal for numerical methods in biomedical engineering.

[3]  C. Farhat,et al.  Interpolation Method for Adapting Reduced-Order Models and Application to Aeroelasticity , 2008 .

[4]  Theo Arts,et al.  Three-Wall Segment (TriSeg) Model Describing Mechanics and Hemodynamics of Ventricular Interaction , 2009, Annals of Biomedical Engineering.

[5]  Gianluigi Rozza,et al.  Reduced basis method for linear elasticity problems with many parameters , 2008 .

[6]  P. Moireau,et al.  Sequential parameter estimation for fluid–structure problems: Application to hemodynamics , 2012, International journal for numerical methods in biomedical engineering.

[7]  Wil H. A. Schilders,et al.  A Novel Approach to Model Order Reduction for Coupled Multiphysics Problems , 2014 .

[8]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947 .

[9]  Roy C. P. Kerckhoffs,et al.  Coupling of a 3D Finite Element Model of Cardiac Ventricular Mechanics to Lumped Systems Models of the Systemic and Pulmonic Circulation , 2006, Annals of Biomedical Engineering.

[10]  K. Kunisch,et al.  Control of the Burgers Equation by a Reduced-Order Approach Using Proper Orthogonal Decomposition , 1999 .

[11]  Alfio Quarteroni,et al.  Numerical approximation of parametrized problems in cardiac electrophysiology by a local reduced basis method , 2018, Computer Methods in Applied Mechanics and Engineering.

[12]  Stefan Volkwein,et al.  Galerkin Proper Orthogonal Decomposition Methods for a General Equation in Fluid Dynamics , 2002, SIAM J. Numer. Anal..

[13]  C. A. Figueroa,et al.  Sequential identification of boundary support parameters in a fluid-structure vascular model using patient image data , 2012, Biomechanics and Modeling in Mechanobiology.

[14]  Dominique Chapelle,et al.  A Galerkin strategy with Proper Orthogonal Decomposition for parameter-dependent problems – Analysis, assessments and applications to parameter estimation , 2013 .

[15]  Danny C. Sorensen,et al.  Nonlinear Model Reduction via Discrete Empirical Interpolation , 2010, SIAM J. Sci. Comput..

[16]  Dragan Stojic,et al.  Advances and Difficulties in Serbia’s Reindustrialization , 2016 .

[17]  Randall J. Allemang,et al.  THE MODAL ASSURANCE CRITERION–TWENTY YEARS OF USE AND ABUSE , 2003 .

[18]  W. Wall,et al.  Towards efficient uncertainty quantification in complex and large-scale biomechanical problems based on a Bayesian multi-fidelity scheme , 2014, Biomechanics and Modeling in Mechanobiology.

[19]  Karen Willcox,et al.  A Survey of Projection-Based Model Reduction Methods for Parametric Dynamical Systems , 2015, SIAM Rev..

[20]  Randall J. Allemang,et al.  A Correlation Coefficient for Modal Vector Analysis , 1982 .

[21]  G. Kerschen,et al.  The Method of Proper Orthogonal Decomposition for Dynamical Characterization and Order Reduction of Mechanical Systems: An Overview , 2005 .

[22]  Maxime Sermesant,et al.  Cardiac Function Estimation from MRI Using a Heart Model and Data Assimilation: Advances and Difficulties , 2005, FIMH.

[23]  Liang Zhong,et al.  Efficient estimation of personalized biventricular mechanical function employing gradient‐based optimization , 2018, International journal for numerical methods in biomedical engineering.

[24]  Nathan M. Newmark,et al.  A Method of Computation for Structural Dynamics , 1959 .

[25]  Brian D. Hong,et al.  Simulation of Left Ventricular Dynamics Using a Low-Order Mathematical Model , 2017, Cardiovascular Engineering and Technology.

[26]  Jack Lee,et al.  Multiphysics and multiscale modelling, data–model fusion and integration of organ physiology in the clinic: ventricular cardiac mechanics , 2016, Interface Focus.

[27]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[28]  M Caruel,et al.  Dimensional reductions of a cardiac model for effective validation and calibration , 2014, Biomechanics and modeling in mechanobiology.

[29]  Michael W Gee,et al.  A monolithic 3D‐0D coupled closed‐loop model of the heart and the vascular system: Experiment‐based parameter estimation for patient‐specific cardiac mechanics , 2017, International journal for numerical methods in biomedical engineering.

[30]  Jean-Frédéric Gerbeau,et al.  Approximated Lax pairs for the reduced order integration of nonlinear evolution equations , 2014, J. Comput. Phys..

[31]  A. Chatterjee An introduction to the proper orthogonal decomposition , 2000 .

[32]  Myrianthi Hadjicharalambous,et al.  Estimation of passive and active properties in the human heart using 3D tagged MRI , 2015, Biomechanics and Modeling in Mechanobiology.

[33]  Radomír Chabiniok,et al.  The importance of the pericardium for cardiac biomechanics: from physiology to computational modeling , 2018, Biomechanics and Modeling in Mechanobiology.

[34]  Jintai Chung,et al.  A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-α Method , 1993 .

[35]  Berend E. Westerhof,et al.  The arterial Windkessel , 2009, Medical & Biological Engineering & Computing.

[36]  C. Farhat,et al.  Structure‐preserving, stability, and accuracy properties of the energy‐conserving sampling and weighting method for the hyper reduction of nonlinear finite element dynamic models , 2015 .

[37]  Theo Arts,et al.  A model of the mechanics of the left ventricle. , 1979, Annals of biomedical engineering.

[38]  J. Marsden,et al.  Dimensional model reduction in non‐linear finite element dynamics of solids and structures , 2001 .

[39]  H. P. Lee,et al.  PROPER ORTHOGONAL DECOMPOSITION AND ITS APPLICATIONS—PART I: THEORY , 2002 .

[40]  Tamara G. Kolda,et al.  An overview of the Trilinos project , 2005, TOMS.

[41]  P Moireau,et al.  Estimation of tissue contractility from cardiac cine-MRI using a biomechanical heart model , 2012, Biomechanics and modeling in mechanobiology.

[42]  Jean-Frédéric Gerbeau,et al.  Reduced order model in cardiac electrophysiology with approximated Lax pairs , 2015, Adv. Comput. Math..

[43]  Hervé Delingette,et al.  Multifidelity-CMA: a multifidelity approach for efficient personalisation of 3D cardiac electromechanical models , 2017, Biomechanics and Modeling in Mechanobiology.

[44]  Kenneth Levenberg A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[45]  Alfio Quarteroni,et al.  A matrix DEIM technique for model reduction of nonlinear parametrized problems in cardiac mechanics , 2017 .

[46]  D. Marquardt An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .