|
[Sponsors] |
January 25, 2022, 07:15 |
Velocity jump near porous media
|
#1 |
Member
Chris
Join Date: Dec 2020
Posts: 45
Rep Power: 5 |
Hey,
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. I´m 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): Code:
while(simple.loop()) { while(i<1000) { //relaxU = 0.1 #include "UEqn.H" #include "pEqn.H" i++ } while(j<500) { //relaxU = 0.001 #include "UEqn.H" #include "pEqn.H" j++ } while(k<100) { //relaxU = 0.001 #include "UEqn.H" #include "pEqn.H" k++ } #include "XEqn.H" runTime.write(); runTime.printExecutionTime(Info); } Thanks in advance Pyro Last edited by Pyrokrates; January 26, 2022 at 04:35. |
|
Tags |
porous media, relaxation factor, velocity jump |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |