BrineCO2 Fluidsystems uses wrong/misleading gasDensity_ calculation
brineco2.hh fluidystems uses an incorrect gasDensity calculation in the private
gasDensity_ function. It returns the unchanged CO2 gas density, while the density function implies a lot more complex density calculation similar to the
This seems to be like this from the very first commit on git, but e.g. in a
biofluidsystem.hh that was derived from the brineco2.hh a long time ago by @hommel on dumux-devel there is a more complex gasDensity calculation using Dalton's law available. (see below)
I suggest we change/fix this once most tests pass, because this might change some results for the
I am not sure, what we should do about the master (or 2.12).
Additionally the normalization of some moleFractions etc. seems a little sketchy to me, but that might be necessary. Additionally the second the brineco2.hh still uses the deprecated xl and xg naming. This should be changed, while fixing this.
static Scalar gasDensity_(Scalar T,
Scalar pH2O = xgH2O*pg; //Dalton' Law
Scalar pCO2 = pg - pH2O;
Scalar gasDensityCO2 = CO2::gasDensity(T, pCO2);
Scalar gasDensityH2O = H2O::gasDensity(T, pH2O);
Scalar gasDensity = gasDensityCO2 + gasDensityH2O;