open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,
I'm trying to simulate the run-up induced by waves on the coast of a lake.
My geometry files includes the bathymetry and the topography.
I tried to use Telemac2D imposing an incident wave with a initial_wave.prn file describing the evolution of the wave crest height according to time. This wave is imposed to one of the boundaries located in the lake.

And I would like to set up free boundary conditions for all the other boundary points. When I give the number 444 in .cli file for the other points, the simulation says
THERE IS 3 LIQUID BOUNDARIES:

Then after the second boundary, the code did not print anything else, but did not crash.


So my questions are :
1- Is the approach with Telemac2d valid to simulate the impact of the wave run-up on the coast ?
2- Could you help me to parameterize my boundary conditions properly ?

Thanks a lot.
Marina

Hello,

The incident wave boundary conditions (1 1 1 in the *.cli file) is no longer in Telemac-2D, it is replaced by the Thompson's boundary conditions, please refer to this condition in the documentation. If you found that incident wave is still mentioned in the documentation please tell us and we'll cancel it. This is for entering a wave, for the run up itself I would not do it like this, I would treat the coast as a dry land that can be flooded, thus the highest point that can possibly reached by water must be in the domain and it will be dry in the beginning. Consequently the corresponding boundary of the domain can be considered as solid.

With best regards,

Jean-Michel Hervouet

Thank you Jean-Michel. It helps me a lot, but I still need some help.

If I'm not mistaken, my case requires 2 different boundary conditions:
* Thompson's BC for the boundaries that are in the lake
* a solid wall for the dry part

For the Thompson's BC, I have two different case :
* for the north-west boundary, I wish to impose the wave: can I still use a .prn file imposing the elevation according to time ?
In this case, I would indicate a 5 4 4 BC. But is it correct ?

* for the north and south boundaries, theyr are in the lake, but the wave is not entering from here. I just want to allow the flow to get out from these boudaries, without imposing any condition (transient flow?). Is it possible?
Can I put then a 4 4 4 BC?

Thanks again.
Marina

in the

Hello,

This is all correct, I understand that there are open boundaries in the lake itself... because you are not representing the whole lake ? The best would be representing the whole lake, even with large elements in areas of lesser interest. You can try 4 4 4 and Thompson, though I am not quite sure that the program will consider 3 different liquid boundaries, it may consider that there is only one big liquid boundary and one solid boundary.

With best regards,

Jean-Michel Hervouet

Hello,

This is as you said: I have two open BC within the lake because I just represent a part of the lake.

The problem with simulating all the lake is that I just know the characteristics of the waves in the area close to the coast I want to study, at the level of the actual 5 4 4 BC.

I tried to set up the open boundary conditions in the lake as 444 BCs, but I think I'm doing something wrong with the OPTION FOR LIQUID BOUNDARIES.
Because now, the simulation is running, but nothing is happening in the domain. I join you the cas file and the cli file, so that you can have a look to the BC and the keywords.

Thank you very much.
Marina

Sorry, I realized that the cli file is not the correct one (the incident wave is still indicated by 1 1 1 in the file I sent you).
When I changed the BC for a 5 4 4, the simulation stopped asking for more prescribed elevations in the liquid boundaries file. So I just added constant SL for the 2 BC in the lake in the .prn file, and restarted the simulation. It seems to be working now! I'll let you know.

Thanks a lot.
Marina