Hi Palabos community,

I’m trying to calculate the absolute permeability on my sample following the instructions on the permeability.cpp tutorial. Everything runs well at the beginning, but the simulation never converges (I even tried to set-up my own convergence criterion).

Taking a closer look at the .vtks I noticed that there is very high velocity ‘currents’ where the outlet of the boundary condition is placed (Figure 1). I tried reducing the pressure gradient but the results look the same:

This is after 1000 iterations.

If I calculate the permeability every 250 iterations, the relative difference with the previous one is around 10% every time (I tried 1M iterations).

Any suggestions are welcome. The code (cpp) can be found here:

And the example is this one:

Thanks!!