Application of reactive fluid transport modeling to hydrothermal systems

A one-dimensional multi-component reactive fluid transport algorithm, 1DREACT (Steefel, 1993) was used to investigate different fluid-rock interaction systems. A major short coming of mass transport calculations which include mineral reactions is that solid solutions occurring in many minerals are not treated adequately. Since many thermodynamic models of solid solutions are highly non-linear, this can seriously impact on the stability and efficiency of the solution algorithms used. Phase petrology community saw itself faced with a similar predicament 10 years ago. To improve performance and reliability, phase equilibrium calculations have been using pseudo compounds. The same approach is used here in the first, using the complex plagioclase solid solution as an example. Thermodynamic properties of a varying number of intermediate plagioclase phases were calculated using ideal molecular, Al-avoidance, and non-ideal mixing models. These different mixing models can easily be incorporated into the simulations without modification of the transport code. Simulation results show that as few as nine intermediate compositions are sufficient to characterize the diffusional profile between albite and anorthite. Hence this approach is very efficient, and can be used with little effort. A subsequent chapter reports the results of reactive fluid transport modeling designed to constrain the hydrothermal alteration of Paleoproterozoic sediments of the Southern Lake Superior region. Field observations reveal that quartz-pyrophyllite (or kaolinite) bearing assemblages have been transformed into muscovite-pyrophyllite-diaspore bearing assemblages due to action of fluids migrating along permeable flow channels. Fluid-rock interaction modeling with an initial qtz-prl assemblage and a K-rich fluid simulates the formation of observed mineralogical transformation. The bulk composition of the system evolves from an SiO2-rich one to an Al2O3+K2O-rich one. Simulations show that the fluid flow was up-temperature (e.g. recharge) and that fluid was K-rich. Pseudo compound approach to include solid solutions in reactive transport models was tested in modeling hydrothermal alteration of Icelandic basalts. Solid solutions of chlorites, amphiboles and plagioclase were included as the secondary mineral phases. Saline and fresh water compositions of geothermal fluids were used to investigate the effect of salinity on alteration. Fluid-rock interaction simulations produce the observed mineral transformations. They show that roughly the same alteration minerals are formed due to reactions with both types of fluid which is in agreement with the field observations. A final application is directed towards the remediation of nitrate rich groundwaters. Removal of excess nitrate from groundwater by pyrite oxidation was modeled using the reactive fluid transport algorithm. Model results show that, when a pyrite-bearing, permeable zone is placed in the flow path, nitrate concentration in infiltrating water can be significantly lowered, in agreement with proposals from the literature. This is due to nitrogen reduction. Several simulations investigate the efficiency of systems with different mineral reactive surface areas, reactive barrier zone widths, and flow rates to identify the optimum setup.

Computer simulations of microstructures related to the olivine spinel transition

Tiefherd-Beben, die im oberen Erdmantel in einer Tiefe von ca. 400 km auftreten, werden gewöhnlich mit dem in gleicher Tiefe auftretenden druckabhängigen, polymorphen Phasenübergang von Olivine (α-Phase) zu Spinel (β-Phase) in Verbindung gebracht. Es ist jedoch nach wie vor unklar, wie der Phasenübergang mit dem mechanischen Versagen des Mantelmaterials zusammenhängt. Zur Zeit werden im Wesentlichen zwei Modelle diskutiert, die entweder Mikrostrukturen, die durch den Phasenübergang entstehen, oder aber die rheologischen Veränderungen des Mantelgesteins durch den Phasenübergang dafür verantwortlich machen. Dabei sind Untersuchungen der Olivin→Spinel Umwandlung durch die Unzugänglichkeit des natürlichen Materials vollständig auf theoretische Überlegungen sowie Hochdruck-Experimente und Numerische Simulationen beschränkt. Das zentrale Thema dieser Dissertation war es, ein funktionierendes Computermodell zur Simulation der Mikrostrukturen zu entwickeln, die durch den Phasenübergang entstehen. Des Weiteren wurde das Computer Modell angewandt um die mikrostrukturelle Entwicklung von Spinelkörnern und die Kontrollparameter zu untersuchen. Die Arbeit ist daher in zwei Teile unterteilt: Der erste Teil (Kap. 2 und 3) behandelt die physikalischen Gesetzmäßigkeiten und die prinzipielle Funktionsweise des Computer Modells, das auf der Kombination von Gleichungen zur Errechnung der kinetischen Reaktionsgeschwindigkeit mit Gesetzen der Nichtgleichgewichtsthermodynamik unter nicht-hydostatischen Bedingungen beruht. Das Computermodell erweitert ein Federnetzwerk der Software latte aus dem Programmpaket elle. Der wichtigste Parameter ist dabei die Normalspannung auf der Kornoberfläche von Spinel. Darüber hinaus berücksichtigt das Programm die Latenzwärme der Reaktion, die Oberflächenenergie und die geringe Viskosität von Mantelmaterial als weitere wesentliche Parameter in der Berechnung der Reaktionskinetic. Das Wachstumsverhalten und die fraktale Dimension von errechneten Spinelkörnern ist dabei in guter Übereinstimmung mit Spinelstrukturen aus Hochdruckexperimenten. Im zweiten Teil der Arbeit wird das Computermodell angewandt, um die Entwicklung der Oberflächenstruktur von Spinelkörnern unter verschiedenen Bedigungen zu eruieren. Die sogenannte ’anticrack theory of faulting’, die den katastrophalen Verlauf der Olivine→Spinel Umwandlung in olivinhaltigem Material unter differentieller Spannung durch Spannungskonzentrationen erklärt, wurde anhand des Computermodells untersucht. Der entsprechende Mechanismus konnte dabei nicht bestätigt werden. Stattdessen können Oberflächenstrukturen, die Ähnlichkeiten zu Anticracks aufweisen, durch Unreinheiten des Materials erklärt werden (Kap. 4). Eine Reihe von Simulationen wurde der Herleitung der wichtigsten Kontrollparameter der Reaktion in monomineralischem Olivin gewidmet (Kap. 5 and Kap. 6). Als wichtigste Einflüsse auf die Kornform von Spinel stellten sich dabei die Hauptnormalspannungen auf dem System sowie Heterogenitäten im Wirtsminerals und die Viskosität heraus. Im weiteren Verlauf wurden die Nukleierung und das Wachstum von Spinel in polymineralischen Mineralparagenesen untersucht (Kap. 7). Die Reaktionsgeschwindigkeit der Olivine→Spinel Umwandlung und die Entwicklung von Spinelnetzwerken und Clustern wird durch die Gegenwart nicht-reaktiver Minerale wie Granat oder Pyroxen erheblich beschleunigt. Die Bildung von Spinelnetzwerken hat das Potential, die mechanischen Eigenschaften von Mantelgestein erheblich zu beeinflussen, sei es durch die Bildung potentieller Scherzonen oder durch Gerüstbildung. Dieser Lokalisierungprozess des Spinelwachstums in Mantelgesteinen kann daher ein neues Erklärungsmuster für Tiefbeben darstellen.

Factorization method and inclusions of mixed type in an inverse elliptic boundary value problem

In various imaging problems the task is to use the Cauchy data of the solutions to an elliptic boundary value problem to reconstruct the coefficients of the corresponding partial differential equation. Often the examined object has known background properties but is contaminated by inhomogeneities that cause perturbations of the coefficient functions. The factorization method of Kirsch provides a tool for locating such inclusions. In this paper, the factorization technique is studied in the framework of coercive elliptic partial differential equations of the divergence type: Earlier it has been demonstrated that the factorization algorithm can reconstruct the support of a strictly positive (or negative) definite perturbation of the leading order coefficient, or if that remains unperturbed, the support of a strictly positive (or negative) perturbation of the zeroth order coefficient. In this work we show that these two types of inhomogeneities can, in fact, be located simultaneously. Unlike in the earlier articles on the factorization method, our inclusions may have disconnected complements and we also weaken some other a priori assumptions of the method. Our theoretical findings are complemented by two-dimensional numerical experiments that are presented in the framework of the diffusion approximation of optical tomography.

Sensitivity analysis of a parabolic-elliptic problem

We consider the heat flux through a domain with subregions in which the thermal capacity approaches zero. In these subregions the parabolic heat equation degenerates to an elliptic one. We show the well-posedness of such parabolic-elliptic differential equations for general non-negative L-infinity-capacities and study the continuity of the solutions with respect to the capacity, thus giving a rigorous justification for modeling a small thermal capacity by setting it to zero. We also characterize weak directional derivatives of the temperature with respect to capacity as solutions of related parabolic-elliptic problems.

Kinetic modeling and experiments on gas uptake and chemical transformation of organic aerosol in the atmosphere

Aerosolpartikel beeinflussen das Klima durch Streuung und Absorption von Strahlung sowie als Nukleations-Kerne für Wolkentröpfchen und Eiskristalle. Darüber hinaus haben Aerosole einen starken Einfluss auf die Luftverschmutzung und die öffentliche Gesundheit. Gas-Partikel-Wechselwirkunge sind wichtige Prozesse, weil sie die physikalischen und chemischen Eigenschaften von Aerosolen wie Toxizität, Reaktivität, Hygroskopizität und optische Eigenschaften beeinflussen. Durch einen Mangel an experimentellen Daten und universellen Modellformalismen sind jedoch die Mechanismen und die Kinetik der Gasaufnahme und der chemischen Transformation organischer Aerosolpartikel unzureichend erfasst. Sowohl die chemische Transformation als auch die negativen gesundheitlichen Auswirkungen von toxischen und allergenen Aerosolpartikeln, wie Ruß, polyzyklische aromatische Kohlenwasserstoffe (PAK) und Proteine, sind bislang nicht gut verstanden.rn Kinetische Fluss-Modelle für Aerosoloberflächen- und Partikelbulk-Chemie wurden auf Basis des Pöschl-Rudich-Ammann-Formalismus für Gas-Partikel-Wechselwirkungen entwickelt. Zunächst wurde das kinetische Doppelschicht-Oberflächenmodell K2-SURF entwickelt, welches den Abbau von PAK auf Aerosolpartikeln in Gegenwart von Ozon, Stickstoffdioxid, Wasserdampf, Hydroxyl- und Nitrat-Radikalen beschreibt. Kompetitive Adsorption und chemische Transformation der Oberfläche führen zu einer stark nicht-linearen Abhängigkeit der Ozon-Aufnahme bezüglich Gaszusammensetzung. Unter atmosphärischen Bedingungen reicht die chemische Lebensdauer von PAK von wenigen Minuten auf Ruß, über mehrere Stunden auf organischen und anorganischen Feststoffen bis hin zu Tagen auf flüssigen Partikeln. rn Anschließend wurde das kinetische Mehrschichtenmodell KM-SUB entwickelt um die chemische Transformation organischer Aerosolpartikel zu beschreiben. KM-SUB ist in der Lage, Transportprozesse und chemische Reaktionen an der Oberfläche und im Bulk von Aerosol-partikeln explizit aufzulösen. Es erforder im Gegensatz zu früheren Modellen keine vereinfachenden Annahmen über stationäre Zustände und radiale Durchmischung. In Kombination mit Literaturdaten und neuen experimentellen Ergebnissen wurde KM-SUB eingesetzt, um die Effekte von Grenzflächen- und Bulk-Transportprozessen auf die Ozonolyse und Nitrierung von Protein-Makromolekülen, Ölsäure, und verwandten organischen Ver¬bin-dungen aufzuklären. Die in dieser Studie entwickelten kinetischen Modelle sollen als Basis für die Entwicklung eines detaillierten Mechanismus für Aerosolchemie dienen sowie für das Herleiten von vereinfachten, jedoch realistischen Parametrisierungen für großskalige globale Atmosphären- und Klima-Modelle. rn Die in dieser Studie durchgeführten Experimente und Modellrechnungen liefern Beweise für die Bildung langlebiger reaktiver Sauerstoff-Intermediate (ROI) in der heterogenen Reaktion von Ozon mit Aerosolpartikeln. Die chemische Lebensdauer dieser Zwischenformen beträgt mehr als 100 s, deutlich länger als die Oberflächen-Verweilzeit von molekularem O3 (~10-9 s). Die ROIs erklären scheinbare Diskrepanzen zwischen früheren quantenmechanischen Berechnungen und kinetischen Experimenten. Sie spielen eine Schlüsselrolle in der chemischen Transformation sowie in den negativen Gesundheitseffekten von toxischen und allergenen Feinstaubkomponenten, wie Ruß, PAK und Proteine. ROIs sind vermutlich auch an der Zersetzung von Ozon auf mineralischem Staub und an der Bildung sowie am Wachstum von sekundären organischen Aerosolen beteiligt. Darüber hinaus bilden ROIs eine Verbindung zwischen atmosphärischen und biosphärischen Mehrphasenprozessen (chemische und biologische Alterung).rn Organische Verbindungen können als amorpher Feststoff oder in einem halbfesten Zustand vorliegen, der die Geschwindigkeit von heterogenen Reaktionenen und Mehrphasenprozessen in Aerosolen beeinflusst. Strömungsrohr-Experimente zeigen, dass die Ozonaufnahme und die oxidative Alterung von amorphen Proteinen durch Bulk-Diffusion kinetisch limitiert sind. Die reaktive Gasaufnahme zeigt eine deutliche Zunahme mit zunehmender Luftfeuchte, was durch eine Verringerung der Viskosität zu erklären ist, bedingt durch einen Phasenübergang der amorphen organischen Matrix von einem glasartigen zu einem halbfesten Zustand (feuchtigkeitsinduzierter Phasenübergang). Die chemische Lebensdauer reaktiver Verbindungen in organischen Partikeln kann von Sekunden bis zu Tagen ansteigen, da die Diffusionsrate in der halbfesten Phase bei niedriger Temperatur oder geringer Luftfeuchte um Größenordnungen absinken kann. Die Ergebnisse dieser Studie zeigen wie halbfeste Phasen die Auswirkung organischeer Aerosole auf Luftqualität, Gesundheit und Klima beeinflussen können. rn

Drainage of thin aqueous films between solid surfaces measured with the colloidal probe technique

Understanding liquid flow at the vicinity of solid surfaces is crucial to the developmentrnof technologies to reduce drag. One possibility to infer flow properties at the liquid-solid interface is to compare the experimental results to solutions of the Navier-Stokes equations assuming the no-slip boundary condition (BC) or the slip BC. There is no consensus in the literature about which BC should be used to model the flow of aqueous solutions over hydrophilic surfaces. Here, the colloidal probe technique is used to systematically address this issue, measuring forces acting during drainage of water over a surface. Results show that experimental variables, especially the cantilever spring constant, lead to the discrepancy observed in the literature. Two different parameters, calculated from experimental variables, could be used to separate the data obtained in this work and those reported in the literature in two groups: one explained with the no-slip BC, and another with the slip BC. The observed residual slippage is a function of instrumental variables, showing a trend incompatible with the available physical justifications. As a result, the no-slip is the more appropriate BC. The parameters can be used to avoid situations where the no-slip BC is not satisfied.

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.

Computer simulation studies of vector spin glasses

Die vorliegende Doktorarbeit befasst sich mit klassischen Vektor-Spingläsern eine Art von ungeordneten Magneten - auf verschiedenen Gittertypen. Da siernbedeutsam für eine experimentelle Realisierung sind, ist ein theoretisches Verständnis von Spinglas-Modellen mit wenigen Spinkomponenten und niedriger Gitterdimension von großer Bedeutung. Da sich dies jedoch als sehr schwierigrnerweist, sind neue, aussichtsreiche Ansätze nötig. Diese Arbeit betrachtet daher den Limesrnunendlich vieler Spindimensionen. Darin entstehen mehrere Vereinfachungen im Vergleichrnzu Modellen niedriger Spindimension, so dass für dieses bedeutsame Problem Eigenschaften sowohl bei Temperatur Null als auch bei endlichen Temperaturenrnüberwiegend mit numerischen Methoden ermittelt werden. Sowohl hyperkubische Gitter als auch ein vielseitiges 1d-Modell werden betrachtet. Letzteres erlaubt es, unterschiedliche Universalitätsklassen durch bloßes Abstimmen eines einzigen Parameters zu untersuchen. "Finite-size scaling''-Formen, kritische Exponenten, Quotienten kritischer Exponenten und andere kritische Größen werden nahegelegt und mit numerischen Ergebnissen verglichen. Eine detaillierte Beschreibung der Herleitungen aller numerisch ausgewerteter Gleichungen wird ebenso angegeben. Bei Temperatur Null wird eine gründliche Untersuchung der Grundzustände und Defektenergien gemacht. Eine Reihe interessanter Größen wird analysiert und insbesondere die untere kritische Dimension bestimmt. Bei endlicher Temperatur sind der Ordnungsparameter und die Spinglas-Suszeptibilität über die numerisch berechnete Korrelationsmatrix zugänglich. Das Spinglas-Modell im Limes unendlich vieler Spinkomponenten kann man als Ausgangspunkt zur Untersuchung der natürlicheren Modelle mit niedriger Spindimension betrachten. Wünschenswert wäre natürlich ein Modell, das die Vorteile des ersten mit den Eigenschaften des zweiten verbände. Daher wird in Modell mit Anisotropie vorgeschlagen und getestet, mit welchem versucht wird, dieses Ziel zu erreichen. Es wird auf reizvolle Wege hingewiesen, das Modell zu nutzen und eine tiefergehende Beschäftigung anzuregen. Zuletzt werden sogenannte "real-space" Renormierungsgruppenrechnungen sowohl analytisch als auch numerisch für endlich-dimensionale Vektor-Spingläser mit endlicher Anzahl von Spinkomponenten durchgeführt. Dies wird mit einer zuvor bestimmten neuen Migdal-Kadanoff Rekursionsrelation geschehen. Neben anderen Größen wird die untere kritische Dimension bestimmt.

Dynamic wetting of complex liquids

A simple dependency between contact angle θ and velocity or surface tension has been predicted for the wetting and dewetting behavior of simple liquids. According to the hydrodynamic theory, this dependency was described by Cox and Voinov as θ ∼ Ca^(1/3) (Ca: Capillary number). For more complex liquids like surfactant solutions, this prediction is not directly given.rnHere I present a rotating drum setup for studying wetting/dewetting processes of surfactant solutions on the basis of velocity-dependent contact angle measurements. With this new setup I showed that surfactant solutions do not follow the predicted Cox-Voinov relation, but showed a stronger contact angle dependency on surface tension. All surfactants independent of their charge showed this difference from the prediction so that electrostatic interactions as a reason could be excluded. Instead, I propose the formation of a surface tension gradient close to the three-phase contact line as the main reason for the strong contact angle decrease with increasing surfactant concentration. Surface tension gradients are not only formed locally close to the three-phase contact line, but also globally along the air-liquid interface due to the continuous creation/destruction of the interface by the drum moving out of/into the liquid. By systematically hindering the equilibration routes of the global gradient along the interface and/or through the bulk, I was able to show that the setup geometry is also important for the wetting/dewetting of surfactant solutions. Further, surface properties like roughness or chemical homogeneity of the wetted/dewetted substrate influence the wetting/dewetting behavior of the liquid, i. e. the three-phase contact line is differently pinned on rough/smooth or homogeneous/inhomogeneous surfaces. Altogether I showed that the wetting/dewetting of surfactant solutions did not depend on the surfactant type (anionic, cationic, or non-ionic) but on the surfactant concentration and strength, the setup geometry, and the surface properties.rnSurfactants do not only influence the wetting/dewetting behavior of liquids, but also the impact behavior of drops on free-standing films or solutions. In a further part of this work, I dealt with the stability of the air cushion between drop and film/solution. To allow coalescence between drop and substrate, the air cushion has to vanish. In the presence of surfactants, the vanishing of the air is slowed down due to a change in the boundary condition from slip to no-slip, i. e. coalescence is suppressed or slowed down in the presence of surfactant.

Molecular plating of thin lanthanide layers with improved material properties for nuclear applications

In der vorliegenden Arbeit werden Experimente beschrieben, die zu einem vertieften Verständnis fundamentaler Prozesse bei der elektrochemischen Herstellung von Dünnschichten, sog. Targets, für kernphysikalische und -chemische Studien führten. Targets wurden mittels 'Molecular Plating' (MP) hergestellt, indem eine Elektrodeposition aus organischem Medium in der Regel bei konstantem Strom in Zwei-Elektroden-Zellen. Die Resultate erlaubten, optimierte Herstellungs-bedingungen zu ermitteln, welche die Produktion deutlich verbesserter Targets erlaubten. MP bei konstantem Strom ist ein massentransportkontrollierter Prozess. Der angelegte Strom wird durch einen konstanten Fluss elektroaktiver Spezies zur Kathode – auf der die Schicht wächst – und Anode aufrechterhalten. Die Untersuchungen zeigten, dass das Zellenpotential des Elektrodepositionsystems immer durch den Ohm'schen Spannungsabfall auf Grund des Widerstandes der verwendeten Lösung dominiert wurde. Dies erlaubte die Herleitung einer Beziehung zwischen dem Zellenpotential und der Konzentration der elektroaktiven Spezies. Die Beziehung erlaubt die Erklärung des gemessenen zeitlichen Verlaufs des Zellenpotentials während der Abscheidung als Funktion der Elektrolytkonzentration. Dies dient als Basis, auf der nun ein umfassenderes Bild der Prozesse, die für die charakteristischen Minima im Potentialverlauf einer Abscheidung verantwortlich sind, gewonnen werden kann. Es konnte gezeigt werden, dass die Minima mit der fast vollständigen Entfernung (durch Abscheidung) der aus einem gelösten Salz erzeugten Nd-Ionen korrespondieren. Die abgeschiedene Spezies wurde als Nd3+ identifiziert, vermutlich als Carboxylat, Oxid oder Hydroxid, was auf Grund der hohen negative Werte des Standardredoxpotentials der Lanthanide verständlich erscheint. Von den vorliegenden elektroaktiven Spezies tragen die Nd3+ Ionen nur zu knapp 20% zum Gesamtstrom bei. Durch Elektrolyse tragen auch die Lösungsmittelkomponenten zu diese Strom bei. Die Gegenwart von elektrolysiertem Lösungsmittel wurde in Analysen der Dünnschichten bestätigt. Diese waren immer mit chemi- und physisorbierten Lösungsmittelmolekülen bedeckt. Die Analyse der Dünnschichten zeigte, dass die Oberflächen von einem furchenartiges Netz durchzogen waren, und dass diese während des Trocknen der Schichten nach dem MP entstanden. Ob die Schichten an Luft oder in inerter Atmosphäre trockneten, hatte keinen Einfluss. Es wurden Experimente mit mehreren Lösungsmitteln durchgeführt, die sich deutlich in ihren physikalischen Eigenschaften, v.a. dem Siedepunkt, unterschieden. Furchenfreie Dünnschichten konnten insbesondere bei MP in N,N-dimethylformamide (DMF) erzeugt werden. Die Verwendung von DMF in Kombination mit einer Abscheidung auf sehr glatten Substraten erlaubte die Produktion von sehr homogenen, glatten und defektfreien Schichten. Diese waren vermutlich geringeren inneren Spannungen während des Trocknens ausgesetzt, als Schichten auf raueren Substraten oder solche, die aus flüchtigeren Lösungsmitteln hergestellt wurden. Die Oberflächenrauigkeit des Substrats und das gewählte Lösungsmittel wurden so als Schlüsselfaktoren für die Produktion hochqualitativer Schichten identifiziert. Es konnte gezeigt werden, dass mit MP eine sehr effiziente Methode zur Herstellung homogener Schichten mit exzellenter Ausbeute ist. In weiteren Experimenten mit dem primordialen Alpha-Emitter 147Sm als Modellisotop wurde die Eignung solcher Schichten als Alpha-Quelle untersucht. Sowohl die Energieauflösung als auch der Anteil der Alpha-Teilchen, die den Detektor erreichten, waren von den Quelleneigenschaften abhängig. Die Effekte wurden verschiedenen Variablen der Dünnschicht zugeordnet, welche die Alpha-Spektren beeinflussten. Dominant war die Wahl des Lösungsmittels und die Rauigkeit des Substrats. Dies beeinflusste Schichtdicke und -morphologie sowie die Art des Schichtwachstums und veränderte die Detektionseffizienz in Alpha-Messungen bis zu 15%. Nur homogene, ebene Schichten, die aus DMF auf glatten Substraten abgeschieden wurden, eignen sich optimal als Alpha-Quelle. Die gewonnenen Ergebnisse erlauben die optimierte Herstellung nuklearer Targets durch MP. Künftige Anwendungen beinhalten insbesondere die Herstellung von Targets für neutroneninduzierte Spaltexperimente und untergrundarmeAlpha-Messungen sehr kleiner Aktivitäten.

Design of multifunctional magnetic nanomaterials for biomedical applications

Diese Arbeit ist ein Beitrag zu den schnell wachsenden Forschungsgebieten der Nano-Biotechnologie und Nanomedizin. Sie behandelt die spezifische Gestaltung magnetischer Nanomaterialien für verschiedene biomedizinische Anwendungsgebiete, wie beispielsweise Kontrastmittel für die magnetische Resonanztomographie (MRT) oder "theragnostische" Agenzien für simultane optische/MR Detektion und Behandlung mittels photodynamischer Therapie (PDT).rnEine Vielzahl magnetischer Nanopartikel (NP) mit unterschiedlichsten magnetischen Eigenschaften wurden im Rahmen dieser Arbeit synthetisiert und erschöpfend charakterisiert. Darüber hinaus wurde eine ganze Reihe von Oberflächenmodifizierungsstrategien entwickelt, um sowohl die kolloidale als auch die chemische Stabilität der Partikel zu verbessern, und dadurch den hohen Anforderungen der in vitro und in vivo Applikation gerecht zu werden. Diese Strategien beinhalteten nicht nur die Verwendung bi-funktionaler und multifunktioneller Polymerliganden, sondern auch die Kondensation geeigneter Silanverbindungen, um eine robuste, chemisch inerte und hydrophile Siliziumdioxid- (SiO2) Schale um die magnetischen NP auszubilden.rnGenauer gesagt, der Bildungsmechanismus und die magnetischen Eigenschaften monodisperser MnO NPs wurden ausgiebig untersucht. Aufgrund ihres einzigartigen magnetischen Verhaltens eignen sich diese NPs besonders als (positive) Kontrastmittel zur Verkürzung der longitudinalen Relaxationszeit T1, was zu einer Aufhellung im entsprechenden MRT-Bild führt. Tatsächlich wurde dieses kontrastverbessernde Potential in mehreren Studien mit unterschiedlichen Oberflächenliganden bestätigt. Au@MnO „Nanoblumen“, auf der anderen Seite, sind Vertreter einer weiteren Klasse von Nanomaterialien, die in den vergangenen Jahren erhebliches Interesse in der wissenschaftlichen Welt geweckt hat und oft „Nano-hetero-Materialien“ genannt wird. Solche Nano-hetero-partikel vereinen die individuellen physikalischen und chemischen Eigenschaften der jeweiligen Komponenten in einem nanopartikulärem System und erhöhen dadurch die Vielseitigkeit der möglichen Anwendungen. Sowohl die magnetischen Merkmale von MnO, als auch die optischen Eigenschaften von Au bieten die Möglichkeit, diese „Nanoblumen“ für die kombinierte MRT und optische Bildgebung zu verwenden. Darüber hinaus erlaubt das Vorliegen zweier chemisch unterschiedlicher Oberflächen die gleichzeitige selektive Anbindung von Katecholliganden (auf MnO) und Thiolliganden (auf Au). Außerdem wurde das therapeutische Potential von magnetischen NPs anhand von MnO NPs demonstriert, die mit dem Photosensibilisator Protoporhyrin IX (PP) funktionalisiert waren. Bei Bestrahlung mit sichtbarem Licht initiiert PP die Produktion von zytotoxisch-reaktivem Sauerstoff. Wir zeigen, dass Nierenkrebszellen, die mit PP-funktionalisierten MnO NPs inkubiert wurden nach Bestrahlung mit Laserlicht verenden, während sie ohne Bestrahlung unverändert bleiben. In einem ähnlichen Experiment untersuchten wir die Eigenschaften von SiO2 beschichteten MnO NPs. Dafür wurde eigens eine neuartige SiO2-Beschichtungsmethode entwickelt, die einer nachfolgende weitere Anbindung verschiedenster Liganden und die Einlagerung von Fluoreszenzfarbstoffen durch herkömmliche Silan- Sol-Gel Chemie erlaubt. Die Partikel zeigten eine ausgezeichnete Stabilität in einer ganzen Reihe wässriger Lösungen, darunter auch physiologische Kochsalzlösung, Pufferlösungen und humanes Blutserum, und waren weniger anfällig gegenüber Mn-Ionenauswaschung als einfache PEGylierte MnO NPs. Des Weiteren konnte bewiesen werden, dass die dünne SiO2 Schicht nur einen geringen Einfluss auf das magnetische Verhalten der NPs hatte, so dass sie weiterhin als T1-Kontrastmittel verwendet werden können. Schließlich konnten zusätzlich FePt@MnO NPs hergestellt werden, welche die individuellen magnetischen Merkmale eines ferromagnetischen (FePt) und eines antiferromagnetischen (MnO) Materials vereinen. Wir zeigen, dass wir die jeweiligen Partikelgrößen, und damit das resultierende magnetische Verhalten, durch Veränderung der experimentellen Parameter variieren können. Die magnetische Wechselwirkung zwischen beiden Materialien kann dabei auf Spinkommunikation an der Grenzfläche zwischen beiden NP-Sorten zurückgeführt werden.rn

Numerical simulation of some viscoelastic fluids

Numerical simulation of the Oldroyd-B type viscoelastic fluids is a very challenging problem. rnThe well-known High Weissenberg Number Problem" has haunted the mathematicians, computer scientists, and rnengineers for more than 40 years. rnWhen the Weissenberg number, which represents the ratio of elasticity to viscosity, rnexceeds some limits, simulations done by standard methods break down exponentially fast in time. rnHowever, some approaches, such as the logarithm transformation technique can significantly improve rnthe limits of the Weissenberg number until which the simulations stay stable. rnrnWe should point out that the global existence of weak solutions for the Oldroyd-B model is still open. rnLet us note that in the evolution equation of the elastic stress tensor the terms describing diffusive rneffects are typically neglected in the modelling due to their smallness. However, when keeping rnthese diffusive terms in the constitutive law the global existence of weak solutions in two-space dimension rncan been shown. rnrnThis main part of the thesis is devoted to the stability study of the Oldroyd-B viscoelastic model. rnFirstly, we show that the free energy of the diffusive Oldroyd-B model as well as its rnlogarithm transformation are dissipative in time. rnFurther, we have developed free energy dissipative schemes based on the characteristic finite element and finite difference framework. rnIn addition, the global linear stability analysis of the diffusive Oldroyd-B model has also be discussed. rnThe next part of the thesis deals with the error estimates of the combined finite element rnand finite volume discretization of a special Oldroyd-B model which covers the limiting rncase of Weissenberg number going to infinity. Theoretical results are confirmed by a series of numerical rnexperiments, which are presented in the thesis, too.

Analysis and numerical solution of the Peterlin viscoelastic model

Liquids and gasses form a vital part of nature. Many of these are complex fluids with non-Newtonian behaviour. We introduce a mathematical model describing the unsteady motion of an incompressible polymeric fluid. Each polymer molecule is treated as two beads connected by a spring. For the nonlinear spring force it is not possible to obtain a closed system of equations, unless we approximate the force law. The Peterlin approximation replaces the length of the spring by the length of the average spring. Consequently, the macroscopic dumbbell-based model for dilute polymer solutions is obtained. The model consists of the conservation of mass and momentum and time evolution of the symmetric positive definite conformation tensor, where the diffusive effects are taken into account. In two space dimensions we prove global in time existence of weak solutions. Assuming more regular data we show higher regularity and consequently uniqueness of the weak solution. For the Oseen-type Peterlin model we propose a linear pressure-stabilized characteristics finite element scheme. We derive the corresponding error estimates and we prove, for linear finite elements, the optimal first order accuracy. Theoretical error of the pressure-stabilized characteristic finite element scheme is confirmed by a series of numerical experiments.