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 bizar! 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

[membreComplexe12]

salut,

voici des pistes de réflexion:

1°) les NaN où complexes peuvent apparaitre lorsque:
- tu fais une division par zero
- tu prends la racine d'un nombre négatif
ou d'autres choses dans ce genre...

peut etre devrais tu poser la question sur le forum matlab:
http://www.developpez.net/forums/f14...pement/matlab/

2°) essai de mettre des points d'arrêt dans matlab et executer ton programme pas par pas afin de voir à chaque étape du calcul comment évolue des grandeurs.
Ensuite, une fois que tu auras fais ceci tu arriveras à localiser la fonction qui pose problème.

3°) sinon pour vérifier ton problème je te conseil de tester des cas simple
à chaque fois. tu te fais par exemple un maillage avec un seul élément et tu vérifie le résultat.
Après d'une mais deux éléments dans ton maillage et tu recommence...Etc

ensuite tu fais des calculs sans température et tu regardes si ça fonctionne...etc

bref, tu fais des petits test et tu fais un bilan de ce qui fonctionne et de ce qui ne fonctionne pas et ensuite tu pourras mieux cibler d'où vient ton problème.

Bon courage.

ps: au fait je ne connais pas trop les méthodes de volumes finis, j'ai entendu dire qu'il fallait nécessairement des maillages admissibles pour pouvoir faire ce genre de calcul, c'est bien ça ?

c'est un peu contraignant ceci...comment générer ce type de maillage ? tu fais des cellules de voronoi ?

[Aleph69]

Bonjour,

il y a deux approches possibles :
1. prendre le temps de faire les tests unitaires qui aurait dû être conçu avant ou pendant la phase de développement de ton code de calcul;
2. se placer dans un cas particulier pour lequel une solution analytique est connue puis voir pas à pas ce qui cloche.

Bon courage.

[elhouyoun]

Merci, pour vos conseils, je vais les essayer et voir ce que ça donne. oui j'utilise des maillages admissibles, ce n'est facile oui, moi ce que je fais j'essaye de les faire dans mon bloc note d'abord, afin de voir comment élaborer le bon algorithme, je vais des petits essais et puis je fais un petit programme je fais des essais en considérant différents cas, si ça marche, j'élabore le programme.
Merci à vous et si vous avez d'autre idée ça me servira toujours