Les conditions initiales (ou conditions aux limites) te permet de résoudre ces équations (avoir une solution très approchée) par une méthode numérique : Différences finis, éléments finis ou comme la dit juju12 les volumes finis. Si tu n'est pas habitué à coder l'une de ces méthodes il vaut mieux utiliser les différences finis. Cependant, même si elle donne de bon résultats, elle est moins optimisée que les deux autre et surtout elle peut être instable (elle diverge).
En gros la philosophie de la méthode c'est de remplacer une dérivée du/dx par une différence (u2-u1)/(x2-x1). (On voit bien l'importance des conditions aux limites par exemple, x2 et x1 sont connues c'est ton maillage pour avoir u2 il faut connaitre u1 qui est donné par les CL)
voir
ICI pour plus de détailsSinon pour les éléments et volumes finis c'est plus complex et faut bien gérer les matrices et la génération du maillage peut être très complexe tout dépend de ta géométrie.
Elements finis
ICIVolumes finis
ICIPour du code en ce qui concerne les volumes et éléments finis je n'ai rien car à l'école on utilisai la toolbox PDETool de matlab. Sinon j'ai des script matlab pour les différences finies
Pour des infos sur l'équation j'imagine que tu as déja vu cet article de
Wikipedia.Et je me répete, il vas falloir poser des hypothèses pour pouvoir résoudre l'équation
Extrait de l'article Wikipedia : La résolution de l'équation de Navier-Stokes est extrêmement difficile.
À la complexité inhérente aux équations aux dérivées partielles
s'ajoutent celle de la non-linéarité introduite par le terme d'advection
de l'accélération. La plupart du temps, on essaie de résoudre une
version simplifiée de l'équation en éliminant l'un de ces termes (En posant des hypothèses). Par
exemple, à faible nombre de Reynolds, on peut négliger le terme advectif (écoulement de Stokes) et à fort nombre de Reynolds, on s'affranchit de la viscosité (équation d'Euler).Voila
Damien