@@ -476,37 +476,6 @@ function update_effective_angle_of_attack_if_VSM(body_aero::BodyAerodynamics,
476476 return alpha_array
477477end
478478
479- # function update_effective_angle_of_attack_if_VSM(body_aero::BodyAerodynamics,
480- # gamma::Vector{Float64},
481- # core_radius_fraction::Float64,
482- # z_airf_array::Matrix{Float64},
483- # x_airf_array::Matrix{Float64},
484- # va_array::Matrix{Float64},
485- # va_norm_array::Vector{Float64},
486- # va_unit_array::Matrix{Float64})
487-
488- # # Calculate AIC matrices at aerodynamic center using LLT method
489- # calculate_AIC_matrices!(
490- # body_aero, LLT, core_radius_fraction, va_norm_array, va_unit_array
491- # )
492- # AIC_x, AIC_y, AIC_z = @views body_aero.AIC[1, :, :], body_aero.AIC[2, :, :], body_aero.AIC[3, :, :]
493-
494- # # Calculate induced velocities
495- # induced_velocity = [
496- # AIC_x * gamma,
497- # AIC_y * gamma,
498- # AIC_z * gamma
499- # ]
500- # induced_velocity = hcat(induced_velocity...)
501-
502- # # Calculate relative velocities and angles
503- # relative_velocity = va_array + induced_velocity
504- # v_normal = sum(z_airf_array .* relative_velocity, dims=2)
505- # v_tangential = sum(x_airf_array .* relative_velocity, dims=2)
506- # alpha_array = atan.(v_normal ./ v_tangential)
507- # return alpha_array
508- # end
509-
510479"""
511480 calculate_results(body_aero::BodyAerodynamics, gamma_new::Vector{Float64},
512481 density::Float64, aerodynamic_model_type::Model,
0 commit comments