reduce range

Author: 1-Bart-1

examples/stall_model.jl

@@ -107,7 +107,7 @@ plot_polars(
         "CFD_Lebesque Rey 30e5"
     ];
     literature_path_list=[path_cfd_lebesque],
-    angle_range=range(0, 25, length=25),
+    angle_range=range(0, 20, length=20),
     angle_type="angle_of_attack",
     angle_of_attack=0,
     side_slip=0,

scripts/polars.jl

@@ -20,10 +20,9 @@ try
 
     @eval @everywhere using VortexStepMethod, Xfoil, Statistics, SharedArrays
 
-    alphas = -deg2rad(5):deg2rad(1.0):deg2rad(50)
-    # alphas = -1.0:1.0:1.0
-    d_trailing_edge_angles = -deg2rad(5):deg2rad(1.0):deg2rad(50)
-    # d_trailing_edge_angles = -1.0:1.0:1.0
+    alphas = -deg2rad(5):deg2rad(1.0):deg2rad(20)
+    d_trailing_edge_angles = -deg2rad(5):deg2rad(1.0):deg2rad(20)
+    
     cl_matrix = SharedArray{Float64}((length(alphas), length(d_trailing_edge_angles)))
     cd_matrix = SharedArray{Float64}((length(alphas), length(d_trailing_edge_angles)))
     cm_matrix = SharedArray{Float64}((length(alphas), length(d_trailing_edge_angles)))
@@ -248,6 +247,30 @@ try
         end
         return nothing
     end
+
+    function plot_values(alphas, d_trailing_edge_angles, matrix, interp, name)
+        fig = plt.figure()
+        ax = fig.add_subplot(projection="3d")
+    
+        X_data = collect(d_trailing_edge_angles) .+ zeros(length(alphas))'
+        Y_data = collect(alphas)' .+ zeros(length(d_trailing_edge_angles))
+    
+        interp_matrix = similar(matrix)
+        int_alphas, int_d_trailing_edge_angles = alphas .+ deg2rad(0.5), d_trailing_edge_angles .+ deg2rad(0.5)
+        interp_matrix .= [interp(alpha, d_trailing_edge_angle) for alpha in int_alphas, d_trailing_edge_angle in int_d_trailing_edge_angles]
+        X_int = collect(int_d_trailing_edge_angles) .+ zeros(length(int_alphas))'
+        Y_int = collect(int_alphas)' .+ zeros(length(int_d_trailing_edge_angles))
+    
+        ax.plot_wireframe(X_data, Y_data, matrix, edgecolor="royalblue", lw=0.5, rstride=5, cstride=5, alpha=0.6)
+        ax.plot_wireframe(X_int, Y_int, interp_matrix, edgecolor="orange", lw=0.5, rstride=5, cstride=5, alpha=0.6)
+        plt.xlabel("Alpha")
+        plt.ylabel("Flap angle")
+        plt.zlabel("$name values")
+        plt.title("$name for different d_flap and angle")
+        plt.legend()
+        plt.grid(true)
+        plt.show()
+    end
 
     replace_nan!(cl_matrix)
     replace_nan!(cd_matrix)
@@ -256,8 +279,13 @@ try
     cl_interp = extrapolate(scale(interpolate(cl_matrix, BSpline(Linear())), alphas, d_trailing_edge_angles), NaN)
     cd_interp = extrapolate(scale(interpolate(cd_matrix, BSpline(Linear())), alphas, d_trailing_edge_angles), NaN)
     cm_interp = extrapolate(scale(interpolate(cm_matrix, BSpline(Linear())), alphas, d_trailing_edge_angles), NaN)
+
+    plot_values(alphas, d_trailing_edge_angles, cl_matrix, cl_interp, "Cl")
+    plot_values(alphas, d_trailing_edge_angles, cd_matrix, cd_interp, "Cd")
+    plot_values(alphas, d_trailing_edge_angles, cm_matrix, cm_interp, "Cm")
+    display(plot(alphas, cd_interp.(alphas, 0.0)))
 
-    @info "Relative trailing_edge height: $upper - $lower"
+    @info "Relative trailing_edge height: $(upper - lower)"
     @info "Reynolds number for flying speed of $kite_speed is $reynolds_number"
 
     toc()