[sangeet]

Hello,
I am currently trying to generate a numerical solution to a converging-diverging, subsonic-supersonic 1D nozzle flow problem. I am using Maccormack's technique to solve the Euler equations. I have already done it using the non-conservation form of the governing equations. I am now trying to solve it using the strong conservation form of the equations, but the solution keeps blowing up after a few time steps. I have checked for syntactical errors or logical errors in my C++ code but i am unable to pin down exactly what i am doing wrong.
I have attached my C++ code and a few photos of the output.
I will try to describe my code here as best as i can.
I have 5 functions in my C++ code:
1-initial_values()=assigns initial values to the primitive variables and the area value at all grid points
2-conservation_form_initialization()=calculates the solution vector and flux vector terms at all the grid points from the initial values of the primitive variables.
3-time_step()=Calculates the time step satisfying the CFL criterion for linear PDEs at each grid point and chooses the lowest.
4-time_march()=Using Mccormack's technique marches one step forward in time
5-display_values()=just displays values with some formatting
All variables are non-dimensionalized.
I believe i have included enough comments in the code to help you understand whats going on in the code.

If it helps, i am currently following the anderson book on CFD (CFD-The Basics with Applications) chapter 7.