open TELEMAC-MASCARET The mathematically superior suite of solvers

bonjour, nous sommes en train de simuler un écoulement de rivière avec charriage du fond, dont la section es relativement simple, quasiment trapezoidal. La pente du fond est de l'ordre de 1.7% donc le régime d'écoulement est évidemment supercritique.

Le nombre de Froude est 1.5 ou environ, mais de toutes manières supérieur a l'unité. Nous sommes en train de fixer les conditions de frontière telle que la vitesse et l'hauteur de l'écoulement (LIHBOR=5 , LIUBOR=6 , LIVBOR=6) mais le modèle n'atteint pas une condition numérique stable. Notre question donc vise a savoir quelle est la méthode numérique le plus pertinent pour résoudre la composante advective de la vitesse de l'écoulement (characteristics, implicit, SUPG, etc?)

Est-il toujours possible de fixer une condition supercritique a l'entrée du modelé sans perdre de stabilité numérique?

Hello,

In your case, if you set the right conditions at the entrance, the velocities should be unchanged by the advection step, so the issue is probably not in the choice of the advection scheme. Your conditions at the entrance should correspond to a balance between the slope and the friction, as if there were upstream an infinite channel, is it what you did ? Your initial conditions may also have an influence.

With best regards,

Jean-Michel Hervouet