Merge pull request #1015 from CombustionToolbox/develop

Solve: error calculation in subsonic region for rocketFAC method
CombustionToolbox · Feb 5, 2025 · e00bb03 · e00bb03
2 parents efecd5f + 8d04ee9
commit e00bb03
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/+combustiontoolbox/+rocket/@RocketSolver/rocketFAC.m b/+combustiontoolbox/+rocket/@RocketSolver/rocketFAC.m
@@ -30,6 +30,7 @@
     mix2_inj = copy(mix1);
     areaRatio = mix1.areaRatio;
     areaRatioChamber = mix1.areaRatioChamber;
+    temp_FLAG_SUBSONIC = obj.FLAG_SUBSONIC;
     obj.FLAG_SUBSONIC = true;
 
     % Unpack additional input parameters
@@ -125,7 +126,7 @@
     mix2_c.areaRatioChamber = areaRatioChamber; % [-]
 
     % Compute chemical equilibria at the given exit points
-    obj.FLAG_SUBSONIC = false;
+    obj.FLAG_SUBSONIC = temp_FLAG_SUBSONIC;
     mix4 = rocketExit(obj, mix2_c, mix3, mix4_guess, areaRatio, mix2_inj);
 
     % Compute rocket parameters