[poeterish]

Hello everyone,


I have been trying to find a way to implement variable permeability in OpenFoam, but with no success.


As I understand it, the fvOptions file contains Darcy coefficients that are constant during all the simulation. However, in my case the mesh is compressed over time therefore the permeability of each cell also changes over time.


Help is much appreciated. If you need more explanation of my problem please let me know.


Thank you in advance.

[olesen]

See how far you get with prototyping of solving with a coded fvOption to make your own version.

[poeterish]

Dear Mark,


Thank you for your response.


I managed to create a copy of the DarcyForchheimer class so that I can modify it according to what I want to implement.


I successfully compiled and linked it to the interFoam solver I'm using by adding "mylibFiniteVolume.so" to the controlDict File.


However, I still have a hard time understanding how some methods beloging to the DarcyForchheimer class are called at each time step. For example, the method "Foam:orosityModels::myDarcyForchheimer::correct " computes the contribution of the Darcy source term to the diagonal and off-diagonal matrices of U, therefore it should be called at each time step to do its job. My problem is: in which file is this method called at every time step?

If I know the answer to this question, I think the rest should be quite straightforward:
- create a method that updates the D coefficients based on a model of my choosing,
- call it a each time step before Foam::..::correct.


I hope my description of this hurdle is clear enough for some suggestions.
Feel free to ask me for more information.


Thank you again.


Imad