diff --git a/interfaces/cython/cantera/liquidvapor.py b/interfaces/cython/cantera/liquidvapor.py index 544947c7bb..96021ae683 100644 --- a/interfaces/cython/cantera/liquidvapor.py +++ b/interfaces/cython/cantera/liquidvapor.py @@ -1,11 +1,19 @@ # This file is part of Cantera. See License.txt in the top-level directory or # at http://www.cantera.org/license.txt for license and copyright information. -from . import PureFluid +from . import PureFluid, _cantera + + def Water(): - """Create a `PureFluid` object using the equation of state for water.""" - return PureFluid('liquidvapor.xml','water') + """ + Create a `PureFluid` object using the equation of state for water and the + `WaterTransport` class for viscosity and thermal conductivity.""" + class WaterWithTransport(PureFluid, _cantera.Transport): + __slots__ = () + + + return WaterWithTransport('liquidvapor.xml', 'water', transport_model='Water') def Nitrogen(): """Create a `PureFluid` object using the equation of state for nitrogen.""" diff --git a/interfaces/cython/cantera/test/test_transport.py b/interfaces/cython/cantera/test/test_transport.py index 73614ce43f..207279b02c 100644 --- a/interfaces/cython/cantera/test/test_transport.py +++ b/interfaces/cython/cantera/test/test_transport.py @@ -268,6 +268,59 @@ def test_molar_fluxes(self): # self.assertNear(sum(fluxes1) / sum(abs(fluxes1)), 0.0) +class TestWaterTransport(utilities.CanteraTest): + """ + Comparison values are taken from the NIST Chemistry WebBook. Agreement is + limited by the use of a different equation of state here (Reynolds) than + in the Webbook (IAPWS95), as well as a different set of coefficients for + the transport property model. Differences are largest in the region near + the critical point. + """ + @classmethod + def setUpClass(cls): + cls.water = ct.Water() + + def check_viscosity(self, T, P, mu, rtol): + self.water.TP = T, P + self.assertNear(self.water.viscosity, mu, rtol) + + def check_thermal_conductivity(self, T, P, k, rtol): + self.water.TP = T, P + self.assertNear(self.water.thermal_conductivity, k, rtol) + + def test_viscosity_liquid(self): + self.check_viscosity(400, 1e6, 2.1880e-4, 1e-3) + self.check_viscosity(400, 8e6, 2.2061e-4, 1e-3) + self.check_viscosity(620, 1.6e7, 6.7489e-5, 2e-3) + self.check_viscosity(620, 2.8e7, 7.5684e-5, 2e-3) + + def test_thermal_conductivity_liquid(self): + self.check_thermal_conductivity(400, 1e6, 0.68410, 1e-3) + self.check_thermal_conductivity(400, 8e6, 0.68836, 1e-3) + self.check_thermal_conductivity(620, 1.6e7, 0.45458, 2e-3) + self.check_thermal_conductivity(620, 2.8e7, 0.49705, 2e-3) + + def test_viscosity_vapor(self): + self.check_viscosity(600, 1e6, 2.1329e-5, 1e-3) + self.check_viscosity(620, 5e6, 2.1983e-5, 1e-3) + self.check_viscosity(620, 1.5e7, 2.2858e-5, 2e-3) + + def test_thermal_conductivity_vapor(self): + self.check_thermal_conductivity(600, 1e6, 0.047636, 1e-3) + self.check_thermal_conductivity(620, 5e6, 0.055781, 1e-3) + self.check_thermal_conductivity(620, 1.5e7, 0.10524, 2e-3) + + def test_viscosity_supercritical(self): + self.check_viscosity(660, 2.2e7, 2.7129e-5, 2e-3) + self.check_viscosity(660, 2.54e7, 3.8212e-5, 1e-2) + self.check_viscosity(660, 2.8e7, 5.3159e-5, 1e-2) + + def test_thermal_conductivity_supercritical(self): + self.check_thermal_conductivity(660, 2.2e7, 0.14872, 1e-2) + self.check_thermal_conductivity(660, 2.54e7, 0.35484, 2e-2) + self.check_thermal_conductivity(660, 2.8e7, 0.38479, 1e-2) + + class TestTransportData(utilities.CanteraTest): @classmethod def setUpClass(cls):