[EPIC] permeability models for PorousSolid #3918
Description
Hello,
My name is Ibrahima Djiba, a research engineer currently working with Olga Dufour’s team (https://github.com/ol-n-ga).
I am writing to you regarding the implementation of the Kozeny-Carman model for permeability calculation in GEOS.
I have tested the following examples:
staircase_3d.xml (CompressibleSolid)
PoroElastic_CarmanKozenyPermeability_base.xml and PoroElastic_CarmanKozenyPermeability_fim_smoke.xml (PorousSolid)
Then, I visualized both models in ParaView. For the first example, permeability is correctly computed and displayed. However, for the second example, no permeability variation is observed.
I compared the source files CompressibleSolid.cpp and CompressibleSolid.hpp with PorousSolid.cpp and PorousSolid.hpp and their dependencies. I noticed that the CompressibleSolid files (used in staircase_3d.xml) include routines that define variations in porosity and permeability, whereas the PorousSolid files only seem to handle porosity changes without updating permeability.
Additionally, I extracted porosity values from the VTK files and then calculated and plotted permeability values using the Kozeny-Carman model in Python. In this case, the permeability variation is clearly visible.
Could you please confirm whether the routines for updating permeability variations are indeed missing in the PorousSolid model?
Best regards,
Ibrahima Djiba