open TELEMAC-MASCARET The mathematically superior suite of solvers

Good morning!

We are trying to modify k-epsilon model in Telemac-3D. By tracing the variable ITURBV, we have identified following subroutines of k-epsilon model:

1. CSTKEP --- Set turbulence constants
2. KEPINI --- Initialize k and epsilon
3. VISCKE --- Compute viscosity
4. KEPICL --- Initialize the boundary conditions for k-epsilon model
5. KEPCL3 --- Compute KBOR, EBOR and AUBOR for the k-epsilon model
6. SOUKEP --- Prepare source terms in the diffusion equation of k and epsilon

In order to do the modification, eddy viscosity, shear production P and the turbulent time scale have to be changed. I found the eddy viscosity in VISCKE, shear production P in SOUKEP, but I couldn't find places to modify other terms like epsilon and turbulence time scale. I guess they have been integrated into the advection-diffusion step. Where can I find them? Thanks!

Best regards,
Qilong

Hello,

Yes exactly. Advection diffusion of k and epsilon is done like velocities and tracers, by a call to cvdf3d.f which is in telemac3d.f.

Regards,

JMH

Good morning,

In our modification of k-epsilon model, the shear velocity at each node was required. However, the variable UETCAR seems only contains NPOIN2 members instead of NPOIN3. So it only covers the bottom layer. Is there a variable which contains shear velocity or shear stress for all vertical layers?

Best regards,
Qilong

Hello,

No, this gradient of velocity with respect to z is computed in soukep.f (with a long series of other derivatves) but not kept.

Regards,

JMH

Thanks!