code_saturne User's Forum

Hello,

I am trying to model the flow in a porous medium. I do not, however, try to model the details (e.g. the pores) but a macro / semi-empirical approach.

Therefore, a Darcy model (grad(p) = mu/alpha . v, as porposed in http://code-saturne.org/forum/viewtopic ... arcy#p4493 would be enough for the momentum.
However, the fluid mass present in the porous model would not be reduced according to the porosity and I am not sure as to the effect on the heat transfer in the porous medium.
Is there a possibility to override the element volume (V*poro) and the element face vectors (F*poro) and if so, would this solve the issue?
al

Hello,

Since the thread you mention (dating back to 2010), a porosity modeling has been implemented, and may be activated setting iporos = 1 in usipsu (in cs_user_parameters.f90) and defined using the usporo.f90 user subroutine (defining the appropriate head loss in the same region is probably also necessary).

This functionnality's documentation is not integrated in the user documenatation, and I do not know if it exists yet or not (I'll check), so you may need to experiment a bit...

Regards,

Yvan