TELEMAC-MASCARET Forum :: Topic: Porosity

Porosity - Not working 1 year 7 months ago #41220

ssiddi7 OFFLINE Senior Boarder Posts: 66	yes, a new polygon where 55 points of coordinates. I hanged the fortran file based on the new polygon points. I mean I changed NSOM and XSOM, YSOM values. What else I need to modify here?
	The administrator has disabled public write access.

Porosity - Not working 1 year 7 months ago #41221

c.coulet
OFFLINE
Moderator
Posts: 3632
Thank you received: 1010

Try to write a consistent test to modify the porosity values ...

IF(INPOLY(XX1,YY1,XSOM,YSOM,NSOM)) THEN
          IF(XX1.GE.-10.D0.AND.XX1.LE.10.D0) THEN
            POROS%R(IELEM) = POR
          ELSEIF(XX1.LT.-10.D0) THEN
            POROS%R(IELEM) = POR - (1.D0-POR) * (XX1+10.D0) / 10.D0
          ELSEIF(XX1.GT.10.D0) THEN
            POROS%R(IELEM) = POR + (1.D0-POR) * (XX1-10.D0) / 10.D0
          ENDIF
        ENDIF

When element XX1, YY1 is in your polygon, the following conditions are probably not satisfied except maybe the last one ...
Then you probably have porosity = 0 in some elements...

Christophe

The administrator has disabled public write access.

The following user(s) said Thank You: ssiddi7

Porosity - Not working 1 year 7 months ago #41292

ssiddi7
OFFLINE
Senior Boarder
Posts: 66

Hi,

Thanks for your last reply. I am still getting the same error. Which is "DRAGFO: AREA OF ZONE EQUAL TO ZERO" after using one condition in the drag force Fortran subroutine. Now I am wondering if I am using the correct Fortran file even. I added my Fortran user_corpor and output file of the simulation here. I input one polygon for changing the drag coefficient. Can you please have a look and let me know what can be the issue?

The administrator has disabled public write access.

Porosity - Not working 1 year 6 months ago #41347

ssiddi7
OFFLINE
Senior Boarder
Posts: 66

Hi Pham,

Have you had any experience with drag force modification? I am struggling with this. I am trying to change the drag coefficient of the polygon, where I put all the coordinates in the drag force Fortran subroutine. I attached my Fortran file and the output error so that you can have a look. The error says DRAGFO: AREA OF ZONE EQUAL TO ZERO.

Let me know if you have any ideas.

File Attachment:

File Name: user_corpor_2022-10-22.f
File Size: 6 KB

The administrator has disabled public write access.

Porosity - Not working 1 year 6 months ago #41348

ssiddi7
OFFLINE
Senior Boarder
Posts: 66

Hi Pham,

Have you had any experience with drag force modification? I am struggling with this. I am trying to change the drag coefficient of the polygon, where I put all the coordinates in the drag force Fortran subroutine. I attached my Fortran file and the output error so that you can have a look. The error says DRAGFO: AREA OF ZONE EQUAL TO ZERO.

Let me know if you have any ideas.

File Attachment:

File Name: user_corpor_2022-10-22-2.f
File Size: 6 KB

The administrator has disabled public write access.

Porosity - Not working 1 year 6 months ago #41363

pham OFFLINE Administrator Posts: 1467 Thank you received: 563	Hello, I think Christophe gave you good tips in this topic and the other one: opentelemac.org/index.php/kunena/16-tele...-of-h?start=30#41351 If you do not succeed, upload a .zip with the every files needed to run your simulation and get your error (in particular, steering file, geometry file, boundary condition file, Fortran files etc.). Chi-Tuan
	The administrator has disabled public write access.

Porosity - Not working 11 months 3 weeks ago #42646

ssiddi7
OFFLINE
Senior Boarder
Posts: 66

Hi,
I have assigned a polygon with drag coefficient value and porosity 0.25 values to simulate a wall. Now I checked the velocity in the cross-section and also water surface elevation. It shows where the wall is, the velocity becomes low as expected, but water surface shows upward or downward spike. I cannot make sense why this is happening. I have attached two cross section in the below. Also, if I try to see the water surface map, the variation looks very abrupt. I sent a screenshot here how it looks like. Do you have any idea if the model is okay and if it makes sense?