CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Velocity jump near porous media

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   January 25, 2022, 07:15
Default Velocity jump near porous media
Join Date: Dec 2020
Posts: 41
Rep Power: 4
Pyrokrates is on a distinguished road

I wrote a program (OF version: v2006) to do some simulations with porous media. To reduce the jump and oscillations nearby the "inlet" and "outlet" of the porous media, I reduce the relaxationFactor of U (from 0.1 to around 0.0001) which also helps to correct the calculated pressure drop. Now I have some problems simulating loading processes. When Darcy coefficient inside my porous media changes locally, the velocity magnitude should change too (some regions get increased velocity other reduced...). Problem here is that with very small relaxationFactor, U can not compensate this change which is quite logical.

Do you have any suggestions how to be able to reduce/remove the velocity jump for larger relaxationFactors to make changes of U possible?

1) I thought about some manual correction (structured hex mesh):
Loop over the cell layer right before my porous media. Find all cells in thickness direction right behind the specific cell. Set U of the cell in front of the porous media to the mean value of the cells behind.

Im not sure, if this is possible and converge since it changes the result of the regular calculation of field U. Another problem is following: How can I get the cells after the specific cell in front of the porous media... For flat porous media I could loop over my filter cells and find the cells with matching x and y values (flow direction is z) but for e.g. more complex filters like pleated ones this is kind of hard to achieve.

2) A second idea would be a split of solving my equations over different time steps. I have UEqn, pEqn and the new one for the loading process XEqn. My idea would be following: Calculate Ueqn and pEqn inside a inner loop where I can use e.g. three relaxation factors: 0.1 for the velocity change, 0.01 for lower osciallation and 0.0001 for a stable result of velocity U paired with the correct pressure p. After this inner loop I can do one single calculation for updating XEqn.

So what I imagine is something like this (just simple example):
        //relaxU = 0.1
        #include "UEqn.H"
         #include "pEqn.H"


         //relaxU = 0.001
        #include "UEqn.H"
         #include "pEqn.H"

         //relaxU = 0.001
        #include "UEqn.H"
         #include "pEqn.H"

    #include "XEqn.H"


Do you have any suggestions about my ideas or maybe another plan?

Thanks in advance


Last edited by Pyrokrates; January 26, 2022 at 04:35.
Pyrokrates is offline   Reply With Quote


porous media, relaxation factor, velocity jump

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
SU2 - Porous media / porous jump model jw-89 SU2 1 January 9, 2022 17:35
FLUENT + porous media and direction angle zog FLUENT 2 August 11, 2019 10:47
inlet velocity in pipe partially filled with porous media Wi A Sa FLUENT 2 October 10, 2017 05:36
2Phase flow simulation in Porous Media Ansys Fluent Yaqub FLUENT 0 June 22, 2017 04:46
Testing the integrity of POROUS media and por jump Azman FLUENT 0 July 31, 2006 12:11

All times are GMT -4. The time now is 09:34.