BUG: Corrected calculation of beta angle in aerodynamic computations

Co-authored-by: kevin-alcaniz <[email protected]>

Author: MateusStano

rocketpy/rocket/aero_surface/generic_surface.py

@@ -273,7 +273,7 @@ def compute_forces_and_moments(
 
         # Angles of attack and sideslip
         alpha = np.arctan2(stream_velocity[1], stream_velocity[2])
-        beta = np.arctan2(-stream_velocity[0], stream_velocity[2])
+        beta = np.arctan2(stream_velocity[0], stream_velocity[2])
 
         # Compute aerodynamic forces and moments
         lift, side, drag, pitch, yaw, roll = self._compute_from_coefficients(
@@ -283,9 +283,9 @@ def compute_forces_and_moments(
             beta,
             stream_mach,
             reynolds,
-            omega[0],
-            omega[1],
-            omega[2],
+            omega[0],  # q
+            omega[1],  # r
+            omega[2],  # p
         )
 
         # Conversion from aerodynamic frame to body frame

rocketpy/simulation/flight.py

@@ -2699,7 +2699,7 @@ def angle_of_sideslip(self):
         # Stream velocity in standard aerodynamic frame
         stream_velocity = -self.stream_velocity_body_frame
         beta = np.arctan2(
-            -stream_velocity[:, 0],
+            stream_velocity[:, 0],
             stream_velocity[:, 2],
         )  # x-z plane
         return np.column_stack([self.time, np.rad2deg(beta)])