open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I was trying to use k-epsilon model in a 3D hydrodynamics model. Because the domain of our mesh is very large, in the steering file I set

HORIZONTAL TURBULENCE MODEL = 1
VERTICAL TURBULENCE MODEL = 3
COEFFICIENT FOR HORIZONTAL DIFFUSION OF VELOCITIES = 1.0E-0
COEFFICIENT FOR VERTICAL DIFFUSION OF VELOCITIES = 1.0E-6

I think the turbulence model cannot solve k-epsilon in horizontal because the scale is too large compared with vertical. But the program always set HORIZONTAL TURBULENCE MODEL to 3 automatically while it's not coupled with Sisyphe. If I couple Telemac-3D with Sisyphe, the program will not give such message.

I'm wondering how does k-epsilon model work in Telemac?

Thanks in advance!

Hello,

In lecdon_telemac3d.f (the subroutine that reads keywords for Telemac-3D), you can see that indeed, if a k-epsilon model is chosen for horizontal or vertical, the other (vertical or horizontal) one must be the same (otherwise, it is changed).

I am surprised that you do not have such a message when coupling Telemac-3D with Sisyphe. When coupling these two codes, you have first a listing for Telemac-3D then a listing for Sisyphe. Can you read the full listing again?

By the way, for k-epsilon model in Telemac, you can first read Jean-Michel's book.

Regards,

Chi-Tuan