@@ -61,12 +61,15 @@ plt.ioff()
6161 rm (" /tmp/Rectangular_wing_geometry_side_view.pdf" )
6262 @test isfile (" /tmp/Rectangular_wing_geometry_top_view.pdf" )
6363 rm (" /tmp/Rectangular_wing_geometry_top_view.pdf" )
64+
6465 # Step 5: Initialize the solvers
6566 vsm_solver = Solver (aerodynamic_model_type= " VSM" )
6667 llt_solver = Solver (aerodynamic_model_type= " LLT" )
68+
6769 # Step 6: Solve the VSM and LLT
6870 results_vsm = solve (vsm_solver, wa)
6971 results_llt = solve (llt_solver, wa)
72+
7073 # Step 7: Plot spanwise distributions
7174 y_coordinates = [panel. aerodynamic_center[2 ] for panel in wa. panels]
7275
@@ -77,6 +80,26 @@ plt.ioff()
7780 title= " Spanwise Distributions"
7881 )
7982 @test fig isa plt. PyPlot. Figure
83+
84+ # Step 8: Plot polar curves
85+ Umag = 20.0 # Magnitude of inflow velocity [m/s]
86+ angle_range = range (0 , 20 , 20 )
87+ fig = plot_polars (
88+ [llt_solver, vsm_solver],
89+ [wa, wa],
90+ [" LLT" , " VSM" ],
91+ angle_range= angle_range,
92+ angle_type= " angle_of_attack" ,
93+ Umag= Umag,
94+ title= " Rectangular Wing Polars" ,
95+ data_type= " .pdf" ,
96+ save_path= " /tmp" ,
97+ is_save= true ,
98+ is_show= false
99+ )
100+ @test fig isa plt. PyPlot. Figure
101+ @test isfile (" /tmp/Rectangular_Wing_Polars.pdf" )
102+ rm (" /tmp/Rectangular_Wing_Polars.pdf" )
80103 end
81104end
82105plt. ion ()
0 commit comments