Universität Stuttgart
Permanent URI for this communityhttps://elib.uni-stuttgart.de/handle/11682/1
Browse
29 results
Search Results
Item Open Access Nichtlineare Halbgruppen und Hysteresemodelle für nichtmonotone Sättigungsprofile(2019) Steinle, Rouven; Hilfer, Rudolf (Prof. Dr. Dr.)Die experimentellen Beobachtungen von Sättigungsüberschüssen und Sättigungsfingern bei gravitationsgetriebenen Infiltrationsprozessen beschäftigt schon seit Jahrzehnten die Forschung von porösen Medien in Mathematik, Physik, Geowissenschaft und Ingenieurwesen. Sowohl die Vorhersage von Sättigungsüberschüssen bei Experimenten als auch das allgemeine theoretische Verständnis von Sättigungsüberschüssen ist noch nicht vollständig geklärt. Des Weiteren ist die Frage, ob Sättigungsüberschüsse und Sättigungsfinger eine Verbindung zueinander haben, weiterhin offen. Zur theoretischen Modellierung von nichtmonotonen Sättigungsprofilen und den Mechanismen, welche zur Entstehung und Ausbreitung dieser Profile nötig sind, gibt es eine Vielzahl von Ansätzen, beschrieben in Beliaev et al (Transport in Porous Media 43 (2001), S.487), Cuesta et al (European Journal of Applied Mathematics 11 (2000), S.381), Cueto-Felgueroso et al (Physical Review Letters 101 (2008), S.244504), DiCarlo et al (Transport in Porous Media 91 (2012), S.955), Doster et al (Physical Review E 86 (2012), S.016317) und Eliassi et al (Water Resources Research 38 (2002), S.1234). Ein Ansatz zur Ausbreitung von laufenden nichtmonotonen Sättigungsprofilen wurde von Hilfer und Steinle (The European Physical Journal Special Topics 223 (2014), S.2323) vorgestellt. Dieser Ansatz benutzt ausschließlich hysteretische relative Permeabilitäten zur Modellierung von Be- und Entwässerungsprozessen bei sich ausbreitenden Sättigungsüberschüssen. Sonstige Erweiterungen der Darcy-Theorie werden dabei nicht benutzt. Die Ausbreitungsgeschwindigkeiten der Sättigungsfronten werden mit einer Rankine-Hugoniot-Bedingung, welche die fraktionalen Flüsse verwendet, berechnet. Bei einer optimal aufeinander abgestimmten Wahl der relativen Permeabilitäten, der Material- und Fluidparameter sowie der Anfangs- und Randbedingungen, existieren Sättigungsüberschüsse, die sich mit konstanter Geschwindigkeit und konstanter Überschussbreite ausbreiten. In der vorliegenden Arbeit wird der Ansatz von Hilfer und Steinle (The European Physical Journal Special Topics 223 (2014), S.2323) auf die Darcy-Theorie in Verbindung mit einem allgemeinen Hysteresemodell mit graphenabhängigen relativen Permeabilitäten und Kapillardrücken angewendet. Das entstandene Anfangs- und Randwertproblem wird analytisch und numerisch untersucht. Im analytischen Teil wird sowohl das traditionelle als auch das hysteretische Anfangs- und Randwertproblem der Darcy-Theorie mit den Methoden der nichtlinearen Halbgruppentheorie diskutiert. Für die traditionelle Darcy-Theorie mit einer geeigneten Anfangsbedingung kann Existenz, Eindeutigkeit und Beschränktheit von Lösungen gezeigt werden. Darauf aufbauend werden diese Ergebnisse auf die hysteretische Darcy-Theorie angewendet. Im Rahmen einer vereinfachten Hystereseformulierung und mit geeigneten Bedingungen an die Anfangsbedingung und die Parameter der Gleichung existieren auch hier eindeutig Lösungen. Bei einer nichtmonotonen Anfangsbedingung sind diese Lösungen auch nichtmonoton. Des Weiteren wird die Ausbreitung von nichtmonotonen Sättigungsprofile mit Hilfe der Rankine-Hugoniot-Bedingung untersucht. Die Wahl der Parameter und Anfangs- und Randbedingungen erzeugt nicht nur laufende Sättigungsüberschüsse, sondern auch Überschüsse mit wachsender oder schrumpfender Plateaubreite. Nicht nur die Wahl der Parameter hat einen Einfluss auf die Ausbreitung der Lösungen, sondern auch die Wahl des Hysteresemodells. Abhängig vom Hysteresemodell, kann sich ein nichtmonotones Sättigungsprofil als laufende Welle oder mit wachsender oder schrumpfender Plateaubreite durch das poröse Medium bewegen. Außerdem kann man in der theoretischen Betrachtung nicht nur die Ausbreitung von Sättigungsüberschüssen bei schwerkraftgetriebenen, sondern auch bei horizontalen Strömungen erklären. Im numerischen Teil wird die hysteretische Darcy-Theorie für variierende Material- und Fluidparameter, Anfangs- und Randbedingungen und Hysteresemodelle untersucht. Die numerischen Simulationen bestätigen die analytischen Resultate und zeigen, dass die Lösungen unabhängig von der Hysterese im Kapillardruck und im kapillaren Fluss sind. Jedoch sind die Lösungen stark abhängig von den Anfangs- und Randwerten. In Folge dieser Abhängigkeit gibt es eine große Variation an qualitativ unterschiedlichen zeitlichen Entwicklungen der Sättigungsüberschüsse. Darüber hinaus werden existierende Experimente modelliert. Es kann gezeigt werden, dass die numerischen nichtmonotonen Sättigungsprofile mit den experimentellen Daten in Bezug auf Ausbreitungsgeschwindigkeiten, Überschussbreite und Überschusshöhe übereinstimmen. Des Weiteren wird numerisch untersucht, ob die Formulierung der Darcy-Theorie einen Einfluss auf die Lösungen hat. Bei einer geeigneten Wahl der Randbedingungen sind die numerischen Ergebnisse nahezu identischen und somit sind die Lösungen unabhängig von der Formulierung. Abschließend ist eine erste numerische Einordnung von Sättigungsfingern bei zweidimensionalen Strömungen möglich. Die numerischen Lösungen zeigen, dass nur ein Sättigungsfinger in Verbindung mit einem Sättigungsüberschuss in der Fingerspitze für lange Zeiten existiert. Diese Arbeit zeigt, dass die Beschreibung von laufenden Sättigungsüberschüssen im Rahmen der Darcy-Theorie mit hysteretischen relativen Permeabilitäten möglich ist. Zusätzlich zu laufenden Sättigungsüberschüssen existieren in der Theorie auch nichtlaufende nichtmonotone Sättigungsprofile mit wachsender oder schrumpfender Plateaubreite. Einerseits sind in der theoretischen Modellierung weitere Untersuchungen und numerische Simulationen für zweidimensionale Strömungen notwendig, um Aussagen über die Ausbreitung von Sättigungsfingern treffen zu können. Auf der anderen Seite sind auch weitere Forschungen mit längeren quasi-eindimensionalen porösen Medien notwendig, um herauszufinden, ob die nichtlaufenden Überschüsse aus der Theorie auch experimentell beobachtbar sind. Darüber hinaus wurden bei einer Kooperation mit dem Institut für Wasser- und Umweltsystemmodellierung (IWS) der Universität Stuttgart nichtmonotone Sättigungsprofile mit unterschiedlichen Implementierungen berechnet und verglichen, Schneider et al (Transport in Porous Media 121 (2018), Nr.3). Die numerischen Lösungen unterscheiden sich dabei nur minimal und damit konnte gezeigt werden, dass die Ausbreitung von Sättigungsüberschüssen nicht von der Wahl der numerischen Lösungsmethode abhängig ist.Item Open Access Simulation of novel magnetic materials in the field of soft matter(2014) Weeber, Rudolf; Holm, Christian (Prof. Dr.)This thesis has dealt with the tailoring of magnetic soft matter. Two strategies are available to achieve this goal. First, it is possible to alter the magnetic nanoparticles, in order to change their interactions. Second, it is possible to exchange the carrier fluid into which the magnetic particles are embedded by a more complex matrix. For each of these two possibilities, an example was studied, namely shifted-dipole particles and magnetic gels. Shifted-dipole particles (sd-particles) are a special kind of model magnetic particles which can be used to explain findings for particles with magnetic caps as well as for particles with magnetic inclusions. Magnetic gels, on the other hand, derive their particular properties from an interplay of the magnetic properties of the nanoparticles and the elastic behavior of the polymer matrix. In the thesis, the findings for these two systems will be discussed, and further research questions will be identified.Item Open Access Lattice Boltzmann simulations of fluid flow in the vicinity of rough and hydrophobic boundaries(2010) Kunert, Christian; Harting, Jens (P.D. Dr.)In recent years, it became possible to perform very well controlled experiments that have shown a violation of the no-slip boundary condition in sub-micron sized geometries. Since then, mostly experimental, but also theoretical works, as well as computer simulations, have been performed to improve our understanding of boundary slip. The topic is of fundamental interest because it has practical consequences in the physical and engineering sciences as well as for medical and industrial applications. This work focuses on numerical investigations of the slip phenomenon by means of lattice Boltzmann simulations with a strong focus on roughness and the interplay between roughness and wetting phenomena. To do so, two different slip measurement methods are simulated. One is to apply a Poiseuille flow between two patterned boundaries, and to record the flow profile. Then, the profile can be compared to the theoretical one, which assumes a slip boundary condition. The second method records the drag force that is acting on a sphere which is moved with a constant velocity towards the observed surface. Due to the influence of the boundary, the drag force acting on the sphere is disturbed and a correction function is needed to describe the measured force.Item Open Access Untersuchung der Pfropfendynamik und -stabilität bei der vertikalen und horizontalen Pfropfenförderung(2005) Strauß, Martin; Herrman, Hans J. (Prof. Dr.)Die Pfropfenförderung stellt ein produktschonendes Transportverfahren für granulare Medien dar, dass in der Industrie breite Verwendung findet. Eine zuverlässige Voraussage des Massentransportes bei der Auslegung industrieller Anlagen ist beim derzeitigen Wissensstand jedoch nicht möglich, da Grundkenntnisse über diesen Fördertyp fehlen. Ursache ist, dass sich die Untersuchung der bei der Pfropfenförderung mittels Druckluft durch Transportröhren getriebenen Granulatpfropfen sich weitestgehend experimentellen Messmethoden entzieht. Um Zugang zu den Prozessen während des Transportes zu erhalten, wurde für diese Arbeit die Förderung mittels einer Kombination aus Molekulardynamiksimulation und einem Löser für den Druckverlust am System auf dem Computer nachgebildet. Die Gültigkeit des Ansatzes wurde mittels experimenteller Messdaten an einer Anlage im Maßstab 1 zu 1 nachgewiesen. Es wurden sowohl Parameterstudien für die vertikale als auch die horizontale Förderung durchgeführt. Die Ergebnisse zeigen unter anderem, dass die Stabilität der Pfropfen durch entgegengesetzt auf den Pfropfen wirkende Kräfte bewirkt wird. Maßgeblich sind dabei die vorantreibende Druckluft und die auf den Pfropfen treffenden bremsenden Granulatteilchen zwischen den Pfropfen. Der Effekt ist unabhängig von der Transportrichtung. Unterschiede im Längenwachstum und der inneren Dynamik der Pfropfen zwischen vertikaler und horizontaler Förderung zeigen, dass eine einheitliche Abhandlung der Förderungsrichtungen bei der Auslegung von industriellen Anlagen nicht sinnvoll ist.Item Open Access The electrophoretic mobility of bare and soft spherical colloids : a molecular dynamics study(2014) Raafatnia, Shervin; Holm, Christian (Prof. Dr.)In this work, the electrophoretic mobility of colloids in salt solutions are studied by means of coarse-grained Molecular Dynamics simulations. Two different types of colloids are considered; bare colloids and polyelectrolyte-grafted colloids. A novel model for simulation of large bare colloids in the presence of explicit ions is developed. Comparison of the results with experimental data helps gain a better understanding of the mechanisms responsible for the interesting phenomenon of mobility reversal. Furthermore, a hitherto unknown electrokinetic behavior of polyelectrolyte-grafted colloids is found from simulations including full hydrodynamic interactions. The validity of the existing theories is verified via comparison with simulation results.Item Open Access Simulations on evolutionary phenomena with ageing models(2006) Schwämmle, Veit; Herrmann, Hans (Prof.)The idea of reducing the characteristics of highly complex systems to some basic features without losing the essential informations has been successful for the comprehension of not only physical systems. The experience in treating systems which consist of many elements that interact over different time scales as well as over wide distances, has opened many interdisciplinary fields to statistical physicists. Although animals and human beings behave in a much more complex way than for instance atoms or molecules, simple statistical models are able to reproduce many biological phenomena. The goal is to reproduce the collective behaviour of a large number of organisms: their individual properties are not of crucial importance. Due to the strongly increasing computer power along the last fifty years, computational models have gained more and more importance for research in many areas. The limitation of analytical models to describe accurately non-linear systems with critical behavior makes computational models nearly irreplaceable.Item Open Access Sand dunes on Mars and on Earth(2007) Ribeiro Parteli, Eric Josef; Herrmann, Hans Jürgen (Prof. Dr.)In this work the dune model introduced by Sauermann et al. (2001) is extended and applied to investigate the formation of different dune shapes on Mars and on Earth as function of wind directionality and sand availability. The formation of sand dunes on Mars under the present atmospheric conditions of the red planet is studied and conclusions about wind speed, migration velocity of dunes and changing wind regimes on Mars are presented. Field measurements of the shape of coastal transverse dunes are presented and the formation of coastal dune fields is explained. Finally, the formation of linear dunes by bimodal wind regimes is calculated. The simulations explain the appearance of exotic bimodal dune shapes in areas of low sand availability on Mars and on Earth.Item Open Access Entropic segregation of polymers under confinement(2016) Minina, Elena; Holm, Christan (Prof. Dr.)Overlapping polymers confined in a cylinder experience strong repulsion that drives them towards segregation. This has biological relevance to chromosome segregation in single-celled elongated bacteria such as Escherichia coli because in principle, chromosomes can segregate for purely entropic reasons without any help from active mechanisms. In this thesis, we investigated entropic segregation of polymers under cylindrical confinement of infinite length where the confining cylinder is so narrow that its diameter is significantly smaller than the radius of gyration of the unconfined polymers.Item Open Access Simulation of peloids(2007) Hecht, Martin; Herrmann, Hans J. (Prof. Dr.)In this work we investigate dense colloidal suspensions of alumina particles, which we regard as a model system for clay-like soils (peloids). Beyond soil mechanics these suspensions are important for wet processing of work-pieces in ceramics. In the present work we investigate the rheological properties of the suspensions and connect them to the underlying microstructure. Experimental findings depending on the pH value and the salt concentration in the sample (expressed as ionic strength) are reviewed (viscosity versus shear rate, shear thinning, oedometer data, sedimentation experiments, cyclic loading...) and computer simulations are performed to carry out further research work. For this purposes a coupled Molecular Dynamics (MD) and Stochastic Rotation Dynamics (SRD) code has been developed and the simulation results are compared to the experimental data. To describe the surface charge of the particles a charge regulation model, which describes adsorption and desorption of the charge determining ions on the particle surface, has been developed within Debye Huckel theory. This model has been calibrated to measurements of the zeta potential. Starting from the known values, the model allows us to extrapolate to different experimental conditions. Based on characteristic time scales and dimensionless numbers like the Reynolds number and the Peclet number, we apply a scaling scheme to determine the simulation parameters, so that we can achieve a quantitative comparability of simulation and experiment. Using the shear viscosity, shear force, pair correlation function, density fluctuations and structure factor we can identify three different regimes: a clustered regime, a stable suspension, a repulsive structure similar to the structure known from glassy systems. The microstructures are plotted in a stability diagram depending on pH value an ionic strength. The microstructures found in the simulations provide a possibility to explain the relations found in the experiments.Item Open Access Laufende Wellenlösungen von Systemen nichtlinearer partieller Differentialgleichungen am Beispiel von Mehrphasenströmungen in porösen Medien(2012) Hönig, Oliver; Hilfer, Rudolf (Prof. Dr. Dr.)Obwohl die Anwendungen der Theorie der Mehrphasenströmungen in porösen Medien äußerst zahlreich und in den letzten hundert Jahren intensiv erforscht worden sind, ist das Verhalten auf makroskopischer Skala bis heute nicht verstanden. Darunter sind insbesondere die Hysterese und die residualen Fluidkonfigurationen zu zählen. Hilfer stellt in Physical Review E, Vol. 73, 016307, 2006 ein Modell vor, das die zwei Phasen in jeweils perkolierende und nichtperkolierende Phasen unterteilt. Seitdem konnten Doster und Hilfer erste vielversprechende Ergebnisse für das Modell finden. Darunter waren z.B. hyperbolische Ähnlichkeitslösungen, wie das Buckley-Leverett-Problem für immobile nichtperkolierende Phasen. Eine andere Ähnlichkeitslösung, nämlich die der laufenden Wellen, soll hier in dieser Arbeit auf das Modell angewandt werden. Laufende Wellen zeichnen sich dadurch aus, dass die Sättigungsprofile sich in der Form nicht ändern und sich mit konstanter Geschwindigkeit fortbewegen. Sie bedient sich der Theorie der dynamischen Systeme, um eine qualitative Beschreibung zu ermöglichen. Laufende Wellen sind zum einen interessant, weil sie in verschiedenen Experimenten beobachtet wurden, zum anderen aber auch, weil sie gewisse Schwierigkeiten im Buckley-Leverett-Problem beheben können. Das traditionelle und das in Physical Review E, Vol. 73, 016307, 2006 vorgeschlagene Modell wird in einer Raumdimension vorgestellt. Für das neue Modell werden dabei einige Näherungen besprochen, um Systeme einer, zweier oder dreier dimensionsloser fraktionaler Flussgleichungen zu erhalten. Die aus einem laufenden Wellenansatz dieser Systeme gewonnen Gleichungen werden mit einer dynamischen Systembeschreibung besprochen. Zuletzt werden Lösungen gezeigt, die direkt aus den partiellen Differentialgleichungen durch einen numerischen Algorithmus berechnet werden können. Hierbei werden geeignete Anfangs- und Randwerte gesucht, die stabile laufende Wellen ermöglichen. In Transport in Porous Media, Vol. 44, pp. 507-537, 2001 wurden von Brevdo et al. vier topologisch unterschiedliche laufende Wellenprofilklassen für das traditionelle Modell gefunden. In dieser Dissertation wird nun gezeigt, dass die Annahme, dass der Wasserfluss überall zu jeder Zeit stetig ist, zu Bedingungen führt, die die nichtdifferenzierbaren Wellenprofilklassen entweder auf ein bezüglich der Randwertsättigungen und der Geschwindigkeit eindeutiges Profil reduziert oder aber komplett ausschließt. Dies hängt vom Verhalten der Flussfunktionen und der Kapillarfunktionen ab. Bei den Parametrisierungsmodellen nach Brooks und Corey und nach Van Genuchten werden die nichtdifferenzierbaren Lösungsklassen ausgeschlossen. Die überall stetig differenzierbare Klasse besitzt die Eigenschaft, dass ihre Profile eindeutig durch die Randwertsättigungen und die Geschindigkeiten bestimmt werden. Damit ist es in dieser Arbeit gelungen, alle qualitativen Informationen der laufenden Wellenlösung in nur einer Abbildung darzustellen. Die Ergebnisse des Perkolationsmodell mit immobilen nichtperkolierenden Sättigungen unterscheiden sich im Vergleich zum traditionellen Modell darin, dass primäre und sekundäre Bewässerungen unterschieden werden können, was zu verschiedenen Geschwindigkeiten führt. Desweiteren ist es nun möglich, ein komplett trockenes Medium mit einer überall stetig differenzierbaren Welle zu bewässern. Bei mobilen nichtperkolierenden Phasen müssen zusätzliche Näherungen getroffen werden, um eine dynamische Systembeschreibung des laufenden Wellenansatzes zu ermöglichen. So muss die Kapillarität vereinfacht werden. Es kann nun zu nichtmonotonem Verhalten in allen vier Phasen kommen. Außerdem können sich Überschusslösungen in der Wassersättigung ergeben. Es werden für verschieden starke Näherungen Lösungen gezeigt, die die Auswirkungen einzelner Aspekte, wie die der gekoppelten Flüsse oder die des Massenaustausches, beleuchten. Die gekoppelten Flüsse spielen gerade bei den Überschusslösungen in der Wassersättigung eine herausragende Rolle. Die Massenaustauschterme können je nachdem, ob sie als Quellen oder Senken dienen, Überschüsse oder Unterschüsse in den einzelnen Phasen zulassen. Da sie jedoch nur zwischen den perkolierenden und nichtperkolierenden Phasen einer Flüssigkeit vorkommen, lassen sie die Wassersättigungsprofile monoton. Mit Hilfe eines numerischen Algorithmus' lassen sich geeignete Anfangs- und Randwerte finden, die in stabilen laufenden Wellenlösungen enden. Bei weniger restriktiven Annahmen kommen entweder ein oder zwei Freiheitsgrade hinzu. Dabei ist die Eindeutigkeit für gegebene Randwerte im Gegensatz zu den Lösungen für das traditionelle Modell nicht mehr gegeben. Dies geht auf die zusätzlichen Freiheitsgrade in den Anfangsbedingungen dreier oder vierer Phasen zurück, die entstehen, indem man die Wellen der einzelnen Phasen relativ zueinander verschieben kann.
- «
- 1 (current)
- 2
- 3
- »