open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,there.

In a riverbank break case, I would like to definite the gate as an open boundary.

However, the width and height of the gate vary with time.

How can I prescribe the open boundary condition?

Should I resort to the subroutine bord.f or something else?

Appreciated could some one give me some suggestions.

Please!

Best wishes

Actually, I have no idea of how to re-program in the bord.f.

Hello

If I have well understood, you want to change a liquid boundary condition in time and in space. This is a little bit tricky!
You need to change boundary condition type for the nodes of interest (for instance from solid to liquid boundary condition) in subroutine propin_telemac2d.f. Moreover, you need to edit bord.f subroutine in order to introduce time variation.
If you are not familiar with fortran programming (and it seems to be the case), this solution is not preferable. If you give more details about your model and what you want to do with riverbank breaking, we may give you easier solutions.
Finally, bank failure can be tackled easily and properly when considering intern nodes.

with my best regards
Riadh

Thank you so much, riadh.

I'm so glad that you understood what I meant.

Yes,my original thought was reprogramming. But I found it seemed quite difficult for me at present because of time limit.

Two days earlier, I studied the validation case 019_dambreak, and I think that maybe I can change the river bank nodes into internal nodes and drive the flow by the water depth difference between the two sides of the bank as the case 019 did. Is this practice what you meant "considering intern nodes"?

Best wishes
Hong

Hello Hong

Mostly yes. It depends on which scenario you want to simulate. People use to consider dykes on the river bank as internal nodes (i.e. the domain is wide enough to consider boundary nodes far from the immediate river bank). Then, you can simulate the scenarii that you want (breaching, bottom change in time, ...) and I think that this is close to what you want to do.
The breaching process, for instance, can be monitored by considering water depth that overtops the dyke (see Breaching section in the user manual).
I hope that this helps

Riadh ATA