AUSM_scheme.mclear allglobal PRL CRL MACHLEFT gamma pleft pright rholeft rhoright uleft... uright tend lambda % lambda = dt/dx % ..........
albada.mfigure(1), clf, hold onx = 0.0:0.01:4.0; y = (x.*x + x)./(1 + x.*x); % Graph of van Albada limiterplot(x,y)x = 0.0:0.01:1.1; y = 0.5*(1+sqrt(2))...
exact_solution.mfunction ye = exact_solution(t,x,c) % Function called: profile yyy = x;for j = 1:length(x), yyy(j) = profile(x(j) - c*t); endye = yy...
stability_domain.momega = 1 - sqrt(0.5); % Parameters in multistage omega schemealpha = (1 - 2*omega)/(1 - omega);alphap = 1 - alpha;omegap = 1 - 2*ome...
exact_solution.mfunction ye = exact_solution(x,y,D)z = (1/sqrt(2-x))*(exp(-y*y/(4*D*(2-x))) + exp(-(2-y)*(2-y)/(4*D*(2-x))));ye = z;fL.mfunction ye = fL(...
burgers_exact.mfunction burgers_exact% Plot the exact solution of Burgers' equation% u_t + u*u_x = 0, u(x,0)=f(x)u = 0.01:0.01:0.999;xp = @(u,t) u.*t + ...