IMEX finite volume methods for the shallow water equations

Die Flachwassergleichungen (SWE) sind ein hyperbolisches System von Bilanzgleichungen, die adäquate Approximationen an groß-skalige Strömungen der Ozeane, Flüsse und der Atmosphäre liefern. Dabei werden Masse und Impuls erhalten. Wir unterscheiden zwei charakteristische Geschwindigkeiten: die Advektionsgeschwindigkeit, d.h. die Geschwindigkeit des Massentransports, und die Geschwindigkeit von Schwerewellen, d.h. die Geschwindigkeit der Oberflächenwellen, die Energie und Impuls tragen. Die Froude-Zahl ist eine Kennzahl und ist durch das Verhältnis der Referenzadvektionsgeschwindigkeit zu der Referenzgeschwindigkeit der Schwerewellen gegeben. Für die oben genannten Anwendungen ist sie typischerweise sehr klein, z.B. 0.01. Zeit-explizite Finite-Volume-Verfahren werden am öftersten zur numerischen Berechnung hyperbolischer Bilanzgleichungen benutzt. Daher muss die CFL-Stabilitätsbedingung eingehalten werden und das Zeitinkrement ist ungefähr proportional zu der Froude-Zahl. Deswegen entsteht bei kleinen Froude-Zahlen, etwa kleiner als 0.2, ein hoher Rechenaufwand. Ferner sind die numerischen Lösungen dissipativ. Es ist allgemein bekannt, dass die Lösungen der SWE gegen die Lösungen der Seegleichungen/ Froude-Zahl Null SWE für Froude-Zahl gegen Null konvergieren, falls adäquate Bedingungen erfüllt sind. In diesem Grenzwertprozess ändern die Gleichungen ihren Typ von hyperbolisch zu hyperbolisch.-elliptisch. Ferner kann bei kleinen Froude-Zahlen die Konvergenzordnung sinken oder das numerische Verfahren zusammenbrechen. Insbesondere wurde bei zeit-expliziten Verfahren falsches asymptotisches Verhalten (bzgl. der Froude-Zahl) beobachtet, das diese Effekte verursachen könnte.Ozeanographische und atmosphärische Strömungen sind typischerweise kleine Störungen eines unterliegenden Equilibriumzustandes. Wir möchten, dass numerische Verfahren für Bilanzgleichungen gewisse Equilibriumzustände exakt erhalten, sonst können künstliche Strömungen vom Verfahren erzeugt werden. Daher ist die Quelltermapproximation essentiell. Numerische Verfahren die Equilibriumzustände erhalten heißen ausbalanciert.rnrnIn der vorliegenden Arbeit spalten wir die SWE in einen steifen, linearen und einen nicht-steifen Teil, um die starke Einschränkung der Zeitschritte durch die CFL-Bedingung zu umgehen. Der steife Teil wird implizit und der nicht-steife explizit approximiert. Dazu verwenden wir IMEX (implicit-explicit) Runge-Kutta und IMEX Mehrschritt-Zeitdiskretisierungen. Die Raumdiskretisierung erfolgt mittels der Finite-Volumen-Methode. Der steife Teil wird mit Hilfe von finiter Differenzen oder au eine acht mehrdimensional Art und Weise approximniert. Zur mehrdimensionalen Approximation verwenden wir approximative Evolutionsoperatoren, die alle unendlich viele Informationsausbreitungsrichtungen berücksichtigen. Die expliziten Terme werden mit gewöhnlichen numerischen Flüssen approximiert. Daher erhalten wir eine Stabilitätsbedingung analog zu einer rein advektiven Strömung, d.h. das Zeitinkrement vergrößert um den Faktor Kehrwert der Froude-Zahl. Die in dieser Arbeit hergeleiteten Verfahren sind asymptotisch erhaltend und ausbalanciert. Die asymptotischer Erhaltung stellt sicher, dass numerische Lösung das "korrekte" asymptotische Verhalten bezüglich kleiner Froude-Zahlen besitzt. Wir präsentieren Verfahren erster und zweiter Ordnung. Numerische Resultate bestätigen die Konvergenzordnung, so wie Stabilität, Ausbalanciertheit und die asymptotische Erhaltung. Insbesondere beobachten wir bei machen Verfahren, dass die Konvergenzordnung fast unabhängig von der Froude-Zahl ist.

Methoden und Modelle zur Phytoplanktonanalyse am Beispiel einer Restaurierungsmaßnahme an der Bleiloch Talsperre

Methoden und Modelle zur Phytoplanktonanalyse am Beispiel einer Restaurierungsmaßnahme an der Bleiloch Talsperre Zusammenfassung An der hypertrophen Bleiloch Talsperre (Thüringen, BRD) werden die Auswirkungen der 1997 in Betrieb genommenen partiellen Destratifikation (0- 20m) untersucht. Ziel der Maßnahme war die Hemmung des Phytoplanktonwachstums durch Lichtlimitation.Da es nicht gelang, auch die Oberflächenlammelle in die Durchmischung einzubeziehen, wurde auch keine Reduktion des Phytoplanktons beobachtet. Seine zeitliche Dynamik nahm deutlich zu und es traten neben typischen Längsgradienten auch extreme, kleinräumige Konzentrationsunterschiede auf. Die vorher häufigen Cyanobakterien spielten jedoch keine quantitativ bedeutsame Rolle mehr. Es wird diskutiert, wie die Belüftung die Saisonalität des Planktons stört und zu Verschiebungen der Zusammensetzung führt. Der Zusammenhang zwischen der Chlorophyll a Menge und der Stabilität der Wassersäule sowie der Einfluß wachstumsbestimmender Größen werden modelliert. Methodische Schwerpunkte umfassen die Auswertung von Phytoplankton-Pigmentdaten (HPLC-Analyse) und die Entwicklung einer Durchflußzytometer-Methode. Bei letzerer werden an lebenden Phytoplanktern die Vorwärtsstreuung sowie vier Fluoreszenzen (2 Phycoerythrine, Phycocyanin und Chlorophyll a, 514nm Anregungswellenlänge) gemessen. Zur Klassifizierung von Algen anhand von Referenzen werden zwei Methoden entwickelt, Möglichkeiten des Sortierens aufgezeigt und es wird ein Vergleich mit Mikroskopie und HPLC- Pigmentanalyse durchgeführt. Anhand von Modellrechnungen wird demonstriert, wie sich systematische bzw. zufällige Abweichungen der Pigmentstöchiometrien auf die Berechnung von Algengruppenanteilen auswirken. Strategien zur Auswertung von Pigmentdaten werden empfohlen.

Phasenverhalten, Kristallisation und Morphologie ladungsstabilisierter kolloidaler Mischungen

Die vorliegende Arbeit möchte die Anwendbarkeit ladungsstabilisierter kolloidaler Systeme als Modellsysteme für fundamentale Fragen der Festkörperphysik und Thermodynamik auf binäre Mischungen erweitern. In diesem Kontext untersucht sie das Phasenverhalten und mit ihm im Zusammenhang stehende Eigenschaften von binären Mischungen ladungsstabilisierter, sphärischer kolloidaler Partikel in wässriger Suspension. Da das Verhalten hartkugelähnlicher Systeme durch hohe Fremdionenkonzentrationen bereits gut bekannt ist, konzentriert sich diese Arbeit auf sehr langreichweitig repulsive Systeme unter deionisierten Bedingungen. Neben etablierten Methoden der Mikroskopie und statischer Lichtstreuung zur Phasendiagrammsbestimmung wird auch die Beobachtung der zeitabhängigen Entwicklung des Schermoduls verwendet, um eine langsame Erstarrungskinetik zu studieren. Es werden insbesondere Mischungen aus Komponenten unterschiedlicher Größe und Ladung der Größenverhältnisse 0,9, 0,82, 0,57, 0,39 und 0,37 untersucht. Diese zeigen in dieser Reihenfolge Phasendiagramme mit spindelförmigem fluid/kristallinen Koexistenzbereich wie auch azeotrope und eutektische Phasendiagramme. Die Strukturuntersuchungen aus der statischen Lichtstreuung stehen in praktisch allen Fällen im Einklang mit ungeordneten bcc- Substitutionskristallen, was über Modelle zu Schermodulmessungen bestätigt wird. Für das spindelförmige System wird ein überraschend weiter Koexistenzbereich beobachtet, wie er nicht von der Theorie erwartet wird. Die Lage, aber nicht die Form des Solidus stimmt quantitativ mit Simulationsvorhersagen zu einkomponentigen Systemen überein. Für das eutektische System bei einem Radienverhältnis von 0,57 wird der Einfluss der Schwerkraft auf das Phasenverhalten und die Erstarrungskinetik untersucht. Die der Kristallisation der kleineren Majoritätskomponente vorgelagerte gravitativ unterstützte Entmischung begünstigt hier die Verfestigung. Beobachtet werden Morphologien, die aus anderen Systemen bekannt sind (Facetten, Dendriten), wie auch erstmals eine kollumnare eutektische Morphologie. Aus den Ergebnissen wird der erste umfassende Überblick über das Phasenverhalten deionisierter Mischungen ladungsstabilisierter, sphärischer Partikel erstellt, die eine Diskussion der Daten anderer Autoren und unserer Gruppe über fluid-fluider Phasenseparation und einem System mit oberem azeotropen Punkt mit einschließt. Die meisten metallspezifischen Phasendiagrammtypen können mit ladungsstabilisierten kolloidalen Partikeln reproduziert werden. Die langreichweitig wechselwirkenden Partikel zeigt eine wesentlich verbesserte substitutionelle Mischbarkeit im Vergleich mit Hartkugel- und Metallsystemen. Das Größenverhältnis der sphärischen Partikel nimmt dabei die bestimmende Rolle für den Phasendiagrammtyp ein.

Boron determination in biological samples - Intercomparison of three analytical methods to assist development of a treatment protocol for neoplastic diseases of the liver with Boron Neutron Capture Therapy

Die Bor-Neuroneneinfang-Therapie (engl.: Boron Neutron Capture Therapy, BNCT) ist eine indirekte Strahlentherapie, welche durch die gezielte Freisetzung von dicht ionisierender Strahlung Tumorzellen zerstört. Die freigesetzten Ionen sind Spaltfragmente einer Kernreaktion, bei welcher das Isotop 10B ein niederenergetisches (thermisches) Neutron einfängt. Das 10B wird durch ein spezielles Borpräparat in den Tumorzellen angereichert, welches selbst nicht radioaktiv ist. rnAn der Johannes Gutenberg-Universität Mainz wurde die Forschung für die Anwendung eines klinischen Behandlungsprotokolls durch zwei Heilversuche bei Patienten mit kolorektalen Lebermetastasen an der Universität Pavia, Italien, angeregt, bei denen die Leber außerhalb des Körpers in einem Forschungsreaktor bestrahlt wurde. Als erster Schritt wurde in Kooperation verschiedener universitärer Institute eine klinische Studie zur Bestimmung klinisch relevanter Parameter wie der Borverteilung in verschiedenen Geweben und dem pharmakokinetischen Aufnahmeverhalten des Borpräparates initiiert.rnDie Borkonzentration in den Gewebeproben wurde hinsichtlich ihrer räumlichen Verteilung in verschiedenen Zellarealen bestimmt, um mehr über das Aufnahmeverhalten der Zellen für das BPA im Hinblick auf ihre biologischen Charakteristika zu erfahren. Die Borbestimung wurde per Quantitative Neutron Capture Radiography, Prompt Gamma Activation Analysis und Inductively Coupled Plasma Mass Spectroscopy parallel zur histologischen Analyse des Gewebes durchgeführt. Es war möglich zu zeigen, dass in Proben aus Tumorgewebe und aus tumorfreiem Gewebe mit unterschiedlichen morphologischen Eigenschaften eine sehr heterogene Borverteilung vorliegt. Die Ergebnisse der Blutproben werden für die Erstellung eines pharmakokinetischen Modells verwendet und sind in Übereinstimmung mit existierenden pharmakokinetische Modellen. Zusätzlich wurden die Methoden zur Borbestimmung über speziell hergestellte Referenzstandards untereinander verglichen. Dabei wurde eine gute Übereinstimmung der Ergebnisse festgestellt, ferner wurde für alle biologischen Proben Standardanalyseprotokolle erstellt.rnDie bisher erhaltenen Ergebnisse der klinischen Studie sind vielversprechend, lassen aber noch keine endgültigen Schlussfolgerungen hinsichtlich der Wirksamkeit von BNCT für maligne Lebererkrankungen zu. rn

Development of models and methods to simulate peptide-assisted nucleation and growth of calcium-minerals

In den vergangenen Jahren wurden einige bislang unbekannte Phänomene experimentell beobachtet, wie etwa die Existenz unterschiedlicher Prä-Nukleations-Strukturen. Diese haben zu einem neuen Verständnis von Prozessen, die auf molekularer Ebene während der Nukleation und dem Wachstum von Kristallen auftreten, beigetragen. Die Auswirkungen solcher Prä-Nukleations-Strukturen auf den Prozess der Biomineralisation sind noch nicht hinreichend verstanden. Die Mechanismen, mittels derer biomolekulare Modifikatoren, wie Peptide, mit Prä-Nukleations-Strukturen interagieren und somit den Nukleationsprozess von Mineralen beeinflussen könnten, sind vielfältig. Molekulare Simulationen sind zur Analyse der Formation von Prä-Nukleations-Strukturen in Anwesenheit von Modifikatoren gut geeignet. Die vorliegende Arbeit beschreibt einen Ansatz zur Analyse der Interaktion von Peptiden mit den in Lösung befindlichen Bestandteilen der entstehenden Kristalle mit Hilfe von Molekular-Dynamik Simulationen.rnUm informative Simulationen zu ermöglichen, wurde in einem ersten Schritt die Qualität bestehender Kraftfelder im Hinblick auf die Beschreibung von mit Calciumionen interagierenden Oligoglutamaten in wässrigen Lösungen untersucht. Es zeigte sich, dass große Unstimmigkeiten zwischen etablierten Kraftfeldern bestehen, und dass keines der untersuchten Kraftfelder eine realistische Beschreibung der Ionen-Paarung dieser komplexen Ionen widerspiegelte. Daher wurde eine Strategie zur Optimierung bestehender biomolekularer Kraftfelder in dieser Hinsicht entwickelt. Relativ geringe Veränderungen der auf die Ionen–Peptid van-der-Waals-Wechselwirkungen bezogenen Parameter reichten aus, um ein verlässliches Modell für das untersuchte System zu erzielen. rnDas umfassende Sampling des Phasenraumes der Systeme stellt aufgrund der zahlreichen Freiheitsgrade und der starken Interaktionen zwischen Calciumionen und Glutamat in Lösung eine besondere Herausforderung dar. Daher wurde die Methode der Biasing Potential Replica Exchange Molekular-Dynamik Simulationen im Hinblick auf das Sampling von Oligoglutamaten justiert und es erfolgte die Simulation von Peptiden verschiedener Kettenlängen in Anwesenheit von Calciumionen. Mit Hilfe der Sketch-Map Analyse konnten im Rahmen der Simulationen zahlreiche stabile Ionen-Peptid-Komplexe identifiziert werden, welche die Formation von Prä-Nukleations-Strukturen beeinflussen könnten. Abhängig von der Kettenlänge des Peptids weisen diese Komplexe charakteristische Abstände zwischen den Calciumionen auf. Diese ähneln einigen Abständen zwischen den Calciumionen in jenen Phasen von Calcium-Oxalat Kristallen, die in Anwesenheit von Oligoglutamaten gewachsen sind. Die Analogie der Abstände zwischen Calciumionen in gelösten Ionen-Peptid-Komplexen und in Calcium-Oxalat Kristallen könnte auf die Bedeutung von Ionen-Peptid-Komplexen im Prozess der Nukleation und des Wachstums von Biomineralen hindeuten und stellt einen möglichen Erklärungsansatz für die Fähigkeit von Oligoglutamaten zur Beeinflussung der Phase des sich formierenden Kristalls dar, die experimentell beobachtet wurde.