porous medium

Questions and remarks about code_saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Post Reply
alberich

porous medium

Post by alberich »

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
Yvan Fournier
Posts: 4208
Joined: Mon Feb 20, 2012 3:25 pm

Re: porous medium

Post by Yvan Fournier »

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
Post Reply