接入test_mql pn_node3_13固定气室支路

This commit is contained in:
huojiarong committed 2026-07-19 09:10:00 +00:00
1 parent 492562a517
commit 4242057bd8
3 files changed
+301 -13

No files matched your search

+65
View File
@@ -4624,6 +4624,15 @@ class TestMqlSystem:
return chamber_spec
raise KeyError(f"No fixed chamber segment starts at {node_alias}")
def chamber_segment_for_outlet_node(node_alias: str):
chamber_specs = (
self.discover_pneumatic_branch_topology().chamber_segment_specs
)
for chamber_spec in chamber_specs:
if chamber_spec.outlet_node_alias == node_alias:
return chamber_spec
raise KeyError(f"No fixed chamber segment ends at {node_alias}")
inlet_line = self.pnl0001_assembly.lines[spec.inlet_line_alias]
outlet_line = self.pnl0001_assembly.lines[spec.outlet_line_alias]
node_line_spec = self._pnl0003_spec_for_node_port(
@@ -4878,6 +4887,28 @@ class TestMqlSystem:
p4_port3_next_resistance_output_line = self.pnl0001_assembly.lines[
p4_port3_next_resistance_output_spec.alias
]
p4_port3_next_resistance_chamber_spec = chamber_segment_for_outlet_node(
p4_port3_next_resistance_node.alias
)
p4_port3_next_resistance_chamber_inlet_line = self.pnl0001_assembly.lines[
p4_port3_next_resistance_chamber_spec.inlet_line_alias
]
p4_port3_next_resistance_chamber = self.network.components[
p4_port3_next_resistance_chamber_spec.volume_alias
]
p4_port3_next_resistance_chamber_outlet_line = self.pnl0001_assembly.lines[
p4_port3_next_resistance_chamber_spec.outlet_line_alias
]
p4_port3_next_resistance_chamber_inlet_orifice = (
self.pneumatic_assembly.fixed_orifices[
p4_port3_next_resistance_chamber_spec.inlet_orifice_alias
]
)
p4_port3_next_resistance_chamber_outlet_orifice = (
self.pneumatic_assembly.fixed_orifices[
p4_port3_next_resistance_chamber_spec.outlet_orifice_alias
]
)
p4_port3_remote_orifice = self.pneumatic_assembly.variable_orifices[
p4_port3_remote_neighborhood.port_4.orifice_alias
]
@@ -4979,6 +5010,24 @@ class TestMqlSystem:
raise TypeError(
f"{p4_port3_next_resistance_orifice_alias} is not an AMESim pneumatic orifice"
)
if not isinstance(p4_port3_next_resistance_chamber, AmesimPneumaticVolume):
raise TypeError(
f"{p4_port3_next_resistance_chamber_spec.volume_alias} is not an AMESim pneumatic volume"
)
if not isinstance(
p4_port3_next_resistance_chamber_inlet_orifice,
AmesimPneumaticOrifice,
):
raise TypeError(
f"{p4_port3_next_resistance_chamber_spec.inlet_orifice_alias} is not an AMESim pneumatic orifice"
)
if not isinstance(
p4_port3_next_resistance_chamber_outlet_orifice,
AmesimPneumaticOrifice,
):
raise TypeError(
f"{p4_port3_next_resistance_chamber_spec.outlet_orifice_alias} is not an AMESim pneumatic orifice"
)
if not isinstance(p4_bridge_primary_chamber, AmesimPneumaticVolume):
raise TypeError(
f"{p4_bridge_primary_chamber_alias} is not an AMESim pneumatic volume"
@@ -5053,6 +5102,22 @@ class TestMqlSystem:
p4_port3_next_resistance_node_line=p4_port3_next_resistance_node_line,
p4_port3_next_resistance_orifice=p4_port3_next_resistance_orifice,
p4_port3_next_resistance_output_line=p4_port3_next_resistance_output_line,
p4_port3_next_resistance_chamber_inlet_line=(
p4_port3_next_resistance_chamber_inlet_line
),
p4_port3_next_resistance_chamber=p4_port3_next_resistance_chamber,
p4_port3_next_resistance_chamber_outlet_line=(
p4_port3_next_resistance_chamber_outlet_line
),
p4_port3_next_resistance_chamber_inlet_orifice=(
p4_port3_next_resistance_chamber_inlet_orifice
),
p4_port3_next_resistance_chamber_outlet_orifice=(
p4_port3_next_resistance_chamber_outlet_orifice
),
p4_port3_next_resistance_chamber_spec=(
p4_port3_next_resistance_chamber_spec
),
p4_port3_remote_orifice=p4_port3_remote_orifice,
p4_port3_remote_orifice_line=p4_port3_remote_orifice_line,
p4_port3_remote_orifice_node=p4_port3_remote_orifice_node,
+175 -5
View File
@@ -246,6 +246,12 @@ class TestMqlPn3P4NodeChamberSegmentComponents:
p4_port3_next_resistance_node_line: AmesimPnl0003Pipe
p4_port3_next_resistance_orifice: AmesimPneumaticOrifice
p4_port3_next_resistance_output_line: AmesimPnl0001Pipe
p4_port3_next_resistance_chamber_inlet_line: AmesimPnl0001Pipe
p4_port3_next_resistance_chamber: AmesimPneumaticVolume
p4_port3_next_resistance_chamber_outlet_line: AmesimPnl0001Pipe
p4_port3_next_resistance_chamber_inlet_orifice: AmesimPneumaticOrifice
p4_port3_next_resistance_chamber_outlet_orifice: AmesimPneumaticOrifice
p4_port3_next_resistance_chamber_spec: TestMqlPneumaticChamberSegmentSpec
p4_port3_remote_orifice: AmesimPneumaticOrifice
p4_port3_remote_orifice_line: AmesimPnl0003Pipe
p4_port3_remote_orifice_node: TestMqlPneumaticNode3
@@ -303,6 +309,9 @@ class TestMqlPn3P4NodeChamberSegmentSnapshot:
p4_port3_next_resistance_node_line_port_1: ThermodynamicProperties
p4_port3_next_resistance_node_line_port_2: ThermodynamicProperties
p4_port3_next_resistance_output_line: ThermodynamicProperties
p4_port3_next_resistance_chamber_inlet_line: ThermodynamicProperties
p4_port3_next_resistance_chamber: ThermodynamicProperties
p4_port3_next_resistance_chamber_outlet_line: ThermodynamicProperties
p4_port3_remote_orifice_line_port_1: ThermodynamicProperties
p4_port3_remote_orifice_line_port_2: ThermodynamicProperties
p4_port3_remote_resistance_node_line_port_1: ThermodynamicProperties
@@ -380,6 +389,10 @@ class TestMqlPn3P4NodeChamberSegmentSnapshot:
p4_port3_next_resistance_line_center_flow: float
p4_port3_next_resistance_to_output_line_flow: float
p4_port3_next_resistance_output_to_p4_flow: float
p4_port3_next_resistance_chamber_inlet_boundary_to_line_flow: float
p4_port3_next_resistance_chamber_inlet_line_to_chamber_flow: float
p4_port3_next_resistance_chamber_outlet_node_to_line_flow: float
p4_port3_next_resistance_chamber_outlet_line_to_chamber_flow: float
p4_port3_remote_node_to_orifice_line_flow: float
p4_port3_remote_node_to_resistance_flow: float
p4_port3_remote_resistance_node_to_line_flow: float
@@ -450,11 +463,14 @@ class TestMqlPn3P4NodeChamberSegmentClosure:
*self.components.p4_port3_next_chamber_outlet_line.get_state_vector(),
*self.components.p4_port3_next_resistance_node_line.get_state_vector(),
*self.components.p4_port3_next_resistance_output_line.get_state_vector(),
*self.components.p4_port3_next_resistance_chamber_inlet_line.get_state_vector(),
*self.components.p4_port3_next_resistance_chamber.get_state_vector(),
*self.components.p4_port3_next_resistance_chamber_outlet_line.get_state_vector(),
]
def apply_state_vector(self, values: list[float]) -> None:
if len(values) != 90:
raise ValueError("PN3/P4 node chamber segment state vector requires ninety values")
if len(values) != 96:
raise ValueError("PN3/P4 node chamber segment state vector requires ninety-six values")
self.components.inlet_line.set_state_vector(values[:2])
self.components.volume.set_state_vector(values[2:4])
self.components.outlet_line.set_state_vector(values[4:6])
@@ -493,7 +509,10 @@ class TestMqlPn3P4NodeChamberSegmentClosure:
self.components.p4_port3_next_chamber.set_state_vector(values[80:82])
self.components.p4_port3_next_chamber_outlet_line.set_state_vector(values[82:84])
self.components.p4_port3_next_resistance_node_line.set_state_vector(values[84:88])
self.components.p4_port3_next_resistance_output_line.set_state_vector(values[88:])
self.components.p4_port3_next_resistance_output_line.set_state_vector(values[88:90])
self.components.p4_port3_next_resistance_chamber_inlet_line.set_state_vector(values[90:92])
self.components.p4_port3_next_resistance_chamber.set_state_vector(values[92:94])
self.components.p4_port3_next_resistance_chamber_outlet_line.set_state_vector(values[94:])
def snapshot(
self,
@@ -581,6 +600,15 @@ class TestMqlPn3P4NodeChamberSegmentClosure:
p4_port3_next_resistance_output_line = (
self.components.p4_port3_next_resistance_output_line.properties()
)
p4_port3_next_resistance_chamber_inlet_line = (
self.components.p4_port3_next_resistance_chamber_inlet_line.properties()
)
p4_port3_next_resistance_chamber = (
self.components.p4_port3_next_resistance_chamber.properties()
)
p4_port3_next_resistance_chamber_outlet_line = (
self.components.p4_port3_next_resistance_chamber_outlet_line.properties()
)
p4_port3_remote_orifice_line_port_1 = (
self.components.p4_port3_remote_orifice_line.properties_1()
)
@@ -870,6 +898,32 @@ class TestMqlPn3P4NodeChamberSegmentClosure:
p4_port3_next_orifice_output_to_p4_flow = (
p4_port3_next_orifice_to_output_line_flow
)
p4_port3_next_resistance_chamber_inlet_boundary_to_line_flow = (
self.components.p4_port3_next_resistance_chamber_inlet_line.resistance_mass_flow(
port_1_pressure_pa=resistance_boundary.p,
port_1_temperature_k=resistance_boundary.T,
)
)
p4_port3_next_resistance_chamber_inlet_line_to_chamber_flow = (
self._line_to_chamber_flow(
line=p4_port3_next_resistance_chamber_inlet_line,
chamber=p4_port3_next_resistance_chamber,
orifice=self.components.p4_port3_next_resistance_chamber_inlet_orifice,
)
)
p4_port3_next_resistance_chamber_outlet_node_to_line_flow = (
self.components.p4_port3_next_resistance_chamber_outlet_line.resistance_mass_flow(
port_1_pressure_pa=p4_port3_next_resistance_node_line_port_1.p,
port_1_temperature_k=p4_port3_next_resistance_node_line_port_1.T,
)
)
p4_port3_next_resistance_chamber_outlet_line_to_chamber_flow = (
self._line_to_chamber_flow(
line=p4_port3_next_resistance_chamber_outlet_line,
chamber=p4_port3_next_resistance_chamber,
orifice=self.components.p4_port3_next_resistance_chamber_outlet_orifice,
)
)
p4_port3_next_resistance_node_balance = (
self.components.p4_port3_next_resistance_node.balance(
port_2_temperature_k=p4_port3_next_resistance_node_line_port_1.T,
@@ -880,8 +934,14 @@ class TestMqlPn3P4NodeChamberSegmentClosure:
connected_h=p4_port3_next_orifice_line_port_1.h,
),
port_1_mass_flow_g_s=p4_port3_next_resistance_node_to_next_flow * 1.0e3,
port_3_enthalpy_flow_w=0.0,
port_3_mass_flow_g_s=0.0,
port_3_enthalpy_flow_w=self._enthalpy_flow_from_node(
flow_kg_s=p4_port3_next_resistance_chamber_outlet_node_to_line_flow,
node_h=p4_port3_next_resistance_node_line_port_1.h,
connected_h=p4_port3_next_resistance_chamber_outlet_line.h,
),
port_3_mass_flow_g_s=(
p4_port3_next_resistance_chamber_outlet_node_to_line_flow * 1.0e3
),
)
)
p4_port3_next_resistance_node_to_line_flow = (
@@ -1218,6 +1278,13 @@ class TestMqlPn3P4NodeChamberSegmentClosure:
p4_port3_next_resistance_node_line_port_2
),
p4_port3_next_resistance_output_line=p4_port3_next_resistance_output_line,
p4_port3_next_resistance_chamber_inlet_line=(
p4_port3_next_resistance_chamber_inlet_line
),
p4_port3_next_resistance_chamber=p4_port3_next_resistance_chamber,
p4_port3_next_resistance_chamber_outlet_line=(
p4_port3_next_resistance_chamber_outlet_line
),
p4_port3_remote_orifice_line_port_1=p4_port3_remote_orifice_line_port_1,
p4_port3_remote_orifice_line_port_2=p4_port3_remote_orifice_line_port_2,
p4_port3_remote_resistance_node_line_port_1=(
@@ -1345,6 +1412,18 @@ class TestMqlPn3P4NodeChamberSegmentClosure:
p4_port3_next_resistance_output_to_p4_flow=(
p4_port3_next_resistance_output_to_p4_flow
),
p4_port3_next_resistance_chamber_inlet_boundary_to_line_flow=(
p4_port3_next_resistance_chamber_inlet_boundary_to_line_flow
),
p4_port3_next_resistance_chamber_inlet_line_to_chamber_flow=(
p4_port3_next_resistance_chamber_inlet_line_to_chamber_flow
),
p4_port3_next_resistance_chamber_outlet_node_to_line_flow=(
p4_port3_next_resistance_chamber_outlet_node_to_line_flow
),
p4_port3_next_resistance_chamber_outlet_line_to_chamber_flow=(
p4_port3_next_resistance_chamber_outlet_line_to_chamber_flow
),
p4_port3_remote_node_to_orifice_line_flow=(
p4_port3_remote_node_to_orifice_line_flow
),
@@ -1800,6 +1879,51 @@ class TestMqlPn3P4NodeChamberSegmentClosure:
p4_port3_next_resistance_output_line.port_2.m_flow = -snapshot.p4_port3_next_resistance_output_to_p4_flow
p4_port3_next_resistance_output_line.port_2.h_outflow = snapshot.p4_port3_next_resistance_output_line.h
resistance_chamber_spec = self.components.p4_port3_next_resistance_chamber_spec
p4_port3_next_resistance_chamber = self.components.p4_port3_next_resistance_chamber
p4_port3_next_resistance_chamber_inlet_port = self._port(
p4_port3_next_resistance_chamber,
resistance_chamber_spec.volume_inlet_port,
)
p4_port3_next_resistance_chamber_outlet_port = self._port(
p4_port3_next_resistance_chamber,
resistance_chamber_spec.volume_outlet_port,
)
TestMqlPnl0001PairChamberSegmentClosure._write_line_side(
self,
line=self.components.p4_port3_next_resistance_chamber_inlet_line,
line_properties=snapshot.p4_port3_next_resistance_chamber_inlet_line,
node_properties=snapshot.resistance_boundary,
node_to_line_flow=(
snapshot.p4_port3_next_resistance_chamber_inlet_boundary_to_line_flow
),
line_to_chamber_flow=(
snapshot.p4_port3_next_resistance_chamber_inlet_line_to_chamber_flow
),
orifice=self.components.p4_port3_next_resistance_chamber_inlet_orifice,
orifice_boundary_port=resistance_chamber_spec.inlet_orifice_boundary_port,
orifice_volume_port=resistance_chamber_spec.inlet_orifice_volume_port,
chamber_port=p4_port3_next_resistance_chamber_inlet_port,
chamber_properties=snapshot.p4_port3_next_resistance_chamber,
)
TestMqlPnl0001PairChamberSegmentClosure._write_line_side(
self,
line=self.components.p4_port3_next_resistance_chamber_outlet_line,
line_properties=snapshot.p4_port3_next_resistance_chamber_outlet_line,
node_properties=snapshot.p4_port3_next_resistance_node_line_port_1,
node_to_line_flow=(
snapshot.p4_port3_next_resistance_chamber_outlet_node_to_line_flow
),
line_to_chamber_flow=(
snapshot.p4_port3_next_resistance_chamber_outlet_line_to_chamber_flow
),
orifice=self.components.p4_port3_next_resistance_chamber_outlet_orifice,
orifice_boundary_port=resistance_chamber_spec.outlet_orifice_boundary_port,
orifice_volume_port=resistance_chamber_spec.outlet_orifice_volume_port,
chamber_port=p4_port3_next_resistance_chamber_outlet_port,
chamber_properties=snapshot.p4_port3_next_resistance_chamber,
)
p4_port3_remote_orifice_line = self.components.p4_port3_remote_orifice_line
p4_port3_remote_orifice_line.port_1.p = snapshot.p4_port3_remote_orifice_line_port_1.p
p4_port3_remote_orifice_line.port_1.m_flow = snapshot.p4_port3_remote_node_to_orifice_line_flow
@@ -2184,6 +2308,49 @@ class TestMqlPn3P4NodeChamberSegmentClosure:
connected_h_2=snapshot.p4_bridge_primary_line.h,
)
)
p4_port3_next_resistance_chamber_inlet_derivative = (
self.components.p4_port3_next_resistance_chamber_inlet_line.derivatives_from_connections(
port_1_m_flow=snapshot.p4_port3_next_resistance_chamber_inlet_boundary_to_line_flow,
connected_h_1=snapshot.resistance_boundary.h,
port_2_m_flow=(
-snapshot.p4_port3_next_resistance_chamber_inlet_line_to_chamber_flow
),
connected_h_2=snapshot.p4_port3_next_resistance_chamber.h,
)
)
p4_port3_next_resistance_chamber_spec = (
self.components.p4_port3_next_resistance_chamber_spec
)
p4_port3_next_resistance_chamber = (
self.components.p4_port3_next_resistance_chamber
)
p4_port3_next_resistance_fixed_chamber_derivative = (
p4_port3_next_resistance_chamber.derivatives_from_two_connections(
port_a_m_flow=self._port(p4_port3_next_resistance_chamber, "port_1").m_flow,
connected_h_a=(
snapshot.p4_port3_next_resistance_chamber_inlet_line.h
if p4_port3_next_resistance_chamber_spec.volume_inlet_port == "port_1"
else snapshot.p4_port3_next_resistance_chamber_outlet_line.h
),
port_b_m_flow=self._port(p4_port3_next_resistance_chamber, "port_2").m_flow,
connected_h_b=(
snapshot.p4_port3_next_resistance_chamber_inlet_line.h
if p4_port3_next_resistance_chamber_spec.volume_inlet_port == "port_2"
else snapshot.p4_port3_next_resistance_chamber_outlet_line.h
),
internal_h=snapshot.p4_port3_next_resistance_chamber.h,
)
)
p4_port3_next_resistance_chamber_outlet_derivative = (
self.components.p4_port3_next_resistance_chamber_outlet_line.derivatives_from_connections(
port_1_m_flow=snapshot.p4_port3_next_resistance_chamber_outlet_node_to_line_flow,
connected_h_1=snapshot.p4_port3_next_resistance_node_line_port_1.h,
port_2_m_flow=(
-snapshot.p4_port3_next_resistance_chamber_outlet_line_to_chamber_flow
),
connected_h_2=snapshot.p4_port3_next_resistance_chamber.h,
)
)
return [
*inlet_derivative.as_vector(),
*chamber_derivative.as_vector(),
@@ -2230,6 +2397,9 @@ class TestMqlPn3P4NodeChamberSegmentClosure:
*p4_port3_next_resistance_derivative_1.as_vector(),
*p4_port3_next_resistance_derivative_2.as_vector(),
*p4_port3_next_resistance_output_derivative.as_vector(),
*p4_port3_next_resistance_chamber_inlet_derivative.as_vector(),
*p4_port3_next_resistance_fixed_chamber_derivative.as_vector(),
*p4_port3_next_resistance_chamber_outlet_derivative.as_vector(),
]
+61 -8
View File
@@ -436,8 +436,28 @@ class TestMqlPn3P4NodeChamberSegmentTests(unittest.TestCase):
closure.components.p4_port3_next_resistance_output_line.name,
"pneumatic_70",
)
self.assertEqual(len(closure.initial_state_vector()), 90)
self.assertEqual(len(rhs), 90)
self.assertEqual(
closure.components.p4_port3_next_resistance_chamber_inlet_line.name,
"pneumatic_101",
)
self.assertEqual(
closure.components.p4_port3_next_resistance_chamber.name,
"pn_general_chamber_4",
)
self.assertEqual(
closure.components.p4_port3_next_resistance_chamber_outlet_line.name,
"pneumatic_102",
)
self.assertEqual(
closure.components.p4_port3_next_resistance_chamber_inlet_orifice.name,
"pn_orifice_22",
)
self.assertEqual(
closure.components.p4_port3_next_resistance_chamber_outlet_orifice.name,
"pn_orifice_23",
)
self.assertEqual(len(closure.initial_state_vector()), 96)
self.assertEqual(len(rhs), 96)
self.assertAlmostEqual(
snapshot.p4_balance.pressure_pa,
snapshot.p4_primary_line.p,
@@ -648,6 +668,29 @@ class TestMqlPn3P4NodeChamberSegmentTests(unittest.TestCase):
"port_4",
).name,
)
resistance_chamber_segment = (
self.system.discover_pneumatic_branch_topology().chamber_segment_specs[2]
)
self.assertEqual(
closure.components.p4_port3_next_resistance_chamber_spec,
resistance_chamber_segment,
)
self.assertEqual(
closure.components.p4_port3_next_resistance_chamber_inlet_line.name,
resistance_chamber_segment.inlet_line_alias,
)
self.assertEqual(
closure.components.p4_port3_next_resistance_chamber_outlet_line.name,
resistance_chamber_segment.outlet_line_alias,
)
self.assertEqual(
closure.components.p4_port3_next_resistance_chamber_inlet_orifice.name,
resistance_chamber_segment.inlet_orifice_alias,
)
self.assertEqual(
closure.components.p4_port3_next_resistance_chamber_outlet_orifice.name,
resistance_chamber_segment.outlet_orifice_alias,
)
self.assertEqual(
closure.components.p4_port3_remote_orifice.name,
port3_remote_neighborhood.port_4.orifice_alias,
@@ -812,7 +855,8 @@ class TestMqlPn3P4NodeChamberSegmentTests(unittest.TestCase):
)
self.assertAlmostEqual(
snapshot.p4_port3_next_resistance_node_to_line_flow
+ snapshot.p4_port3_next_resistance_node_to_next_flow,
+ snapshot.p4_port3_next_resistance_node_to_next_flow
+ snapshot.p4_port3_next_resistance_chamber_outlet_node_to_line_flow,
0.0,
delta=1.0e-12,
)
@@ -830,14 +874,23 @@ class TestMqlPn3P4NodeChamberSegmentTests(unittest.TestCase):
),
delta=1.0e-12,
)
mass_derivative_sum = sum(rhs[index] for index in range(0, 90, 2))
resistance_chamber_mass_sum = rhs[90] + rhs[92] + rhs[94]
self.assertAlmostEqual(
resistance_chamber_mass_sum,
(
snapshot.p4_port3_next_resistance_chamber_inlet_boundary_to_line_flow
+ snapshot.p4_port3_next_resistance_chamber_outlet_node_to_line_flow
),
delta=1.0e-12,
)
mass_derivative_sum = sum(rhs[index] for index in range(0, 96, 2))
self.assertAlmostEqual(
mass_derivative_sum,
0.0,
snapshot.p4_port3_next_resistance_chamber_inlet_boundary_to_line_flow,
delta=1.0e-12,
)
def test_simulates_ninety_state_pn3_p4_node_chamber_segment(self) -> None:
def test_simulates_ninety_six_state_pn3_p4_node_chamber_segment(self) -> None:
solution = self.system.simulate_pn3_p4_node_chamber_segment_from_spec(
self.spec,
inlet_node_pressure_pa=15.31e6,
@@ -852,8 +905,8 @@ class TestMqlPn3P4NodeChamberSegmentTests(unittest.TestCase):
self.assertTrue(solution.success)
self.assertEqual(len(solution.t), 3)
self.assertEqual(len(solution.y), 90)
self.assertEqual([len(row) for row in solution.y], [3] * 90)
self.assertEqual(len(solution.y), 96)
self.assertEqual([len(row) for row in solution.y], [3] * 96)
if __name__ == "__main__":