def test_adpac(self):
# Verify correct metric values for data from real scenario.
logging.debug("")
logging.debug("test_adpac")
domain = restart.read("lpc-test", logging.getLogger())
surfaces = [("zone_1", 2, 2, 0, -1, 0, -1), ("zone_2", 2, 2, 0, -1, 0, -1)]
variables = [
("area", "inch**2"),
("pressure_stagnation", "psi"),
("pressure", "psi"),
("temperature_stagnation", "degR"),
("temperature", "degR"),
("mass_flow", "lbm/s"),
("corrected_mass_flow", "lbm/s"),
]
metrics = surface_probe(domain, surfaces, variables, "mass")
logging.debug("lpc-test I face data:")
for i, (name, units) in enumerate(variables):
logging.debug(" %s = %g %s" % (name, metrics[i], units))
# These values have been verified with ADSPIN.
assert_rel_error(self, metrics[0], 830.494, 0.00001)
assert_rel_error(self, metrics[1], 8.95658, 0.00001)
assert_rel_error(self, metrics[2], 6.97191, 0.00001)
assert_rel_error(self, metrics[3], 547.784, 0.00001)
assert_rel_error(self, metrics[4], 509.909, 0.00001)
assert_rel_error(self, metrics[5], 120.092, 0.00001)
assert_rel_error(self, metrics[6], 202.573, 0.00001)
surfaces = [("zone_1", 0, -1, 2, 2, 0, -1)]
metrics = surface_probe(domain, surfaces, variables, "mass")
logging.debug("lpc-test J face data:")
for i, (name, units) in enumerate(variables):
logging.debug(" %s = %g %s" % (name, metrics[i], units))
assert_rel_error(self, metrics[0], 1089.77, 0.00001)
assert_rel_error(self, metrics[1], 7.08859, 0.00001)
assert_rel_error(self, metrics[2], 6.52045, 0.00001)
assert_rel_error(self, metrics[3], 553.772, 0.00001)
assert_rel_error(self, metrics[4], 540.711, 0.00001)
assert_rel_error(self, metrics[5], 0.421583, 0.00001)
assert_rel_error(self, metrics[6], 0.903318, 0.00001)
surfaces = [("zone_1", 0, -1, 0, -1, 2, 2)]
metrics = surface_probe(domain, surfaces, variables, "mass")
logging.debug("lpc-test K face data:")
for i, (name, units) in enumerate(variables):
logging.debug(" %s = %g %s" % (name, metrics[i], units))
assert_rel_error(self, metrics[0], 2870.64, 0.00001)
assert_rel_error(self, metrics[1], 8.69595, 0.00001)
assert_rel_error(self, metrics[2], 6.89906, 0.00001)
assert_rel_error(self, metrics[3], 546.863, 0.00001)
assert_rel_error(self, metrics[4], 511.839, 0.00001)
assert_rel_error(self, metrics[5], -156.175, 0.00001)
assert_rel_error(self, metrics[6], -271.477, 0.00001)
def test_adpac(self):
# Verify correct metric values for data from real scenario.
logging.debug('')
logging.debug('test_adpac')
domain = restart.read('lpc-test', logging.getLogger())
regions = [('zone_1', 2, 2, 0, -1, 0, -1),
('zone_2', 2, 2, 0, -1, 0, -1)]
variables = [('area', 'inch**2'),
('pressure_stagnation', 'psi'),
('pressure', 'psi'),
('temperature_stagnation', 'degR'),
('temperature', 'degR'),
('mass_flow', 'lbm/s'),
('corrected_mass_flow', 'lbm/s')]
metrics = mesh_probe(domain, regions, variables, 'mass')
logging.debug('lpc-test I face data:')
for i, (name, units) in enumerate(variables):
logging.debug(' %s = %g %s' % (name, metrics[i], units))
assert_rel_error(self, metrics[0], 830.494, 0.00001)
assert_rel_error(self, metrics[1], 8.95673, 0.00001)
assert_rel_error(self, metrics[2], 6.97191, 0.00001)
assert_rel_error(self, metrics[3], 547.784, 0.00001)
assert_rel_error(self, metrics[4], 509.909, 0.00001)
assert_rel_error(self, metrics[5], 120.092, 0.00001)
assert_rel_error(self, metrics[6], 202.570, 0.00001)
regions = [('zone_1', 0, -1, 2, 2, 0, -1)]
metrics = mesh_probe(domain, regions, variables, 'mass')
logging.debug('lpc-test J face data:')
for i, (name, units) in enumerate(variables):
logging.debug(' %s = %g %s' % (name, metrics[i], units))
assert_rel_error(self, metrics[0], 1089.77, 0.00001)
assert_rel_error(self, metrics[1], 7.09975, 0.00001)
assert_rel_error(self, metrics[2], 6.52045, 0.00001)
assert_rel_error(self, metrics[3], 554.000, 0.00001)
assert_rel_error(self, metrics[4], 540.770, 0.00001)
assert_rel_error(self, metrics[5], 0.421583, 0.00001)
assert_rel_error(self, metrics[6], 0.898204, 0.00001)
regions = [('zone_1', 0, -1, 0, -1, 2, 2)]
metrics = mesh_probe(domain, regions, variables, 'mass')
logging.debug('lpc-test K face data:')
for i, (name, units) in enumerate(variables):
logging.debug(' %s = %g %s' % (name, metrics[i], units))
assert_rel_error(self, metrics[0], 2870.64, 0.00001)
assert_rel_error(self, metrics[1], 8.70534, 0.00001)
assert_rel_error(self, metrics[2], 6.89906, 0.00001)
assert_rel_error(self, metrics[3], 546.989, 0.00001)
assert_rel_error(self, metrics[4], 511.875, 0.00001)
assert_rel_error(self, metrics[5], -156.175, 0.00001)
assert_rel_error(self, metrics[6], -270.858, 0.00001)
surface = domain.extract([(0, -1, 0, -1, 2, 2)])
surface.demote()
regions = [('zone_1', 0, -1, 0, -1)]
metrics = mesh_probe(surface, regions, variables, 'mass')
# These are different than the above metrics since the extraction
# doesn't average cell values on either side of the surface.
assert_rel_error(self, metrics[0], 2870.64, 0.00001)
assert_rel_error(self, metrics[1], 8.56340, 0.00001)
assert_rel_error(self, metrics[2], 6.92024, 0.00001)
assert_rel_error(self, metrics[3], 541.789, 0.00001)
assert_rel_error(self, metrics[4], 509.850, 0.00001)
assert_rel_error(self, metrics[5], -149.525, 0.00001)
assert_rel_error(self, metrics[6], -262.976, 0.00001)
