Simulation d'un fluide avec la méthode SPH

• Objectif : Implémentation de la simulation d'un fluide représenté par des particules

• Modélisation du fluide : Les paramètres du fluide sont mis dans les fichiers de données de l'application qui se trouvent dans le répertoire gkit2light/src/master_MecaSim/exec.

• Travail à effectuer :
  • Complétez le code des différentes méthodes appelées dans ObjetSimuleSPH::Simulation(..) : méthodes se trouvant dans le fichier CalculsSPH.cpp

• Il faut récupérer les modifications de la base Git - ajout des fichiers relatifs à une simulation de fluide basée SPH (ObjetSimuleSPH.h, ObjetSimuleSPH.cpp, CalculsSPH.cpp) et modifications de certains fichiers (Viewer.cpp, Viewer-init.cpp, Param.cpp, etc.)

git pull

• Complétez la méthode suivante permettant de calculer la densité de chaque particule du fluide.

void ObjetSimuleSPH::CalculDensite(){}

• Complétez la méthode suivante permettant de calculer les interactions qu'il y a entre deux particules i et j, dans le cas où celles-ci rentrent en interaction l'une avec l'autre.
• Nous ferons directement le calcul de fij / rho_i pour la particule i.

void ObjetSimuleSPH::CalculInteraction(float visco)

• Comme nous faisons directement le calcul des forces d'interaction entre les particules i et j divisées par la densité de la particule i (soit fij / rho_i), nous devons modifier légèrement le solveur dans le cas où nous simulons un fluide par la méthode SPH. Ci-dessous la modification à apporter :

void SolveurExpl::CalculAccel_ForceGravite(...){
// On a calcule dans Force[i] : fij / rho_i
// Il ne reste qu'à ajouter le vecteur g de la gravité
A[i] = Force[i] + g;
}

• On considère que le fluide se trouve dans une boite
• Il faut donc gérer les collisions entre les différentes frontières de la boite que vous définirez

void ObjetSimuleSPH::CollisionPlan()