“…The volumes, pressures, and temperatures ∆V 4 , ∆P 4 , ∆T 4 and ∆V 1 , ∆P 1 , ∆T 1 from the expansion and exhaust (displacement phase) processes 3→4 and 5→1 are computed by numerically integrating Eqs. [7][8][9][10] from the initial conditions ∆V 3 , ∆P 3 , ∆T 3 , ∆V 3 and ∆V 5 , ∆P 5 , ∆T 5 , ∆V 5 = 0 for timest 34 andt 51 , respectively such that the piston velocity ∆V 4 = ∆V 1 = 0. It is worth nothing that the timest 12 ,t 34 ,t 51 are 'a priori ' and is determined during the integration.…”