Discussion :

[elhouyoun]

Bonjour!
je suis doctorant en Physique, et ma thèse consiste à développer un code volume finis pour les écoulement polyphasiques non isotherme en milieu poreux. Pratiquement j'ai terminé le code mais j'ai des problèmes sur les resultats qui sont très bizzar! par exemple j’obtiens des pressions ou des températures complexes, ou qui affiche NaN avec matlab, des saturations complexes ou infinis, etc. je ne comprends pas d'où viennent ces erreurs car j'ai vérifié et revérifié tous mes programmes (je programme sous matlab), j'ai vérifié toutes mes équations, les maillages, ...etc. Tous me semble correcte, donc je ne comprends rien sur ce qui cloche. Ainsi, s'il y a des personnes qui peuvent m'aider je vous serais très reconnaissant.Je voudrai juste qu'on me dise les causes possibles qui peuvent donner ces types de resultats qui n'ont aucun sens physique. qu'est ce que je dois vérifié, quelle suggestion vous me proposez, etc. toute sorte de remarque me sera utile
Je vous remercie infiniment

[bretus]

Bonjour,

Plusieurs pistes :
- Ecrire des tests unitaires sur chaque fonction (ça fera sauter les coquilles dans le code et vous penserez peut-être à certains cas qui vous ont échappé : des divisions par zéro par exemple)
- Ecrire des journaux de calcul contenant les résultats intermédiaire (pour détecter à partir de quel moment vous avez des résultats étranges)