Nonlinear object-oriented modeling based optimal control of the heart: Performing precise preload manipulation maneuvers using a ventricular assist device

This paper presents a simulation study of an optimal control approach being applied to control the volume trajectory of a heart's ventricle using a ventricular assist device (VAD). Usually VADs are used for heart failure therapy. However, in this study a novel VAD application within the context of animal trials will be investigated. Preload (the volume of blood within the ventricle just before the contraction) is an important determinant of heart function and therefore an interesting quantity to be manipulated in trials. To overcome limitations of classical preload manipulation techniques we propose to use a VAD to manipulate preload. The aim of this study will be to find the input trajectories for the VAD which are required to precisely hit specific preload values. Two different test scenarios were defined: the first scenario aims at resembling a traditional preload reduction technique. The second scenario defines a preload progression which is not possible to realize with traditional preload manipulation techniques. Besides the test scenarios, a comparison to experimentally obtained data is drawn. Object-oriented Modélica models of a simplified cardiovascular system and the VAD were developed to formulate a nonlinear optimal control problem. A fully open-source software tool-chain was used to automatically discretize the problem into a collocation scheme and to generate the necessary inputs for the nonlinear optimization solver. This tool-chain includes the software packages JModelica, CasADi and Ipopt. Object-oriented modeling in combination with the automated tool-chain allow for easy and fast modifications of the underlying models as well as the introduction of new objectives. Hence, the presented approach can easily be modified to be applied to other scopes of VAD and cardiovascular research.

[1]  Karl-Erik Årzén,et al.  Modeling and optimization with Optimica and JModelica.org - Languages and tools for solving large-scale dynamic optimization problems , 2010, Comput. Chem. Eng..

[2]  Gregor Ochsner,et al.  Numerical Optimal Control of Turbo Dynamic Ventricular Assist Devices , 2013 .

[3]  Patrick Amestoy,et al.  MUMPS : A General Purpose Distributed Memory Sparse Solver , 2000, PARA.

[4]  J. Clark,et al.  A dynamic model of ventricular interaction and pericardial influence. , 1997, The American journal of physiology.

[5]  Pascal Verdonck,et al.  Modeling Ventricular Function during Cardiac Assist: Does Time-Varying Elastance Work? , 2006, ASAIO journal.

[6]  Lorenz T. Biegler,et al.  On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming , 2006, Math. Program..

[7]  Gregor Ochsner,et al.  A physiological controller for turbodynamic ventricular assist devices based on a measurement of the left ventricular volume. , 2014, Artificial organs.

[8]  F. Colacino,et al.  Left Ventricle Load Impedance Control by Apical VAD Can Help Heart Recovery and Patient Perfusion: A Numerical Study , 2007, ASAIO journal.

[9]  Johan Åkesson Optimica—An Extension of Modelica Supporting Dynamic Optimization , 2008 .

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

[11]  R Rossaint,et al.  Xenon and isoflurane improved biventricular function during right ventricular ischemia and reperfusion , 2010, Acta anaesthesiologica Scandinavica.

[12]  Dirk Abel,et al.  Object-oriented Model Library of the Cardiovascular System Including Physiological Control Loops , 2009 .

[13]  Moritz Diehl,et al.  CasADi -- A symbolic package for automatic differentiation and optimal control , 2012 .

[14]  Dirk Abel,et al.  Cardiac Modeling: Identification of Subject Specific Left-Ventricular Isovolumic Pressure Curves , 2015 .

[15]  Francesco Casella,et al.  Integration of CasADi and JModelica.org , 2011 .