Theory behind porousSimpleFoam - where to find citations?
 

Hi,

I was wondering if anyone has any references explaining the theory behind porousSimpleFoam? I can see from the source code and this tutorial ( http://www.tfd.chalmers.se/~hani/kur...ukurReport.pdf ) that there is a Darcy term added to the momentum equation. However, I can't find any papers or books to cite that explain why this is a model for porous zones. The papers I have found using porousSimpleFoam don't give a detailed derivation / explanation of the model, or any citations.

Sorry if this thread is out of place here - I couldn't find a better place for theory questions.

did u find any citation or book?
Thanks

The momentum extraction in the openFOAM implementation seems to be just Darcy-Forchheimer. The coefs are in Tensor form though, so search the openFoam user manual for the tensor syntax on that. It seems that these days you can do the same thing with FV options.

Thanks, but I do not see the equations used by porous SimpleFoam.
how you go from

U_Darcy = viscous terms + inertial terms

to

N-S = Source (inertial + viscous)

For me it is not clear at all.

Yes, I saw, and I have tried both approaches, but I have some doubts:

a) Why there is a porousSimpleFoam solver if the porous region can be implemented with fvOptions using any solver?

b) How I should modify codeAddSup in order to allow D and F variable coefficients? in particular, i do not fully understand the coordinateSystem set with regards to the source term: https://openfoamwiki.net/index.php/DarcyForchheimer

thanks in advance

Book by de Lemos is a general reference.

Cyprien Soulaine has done a lot of OpenFoam specific stuff.

https://link.springer.com/article/10...242-016-0701-3

That paper looks interesting. With regards to the questions:


a) I am not sure. My guess is that the solvers existed before the fvoptions framework was implemented, and they have stuck around so industry people can use them without doing extra work.



b) I would like to hear the answer to this question as well! I think that #codestream should work for any dictionary entry, so maybe that is something you can try. You might be able to modify the example for setfields to work in this situation:
https://wiki.openfoam.com/Programming1