Description
Summary:In this thesis, we study mathematical models that describe the morphology of a generalized biological cell in equilibrium or under the influence of external forces. Within these models, the cell is considered as a thermodynamic system, where streaming effects in the cell bulk and the surrounding are coupled with a Helfrich-type model for the cell membrane. The governing evolution equations for the cell given in a continuum formulation are derived using an energy variation approach. Such two-phase flow problems that combine streaming effects with a free boundary problem that accounts for bending and surface tension can be described effectively by a diffuse interface approach. An advantage of the diffuse interface approach is that models for e.g. different biophysical processes can easily be combined. That makes this method suitable to describe complex phenomena such as cell motility and multi-cell dynamics. Within the first model for cell motility, we combine a biological network for GTPases with the hydrodynamic Helfrich-type model. This model allows to account for cell motility driven by membrane protrusion as a result of actin polymerization. Within the second model, we moreover extend the Helfrich-type model by an active gel theory to account for the actin filaments in the cell bulk. Caused by contractile stress within the actin-myosin solution, a spontaneous symmetry breaking event occurs that lead to cell motility. In this thesis, we further study the dynamics of multiple cells which is of wide interest since it reveals rich non-linear behavior. To apply the diffuse interface framework, we introduce several phase field variables to account for several cells that are coupled by a local interaction potential. In a first application, we study white blood cell margination, a biological phenomenon that results from the complex relation between collisions, different mechanical properties and lift forces of red blood cells and white blood cells within the vascular system. Here, it is shown that inertial effects, which can become of relevance in various parts of the cardiovascular system, lead to a decreasing tendency for margination with increasing Reynolds number. Finally, we combine the active polar gel theory and the multi-cell approach that is capable of studying collective migration of cells. This hydrodynamic approach predicts that collective migration emerges spontaneously forming coherently-moving clusters as a result of the mutual alignment of the velocity vectors during inelastic collisions. We further observe that hydrodynamics heavily influence those systems. However, a complete suppression of the onset of collective migration cannot be confirmed. Moreover, we give a brief insight how such highly coupled systems can be treated numerically using finite elements and how the numerical costs can be limited using operator splitting approaches and problem parallelization with OPENMP. === Diese Dissertation beschäftigt sich mit mathematischen Modellen zur Beschreibung von Gleichgewichts- und dynamischen Zuständen von verallgemeinerten biologischen Zellen. Die Zellen werden dabei als thermodynamisches System aufgefasst, bei dem Strömungseffekte innerhalb und außerhalb der Zelle zusammen mit einem Helfrich-Modell für Zellmembranen kombiniert werden. Schließlich werden durch einen Energie-Variations-Ansatz die Evolutionsgleichungen für die Zelle hergeleitet. Es ergeben sie dabei Mehrphasen-Systeme, die Strömungseffekte mit einem freien Randwertproblem, das zusätzlich physikalischen Einflüssen wie Biegung und Oberflächenspannung unterliegt, vereinen. Um solche Probleme effizient zu lösen, wird in dieser Arbeit die Diffuse-Interface-Methode verwendet. Ein Vorteil dieser Methode ist, dass es sehr einfach möglich ist, Modelle, die verschiedenste Prozesse beschreiben, miteinander zu vereinen. Dies erlaubt es, komplexe biologische Phänomene, wie zum Beispiel Zellmotilität oder auch die kollektive Bewegung von Zellen, zu beschreiben. In den Modellen für Zellmotilität wird ein biologisches Netzwerk-Modell für GTPasen oder auch ein Active-Polar-Gel-Modell, das die Aktinfilamente im Inneren der Zellen als Flüssigkristall auffasst, mit dem Multi-Phasen-Modell kombiniert. Beide Modelle erlauben es, komplexe Vorgänge bei der selbst hervorgerufenen Bewegung von Zellen, wie das Vorantreiben der Zellmembran durch Aktinpolymerisierung oder auch die Kontraktionsbewegung des Zellkörpers durch kontraktile Spannungen innerhalb des Zytoskelets der Zelle, zu verstehen. Weiterhin ist die kollektive Bewegung von vielen Zellen von großem Interesse, da sich hier viele nichtlineare Phänomene zeigen. Um das Diffuse-Interface-Modell für eine Zelle auf die Beschreibung mehrerer Zellen zu übertragen, werden mehrere Phasenfelder eingeführt, die die Zellen jeweils kennzeichnen. Schließlich werden die Zellen durch ein lokales Abstoßungspotential gekoppelt. Das Modell wird angewendet, um White blood cell margination, das die Annäherung von Leukozyten an die Blutgefäßwand bezeichnet, zu verstehen. Dieser Prozess wird dabei bestimmt durch den komplexen Zusammenhang zwischen Kollisionen, den jeweiligen mechanischen Eigenschaften der Zellen, sowie deren Auftriebskraft innerhalb der Adern. Die Simulationen zeigen, dass diese Annäherung sich in bestimmten Gebieten des kardiovaskulären Systems stark vermindert, in denen die Blutströmung das Stokes-Regime verlässt. Schließlich wird das Active-Polar-Gel-Modell mit dem Modell für die kollektive Bewegung vom Zellen kombiniert. Dies macht es möglich, die kollektive Bewegung der Zellen und den Einfluss von Hydrodynamik auf diese Bewegung zu untersuchen. Es zeigt sich dabei, dass der Zustand der kollektiven gerichteten Bewegung sich spontan aus der Neuausrichtung der jeweiligen Zellen durch inelastische Kollisionen ergibt. Obwohl die Hydrodynamik einen großen Einfluss auf solche Systeme hat, deuten die Simulationen nicht daraufhin, dass Hydrodynamik die kollektive Bewegung vollständig unterdrückt. Weiterhin wird in dieser Arbeit gezeigt, wie die stark gekoppelten Systeme numerisch gelöst werden können mit Hilfe der Finiten-Elemente-Methode und wie die Effizienz der Methode gesteigert werden kann durch die Anwendung von Operator-Splitting-Techniken und Problemparallelisierung mittels OPENMP.