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.

Petrographie und Geochemie des Eckfelder Maares und von Hillscheid und die Einbindung in den Hocheifelvulkanismus

Tertiäre Vulkanite aus dem Eckfelder Maar, dem Hillscheider Diatrem und dem Hillscheid Basalt (Schlot) wurden petrologisch und geochemisch untersucht. Bis auf tonige Klasten aus dem Bohrkern des Eckfelder Maares handelt es sich bei allen weiteren Proben um undifferenzierte basische Vulkanite. Die tonigen Klasten aus dem Bohrkern müssen der ehemaligen Landoberfläche vor der Eruption des Eckfelder Maares zugerechnet werden, in dessen Krater sie während der Eruption hineingefallen sind. Bis auf die Proben des Hillscheid Basaltes sind die Proben alteriert. Die Alteration zeigt sich an der Bildung von Zeolithen und Calcitmineralisationen, die primär und sekundär gebildete Hohlräume aufgefüllt haben oder an einer vertonten Grundmasse der Proben, die daneben Mineraleinschlüsse (Spinell) und kantige Fremdgesteinsbruchstücke enthalten können. Bei den Proben mit vertonter Grundmasse handelt es sich um Palagonite, Umwandlungsprodukte aus Sideromelan (basaltischem Glas). Geochemische Analysen an Grundmassepräparaten der alterierten bis vertonten Proben zeigen, dass außer den immobilen Elementen Ti, Nb, Zr, Y alle weiteren Elemente teilweise bis vollständig abgereichert worden sind. Eine Ausnahme bildet Barium (Ba), welches z.T. in beträchtlichen Mengen in Zeolithen (Harmotom) angereichert wurde. Bei den Proben aus dem Eckfelder Maar kann die Alteration bis Vertonung der Proben alleine mit der Palagonitisierung und Verwitterung erklärt werden. Es gibt keine Hinweise auf Materialzufuhr und damit für sich anschließende hydrothermale Prozesse. Die Proben des Hillscheider Diatrem sind wesentlich geringer alteriert (glasige Grundmasse). Neuste Erkenntnisse aus einer Bohrung im Sommer 1999 im vermuteten Zentrum des Hillscheider Diatrems beschränken das Diatrem maximal auf einen kleineren Bereich im Nordosten der bisherigen Lokation. Bei der Bohrung stieß man nach 20 Meter auf Anstehendes. Im Hangschutt darüber fand man Blöcke des Hillscheid Basaltes. Eine geringere Größe der Lokation zusammen mit der geringen Alteration könnten auf deren Entstehung mit einer initialen Maarphase gefolgt von Schlackentätigkeit hinweisen. Die Schlacken könnten die ersten Ablagerungen vor Verwitterung geschützt haben. Allerdings gibt es keine Funde die eine Schlackentätigkeit belegen. Beim sogenannten 'Hillscheider Diatrem' könnte es sich aber auch um Hangschutt aus der Randbreccie des Hillscheider Basaltes handeln. Zusammen mit Bruchstücken aus dem Schlot des Hillscheid Basaltes wären die Palagonite des sogenannten 'Hillscheider Diatrem' erst in jüngster Zeit im Bereich einer Uferböschung zur Ablagerung gekommen. Dies würde allerdings das sogenannte 'Hillscheider Diatrem' in seiner Existenz in Frage stellen. Vergleiche der Proben des Hillscheid Basaltes mit basischen Hocheifelvulkaniten deuten auf kogenetische Beziehung aller Proben untereinander und ordnen die Proben des Hillscheid Basaltes geochemisch dem Hocheifelvulkanismus zu. REE- und weitere Spurenelementgehalte und auch deren Elementverhältnisse weisen für alle tertiären Eifelvulkanite gemeinsam auf Mantelschmelzen aus dem Bereich eines Granatperidotits mit niedrigen Aufschmelzgraden um ein Prozent hin. Vergleich der Elementverhältnisse hochinkompatibler Elemente im Bezug auf die Bildung mafischer Schmelzen mit primitivem Mantel deuten darauf hin, dass der Mantel im Bereich der Hocheifel verarmt ist an K, Rb, Sr und angereichert an Ba und eventuell an Nb. Ursachen für diese von typisch primären Mantelzusammensetzung abweichenden Verhältnisse könnten durch Mischungen von Mantelschmelzen mit lithosphärischem Mantel (K-Anomalie) und durch Anreicherungen mit fluiden Phasen (Ba-Anomalie) oder auch Schmelzen aus einem tieferliegenden Plume (Kelberger Hoch) verursacht worden sein. Englischer Zusammenfassung:

Computer simulations of colloidal fluids in confinement

Computer-Simulationen von Kolloidalen Fluiden in Beschränkten Geometrien Kolloidale Suspensionen, die einen Phasenübergang aufweisen, zeigen eine Vielfalt an interessanten Effekten, sobald sie auf eine bestimmte Geometrie beschränkt werden, wie zum Beispiel auf zylindrische Poren, sphärische Hohlräume oder auf einen Spalt mit ebenen Wänden. Der Einfluss dieser verschiedenen Geometrietypen sowohl auf das Phasenverhalten als auch auf die Dynamik von Kolloid-Polymer-Mischungen wird mit Hilfe von Computer-Simulationen unter Verwendung des Asakura-Oosawa- Modells, für welches auf Grund der “Depletion”-Kräfte ein Phasenübergang existiert, untersucht. Im Fall von zylindrischen Poren sieht man ein interessantes Phasenverhalten, welches vom eindimensionalen Charakter des Systems hervorgerufen wird. In einer kurzen Pore findet man im Bereich des Phasendiagramms, in dem das System typischerweise entmischt, entweder eine polymerreiche oder eine kolloidreiche Phase vor. Sobald aber die Länge der zylindrischen Pore die typische Korrelationslänge entlang der Zylinderachse überschreitet, bilden sich mehrere quasi-eindimensionale Bereiche der polymerreichen und der kolloidreichen Phase, welche von nun an koexistieren. Diese Untersuchungen helfen das Verhalten von Adsorptionshysteresekurven in entsprechenden Experimenten zu erklären. Wenn das Kolloid-Polymer-Modellsystem auf einen sphärischen Hohlraum eingeschränkt wird, verschiebt sich der Punkt des Phasenübergangs von der polymerreichen zur kolloidreichen Phase. Es wird gezeigt, dass diese Verschiebung direkt von den Benetzungseigenschaften des Systems abhängt, was die Beobachtung von zwei verschiedenen Morphologien bei Phasenkoexistenz ermöglicht – Schalenstrukturen und Strukturen des Janustyps. Im Rahmen der Untersuchung von heterogener Keimbildung von Kristallen innerhalb einer Flüssigkeit wird eine neue Simulationsmethode zur Berechnung von Freien Energien der Grenzfläche zwischen Kristall- bzw. Flüssigkeitsphase undWand präsentiert. Die Resultate für ein System von harten Kugeln und ein System einer Kolloid- Polymer-Mischung werden anschließend zur Bestimmung von Kontaktwinkeln von Kristallkeimen an Wänden verwendet. Die Dynamik der Phasenseparation eines quasi-zweidimensionalen Systems, welche sich nach einem Quench des Systems aus dem homogenen Zustand in den entmischten Zustand ausbildet, wird mit Hilfe von einer mesoskaligen Simulationsmethode (“Multi Particle Collision Dynamics”) untersucht, die sich für eine detaillierte Untersuchung des Einflusses der hydrodynamischen Wechselwirkung eignet. Die Exponenten universeller Potenzgesetze, die das Wachstum der mittleren Domänengröße beschreiben, welche für rein zwei- bzw. dreidimensionale Systeme bekannt sind, können für bestimmte Parameterbereiche nachgewiesen werden. Die unterschiedliche Dynamik senkrecht bzw. parallel zu den Wänden sowie der Einfluss der Randbedingungen für das Lösungsmittel werden untersucht. Es wird gezeigt, dass die daraus resultierende Abschirmung der hydrodynamischen Wechselwirkungsreichweite starke Auswirkungen auf das Wachstum der mittleren Domänengröße hat.

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.