Auxiliary conditions (6, 6)
(1): P2 = 0.5 * (T2_P1 + T2_P2)
(1): P1 = 0.5 * (T1_P1 + T1_P2)
(1): V_h = cp * T + b * P
(1): T2_h = cp * T2 + b * P2
(1): T1_h = cp * T1 + b * P1
(1): T1_Q2 = Q1