OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Merge 1b2c537e2f into 9a0ccfd853

Browse Source
This commit is contained in:
Tor Harald Sandve 2023-11-26 12:18:38 +00:00 committed by GitHub
commit 379d8af756
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 28 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

41
opm/material/components/Brine.hpp
View File

@ -30,7 +30,6 @@
#include <opm/material/components/Component.hpp>
#include <opm/material/common/MathToolbox.hpp>
namespace Opm {
/*!
@ -329,23 +328,39 @@ public:
/*!
* \copydoc H2O::liquidViscosity
*
* Equation given in:
* - Batzle & Wang (1992)
* - cited by: Bachu & Adams (2002)
* "Equations of State for basin geofluids"
* Modified Jones - Doles equation to account for temperature effects
* Coeficients from McBride-Wright, 2013
* Viscosity and density of aqueous fluids with dissolved CO2"
*/
template <class Evaluation>
static Evaluation liquidViscosity(const Evaluation& temperature, const Evaluation& /*pressure*/)
static Evaluation liquidViscosity(const Evaluation& temperature, const Evaluation& pressure)
{
Evaluation T_C = temperature - 273.15;
if(temperature <= 275.) // regularization
T_C = 275.0;
// Modified Jones - Doles equation to account for temperature effects
// Coefficients from McBride-Wright, 2013
// "Viscosity and density of aqueous fluids with dissolved CO2"
std::vector<double> c = {0.291015724, -0.119318366, 0.190911976};
const Scalar MmNaCl = 58e-3; // molar mass of NaCl [kg/mol]
// compute molality mol/kg
Scalar molality = salinity / (MmNaCl*(1 - salinity));
// compute molar concentration mol/l
Evaluation b = molality/H2O::liquidDensity(temperature, pressure, true)*1000;
Evaluation mu_corr = 1 + (c[0] + c[1] * 1000/temperature) * pow(b, 0.5) + c[2] * b;
return H2O::liquidViscosity(temperature, pressure, true) * mu_corr;
Evaluation A = (0.42*std::pow((std::pow(salinity, 0.8)-0.17), 2) + 0.045)*pow(T_C, 0.8);
Evaluation mu_brine = 0.1 + 0.333*salinity + (1.65+91.9*salinity*salinity*salinity)*exp(-A);
// old formulation (does not depend on pressure)
// Equation given in:
// - Batzle & Wang (1992)
// - cited by: Bachu & Adams (2002)
// "Equations of State for basin geofluids"
//Evaluation T_C = temperature - 273.15;
//if(temperature <= 275.) // regularization
// T_C = 275.0;
return mu_brine/1000.0; // convert to [Pa s] (todo: check if correct cP->Pa s is times 10...)
//Evaluation A = (0.42*std::pow((std::pow(salinity, 0.8)-0.17), 2) + 0.045)*pow(T_C, 0.8);
//Evaluation mu_brine = 0.1 + 0.333*salinity + (1.65+91.9*salinity*salinity*salinity)*exp(-A);
//return mu_brine/1000.0; // convert to [Pa s] (todo: check if correct cP->Pa s is times 10...)
//}
}
};