research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

[Test] Add tests for elemental mass/mole fraction and mixture fraction

Browse source
This commit is contained in:
Ray Speth 2015-02-17 23:24:00 +00:00
parent faf222ca57
commit 10b73473f9
2 changed files with 42 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
interfaces/cython/cantera/test/test_onedim.py
View file

@ -496,6 +496,14 @@ class TestDiffusionFlame(utilities.CanteraTest):
with self.assertRaises(KeyError): # missing 'stoich'
self.sim.strain_rate('stoichiometric', fuel='H2', oxidizer='H2O2')
def test_mixture_fraction(self):
self.create_sim(p=ct.one_atm)
Z = self.sim.mixture_fraction('H')
self.assertNear(Z[0], 1.0)
self.assertNear(Z[-1], 0.0)
self.assertTrue(all(Z >= 0))
self.assertTrue(all(Z <= 1.0))
class TestCounterflowPremixedFlame(utilities.CanteraTest):
referenceFile = '../data/CounterflowPremixedFlame-h2-mix.csv'

34
interfaces/cython/cantera/test/test_thermo.py
View file

@ -40,6 +40,40 @@ class TestThermoPhase(utilities.CanteraTest):
with self.assertRaises(ValueError):
self.phase.n_atoms('H', 'CH4')
def test_elemental_mass_fraction(self):
self.phase.Y = 'H2O:0.5, O2:0.5'
Zo = self.phase.elemental_mass_fraction('O')
Zh = self.phase.elemental_mass_fraction('H')
Zar = self.phase.elemental_mass_fraction('Ar')
mO = self.phase.element_index('O')
self.assertEqual(Zo, self.phase.elemental_mass_fraction(mO))
self.assertNear(Zo, 0.5 + 0.5 * (15.9994 / 18.01528))
self.assertNear(Zh, 0.5 * (2.01588 / 18.01528))
self.assertEqual(Zar, 0.0)
with self.assertRaises(ValueError):
self.phase.elemental_mass_fraction('C')
with self.assertRaises(ValueError):
self.phase.elemental_mass_fraction(5)
def test_elemental_mole_fraction(self):
self.phase.X = 'H2O:0.5, O2:0.5'
Zo = self.phase.elemental_mole_fraction('O')
Zh = self.phase.elemental_mole_fraction('H')
Zar = self.phase.elemental_mole_fraction('Ar')
mO = self.phase.element_index('O')
self.assertEqual(Zo, self.phase.elemental_mole_fraction(mO))
self.assertNear(Zo, 0.5/3 + 0.5)
self.assertNear(Zh, 0.5*2/3)
self.assertEqual(Zar, 0.0)
with self.assertRaises(ValueError):
self.phase.elemental_mole_fraction('C')
with self.assertRaises(ValueError):
self.phase.elemental_mole_fraction(5)
def test_weights(self):
atomic_weights = self.phase.atomic_weights
molecular_weights = self.phase.molecular_weights