add second benchmark

Author: ufechner7

examples/bench.jl

@@ -2,7 +2,12 @@ using LinearAlgebra
 using ControlPlots
 using VortexStepMethod
 
-plot = true
+using Pkg
+
+if !("CSV" ∈ keys(Pkg.project().dependencies))
+    using TestEnv
+    TestEnv.activate()
+end
 
 # Step 1: Define wing parameters
 n_panels = 20          # Number of panels
@@ -39,9 +44,38 @@ vsm_solver = Solver(aerodynamic_model_type=VSM)
 
 # Step 5: Solve using both methods
 results_vsm = solve(vsm_solver, wa)
+println("Rectangular wing:")
 @time results_vsm = solve(vsm_solver, wa)
 # time Python: 32.0 ms Ryzen 7950x
 # time Julia:   0.6 ms Ryzen 7950x
 #               0.8 ms laptop, performance mode, battery 
 
+# Create wing geometry
+wing = KiteWing("data/ram_air_kite_body.obj", "data/ram_air_kite_foil.dat")
+body_aero = BodyAerodynamics([wing])
+
+# Create solvers
+VSM_solver = Solver(
+    aerodynamic_model_type=VSM,
+    is_with_artificial_damping=false
+)
+
+# Setting velocity conditions
+v_a = 15.0
+aoa = 15.0
+side_slip = 0.0
+yaw_rate = 0.0
+aoa_rad = deg2rad(aoa)
+vel_app = [
+    cos(aoa_rad) * cos(side_slip),
+    sin(side_slip),
+    sin(aoa_rad)
+] * v_a
+set_va!(body_aero, vel_app)
+
+# Solving and plotting distributions
+results = solve(VSM_solver, body_aero)
+println("RAM-air kite:")
+@time results = solve(VSM_solver, body_aero)
+
 nothing