open TELEMAC-MASCARET The mathematically superior suite of solvers

Bonjour,

Je suis en train de simuler sur Télémac2D v7p1 un écoulement sur le littoral. J'utilise un modèle de transport des algues. Sur le manuel, il est indiqué que l'utilisation de ce modèle exige comme modèle de turbulence le modèle k-epsilon et que le schéma d'advection convenable pour k-epsilon est 14. Sauf que ce dernier n'est pas implémenté.

En utilisant la méthode des caractéristiques pour le transport de k-epsilon, j'ai un dépassement du maximum des itérations pour le terme diffusive en k-epsilon. J'ai essayé les consignes données dans les autres topic notamment:

- Changer le solveur de k-epsilon et prendre 7.
- Ajouter du préconditionnement.
- Augmenter le nombre d'itération pour k-epsilon.
- Réduire la précision pour k-epsilon (1e-6).

Mais en vain.

Si vous avez d'autres pistes je serai reconnaissante.

Hello,

Generally speaking the k-epsilon model with tidal flats is a difficult topic. Do you have tidal flats and do you really nead them in your application ? Anyway the method of characteristics should be OK as well (scheme 14 will be in version 7.2). Maybe posting your steering file would allow us to give you more hints.

With best regards,

Jean-Michel Hervouet

Bonjour monsieur Jean Michel,

En effet, pour mon cas j'ai activé la prise en compte de bancs découvrant.
Je vous mets en pièce jointe le fichier .cas. Il y a dessus les différents paramètres testés pour k-epsilon, pas toujours tous en même temps.

En vous remerciant,

Kawtar

Hello,

OK. The only improvement I see would be :

TRAITEMENT DES HAUTEURS NEGATIVES : 2

to have strictly positive depths, this could help k-epsilon.

Also you do not need to have :

DISCRETISATIONS EN ESPACE : 12 ; 11 just try 11 ; 11 instead.

You can add the diffusion coefficient in the variables for graphic printout, so that we see what kind of diffusion gives the k-epsilon model. You can also try to change EMAX and KMAX (maximum epsilon and k) in subroutine cosake.f, trying 100. instead of 1.D10.

If it is not already done, make sure that your model works without algae and with a constant diffusion, let's say = 1 m2/s.

With best regards,

Jean-Michel Hervouet

Rebonjour,

Pour votre questionnement monsieur Jean Michel, le modèle fonctionne bien sans ALGUE avec un modèle de turbulence à coéfficient de diffusion constante.

Il fonctionne bien aussi pour le cas d'un flotteur ponctuel.

Pour ce qui est de vos propostions, j'ai remis la discrétisation à 11;11. ( Ceci n'a pas d'influence).

Par contre en clipant K et epsilon à 100, les premières itérations convergent mais le calcul finit par diverger et cette fois ci c'est plutot le calcul hydrodynamique qui diverge.

Cordialement,

Kawtar

Hello,

OK, maybe we should add a clipping of diffusion coefficient on dry zones, if it is the problem. We would then need your model to investigate the case.

With best regards,

Jean-Michel