Discussion :

[Shadiest]

Salut à tous!

Je suis actuellement en train d'étudier l'application des équations de Navier-Stokes afin de simuler des fluides dynamiques.

J'ai eu l'occasion de lire l'article de Millie à ce propos, il m'a grandement aidé, je tiens tout d'abord à remercier l'auteur! (s'il passe par ici...)

J'ai donc suivi ce cours sur la simulation de fluide et ai probablement commis une erreur, voici mon résultat (le vecteur vitesse de départ est en rouge):

Etat Initial : dt=0.01 et nu=0.1


Premiere Iteration :


Seconde :


Troisième :

Edit : avec dt=0.01 et nu=0.1




Après, le système est figé

Comme vous pouvez le constater, le vecteur d'origine ne ralentit pas... A quelle étape y a-t-il problème ?

Merci d'avance

Lien vers mes sources (loin d'être optimisé)

[Shadiest]

Je vous jure, j'en fais des bonnes, j'ai oublié de remettre mon image à blanc entre chaque itération, du coup, il faut faire la soustraction entre deux images, je reposterai le résultat parce qu'il n'empêche que la velocité augmente trop vite...

Mais ce problème me parait maintenant plus clair !

[millie]

Je te conseille de jetter également un oeil à mes sources C++ (d'un piètre qualité...) que aux algorithmes. On est pas à l'abris d'une erreur dans ceux ci, alors que normalement, le code est juste

[Shadiest]

merci, c'est justement ce que je suis en train de faire, elle ne sont pas si mal que ça !

ça manque un peu de classe vecteur avec des surdefinitions d'opérateurs à mon goût... Enfin, tu dois avoir tes raisons (surement optimisation ?).

En tout cas, merci bien pour ce conseil !

P.S : C'est aparament sur le calcul de la viscocité que ça bloque...

[millie]

Dans le source, il y a une ligne décommenté pour le calcul de la viscosité car ça plombe à moitié l'application (ça rend 3 4 fois plus lent). Sinon, j'ai trop mal optimisé (en plus, on pourrait répartir facilement les calculs sur plusieurs threads), il faudra que je le refasse... mais il y a 9 chances sur 10 pour se planter dans un calcul et perdre des heures à débugger...

[Shadiest]

la solution est surement dans le GPU, c'est ce que je compte faire ensuite (ça tombe bien, c'est ma formation)