Rheology of Brownian Discs
Dateien
Datum
Autor:innen
Herausgeber:innen
ISSN der Zeitschrift
Electronic ISSN
ISBN
Bibliografische Daten
Verlag
Schriftenreihe
Auflagebezeichnung
URI (zitierfähiger Link)
Internationale Patentnummer
Link zur Lizenz
Angaben zur Forschungsförderung
Projekt
Open Access-Veröffentlichung
Sammlungen
Core Facility der Universität Konstanz
Titel in einer weiteren Sprache
Publikationstyp
Publikationsstatus
Erschienen in
Zusammenfassung
In this thesis the properties of binary mixtures of hard discs, undergoing Brownian motion, have been studied. Two major cases have been considered: First the approach of the glass transition for a binary mixture along with the effect of changes in the composition within that mixture. For one selected mixture a detailed discussion of its transition followed. Second the rheological properties of this selected mixture were investigated by focusing on the distorted microstructure and the dynamical correlation functions of mainly tagged quantities. Finally time dependent shear flows were investigated. Performing simulations for different binary mixtures at the glass transition was motivated by recent results of MCT for two dimensions. With the simulations it was possible to confirm four mixing effects, predicted by MCT for these binary mixtures, by David Hajnal. It could be shown that for big radius ratios when we keep the total packing fraction constant, the increase in the concentration of the small particles leads to a melting of the glass, the plasticization. For small radius ratios the increase in the concentration of small particles at constant volume fraction leads to an even stronger glass. Both remaining effects, the increase of non-ergodicity parameters, and the slowing down of the relaxation towards the non-ergodicity parameters on increasing the small particles concentrations could be confirmed. These qualitative effects, especially the first two ones, are not only of theoretical interest, as industrial applications of so called plasticizers show. Going into detail for a selected mixture made it possible to determine the ideal MCT glass transition point. Along with the verification of the factorization close to the critical point, the α-scaling and even the increase of the plateaus according to a square root law could be validated. The simulation shows additional decay processes, leading to a decay of the correlation functions even above the ideal glass transition. This is not in accordance with MCT, however, the rise of the plateau values can be explained. This analysis closes a gap in the field of the glass transition as: First, a quantitative test for exactly the model system MCT uses was performed, and second to the knowledge of the author, no such analysis for a this two dimensional glass former has been performed before. With the characterization of the considered system the foundation for the rheological part was laid. That the simulation algorithm yields Brownian motion for the quiescent system was known before from the work of Erik Lange. This thesis goes one step further and shows, with the same arguments, that even for the sheared case the algorithm is still in accordance with short time expansions of the Smoluchowski equation for the shear modulus. The significance here is, that even on the two particle level used in the theoretical derivation, the simulation agrees in its short time asymptote with the theoretical findings. Having ensured that the simulation actually solves the Smoluchowski equation for Brownian motion under shear, and having characterized the glass transition for one special mixture in the preceding chapters, the rheological properties of exactly that mixture could be probed. In the framework of an extended MCT, MCT-ITT, which was developed by Matthias Fuchs and Michael Cates, the microscopic structure was investigated. The appearance of a yield stress, shear thinning and the distortion of the structure factor in the simulation can be qualitatively explained by MCT-ITT and the interplay of different relaxation time scales involved in the glass and in the liquid. Quantitatively we found that MCT-ITT overestimates the distortion of the structure and the resulting shear stresses and viscosities by about a factor of 10. However, it should be remarked that a simplified monodisperse MCT-ITT calculation was used, which could explain a part of the deviation. Future works should consider using a multi-component approach, as used for the analysis of the mixture effects in this thesis. Scrutinizing the tagged particle correlation functions, the interplay between purely structural relaxation and shear induced relaxation could be observed in detail in the simulation: For correlators in the liquid, the shear induced decay competes with the intrinsic decay, whereas in the glassy system the time scale is always set by the shear field. The relaxation in the shear melted glass thus follows a master function which could be found for very low shear rates in the simulation. The famous Taylor dispersion, expressing itself, for example, through a cubic growth in the mean squared displacement in the shear direction could be found in the simulations. The MCT-ITT prediction, that actually the shear and the gradient direction are connected via one single long time diffusion coefficient, could be confirmed. Furthermore the MCT-ITT scaling properties of these coefficients, stating that they scale linearly with the shear rate in the limit of vanishing shear, could be worked out with the simulation. Yet, the Taylor dispersion still leaves a riddle, as next-to-leading order asymptotes in MCT-ITT and the simulation seem to differ: A term proportional to the square of time couldn’t be found as proposed by MCT-ITT in the simulations, which propose a next-to-leading order dependence that is linear in time. Along with the discussion of the mean squared displacement, a connection between the superdiffusivity and the overshoot in the transient stress could be pointed out. Also the stress overshoot is connected to the anisotropy in the correlation functions for certain directions. Coherent correlators for the stationary state have been shown. In principle they show that simulating this quantities, although it requires a lot more numerical tricks and computational power, is possible. Calculating the transient correlators, which are the main quantity in MCT-ITT, will be a task for the near future, and will make further tests possible. Finally in the last chapter with a method, developed by Matthias Krüger, it was possible to close the gap between stationary and transient correlation functions. The simulation showed that for small waiting times the connection given between transient, stationary and equilibrium correlator is perfectly reproduced by Krüger’s relation. Furthermore, it could be shown, that even for large waiting times reasonable results can be achieved, although the shear modulus entered as a fitting parameter. Extending the simulation to oscillatory shear a connection between the yield stress of a glassy system and its non-linear reaction upon oscillatory shear for low frequencies could be inferred. The results support the theoretical results from a schematic MCT model invented by Joseph Brader. For time dependent stresses below the yield stress the system shows the elastic behavior of a solid state body. Above the yield stress the system shows the dissipative properties of a liquid accompanied by higher harmonics in the stress response. An upper bound for the stress can be estimated from the flow curves.
Zusammenfassung in einer weiteren Sprache
In dieser Arbeit wurden die Eigenschaften von Mischungen harter Scheiben mit Brownscher Molekularbewegung untersucht. Dabei wurden zwei prinzipielle Fälle betrachtet: Erstens wurde der Glasübergang einer binären Mischung unter Berücksichtigung der Zusammensetzung dieser Mischung untersucht. Detailliert wurde dabei der Glasübergang für eine ausgewählte Mischung diskutiert. Zweitens wurden die Fließeigenschaften dieser ausgewählten Mischung untersucht, wobei der Schwerpunkt auf der verzerrten Mikrostrukutur und den dynamischen Korrelationsfunktionen, hierbei insbesondere bei den inkohärenten Korrelatoren, lag. Kürzlich erzielte Ergebnisse aus der Modenkopplungstheorie haben den Anstoß zu den Simulationen für Mischungen binärer harter Scheiben gegeben. Mit den Simulationen konnten vier Mischungseffekte, die von den Modenkopplungsergebnissen von David Hajnal vorhergesagt wurden, bestätigt werden. Es wurde gezeigt, dass bei konstantem totalen Packungsbruch, für große Radienverhältnisse, das Erhöhen der Konzentration der kleinen Partikel, zum Schmelzen des Glases, zur sogenannten Plastifizierung, führt. Für kleine Radienverhältnisse führt das Erhöhen der Konzentration der kleinen Partikel bei festem Packungsbruch zu einer Verstärkung des Glases. Die beiden weiteren verbleibenden Effekte, das Ansteigen der Nichtergodizitätsparameter und das Verlangsamen der Relaxation auf die Nichtergodizitätsparameter beim Erhöhen der Konzentration der kleinen Partikel konnten bestätigt werden. Diese qualitativen Effekte, speziell die ersten beiden, sind nicht nur von Interesse für die Theorie, wie die industrielle Anwendung sogenannter Plastifizierer oder Weichmacher zeigt. Die detaillierte Untersuchung einer ausgewählten Mischung ermöglichte es, für diese den idealen Glasübergangspunkt innerhalb der Modenkopplungstheorie zu bestimmen. Zusammen mit der Bestätigung der Faktorisierung, dem α-Skalengesetz konnte sogar das Ansteigen der Plateaus gemäß eines Wurzelgesetzes gezeigt werden. Die Simulation enthält weitere Zerfallsprozesse, die zu einem Zerfall der Korrelatoren oberhalb des Glasübergangspunkts führen. Dies ist nicht im Einklang mit der Modenkoppungstheorie, jedoch kann der Anstieg der Plateaus erklärt werden. Diese Analyse schließt eine Lücke auf dem Gebiet des Glasübergangs: Erstens konnte genau das Modellsystem, das die Modenkopplungstheorie benutzt, einer quantitativen Prüfung unterzogen werden und zweitens wurde, nach Kenntnis des Autors, bis jetzt keine derartige Analyse für einen solchen zweidimensionalen Glasbildner durchgeführt. Mit der Charakterisierung dieses Systems wurde der Grundstein für eine Analyse des Fließverhaltens gelegt. Dass die Simulation Brownsche Dynamik für das ungescherte System erzeugt, wurde bereits in der Arbeit von Erik Lange gezeigt. Diese Arbeit geht einen Schritt weiter, und zeigt mit den selben Argumenten, dass der Algorithmus sogar im gescherten Fall mit der Kurzzeitentwicklungen der Smoluchowski-Gleichung für den Schermodul übereinstimmt. Bedeutend ist dabei, dass bei der theoretischen Betrachtung, die auch Zweipartikelwechselwirkungen mit einbezieht, die Kurzzeitasymptote mit denen der Simulationen übereinstimmt. Nachdem sichergestellt wurde, dass die Simulation tatsächlich die Smoluchowski-Gleichung der Brownschen Bewegung unter Scherung löst, und der Glasübergang der speziellen Mischung der vorangegangenen Kapitel charakterisiert wurde, konnte das Fließverhalten genau dieser Mischung untersucht werden. Die mikroskopische Struktur wurde im Rahmen einer erweiterten Modenkopplungstheorie, die von Matthias Fuchs und Michael Cates entwickelt wurde, untersucht. Das Auftreten einer Fließspannung, der Scherverdünnung und die Verzerrung des Strukturfaktors können qualitativ mit der erweiterten Modenkopplungstheorie und den verschiedenen, unter Scherung auftretenden Zeitskalen erklärt werden. Die erweiterte Modenkopplungstheorie überschätzt die Verzerrung der Struktur, woraus ungefähr um den Faktor zehn zu große Scherviskositäten und Scherspannungen folgen. Es sollte hinzugefügt werden, dass in diesem Fall jedoch ein vereinfachtes einkomponentiges Modell für die Modenkopplungsrechnungen verwendet wurde, was einen Teil der Abweichungen erklären kann. Zukünftige Arbeiten sollten in Erwägung ziehen, Vielkomponentenrechnungen für die erweiterte Modenkopplungstheorie durchzuführen, so wie für die Mischungseffekte in dieser Arbeit. In der Simulation konnte mit den inkohärenten Korrelatoren das Zusammenspiel reiner struktureller Relaxation und der scherinduzierten Relaxation detailliert untersucht werden: Bei Korrelatoren in der Flüssigkeit konkurriert der scherinduzierte Zerfall mit dem intrinsischen Zerfall. Tief genug im Glas wird der Zerfall immer von der Scherung eingeleitet. Die Korrelatoren im schergeschmolzenen Glas fallen deshalb auf eine Masterfunktion, die auch in der Simulation für kleine Scherraten gefunden wurde. Die bekannte Taylordispersion, die sich z.B. durch ein kubisches Anwachsen im mittleren Verschiebungsquadrat in Scherrichtung bemerkbar macht, konnte in der Simulation nachgewiesen werden. Die Voraussage, der erweiterten Modenkopplungstheorie, dass die Scher-und Gradientenrichtung über den selben Diffusionskoeffizienten verknüpft sind, konnte bestätigt werden. Außerdem wurde das von der Modenkopplungstheorie vorhergesagte Skalenverhalten, das besagt, dass die Diffusionskoeffizienten für kleine Scherraten linear mit der Scherrate skalieren, herausgearbeitet. Dennoch hinterlässt die Taylordispersion ein letztes Rätsel, da die Ordnung nach der führenden Ordnung in der Simulation und der Theorie unterschiedlich sind: Ein quadratischer Term, der in den Modenkopplungsergebnissen vorkommt, konnte in der Simulation nicht gefunden werden. In Verbindung mit der Diskussion der mittleren Verschiebungsquadrate konnte ein Zusammenhang zwischen dem Überschwingen der Spannungen und der Superdiffusivität hergestellt werden. Das Überschwingen der Spannungen steht auch in Verbindung mit der Anisotropie der Korrelatoren in bestimmte Richtungen. Weiter wurden Kohärente Korrelatoren für den stationären Fall berechnet. Damit wird gezeigt, dass es möglich ist, diese Größen aus der Simulation zu berechnen, obwohl es eines erheblich größeren numerischen Aufwands und einiger numerischer Tricks bedarf. Die Berechnung der transienten Korrelatoren, die die Hauptgrößen in der erweiterten Modenkopplungstheorie darstellen, ist eine Aufgabe, die in Bälde gelöst werden kann, was einen Vergleich mit der Theoriemöglich machen wird. Schlussendlich wurde im letzten Kapitel mit einer Methode, die von Matthias Krüger entwickelt wurde, die Lücke zwischen stationären und transienten Korrelatoren geschlossen. Die Simulation konnte zeigen, dass für kurze Wartezeiten der Zusammenhang, der zwischen den transienten, stationären und Korrelatoren des Ruhesystems besteht, in perfektem Einklang mit der Gleichung von Krüger steht. Außerdem konnte sogar für große Wartezeiten ein gutes Ergebnis erzielt werden, obwohl hierbei der Schermodul als Fitparameter einging. Eine Erweiterung der bestehenden Simulation auf oszillatorische Scherung zeigt eine Verbin- dung zwischen der Fließspannung und der nichtlinearen Antwort auf oszillatorische Scherung bei kleinen Frequenzen auf. Dies unterstützt die theoretischen Vorhersagen eines schematischen Modenkopplungsmodells, das von Joseph Brader entwickelt wurde. Für dynamische Scherspannungen unterhalb der Fließspannung reagiert das System mit linearer Antwort, oberhalb der Fließspannung reagiert das System als dissipative Flüssigkeit und zeigt höhere harmonische Anteile in der Spannungsantwort. Eine obere Schranke für die Spannungen kann durch die Fließkurven abgeschätzt werden.
Fachgebiet (DDC)
Schlagwörter
Konferenz
Rezension
Zitieren
ISO 690
WEYSSER, Fabian, 2011. Rheology of Brownian Discs [Dissertation]. Konstanz: University of Konstanz. Ersch. in: München : Dr. Hut. ISBN 978-3-86853-953-0BibTex
@phdthesis{Weysser2011Rheol-13782, year={2011}, publisher={Ersch. in: München : Dr. Hut}, title={Rheology of Brownian Discs}, author={Weysser, Fabian}, address={Konstanz}, school={Universität Konstanz} }
RDF
<rdf:RDF xmlns:dcterms="http://purl.org/dc/terms/" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bibo="http://purl.org/ontology/bibo/" xmlns:dspace="http://digital-repositories.org/ontologies/dspace/0.1.0#" xmlns:foaf="http://xmlns.com/foaf/0.1/" xmlns:void="http://rdfs.org/ns/void#" xmlns:xsd="http://www.w3.org/2001/XMLSchema#" > <rdf:Description rdf:about="https://kops.uni-konstanz.de/server/rdf/resource/123456789/13782"> <dc:language>eng</dc:language> <dc:publisher>Ersch. in: München : Dr. Hut</dc:publisher> <dcterms:available rdf:datatype="http://www.w3.org/2001/XMLSchema#dateTime">2011-06-20T11:16:17Z</dcterms:available> <dc:creator>Weysser, Fabian</dc:creator> <dcterms:title>Rheology of Brownian Discs</dcterms:title> <dcterms:isPartOf rdf:resource="https://kops.uni-konstanz.de/server/rdf/resource/123456789/41"/> <bibo:uri rdf:resource="http://kops.uni-konstanz.de/handle/123456789/13782"/> <bibo:issn>978-3-86853-953-0</bibo:issn> <dc:date rdf:datatype="http://www.w3.org/2001/XMLSchema#dateTime">2011-06-20T11:16:17Z</dc:date> <void:sparqlEndpoint rdf:resource="http://localhost/fuseki/dspace/sparql"/> <dc:contributor>Weysser, Fabian</dc:contributor> <dspace:hasBitstream rdf:resource="https://kops.uni-konstanz.de/bitstream/123456789/13782/2/Dissertation_Weysser.pdf"/> <dcterms:issued>2011</dcterms:issued> <dspace:isPartOfCollection rdf:resource="https://kops.uni-konstanz.de/server/rdf/resource/123456789/41"/> <dcterms:hasPart rdf:resource="https://kops.uni-konstanz.de/bitstream/123456789/13782/2/Dissertation_Weysser.pdf"/> <dcterms:rights rdf:resource="https://rightsstatements.org/page/InC/1.0/"/> <dcterms:abstract xml:lang="eng">In this thesis the properties of binary mixtures of hard discs, undergoing Brownian motion, have been studied. Two major cases have been considered: First the approach of the glass transition for a binary mixture along with the effect of changes in the composition within that mixture. For one selected mixture a detailed discussion of its transition followed. Second the rheological properties of this selected mixture were investigated by focusing on the distorted microstructure and the dynamical correlation functions of mainly tagged quantities. Finally time dependent shear flows were investigated. Performing simulations for different binary mixtures at the glass transition was motivated by recent results of MCT for two dimensions. With the simulations it was possible to confirm four mixing effects, predicted by MCT for these binary mixtures, by David Hajnal. It could be shown that for big radius ratios when we keep the total packing fraction constant, the increase in the concentration of the small particles leads to a melting of the glass, the plasticization. For small radius ratios the increase in the concentration of small particles at constant volume fraction leads to an even stronger glass. Both remaining effects, the increase of non-ergodicity parameters, and the slowing down of the relaxation towards the non-ergodicity parameters on increasing the small particles concentrations could be confirmed. These qualitative effects, especially the first two ones, are not only of theoretical interest, as industrial applications of so called plasticizers show. Going into detail for a selected mixture made it possible to determine the ideal MCT glass transition point. Along with the verification of the factorization close to the critical point, the α-scaling and even the increase of the plateaus according to a square root law could be validated. The simulation shows additional decay processes, leading to a decay of the correlation functions even above the ideal glass transition. This is not in accordance with MCT, however, the rise of the plateau values can be explained. This analysis closes a gap in the field of the glass transition as: First, a quantitative test for exactly the model system MCT uses was performed, and second to the knowledge of the author, no such analysis for a this two dimensional glass former has been performed before. With the characterization of the considered system the foundation for the rheological part was laid. That the simulation algorithm yields Brownian motion for the quiescent system was known before from the work of Erik Lange. This thesis goes one step further and shows, with the same arguments, that even for the sheared case the algorithm is still in accordance with short time expansions of the Smoluchowski equation for the shear modulus. The significance here is, that even on the two particle level used in the theoretical derivation, the simulation agrees in its short time asymptote with the theoretical findings. Having ensured that the simulation actually solves the Smoluchowski equation for Brownian motion under shear, and having characterized the glass transition for one special mixture in the preceding chapters, the rheological properties of exactly that mixture could be probed. In the framework of an extended MCT, MCT-ITT, which was developed by Matthias Fuchs and Michael Cates, the microscopic structure was investigated. The appearance of a yield stress, shear thinning and the distortion of the structure factor in the simulation can be qualitatively explained by MCT-ITT and the interplay of different relaxation time scales involved in the glass and in the liquid. Quantitatively we found that MCT-ITT overestimates the distortion of the structure and the resulting shear stresses and viscosities by about a factor of 10. However, it should be remarked that a simplified monodisperse MCT-ITT calculation was used, which could explain a part of the deviation. Future works should consider using a multi-component approach, as used for the analysis of the mixture effects in this thesis. Scrutinizing the tagged particle correlation functions, the interplay between purely structural relaxation and shear induced relaxation could be observed in detail in the simulation: For correlators in the liquid, the shear induced decay competes with the intrinsic decay, whereas in the glassy system the time scale is always set by the shear field. The relaxation in the shear melted glass thus follows a master function which could be found for very low shear rates in the simulation. The famous Taylor dispersion, expressing itself, for example, through a cubic growth in the mean squared displacement in the shear direction could be found in the simulations. The MCT-ITT prediction, that actually the shear and the gradient direction are connected via one single long time diffusion coefficient, could be confirmed. Furthermore the MCT-ITT scaling properties of these coefficients, stating that they scale linearly with the shear rate in the limit of vanishing shear, could be worked out with the simulation. Yet, the Taylor dispersion still leaves a riddle, as next-to-leading order asymptotes in MCT-ITT and the simulation seem to differ: A term proportional to the square of time couldn’t be found as proposed by MCT-ITT in the simulations, which propose a next-to-leading order dependence that is linear in time. Along with the discussion of the mean squared displacement, a connection between the superdiffusivity and the overshoot in the transient stress could be pointed out. Also the stress overshoot is connected to the anisotropy in the correlation functions for certain directions. Coherent correlators for the stationary state have been shown. In principle they show that simulating this quantities, although it requires a lot more numerical tricks and computational power, is possible. Calculating the transient correlators, which are the main quantity in MCT-ITT, will be a task for the near future, and will make further tests possible. Finally in the last chapter with a method, developed by Matthias Krüger, it was possible to close the gap between stationary and transient correlation functions. The simulation showed that for small waiting times the connection given between transient, stationary and equilibrium correlator is perfectly reproduced by Krüger’s relation. Furthermore, it could be shown, that even for large waiting times reasonable results can be achieved, although the shear modulus entered as a fitting parameter. Extending the simulation to oscillatory shear a connection between the yield stress of a glassy system and its non-linear reaction upon oscillatory shear for low frequencies could be inferred. The results support the theoretical results from a schematic MCT model invented by Joseph Brader. For time dependent stresses below the yield stress the system shows the elastic behavior of a solid state body. Above the yield stress the system shows the dissipative properties of a liquid accompanied by higher harmonics in the stress response. An upper bound for the stress can be estimated from the flow curves.</dcterms:abstract> <dc:rights>terms-of-use</dc:rights> <foaf:homepage rdf:resource="http://localhost:8080/"/> </rdf:Description> </rdf:RDF>