Mercurial > repos > public > sbplib

comparison +scheme/Euler1d.m @ 58:24103284e09d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed coding errors in Euler1d.
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 13 Nov 2015 16:39:45 -0800
parents 9a647dcccbdd
children e431c1260f52
comparison
equal deleted inserted replaced
57:9a647dcccbdd 58:24103284e09d
77 % p=(gamma-1)*(e-rho*u^2/2); 77 % p=(gamma-1)*(e-rho*u^2/2);
78 78
79 %Solving on form q_t + F_x = 0 79 %Solving on form q_t + F_x = 0
80 function o = F(Q) 80 function o = F(Q)
81 % Flux: f = [q2; q2.^2/q1 + p(q); (q3+p(q))*q2/q1]; 81 % Flux: f = [q2; q2.^2/q1 + p(q); (q3+p(q))*q2/q1];
82 o = [Q(2,:); Q(2,:).^2/Q(1,:) + p(Q); (Q(3,:)+p(Q)).*Q(2,:)./Q(1,:)]; 82 o = [Q(2,:); Q(2,:).^2./Q(1,:) + p(Q); (Q(3,:)+p(Q)).*Q(2,:)./Q(1,:)];
83 end 83 end
84 84
85 % Equation of state 85 % Equation of state
86 function o = p(Q) 86 function o = p(Q)
87 % Pressure p = (gamma-1)*(q3-q2.^2/q1/2) 87 % Pressure p = (gamma-1)*(q3-q2.^2/q1/2)
88 o = (gamma-1)*(Q(3,:)-Q(2,:).^2/Q(1,:)/2); 88 o = (gamma-1)*( Q(3,:)-1/2*Q(2,:).^2./Q(1,:) );
89 end 89 end
90 90
91 function o = c(Q) 91 function o = c(Q)
92 % Speed of light c = sqrt(obj.gamma*p/rho); 92 % Speed of light c = sqrt(obj.gamma*p/rho);
93 o = sqrt(gamma*p(Q)./Q(1,:)); 93 o = sqrt(gamma*p(Q)./Q(1,:));