open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello!

I just want to clarify how Telemac handles porosity. From reading Hervouet (2007) the porosity terms are added into the mass conservation and momentum equations after being previously defined. From what I've read these terms defined by the user in corpor.f, but in which subroutines are the mass and momentum equations for these terms?

Appreciate any information.
Many thanks.

Hello,

There is no specific subroutines for porosity, the trick is that a mask is defined for every element and the integral of all terms computed at element level are multiplied by the mask, and the mask is the porosity. So the treatment of porosity will be triggered only by the fact that the logical MASK will be set to true and that the value of MASKEL will contain the porosity or a function of porosity.

With best regards,

Jean-Michel Hervouet