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)

[jeroman]

Bonjour.

Je suis désolé de dépoussiérer un aussi vieux post.

J'ai beaucoup aimé l'article de Millie ( http://humbert-florent.developpez.co...hmique/fluide/ ). Le soucis, c'est que le code source du projet n'est hélas plus disponible. Serait-il possible de le remettre en ligne ? Merci.

[FR119492]

Salut!

La résolution des équations de Navier-Stokes constitue l'un des problèmes du prix du millénaire.

Alors, si tu arrives à résoudre parfaitement ce problème, tu gagneras un million de dollars américains.
Bonne chance.
Jean-Marc Blanc

[millie]

Je suis désolé de dépoussiérer un aussi vieux post.

J'ai beaucoup aimé l'article de Millie ( http://humbert-florent.developpez.co...hmique/fluide/ ). Le soucis, c'est que le code source du projet n'est hélas plus disponible. Serait-il possible de le remettre en ligne ? Merci.

Je réponds, à peine... 10 ans après. "Vaut mieux tard que jamais", mais c'est quand même très à la bourre.

Je n'ai malheureusement plus les sources (et ça devait pas être ouf).

[f-leb]

Je récupère normalement les sources en cliquant sur le lien proposé, le lien a été réparé depuis je suppose. Il y a aussi un exécutable navier.exe qui fonctionne très bien, même 15 ans après

[Mat.M]

Grand merci Millie pour le code source.
Je vais m'en inspirer