1-D nozzle flow equations using MacCormack method
 

Hello,

I hope someone can help me out.

This may seem like a reasonably simple question. I need to modify MacCormack's method to include the Area and Source terms in order to solve the governing equations for quasi-1-D inviscid nozzle flow:

dA*rho/dt + dA*rho*u/dx = 0

dA*rho*u/dt + dA(P + rho*u^2)/dx = P*dA/dx

dAE/dt + dA(P + E)u/dx = 0

I would be grateful if somebody could tell me how I should modify the Predictor and Corrector stages to include the additional terms.

Thanks very much.