remove fuselage influences for better numerical stability and more ac…

…curate drag
peterdsharpe · Oct 5, 2023 · 9e04a1e · 9e04a1e
1 parent 7ffd0c2
commit 9e04a1e
Show file tree

Hide file tree

Showing 2 changed files with 16 additions and 14 deletions.
diff --git a/aerosandbox/aerodynamics/aero_3D/lifting_line.py b/aerosandbox/aerodynamics/aero_3D/lifting_line.py
@@ -903,13 +903,14 @@ def get_velocity_at_points(self,
 
         freestream_velocities = np.add(wide(self.steady_freestream_velocity), rotation_freestream_velocities)
 
-        if self.airplane.fuselages:
-            V_induced_fuselage = self.calculate_fuselage_influences(
-                points=points
-            )
-            V = V_induced + V_induced_fuselage + freestream_velocities
-        else:
-            V = V_induced + freestream_velocities
+        V = V_induced + freestream_velocities
+
+        # if self.airplane.fuselages:
+        #     V_induced_fuselage = self.calculate_fuselage_influences(
+        #         points=points
+        #     )
+        #     V = V + V_induced_fuselage
+
         return V
 
     def calculate_fuselage_influences(self, points: np.ndarray) -> np.ndarray:

diff --git a/aerosandbox/aerodynamics/aero_3D/nonlinear_lifting_line.py b/aerosandbox/aerodynamics/aero_3D/nonlinear_lifting_line.py
@@ -575,13 +575,14 @@ def get_velocity_at_points(self,
 
         freestream_velocities = np.add(wide(self.steady_freestream_velocity), rotation_freestream_velocities)
 
-        if self.airplane.fuselages:
-            V_induced_fuselage = self.calculate_fuselage_influences(
-                points=points
-            )
-            V = V_induced + V_induced_fuselage + freestream_velocities
-        else:
-            V = V_induced + freestream_velocities
+        V = V_induced + freestream_velocities
+
+        # if self.airplane.fuselages:
+        #     V_induced_fuselage = self.calculate_fuselage_influences(
+        #         points=points
+        #     )
+        #     V = V + V_induced_fuselage
+
         return V
 
     def calculate_fuselage_influences(self, points: np.ndarray) -> np.ndarray: