def setUp(self):
alpha = [-10.1, -8.2, -6.1, -4.1, -2.1, 0.1, 2, 4.1, 6.2, 8.1, 10.2,
11.3, 12.1, 13.2, 14.2, 15.3, 16.3, 17.1, 18.1, 19.1, 20.1]
cl = [-0.6300, -0.5600, -0.6400, -0.4200, -0.2100, 0.0500, 0.3000,
0.5400, 0.7900, 0.9000, 0.9300, 0.9200, 0.9500, 0.9900, 1.0100,
1.0200, 1.0000, 0.9400, 0.8500, 0.7000, 0.6600]
cd = [0.0390, 0.0233, 0.0131, 0.0134, 0.0119, 0.0122, 0.0116, 0.0144,
0.0146, 0.0162, 0.0274, 0.0303, 0.0369, 0.0509, 0.0648, 0.0776,
0.0917, 0.0994, 0.2306, 0.3142, 0.3186]
cm = [-0.0044, -0.0051, 0.0018, -0.0216, -0.0282, -0.0346, -0.0405,
-0.0455, -0.0507, -0.0404, -0.0321, -0.0281, -0.0284, -0.0322,
-0.0361, -0.0363, -0.0393, -0.0398, -0.0983, -0.1242, -0.1155]
cm_zeros = np.zeros(len(cm))
Re = 1
self.polar = Polar(Re, alpha, cl, cd, cm)
self.polar2 = Polar(Re, alpha, cl, cd, cm_zeros)
def setUp(self):
alpha = [-3.04, -2.03, -1.01, 0.01, 1.03, 2.05, 3.07, 4.09, 5.11,
6.13, 7.14, 8.16, 9.17, 10.18, 11.18, 12.19, 13.18,
14.18, 15.18, 16.17, 17.14, 18.06, 19.06, 20.07, 25]
cl = [-0.071, 0.044, 0.144, 0.241, 0.338, 0.435, 0.535, 0.632,
0.728, 0.813, 0.883, 0.946, 1.001, 1.054, 1.056, 1.095,
1.138, 1.114, 1.073, 1.008, 0.95, 0.902, 0.795, 0.797, 0.8]
cd = [0.0122, 0.0106, 0.0114, 0.0134, 0.0136, 0.014, 0.0147,
0.0156, 0.0162, 0.0173, 0.0191, 0.0215, 0.0248, 0.0339,
0.0544, 0.0452, 0.0445, 0.067, 0.0748, 0.1028, 0.1473,
0.2819, 0.2819, 0.2819, 0.3]
cm = [-0.0044, -0.0051, 0.0018, -0.0216, -0.0282, -0.0346, -0.0405,
-0.0455, -0.0507, -0.0404, -0.0321, -0.0281, -0.0284, -0.0322,
-0.0361, -0.0363, -0.0393, -0.0398, -0.0983, -0.1242, -0.1155,
-0.1068, -0.0981, -0.0894, -0.0807]
Re = 1
self.polar1 = Polar(Re, alpha, cl, cd, cm)
alpha = [-3.04, -2.03, -1.01, 0.01, 1.03, 2.05, 3.07, 4.09, 5.11,
6.13, 7.14, 8.16, 9.17, 10.18, 11.18, 12.19, 13.18, 14.18,
15.189, 16.17, 17.14, 18.06, 19.06, 20.07, 21.08, 22.09,
23.1, 25]
cl = [-0.0852, 0.0528, 0.1728, 0.2892, 0.4056, 0.522, 0.642, 0.7584,
0.8736, 0.9756, 1.0596, 1.1352, 1.2012, 1.2648, 1.2672, 1.314,
1.3656, 1.3368, 1.2876, 1.2096, 1.14, 1.0824, 0.954, 0.9564, 1,
1.2, 1.4, 1.6]
cd = [0.01464, 0.01272, 0.01368, 0.01608, 0.01632, 0.0168, 0.01764,
0.01872, 0.01944, 0.02076, 0.02292, 0.0258, 0.02976, 0.04068,
0.06528, 0.05424, 0.0534, 0.0804, 0.08976, 0.12336, 0.17676,
0.33828, 0.33828, 0.33828, 0.35, 0.4, 0.45, 0.5]
cm = [-0.0037, -0.0044, -0.0051, 0.0018, -0.0216, -0.0282, -0.0346,
-0.0405, -0.0455, -0.0507, -0.0404, -0.0321, -0.0281, -0.0284,
-0.0322, -0.0361, -0.0363, -0.0393, -0.0398, -0.0983, -0.1242,
-0.1155, -0.1068, -0.0981, -0.0894, -0.0807, -0.072, -0.0633]
self.polar2 = Polar(Re, alpha, cl, cd, cm)
def setUp(self):
alpha = [-9.000, -8.000, -7.000, -6.000, -5.000, -4.000, -3.000,
-2.000, -1.000, 0.000, 1.000, 2.000, 3.000, 4.000, 5.000,
6.000, 7.000, 8.000, 9.000, 10.000, 11.000, 12.000, 13.000,
14.000, 15.000, 16.000, 17.000, 18.000, 19.000, 20.000,
30.000, 40.000, 50.000]
cl = [-0.802, -0.721, -0.611, -0.506, -0.408, -0.313, -0.220, -0.133,
-0.060, 0.036, 0.227, 0.342, 0.436, 0.556, 0.692, 0.715, 0.761,
0.830, 0.893, 0.954, 1.013, 1.042, 1.061, 1.083, 1.078, 0.882,
0.811, 0.793, 0.793, 0.798, 0.772, 0.757, 0.700]
cd = [0.027, 0.025, 0.024, 0.023, 0.022, 0.022, 0.023, 0.025, 0.027,
0.028, 0.024, 0.019, 0.017, 0.015, 0.017, 0.019, 0.021, 0.024,
0.027, 0.031, 0.037, 0.046, 0.058, 0.074, 0.088, 0.101, 0.114,
0.128, 0.142, 0.155, 0.321, 0.525, 0.742]
cm = [-0.0037, -0.0044, -0.0051, 0.0018, -0.0216, -0.0282, -0.0346,
-0.0405, -0.0455, -0.0507, -0.0404, -0.0321, -0.0281, -0.0284,
-0.0322, -0.0361, -0.0363, -0.0393, -0.0398, -0.0983, -0.1242,
-0.1155, -0.1068, -0.0981, -0.0894, -0.0807, -0.072, -0.0633,
-0.054, -0.045, -0.036, -0.22, -0.13]
cm_zeros = np.zeros(len(cm))
Re = 1
self.polar = Polar(Re, alpha, cl, cd, cm)
self.polar2 = Polar(Re, alpha, cl, cd, cm_zeros)
class TestMisc(unittest.TestCase):
def setUp(self):
alpha = [-10.1, -8.2, -6.1, -4.1, -2.1, 0.1, 2, 4.1, 6.2, 8.1, 10.2,
11.3, 12.1, 13.2, 14.2, 15.3, 16.3, 17.1, 18.1, 19.1, 20.1]
cl = [-0.6300, -0.5600, -0.6400, -0.4200, -0.2100, 0.0500, 0.3000,
0.5400, 0.7900, 0.9000, 0.9300, 0.9200, 0.9500, 0.9900, 1.0100,
1.0200, 1.0000, 0.9400, 0.8500, 0.7000, 0.6600]
cd = [0.0390, 0.0233, 0.0131, 0.0134, 0.0119, 0.0122, 0.0116, 0.0144,
0.0146, 0.0162, 0.0274, 0.0303, 0.0369, 0.0509, 0.0648, 0.0776,
0.0917, 0.0994, 0.2306, 0.3142, 0.3186]
cm = [-0.0044, -0.0051, 0.0018, -0.0216, -0.0282, -0.0346, -0.0405,
-0.0455, -0.0507, -0.0404, -0.0321, -0.0281, -0.0284, -0.0322,
-0.0361, -0.0363, -0.0393, -0.0398, -0.0983, -0.1242, -0.1155]
cm_zeros = np.zeros(len(cm))
Re = 1
self.polar = Polar(Re, alpha, cl, cd, cm)
self.polar2 = Polar(Re, alpha, cl, cd, cm_zeros)
def test_unsteady(self):
alpha0, alpha1, alpha2, cnSlope, cn1, cn2, cd0, cm0 = self.polar.unsteadyParams()
np.testing.assert_allclose(alpha0, -0.32307692307692304)
np.testing.assert_allclose(alpha1, 9.260783831245934)
np.testing.assert_allclose(alpha2, -6.779334979177289)
np.testing.assert_allclose(cnSlope, 6.4380618436681765)
np.testing.assert_allclose(cn1, 0.9201540372961516)
np.testing.assert_allclose(cn2, -0.6377683435797556)
np.testing.assert_allclose(cd0, 0.012142307692307694)
np.testing.assert_allclose(cm0, -0.03336923076923077)
def test_fully_separated(self):
cl_fs, f_st = self.polar.cl_fully_separated()
cl_fs_ref = np.array([-0.63 , -0.42017185, -0.35815607, -0.23440711, -0.11213462,
0.02669872, 0.15 , 0.2815297 , 0.41432191, 0.51685242,
0.60852946, 0.6464375 , 0.68202361, 0.7299095 , 0.76769179,
0.8037866 , 0.82370687, 0.81723832, 0.78926905, 0.69419819,
0.65999953])
f_st_ref = np.array([0.00000000e+00, 2.34199688e-01, 7.26644559e-01, 7.32580663e-01,
8.34063987e-01, 8.34063987e-01, 1.00000000e+00, 8.92315821e-01,
8.77625013e-01, 6.71133852e-01, 4.28392660e-01, 3.20122429e-01,
2.90558283e-01, 2.55881726e-01, 2.18728235e-01, 1.78134763e-01,
1.33254382e-01, 8.56818538e-02, 3.81986876e-02, 3.19820908e-03,
2.39632149e-07])
np.testing.assert_allclose(cl_fs, cl_fs_ref)
np.testing.assert_allclose(f_st, f_st_ref)
def test_cl_max(self):
cl_max, alpha_cl_max = self.polar.cl_max()
np.testing.assert_allclose(cl_max, 1.02)
np.testing.assert_allclose(alpha_cl_max, 15.3)
def test_linear_region(self):
alpha_linear_region, cl_linear_region, slope, alpha0 = self.polar.linear_region()
np.testing.assert_allclose(alpha_linear_region, np.array([-6.17381944, 7.43986639]))
np.testing.assert_allclose(cl_linear_region, np.array([-0.68718783, 0.91178174]))
np.testing.assert_allclose(slope, 0.11745309755638363)
np.testing.assert_allclose(alpha0, -0.32307692307692304)
class TestExtrap(unittest.TestCase):
def setUp(self):
alpha = [-10.1, -8.2, -6.1, -4.1, -2.1, 0.1, 2, 4.1, 6.2, 8.1, 10.2,
11.3, 12.1, 13.2, 14.2, 15.3, 16.3, 17.1, 18.1, 19.1, 20.1]
cl = [-0.6300, -0.5600, -0.6400, -0.4200, -0.2100, 0.0500, 0.3000,
0.5400, 0.7900, 0.9000, 0.9300, 0.9200, 0.9500, 0.9900, 1.0100,
1.0200, 1.0000, 0.9400, 0.8500, 0.7000, 0.6600]
cd = [0.0390, 0.0233, 0.0131, 0.0134, 0.0119, 0.0122, 0.0116, 0.0144,
0.0146, 0.0162, 0.0274, 0.0303, 0.0369, 0.0509, 0.0648, 0.0776,
0.0917, 0.0994, 0.2306, 0.3142, 0.3186]
cm = [-0.0044, -0.0051, 0.0018, -0.0216, -0.0282, -0.0346, -0.0405,
-0.0455, -0.0507, -0.0404, -0.0321, -0.0281, -0.0284, -0.0322,
-0.0361, -0.0363, -0.0393, -0.0398, -0.0983, -0.1242, -0.1155]
cm_zeros = np.zeros(len(cm))
Re = 1
self.polar = Polar(Re, alpha, cl, cd, cm)
self.polar2 = Polar(Re, alpha, cl, cd, cm_zeros)
def test_extrap1(self):
cdmax = 1.29
newpolar = self.polar.extrapolate(cdmax=cdmax)
alpha_extrap = [-180, -170, -160, -150, -140, -130, -120, -110, -100,
-90, -80, -70, -60, -50, -40, -30, -20, -10.1, -8.2,
-6.1, -4.1, -2.1, 0.1, 2, 4.1, 6.2, 8.1, 10.2, 11.3,
12.1, 13.2, 14.2, 15.3, 16.3, 17.1, 18.1, 19.1, 20.1,
30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140,
150, 160, 170, 180]
cl_extrap = [0.0000, 0.2299, 0.4597, 0.4907, 0.5053, 0.4805, 0.4102,
0.2985, 0.1565, 0.0000, -0.1565, -0.2985, -0.4102,
-0.4805, -0.5053, -0.4907, -0.4637, -0.6300, -0.5600,
-0.6400, -0.4200, -0.2100, 0.0500, 0.3000, 0.5400,
0.7900, 0.9000, 0.9300, 0.9200, 0.9500, 0.9900, 1.0100,
1.0200, 1.0000, 0.9400, 0.8500, 0.7000, 0.6600, 0.7010,
0.7219, 0.6864, 0.5860, 0.4264, 0.2235, 0.0000, -0.1565,
-0.2985, -0.4102, -0.4805, -0.5053, -0.4907, -0.4597,
-0.2299, 0.0000]
cd_extrap = [0.1770, 0.2132, 0.3173, 0.4758, 0.6686, 0.8708, 1.0560,
1.1996, 1.2818, 1.2900, 1.2818, 1.1996, 1.0560, 0.8708,
0.6686, 0.4758, 0.3158, 0.0390, 0.0233, 0.0131, 0.0134,
0.0119, 0.0122, 0.0116, 0.0144, 0.0146, 0.0162, 0.0274,
0.0303, 0.0369, 0.0509, 0.0648, 0.0776, 0.0917, 0.0994,
0.2306, 0.3142, 0.3186, 0.4758, 0.6686, 0.8708, 1.0560,
1.1996, 1.2818, 1.2900, 1.2818, 1.1996, 1.0560, 0.8708,
0.6686, 0.4758, 0.3173, 0.2132, 0.1770]
cm_extrap = [0.0000, 0.4000, 0.2431, 0.2568, 0.2865, 0.3185, 0.3458,
0.3632, 0.3672, 0.3559, 0.3443, 0.3182, 0.2808, 0.2362,
0.1886, 0.1414, 0.0942, -0.0044, -0.0051, 0.0018, -0.0216,
-0.0282, -0.0346, -0.0405, -0.0455, -0.0507, -0.0404, -0.0321,
-0.0281, -0.0284, -0.0322, -0.0361, -0.0363, -0.0393, -0.0398,
-0.0983, -0.1242, -0.1155, -0.1710, -0.2202, -0.2637, -0.3002,
-0.3284, -0.3471, -0.3559, -0.3672, -0.3632, -0.3458, -0.3185,
-0.2865, -0.2568, -0.2431, -0.5000, 0.0000]
# re-interpolate b/c angles of attack are different
cl = np.interp(alpha_extrap, newpolar.alpha, newpolar.cl)
cd = np.interp(alpha_extrap, newpolar.alpha, newpolar.cd)
cm = np.interp(alpha_extrap, newpolar.alpha, newpolar.cm)
# test equality
np.testing.assert_allclose(cl, cl_extrap, atol=1.5e-4)
np.testing.assert_allclose(cd, cd_extrap, atol=1.5e-4)
np.testing.assert_allclose(cm, cm_extrap, atol=5e-3)
def test_extrap2(self):
cdmax = 1.0
newpolar = self.polar.extrapolate(cdmax=cdmax)
alpha_extrap = [-180, -170, -160, -150, -140, -130, -120, -110, -100,
-90, -80, -70, -60, -50, -40, -30, -20, -10.1, -8.2,
-6.1, -4.1, -2.1, 0.1, 2, 4.1, 6.2, 8.1, 10.2, 11.3,
12.1, 13.2, 14.2, 15.3, 16.3, 17.1, 18.1, 19.1, 20.1,
30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140,
150, 160, 170, 180]
cl_extrap = [0.0000, 0.2299, 0.4597, 0.4411, 0.4287, 0.3943, 0.3297,
0.2364, 0.1225, 0.0000, -0.1225, -0.2364, -0.3297,
-0.3943, -0.4287, -0.4411, -0.4637, -0.6300, -0.5600,
-0.6400, -0.4200, -0.2100, 0.0500, 0.3000, 0.5400,
0.7900, 0.9000, 0.9300, 0.9200, 0.9500, 0.9900, 1.0100,
1.0200, 1.0000, 0.9400, 0.8500, 0.7000, 0.6600, 0.6302,
0.6124, 0.5633, 0.4710, 0.3378, 0.1750, 0.0000, -0.1225,
-0.2364, -0.3297, -0.3943, -0.4287, -0.4411, -0.4597,
-0.2299, 0.0000]
cd_extrap = [0.2135, 0.2404, 0.3176, 0.4349, 0.5767, 0.7241, 0.8568,
0.9560, 1.0069, 1.0000, 1.0069, 0.9560, 0.8568, 0.7241,
0.5767, 0.4349, 0.3158, 0.0390, 0.0233, 0.0131, 0.0134,
0.0119, 0.0122, 0.0116, 0.0144, 0.0146, 0.0162, 0.0274,
0.0303, 0.0369, 0.0509, 0.0648, 0.0776, 0.0917, 0.0994,
0.2306, 0.3142, 0.3186, 0.4349, 0.5767, 0.7241, 0.8568,
0.9560, 1.0069, 1.0000, 1.0069, 0.9560, 0.8568, 0.7241,
0.5767, 0.4349, 0.3176, 0.2404, 0.2135]
cm_extrap = [0.0000, 0.4000, 0.2432, 0.2354, 0.2500, 0.2695, 0.2864,
0.2961, 0.2956, 0.2834, 0.2776, 0.2603, 0.2337, 0.2013,
0.1663, 0.1310, 0.0942, -0.0044, -0.0051, 0.0018, -0.0216,
-0.0282, -0.0346, -0.0405, -0.0455, -0.0507, -0.0404, -0.0321,
-0.0281, -0.0284, -0.0322, -0.0361, -0.0363, -0.0393, -0.0398,
-0.0983, -0.1242, -0.1155, -0.1577, -0.1930, -0.2239, -0.2494,
-0.2683, -0.2798, -0.2834, -0.2956, -0.2961, -0.2864, -0.2695,
-0.2500, -0.2354, -0.2432, -0.5000, 0.0000]
# re-interpolate b/c angles of attack are different
cl = np.interp(alpha_extrap, newpolar.alpha, newpolar.cl)
cd = np.interp(alpha_extrap, newpolar.alpha, newpolar.cd)
cm = np.interp(alpha_extrap, newpolar.alpha, newpolar.cm)
# test equality
np.testing.assert_allclose(cl, cl_extrap, atol=1.5e-4)
np.testing.assert_allclose(cd, cd_extrap, atol=1.5e-4)
np.testing.assert_allclose(cm, cm_extrap, atol=5e-3)
def test_extrap3(self):
cdmax = 1.5
newpolar = self.polar.extrapolate(cdmax)
alpha_extrap = [-180, -170, -160, -150, -140, -130, -120, -110, -100,
-90, -80, -70, -60, -50, -40, -30, -20, -10.1, -8.2,
-6.1, -4.1, -2.1, 0.1, 2, 4.1, 6.2, 8.1, 10.2, 11.3,
12.1, 13.2, 14.2, 15.3, 16.3, 17.1, 18.1, 19.1, 20.1,
30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140,
150, 160, 170, 180]
cl_extrap = [0.0000, 0.2299, 0.4597, 0.5266, 0.5608, 0.5429, 0.4685,
0.3434, 0.1810, 0.0000, -0.1810, -0.3434, -0.4685,
-0.5429, -0.5608, -0.5266, -0.4637, -0.6300, -0.5600,
-0.6400, -0.4200, -0.2100, 0.0500, 0.3000, 0.5400, 0.7900,
0.9000, 0.9300, 0.9200, 0.9500, 0.9900, 1.0100, 1.0200,
1.0000, 0.9400, 0.8500, 0.7000, 0.6600, 0.7523, 0.8012,
0.7756, 0.6693, 0.4906, 0.2586, 0.0000, -0.1810, -0.3434,
-0.4685, -0.5429, -0.5608, -0.5266, -0.4597, -0.2299,
0.0000]
cd_extrap = [0.1506, 0.1936, 0.3170, 0.5054, 0.7351, 0.9771, 1.2003,
1.3760, 1.4809, 1.5000, 1.4809, 1.3760, 1.2003, 0.9771,
0.7351, 0.5054, 0.3158, 0.0390, 0.0233, 0.0131, 0.0134,
0.0119, 0.0122, 0.0116, 0.0144, 0.0146, 0.0162, 0.0274,
0.0303, 0.0369, 0.0509, 0.0648, 0.0776, 0.0917, 0.0994,
0.2306, 0.3142, 0.3186, 0.5054, 0.7351, 0.9771, 1.2003,
1.3760, 1.4809, 1.5000, 1.4809, 1.3760, 1.2003, 0.9771,
0.7351, 0.5054, 0.3170, 0.1936, 0.1506]
cm_extrap = [0.0000, 0.4000, 0.2431, 0.2723, 0.3130, 0.3540, 0.3888,
0.4118, 0.4190, 0.4084, 0.3926, 0.3602, 0.3148, 0.2614,
0.2049, 0.1488, 0.0942, -0.0044, -0.0051, 0.0018, -0.0216,
-0.0282, -0.0346, -0.0405, -0.0455, -0.0507, -0.0404, -0.0321,
-0.0281, -0.0284, -0.0322, -0.0361, -0.0363, -0.0393, -0.0398,
-0.0983, -0.1242, -0.1155, -0.1807, -0.2399, -0.2925, -0.3370,
-0.3719, -0.3959, -0.4084, -0.4190, -0.4118, -0.3888, -0.3540,
-0.3130, -0.2723, -0.2431, -0.5000, 0.0000]
# re-interpolate b/c angles of attack are different
cl = np.interp(alpha_extrap, newpolar.alpha, newpolar.cl)
cd = np.interp(alpha_extrap, newpolar.alpha, newpolar.cd)
cm = np.interp(alpha_extrap, newpolar.alpha, newpolar.cm)
# test equality
np.testing.assert_allclose(cl, cl_extrap, atol=1.5e-4)
np.testing.assert_allclose(cd, cd_extrap, atol=1.5e-4)
np.testing.assert_allclose(cm, cm_extrap, atol=5e-3)
class Test3DStall(unittest.TestCase):
def setUp(self):
alpha = [-9.000, -8.000, -7.000, -6.000, -5.000, -4.000, -3.000,
-2.000, -1.000, 0.000, 1.000, 2.000, 3.000, 4.000, 5.000,
6.000, 7.000, 8.000, 9.000, 10.000, 11.000, 12.000, 13.000,
14.000, 15.000, 16.000, 17.000, 18.000, 19.000, 20.000,
30.000, 40.000, 50.000]
cl = [-0.802, -0.721, -0.611, -0.506, -0.408, -0.313, -0.220, -0.133,
-0.060, 0.036, 0.227, 0.342, 0.436, 0.556, 0.692, 0.715, 0.761,
0.830, 0.893, 0.954, 1.013, 1.042, 1.061, 1.083, 1.078, 0.882,
0.811, 0.793, 0.793, 0.798, 0.772, 0.757, 0.700]
cd = [0.027, 0.025, 0.024, 0.023, 0.022, 0.022, 0.023, 0.025, 0.027,
0.028, 0.024, 0.019, 0.017, 0.015, 0.017, 0.019, 0.021, 0.024,
0.027, 0.031, 0.037, 0.046, 0.058, 0.074, 0.088, 0.101, 0.114,
0.128, 0.142, 0.155, 0.321, 0.525, 0.742]
cm = [-0.0037, -0.0044, -0.0051, 0.0018, -0.0216, -0.0282, -0.0346,
-0.0405, -0.0455, -0.0507, -0.0404, -0.0321, -0.0281, -0.0284,
-0.0322, -0.0361, -0.0363, -0.0393, -0.0398, -0.0983, -0.1242,
-0.1155, -0.1068, -0.0981, -0.0894, -0.0807, -0.072, -0.0633,
-0.054, -0.045, -0.036, -0.22, -0.13]
cm_zeros = np.zeros(len(cm))
Re = 1
self.polar = Polar(Re, alpha, cl, cd, cm)
self.polar2 = Polar(Re, alpha, cl, cd, cm_zeros)
def test_stall1(self):
R = 2.4
r = 0.25*R
chord = 0.18
Omega = 200*pi/30
Uinf = 10.0
tsr = Omega*R/Uinf
newpolar = self.polar.correction3D(r/R, chord/r, tsr,
alpha_max_corr=30,
alpha_linear_min=-4,
alpha_linear_max=4)
cl_3d = [-0.8466, -0.7523, -0.6420, -0.5342, -0.4302, -0.3284,
-0.2276, -0.1303, -0.0404, 0.0618, 0.2191, 0.3321, 0.4336,
0.5501, 0.6755, 0.7363, 0.8101, 0.8973, 0.9810, 1.0640,
1.1450, 1.2098, 1.2682, 1.3281, 1.3731, 1.3088, 1.3159,
1.3534, 1.4010, 1.4515, 1.9140, 1.8857, 1.6451]
# Eggers method
cd_3d = [0.0399, 0.0334, 0.0316, 0.0293, 0.0269, 0.0254, 0.0246,
0.0246, 0.0246, 0.0252, 0.0249, 0.0200, 0.0167, 0.0157,
0.0174, 0.0183, 0.0212, 0.0255, 0.0303, 0.0367, 0.0465,
0.0615, 0.0800, 0.1047, 0.1301, 0.1695, 0.2047, 0.2384,
0.2728, 0.3081, 0.8097, 1.2625, 1.6280]
# # Du Selig method
# cd_3d = [0.027, 0.024, 0.023, 0.021, 0.02, 0.02, 0.021, 0.024, 0.027, 0.028, 0.023, 0.016, 0.013, 0.011, 0.013, 0.016, 0.019, 0.023, 0.027, 0.032, 0.04, 0.052, 0.068, 0.09, 0.109, 0.126, 0.144, 0.162, 0.181, 0.199, 0.422, 0.696, 0.987]
# test equality
np.testing.assert_allclose(newpolar.cl, cl_3d, atol=1e-3, rtol=1e-3)
np.testing.assert_allclose(newpolar.cd, cd_3d, atol=1e-3, rtol=1e-3)
def test_stall2(self):
R = 2.4
r = 0.75*R
chord = 0.28
Omega = 200*pi/30
Uinf = 14.0
tsr = Omega*R/Uinf
newpolar = self.polar.correction3D(r/R, chord/r, tsr,
alpha_max_corr=30,
alpha_linear_min=-4,
alpha_linear_max=4)
cl_3d = [-0.81340155, -0.72876051, -0.61903798, -0.51322348,
-0.41336822, -0.31696485, -0.22214149, -0.13269893,
-0.05485453, 0.04222704, 0.22525537, 0.33917483,
0.43518608, 0.55464051, 0.68785835, 0.72023796,
0.77302335, 0.84665343, 0.91485674, 0.98191931, 1.04592758,
1.08446883, 1.11313747, 1.14423161, 1.15194066, 0.98921407,
0.93776667, 0.93384528, 0.94558296, 0.96199091, 1.05910388,
1.04054486, 0.93735382]
# Eggers method
cd_3d = [0.03050922, 0.02712935, 0.02589588, 0.02453937, 0.02341344,
0.02320787, 0.02359745, 0.02497252, 0.02653913, 0.02751806,
0.02430795, 0.01935093, 0.01663156, 0.01552516, 0.01698944,
0.01853615, 0.02107760, 0.02443710, 0.02784230, 0.03217433,
0.03929881, 0.05021192, 0.06322801, 0.08159739, 0.09837902,
0.11798276, 0.13692472, 0.15565820, 0.17470667, 0.19368328,
0.44408310, 0.71034295, 0.96437541]
# # Du Selig method
# cd_3d = [0.027, 0.025, 0.024, 0.023, 0.022, 0.022, 0.023, 0.025, 0.027, 0.028, 0.024, 0.019, 0.017, 0.015, 0.017, 0.019, 0.021, 0.024, 0.027, 0.031, 0.037, 0.046, 0.059, 0.075, 0.089, 0.102, 0.116, 0.13, 0.144, 0.157, 0.326, 0.534, 0.755]
# test equality
np.testing.assert_allclose(newpolar.cl, cl_3d, atol=1e-3)
np.testing.assert_allclose(newpolar.cd, cd_3d, atol=1e-3)
def test_stall3(self):
R = 5.0
r = 0.5*R
chord = 0.5
Omega = 100*pi/30
Uinf = 10.0
tsr = Omega*R/Uinf
newpolar = self.polar.correction3D(r/R, chord/r, tsr,
alpha_max_corr=30,
alpha_linear_min=-4,
alpha_linear_max=4)
cl_3d = [-0.8240, -0.7363, -0.6264, -0.5199, -0.4188, -0.3206, -0.2239,
-0.1319, -0.0502, 0.0485, 0.2233, 0.3369, 0.4347, 0.5532,
0.6839, 0.7254, 0.7849, 0.8629, 0.9361, 1.0082, 1.0777,
1.1246, 1.1628, 1.2031, 1.2228, 1.0916, 1.0589, 1.0682,
1.0914, 1.1188, 1.3329, 1.3112, 1.1640]
# Eggers method
cd_3d = [0.0335, 0.0291, 0.0277, 0.0261, 0.0245, 0.0239, 0.0239,
0.0249, 0.0259, 0.0268, 0.0245, 0.0195, 0.0167, 0.0156,
0.0171, 0.0185, 0.0211, 0.0248, 0.0286, 0.0336, 0.0416,
0.0538, 0.0686, 0.0890, 0.1085, 0.1345, 0.1586, 0.1822,
0.2061, 0.2303, 0.5612, 0.8872, 1.1769]
# # Du Selig method
# cd_3d = [0.027, 0.025, 0.024, 0.022, 0.021, 0.021, 0.022, 0.025, 0.027, 0.028, 0.024, 0.018, 0.016, 0.014, 0.016, 0.018, 0.02, 0.024, 0.027, 0.031, 0.038, 0.048, 0.061, 0.079, 0.095, 0.109, 0.123, 0.139, 0.155, 0.169, 0.353, 0.58, 0.821]
# test equality
np.testing.assert_allclose(newpolar.cl, cl_3d, atol=1e-3)
np.testing.assert_allclose(newpolar.cd, cd_3d, atol=1e-3)
def test_stall4_cm(self):
R = 5.0
r = 0.5*R
chord = 0.5
Omega = 100*pi/30
Uinf = 10.0
tsr = Omega*R/Uinf
newpolar = self.polar2.correction3D(r/R, chord/r, tsr,
alpha_max_corr=30,
alpha_linear_min=-4,
alpha_linear_max=4)
cl_3d = [-0.8240, -0.7363, -0.6264, -0.5199, -0.4188, -0.3206, -0.2239,
-0.1319, -0.0502, 0.0485, 0.2233, 0.3369, 0.4347, 0.5532,
0.6839, 0.7254, 0.7849, 0.8629, 0.9361, 1.0082, 1.0777,
1.1246, 1.1628, 1.2031, 1.2228, 1.0916, 1.0589, 1.0682,
1.0914, 1.1188, 1.3329, 1.3112, 1.1640]
# Eggers method
cd_3d = [0.0335, 0.0291, 0.0277, 0.0261, 0.0245, 0.0239, 0.0239,
0.0249, 0.0259, 0.0268, 0.0245, 0.0195, 0.0167, 0.0156,
0.0171, 0.0185, 0.0211, 0.0248, 0.0286, 0.0336, 0.0416,
0.0538, 0.0686, 0.0890, 0.1085, 0.1345, 0.1586, 0.1822,
0.2061, 0.2303, 0.5612, 0.8872, 1.1769]
# # Du Selig method
# cd_3d = [0.027, 0.025, 0.024, 0.022, 0.021, 0.021, 0.022, 0.025, 0.027, 0.028, 0.024, 0.018, 0.016, 0.014, 0.016, 0.018, 0.02, 0.024, 0.027, 0.031, 0.038, 0.048, 0.061, 0.079, 0.095, 0.109, 0.123, 0.139, 0.155, 0.169, 0.353, 0.58, 0.821]
# cm = [-0.0037, -0.0044, -0.0051, 0.0018, -0.0216, -0.0282, -0.0346,
# -0.0405, -0.0455, -0.0507, -0.0404, -0.0321, -0.0281, -0.0284,
# -0.0322, -0.0361, -0.0363, -0.0393, -0.0398, -0.0983, -0.1242,
# -0.1155, -0.1068, -0.0981, -0.0894, -0.0807, -0.072, -0.0633,
# -0.054, -0.045, -0.036, -0.22, -0.13]
cm_zeros = np.zeros(len(cd_3d))
# test equality
np.testing.assert_allclose(newpolar.cl, cl_3d, atol=1e-3)
np.testing.assert_allclose(newpolar.cd, cd_3d, atol=1e-3)
np.testing.assert_allclose(newpolar.cm, cm_zeros, atol=1e-3)