Numerical modelling of fluid-structure interaction with application in hemodynamics

The thesis deals with numerical algorithms for fluid-structure interaction problems with application in blood flow modelling. It starts with a short introduction on the mathematical description of incompressible viscous flow with non-Newtonian viscosity and a moving linear viscoelastic structure. The mathematical model consists of the generalized Navier-Stokes equation used for the description of fluid flow and the generalized string model for structure movement. The arbitrary Lagrangian-Eulerian approach is used in order to take into account moving computational domain. A part of the thesis is devoted to the discussion on the non-Newtonian behaviour of shear-thinning fluids, which is in our case blood, and derivation of two non-Newtonian models frequently used in the blood flow modelling. Further we give a brief overview on recent fluid-structure interaction schemes with discussion about the difficulties arising in numerical modelling of blood flow. Our main contribution lies in numerical and experimental study of a new loosely-coupled partitioned scheme called the kinematic splitting fluid-structure interaction algorithm. We present stability analysis for a coupled problem of non-Newtonian shear-dependent fluids in moving domains with viscoelastic boundaries. Here, we assume both, the nonlinearity in convective as well is diffusive term. We analyse the convergence of proposed numerical scheme for a simplified fluid model of the Oseen type. Moreover, we present series of experiments including numerical error analysis, comparison of hemodynamic parameters for the Newtonian and non-Newtonian fluids and comparison of several physiologically relevant computational geometries in terms of wall displacement and wall shear stress. Numerical analysis and extensive experimental study for several standard geometries confirm reliability and accuracy of the proposed kinematic splitting scheme in order to approximate fluid-structure interaction problems.

Fourier-transform rheology on anionically synthesised polymer melts and solutions of various topology

Aim of this thesis was to further extend the applicability of the Fourier-transform (FT) rheology technique especially for non-linear mechanical characterisation of polymeric materials on the one hand and to investigated the influence of the degree of branching on the linear and non-linear relaxation behaviour of polymeric materials on the other hand. The latter was achieved by employing in particular FT-rheology and other rheological techniques to variously branched polymer melts and solutions. For these purposes, narrowly distributed linear and star-shaped polystyrene and polybutadiene homo-polymers with varying molecular weights were anionically synthesised using both high-vacuum and inert atmosphere techniques. Furthermore, differently entangled solutions of linear and star-shaped polystyrenes in di-sec-octyl phthalate (DOP) were prepared. The several linear polystyrene solutions were measured under large amplitude oscillatory shear (LAOS) conditions and the non-linear torque response was analysed in the Fourier space. Experimental results were compared with numerical predictions performed by Dr. B. Debbaut using a multi-mode differential viscoelastic fluid model obeying the Giesekus constitutive equation. Apart from the analysis of the relative intensities of the harmonics, a detailed examination of the phase information content was developed. Further on, FT-rheology allowed to distinguish polystyrene melts and solutions due to their different topologies where other rheological measurements failed. Significant differences occurred under LAOS conditions as particularly reflected in the phase difference of the third harmonic, Ħ3, which could be related to shear thinning and shear thickening behaviour.

Influences of deep convective cloud systems on tropospheric trace gases and photochemistry over the tropical West Pacific

A numerical model for studying the influences of deep convective cloud systems on photochemistry was developed based on a non-hydrostatic meteorological model and chemistry from a global chemistry transport model. The transport of trace gases, the scavenging of soluble trace gases, and the influences of lightning produced nitrogen oxides (NOx=NO+NO2) on the local ozone-related photochemistry were investigated in a multi-day case study for an oceanic region located in the tropical western Pacific. Model runs considering influences of large scale flows, previously neglected in multi-day cloud resolving and single column model studies of tracer transport, yielded that the influence of the mesoscale subsidence (between clouds) on trace gas transport was considerably overestimated in these studies. The simulated vertical transport and scavenging of highly soluble tracers were found to depend on the initial profiles, reconciling contrasting results from two previous studies. Influences of the modeled uptake of trace gases by hydrometeors in the liquid and the ice phase were studied in some detail for a small number of atmospheric trace gases and novel aspects concerning the role of the retention coefficient (i.e. the fraction of a dissolved trace gas that is retained in the ice phase upon freezing) on the vertical transport of highly soluble gases were illuminated. Including lightning NOx production inside a 500 km 2-D model domain was found to be important for the NOx budget and caused small to moderate changes in the domain averaged ozone concentrations. A number of sensitivity studies yielded that the fraction of lightning associated NOx which was lost through photochemical reactions in the vicinity of the lightning source was considerable, but strongly depended on assumptions about the magnitude and the altitude of the lightning NOx source. In contrast to a suggestion from an earlier study, it was argued that the near zero upper tropospheric ozone mixing ratios which were observed close to the study region were most probably not caused by the formation of NO associated with lightning. Instead, it was argued in agreement with suggestions from other studies that the deep convective transport of ozone-poor air masses from the relatively unpolluted marine boundary layer, which have most likely been advected horizontally over relatively large distances (both before and after encountering deep convection) probably played a role. In particular, it was suggested that the ozone profiles observed during CEPEX (Central Equatorial Pacific Experiment) were strongly influenced by the deep convection and the larger scale flow which are associated with the intra-seasonal oscillation.

Theoretical studies of fluid membrane mechanics

Biologische Membranen sind Fettmolekül-Doppelschichten, die sich wie zweidimensionale Flüssigkeiten verhalten. Die Energie einer solchen fluiden Oberfläche kann häufig mit Hilfe eines Hamiltonians beschrieben werden, der invariant unter Reparametrisierungen der Oberfläche ist und nur von ihrer Geometrie abhängt. Beiträge innerer Freiheitsgrade und der Umgebung können in den Formalismus mit einbezogen werden. Dieser Ansatz wird in der vorliegenden Arbeit dazu verwendet, die Mechanik fluider Membranen und ähnlicher Oberflächen zu untersuchen. Spannungen und Drehmomente in der Oberfläche lassen sich durch kovariante Tensoren ausdrücken. Diese können dann z. B. dazu verwendet werden, die Gleichgewichtsposition der Kontaktlinie zu bestimmen, an der sich zwei aneinander haftende Oberflächen voneinander trennen. Mit Ausnahme von Kapillarphänomenen ist die Oberflächenenergie nicht nur abhängig von Translationen der Kontaktlinie, sondern auch von Änderungen in der Steigung oder sogar Krümmung. Die sich ergebenden Randbedingungen entsprechen den Gleichgewichtsbedingungen an Kräfte und Drehmomente, falls sich die Kontaktlinie frei bewegen kann. Wenn eine der Oberflächen starr ist, muss die Variation lokal dieser Fläche folgen. Spannungen und Drehmomente tragen dann zu einer einzigen Gleichgewichtsbedingung bei; ihre Beiträge können nicht mehr einzeln identifiziert werden. Um quantitative Aussagen über das Verhalten einer fluiden Oberfläche zu machen, müssen ihre elastischen Eigenschaften bekannt sein. Der "Nanotrommel"-Versuchsaufbau ermöglicht es, Membraneigenschaften lokal zu untersuchen: Er besteht aus einer porenüberspannenden Membran, die während des Experiments durch die Spitze eines Rasterkraftmikroskops in die Pore gedrückt wird. Der lineare Verlauf der resultierenden Kraft-Abstands-Kurven kann mit Hilfe der in dieser Arbeit entwickelten Theorie reproduziert werden, wenn der Einfluss von Adhäsion zwischen Spitze und Membran vernachlässigt wird. Bezieht man diesen Effekt in die Rechnungen mit ein, ändert sich das Resultat erheblich: Kraft-Abstands-Kurven sind nicht länger linear, Hysterese und nichtverschwindende Trennkräfte treten auf. Die Voraussagen der Rechnungen könnten in zukünftigen Experimenten dazu verwendet werden, Parameter wie die Biegesteifigkeit der Membran mit einer Auflösung im Nanometerbereich zu bestimmen. Wenn die Materialeigenschaften bekannt sind, können Probleme der Membranmechanik genauer betrachtet werden. Oberflächenvermittelte Wechselwirkungen sind in diesem Zusammenhang ein interessantes Beispiel. Mit Hilfe des oben erwähnten Spannungstensors können analytische Ausdrücke für die krümmungsvermittelte Kraft zwischen zwei Teilchen, die z. B. Proteine repräsentieren, hergeleitet werden. Zusätzlich wird das Gleichgewicht der Kräfte und Drehmomente genutzt, um mehrere Bedingungen an die Geometrie der Membran abzuleiten. Für den Fall zweier unendlich langer Zylinder auf der Membran werden diese Bedingungen zusammen mit Profilberechnungen kombiniert, um quantitative Aussagen über die Wechselwirkung zu treffen. Theorie und Experiment stoßen an ihre Grenzen, wenn es darum geht, die Relevanz von krümmungsvermittelten Wechselwirkungen in der biologischen Zelle korrekt zu beurteilen. In einem solchen Fall bieten Computersimulationen einen alternativen Ansatz: Die hier präsentierten Simulationen sagen voraus, dass Proteine zusammenfinden und Membranbläschen (Vesikel) bilden können, sobald jedes der Proteine eine Mindestkrümmung in der Membran induziert. Der Radius der Vesikel hängt dabei stark von der lokal aufgeprägten Krümmung ab. Das Resultat der Simulationen wird in dieser Arbeit durch ein approximatives theoretisches Modell qualitativ bestätigt.

Magnetic tunneling junctions with the Heusler compound Co 2 Cr 0.6 Fe 0.4 Al

Materialen mit sehr hoher Spinpolarisation werden für Anwendungen im Bereich der Spin-Elektronik benötigt. Deshalb werden große Forschungsanstrengungen zur Untersuchung der Eigenschaften von Verbindungen mit potentiell halbmetallischem Charakter, d. h.mit 100% Spinpolarisation, unternommen. In halbmetallischen Verbindungen, erwartet man eine Lücke in der Zustandsdichte an der Fermi Energie für Ladungsträger einer Spinrichtung, wahrend die Ladungsträger mit der anderen Spinrichtung sich metallisch verhalten. Eine Konsequenz davon ist, dass ein Strom, der durch solche Verbindung fließt, voll spinpolarisiert ist. Die hohe Curie-Temperatur Tc (800 K) und der theoretisch vorhergesagte halbmetallische Charakter machen Co2Cr0.6Fe0.4Al (CCFA) zu einem guten Kandidaten für Spintronik-Anwendungen wie magnetische Tunnelkontakte (MTJs = Magnetic Tunneling Junctions). In dieser Arbeit werden die Ergebnisse der Untersuchung der elektronischen und strukturellen Eigenschaften von dünnen CCFA Schichten dargestellt. Diese Schichten wurden in MTJs integriert und der Tunnel-Magnetowiderstands-Effekt untersucht. Hauptziele waren die Messung der Spinpolarisation und Untersuchungen der elektronischen Struktur von CCFA. Der Einfluss verschiedener Depositionsparameter auf die Eigenschaften der Schichten, speziell auf der Oberflächenordnung und damit letztlich auf den Tunnel-Magnetowiderstand (TMR), wurde bestimmt. Epitaktische d¨unne CCFA Schichten mit zwei verschiedenen Wachstumsrichtungen wurden auf verschiedene Substrate und Pufferschichten deponiert. Ein Temperverfahren wurde eingesetzt um die strukturelle Eigenschaften der dünnen Schichten zu verbessern. Für die MTJs wurde Al2O3 als Barrierenmaterial verwendet und Co als Gegenelektrode gewählt. Die Mehrschicht-Systeme wurden in Mesa-Geometrie mit lithographischen Methoden strukturiert. Eine maximal Jullière Spinpolarisation von 54% wurde an Tunnelkontakte mit epitaktischen CCFA Schichten gemessen. Ein starker Einfluss der Tempernbedingungen auf dem TMR wurde festgestellt. Eine Erhörung des TMR wurde mit einer Verbesserung der Oberflächenordung der CCFA Schichten korreliert. Spektroskopische Messungen wurden an den MTJs durchgeführt. Diesen Messungen liefern Hinweise auf inelastische Elektron-Magnon und Elektron-Phonon Stossprozesse an den Grenzflächen. Einige der beobachteten Strukturen konnten mit der berechneten elektronischen Struktur von CCFA korreliert worden.

Palaeo wind system reconstruction of the last glacial period over Europe, using high resolution proxy data and model-data-comparison

The atmosphere is a global influence on the movement of heat and humidity between the continents, and thus significantly affects climate variability. Information about atmospheric circulation are of major importance for the understanding of different climatic conditions. Dust deposits from maar lakes and dry maars from the Eifel Volcanic Field (Germany) are therefore used as proxy data for the reconstruction of past aeolian dynamics.rnrnIn this thesis past two sediment cores from the Eifel region are examined: the core SM3 from Lake Schalkenmehren and the core DE3 from the Dehner dry maar. Both cores contain the tephra of the Laacher See eruption, which is dated to 12,900 before present. Taken together the cores cover the last 60,000 years: SM3 the Holocene and DE3 the marine isotope stages MIS-3 and MIS-2, respectively. The frequencies of glacial dust storm events and their paleo wind direction are detected by high resolution grain size and provenance analysis of the lake sediments. Therefore two different methods are applied: geochemical measurements of the sediment using µXRF-scanning and the particle analysis method RADIUS (rapid particle analysis of digital images by ultra-high-resolution scanning of thin sections).rnIt is shown that single dust layers in the lake sediment are characterized by an increased content of aeolian transported carbonate particles. The limestone-bearing Eifel-North-South zone is the most likely source for the carbonate rich aeolian dust in the lake sediments of the Dehner dry maar. The dry maar is located on the western side of the Eifel-North-South zone. Thus, carbonate rich aeolian sediment is most likely to be transported towards the Dehner dry maar within easterly winds. A methodology is developed which limits the detection to the aeolian transported carbonate particles in the sediment, the RADIUS-carbonate module.rnrnIn summary, during the marine isotope stage MIS-3 the storm frequency and the east wind frequency are both increased in comparison to MIS-2. These results leads to the suggestion that atmospheric circulation was affected by more turbulent conditions during MIS-3 in comparison to the more stable atmospheric circulation during the full glacial conditions of MIS-2.rnThe results of the investigations of the dust records are finally evaluated in relation a study of atmospheric general circulation models for a comprehensive interpretation. Here, AGCM experiments (ECHAM3 and ECHAM4) with different prescribed SST patterns are used to develop a synoptic interpretation of long-persisting east wind conditions and of east wind storm events, which are suggested to lead to an enhanced accumulation of sediment being transported by easterly winds to the proxy site of the Dehner dry maar.rnrnThe basic observations made on the proxy record are also illustrated in the 10 m-wind vectors in the different model experiments under glacial conditions with different prescribed sea surface temperature patterns. Furthermore, the analysis of long-persisting east wind conditions in the AGCM data shows a stronger seasonality under glacial conditions: all the different experiments are characterized by an increase of the relative importance of the LEWIC during spring and summer. The different glacial experiments consistently show a shift from a long-lasting high over the Baltic Sea towards the NW, directly above the Scandinavian Ice Sheet, together with contemporary enhanced westerly circulation over the North Atlantic.rnrnThis thesis is a comprehensive analysis of atmospheric circulation patterns during the last glacial period. It has been possible to reconstruct important elements of the glacial paleo climate in Central Europe. While the proxy data from sediment cores lead to a binary signal of the wind direction changes (east versus west wind), a synoptic interpretation using atmospheric circulation models is successful. This shows a possible distribution of high and low pressure areas and thus the direction and strength of wind fields which have the capacity to transport dust. In conclusion, the combination of numerical models, to enhance understanding of processes in the climate system, with proxy data from the environmental record is the key to a comprehensive approach to paleo climatic reconstruction.rn

The characterization of the Atxn2-CAG42-knock-in mouse as a model for Spinocerebellar Ataxia Type 2

Die Ursache der neurodegenerativen Erkrankung Spinozerebelläre Ataxie Typ 2 (SCA2) ist eine expandierte Polyglutamin-Domäne im humanen ATXN2-Gen von normalerweise 22 auf über 31 CAGs. Von der Degeneration sind vorwiegend die zerebellären Purkinje Neuronen betroffen, in denen zunehmend zytoplasmatische Aggregate sichtbar werden. Auch wenn die genaue Funktion von ATXN2 und die zugrunde liegenden molekularen Mechanismen noch immer ungeklärt sind, werden ein toxischer Funktionsgewinn sowie der Verlust der normalen Proteinfunktion als mögliche Ursachen diskutiert.rnUm ein wirklichkeitsgetreues Tiermodell für die SCA2 zu haben, wurde eine knock-in Maus generiert, deren einzelnes CAG im Atxn2-Gen durch 42 CAGs ersetzt wurde. Dieses Mausmodell ist durch eine stabile Vererbung der Expansion charakterisiert. Weiterhin zeigt sie ein verringertes Körpergewicht sowie eine spät beginnende motorische Inkoordination, was dem Krankheitsbild von SCA2 entspricht. rnIm Weiteren konnte gezeigt werden, dass, obwohl die Atxn2 mRNA-Spiegel in Großhirn und Kleinhirn erhöht waren, die Menge an löslichem ATXN2 im Laufe der Zeit abnahm und dies mit einem Auftreten an unlöslichem ATXN2 korrelierte. Dieser im Kleinhirn progressive Prozess resultierte schließlich in zytoplasmatischen Aggregaten innerhalb der Purkinje Neuronen alter Mäuse. Der Verlust an löslichem ATXN2 könnte Effekte erklären, die auf einen partiellen Funktionsverlust von ATXN2 zurückzuführen sind, wobei die Aggregatbildung einen toxischen Funktionsgewinn wiederspiegeln könnte. Neben ATXN2 wurde auch sein Interaktor PABPC1 zunehmend unlöslich. Während dies im Großhirn eine Erhöhung der PABPC1 mRNA- und löslichen Proteinspiegel zur Folge hatte, konnte keine kompensatorische Veränderung seiner mRNA und zudem eine Verminderung an löslichem PABPC1 im Kleinhirn beobachtet werden. Auch PABPC1 wurde in Aggregate sequestriert. Diese Unterschiede zwischen Großhirn und Kleinhirn könnten zu der spezifischen Vulnerabilität des Kleinhirns beitragen.rnUm die Folgen auf mRNA-Prozessierung zu untersuchen, wurde ein Transkriptomprofil im mittleren sowie fortgeschrittenen Alter der Mäuse erstellt. Hierbei war eine erhöhte Expression von Fbxw8 im Kleinhirn alter Mäuse auffällig. Als Komponente eines Ubiquitin-E3-Ligase-Komplexes, hilft FBXW8 in der Degradierung von Zielproteinen und könnte somit die Toxizität des expandieren ATXN2 verringern. rnZur näheren Beschreibung der physiologischen Funktion von ATXN2, konnte in ATXN2-knock-out Mäusen gezeigt werden, dass das Fehlen von ATXN2 zu einer reduzierten globalen Proteinsyntheserate führte und somit eine Rolle als Translationsaktivator möglich erscheint. Kompensatorisch wurde eine erhöhte S6-Phosphorylierung gemessen.