Hi,
I want to model evaporation of brine in porous media.
after checking fugacityCoefficient of brineair.hh file,I found that we use H2O
fugacity for liquid phase.
for this reason, evaporation rate for brine and water is equal.but It is not
correct.
I found TODO for changing fugacity of water to brine,but I don't understand why
we can not add bellow code for solving this problem:
change vapor pressure of brine component to:
static Scalar vaporPressure(Scalar temperature, Scalar salinityy)
{
Scalar ps = H2O::vaporPressure(temperature); //Saturation vapor
pressure for pure water
Scalar pi = 0;
using std::log;
if (ThisType::salinity() < 0.26) // here we have hard coded the
solubility limit for NaCl
pi = (R * temperature * log(1- salinityy)); // simplified version
of Eq 2.29 in Vishal Jambhekar's Promo
else
pi = (R * temperature * log(0.74));
using std::exp;
ps *= 0.9;// Kelvin's law for reduction in saturation vapor pressure
due to osmotic potential
return ps;
}
then,change fugacitycoefficent of brineair.hh to :
--
This email was Anti Virus checked by Security Gateway.
_______________________________________________
DuMux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
