Just for clarity, correcting post #11 using the proper way to capture the accumulation of momentum in the ejecta from
@anuttarasammyak:
$$\sum F_{ext} = \frac{d}{dt} \left( p_{\rm{rocket}} + p_{\rm{fuel}} + p_{\rm{ejecta}} \right)$$
Mass of the rocket is constant ##M_r##
The mass of fuel lost from the rocket is gained by the ejecta ## \frac{dM_e}{dt} = -\frac{dM_f}{dt}##
The velocity of the ejecta in the rest frame is ##v_{e/O} = v_r - v_{e/r}## and is constant after being ejected, but the momentum of the ejecta is accumulating, hence:
$$\sum F_{ext} = \frac{d}{dt} \left( M_r v_r + M_f v_r + \int \frac{dM_e}{dt}( v_r - v_{e/r}) ~dt \right)$$
Applying the derivative:
$$\sum F_{ext} = \overbrace{\left( M_r\frac{dv_r}{dt} \right)}^{\rm {rocket}} +\overbrace{\left(M_f\frac{dv_r}{dt} + \frac{dM_f}{dt}v_r\right)}^{\rm{fuel}}+ \overbrace{ \frac{dM_e}{dt}( v_r - v_{e/r})}^{\rm{ejecta}} $$
$$\sum F_{ext} = \left( M_r + M_f \right) \frac{dv_r}{dt}+ \cancel{ \left( \frac{dM_f}{dt}v_r + \frac{dM_e}{dt}v_r \right)}^0 -\frac{dM_e}{dt} v_{e/r}$$
The Rocket Equation:
$$ \sum F_{ext} = \left( M_r + M_f \right) \frac{dv_r}{dt} + \frac{dM_f}{dt} v_{e/r}$$