Grenzflächenspannung und Morphologie fließender zweiphasiger Polymermischungen

In der vorliegenden Arbeit wurde die Morphologie von zweiphasigen Polymermischungen unter Scherung in situ mit Hilfe einer Kombination aus optischer Scherzelle, Durchlichtmikroskop und computergestützten CCD-Kamera untersucht. Als Modellblends dienten die unverträglichen, bei Raumtemperatur flüssigen Polymersysteme Polyisobutylen (PIB)/Polydimethylsiloxan (PDMS) (I) und Poly(dimethyl-co-methylphenyl)siloxan/PDMS (II). Alle Komponenten verhalten sich bei den verwendeten Scherraten newtonisch.Eine der wichtigsten Einflussgrößen für die Blendmorphologie ist die Grenzflächenspannung gamma 12. Sie wurde für I und II mit Hilfe der Methode der Tropfenrelaxation (dynamisch) als Funktion der Zeit bestimmt. Diese Methode erlaubt die Messung von gamma 12 für Tropfen der Phase A in B sowie von Tropfen B in A. Bei der Methode des hängenden Tropfens (statisch) muss der Tropfen aus der Phase mit der höheren Dichte bestehen. Wo der Vergleich der beiden Methoden möglich ist, stimmen die Ergebnisse für beide Systeme sehr gut überein. Bei II sind die aus der Tropfenrelaxation erhaltenen gamma 12-Werte der beiden komplementären Zusammensetzungen im Rahmen des Fehlers gleich, bei I zeigt ein PIB-Tropfen in PDMS einen um 40 % niedrigeren Wert als ein PDMS-Tropfen in PIB, dies wird auf die Diffusion von kurzkettigen Anteilen des PDMS in die Grenzschicht zurückgeführt. Die Grenzflächenspannung hängt also unter Umständen auch bei binären Systemen deutlich von der Zusammensetzung ab.Für II wurde die Blendmorphologie über den gesamten Zusammensetzungsbereich untersucht. Die häufig beobachteten cokontinuierlichen Strukturen treten bei keiner Zusammensetzung auf. Die Phaseninversion erkennt man in einer sprunghaften Änderung der Tropfengröße zwischen phiPDMS <= 0,400 und 0,500; zudem lässt sich die Zeitabhängigkeit der Radien durch Auftragung gegen das Produkt aus der Deformation und dem Quadrat des Volumenbruchs der Tropfenphase für 0 <= phiPDMS <= 0,400 sowie 0,500 <= phiPDMS <= 1 normieren. Für I und II wurde die Morphologieentwicklung bei 25 °C nach Vorscherung bei 100 bzw. 50 s-1 und anschließendem Sprung der Scherrate auf deutlich niedrigere Werte als Funktion der Zeit verfolgt. Hierbei erhält man bei genügend langer Messdauer (mindestens 200 000-300 000 Schereinheiten) konstante Tropfengrößen. Zum einen handelt es sich dabei um pseudo-stationäre Werte, die nur durch Koaleszenz bestimmt sind, zum anderen um echte stationäre Radien, die durch gleichzeitig ablaufende Koaleszenz und Zerteilung entstehen. Für I liegen die stationären Mittelwerte auf der Zerteilungskurve, für II hingegen auf der Koaleszenzkurve.Der Einfluss eines grenzflächenwirksamen Additivs wurde anhand von I durch Zugabe des Blockcopolymer PIB-b-PDMS zu PIB untersucht. Der Vergleich des zeitlichen Verlaufs von gamma 12 mit der Morphologieentwicklung zeigt, dass das Additiv eine Stabilisierung der feinen Tropfen/Matrix-Struktur des Blends durch Hinderung der Koaleszenz und nicht durch Reduktion der Grenzflächenspannung bewirkt.

Efficient parallel proximity queries and an application to highly complex motion planning problems with many narrow passages

In vielen Bereichen der industriellen Fertigung, wie zum Beispiel in der Automobilindustrie, wer- den digitale Versuchsmodelle (sog. digital mock-ups) eingesetzt, um die Entwicklung komplexer Maschinen m ̈oglichst gut durch Computersysteme unterstu ̈tzen zu k ̈onnen. Hierbei spielen Be- wegungsplanungsalgorithmen eine wichtige Rolle, um zu gew ̈ahrleisten, dass diese digitalen Pro- totypen auch kollisionsfrei zusammengesetzt werden k ̈onnen. In den letzten Jahrzehnten haben sich hier sampling-basierte Verfahren besonders bew ̈ahrt. Diese erzeugen eine große Anzahl von zuf ̈alligen Lagen fu ̈r das ein-/auszubauende Objekt und verwenden einen Kollisionserken- nungsmechanismus, um die einzelnen Lagen auf Gu ̈ltigkeit zu u ̈berpru ̈fen. Daher spielt die Kollisionserkennung eine wesentliche Rolle beim Design effizienter Bewegungsplanungsalgorith- men. Eine Schwierigkeit fu ̈r diese Klasse von Planern stellen sogenannte “narrow passages” dar, schmale Passagen also, die immer dort auftreten, wo die Bewegungsfreiheit der zu planenden Objekte stark eingeschr ̈ankt ist. An solchen Stellen kann es schwierig sein, eine ausreichende Anzahl von kollisionsfreien Samples zu finden. Es ist dann m ̈oglicherweise n ̈otig, ausgeklu ̈geltere Techniken einzusetzen, um eine gute Performance der Algorithmen zu erreichen.rnDie vorliegende Arbeit gliedert sich in zwei Teile: Im ersten Teil untersuchen wir parallele Kollisionserkennungsalgorithmen. Da wir auf eine Anwendung bei sampling-basierten Bewe- gungsplanern abzielen, w ̈ahlen wir hier eine Problemstellung, bei der wir stets die selben zwei Objekte, aber in einer großen Anzahl von unterschiedlichen Lagen auf Kollision testen. Wir im- plementieren und vergleichen verschiedene Verfahren, die auf Hu ̈llk ̈operhierarchien (BVHs) und hierarchische Grids als Beschleunigungsstrukturen zuru ̈ckgreifen. Alle beschriebenen Verfahren wurden auf mehreren CPU-Kernen parallelisiert. Daru ̈ber hinaus vergleichen wir verschiedene CUDA Kernels zur Durchfu ̈hrung BVH-basierter Kollisionstests auf der GPU. Neben einer un- terschiedlichen Verteilung der Arbeit auf die parallelen GPU Threads untersuchen wir hier die Auswirkung verschiedener Speicherzugriffsmuster auf die Performance der resultierenden Algo- rithmen. Weiter stellen wir eine Reihe von approximativen Kollisionstests vor, die auf den beschriebenen Verfahren basieren. Wenn eine geringere Genauigkeit der Tests tolerierbar ist, kann so eine weitere Verbesserung der Performance erzielt werden.rnIm zweiten Teil der Arbeit beschreiben wir einen von uns entworfenen parallelen, sampling- basierten Bewegungsplaner zur Behandlung hochkomplexer Probleme mit mehreren “narrow passages”. Das Verfahren arbeitet in zwei Phasen. Die grundlegende Idee ist hierbei, in der er- sten Planungsphase konzeptionell kleinere Fehler zuzulassen, um die Planungseffizienz zu erh ̈ohen und den resultierenden Pfad dann in einer zweiten Phase zu reparieren. Der hierzu in Phase I eingesetzte Planer basiert auf sogenannten Expansive Space Trees. Zus ̈atzlich haben wir den Planer mit einer Freidru ̈ckoperation ausgestattet, die es erlaubt, kleinere Kollisionen aufzul ̈osen und so die Effizienz in Bereichen mit eingeschr ̈ankter Bewegungsfreiheit zu erh ̈ohen. Optional erlaubt unsere Implementierung den Einsatz von approximativen Kollisionstests. Dies setzt die Genauigkeit der ersten Planungsphase weiter herab, fu ̈hrt aber auch zu einer weiteren Perfor- mancesteigerung. Die aus Phase I resultierenden Bewegungspfade sind dann unter Umst ̈anden nicht komplett kollisionsfrei. Um diese Pfade zu reparieren, haben wir einen neuartigen Pla- nungsalgorithmus entworfen, der lokal beschr ̈ankt auf eine kleine Umgebung um den bestehenden Pfad einen neuen, kollisionsfreien Bewegungspfad plant.rnWir haben den beschriebenen Algorithmus mit einer Klasse von neuen, schwierigen Metall- Puzzlen getestet, die zum Teil mehrere “narrow passages” aufweisen. Unseres Wissens nach ist eine Sammlung vergleichbar komplexer Benchmarks nicht ̈offentlich zug ̈anglich und wir fan- den auch keine Beschreibung von vergleichbar komplexen Benchmarks in der Motion-Planning Literatur.

Computer simulations of charged systems in partially periodic geometries

This work presents algorithms for the calculation of the electrostatic interaction in partially periodic systems. The framework for these algorithms is provided by the simulation package ESPResSo, of which the author was one of the main developers. The prominent features of the program are listed and the internal structure is described. In the following, algorithms for the calculation of the Coulomb sum in three dimensionally periodic systems are described. These methods are the foundations for the algorithms for partially periodic systems presented in this work. Starting from the MMM2D method for systems with one non-periodic coordinate, the ELC method for these systems is developed. This method consists of a correction term which allows to use methods for three dimensional periodicity also for the case of two periodic coordinates. The computation time of this correction term is neglible for large numbers of particles. The performance of MMM2D and ELC are demonstrated by results from the implementations contained in ESPResSo. It is also discussed, how different dielectric constants inside and outside of the simulation box can be realized. For systems with one periodic coordinate, the MMM1D method is derived from the MMM2D method. This method is applied to the problem of the attraction of like-charged rods in the presence of counterions, and results of the strong coupling theory for the equilibrium distance of the rods at infinite counterion-coupling are checked against results from computer simulations. The degree of agreement between the simulations at finite coupling and the theory can be characterized by a single parameter gamma_RB. In the special case of T=0, one finds under certain circumstances flat configurations, in which all charges are located in the rod-rod plane. The energetically optimal configuration and its stability are determined analytically, which depends on only one parameter gamma_z, similar to gamma_RB. These findings are in good agreement with results from computer simulations.

Computer simulations of slowly relaxing systems in external fields

Eine der offenen Fragen der aktuellen Physik ist das Verständnis von Systemen im Nichtgleichgewicht. Im Gegensatz zu der Gleichgewichtsphysik ist in diesem Bereich aktuell kein Formalismus bekannt der ein systematisches Beschreiben der unterschiedlichen Systeme ermöglicht. Um das Verständnis über diese Systeme zu vergrößern werden in dieser Arbeit zwei unterschiedliche Systeme studiert, die unter einem externen Feld ein starkes nichtlineares Verhalten zeigen. Hierbei handelt es sich zum einen um das Verhalten von Teilchen unter dem Einfluss einer extern angelegten Kraft und zum anderen um das Verhalten eines Systems in der Nähe des kritischen Punktes unter Scherung. Das Modellsystem in dem ersten Teil der Arbeit ist eine binäre Yukawa Mischung, die bei tiefen Temperaturen einen Glassübergang zeigt. Dies führt zu einer stark ansteigenden Relaxationszeit des Systems, so dass man auch bei kleinen Kräften relativ schnell ein nichtlineares Verhalten beobachtet. In Abhängigkeit der angelegten konstanten Kraft können in dieser Arbeit drei Regime, mit stark unterschiedlichem Teilchenverhalten, identifiziert werden. In dem zweiten Teil der Arbeit wird das Ising-Modell unter Scherung betrachtet. In der Nähe des kritischen Punkts kommt es in diesem Modell zu einer Beeinflussung der Fluktuationen in dem System durch das angelegte Scherfeld. Dies hat zur Folge, dass das System stark anisotrop wird und man zwei unterschiedliche Korrelationslängen vorfindet, die mit unterschiedlichen Exponenten divergieren. Infolgedessen lässt sich der normale isotrope Formalismus des "finite-size scaling" nicht mehr auf dieses System anwenden. In dieser Arbeit wird gezeigt, wie dieser auf den anisotropen Fall zu verallgemeinern ist und wie damit die kritischen Punkte, sowie die dazu gehörenden kritischen Exponenten berechnet werden können.