CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

cell volume-dependent permeability

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   April 8, 2021, 06:54
Default cell volume-dependent permeability
New Member
Join Date: Apr 2021
Posts: 2
Rep Power: 0
poeterish is on a distinguished road
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.
poeterish is offline   Reply With Quote

Old   April 18, 2021, 06:18
Senior Member
Mark Olesen
Join Date: Mar 2009
Posts: 1,503
Rep Power: 35
olesen will become famous soon enougholesen will become famous soon enough
See how far you get with prototyping of solving with a coded fvOption to make your own version.
olesen is offline   Reply With Quote

Old   April 21, 2021, 10:02
New Member
Join Date: Apr 2021
Posts: 2
Rep Power: 0
poeterish is on a distinguished road
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 "" 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.

poeterish is offline   Reply With Quote


darcy, permeability

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
[Other] mesh airfoil NACA0012 anand_30 OpenFOAM Meshing & Mesh Conversion 13 March 7, 2022 18:22
[Other] refineWallLayer Error Yuby OpenFOAM Meshing & Mesh Conversion 2 November 11, 2021 12:04
how to set periodic boundary conditions Ganesh FLUENT 15 November 18, 2020 07:09
SU2 7.0.7 Built on CentOS 7, parallel computation pyscript mpi exit error? EternalSeekerX SU2 3 October 9, 2020 19:28
Negative cell volume problem in-cylinder ziani.lotfi FLUENT 0 April 11, 2012 08:54

All times are GMT -4. The time now is 15:29.