La physique
Hypothèses de la simulation
Afin de bien situer le cadre de la simulation, il faut expliquer l'environnement physique de l'expérience :
- On se place dans un cube de côté d'unité unité
. - Chaque particule
est modélisée par une sphère de rayon
et de masse
qui a une position Position
et une vitesse vitesse
initiale. - On suppose que seules les trois loi de newton (Principe de inertie, Principe fondamentale de la dynamique et Principe des actions réciproques) sont valables. Les forces de frotement et d'interraction (gravité et électromagnétique) ne sont donc pas prises en compte.
Théorie de la physique sur les Chocs
On prend deux particules
et
qui entrent en collision entre deux instant
(avant la collision) et
(après la collision)

Lors d'une collision entre deux particules, des lois énergétiques doivent être vérifiées afin de satisfaire le modèle de la simulation :
-
La loi de la conservation de la quantité de mouvement
La quantité de mouvement
des deux particules est définie par
Par la loi de la conservation de la quantité de mouvement, on a :

-
La loi de la conservation de l'energie
L'energie
des deux particules est définie par
Par la loi de la conservation de l'énergie, on a :

Calcul des nouvelles vitesses après un choc
Dans ce paragraphe, nous allons calculer les nouveaux vecteurs vitesses après une collision.
Définitions et contexte
On note :
-
le vecteur position,
le vecteur vitesse avant la collision et
le vecteur vitesse après la collision de la particule 
-
le vecteur position,
le vecteur vitesse avant la collision et
le vecteur vitesse après la collision de la particule 
On considère alors le schéma suivant à l'instant précis où deux particules rentrent en collision.
On pose
, un vecteur unitaire de direction
.

Calculs théorique d'une collision entre deux particules
et 
On peut réécrire les équations énergétique
et
qui deviennent



Dans une collision, la quantité de mouvement pour chaque particule est dans la direction du vecteur
au moment du contact. La quantité de mouvement est alors donnée par

On peut alors réécrire l'équation de conservation de l'énergie

Or nous savons que par
, nous avons


En remplacant
et
dans l'équation
, on obtient

Or, comme
,
devient

Ainsi, nous connaissons maintenant la valeur de
qui est

En remplacant
dans
et
, nous pouvons donc connaître les vitesses des deux particules après la collision qui sont


Calculs théorique d'une collision entre une particule
et un plan 
Nous n'allons pas refaire tous les calculs, mais nous servir des vitesses
et
trouvées.
Dans le cas présent, nous avons une particule
de vecteur vitesse
avant la collision et le vecteur vitesse
après la collision. Pour modéliser le plan
, nous considérons qu'il a un vecteur vitesse vecteur
et une masse
. On peut aussi considérer que
, où
est le vecteur normal au plan
. Par passage à la limite sur
, on a


Nous pouvons donc connaître les vitesses de la particule
et du plan
après la collision qui sont


Recherche de l'instant du choc
Définition de l'instant t' de collision
Avant de calculer des vitesses après une collision, nous devons d'abord trouver l'instant
qui va produire une collision. Pour cela, nous recherchons le premier instant de collision
qui aura lieu avant toutes les autres collisions

avec


est l'ensemble des particules
est l'ensemble des plans (les 6 facettes du cube)
Nous allons par la suite indiquer comment on calcule
et 
Instant t' de collision entre une particule
et un plan 
Nous allons expliquer comment la fonction
calcule l'instant
de collision entre une particule
et un plan 

La position de la particule
à l'instant
de collision est 
L'équation du plan
est 
Si la particule
intersectionne le plan
en
, alors,
doit satisfaire l'équation du plan

ce qui donne

On trouve alors la valeur exacte de 

Notez que la fonction
ne calcule pas la valeur de
si
(le cas où la particule ne se dirige pas en direction du plan). On aura toujours
car les 6 facettes du cube ont leur normale orienté vers l'intérieur du cube.
Instant t' de collision entre deux particules
et 
Nous allons expliquer comment la fonction
calcule l'instant
de collision entre deux particules
et 

La position de la particule
à un instant
ultérieur est 
La position de la particule
à un instant
ultérieur est 
Pour calculer l'instant
de collision, il faut résoudre l'équation

Afin d'écrire le plus simplement possible les calculs qui vont suivre, on pose




L'équation
devient

Il faut isoler la variable
. Pour cela, on utilise un programme tel que Maple qui va résoudre l'équation. On obtient alors deux solutions de la forme

avec

La fonction
ne met pas à jour la valeur de
si
(le cas où les particules ne se choquent pas) ou si
(le cas où l'on a bien une collision, mais à un instant t prime inférieur strict à 0). Dans le cas contraire, la fonction
met a jour
et 
