ufz / ogs5

OpenGeoSys THM/C simulator version 5
https://www.opengeosys.org/ogs-5
Other
39 stars 97 forks source link

Mass coefficient if $DENSITY model 2 used #147

Open lurpi opened 5 years ago

lurpi commented 5 years ago

Hi all, I think there is a missing condition in the if here: https://github.com/ufz/ogs5/blob/8b235b2dc7ee43b3df6ea9e29ebe95b7a0e97b70/FEM/fem_ele_std.cpp#L1581

If fluid compressibility is taken into account via $DENSITY (without $COMPRESSIBILITY) then the resulting mass coefficient lacks the "divided density" term. Best,

Luca