3.3. Delta wing
Codes: delta-wing.py
3.3.1. Diamond wing
Generate a diamond wing jet.
Figure 3.3.1 Diamond wing (1)
Figure 3.3.2 Diamond wing (2)
1def diamond_wing():
2 '''
3 Note: the surface array `[surf_x, surf_y, surf_z]` have shape `[nz, nx]`.
4
5 - `nz` is number of points in Z direction (span-wise direction).
6 - `nx` is number of points in X direction (chord-wise direction).
7 '''
8 secs, n_point = diamond_wing_parameters()
9
10 diamond = DiamondWing(secs, ns=21, projection=True)
11 diamond.add_sec([0.17468, 0.19964, 0.34936])
12 diamond.geo()
13
14 # Split the surface at leading edge
15 diamond.split([n_point])
16
17 # For inner wing lower surface: nx=0 is trailing edge, nz=0 is symmetry plan, so dyn0=0
18 diamond.smooth(0, diamond.n_sec-2, dyn0=0, ratio_end=10)
19
20 # For outer wing lower surface: nx=0 is trailing edge, nz=0 connects inner wing, so smooth0
21 diamond.smooth(diamond.n_sec-2, diamond.n_sec-1, smooth0=True, ratio_end=[10, 1, 10])
22
23 # For inner wing upper surface: nx=0 is leading edge, nz=0 is symmetry plan, so dyn0=0
24 diamond.smooth(diamond.n_sec-1, 2*diamond.n_sec-3, dyn0=0, ratio_end=10)
25
26 # For outer wing lower surface: nx=0 is leading edge, nz=0 connects inner wing, so smooth0
27 diamond.smooth(2*diamond.n_sec-3, 2*diamond.n_sec-2, smooth0=True, ratio_end=[1, 10, 10])
28
29 ax = diamond.plot(type='surface', show=False)
30 ax.view_init(elev=150, azim=-90)
31
32 plt.savefig('figures/diamond_wing.jpg', dpi=300)
33 plt.close()
34
35 # Scale by 1000
36 diamond.scale(scale=1000)
37 diamond.output_tecplot(fname=file_dump+'diamond_wing.dat', one_piece=False)
3.3.2. VFE-2
Experimental Surface Pressure Data Obtained on 65° Delta Wing Across Reynolds Number and Mach Number Ranges (Volume 2—Small-Radius Leading Edges)
https://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19960025648.pdf
Figure 3.3.3 VFE-2