def test_cube(self): logging.debug("") logging.debug("test_cube") cube = create_cube((30, 20, 10), 5.0, 4.0, 3.0) # I face. surfaces = (("xyzzy", 2, 2, 0, -1, 0, -1),) variables = (("area", "inch**2"),) metrics = surface_probe(cube, surfaces, variables) area = metrics[0] logging.debug("area = %g (%g ft)", area, area / 144.0) self.assertEqual(area, 4.0 * 3.0 * 144.0) # J face. surfaces = (("xyzzy", 0, -1, 2, 2, 0, -1),) variables = (("area", "inch**2"),) metrics = surface_probe(cube, surfaces, variables) area = metrics[0] logging.debug("area = %g (%g ft)", area, area / 144.0) self.assertEqual(area, 5.0 * 3.0 * 144.0) # K face. surfaces = (("xyzzy", 0, -1, 0, -1, 2, 2),) variables = (("area", "inch**2"),) metrics = surface_probe(cube, surfaces, variables) area = metrics[0] logging.debug("area = %g (%g ft)", area, area / 144.0) self.assertEqual(area, 5.0 * 4.0 * 144.0)
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_wedge(self): logging.debug("") logging.debug("test_wedge") wedge = create_wedge_3d((30, 20, 100), 5.0, 0.5, 2.0, 30.0) # I face. surfaces = (("xyzzy", 2, 2, 0, -1, 0, -1),) variables = (("area", "inch**2"),) metrics = surface_probe(wedge, surfaces, variables) area = metrics[0] expected = (((pi * 2.0 ** 2.0) - (pi * 0.5 ** 2.0)) * 30.0 / 360.0) * 144.0 logging.debug("area = %g (%g ft), expected %g", area, area / 144.0, expected) assert_rel_error(self, area, expected, 0.00001) # J face. surfaces = (("xyzzy", 0, -1, 0, 0, 0, -1),) variables = (("area", "inch**2"),) metrics = surface_probe(wedge, surfaces, variables) area = metrics[0] expected = (2.0 * pi * 0.5 * 30.0 / 360.0) * 5.0 * 144.0 logging.debug("area = %g (%g ft), expected %g", area, area / 144.0, expected) assert_rel_error(self, area, expected, 0.00001) surfaces = (("xyzzy", 0, -1, -1, -1, 0, -1),) variables = (("area", "inch**2"),) metrics = surface_probe(wedge, surfaces, variables) area = metrics[0] expected = (2.0 * pi * 2.0 * 30.0 / 360.0) * 5.0 * 144.0 logging.debug("area = %g (%g ft), expected %g", area, area / 144.0, expected) assert_rel_error(self, area, expected, 0.00001) # K face. surfaces = (("xyzzy", 0, -1, 0, -1, 2, 2),) variables = (("area", "inch**2"),) metrics = surface_probe(wedge, surfaces, variables) area = metrics[0] expected = 5.0 * (2.0 - 0.5) * 144.0 logging.debug("area = %g (%g ft), expected %g", area, area / 144.0, expected) assert_rel_error(self, area, expected, 0.000001)