open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello.
I am simulating a 2d flow case with a rectangular domain of constant bathy (-20m).

The flow enters in the lhs of domain and leaves at rhs of the domain.

I want to keep a strip at the mid of the domain so that the water passing through the strip shoul be restricted 50%.

I tried to simulate the flow using CORPOR subroutine with defined 0.5 POROSITY for the elements that contains the strip


There seems to be no flow velocity suppression due to this POROSITY. I need anyones suggestions to carry out the simulation with the velocity suppression along the strip.

Thanks
Vks

Hi,

As explained in the user manual, porosity correction is only used for tidal flats treatment, on half-dry elements with the keyword OPTION FOR THE TREATMENT OF TIDAL FLATS = 3.
As you have 20m depth in your domain, I guess there is no dry elements and then this option is never activated.
I don't really know what you're trying to represent with this strip at the mid of the domain, but I think you should use another methodology, by increasing the bottom friction for instance or something else...

Regards,
Laurent

Hi Leballeur,

I have a similar situation. I am trying to simulate a wall in the domain with 0.5 porosity. There are two parameters I think I need to work with, porosity and dragforce. I first tried poros fortran subroutine. But later I saw there is another fortran file for dragforce. Now my question is if the dragforce fortran works for both porosity and dragforce changes in simulation. I attached the drag force fortran file here. Can you see and let me know if this file works for porosity changes as well? What if I only use poros subroutine and don't use drag force subroutine, will it represent a wall in the domain?

Hello,

What do you really mean when you says you want to "simulate a wall". Can you be clearer?

Chi-Tuan

Hi, I am also trying to change porosity for the cells of a training dam. Could you figure out how to do it in telemac? Can I define the porosity for the elements in Bluekenue? Or I need to do it in case file?

Hello,

As suggested by user Leballeur, what about reading the TELEMAC-2D user manual + TELEMAC-2D validation manual, in particular porosity and dragforce examples, where you can find an example of Fortran file, thanks to what you can define porosity (where and how).

If it is a Blue Kenue question, it should be written in the dedicated section.

Hope this helps,

Chi-Tuan

Hello Pham,

Thanks for your reply, I found the fortran file, where the example shows a porosity change in a polygon defining the polygon points. If I need to change porosity of the cells of multiple polygon, then what I need to do? What things in the fortran file I need to change to change porosity of three polygons?

Hello,

In user_corpor user subroutine, you can see that 1 polygon of NSOM vertices is defined with XSOM.

If you want 3 polygons, I would suggest to define 3 pairs of arrays XSOM1, YSOM1, XSOM2, YSOM2, XSOM3, YSOM3 as done before (with potentially a different number of vertices NSOM1, NSOM2 and NSOM3).

Then, rather than previous test:I would write something likewhich means you check if points with coordinates XX1, YY1 is in polygon 1, or 2 or 3.

Hope this helps,

Chi-Tuan

Hi Pham,

Thanks for your answer and it helped me. I am trying to simulate a training wall changing the porosity of the cells. Now my question is- if I need to use another Fortran routine to represent the hydraulic behavior. For example: drag coefficient for the polygons?

Hello,

I am not sure to understand your question.
If you want the hydraulic behavior you want to model can be totally steered by keywords in the steering file, no need for extra Fortran subroutines.
Otherwise, in particular if a specific treatment is required, yes you need to implement additional Fortran subroutines, e.g. for drag coefficient.
I think you can have a look at the dragforce example of the TELEMAC-2D module ($HOMETEL/examples/telemac2d/dragforce).

Please read the TELEMAC-2D user manual to better know what is already available:
wiki.opentelemac.org/doku.php?id=documentation_latest

Hope this helps,

Chi-Tuan