Simulation of a 1D Super-sonic nozzle flow

Objective: CONSERVATIVE GOVERNING EQUATIONS clear all close all clc n = 31; x = linspace(0,3,n); dx = x(2)-x(1); gamma = 1.4; % calculate initial profiles rho = 1-0.3146*x; t = 1-0.2314*x; % t = temperature v = (0.1+1.09*x).*t.^0.5; a = 1+2.2*(x-1.5).^2; s = ((t./(gamma-1) ) + ((gamma/2).*v.^2)); % time steps nt = 200; dt = 0.01;…