Stabilité numérique modèle 1D

[Victor987]

Bonjour à tous,

Je me posais la question de savoir pourquoi mon modèle 1D d'un circuit hydraulique sous AMESim en transitoire ne converge pas lorsque le débit circulant est nul. Je sais que c'est un problème de stabilité numérique du logiciel dû à la résolution des équations de Navier-Stokes et de l'algorithme de résolution. Cependant, sur un autre logiciel sur lequel je travaille, ce problème ne se présente pas. D'où peut venir l'astuce de résolution ? Mathématiquement parlant, d'où vient ce problème ?

Cordialement,
-----

[obi76]

Bonjour,

ça peut venir de l'estimation du CFL, qui lorsque la vitesse est nulle (et s'il n'y a pas de pas de temps maximum) tend vers l'infini.

Ca plante dès la première itération ?

[mAx6010]

AMESim n est pas un soft CFD, donc il ne resoud pas les equations de NS.
On parle d un soft multiphysique avec des bibliotheques hydraulique, mecanique, elec, etc....

[obi76]

Je vois un "débit circulant"...

[A voir en vidéo sur Futura]

[mAx6010]

Il y a 3 elements de la bibliotheque hydraulique qui necessitent une resolution "directe" des equations de NS en 1D.
Si je prends un systeme simple avec une source de debit, un orifice et un reservoir, alors la seule equation mise en jeu est la relation entre le debit et la perte de charge
Q(t)=alpha.A.sqrt(2.dp/rho) avec Q(t) etant une donnée (AMESim crache donc la pression en amont de l orifice)
c'est tout
C'etait juste pour mentionner qu AMESim n est pas un soft cfd

[Victor987]

Bonjour,

Merci pour vos réponses.

[obi76]

D'accord, merci.

Dans ce cas je ne sais pas.