Welcome, Guest
Username: Password: Remember me

TOPIC: Porosity - Not working

Porosity - Not working 2 years 4 weeks ago #41220

  • ssiddi7
  • ssiddi7's Avatar
  • 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 2 years 4 weeks ago #41221

  • c.coulet
  • c.coulet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 3722
  • Thank you received: 1031
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 2 years 2 weeks ago #41292

  • ssiddi7
  • ssiddi7's Avatar
  • 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 2 years 1 week ago #41347

  • ssiddi7
  • ssiddi7's Avatar
  • 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 2 years 1 week ago #41348

  • ssiddi7
  • ssiddi7's Avatar
  • 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 2 years 6 days ago #41363

  • pham
  • pham's Avatar
  • OFFLINE
  • Administrator
  • Posts: 1549
  • Thank you received: 598
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 1 year 5 months ago #42646

  • ssiddi7
  • ssiddi7's Avatar
  • 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?

cros3v.png


cros3ws.png


cross4v.png


cross4ws.png


Capture_2023-05-23.png
The administrator has disabled public write access.
Moderators: pham

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.