Modeling, analysis, and simulation of thrombosis and hemostasis [Elektronische Ressource] / vorgelegt von Frédéric Weller

Sujets

• Gesundheit

Informations

• Publié par: ruprecht-karls-universitat_heidelberg
• Publié le: 01 janvier 2008
• Nombre de lectures: 19
• Langue: English
• Poids de l'ouvrage: 4 Mo

Extrait

Inaugural-Dissertation
zur
Erlangung der Doktorwürde
der
Naturwissenschaftlich-Mathematischen Gesamtfakultät
der
Ruprecht-Karls-Universität
Heidelberg
vorgelegt von
Diplom-Mathematiker Frédéric Weller
aus Stuttgart
Tag der mündlichen Prüfung: 14. Juli 2008Modeling, Analysis, and Simulation of
Thrombosis and Hemostasis
Gutachter: Prof. Dr. Dr. h.c. mult. Willi Jäger
Prof. Dr. Rolf RannacherAbstract
This thesis investigates the influences of shear stress, saturation-dependent changes in surface reactivity,
and thrombus growth on platelet deposition to reactive materials, which is of paramount interest
in bioengineering and clinical practice. For this purpose, two mathematical models based on the
Navier-Stokes equations and on particle conservation are developed. The first model is formulated on a
fixed domain (“FD-model”) and describes the initial phase of platelet adhesion, whereas the second
one is a free boundary problem capturing long-term thrombus growth. Several vessel geometries are
considered: Stagnation point flow, tubular expansion, and t-junction. Model parameters are optimized
to fit the data and their so obtained values are justified on the basis of experimental observations.
The FD-model does not match the experimental data at all, when platelet adhesion is assumed
independent of shear stress. In contrast, when adhesion is assumed shear-dependent, at least qualitative
agreement is achieved. Solely by consideration of both shear stress and saturation-dependent changes
in surface reactivity, good quantitative agreement of FD-model and data is obtainable. Such changes in
surface reactivity are taken into account by coupling platelet flux conditions to ordinary differential
equations (ODEs) for the evolution of surface-bound platelets. The free boundary problem is simulated
by the level set method. Like the FD-model, it shows good qualitative agreement with the experimental
evidence when shear stress is taken into account, whereas negligence of shear leads to completely false
predictions.
Regarding mathematical well-posedness of the FD-model, existence of weak solutions is shown for
generalized parabolic systems having ODE-coupled flux conditions. Uniqueness and positivity of
solutions are also investigated. Regarding the free boundary problem, a detailed proof of classical
solvability in terms of Hölder spaces is presented.
Zusammenfassung
Diese Arbeit untersucht die Einflüsse von Scherkräften, sättigungsbedingten Änderungen der Oberflä-
chenreaktivität und die Auswirkungen des Thrombenwachstums auf die Adhäsion von Blutplättchen
an reaktiven Materialien, was von großem Interesse in Biotechnik und klinischer Praxis ist. Zu diesem
Zweck werden zwei mathematische Modelle basierend auf den Navier-Stokes Gleichungen und der
Teilchenerhaltung entwickelt. Das erste Modell beschreibt die Anfangsphase der Plättchenadhäsion und
nimmt daher ein festes Gebiet an („FD-Modell“), wohingegen das zweite ein freies Randwertproblem
zur Beschreibung des längerfristigen Thrombenwachstums darstellt. Es werden mehrere Gefäßgeometri-
en betrachtet: Staupunkt, Gefäßerweiterung und T-Kreuzung. Die Modellparameter werden anhand
der experimentellen Daten gefittet und ihre so erhaltenen Werte mit experimentellen Beobachtungen
gerechtfertigt.
Beim FD-Modell stellt sich bei Annahme scherkraftunabhängiger Plättchenadhäsion keine Übereinstim-
mung mit den experimentellen Daten ein. Bei scherkraftabhängiger Adhäsion wird immerhin qualitative
Übereinstimmung erzielt. Gute quantitative Übereinstimmung mit den Daten zeigt das FD-Modell
dagegen nur bei gleichzeitiger Berücksichtigung von Scherkräften und sättigungsbedingten Änderungen
der Oberflächenreaktivität. Letzteres wird über eine Kopplung des Plättchenflusses mit gewöhnlichen
Differenzialgleichungen für die zeitliche Entwicklung der Konzentration gebundener Plättchen realisiert.
Das freie Randwertproblem wird mit Hilfe der Level Set Methode numerisch simuliert. Ebenso wie das
FD-Modell zeigt es bei Berücksichtigung der Scherkräfte gute qualitative Übereinstimmung mit dem
Experiment, wohingegen die Vernachlässigung der Scherkräfte zu völlig falschen Voraussagen führt.
Zur Sicherstellung der mathematischen Wohlgestelltheit des FD-Modells wird die Existenz schwacher
Lösungen für allgemeinere parabolische Systeme, die solch gekoppelten Randbedingungen unterworfen
sind, gezeigt. Darüberhinaus wird die Eindeutigkeit und Positivität der Lösung untersucht. Für das
freie Randwertproblem wird dessen klassische Lösbarkeit in Hölderräumen umfassend bewiesen.Contents
1 Introduction 1
2 Biological background 5
2.1 Primary hemostasis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.2 Secondary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3 The fixed domain model 13
3.1 Model development . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3.2 Numerical methods and optimization of parameters . . . . . . . . . . . . . . . . 17
3.3 Numerical results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3.3.1 Stagnation point flow . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3.3.2 Tubular expansion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
3.3.3 T-junction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
3.4 Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
3.5 Mathematical analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3.5.1 Notation and function spaces . . . . . . . . . . . . . . . . . . . . . . . . 27
3.5.2 Problem statement and assumptions . . . . . . . . . . . . . . . . . . . . 28
3.5.3 Existence and uniqueness for the linear problem . . . . . . . . . . . . . 29
3.5.4 and for the nonlinear problem . . . . . . . . . . . 33
3.5.5 Positivity of solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
4 The free boundary problem 39
4.1 Model development . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.2 The level set method and its implementation . . . . . . . . . . . . . . . . . . . 41
4.3 Numerical results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
5 Classical solvability of the free boundary problem 53
5.1 Preliminaries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
5.1.1 Geometry and model equations . . . . . . . . . . . . . . . . . . . . . . . 53
5.1.2 Notation and function spaces . . . . . . . . . . . . . . . . . . . . . . . . 54
5.1.3 Assumptions and main theorem . . . . . . . . . . . . . . . . . . . . . . . 55
5.2 Transformation to fixed domain . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
5.3 The linear problem for the platelets; fixed flow . . . . . . . . . . . . . . . . . . 58
5.3.1 Investigation of auxiliary problems in half space . . . . . . . . . . . . . 59
5.3.2 Solution of the linear problem in Q . . . . . . . . . . . . . . . . . . . . 79T
5.4 Solution of the nonlinear problem for the platelets; fixed flow . . . . . . . . . . 86
5.5 of the full nonlinear problem including flow . . . . . . . . . . . . . . . 94
iContents
6 Conclusions and outlook 97
Acknowledgments 99
List of Figures 101
List of Tables 103
Bibliography 105
ii1 Introduction
Hemostasis is responsible to stem blood loss after injury by platelet plug formation. Although
being life essential, a major part of deaths in the western society is due to thrombotic events
provoked by disorders of the hemostatic system. Therefore, a better understanding of the
underlying mechanisms is needed. It is known that the overall process is governed by Virchow’s
triad [94] which comprises composition of blood, surface reactivity, and flow. Nevertheless,
previous mathematical models, such as [5, 24, 37, 49, 82, 83], concentrated mainly on kinetics
and, if at all, accounted for the in vivo flow situation only in simple terms of transport.
But, since the end of the 1970s there is experimental evidence that shear stress strongly
influences the activation of platelets and their adhesion to injured tissue [88, 95]—both crucial
steps in plug formation. However, there is still confusion about the exact physical quantities
that determine spatial platelet distribution and give rise to the well known fact that sites of
disturbed flow are prone to platelet deposition [32, 90].
This thesis investigates the influences of shear stress, saturation-dependent changes in surface
reactivity, and thrombus growth on the adhesion and aggregation of platelets to reactive
materials. For this purpose, two mathematical models based on fluid dynamic and species
conservation equations are developed and their ability to match given in vitro experimental
data is studied—dependent on whether the above mentioned effects are taken into account
or not. The first model describes the initial phase of platelet deposition, when thrombus
growth can be neglected. Therefore, this model is formulated on a fixed domain and henceforth
referred to as “Fixed Domain model” (FD-model). It accounts for shear stress and changes