II-2.2.1.4. Mécanique des fluides numérique 1D

Pour l'équation des fluides, on résoud le système d'équations d'EULER : Permettant de simuler un écoulement dans un conduit rectiligne mais avec une section variable |S|. Pour cette simulation de type Volumes Finis , on utilise le schéma de LAX-WENDROFF. On part d'une solution stationnaire (ou on l'atteint) obéissant au théorème de BERNOULLI valable en bas subsonique. On injecte une perturbation sur la face d'entrée au temps t=tpa < tmax ; on suit son évolution. Les conditions aux limites sont absorbantes aux deux extrémités du conduit. L'écoulement est adiabatique mais on fait l'hypothèse de l'isothermie aux extrémités avec de plus une section constante (conditions aux limites absorbantes). Les résultats sont les deux profils monodimensionnels de la vitesse et de la pression. Pour le système d'équations d'Euler (mécanique des fluides) les fichiers résultats sont : Le post-processing gère l'exploitation graphique (gnuplot) des résultats.
Modules sources en Fortran 90 : interface , géométrie du conduit , mécanique des fluides 1D , subroutines (linéaire 2), module gnuplot
Procédure IDL associée (visualisation de la géometrie du conduit) : voisgeom.pro

tests disponibles


Sommaire
page precedente home page page suivante