open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

Telemac-3d allows the use of a rating curve (Q-H) as a downstream boundary condition. But is there an existing way, in which a dynamic rating curve can me implemented. By, dynamic rating curve I mean, a rating curve that changes with time.

for example, amount of water passing through a dam with multiple sluice gates depends not only on the WSE difference across the structure but also on number of gates open, which can change with time. Thus the rating curve (relationship) is usually of the form

Q = f(h1, h2, n), where n is the number of gates open.

So, is there a way (without making changes to the fortran file) to implement the aforementioned rating curve in telemac-3d ?


thank you

regards

Som

Hello,

I'm afraid that it is a case that requires a bit of implementation in function Q3 of library telemac3d... we can provide it later as an example but it is unlikely that we can offer something sufficiently automatic considering the number of possible scenarios.

With best regards,

Jean-Michel Hervouet

Thanks a lot Dr. Hervouet for your reply,

I will then go ahead and do the implementation in function Q3. It would be very helpful if you know or have a fortran file that has previously done the implementation of a dynamic rating curve.


I am also listing how I think the implementation should be done, and it would be great if you can provide some hints on it

1) I will need to check how the current version of the static Q-H rating implementation takes place, and replace it with a user defined formula like Q = nA(H)^b

where A and b are constants, n is number of gates open, and H is WSE elevation.

2) I also need to figure out a way to take in the time-series of "n" as an input, similar to the time-series of the liquid boundary conditions.


thank you

regards

Som

Hello,

Your implementation will probably be simpler that reading a rating curve file (see read_fic_curve.f in library telemac2d). Your formula will depend on depth at given points, in Q3.f it will be H%R(point number) or ZF%R(point number)+H%R(point number) for the free surface, the only difficulty being in parallel, where you will have to use the global point number and look for its value in a subdolain, but this can come once everything works in scalar mode.

With best regards,

Jean-Michel Hervouet