[james_minto]

Hi everyone,

I'm using OpenFOAM 2.3.0 and would like to model flow through a porous media such as a sandstone core in which the permeability of the media varies spatially.

I think one approach (mentioned here) would be to use porousSimpleFoam but modify the DarcyForchheimer porosityModel so that the d and f coefficients are volume scalar fields that can be set with setFields (instead of being constants specified in the porosityProperties dictionary).

I am having two problems with this. Firstly the DarcyForchheimer class is located in src/finiteVolume/cfdTools/general/porosityModel and I am having trouble separating DarcyForchheimer from all the other finiteVolume classes so that it can be compiled separately in my user directory without getting a lot of dependency issues. Secondly I can create the new volume scalar fields for d and f (called d_local and f_local) following the adding temperature to icoFoam tutorial, but don't know where to start with changing references to d and f to d_local and f_local.

Another approach may be to modify the icoFoam solver so that d_local and f_local are included within the calculation of p and U in icoFoam.C. This might not be as nice a solution as using the existing porosityModel class but if it is a lot simpler then I would prefer this option as I have been using OpenFOAM for 2 years now but I'm new to modifying solvers.

And advice on how to implement either of these solutions (or any other solutions you may know) would be greatly appreciated. If anyone knows of a beginners guide/tutorial for creating your own solver than could you please send me a link to it.

If I make any progress on my own I'll post an update here.

Thanks a lot

James