open TELEMAC-MASCARET The mathematically superior suite of solvers

Is there a possibility to model a sluice gate within telemac (in matisse)

Hi
The answer is YES but this is not probably as simplis as you want, and it depends also of the sluice gate dimension.

One of the possibility is to manage the both part of the mesh separately and merge them in only one mesh (with two unmerged side) and in the simulation you manage the discharge exchange by a fortran program.
Another possibility is to change the bottom elevation according to the level upstream and downstream to block or let pas the flow

There is probably a lot of other possibilities with sources/sink points or element masking

Hope this helps

As you said it is not easy, actually it is quiet complicated. Can you tell me how the fortran file will look like (for exchanging the discharge), it will be very kind of you if you send me an example of such fortran file

I didn't have any fortran for this (i've made similar things long time ago...)
For a first test, I think you could create a weir between the 2 boundaries which represent the sluice gate.
Then you could adapt the CLSING.F subroutine to compute the discharge assuming the flow at the gate follow a weir law and you manage the discharge as a function of upstream and downstream level

Regards

Hello All,

In the meaning same of topic, could i change the elevation bottom of bed river follow the specified time? I would like to control my sluice follow time.
In my code.f:
if(AT.GT.43200.D0) then
zf%r(795) =2.d0 ! +2.d0
zf%r(797) =2.d0 ! +2.d0
zf%r(808) =2.d0 ! +2.d0
zf%r(818) =2.d0 ! +2.d0
ENDIF
ELSE
ENDIF

I got error report after run cas file.

Thanks for your time.

Best.

Hoa.

Hello,

I see a syntax error and it should be rather:

if(AT.GT.43200.D0) then
zf%r(795) =2.d0 ! +2.d0
zf%r(797) =2.d0 ! +2.d0
zf%r(808) =2.d0 ! +2.d0
zf%r(818) =2.d0 ! +2.d0
ENDIF

The next question is that when you do this the depth remains unchanged and so the free surface is modified also. A possible solution could be to do:

if(AT.GT.43200.D0) then
H%R(795) =H%R(795)+ZF%R(795)-2.D0
H%R(797) =H%R(797)+ZF%R(797)-2.D0
H%R(808) =H%R(808)+ZF%R(808)-2.D0
H%R(818) =H%R(818)+ZF%R(818)-2.D0
zf%r(795) =2.d0 ! +2.d0
zf%r(797) =2.d0 ! +2.d0
zf%r(808) =2.d0 ! +2.d0
zf%r(818) =2.d0 ! +2.d0
ENDIF

provided that it does not give negative depths.

With best regards,

Jean-Michel Hervouet

Hello!

I'm trying a similar approach to simulate the opening of a tidal gate, by changing bottom elevation as function of free surface gradient between the two sides of the gate.

I modified the SUBROUTINE UTIMP to program the bottom change as suggested by Jean-Michel Hervouet in the post #17544.

The problem is that the bottom is changing, but the flux not agree with the new bottom. I think that Telemac is not updating the bottom for flux calculation.

I attached my files for evaluation.





Thanks in advance.

With best regards,

Taoan

Hello!

I found a solution for my problem. The trick is to activate the coupling with SISYPHE in my steering file.

Of course I needed to create a SISYPHE steering file too. I my case, I set up an equilibrium inflow discharge and a big diameter for my bed sediments to prevent some erosion or deposition in SISYPHE.

I believe this is not the best solution... maybe is there some change that can I do in telemac3d source file to activate morphodynamic coupling independent of SISYPHE?

Thanks in advance.

With best regards,

Taoan