Bonsoir,
j'ai "encore" un petit problème,voici le sujet:
On se propose d’observer le comportement de N particules (N=4 pour commencer ….) en 2D, chaque particule interagissant avec chaque autre particule par le biais d’un potentiel de Lenard-Jones 6/12 à l'aide de scilab.
V(r)=4ε((σ/r)^12 – (σ/r)^6))
On rappel e que la distance à l’équilibre vaut r0 = σ* 21/6
a) Montrer qu’un adimensionnement judicieux permet sans perte de généralité de s’intéresser seulement à σ= 1 et ε= 1.
b)Que valent les position initiales des particules*?Pourquoi «*gigotent*» elles*?
c) A partir de la position initiale la plus stable, étudier les positions comme ci-dessus mais avec des vitesses initiales non nulles. Que se passe-t-il sur la distance moyenne entre les particules*?
d) Reprendre l’étude avec 1000 particules.
Voici ce que j'ai essayé de faire:
a) Je n'ai pas trop saisi ce qu'était un adimensionnement, mais voici ce que je fais
Si on injecte R=r/σ alors dans la formule on peut dire que σ/r=1/R soit σ=1
Pour ε je n'ai pas trouver la solution...
b)sur scilab les 4 particules forment un carré séparé de r=1.12 angstrom ce qui est différent de la valeur ro=21/6=1.12
On en conclut que les particules gigotent.
c)La distance moyenne entre les 4 particules à vitesse non nulles sera plus élevée mais je ne sais pas comment le justifier
d)Pour 1000 particules...les gigotements seront plus importants je pense
Merci de m'éclairer car mes réponses sont vraiment pas digne du bon physicien que je voudrais devenir!
-----