Mathematical Modeling and Sensitivity Analysis of Lumped-Parameter Model of the Human Cardiovascular System

Abstract

In this thesis, a lumped-parameter model of the systemic circulation is developed and sensitivity analysis (local and global) is applied to quantify the impact of input parameters on output variables. The thesis is divided into two major parts, (1) cardiovascular (CV) modeling and (2) sensitivity analysis. In the CV modeling part, the major arteries of the systemic circulation are modeled explaining the elastic and the visco-elastic vessel walls behavior. In the sensitivity analysis part of the thesis, two local sensitivity analysis (LSA) and three variance-based global sensitivity analysis (GSA) methods are applied on a full and partial cardiovascular system (CVS). LSA is applied on a linear elastic model of arm arteries (with and without anastomosis). While GSA is applied on MACSim (Major Arterial Cardiovascular Simulator), visco-elastic model of the CVS and carotid bifurcation. The main aim of this research work is to provide a general framework of parameter estimation using sensitivity analysis. However, the objectives of this thesis are summarized as: • Formulate a reliable and clinically relevant model of the systemic circulation using the lumped- parameter approach (linear elastic and visco-elastic). • Identify and rank the most important CV parameters, which contribute most on the output uncertainty. Furthermore, these sensitive parameter can be accurately estimated from measurements. • Find the optimal measurement locations and optimal time regions in pressure and flow waves, which are very helpful in the parameter estimation. • Find a GSA method for multi-compartment lumped-parameter model of the CVS on the basis of computational cost, simplicity and straightforward implementation. • Discuss and explain the combine effect of vasodilation and vasoconstriction in linear elastic model of MACSim (Major Arterial Cardiovascular Simulator) by considering the sensitivity of the boundary resistance. • Study different levels of stenosis and aneurysm in visco-elastic model of the CVS. The study could benefit medical students and doctors to understand the dependence of hemodynamic variables on CV model parameters.

Die grundlegenden physikalischen Eigenschaften des Herz-Kreislaufsystems lassen sich durch dimensionsreduzierte parametrisierte Multi-Kompartiment- Modelle mit Hilfe von computergest¨utzten Rechenverfahren effizient beschreiben. Einblicke in das Verhalten des Herz-Kreislaufsystems bietet dabei vor allem die Ver¨anderung von Modellparametern. In dieser Dissertation werden die systemischen Hauptarterien des menschlichen Kreislaufsystems in 122 Kompartimente (Knoten oder Segmente) unterteilt. Die Zustandsvariablen (Druck und Fluss) eines jeden Segments werden durch 3-elektrische (RCL) bzw. 4-strukturelle Parameter (Eldh) beschrieben. In computergest¨utzten Herz- Kreislaufmodellen stellen die Parameter die Hauptquelle der Unsicherheit dar, weil diese sich direkt auf die Unsicherheit der Zustandsvariablen und damit auf die Verl¨asslichkeit des Modells auswirken. Bei der Erstellung von patienten-spezifischen Herz-Kreislaufmodellen die sich zur Untersuchung von Herz-Kreislauferkrankungen eignen, werden die Parameter aus patienten- spezifischen Messdaten gesch¨atzt. In Multi-Kompartiment Herz- Kreislaufmodellen gibt es eine Vielzahl von elektrischen und strukturellen Parametern, eine genaue Sch¨atzung aller Parameter aus gegebenen Messdaten ist deshalb nur bedingt m¨oglich. Ein m¨oglicher Ausweg aus dieser Situation ergibt sich durch die Quantifizierung der einflussreichsten Modellparameter bez¨uglich der Zustandsvariablen erg¨anzt durch optimaleMessstellen und Zeitr¨aume, um aus den patienten-spezifischen Daten die einflussreichsten Parameter zu sch¨atzen. In dieser Arbeit wurden die Modellparameter mit Hilfe der Sensitivit¨atsanalyse bez¨uglich ihrer Wichtigkeit geordnet, sowie optimale Messstellen und Zeitr¨aume in den Druck- und Flusswellenformen identifiziert. Hierf¨ur wurde die lokale Sensitivit ¨atsanalyse auf die Armarterien und die globale Sensitivit¨atsanalyse auf eine Karotis Bifurkation und den gesamten Herz-Kreislauf angewendet. Die Ergebnisse zeigen eine starke Ortsabh¨angigkeit der Sensitivit¨at bez¨uglich Druck und Fluss, d.h. in jedem Ort des arteriellen Netzwerks ergab sich eine unterschiedliche Ordnung der einflussreichsten Parameter. Grunds¨atzlich erwiesen sich der Durchmesser, die L¨ange des Gef¨aßes und der Widerstand des Blutflusses, speziell an den Gef¨aßenden als einflussreiche Parameter des Herz-Kreislaufsystems. Weiterhin bietet die vorliegende Arbeit eine Orientierung f¨ur Experimentatoren, welche Messgr¨oßen (Druck und Fluss), an welchen Orten (optimale Messstellen) und zu welchen Zeiten bei der Sch¨atzung von einflussreichen Parametern ber¨ucksichtigt werden sollten. Die optimalen Druck-/Flussmessungen, die zur Sch¨atzung von einflussreichen Parametern verwendet werden sollten, h¨angen stark vom Ort im Netzwerk ab, wohingegen sich die optimalen Zeitr¨aume im Druck- und Flussverlauf bez¨uglich der einzelnen Parameter gleichen, d.h. in der fr¨uhen Systole, End Systole und fr¨uhen Diastole.