Air flow through porous media
 

Hello.


I need to simulate an air flow through a porous media, but I do not want to use the porousSimpleFoam or the rhoPorousSimpleFoam. I want to import the equations governing the problem in another solver, such as icoFoam.



I tried to build a new solver, named porousIcoFoam. I used the equation:

fvVectorMatrix UEqn
(
fvm::ddt(gamma,U)
+ fvm::div(phi, U)
- fvm::laplacian(nu, U)
- (nu)*(D&U)*lamda
);


where gamma=porosity at the porous media, and gamma=0 otherwise

lamda=1 at the porous media, and lamda=0 otherwise


In setFeilds i have defined lamda and gamma, and in transport properties the tensor D.


Is this a correct way to solve my problem?