TELEMAC-MASCARET Forum: Combination: k-epsilon + Smagorinski (1/1)

Welcome, Guest

TOPIC: Combination: k-epsilon + Smagorinski

Combination: k-epsilon + Smagorinski 11 years 10 months ago #6891

robinson

Hello,

It seems that in TELEMAC-3D, if one uses the k-epsilon model in the vertical, horizntal turbulence model will be changed to the same automatically regadless of the real horizontal turbulence model in the steering file. I'm trying to use TELEMAC-3D to study tidal hydrodynamics in an estuary, with really large horizontal length scale. So I'm wondering if it is possible to use a combination of k-epsilon model(vertial) and Smagorinski model(horizontal)...

Thanks for any hint.

The administrator has disabled public write access.

Combination: k-epsilon + Smagorinski 11 years 10 months ago #6892

pham
OFFLINE
Administrator
Posts: 1559
Thank you received: 602

Hello Robinson,

You are right. If you read lecdon_telemac3d.f, around line 1262 (search ITURBV or ITURBH), you will find the change of turbulence model for the other direction when asking to use the k-epsilon model in only one direction.

When using the k-epsilon model, you have to solve a transport equation similar to the ones for U and V or a tracer. I wonder if it is really possible to split in different directions.

You can try to include this subroutine in your Fortran file + comment these lines to test your idea, at least informatically. Just try to do it "for fun" if you want, but I am not sure that it is really correct. Maybe other users already have an answer.

Regards,

Chi-Tuan

The administrator has disabled public write access.

The following user(s) said Thank You: robinson

Combination: k-epsilon + Smagorinski 11 years 10 months ago #6895

robinson

Hi, Chi-Tuan

Thanks for the suggestion, I think this is possible because in one of the validation cases of v6p2:'119_stratification', it uses k-epsilon in the vertical and constant diffusivity in the horizontal, maybe some subroutine has been modified but not 't3d_stratification.f' at least. I don't know which one it is ... could you please further elaborate how it is achieved in this validation case. Thanks again.

with best regard.

The administrator has disabled public write access.

Combination: k-epsilon + Smagorinski 11 years 4 months ago #9733

robinson	Hello, Could anyone be so kind to help... I know that TELEMAC-3D will only allow k-epsilon model for horizontal when the vertical turbulence model is set to be k-epsilon model. But I cannot find any literature/paper supporting this. so could anyone recommend one or two... Thank you very much. Best regards.
	The administrator has disabled public write access.