[Pyrokrates]

Hi,


I would like to change the way, I can input porosity for my medium. I know that i can use fvOptions to define porous media and set d and f coeffs. I would like to define this locally later on due to inhomogeneous material and porosity. I know I could define multiple cellZones for each permeability to set d and f local but I need to do this with thousands of cells meaning the better way for me would be to give a field of d and f coeffs out of another program (in my case python calculates K local) instead of hundreds of porostiy zones.



Now I would like to change the input but I dont understand the whole file system.


What I know so far:
1) porousSimpleFoam uses Ueqn.H with fv.Options(U) on right hand side. Thats why I need to use #include "fvOptions.H" at porousSimpleFoam.C.
2) Now the part that confusing me starts.
2.1) fvOptions.C and fvOptions.H are also coupled with fvOption.C and fvOption.H and same with fvOptionList. But I cant finde the connection to DarcyForchheimer.H and DarcyForchheimer.C. This construct seems too confusing for me :/

2.2) To solve Ueqn and peqn I need to import the .H files with the equations. For DarcyForchheimer the equations are at the .C files... I thought .H are header files to define Variables etc. and .C are the scource files to solve equations or am I wrong with this.



If I could understand all the connections I may be able to replace all this and pass a field of K or directly d and f coeffs to the field of DarcyForchheimer equations for calculating the results.


Thanks in advance for any help.




Chris