rename Umag to v_a

Author: 1-Bart-1

data/HL5_ram_air_kite_foil.dat

@@ -2,6 +2,8 @@ using LinearAlgebra
 using ControlPlots
 using VortexStepMethod
 
+plot = false
+
 # Step 1: Define wing parameters
 n_panels = 20          # Number of panels
 span = 20.0            # Wing span [m]
@@ -53,7 +55,7 @@ println("CD = $(round(results_vsm["cd"], digits=4))")
 println("Projected area = $(round(results_vsm["projected_area"], digits=4)) m²")
 
 # Step 6: Plot geometry
-plot_geometry(
+plot && plot_geometry(
       wa,
       "Rectangular_wing_geometry";
       data_type=".pdf",
@@ -66,7 +68,7 @@ nothing
 # Step 7: Plot spanwise distributions
 y_coordinates = [panel.aerodynamic_center[2] for panel in wa.panels]
 
-plot_distribution(
+plot && plot_distribution(
     [y_coordinates, y_coordinates],
     [results_vsm, results_llt],
     ["VSM", "LLT"],
@@ -75,7 +77,7 @@ plot_distribution(
 
 # Step 8: Plot polar curves
 angle_range = range(0, 20, 20)
-plot_polars(
+plot && plot_polars(
     [llt_solver, vsm_solver],
     [wa, wa],
     ["LLT", "VSM"];

examples/rectangular_wing.jl

@@ -199,8 +199,8 @@ function gamma_loop(
     converged = false
     n_panels = wing_aero.n_panels
     alpha_array = wing_aero.alpha_array
-    Umag_array = wing_aero.Umag_array
-    Umagw_array = similar(Umag_array)
+    v_a_array = wing_aero.v_a_array
+    Umagw_array = similar(v_a_array)
 
     gamma = copy(gamma_new)
     abs_gamma_new = copy(gamma_new)
@@ -242,14 +242,14 @@ function gamma_loop(
         alpha_array .= atan.(v_normal_array, v_tangential_array)
 
         for i in 1:n_panels
-            @views Umag_array[i] = norm(relative_velocity_crossz[i, :])
+            @views v_a_array[i] = norm(relative_velocity_crossz[i, :])
             @views Umagw_array[i] = norm(Uinfcrossz_array[i, :])
         end
 
         for (i, (panel, alpha)) in enumerate(zip(panels, alpha_array))
             cl_array[i] = calculate_cl(panel, alpha)
         end
-        gamma_new .= 0.5 .* Umag_array.^2 ./ Umagw_array .* cl_array .* chord_array
+        gamma_new .= 0.5 .* v_a_array.^2 ./ Umagw_array .* cl_array .* chord_array
 
         # Apply damping if needed
         if solver.is_with_artificial_damping