Add: solve Regular Reflection - shock polar for a given theta or beta #…

…309
CombustionToolbox · Mar 24, 2022 · 7d866db · 7d866db
1 parent 37017be
commit 7d866db
Show file tree

Hide file tree

Showing 4 changed files with 9 additions and 101 deletions.
diff --git a/Settings/functions/check_inputs.m b/Settings/functions/check_inputs.m
@@ -27,7 +27,7 @@
                 self = set_prop(self, 'pP', self.PD.pR.value); % Guess
             case {'SHOCK_I', 'SHOCK_R'} % * SHOCK_I and SHOCK_R: CALCULATE PLANAR SHOCK WAVE
                 self = check_inputs_prop(self, 'u1');
-            case 'SHOCK_POLAR' % * SHOCK_OBLIQUE: CALCULATE OBLIQUE SHOCK WAVE
+            case {'SHOCK_POLAR', 'SHOCK_POLAR_R'} % * SHOCK_POLAR: CALCULATE OBLIQUE SHOCK POLARS
                 self = check_inputs_prop(self, 'u1');
             case 'SHOCK_OBLIQUE' % * SHOCK_OBLIQUE: CALCULATE OBLIQUE SHOCK WAVE
                 self = check_inputs_prop(self, 'u1');

diff --git a/Settings/functions/postResults.m b/Settings/functions/postResults.m
@@ -78,6 +78,13 @@ function postResults(self)
 %         beta = cell2vector(mix2{i}.polar, 'beta');
 %         displaysweepresults(self, mix2{i}.polar, beta);
 %     end
+elseif strcmp(ProblemType,{'SHOCK_POLAR_R'})
+    mix3 = mix2;
+    mix2 = self.PS.str2;
+    % Shock polars incident
+    plot_shock_polar(self, mix1, mix2);
+    % Shock polars reflected
+    plot_shock_polar(self, mix2, mix3, self.PS.str2_case, mix1);
 elseif strcmp(ProblemType,{'SHOCK_I'}) && length(phi) > 1
     self.Misc.config.labelx = 'Incident velocity $u_1$ [m/s]';
     self.Misc.config.labely = 'Temperature $T$ [K]';

diff --git a/Solver/Shocks and detonations/shock_polar.m b/Solver/Shocks and detonations/shock_polar.m
@@ -53,6 +53,7 @@
     % Range values for the shock polar
     mix2.polar.p = p2;    % [bar]
     mix2.polar.Mach = M2; % [-]
+    mix2.polar.u = u2;    % [m/s]
     mix2.polar.un = u2n;  % [m/s]
     mix2.polar.ux = u2x;  % [m/s]
     mix2.polar.uy = u2y;  % [m/s]

diff --git a/Solver/Shocks and detonations/shock_polar_reflected.m b/Solver/Shocks and detonations/shock_polar_reflected.m