Startseite
Forum
Fragen
Suchen
Formeleditor
Über Uns
Registrieren
Login
FAQ
Suchen
Foren-Übersicht
->
Elektrik
Antwort schreiben
Benutzername
(du bist
nicht
eingeloggt!)
Titel
Nachrichtentext
Smilies
Weitere Smilies ansehen
Schriftfarbe:
Standard
Dunkelrot
Rot
Orange
Braun
Gelb
Grün
Oliv
Cyan
Blau
Dunkelblau
Indigo
Violett
Weiß
Schwarz
Schriftgröße:
Schriftgröße
Winzig
Klein
Normal
Groß
Riesig
Tags schließen
Schreibt eure Formeln hier im Board am besten mit Latex!
So gehts:
Latex-Kurzbeschreibung
|
Formeleditor
[quote="capstrovor1234"][b]Meine Frage:[/b] Hallo, ich muss numerisch die Differentialgleichung für n-Teilchen im Lennard-Jones Potential [latex]V(r)=(1/r)^{12} - (1/r)^6[/latex] 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). [b]Meine Ideen:[/b] 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?[/quote]
Optionen
HTML ist
aus
BBCode
ist
an
Smilies sind
an
BBCode in diesem Beitrag deaktivieren
Smilies in diesem Beitrag deaktivieren
Spamschutz
Text aus Bild eingeben
Alle Zeiten sind GMT + 1 Stunde
Gehe zu:
Forum auswählen
Themenbereiche
----------------
Mechanik
Elektrik
Quantenphysik
Astronomie
Wärmelehre
Optik
Sonstiges
FAQ
Sonstiges
----------------
Off-Topic
Ankündigungen
Thema-Überblick
Autor
Nachricht
Huggy
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?
capstrovor
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
Huggy
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?
capstrovor1234
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?