 Verfasst am: 29. März 2018 19:18 Titel: |
|
| Physikalisch sollte das dann prinzipiell stimmen, falls die Massen der Teilchen zu 1 normiert sind. In dem Potential fehlen aber noch Konstanten. Zu dem Code selbst kann ich nichts sagen, da ich Python nicht kenne. Hast du das Programm mal für ein einzelnes Teilchen getestet? |
|
| Verfasst am: 29. März 2018 18:12 Titel: Festes Potential |
|
| Hallo, es soll ein festes Potentialfeld existieren. Die Teilchen beeinflussen sich nicht untereinander. LG |
|
| Verfasst am: 29. März 2018 18:07 Titel: |
|
| Sollen sich die n Teilchen in einem festen Lennard-Jones-Potential befinden oder sollen sie gemäß einem solchen Potential paarweise untereinander wechselwirken? | |
| Verfasst am: 29. März 2018 14:00 Titel: Potentielle und kinetische Energie im Lennard-Jones Potentia |
|
| Meine Frage: Hallo, ich muss numerisch die Differentialgleichung für n-Teilchen im Lennard-Jones Potential mit Hilfe von numerischen Methoden in Python lösen und die potentielle, kinetische und die Summe daraus Plotten. Gelöst hab ich die DGL schon mit Hilfe des Velocity-Verlet Verfahren (habe also ein Array mit den Positionen und Geschwindigkeiten zu allen Zeitschritten). Meine Ideen: Das Problem ist, dass die Gesamtenergie nicht konstant bleibt. for i in range(n): Vkin += 0.5*(norm(v[i])**2) Vpot += V(norm(x[i])) print(x[i]) print(" Norm " + str(norm(x[i]))) print(" Pot " + str(V(norm(x[i])))) print(" Kin " + str(0.5*(norm(v[i])**2))) Das ist der Code für die Berechnung der Energien. i iteriert durch alle Teilchen. Die Funktion wird für jeden Zeitschritt aufgerufen. Sieht jemand einen Fehler? |
|