This question is related to gas-dynamics/ thermodynamics of compressible flows. The flow in the Laval nozzle is isentropic. This means that there is no irreversibility in the flow. If we consider this, the expansion in the nozzle causes a drop in the temperature and thus, resulting in a rise in the velocity of the fluid. But this happens only for a particular outlet pressure - inlet pressure ratio, 0.528. if the pressure ratio is more than this, the flow in the diffuser section (after the throat) decelerates. This is because of the Area-velocity relation (which is obtained from compressible continuity equation, and Laplace equation for velocity of sound), in terms of the Mach number. Its implications are:
- In subsonic flow regime: Convergent duct causes acceleration; divergent, deceleration.
- In supersonic regime: Convergent duct causes DECELERATION; AND DIVERGENT, ACCELERATION.
But this also is not applicable for all cases, as already mentioned, the pressure ratio must be below 0.528.
The Area-Mach number relation (which is obtained from compressible continuity equation, and the stagnation density relation to the mach number, again applicable for isentropic flows) gives the area required at the throat to "choke" the flow (i.e. to accelerate it to sonic speed at the throat) at a given pressure ratio.
The area-mach number relation gives a solution which has TWO VALUES for every area-ratio (area-ratio: ratio of the area of section under consideration to the area of the throat when flow is sonic at the throat. NOTE: flow at throat becomes sonic at throat only when Pe/P < 0.528). one value for subsonic; and the other for supersonic flows. So when the flow is supersonic, the flow accelerates; and decelerates otherwise.
So the answer to your question: Supersonic flow is accelerated by the expansion of the flow and reduction in temperature. In a nutshell, acceleration-deceleration with variation in area depends on 2 things:
1. Mach number, and
2. Pressure ratio.
If you need a more detailed explanation, I refer you to Gas Dynamics by E. Rathakrishnan for a more detailed explanation and the derivations of the area-mach number relations and so.