Beiträge zur Modellierung technischer dynamischer Systeme

Die Nichtlineare Dynamik verallgemeinert Aussagen über dynamische Systeme durch Abstraktion von konkreten Systemen. In der Technik sind Maschinen dagegen sehr konkret und die Behandlung auftretender Probleme mit Methoden der theoretischen Physik ist nicht trivial. Diese Arbeit versucht einige der Schwierigkeiten einer technischen Anwendung der nichtlinearen Theorie zu lokalisieren. Am Beispiel von vier Klassen von Modellansätzen, werden Anwendungsschnittstellen beleuchtet und systematisiert. Die Anwendung von Modellen, die explizit auf bekannten physikalischen Gesetzmäßigkeiten aufbauen, findet Grenzen in der Anzahl der Freiheitsgrade und den Nebenbedingungen konkreter Systeme. Solche Modelle liefern jedoch wichtige Hinweise auf die Vielfalt der nichtlinearen Phänomene und tragen zu ihrem Verständnis bei. Daher sind sie für die Konstruktionspraxis wichtig. Es werden typisch nichtlineare Phänomene und ihre zugrundeliegenden Mechanismen vorgestellt und klassifiziert, sowie grundsätzliche Probleme der Berechenbarkeit analytisch formulierter Modelle betrachtet. Eine zweite Schnittstelle bieten die Darstellungen des Systemverhaltens als überlagerung spezieller Funktionen, diez.B. Symmetrieeigenschaften des betrachteten Systems besonders deutlich widerspiegeln. Gegenüber der klassischen Fourierzerlegung nach Frequenz und Phase bringt die Analyse nach Detaillierungsgrad und Position von Waveletfunktionen wichtige Vorteile für die nichtlineare zustandsraumbasierte Datenanalyse. Viele Verfahren der Nichtlinearen Datenanalyse beruhen auf metrischen Eigenschaften der dynamischen Systeme. Als dritte Gruppe werden demgegenüber topologische Methoden beleuchtet. Die Konstruktion von Simplexen aus Zeitreihen mittels der Zeitversatzmethode ist die Grundlage für eine Triangulation der Zustandsräume. Die Methoden, z.B. Templateverfahren, die auf der Einbettung von eindimensionalen Trajektorien in den R^3 basieren, lassen sich hingegen nicht einfach auf hochdimensionale Zustandsmannigfaltigkeiten anwenden. Schließlich werden stochastische Aspekte behandelt. Schwankungen des Systemverhaltens können auf Schwankungen der Anfangswerte und/oder auf Schwankungen der eigentlichen Systemdynamik beruhen. Die Einordnung des konkreten Anwendungsfalles setzt jedoch ein sicheres Verständnis stochastischer Prozesse voraus. Am Beispiel der Rekonstruktion der stochastischen Dynamik über eine eindimensionale Fokker-Planck-Gleichung zeigen sich deutlich die praktischen Grenzen solcher Ansätze.

Modellierung von Porenraumgeometrien und Transport in korngestützten porösen Medien

Poröse Medien spielen in der Hydrosphäre eine wesentliche Rolle bei der Strömung und beim Transport von Stoffen. In diesem Raum finden komplexe Prozesse statt: Advektion, Kon-vektion, Diffusion, hydromechanische Dispersion, Sorption, Komplexierung, Ionenaustausch und Abbau. Die strömungsmechanischen- und die Transportverhältnisse in porösen Medien werden direkt durch die Geometrie des Porenraumes selbst und durch die Eigenschaften der transportierten (oder strömenden) Medien bestimmt. In der Praxis wird eine Vielzahl von empirischen Modellen verwendet, die die Eigenschaften des porösen Mediums in repräsentativen Elementarvolumen wiedergeben. Die Ermittlung der in empirischen Modellen verwendeten Materialparameter erfolgt über Labor- oder Feldbestimmungsmethoden. Im Rahmen dieser Arbeit wurde das Computer-modell PoreFlow entwickelt, welches die hydraulischen Eigenschaften eines korngestützten porösen Mediums aus der mikroskopischen Modellierung des Fluidflusses und Transportes ableitet. Das poröse Modellmedium wird durch ein dreidimensionales Kugelpackungsmodell, zusam-mengesetzt aus einer beliebigen Kornverteilung, dargestellt. Im Modellporenraum wird die Strömung eines Fluids basierend auf einer stationären Lösung der Navier-Stokes-Gleichung simuliert. Die Ergebnisse der Modellsimulationen an verschiedenen Modellmedien werden mit den Ergebnissen von Säulenversuchen verglichen. Es zeigt sich eine deutliche Abhängigkeit der Strömungs- und Transportparameter von der Porenraumgeometrie sowohl in den Modell-simulationen als auch in den Säulenexperimenten.

Beiträge zur Modellierung von Parallelrechnern

Diese Arbeit versteht sich als Beitrag zur Modellierung von Parallelrechnern. Ein solcher Parallelrechner kann als makroskopisches physikalisches dynamisches System mit einer sehr großen Anzahl von Freiheitsgraden, diskretem Zustandsraum und diskreter Zeit aufgefasst werden. Derartige Systeme werden von der Nichtlinearen Dynamik behandelt. Jede modellmäßige Behandlung eines Systems mit derart differenzierten Wechselwirkungen muss sich auf bestimmte, dem Ziel und Zweck der Untersuchung angepasste Aspekte beschränken. Dabei müssen sowohl allgemeine Vorstellungen als auch konkretes Wissen in ein mathematisch behandelbares Modell umgesetzt werden. Die in dieser Arbeit vorgestellten Beiträge zur Modellierung von Parallelrechnern dienen mehreren Zielen. Zum einen wird ein Modell kritisch untersucht und weiterentwickelt, das dazu dienen soll, die Ausführungszeit eines konkreten parallelen Programmes auf einem konkreten Parallelrechner brauchbar vorherzusagen. Zum anderen soll die Untersuchung eines konkreten Problems aus dem Bereich von Computerwissenschaft und -technik dazu genutzt werden, ein tieferes Verständnis für das zu modellierende System zu entwickeln und daraus neue Aspekte für die Modellierung dynamischer Systeme im Allgemeinen zu gewinnen. In dieser Arbeit wird gezeigt, dass es bei der Modellierung von Parallelrechnern notwendig ist, viele technische Konstruktionseigenschaften in das Modell zu integrieren. Diese aber folgen der sehr raschen Entwicklung der Technik auf diesem Gebiet. Damit Formulierung, Test und Validierung des Modells mit der Entwicklung des Objektbereiches Schritt halten können, müssen in Zukunft neue Modellierungsverfahren entwickelt und angewendet werden, die bei Änderungen im Objektbereich eine rasche Anpassung ermöglichen. Diese Untersuchung entspricht einem interdisziplinären Ansatz, in dem einerseits Inhalte der Computerwissenschaften und andererseits Grundmethoden der experimentellen Physik verwendet werden. Dazu wurden die Vorhersagen der abstrakten Modelle mit den experimentell gewonnen Messergebnissen von realen Systemen verglichen. Auf dieser Basis wird gezeigt, dass der hierarchische Aufbau des Speichers Einflüsse von mehreren Größenordnungen auf die Ausführungsgeschwindigkeit einer Anwendung ausüben kann. Das im Rahmen der vorliegenden Arbeit entwickelte Modell der einzelnen Rechenknoten eines Parallelrechners gibt diese Effekte innerhalb eines relativen Vorhersagefehlers von nur wenigen Prozent korrekt wieder.

Prozessorientierte Modellierung des Wassertransports zur Bewertung von Hochwasserschutzmaßnahmen in bewaldeten Entstehungsgebieten

Im der vorliegenden Arbeit wurden die Einflüsse verschiedener wasserwirtschaftlicher und forstwirtschaftlicher Maßnahmen auf das Abflussverhalten von bewaldeten Gebieten untersucht. Der Schwerpunkt lag auf der Erfassung der Veränderungen im natürlichen System durch die jeweilige Maßnahme. Der Wassertransport wurde dabei auf der Prozessebene im Hangmaßstab betrachtet. Wichtige Erkenntnisse hinsichtlich der Hochwasserentstehung in bewaldeten Gebieten wurden gewonnen. Daraus abgeleitet wurde eine Bewertung sowohl der natürlich gegebenen Systemgrößen als auch der einzelnen Maßnahmen bezüglich einer Hochwasservorsorge in Wäldern.Berücksichtigt wurden Änderungen in der Bestockung und in der Feinerschließung der bewaldeten Flächen durch Wege und Wirtschaftstrecken. Die Wassertransportprozesse wurden im Gelände untersucht und in Abhängigkeit der globalen klimatischen Randbedingungen sowie den bodenhydraulischen und landnutzungsspezifischen Parameter mit dem Niederschlag-Abfluss-Modell CATFLOW simuliert.Die Abflussdynamik wird im untersuchten Maßstab im wesentlichen durch die gegebene Bodenhydraulik, die Vorfeuchte des Bodens bei einem Niederschlag und durch die Niederschlagsintensität bestimmt. Die vergleichende Untersuchung von zwei rheinland-pfälzischen Waldstandorten ergab, dass im Bereich Soonwald ein relativ hohes Hochwasserrisiko besteht. Hier wirken sich Eingriffe durch forstwirtschaftliche Bewirtschaftungsmaßnahmen stark aus. Im Bereich des Pfälzerwaldes ist nur ein sehr geringes Hochwasserrisiko gegeben, da hier die natürlichen Gebietseigenschaften die Eingriffe durch Bewirtschaftungsmaßnahmen überprägen.

Modellierung der Oberflächenreaktionen von Vesikeln und Proteinen mit Hilfe dynamischer Monte-Carlo-Simulationen

Die vorliegende Arbeit behandelt die Adsorption von Phospholipidvesikeln, Proteinen und Latex Beads an funktionalisierten festkörpergestützten Oberflächen. Der Schwerpunkt der Abhandlung liegt in der Modellierung der Kinetik mit Hilfe eines Simulationsprogramms unter Berücksichtigung des Massentransports in Staupunktgeometrie, der Adsorption und nachfolgender konformativer Änderungen adhärierter Partikel. Das Simulationsprogramm basiert auf einem RSA (random sequential adsorption) Algorithmus und erlaubt es aufgrund zusätzlich implementierter Optionen, z.B. zur Behandlung elektrostatischer Abstoßung, von Spreitprozessen oder von Desorptionsereignissen, den Adsorptionsprozess unter realistischen physikalischen Bedingungen nachzuempfinden und spezifische biologische Fragestellungen zu beantworten. Aus Anpassungen von Simulationen an experimentelle Daten ergeben sich dynamische Parameter, wie z.B. die Transport-, Adsorptions-, Spreit- und die Desorptionsrate. Die experimentellen Daten wurden mit Hilfe der Quarzmikrowaagetechnik (QCM), der Impedanzspektroskopie (IS) und der Rasterkraftmikroskopie (AFM) erhoben. Zusätzlich zur Kinetik gibt die graphische Ausgabe des Simulationsprogramms Aufschluss über die Oberflächenverteilung der adsorbierten und gespreiteten Partikel. Bei der Untersuchung von Systemen, in denen die Adsorption reversibel ist und infolge dessen Desorptionsprozesse eine wichtige Rolle spielen, wurde mit Hilfe des Simulationsprogramms ein völlig neuartiges Sensorkonzept in Aussicht gestellt. Es wurde gezeigt, dass aus der Analyse des Leistungsspektrums der Fluktuationen in der Belegung die Kinetik von Adsorption und Desorption ermittelt werden kann. Dies ist aus sensorischer Sicht interessant, da eine Vielzahl störender Einflüsse, wie z.B: der Massentransport und die elektronische Drift in diesem Fall nicht berücksichtigt werden müssen.

Methoden zur Analyse, Interpretation und Modellierung des elektrischen Widerstandes von archäologischen Steinstrukturen

In der Archäologie werden elektrische Widerstandsmessungen routinemäßig zur Prospektion von Fundstellen eingesetzt. Die Methode ist kostengünstig, leicht anwendbar und liefert in den meisten Fällen zuverlässige und leicht zu interpretierende Ergebnisse. Dennoch kann die Methode die archäologischen Strukturen in manchen Fällen nur teilweise oder gar nicht abbilden, wenn die bodenphysikalischen und bodenchemischen Eigenschaften des Bodens und der archäologischen Strukturen dies nicht zulassen. Der spezifische elektrische Widerstand wird durch Parameter wie Wassergehalt, Bodenstruktur, Bodenskelett, Bodentextur, Salinität und Bodentemperatur beeinflusst. Manche dieser Parameter, wie z.B. der Wassergehalt und die Bodentemperatur, unterliegen einer saisonalen Veränderung. Die vorliegende Arbeit untersucht den spezifischen elektrischen Widerstand von archäologischen Steinstrukturen und evaluiert die Möglichkeit, auf Grundlage von Geländemessungen und Laboranalysen archäologische Strukturen und Böden als numerische Modelle darzustellen. Dazu wurde eine Kombination von verschiedenen bodenkundlichen, geoarchäologischen und geophysikalischen Methoden verwendet. Um archäologische Strukturen und Bodenprofile als numerische Widerstandsmodelle darstellen zu können, werden Informationen zur Geometrie der Strukturen und ihren elektrischen Widerstandswerten benötigt. Dabei ist die Qualität der Hintergrundinformationen entscheidend für die Genauigkeit des Widerstandsmodells. Die Geometrie der Widerstandsmodelle basiert auf den Ergebnissen von Rammkernsondierungen und archäologische Ausgrabungen. Die an der Ausbildung des elektrischen Widerstands beteiligten Parameter wurden durch die Analyse von Bodenproben gemessen und ermöglichen durch Pedotransfer-Funktion, wie die Rhoades-Formel, die Abschätzung des spezifischen elektrischen Widerstandes des Feinbodens. Um den Einfluss des Bodenskeletts auf den spezifischen elektrischen Widerstand von Bodenprofilen und archäologischen Strukturen zu berechnen, kamen die Perkolationstheorie und die Effective Medium Theory zum Einsatz. Die Genauigkeit und eventuelle Limitierungen der Methoden wurden im Labor durch experimentelle Widerstandsmessungen an ungestörten Bodenproben und synthetischen Materialien überprüft. Die saisonale Veränderung des Wassergehalts im Boden wurde durch numerische Modelle mit der Software HYDRUS simuliert. Die hydraulischen Modelle wurden auf Grundlage der ermittelten bodenkundlichen und archäologischen Stratigraphie erstellt und verwenden die Daten von lokalen Wetterstationen als Eingangsparameter. Durch die Kombination der HYDRUS-Ergebnisse mit den Pedotransfer-Funktionen konnte der Einfluss dieser saisonalen Veränderung auf die Prospektionsergebnisse von elektrischen Widerstandsmethoden berechnet werden. Die Ergebnisse der Modellierungsprozesse wurden mit den Geländemessungen verglichen. Die beste Übereinstimmung zwischen Modellergebnissen und den Prospektionsergebnissen konnte für die Fallstudie bei Katzenbach festgestellt werden. Bei dieser wurden die Modelle auf Grundlage von archäologischen Grabungsergebnissen und detaillierten bodenkundlichen Analysen erstellt. Weitere Fallstudien zeigen, dass elektrische Widerstandsmodelle eingesetzt werden können, um den Einfluss von ungünstigen Prospektionsbedingungen auf die Ergebnisse der elektrischen Widerstandsmessungen abzuschätzen. Diese Informationen unterstützen die Planung und Anwendung der Methoden im Gelände und ermöglichen eine effektivere Interpretation der Prospektionsergebnisse. Die präsentierten Modellierungsansätze benötigen eine weitere Verifizierung durch den Vergleich der Modellierungsergebnisse mit detailliertem geophysikalischem Gelände-Monitoring von archäologischen Fundstellen. Zusätzlich könnten elektrische Widerstandsmessungen an künstlichen Mauerstrukturen unter kontrollierten Bedingungen zur Überprüfung der Modellierungsprozesse genutzt werden.

Hemiklonale Vererbung im hybridogenetischen Rana grafi-Komplex (Anura: Ranidae)

Hemiklonale Vererbung im hybridogenetische Rana grafi-Komplex (Anura: Ranidae)Im Rahmen der vorliegenden Studie wurde eine großflächige Untersuchung an südfranzösischen Wasserfröschen durchgeführt. Es wurden 31 Populationen mit 918 Tieren beprobt und mit Referenzproben genetisch verglichen. Die Genotypen der Tiere wurden mittels Allozymelektrophorese an sieben diagnostischen Loci bestimmt. Für Teilproben wurde das Alter, die Fläche der Erythrocyten und der DNA-Gehalt bestimmt. Die wichtigsten Ergebnisse der unterschiedlichen methodischen Ansätze lassen sich wie folgt zusammenfassen:(1) Es wurden neben einem großen Anteil von R. grafi (452 Tiere), R. perezi (200 Tiere) und R. ridibunda (254 Tiere) auch ein R. esculenta und zwei R. lessonae in der Camargue nachgewiesen.(2) Die Geschlechterverhältnisse der einzelnen Taxa waren stark gestört. Dies gilt besonders für den Hybriden R. grafi mit einem Männchenanteil von nur 7,3%. Das Geschlechterverhältnis von R. ridibunda ist besonders in PGR-Population zu Gunsten von Weibchen verschoben. Eine Modellierung von Populationen zeigt, daß wahrscheinlich R. grafi-Männchen durch eine gestörte Gametogenese für die gestörten Geschlechtsverhältnisse seiner Parentalformen verantwortlich ist.(3) Die genetische Analyse zeigt eine unerwartet hohe genetische Variabilität von R. ridibunda im Vergleich zu R. perezi und eine geographische Varianz der genetischen Variabilität in allen untersuchten Taxa. Die genetische Variabilität zeigte bei allen Wasserfroschformen in der Camargue ein Maximum.(4) Die hohe genetische Variabilität geht hauptsächlich, neben Einführung von verschiedenen genetischen R. ridibunda-Linien, auf Introgression zurück.

Entwicklung und Evaluierung einer Modellhierarchie zur Simulation der aktinischen Strahlung in aerosolbelasteter und bewölkter Atmosphäre

Das Verstaendnis der in der Atmosphaere stattfindendenOzonproduktions- und Ozonabbaumechanismen ist eines derwichtigsten Ziele der aktuellen Umweltforschung. Da Ozon inder Atmosphaere ausschliesslich durch photochemischeProzesse gebildet wird, ist die Kenntnis der aktinischenStrahlung eine grundlegende Voraussetzung fuer die Prognosevon Ozonkonzentrationen. Die vorliegende Arbeit befasst sichdaher mit der Modellierung der aktinischen Strahlung alsBasis zur Untersuchung photochemischer Prozesse. Dazu wurde eine Modellhierarchie entwickelt, welcheausgehend von der Berechnung der optischen Eigenschaften deratmosphaerischen Extingenten den Strahlungstransport in derAtmosphaere bestimmt. Dazu wurden verschiedeneStrahlungstransportmodelle verwendet, welche sich durch dieLoesungsmethode der Strahlungstransportgleichung und dieAnzahl der raeumlichen Dimensionen unterscheiden. Zur Ueberpruefung der entwickelten Modellhierarchie und zurUntersuchung des Einflusses der verschiedenen Extingentenauf das atmosphaerische Strahlungsfeld wurden zahlreicheSensitivitaetsstudien durchgefuehrt. Dabei zeigte sich, dasssich der Einfluss von Ozon fast ausschliesslich auf denAnteil der Strahlung mit Wellenlaengen kleiner 340 nmauswirkt. Bei der Untersuchung des Einflusses vonAerosolpartikeln und Wolkentropfen auf die aktinischeStrahlung wurde festgestellt, dass die chemischeZusammensetzung, das Vertikalprofil, die Groesse und dieGroessenverteilung der Partikel einen entscheidendenEinfluss nehmen. Eine Evaluierung des 1-dimensionalenStrahlungstransportmodells an Hand von Vergleichen mit denErgebnissen anderen Strahlungstransportmodelle undverschiedenen Messungen ergab, dass die entwickelteModellhierarchie sehr zuverlaessige Resultate liefert. Die Effekte durchbrochener Bewoelkung konnten mit einer2-dimensionalen Simulation untersucht werden. DieUntersuchungen zeigen, dass die haeufig durchgefuehrteapproximative Betrachtung von partieller Bewoelkung mit1-dimensionalen Strahlungstransportmodellen im Vergleich zuden Ergebnissen der 2-dimensionalen Simulation grosse Fehlerverursacht, da die Effekte an den seitlichen Wolkenraendernvon den 1-dimensionalen Strahlungstransportmodellen nichtberuecksichtigt werden koennen.

Ternäre polymerhaltige Lösungen - Phasenverhalten und Grenzflächenspannung

Im Rahmen dieser Arbeit wurden experimentelle und theoretische Untersuchungen zum Phasen- und Grenzflächenverhalten von ternären Systemen des Typs Lösungsmittel/Fällungsmittel/Polymer durchgeführt. Diese Art der Mischungen ist vor allem für die Planung und Durchführung der Membranherstellung von Bedeutung, bei der die genaue Kenntnis des Phasendiagramms und der Grenzflächenspannung unabdingbar ist. Als Polymere dienten Polystyrol sowie Polydimethylsiloxan. Im Fall des Polystyrols kam Butanon-2 als Lösungsmittel zum Einsatz, wobei drei niedrigmolekulare lineare Alkohole als Fällungsmittel verwendet wurden. Für Polydimethylsiloxan eignen sich Toluol als Lösungsmittel und Ethanol als Fällungsmittel. Durch Lichtstreumessungen, Dampfdruckbestimmungen mittels Headspace-Gaschromatographie (VLE-Gleichgewichte) sowie Quellungsgleichgewichten lassen sich die thermodynamischen Eigenschaften der binären Subsysteme charakterisieren. Auf Grundlage der Flory-Huggins-Theorie kann das experimentell bestimmte Phasenverhalten (LLE-Gleichgewichte) in guter Übereinstimmung nach der Methode der Direktminimierung der Gibbs'schen Energie modelliert werden. Zieht man die Ergebnisse der Aktivitätsbestimmung von Dreikomponenten-Mischungen mit in Betracht, so ergeben sich systematische Abweichungen zwischen Experiment und Theorie. Sie können auf die Notwendigkeit ternärer Wechselwirkungsparameter zurückgeführt werden, die ebenfalls durch Modellierung zugänglich sind.Durch die aus den VLE- und LLE-Untersuchungen gewonnenen Ergebnissen kann die sog. Hump-Energie berechnet werden, die ein Maß für die Entmischungstendenz darstellt. Diese Größe eignet sich gut zur Beschreibung von Grenzflächenphänomenen mittels Skalengesetzen. Die für binäre Systeme gefundenen theoretisch fundierten Skalenparameter gelten jedoch nur teilweise. Ein neues Skalengesetz lässt erstmals eine Beschreibung über die gesamte Mischungslücke zu, wobei ein Parameter durch eine gemessene Grenzflächenspannung (zwischen Fällungsmittel/Polymer) ersetzt werden kann.

Constraints on mantle melting from major and trace element systematics in residual abyssal peridotites

ZusammenfassungDie Bildung von mittelozeanischen Rückenbasalten (MORB) ist einer der wichtigsten Stoffflüsse der Erde. Jährlich wird entlang der 75.000 km langen mittelozeanischen Rücken mehr als 20 km3 neue magmatische Kruste gebildet, das sind etwa 90 Prozent der globalen Magmenproduktion. Obwohl ozeanische Rücken und MORB zu den am meisten untersuchten geologischen Themenbereichen gehören, existieren weiterhin einige Streit-fragen. Zu den wichtigsten zählt die Rolle von geodynamischen Rahmenbedingungen, wie etwa Divergenzrate oder die Nähe zu Hotspots oder Transformstörungen, sowie der absolute Aufschmelzgrad, oder die Tiefe, in der die Aufschmelzung unter den Rücken beginnt. Diese Dissertation widmet sich diesen Themen auf der Basis von Haupt- und Spurenelementzusammensetzungen in Mineralen ozeanischer Mantelgesteine.Geochemische Charakteristika von MORB deuten darauf hin, dass der ozeanische Mantel im Stabilitätsfeld von Granatperidotit zu schmelzen beginnt. Neuere Experimente zeigen jedoch, dass die schweren Seltenerdelemente (SEE) kompatibel im Klinopyroxen (Cpx) sind. Aufgrund dieser granatähnlichen Eigenschaft von Cpx wird Granat nicht mehr zur Erklärung der MORB Daten benötigt, wodurch sich der Beginn der Aufschmelzung zu geringeren Drucken verschiebt. Aus diesem Grund ist es wichtig zu überprüfen, ob diese Hypothese mit Daten von abyssalen Peridotiten in Einklang zu bringen ist. Diese am Ozeanboden aufgeschlossenen Mantelfragmente stellen die Residuen des Aufschmelz-prozesses dar, und ihr Mineralchemismus enthält Information über die Bildungs-bedingungen der Magmen. Haupt- und Spurenelementzusammensetzungen von Peridotit-proben des Zentralindischen Rückens (CIR) wurden mit Mikrosonde und Ionensonde bestimmt, und mit veröffentlichten Daten verglichen. Cpx der CIR Peridotite weisen niedrige Verhältnisse von mittleren zu schweren SEE und hohe absolute Konzentrationen der schweren SEE auf. Aufschmelzmodelle eines Spinellperidotits unter Anwendung von üblichen, inkompatiblen Verteilungskoeffizienten (Kd's) können die gemessenen Fraktionierungen von mittleren zu schweren SEE nicht reproduzieren. Die Anwendung der neuen Kd's, die kompatibles Verhalten der schweren SEE im Cpx vorhersagen, ergibt zwar bessere Resultate, kann jedoch nicht die am stärksten fraktionierten Proben erklären. Darüber hinaus werden sehr hohe Aufschmelzgrade benötigt, was nicht mit Hauptelementdaten in Einklang zu bringen ist. Niedrige (~3-5%) Aufschmelzgrade im Stabilitätsfeld von Granatperidotit, gefolgt von weiterer Aufschmelzung von Spinellperidotit kann jedoch die Beobachtungen weitgehend erklären. Aus diesem Grund muss Granat weiterhin als wichtige Phase bei der Genese von MORB betrachtet werden (Kapitel 1).Eine weitere Hürde zum quantitativen Verständnis von Aufschmelzprozessen unter mittelozeanischen Rücken ist die fehlende Korrelation zwischen Haupt- und Spuren-elementen in residuellen abyssalen Peridotiten. Das Cr/(Cr+Al) Verhältnis (Cr#) in Spinell wird im Allgemeinen als guter qualitativer Indikator für den Aufschmelzgrad betrachtet. Die Mineralchemie der CIR Peridotite und publizierte Daten von anderen abyssalen Peridotiten zeigen, dass die schweren SEE sehr gut (r2 ~ 0.9) mit Cr# der koexistierenden Spinelle korreliert. Die Auswertung dieser Korrelation ergibt einen quantitativen Aufschmelz-indikator für Residuen, welcher auf dem Spinellchemismus basiert. Damit kann der Schmelzgrad als Funktion von Cr# in Spinell ausgedrückt werden: F = 0.10×ln(Cr#) + 0.24 (Hellebrand et al., Nature, in review; Kapitel 2). Die Anwendung dieses Indikators auf Mantelproben, für die keine Ionensondendaten verfügbar sind, ermöglicht es, geochemische und geophysikalischen Daten zu verbinden. Aus geodynamischer Perspektive ist der Gakkel Rücken im Arktischen Ozean von großer Bedeutung für das Verständnis von Aufschmelzprozessen, da er weltweit die niedrigste Divergenzrate aufweist und große Transformstörungen fehlen. Publizierte Basaltdaten deuten auf einen extrem niedrigen Aufschmelzgrad hin, was mit globalen Korrelationen im Einklang steht. Stark alterierte Mantelperidotite einer Lokalität entlang des kaum beprobten Gakkel Rückens wurden deshalb auf Primärminerale untersucht. Nur in einer Probe sind oxidierte Spinellpseudomorphosen mit Spuren primärer Spinelle erhalten geblieben. Ihre Cr# ist signifikant höher als die einiger Peridotite von schneller divergierenden Rücken und ihr Schmelzgrad ist damit höher als aufgrund der Basaltzusammensetzungen vermutet. Der unter Anwendung des oben erwähnten Indikators ermittelte Schmelzgrad ermöglicht die Berechnung der Krustenmächtigkeit am Gakkel Rücken. Diese ist wesentlich größer als die aus Schweredaten ermittelte Mächtigkeit, oder die aus der globalen Korrelation zwischen Divergenzrate und mittels Seismik erhaltene Krustendicke. Dieses unerwartete Ergebnis kann möglicherweise auf kompositionelle Heterogenitäten bei niedrigen Schmelzgraden, oder auf eine insgesamt größere Verarmung des Mantels unter dem Gakkel Rücken zurückgeführt werden (Hellebrand et al., Chem.Geol., in review; Kapitel 3).Zusätzliche Informationen zur Modellierung und Analytik sind im Anhang A-C aufgeführt

Neue Methoden zur Computersimulation von Polymersystemen auf verschiedenen Längenskalen und ihre Anwendungen

In dieser Dissertation stellen wir einen neuen Ansatz zurModellierungvon Polymersystemen vor. Es werden (von methodischer Seiteher) zweiautomatisierte Iterationschemata dazu eingeführt,Kraftfeldparametermesoskopischer Polymersysteme systematisch zu optimieren:DasSimplex-Verfahren und das Struktur-Differenzen-Verfahren. Sowerdendiejenigen Freiheitsgrade aus Polymersystemen eliminiert,die einehohe Auflösung erfordern, was die Modellierung größerersystemeermöglicht. Nach Tests an einfachen Flüssigkeiten werdenvergröberteModelle von drei prototypischen Polymeren (Polyacrylsäure,Polyvinylalkohol und Polyisopren) in unterschiedlichenUmgebungen(gutes Lösungsmittel und Schmelze) entwickelt und ihrVerhalten aufder Mesoskala ausgiebig geprüft. Die zugehörige Abbildung(vonphysikalischer Seite her) so zu gestalten, daß sie dieunverwechselbaren Charakteristiken jedes systems auf diemesoskopischeLängenskala überträgt, stellt eine entscheidende Anforderungan dieautomatisierten Verfahren dar. Unsere Studien belegen, daß mesoskopische Kraftfeldertemperatur- unddichtespezifisch sind und daher bei geändernden Bedingungennachoptimiert werden müssen. Gleichzeitig läßt sichabschätzen, beiwelchen Umgebungsbedingungen dies noch nicht notwendig wird.In allenFällen reichen effektive Paarpotentiale aus, einrealistischesmesoskopisches Modell zu konstruieren. VergröberteSimulationenwerden im Falle der Polyacrylsäure erfolgreich gegenexperimentelleLichtstreudaten getestet. Wir erzielen für Molmassen bis zu300000g/mol eine hervorragende Übereinstimmung für denhydrodynamischenRadius. Unsere Ergebnisse erklären auch Korrekturen zudessenVerhalten als Funktion der Kettenlänge ('Skalenverhalten'). Im Fallevon Polyisopren untersuchen wir sowohl statische als auchdynamischeGrößen und stellen klare Unterschiede unserer Ergebnisse zudeneneines einfachen semi-flexiblen Mesoskalenmodells fest. InderProteinforschung werden aus Datenbanken gewonnene effektivePaarwechselwirkungen dazu verwendet, die freie Energie einesneuensystems vorherzusagen. Wir belegen in einem Exkurs mittelsGittersimulationen, daß es selbst in einfachsten Fällennicht gelingt,dies auch nur qualitativ korrekt zu bewerkstelligen.

Untersuchung der Strukturdynamik arbeitender F1-ATPase und ATP-Synthase aus Micrococcus luteus in Einzelschußexperimenten mit Synchrotronstrahlung

ZusammenfassungDie ATP-Synthase koppelt im Energiestoffwechsel der Zellen den Protonentransport über die biologische Membran mit der Synthese des energiespeichernden Moleküls ATP aus ADP und Phosphat. ATP-Synthasen bestehen aus 2 Subkomplexen, wobei der katalytische F1-Teil von der membranständigen Domäne abgelöst werden kann und nur zur ATP-Hydrolyse fähig ist. Der hochkooperative Reaktionsmechanismus der dreizentrigen ATP-Synthasen ist weitgehend unklar.Im Rahmen dieser Arbeit wurde der ATP-Synthasekomplex und ihr wasserlösliches katalytisches F1-Fragment aus Micrococcus luteus in präparativem Maßstab mittels chromatographischer Trennmethoden isoliert. Die Überprüfung der Funktionalität beider Enzyme erfolgte mit enzymatischen Methoden. Durch zeitaufgelöste Röntgenkleinwinkelstreuung wurde die Strukturdynamik der arbeitenden ATP-Synthase und ihres F1-Fragmentes aus Micrococcus luteus im Laufe des ATP-Hydrolysezyklus untersucht. Diese Methode diente zum Nachweis weiträumiger Konformationsänderungen innerhalb der arbeitenden Enzyme unter nativen physiologischen Bedingungen. Die zeitaufgelösten Streuexperimente fanden an der ESRF (Europäische Synchrotronstrahlungsquelle) in Grenoble (F) statt. Dort wurden für beide Enzyme im Laufe des ATP-Hydrolysezykus molekulare Bewegungen nachgewiesen. Als Referenz zu den zeitaufgelösten Messungen dienten statische Messungen zur Strukturuntersuchung der Proteine am schwächeren DESY. Anhand dieser Strukturdaten wurden Molekülmodelle der F1-ATPase und ATP-Synthase aus Micrococcus luteus konstruiert. Das Molekülmodell der F1-ATPase war die Grundlage zur Modellierung einzelner Teilschritte des ATP-Hydrolysezyklus bei 20°C. Die experimentellen Daten wurden mit einer Kippbewegung der membranseitigen Domäne der katalytischen b-Untereinheiten der F1-ATPase während des ATP-Hydrolysezyklus interpretiert.

Tropospheric photochemistry of ozone, its precursors and the hydroxyl radical

Zusammenfassung Ein 3-dimensionales globales Modell der unterenAtmosphäre wurde für die Untersuchung derOzonchemie, sowie der Chemie des Hydroxylradikals (OH) undwichtiger Vorläufersubstanzen, wie reaktiverStickstoffverbindungen und Kohlenwasserstoffe, verwendet.Hierfür wurde die Behandlung vonNicht-Methan-Kohlenwasserstoffen (NMKW) hinzugefügt,was auch die Entwicklung einer vereinfachten Beschreibungihrer Chemie, sowie die Erfassung von Depositionsprozessenund Emissionen erforderte. Zur Lösung der steifengewöhnlichen Differentialgleichungen der Chemie wurdeeine schnelles Rosenbrock-Verfahren eingesetzt, das soimplementiert wurde, dass die Modell-Chemie fürzukünftige Studien leicht abgeändert werden kann. Zur Evaluierung des Modells wurde ein umfangreicherVergleich der Modellergebnisse mit Bodenmessungen, sowieFlugzeug-, Sonden- und Satelliten-Daten durchgeführt.Das Modell kann viele Aspekte der Beobachtungen derverschieden Substanzen realistisch wiedergeben. Es wurdenjedoch auch einige Diskrepanzen festgestellt, die Hinweiseauf fehlerhafte Emissionsfelder oder auf die Modell-Dynamikund auch auf fehlende Modell-Chemie liefern. Zur weiteren Untersuchung des Einflusses verschiedenerStoffgruppen wurden drei Läufe mit unterschiedlichkomplexer Chemie analysiert. Durch das Berücksichtigender NMKW wird die Verteilung mehrerer wichtiger Substanzensignifikant beeinflusst, darunter z.B. ein Anstieg desglobalen Ozons. Es wurde gezeigt, dass die biogene SubstanzIsopren etwa die Hälfte des Gesamteffekts der NMKWausmachte (mehr in den Tropen, weniger anderswo). In einer Sensitivitätsstudie wurden die Unsicherheitenbei der Modellierung von Isopren weitergehend untersucht.Dabei konnte gezeigt werden, dass die Unsicherheit beiphysikalischen Aspekten (Deposition und heterogene Prozesse)ebenso groß sein kann, wie die aus dem chemischenGasphasen-Mechanismus stammende, welche zu globalbedeutsamen Abweichungen führte. Lokal können sichnoch größere Abweichungen ergeben. Zusammenfassend kann gesagt werden, dass die numerischenStudien dieser Arbeit neue Einblicke in wichtige Aspekte derPhotochemieder Troposphäre ergaben und in Vorschläge fürweiter Studien mündeten, die die wichtigsten gefundenenUnsicherheiten weiter verringern könnten.

Monte Carlo-Simulationen zum Adsorptionsverhalten von Homo- und Copolymerlösungen

Die vorliegende Arbeit beschäftigt sich mit demAdsorptionsverhalten von Homo- und Blockopolymeren anglatten Wänden im guten Lösungsmittel. Zur Modellierung wirddabei ein vergröbertes Kugel-Feder-Modell benutzt, das diePolymere als Kette von effektiven Monomeren auffasst.Benachbarte Monomere einer Kette werden durch einBindungspotential aneinander gekoppelt; zusätzlich wird eineVolumenausschlusswechselwirkung zwischen allen Monomeren derLösung berücksichtigt. Das Lösungsmittel wird dabei implizitdurch diese effektiven Potentiale berücksichtigt, wobei dieLösungsmittelqualität über die Temperatur geregelt wird.Während das System in z-Richtung durch zwei ebene Wändebegrenzt wird, von denen eine attraktiv mit den Monomerenwechselwirkt, werden in x- und y-Richtung periodischeRandbedingungen vorgegeben. Die Bewegung der Polymereerfolgt durch Monte Carlo-Simulationen im räumlichenKontinuum, welche auf (großkanonischen) Configurational BiasMonte Carlo-Algorithmen sowie lokalen Monomerverschiebungenund Reptationsbewegungen basieren.Über verschiedene Auswertemethoden werden anhand vonEinzelkettensimulationen die (modellabhängige)Adsorptionsschwelle - welche durch Simualtionen derAdsorption aus der Lösung bestätigt wird - sowie der(universelle) Crossoverexponent - dessen gefundener Wertsehr gut mit den neusten analytischen Berechnungenübereinstimmt - bestimmt. Bei der Simulation derBlockcopolymere bilden sich bei niedrigerOberflächenbelegungsdichte pilzartige Konformationen('mushroom'-Bereich), bei hoher OberflächenbelegungsdichtePolymerbürsten ('brush'-Bereich). Die theoretischenVorhersagen für den jeweiligen Bereich werden von denSimulationen im Wesentlichen gut wiedergegeben.

Informationsstrukturen der protokollbasierten Medizin

Die protokollbasierte Medizin stellt einen interdisziplinären Brennpunkt der Informatik dar. Als besonderer Ausschnitt der medizinischen Teilgebiete erlaubt sie die relativ formale Spezifikation von Prozessen in den drei Bereichen der Prävention, Diagnose und Therapie.Letzterer wurde immer besonders fokussiert und gilt seit jeher im Rahmen klinischer Studien als Projektionsfläche für informationstechnologische Konzepte. Die Euphorie der frühen Jahre ernüchtert sich jedoch bei jeder Bilanz. Nur sehr wenige der unzählbaren Projekte haben ihre Routine in der alltäglichen Praxis gefunden. Die meisten Vorhaben sind an der Illusion der vollständigen Berechenbarkeit medizinischer Arbeitsabläufe gescheitert. Die traditionelle Sichtweise der klinischen Praxis beruht auf einer blockorientierten Vorstellung des Therapieausführungsprozesses. Sie entsteht durch seine Zerlegung in einzelne Therapiezweige, welche aus vordefinierten Blöcken zusammengesetzt sind. Diese können sequentiell oder parallel ausgeführt werden und sind selbst zusammengesetzt aus jeweils einer Menge von Elementen,welche die Aktivitäten der untersten Ebene darstellen. Das blockorientierte Aufbaumodell wird ergänzt durch ein regelorientiertes Ablaufmodell. Ein komplexes Regelwerk bestimmt Bedingungen für die zeitlichen und logischen Abhängigkeiten der Blöcke, deren Anordnung durch den Ausführungsprozeß gebildet wird. Die Modellierung der Therapieausführung steht zunächst vor der grundsätzlichen Frage, inwieweit die traditionelle Sichtweise für eine interne Repräsentation geeignet ist. Das übergeordnete Ziel besteht in der Integration der unterschiedlichen Ebenen der Therapiespezifikation. Dazu gehört nicht nur die strukturelle Komponente, sondern vorallem die Ablaufkomponente. Ein geeignetes Regelmodell ist erforderlich, welches den spezifischen Bedürfnissen der Therapieüberwachung gerecht wird. Die zentrale Aufgabe besteht darin, diese unterschiedlichen Ebenen zusammenzuführen. Eine sinnvolle Alternative zur traditionellen Sichtweise liefert das zustandsorientierte Modell des Therapieausführungsprozesses. Das zustandsorientierte Modell beruht auf der Sichtweise, daß der gesamte Therapieausführungsprozeß letztendlich eine lineare Folge von Zuständen beschreibt, wobei jeder Zustandsübergang durch ein Ereignis eingeleitet wird, an bestimmte Bedingungen geknüpft ist und bestimmte Aktionen auslösen kann. Die Parallelität des blockorientierten Modells tritt in den Hintergrund, denn die Menge der durchzuführenden Maßnahmen sind lediglich Eigenschaften der Zustände und keine strukturellen Elemente der Ablaufspezifikation. Zu jedem Zeitpunkt ist genau ein Zustand aktiv, und er repräsentiert eine von endlich vielen klinischen Situationen, mit all ihren spezifischen Aktivitäten und Ausführungsregeln. Die Vorteile des zustandsorientierten Modells liegen in der Integration. Die Grundstruktur verbindet die statische Darstellung der möglichen Phasenanordnungen mit der dynamischen Ausführung aktiver Regeln. Die ursprünglichen Inhalte des blockorientierten Modells werden als gewöhnliche Eigenschaften der Zustände reproduziert und stellen damit nur einen Spezialfall der zustandsbezogenen Sicht dar.Weitere Möglichkeiten für die Anreicherung der Zustände mit zusätzlichen Details sind denkbar wie sinnvoll. Die Grundstruktur bleibt bei jeder Erweiterung jedoch die gleiche. Es ergibt sich ein wiederverwendbares Grundgerüst,ein gemeinsamer Nenner für die Erfüllung der Überwachungsaufgabe.