= 1.4
a(
gama, q,
k) = sqrt( 1 - (
gama - 1) * q * q / 2 )
r(
gama, q,
k) = a^(2/(gama-1))
p(gama, q,
k) = 1/gama * a^(2*gama/(gama-1))
J(gama, q,
k) = 1/a + 1/(3*a^3) + 1/(5*a^5) - 1/2*log((1+a)/(1-a))
% --------------------------------------- %
x (gama, q,
k) = 1/(2*r) * (2/
k^2 - 1/q^2) - J/2
y (gama, q,
k) = 1/(
k*r*q) * sqrt( 1 - (q/
k)^2 )
% -------- gama = 1.4 ------------------- %
x_14 = subs(
x, gama, 1.4)
y_14 = subs(y, gama, 1.4)
% -------- inflow --------%
% q = q_0 = 0.5
figure
hold on
gama
Definition: Ringleb_geom_2D.m:3
syms gama q k
Definition: Ringleb_geom_2D.m:1
x
Definition: Ringleb_geom_2D.m:31