Fix oblique shock rho2_rho1 function

Add missing brackets to oblique shock state relations
kieran-mackle · Oct 30, 2024 · a995d04 · a995d04
1 parent 5b9616f
commit a995d04
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/pysagas/cfd/oblique_prandtl_meyer.py b/src/pysagas/cfd/oblique_prandtl_meyer.py
@@ -507,7 +507,7 @@ def oblique_rho2_rho1(M1: float, beta: float, gamma: float = 1.4):
         Peter Jacobs
         """
         M1n = M1 * abs(np.sin(beta))
-        rho2_rho1 = (gamma + 1) * M1n**2 / 2 + (gamma - 1) * M1n**2
+        rho2_rho1 = (gamma + 1) * M1n**2 / (2 + (gamma - 1) * M1n**2)
         return rho2_rho1
 
     @staticmethod