openpnm.models.phases.viscosity.water(target, temperature='pore.temperature', salinity='pore.salinity')[source]

Calculates viscosity of pure water or seawater at atmospheric pressure using Eq. (22) given by Sharqawy et. al [1]. Values at temperature higher than the normal boiling temperature are calculated at the saturation pressure.

  • target (GenericPhase) – The object for which these values are being calculated. This controls the length of the calculated array, and also provides access to other necessary thermofluid properties.

  • temperature (str) – The dictionary key containing the temperature values. Temperature must be in Kelvin for this emperical equation to work. Can be either a pore or throat array.

  • salinity (str) – The dictionary key containing the salinity values. Salinity must be expressed in g of salt per kg of solution (ppt). Can be either a pore or throat array, but must be consistent with temperature.


value – Array containing viscosity of water/seawater in [kg/m.s]

Return type



T must be in K, and S in g of salt per kg of phase, or ppt (parts per


VALIDITY: 273 < T < 453 K; 0 < S < 150 g/kg; ACCURACY: 1.5 %


[1] Sharqawy M. H., Lienhard J. H., and Zubair, S. M., Desalination and

Water Treatment, 2010.