Discussion :

[arriftou]

bon jours ,d'abord je tien a m'éxcuser je suis nouvo dans le forum
donc si vous voulez bien je vé reprendre ma requette
je veut programmer les equation de navivier stockes en 2d pour la pile pemfc ,et jé besoin de l'algorythme simpler en fortran
quelqu'un peut m'aider s'il vous plais
merci

[FR119492]

Salut!

jé besoin de l'algorythme simpler en fortran

Deux remarques:

1. Il n'y a pas d'algorithme simple pour l'intégration des équations de Navier-Stokes.
2. Les algorithmes ne sont généralement pas liés au langage. Je transfère donc cette discussion sur le forum algo/maths.

Jean-Marc Blanc

[bretus]

Bonsoir,

Je comprend pas tout à ta question. Résoudre Navier-Stokes, c'est résoudre une équation différentielle. Tu auras plusieurs méthodes possibles pour faire résoudre ça à une machine, par ex :
- Méthode des éléments finis
- Méthode des différences finies

Dans tous les cas, comme dit FR119492, peut importe Fortran où autres, tu dois trouver la méthode mathématique et elle dépend de beaucoup de paramètres :
- Il y a des conditions aux limites? Des bords de forme rectangulaire, circulaire?
- Tu as besoin de précision ou de calcul rapide?

Bref, c'est pas chose simple qui s'explique en un post et si tu as compris les principes mathématiques de résolution numérique, tu as l'algorithme... Là je t'avoue que je vois déjà pas où tu retrouves Navier Stokes dans cette pile...

Désolé de pas pouvoir t'aider d'avantage, je t'encourage à rechercher des exemples de résolutions plus simple sur le net d'abord (courant : équation de la chaleur)

Bye

[souviron34]

c'est surout que le nom de l'algorithme est Simplex

une recherche minimale sur Google donne des codes...

[arriftou]

bon jour tou le monde
je m'excuse j'ai pas préciser ,donc: pour l'algorythme c'est simpler ou simplex
pour la methode c'est volume finis ou différence finies(suivant x et y) 2 d
les équations a modéliser c'est l'équation de continuité ,l'équation de mouvement et l'équation d'énérgie ou de chaleur
toute les équation sont sans terme source
les limites sont réctangulaire .
les conditions aux limites on peut les initialisees avec des zéro si on veut l'essentiel c'est que le code tourne .
j'ai toutes les equations pour mon probléme mais je ne sais pas est ce que je peut l'envoyer sur le forum ou quoi faire ...
je vous remercie .

[FR119492]

Salut!
Ce qui reste incompréhensible, c'est pourquoi tu veux utiliser pour intégrer l'équation de Navier-Stokes un algorithme simpler ou somplex qui n'est pas fait pour ça.
Jean-Marc Blanc

[souviron34]

c'est Simplex

http://en.wikipedia.org/wiki/Simplex_method

C'est une méthode d'approximation par "raffinnement d'encadrement", une variante générlaisée du principe dichotomique..

Maintenant pour Navier-Stokes, mes lointains souvenirs ne seraient d'aucun secours, mais il me semble me souvenir qu'effectivement ça ne passe pas par ce genre de résolution.. Eléments finis ? je ne sais plus;..

Tu es le mieux placé, JM, pour répondre...

Je ne parlais que du Simplex

[FR119492]

Salut!

toute les équation sont sans terme source
les limites sont réctangulaire .
les conditions aux limites on peut les initialisees avec des zéro

Si tu as des conditions aux limites nulles et et pas de source, ton problème est homogène et n'admet que la solution triviale identiquement nulle.

Pour ce qui est de la résolutions, essayons de faire aussi simple que possible: comme ton domaine est rectangulaire, tu le discrétises par différences finies. Ensuite, je vois deux possibilités:

1. La discrétisation ayant remplacé ton équation aux dérivées partielles par un système d'équations différentielles ordinaires, tu simules le transitoire qui se produirait en partant de l'état de repos en intégrant pas à pas (RK4 ou autre) ce système; attention, il risque d'être assez raide (stiff)!
2. Tu cherches directement le régime permanent. Malheureusement, le système à résoudre n'est pas linéaire. Personnellement, j'essaierais de le résoudre avec une méthode de Gauss-Seidel modifiée: à la première itération, tu laisses tomber les termes quadratiques. Pour les itérations suivantes, tu recalcules chaque fois ces termes quadratiques avec lesquels tu corriges ton système. Si le nombre de Reynolds n'est pas trop grand, ça devrait aller assez bien.

A vrai dire, je n'ai jamais eu l'occasion (ou le temps) de le faire, mais, faute de mieux, j'essaierais la seconde méthode qui me semble assez simple à implémenter.
Jean-Marc Blanc

[Ehouarn]

Bonjour,

c'est Simplex

Non, vu qu'il est question de volumes finis et de Navier-Stockes, je pencherai plutôt pour SIMPLER (pour "SIMPLE Revisited"; SIMPLE étant l'acronyme pour "Semi-Implicit Method for Pressure-Linked Equations"), un algorithme prédicteur-correcteur pour le calcul de la pression.