open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear All,

I am trying to set "movable (Operational" weirs and sluice gates in a river.
In other words, I want to apply a time series of the gates' height and opening.
Do you have any idea/experience on this issue?

Many thanks in advance.

Joon

Hi

You could start from the buse.f subroutine and manage the movement of the dimension...
If the size of the weirs is large in comparison to the mesh size, it's better to envisage to work on the boundary conditions...

hope this helps

Hello,

Thanks a lot for your help!


You also talked about "working on the boundary conditions". Could you please explain the way of working on the (internal) boundary conditions?

Actually, I prefer to use TELEMAC-3D, which however does not provide a built-in weir function. So, for this TELEMAC-3D simulation, I want to connect the two boundaries of the up- and down-streams of a movable weir. However, I do not know how to treat this internal boundaries in a TELEMAC3D simulation. Do I have to make separate runs for the up- and down-streams of a weir? Or, Do you have a way to make a single run for the up- and down-streams?

Anyway, thank you again!

Joon

Hi
As you are in the topic Telemac-2D, I gave you some possible solution for Telemac2D
You could do more or less the same thing in Telemac3D but it's a little bit more complicated as you have to manage the vertical profile on the boundary...

Hope this helps

Hello Joon,

I haven't tried this but it may give you some idea.

In Telemac-2D it uses the weir data file, in which you can define the height and flow rate coefficient of a weir. I remember you used this file in your simulations. I checked the code and I found this weir data file is just read once by the subroutine LECSNG. You can find it in the Telemac source code. In that subroutine, you can find following defined variables:

NPSING%I(N) - Number of points for 1 side
NDGA1%ADR(N)%P%I(I) - Points side 1 (node number)
NDGB1%ADR(N)%P%I(I) - Points side 2 (node number)
ZDIG%ADR(N)%P%R(I) - level of the dike (weir height at each node)
PHIDIG%ADR(N)%P%R(I) - flowrate coefficients at each node

In this case, it is possible to change ZDIG%ADR(N)%P%R(I) within each time loop with a little modifications in the code.

One way is to create another piece of code that changes weir height over time and call it from the time loop in Telemac-2D; another way is to recall subroutine LECSNG from the time loop but you have to modify it and allow it to use different values for weir height.

Best regards,
Qilong

Thank you, Master.

I will try it, sooner or later.

Best Regards,
Joon