Mathesis Wiki

Diese Seite ist nicht editierbar. Sie können den Quelltext sehen, jedoch nicht verändern. Kontaktieren Sie den Administrator, wenn Sie glauben, dass hier ein Fehler vorliegt.
=== Protokoll 13.03. ===
==geschafft==
Nach großen Änderungen am Programmcode über die letzten Tage und teilweise verschiedenen in Konflikt stehenden Versionen und Ansätzen mussten wir zunächst Zeit darauf verwenden, bisherige Ergebnisse zu sortieren.
Im Anschluss haben wir uns mit Abschnitten aus Perrins "Brownsche Bewegung und die wahre Existenz der Moleküle" und insbesondere seiner Vorgehensweise bei der Betrachtung Brownscher Teilchen auf verscheidenen Höhenstufen innerhalb einer Probe unter dem Mikroskop beschäftigt, mit dem Ziel, daraus Ideen für unsere eigene Simulation mitzunehmen.
Bei Perrin trat das Problem auf, dass er aufgrund der optischen Einschränkungen seines (trotzdem hochpräzisen) Mikroskops nur jeweils eine kleine Horizontalschicht von weniger als zehn Mal des Teilchendurchmessers Dicke beobachten konnte, die die Teilchen unter Umständen auch schnell wieder verließen.
Angelehnt daran wollten wir unser Programm so umschreiben, dass wir überprüfen können, ob sich das Brownsche Teilchen während 50 Schritten (die Schrittzahl, die Perrin und Chaudesaiges in der berüchtigten Zeichnung zu Papier brachten) in einer bestimmten Schicht befindet. Unseren bisherigen Erfahrungen nach schien diese Zahl allerdings sehr hoch zu sein.
Wir haben dafür zunächst in der Umgebung des bestehenden Programms für Kollisionen in 3 Dimensionen(((Die matplotlib-Version, da sich VPython, worauf wir eigentlich umsteigen wollten, dafür sehr wenig eignet. In 2 Dimensionen würde sich vermutlich pygame am besten eignen, allerdings fehlten in der pygame-Realisierung, die wir anfangs erstellt hatten, inzwischen einige wichtige Elemente)) eine Projektion auf die XZ-Ebene (also eine Seitenansicht) geplottet, die Horizontalschicht in Form von einer Ober- und einer Untergrenze darauf angelegt und gezählt, wie oft das Teilchen diese verlässt.
==Offene Fragen, Probleme:==
Das Ergebnis der vorangegangenen Betrachtungen ist eng verknüpft mit der zweiten Hauptbemühung des Tages, der Anpassung der Quasi-Platzhalterwerte in der Simulation an tatsächliche physikalische Werte. Wir wollten zumindest annähernde oder aber immerhin begründet abweichende tatsächliche Teilchenradien, Massen und Geschwindigkeiten, Werte für Gravitation und Behältergröße etc. einführen, auch um uns an einem späteren Zeitpunkt evtl. die Arbeit mit weiteren relevanten Größen wie Temperatur oder z.B. der Boltzmann-Konstante zu erleichtern oder überhaupt zu ermöglichen. Allerdings konkretisierten sich dabei unsere Befürchtungen, dass es uns, da wir durch die Beschränkung der Rechenleistung und damit der drastischen Untertreibung der Teilchenzahl, nahezu unmöglich sein würde, stets "realistische" Werte zu nutzen. Aber auch die Umstellung auf höhere Geschwindigkeiten, längere Zeitintervalle (denn um in einer "kleinen" Box für einen außenstehenden Beobachter sichtbar zu sein, müssen sich Teilchen entweder sehr langsam bewegen, oder aber es müssen sehr kurze Zeitintervalle betrachtet werden- wir tun bislan defacto letzteres) und kleinere Massen führte uns nur zu fehlerhaften Programmen oder, unter der Maßgabe die Eigenschaften der Simulation in irgendeiner Weise beibehalten zu wollen, für uns nur sehr schwer zu entwirrenden Ketten aus Interdependenzen der verschiedenen Größen. Letztendlich haben wir uns dafür entschieden,nur Verhältnisse zu ändern, die detaillierten Anpassungen ans Ende zu verschieben und uns zunächst weiteren Punkten zu widmen. Nichtsdestotrotz sind weiterhin offene Probleme:
  * insgesamt Simulation zu weniger Teilchen, fallen daher unter Einfluss der Gravitation zu schnell aus der Horizontalschicht
  * Brownsches Teilchen wird zu klein dargestellt, //siehe auch vorherige Protokolle//: Die Berechnung des markersize-Attributs, das wir zur Darstellung des Radius der Teilchen verwenden, ist wider erwarten sehr kompliziert. matplotlib berechnet laut Dokumentation die Größe s mit ''s = markersize**2''. Die Einheit dieser Größe sind allerdings Punkte, nicht Pixel, und zudem wird mit der Fenstergröße skaliert. Im Zweidimensionalen konnten wir das nach langem Suchen durch ''(ax2.get_window_extent().width*(BROWNSCHESTEILCHEN_RADIUS/WIDTH)) ** 2)* 72./fig1.dpi'' beheben, im 3D-Plot funktioniert die Methode leider nicht.
Mathesis Wiki

Benutzer-Werkzeuge

Webseiten-Werkzeuge

Seiten-Werkzeuge
