[nuttita]

Hi, I’m simulating porous media flows using porousInterFoam in OF 2.3. I created porous zone in blockMesh for flows through static porous medium and it works fine. However, I need to work with deforming medium. That means the shape of the porous zone changes. Hence, the cellzone needs to be changed during simulation. Does anyone know how to do this? Thank you.

[Santiago]

Is it the domain that changes, or you choose a different porous region everytime?

[nuttita]

Santiago, thank you for your reply.

The domain doesn't change. Only the porous region in the domain changes every time step. The change is based on a computed field variable alpha.

if alpha = 1 -->> porous region
else non-porous region

This alpha is obtained from another solver that is coupled with porousInterFoam.

[Santiago]

Is the porous region adding some source/sink of some kind?

[nuttita]

Yes, it’s the Darcy-Forchheimer sink term added to the Navier-Stokes equation. This term is active only in the porous region.

[hk318i]

Hello,


There are few ways to solve this problem but it depends on your implementation.

1. Of course you can create a cellSet each time step and apply the source term for these specific set.

• It could cause issues in parallel.
• Extra computation time to select those cells.
• An example maybe found in topSet, adaptive dynamic mesh ....

2. Don't use cellSet at all, calculate the source term as a function of alpha for each grid cell. Of you kind of if condition if it needed.

Probably there are other solution but all depend on your current code and what features you need in the future.


Best wishes,
Hassan

[nuttita]

Dear Hassan,

Thank you for your suggestion.

Nuttita