Mit dem Impulssatz erhält man \( \vec F_{ext} = \frac{\partial}{\partial t} \int_{KV} \vec v \rho\, dV + \int_{KF} \vec v \rho (\vec v \vec n)\, dA \), wobei die Geschwindigkeiten um den Einström- bzw. Ausströmwinkel \( \varepsilon \) gedreht sind.
Betrachtet man nur die axiale Richtung erhält man \( F_{Str}=\rho\frac{\partial}{\partial t} ((\vec v)_x V_{KV}) + (\vec v_1)_x \rho (\vec v_1 \vec n_1)A_1 + (\vec v_2)_x \rho (\vec v_2 \vec n_2) A_2 \)
\( =-\rho l \frac{\partial Q}{\partial t}+(-v_{1,x})\rho (-v_{1,r}) A_1 + (-v_{2,x}) \rho (v_{2,r}) A_2 \).
Ersetzt man noch die Geschwindigkeit durch den Volumenstrom erhält man für den stationären Fall \( F_{str,stat}=\frac{\cos\varepsilon_1}{\sin\varepsilon_1 A_1} \rho Q^2 - \frac{\cos\varepsilon_2}{\sin\varepsilon_2 A_2} \rho Q^2 \)