Stability, viscous dissipation and local thermal non-equilibrium in fluid saturated porous media

In this thesis, the field of study related to the stability analysis of fluid saturated porous media is investigated. In particular the contribution of the viscous heating to the onset of convective instability in the flow through ducts is analysed. In order to evaluate the contribution of the viscous dissipation, different geometries, different models describing the balance equations and different boundary conditions are used. Moreover, the local thermal non-equilibrium model is used to study the evolution of the temperature differences between the fluid and the solid matrix in a thermal boundary layer problem. On studying the onset of instability, different techniques for eigenvalue problems has been used. Analytical solutions, asymptotic analyses and numerical solutions by means of original and commercial codes are carried out.

Physikalische Eigenschaften kolloidaler Festkörper

In dieser Arbeit wird das Phasenverhalten fluid-kristallin und kristallin-amorph, die elastischen Eigenschaften, das Nukleationsverhalten und das diffusive Verhalten ladungsstabilisierter Kolloide aus sphärischen Polystyrol- und Polytetrafluorethylenpartikeln in wässerigen Dispersionsmitteln bei sehr geringem Fremdionengehalt systematisch untersucht. Die dazugehörigen Messungen werden an einer neuartigen selbstkonstruierten Kombinationslichtstreuapparatur durchgeführt, die die Meßmethoden der dynamischen Lichtstreuung, statischen Lichtstreuung und Torsionsresonanzspektroskopie in sich vereint. Die drei Meßmethoden sind optimal auf die Untersuchung kolloidaler Festkörper abgestimmt. Das elastische Verhalten der Festkörper kann sehr gut durch die Elastizitätstheorie atomarer Kristallsysteme beschrieben werden, wenn ein Debye-Hückel-Potential im Sinne des Poisson-Boltzmann-Cell Modells als Wechselwirkungspotential verwendet wird. Die ermittelten Phasengrenzen fluid-kristallin stehen erstmalig in guter Übereinstimmung mit Ergebnissen aus molekulardynamischen Simulationen, wenn die in der Torsionsresonanzspektroskopie bestimmte Wechselwirkungsenergie zu Grunde gelegt wird. Neben der Gleichgewichtsstruktur sind Aussagen zur Verfestigungskinetik möglich. Das gefundene Nukleationserhalten kann gut durch die klassische Nukleationstheorie beschrieben werden, wenn bei niedriger Unterkühlung der Schmelze ein Untergrund heterogener Keimung berücksichtigt wird. PTFE-Partikel zeigen auch bei hohen Konzentrationen nur geringfügige Mehrfachstreuung. Durch ihren Einsatz ist erstmals eine systematische Untersuchung des Glasübergangs in hochgeladenen ladungsstabilisierten Systemen möglich. Ladungsstabilisierte Kolloide unterscheiden sich vor allem durch ihre extreme Kristallisationstendenz von früher untersuchten Hartkugelsystemen. Bei hohen Partikelkonzentrationen (Volumenbrüche größer 10 Prozent) kann ein glasartiger Festkörper identifiziert werden, dessen physikalisches Verhalten die Existenz eines Bernalglases nahe legt. Der Glasübergang ist im Vergleich mit den in anderen kolloidalen Systemen und atomaren Systemen beobachteten Übergängen von sehr unterschiedlichem Charakter. Im verwendeten PTFE-System ist auf Grund der langreichweitigen stark repulsiven Wechselwirkung kein direkter Zugang des Glaszustandes aus der übersättigten Schmelze möglich. Der amorphe Festkörper entsteht hier aus einer nanokristallinen Phase. Die Keimrate steigt im zugänglichen Meßbereich annähernd exponentiell mit der Partikelanzahldichte, so daß man feststellen kann, daß der Glaszustand nicht durch Unterdrückung der Nukleation, sondern durch eine Forcierung derselben erreicht wird.

Molekulardynamik-Untersuchungen zur atomistischen Struktur & Dynamik von binären Mischgläsern

Mit Hilfe von Molekulardynamik-Computersimulationen werdenin dieser Arbeit die Struktur und Dynamik von Gläsern undSchmelzen der Systeme 'NSx'=(Na2O)(xSiO2), mit x=2,3,5, und'AS2'=(Al2O3)(2SiO2) untersucht. Zur Beschreibung dermikroskopischen Wechselwirkungen dient ein Modellpotenzial,das auf einem effektiven Paarpotenzial aus der Literaturaufbaut. Simuliert wurden Teilchentrajektorien über mehrereNanosekunden im Bereich 6100 K >= T >= 2100 K sowie dieGlasstruktur bei 100 K (NSx) bzw. 300 K (AS2). Das Aufbrechen der tetraedrischen Netzwerkstruktur durch denNetzwerkmodifikator Na2O führt zur Ausbildung einerzusätzlichen Struktur auf intermediären Längenskalen, diebei allen Systemen NSx etwa dem Abstand übernächster Na-bzw. Si-Nachbarn eines Na-Atoms entspricht. Die diffusiveDynamik ist in allen Systeme NSx bis zu drei Größenordnungenschneller als in SiO2 . Sie nimmt mit wachsenderNa-Konzentration zu. Die Na-Diffusion zeigtArrhenius-Verhalten; hierbei vollführen die Na-Atome einaktiviertes Hüpfen durch eine erstarrte Si-O-Matrix. DieZustandsdichten werden bis ca. 20 THz durch dominanteNa-Moden bestimmt. Bei hohen Frequenzen weichen die für SiO2 typischen intratetraedrischen Si-O-Schwingungsmodenauf.Im Gegensatz zu Na2O wird Al2O3 in die tetraedrischeNetzwerkstruktur eingebunden. AS2 zeigt eine überwiegend ausAlO4- und SiO4-Tetraedern verknüpfte Polyederstruktur, beider sich die AlO4-Tetraeder lokal anders anordnen als dieSiO4-Tetraeder, um Ladungsausgleich zu gewährleisten. Esbilden sich typische 3(Si,Al)O4-Bausteine ('3-Cluster'). Sie führen auf intermediären Längenskalen zur Ausbildungeines Al-reichen perkolierenden Netzwerks, das dieSiO4-Struktur durchdringt. Im Vergleich zu SiO2 erhöht sichdie diffusive Dynamik aller Komponenten in AS2 um ca. zweiGrößenordnungen. Die intratetraedrischenAl-O-Steckschwingungsmoden sind wesentlich weicher als die typischen intratetraedrischen Si-O-Moden des SiO2.

Sintesi di immagini realistiche mediante l'uso del Photon differentials come estensione del Photon mapping

Questa tesi intende presentare una tecnica per la sintesi di immagini realistiche al calcolatore basata sul concetto di particle tracing. Il metodo proposto opera una stima sulla densità locale dei fotoni estendendo il concetto del photon differentials anche alla gestione delle riflessioni diffusive. Si è scelto di implementare il nuovo algoritmo di illuminazione globale all’interno di XCModel come estensione del photon mapping.

Nature and origin of ultramafic lamprophyres and carbonatites from the borders of the Labrador Sea

Die Ränder des Labrador Meeres wurden während des späten Neoproterozoikums intensiv von karbonatreichen silikatischen Schmelzen durchsetzt. Diese Schmelzen bildeted sich bei Drucken zwischen ca. 4-6 GPa (ca. 120-180 km Tiefe) an der Basis der kontinentalen Mantel-Lithosphäre. Diese Magmengenerierung steht in zeitlichem und räumlichem Zusammenhang mit kontinentalen Extensionsprozessen, welche zu beiden Seiten des sich öffnenden Iapetus-Ozeans auftraten.

Eigenschaftskorrelationen in kolloidalen Festkörpern und Fluiden aus optischen Experimenten

Die Arbeit beschäftigt sich mit ein- und zweikomponentigen, geladenen Kolloidsystemen, die in vollentsalzten wässrigen und organischen Dispersionen kristalline Strukturen ausbilden. Im ersten Teil der Arbeit wird die Wechselwirkung der Kolloide mit verschiedenen Methoden charakterisiert. Dabei zeigten sich quantitative Übereinstimmungen zwischen den Resultaten aus Zellenmodellrechnungen und aus elektrokinetischen Messungen einerseits und Messungen des Phasenverhaltens und der Elastizität andererseits. Diese nunmehr gut gesicherten Diskrepanzen und Korrelationen bedürfen des theoretischen Verständnisses. Im zweiten Teil der Arbeit wurde das Erstarrungsverhalten kolloidaler Scherschmelzen in den kristallinen Zustand mit (zeitaufgelöster) statischer Lichtstreuung und mikroskopischen Methoden untersucht. Dies erlaubte zunächst die kritische Überprüfung klassischer Modelle zur Kristallisationskinetik (Wilson- Frenkel- Gesetz, klassische Nukleationstheorie, Kolmogorov- Johnson- Mehl- Avrami (KJMA)- Modell). Es zeigte sich, dass diese Modelle gut geeignet sind die Verfestigung auch geladener kolloidaler Schmelzen zu beschreiben, wenn die diffusive Einteilchendynamik korrekt berücksichtigt wird. Erstmals wurden Oberflächenspannungen zwischen Kristallkeim und Schmelze für geladene Systeme bestimmt, die im Gegensatz zu Hartkugel- Systemen eine lineare Zunahme mit der Partikelkonzentration aufweisen. Der Methodenpark und die Auswerteverfahren wurden sodann auf binäre kolloidale Mischungen übertragen. Entsprechend den Einzelkomponenten kristallisieren alle Mischungen in einer kubischen Struktur. Leitfähigkeitsmessungen und Elastizität stehen meist im Einklang mit der Nukleation zufallsgeordneter Substitutionskristalle. Für mehrere Proben mit unterschiedlichen Größenverhältnissen wurde mit statischer Lichtstreuung der Einfluss der Komposition und der Partikelkonzentration auf das Nukleationsverhalten untersucht. Im Allgemeinen wurde das Nukleationsszenario einkomponentiger Systeme mit einigen unerwarteten, quantitativen Unterschieden reproduziert. Für eine Probe, die eine Kompositionsordnung andeutet, wurden interessante Korrelationen zwischen der Nukleationskinetik und den Eigenschaften des resultierenden Festkörpers gefunden.

Diffusionsgewichtete Helium-3 Magnetresonanztomographie zur Untersuchung der Lunge

Die diffusionsgewichtete Magnetresonanztomographie (MRT) mit dem hyperpolarisierten Edelgas-Isotop 3He ist ein neues Verfahren zur Untersuchung von Erkrankungen der Atem-wege und der Lunge. Die Diffusionsbewegung der 3He-Atome in den Luftwegen der Lunge wird durch deren Wände begrenzt, wobei diese Einschränkung sowohl von den Dimensionen der Atemwege als auch von den Messparametern abhängt. Man misst daher einen scheinbaren Diffusionskoeffizienten (Apparent Diffusion Coefficient, ADC) der kleiner ist als der Diffusionskoeffizient bei freier Diffusion. Der ADC gestattet somit eine qualitative Abschät-zung der Größe der Luftwege und deren krankhafte Veränderung, ohne eine direkte Abbil-dung der Luftwege selbst. Eine dreidimensionale Abbildung der räumlichen Verteilung von Lungenschädigungen wird dadurch möglich. Ziel der vorliegenden Arbeit war es, ein tieferes physikalisch fundiertes Verständnis der 3He-Diffusionsmessung zu ermöglichen und die Methode der diffusionsgewichteten 3He-MRT hin zur Erfassung des kompletten 3He-Diffusionstensors weiterzuentwickeln. Dazu wurde systematisch im Rahmen von Phantom- und tierexperimentellen Studien sowie Patientenmes-sungen untersucht, inwieweit unterschiedliche Einflussfaktoren das Ergebnis der ADC-Messung beeinflussen. So konnte beispielsweise nachgewiesen werden, dass residuale Luftströmungen am Ende der Einatmung keinen Einfluss auf den ADC-Wert haben. Durch Simulationsrechnungen konnte gezeigt werden, in welchem Maße sich die durch den Anregungspuls hervorgerufene Abnah-me der Polarisation des 3He-Gases auf den gemessenen ADC-Wert auswirkt. In einer Studie an lungengesunden Probanden und Patienten konnte die Wiederholbarkeit der ADC-Messung untersucht werden, aber auch der Einfluss von Gravitationseffekten. Diese Ergebnisse ermöglichen genauere Angaben über systematische und statistische Messfehler, sowie über Grenzwerte zwischen normalem und krankhaft verändertem Lungengewebe. Im Rahmen dieser Arbeit wurde die bestehende diffusionsgewichtete Bildgebung methodisch zur Erfassung des kompletten Diffusionstensors von 3He in der Lunge weiterentwickelt. Dies war wichtig, da entlang der Luftwege weitestgehend freie Diffusion vorherrscht, während senkrecht zu den Luftwegen die Diffusion eingeschränkt ist. Mit Hilfe von Simulationsrech-nungen wurde der kritische Einfluss von Rauschen in den MRT-Bildern auf die Qualität der Messergebnisse untersucht. Diese neue Methodik wurde zunächst an einem Phantom beste-hend aus einem Bündel aus Glaskapillaren, deren innerer Durchmesser mit dem des mensch-lichen Azinus übereinstimmt, validiert. Es ergab sich eine gute Übereinstimmung zwischen theoretischen Berechnungen und experimentellen Ergebnissen. In ersten Messungen am Menschen konnten so unterschiedliche Anisotropiewerte zwischen lungengesunden Proban-den und Patienten gefunden werden. Es zeigte sich eine Tendenz zu isotroper Diffusion bei Patienten mit einem Lungenemphysem. Zusammenfassend tragen die Ergebnisse der vorliegenden Arbeit zu einem besseren Ver-ständnis der ADC-Messmethode bei und helfen zukünftige Studien aufgrund des tieferen Verständnisses der die 3He Messung beeinflussenden Faktoren besser zu planen.

Numerical analysis of Dielectric Barrier Discharge

A numerical investigation of dielectric barrier discharge aimed to simulate the electro hydro dynamic interaction is presented. A discharge drift diffusive model according to the Townsend avalanche is described and used to duplicate the plasma kinetics of a DBD actuator. The discharge characteristics dependence upon dielectric material and applied voltage are simulated and the EHD force field according to a simplified approach is presented and discussed. The coupling of DBD results with a fluid dynamic code is also shown. Finally, a new non invasive diagnostic technique for EHD interaction based on Schlieren imaging is computationally validated.

Influence of spatial constraints in non-crystalline regions on the polymer dynamics in semi-crystalline polyethylene: a solid state NMR study

A broad variety of solid state NMR techniques were used to investigate the chain dynamics in several polyethylene (PE) samples, including ultrahigh molecular weight PEs (UHMW-PEs) and low molecular weight PEs (LMW-PEs). Via changing the processing history, i.e. melt/solution crystallization and drawing processes, these samples gain different morphologies, leading to different molecular dynamics. Due to the long chain nature, the molecular dynamics of polyethylene can be distinguished in local fluctuation and long range motion. With the help of NMR these different kinds of molecular dynamics can be monitored separately. In this work the local chain dynamics in non-crystalline regions of polyethylene samples was investigated via measuring 1H-13C heteronuclear dipolar coupling and 13C chemical shift anisotropy (CSA). By analyzing the motionally averaged 1H-13C heteronuclear dipolar coupling and 13C CSA, the information about the local anisotropy and geometry of motion was obtained. Taking advantage of the big difference of the 13C T1 relaxation time in crystalline and non-crystalline regions of PEs, the 1D 13C MAS exchange experiment was used to investigate the cooperative chain motion between these regions. The different chain organizations in non-crystalline regions were used to explain the relationship between the local fluctuation and the long range motion of the samples. In a simple manner the cooperative chain motion between crystalline and non-crystalline regions of PE results in the experimentally observed diffusive behavior of PE chain. The morphological influences on the diffusion motion have been discussed. The morphological factors include lamellar thickness, chain organization in non-crystalline regions and chain entanglements. Thermodynamics of the diffusion motion in melt and solution crystallized UHMW-PEs is discussed, revealing entropy-controlled features of the chain diffusion in PE. This thermodynamic consideration explains the counterintuitive relationship between the local fluctuation and the long range motion of the samples. Using the chain diffusion coefficient, the rates of jump motion in crystals of the melt crystallized PE have been calculated. A concept of "effective" jump motion has been proposed to explain the difference between the values derived from the chain diffusion coefficients and those in literatures. The observations of this thesis give a clear demonstration of the strong relationship between the sample morphology and chain dynamics. The sample morphologies governed by the processing history lead to different spatial constraints for the molecular chains, leading to different features of the local and long range chain dynamics. The knowledge of the morphological influence on the microscopic chain motion has many implications in our understanding of the alpha-relaxation process in PE and the related phenomena such as crystal thickening, drawability of PE, the easy creep of PE fiber, etc.

Development of parallelizable flood inundation models for large scale analysis

Flood disasters are a major cause of fatalities and economic losses, and several studies indicate that global flood risk is currently increasing. In order to reduce and mitigate the impact of river flood disasters, the current trend is to integrate existing structural defences with non structural measures. This calls for a wider application of advanced hydraulic models for flood hazard and risk mapping, engineering design, and flood forecasting systems. Within this framework, two different hydraulic models for large scale analysis of flood events have been developed. The two models, named CA2D and IFD-GGA, adopt an integrated approach based on the diffusive shallow water equations and a simplified finite volume scheme. The models are also designed for massive code parallelization, which has a key importance in reducing run times in large scale and high-detail applications. The two models were first applied to several numerical cases, to test the reliability and accuracy of different model versions. Then, the most effective versions were applied to different real flood events and flood scenarios. The IFD-GGA model showed serious problems that prevented further applications. On the contrary, the CA2D model proved to be fast and robust, and able to reproduce 1D and 2D flow processes in terms of water depth and velocity. In most applications the accuracy of model results was good and adequate to large scale analysis. Where complex flow processes occurred local errors were observed, due to the model approximations. However, they did not compromise the correct representation of overall flow processes. In conclusion, the CA model can be a valuable tool for the simulation of a wide range of flood event types, including lowland and flash flood events.

The role of the West African Monsoon in the tropical to mid-latitudes climate

The interaction between atmosphere–land–ocean–biosphere systems plays a prominent role on the atmospheric dynamics and on the convective rainfall distribution over the West Africa monsoon area during the boreal summer. In particular, the initialization of convective systems in the Sub – Sahelian region has been directly linked to soil moisture heterogeneities identified as the major triggering, development and propagation of convective systems. The present study aims at investigating African monsoon large scale convective dynamics and rainfall diurnal cycle through an exploration of the hypothesis behind the mechanisms of a monsoon phenomenon as an emergence of a collective dynamics of many propagating convective systems. Such hypothesis is based on the existence of an internal self – regulation mechanism among the various components. To achieve these results a multiple analysis was performed based on remote sensed rainfall dataset, and global and regional modelling data for a period of 5 seasons: 2004 - 2008. Satellite rainfall data and convective occurrence variability were studied for assessing typical spatio – temporal signatures and characteristics with an emphasis to the diurnal cycle footprint. A global model and regional model simulation datasets, specifically developed for this analysis and based on Regional Atmospheric Modelling System – RAMS, have been analysed. Results from numerical model datasets highlight the evidence of a synchronization between the destabilization of the convective boundary layer and rainfall occurrence due to the solar radiation forcing through the latent heat release. This supports the conclusion that the studied interacting systems are associated with a process of mutual adjustment of rhythms. Furthermore, this rainfall internal coherence was studied in relation to the West African Heat Low pressure system, which has a prominent role in the large scale summer variability over the Mediterranean area since it is acting as one of dynamic link between sub tropical and midlatitudes variability.

Diffusion of laser polarized gases in MRI

The use of Magnetic Resonance Imaging (MRI) as a diagnostic tool is increasingly employing functional contrast agents to study or contrast entire mechanisms. Contrast agents in MRI can be classified in two categories. One type of contrast agents alters the NMR signal of the protons in its surrounding, e.g. lowers the T1 relaxation time. The other type enhances the Nuclear Magnetic Resonance (NMR) signal of specific nuclei. For hyperpolarized gases the NMR signal is improved up to several orders of magnitude. However, gases have a high diffusivity which strongly influences the NMR signal strength, hence the resolution and appearance of the images. The most interesting question in spatially resolved experiments is of course the achievable resolution and contrast by controlling the diffusivity of the gas. The influence of such diffusive processes scales with the diffusion coefficient, the strength of the magnetic field gradients and the timings used in the experiment. Diffusion may not only limit the MRI resolution, but also distort the line shape of MR images for samples, which contain boundaries or diffusion barriers within the sampled space. In addition, due to the large polarization in gaseous 3He and 129Xe, spin diffusion (different from particle diffusion) could play a role in MRI experiments. It is demonstrated that for low temperatures some corrections to the NMR measured diffusion coefficient have to be done, which depend on quantum exchange effects for indistinguishable particles. Physically, if these effects can not change the spin current, they can do it indirectly by modifying the velocity distribution of the different spin states separately, so that the subsequent collisions between atoms and therefore the diffusion coefficient can eventually be affected. A detailed study of the hyperpolarized gas diffusion coefficient is presented, demonstrating the absence of spin diffusion (different from particle diffusion) influence in MRI at clinical conditions. A novel procedure is proposed to control the diffusion coefficient of gases in MRI by admixture of inert buffer gases. The experimental measured diffusion agrees with theoretical simulations. Therefore, the molecular mass and concentration enter as additional parameters into the equations that describe structural contrast. This allows for setting a structural threshold up to which structures contribute to the image. For MRI of the lung this allows for images of very small structural elements (alveoli) only, or in the other extreme, all airways can be displayed with minimal signal loss due to diffusion.

Continuous geochemical monitoring by mass-spectometer in the Campi Flegrei geothermal area. An application at Pisciarelli-Solfatara (diffuse and fumarolic gases) and at the mud gases during drilling of the CFDDP pilot hole

During this work has been developed an innovative methodology for continuous and in situ gas monitoring (24/24 h) of fumarolic and soil diffusive emissions applied to the geothermal and volcanic area of Pisciarelli near Agnano inside the Campi Flegrei caldera (CFc). In literature there are only scattered and in discrete data of the geochemical gas composition of fumarole at Campi Flegrei; it is only since the early ’80 that exist a systematic record of fumaroles with discrete sampling at Solfatara (Bocca Grande and Bocca Nuova fumaroles) and since 1999, even at the degassing areas of Pisciarelli. This type of sampling has resulted in a time series of geochemical analysis with discontinuous periods of time set (in average 2-3 measurements per month) completely inadequate for the purposes of Civil Defence in such high volcanic risk and densely populated areas. For this purpose, and to remedy this lack of data, during this study was introduced a new methodology of continuous and in situ sampling able to continuously detect data related and from its soil diffusive degassing. Due to its high sampling density (about one measurement per minute therefore producing 1440 data daily) and numerous species detected (CO2, Ar, 36Ar, CH4, He, H2S, N2, O2) allowing a good statistic record and the reconstruction of the gas composition evolution of the investigated area. This methodology is based on continuous sampling of fumaroles gases and soil degassing using an extraction line, which after undergoing a series of condensation processes of the water vapour content - better described hereinafter - is analyzed through using a quadrupole mass spectrometer

Dynamics, rheology and critical properties of colloidal fluid mixtures

Liquids under the influence of external fields exhibit a wide range of intriguing phenomena that can be markedly different from the behaviour of a quiescent system. This work considers two different systems — a glassforming Yukawa system and a colloid-polymer mixture — by Molecular Dynamics (MD) computer simulations coupled to dissipative particle dynamics. The former consists of a 50-50 binary mixture of differently-sized, like-charged colloids interacting via a screened Coulomb (Yukawa) potential. Near the glass transition the influence of an external shear field is studied. In particular, the transition from elastic response to plastic flow is of interest. At first, this model is characterised in equilibrium. Upon decreasing temperature it exhibits the typical dynamics of glassforming liquids, i.e. the structural relaxation time τα grows strongly in a rather small temperature range. This is discussed with respect to the mode-coupling theory of the glass transition (MCT). For the simulation of bulk systems under shear, Lees-Edwards boundary conditions are applied. At constant shear rates γ˙ ≫ 1/τα the relevant time scale is given by 1/γ˙ and the system shows shear thinning behaviour. In order to understand the pronounced differences between a quiescent system and a system under shear, the response to a suddenly commencing or terminating shear flow is studied. After the switch-on of the shear field the shear stress shows an overshoot, marking the transition from elastic to plastic deformation, which is connected to a super-diffusive increase of the mean squared displacement. Since the average static structure only depends on the value of the shear stress, it does not discriminate between those two regimes. The distribution of local stresses, in contrast, becomes broader as soon as the system starts flowing. After a switch-off of the shear field, these additional fluctuations are responsible for the fast decay of stresses, which occurs on a time scale 1/γ˙ . The stress decay after a switch-off in the elastic regime, on the other hand, happens on the much larger time scale of structural relaxation τα. While stresses decrease to zero after a switch-off for temperatures above the glass transition, they decay to a finite value for lower temperatures. The obtained results are important for advancing new theoretical approaches in the framework of mode-coupling theory. Furthermore, they suggest new experimental investigations on colloidal systems. The colloid-polymer mixture is studied in the context of the behaviour near the critical point of phase separation. For the MD simulations a new effective model with soft interaction potentials is introduced and its phase diagram is presented. Here, mainly the equilibrium properties of this model are characterised. While the self-diffusion constants of colloids and polymers do not change strongly when the critical point is approached, critical slowing down of interdiffusion is observed. The order parameter fluctuations can be determined through the long-wavelength limit of static structure factors. For this strongly asymmetric mixture it is shown how the relevant structure factor can be extracted by a diagonalisation of a matrix that contains the partial static structure factors. By presenting first results of this model under shear it is demonstrated that it is suitable for non-equilibrium simulations as well.

Ionenleitfähigkeit in polymeren Matrizes – Synthese, Charakterisierung und Untersuchung dynamischer Prozesse von neuen Lithium- und Protonenleitern

Die vorliegende Dissertation befasst sich mit der Synthese, physikochemischen und polymerspezifischen Charakterisierung und insbesondere der impedanzspektroskopischen Untersuchung von sowohl neuartigen, solvensfreien lithiumionen- als auch protonenleitfähigen Polymermaterialien für potentielle Anwendungen in sekundären Lithiumionenbatterien bzw. in Hochtemperatur-Protonenaustauschmembran-Brennstoffzellen (engl.: proton exchange membrane fuel cell, auch: polymer electrolyte membrane fuel cell, PEMFC). Beiden Typen von ionenleitfähigen Membranen liegt das gängige Prinzip der chemischen Anbindung einer für den Ionentransport verantwortlichen Seitengruppe an eine geeignete Polymerhauptkette zugrunde („Entkopplung“; auch Immobilisierung), welcher hinsichtlich Glasübergangstemperatur (Tg), elektrochemischer und thermischer Stabilität (Td) eine dynamisch entkoppelte, aber nicht minder bedeutsame Rolle zukommt. Die Transportaktivierung erfolgt in beiden Fällen thermisch. Im Falle der Protonenleiter liegt die zusätzliche Intention darin, eine Alternative aufzuzeigen, in der die Polymerhauptkette gekoppelt direkt am Protonentransportmechanismus beteiligt ist, d.h., dass der translatorisch diffusive Ionentransport entlang der Hauptkette stattfindet und nicht zwischen benachbarten Seitenketten. Ein Hauptaugenmerk der Untersuchungen liegt sowohl bei den lithiumionen- als auch den protonenleitfähigen Polymermembranen auf temperaturabhängigen dynamischen Prozessen der jeweiligen Ionenspezies in der polymeren Matrix, was die Ionenleitfähigkeit selbst, Relaxationsphänomene, die translatorische Ionendiffusion und im Falle der Protonenleiter etwaige mesomere Grenzstrukturübergänge umfasst. Lithiumionenleiter: Poly(meth)acrylate mit (2-Oxo-1,3-dioxolan)resten (Cyclocarbonat-) in der Seitenkette unterschiedlicher Spacerlänge wurden synthetisiert und charakterisiert. Die Leitfähigkeit s(,T) erreicht bei Poly(2-oxo-[1,3]dioxolan-4-yl)methylacrylat (PDOA): Lithium-bis-trifluormethansulfonimid (LiTFSI) (10:3) ca. 10^-3,5 S cm^-1 bei 150 °C. Weichmachen (Dotieren) mit äquimolaren Mengen an Propylencarbonat (PC) bewirkt in allen Fällen einen enormen Anstieg der Leitfähigkeit. Die höchsten Leitfähigkeiten von Mischungen dieser Polymere mit LiTFSI (und LiBOB) werden nicht beim System mit der niedrigsten Tg gefunden. Auch dient Tg nicht als Referenztemperatur (Tref) nach Williams-Landel-Ferry (WLF), so dass eine WLF-Anpassung der Leitfähigkeitsdaten nur über einen modifizierten WLF-Algorithmus gelingt. Die ermittelten Tref liegen deutlich unterhalb von Tg bei Temperaturen, die charakteristisch für die Seitenkettenrelaxation sind („Einfrieren“). Dies legt nahe, dass der Relaxation der Seitenketten eine entscheidende Rolle im Li^+-Leitfähigkeitsmechanismus zukommt. Die Li^+-Überführungszahlen tLi^+ in diesen Systemen schwanken zwischen 0,13 (40 °C) und 0,55 (160 °C). Protonenleiter: Polymere mit Barbitursäure- bzw. Hypoxanthinresten in der Seitenkette und Polyalkylenbiguanide unterschiedlicher Spacerlänge wurden synthetisiert und charakterisiert. Die Leitfähigkeit s(,T) erreicht bei Poly(2,4,6(1H,3H,5H)-trioxopyrimidin-5-yl)methacrylat (PTPMA) maximal ca. 10^-4,4 S cm^-1 bei 140 °C. Höhere Leitfähigkeiten sind nur durch Mischen mit aprotischen Lösungsmitteln erreichbar. Die höchste Leitfähigkeit wird im Falle der Polyalkylenbiguanide bei Polyethylenbiguanid (PEB) erzielt. Sie erreicht 10^-2,4 S cm^-1 bei 190 °C. Die Aktivierungsenergien EA der Polyalkylenbiguanide liegen (jeweils unterhalb von Tg) zwischen ca. 3 – 6 kJ mol^-1. In allen beobachteten Fällen dient Tg als Tref, so dass eine konventionelle WLF-Behandlung möglich ist und davon auszugehen ist, dass die Leitfähigkeit mit dem freien Volumen Vf korreliert.