[[1.0423, 1.0412, 1.0393, 1.0367, 1.0333, 1.0292, 1.0234, 1.0151, 1.0043, 1.0000, 1.0000, 1.0000],
[1.0760, 1.0738, 1.0704, 1.0658, 1.0600, 1.0530, 1.0434, 1.0297, 1.0122, 1.0000, 1.0000, 1.0000],
[1.1215, 1.1180, 1.1127, 1.1055, 1.0965, 1.0856, 1.0707, 1.0497, 1.0228, 1.0000, 1.0000, 1.0000],
[1.1789, 1.1738, 1.1660, 1.1555, 1.1423, 1.1266, 1.1052, 1.0753, 1.0374, 1.0000, 1.0000, 1.0000],
[1.2494, 1.2422, 1.2312, 1.2167, 1.1986, 1.1771, 1.1481, 1.1078, 1.0572, 1.0000, 1.0000, 1.0000],
[1.3353, 1.3253, 1.3105, 1.2910, 1.2669, 1.2384, 1.2002, 1.1476, 1.0822, 1.0056, 1.0000, 1.0000],
[1.4411, 1.4282, 1.4088, 1.3830, 1.3512, 1.3136, 1.2632, 1.1942, 1.1088, 1.0101, 1.0000, 1.0000],
[1.5724, 1.5561, 1.5313, 1.4982, 1.4572, 1.4088, 1.3440, 1.2556, 1.1472, 1.0233, 1.0000, 1.0000],
[1.7323, 1.7101, 1.6785, 1.6379, 1.5888, 1.5318, 1.4572, 1.3572, 1.2358, 1.0982, 1.0000, 1.0000],
[1.9360, 1.9056, 1.8662, 1.8184, 1.7625, 1.6991, 1.6190, 1.5142, 1.3887, 1.2471, 1.0944, 1.0000],
[2.0507, 2.0158, 1.9729, 1.9223, 1.8645, 1.8000, 1.7201, 1.6176, 1.4958, 1.3584, 1.2098, 1.0546]]])
#LPCmap.Nc_data, LPCmap.alpha_data, LPCmap.Rline_data = np.meshgrid(LPCmap.Nc_vals, LPCmap.alpha_vals, LPCmap.Rline_vals, sparse=False)
LPCmap.Npts = LPCmap.NcMap.size
LPCmap.units = {}
LPCmap.units['NcMap'] = 'rpm'
LPCmap.units['WcMap'] = 'lbm/s'
# format for new regular grid interpolator:
LPCmap.param_data = []
LPCmap.output_data = []
LPCmap.param_data.append({'name': 'alphaMap', 'values': LPCmap.alphaMap,
'default': 0, 'units': None})
LPCmap.param_data.append({'name': 'NcMap', 'values': LPCmap.NcMap,
'default': 1.1, 'units': 'rpm'})
LPCmap.param_data.append({'name': 'RlineMap', 'values': LPCmap.RlineMap,
'default': 2.2, 'units': None})
2.2764, 2.2510, 2.2248, 2.1978, 2.1701, 2.1416,
2.1118, 2.0801, 2.0464, 2.0108, 1.9735
],
[
2.3771, 2.3664, 2.3557, 2.3448, 2.3339, 2.3229,
2.3118, 2.3004, 2.2887, 2.2768, 2.2646
],
[
2.4343, 2.4365, 2.4387, 2.4409, 2.4430, 2.4452,
2.4473, 2.4495, 2.4516, 2.4538, 2.4559
]]])
#LPCMap.Nc_data, LPCMap.alpha_data, LPCMap.Rline_data = np.meshgrid(LPCMap.Nc_vals, LPCMap.alpha_vals, LPCMap.Rline_vals, sparse=False)
LPCMap.Npts = LPCMap.NcMap.size
LPCMap.units = {}
LPCMap.units['NcMap'] = 'rpm'
LPCMap.units['WcMap'] = 'lbm/s'
# format for new regular grid interpolator:
LPCMap.param_data = []
LPCMap.output_data = []
LPCMap.param_data.append({
'name': 'alphaMap',
'values': LPCMap.alphaMap,
'default': 0,
'units': None
})
LPCMap.param_data.append({
],
[
29.856, 29.948, 30.013, 30.059, 30.091, 30.113, 30.128, 30.139,
30.145, 30.149, 30.150, 30.150, 30.150, 30.150, 30.150, 30.150,
30.150, 30.150, 30.150, 30.150
],
[
29.799, 29.870, 29.920, 29.955, 29.979, 29.997, 30.009, 30.017,
30.023, 30.026, 30.028, 30.029, 30.029, 30.029, 30.029, 30.029,
30.029, 30.029, 30.029, 30.029
]]])
#HPTmap.Np_data, HPTmap.alpha_data, HPTmap.PR_data = np.meshgrid(HPTmap.Nc_vals, HPTmap.alpha_vals, HPTmap.PR_vals, sparse=False)
HPTmap.Npts = HPTmap.NpMap.size
HPTmap.units = {}
HPTmap.units['NpMap'] = 'rpm'
HPTmap.units['WpMap'] = 'lbm/s'
# format for new regular grid interpolator:
HPTmap.param_data = []
HPTmap.output_data = []
HPTmap.param_data.append({
'name': 'alphaMap',
'values': HPTmap.alphaMap,
'default': 1,
'units': None
})
HPTmap.param_data.append({
],
[
43.553, 43.780, 43.977, 44.150, 44.302, 44.437, 44.556, 44.655,
44.744, 44.825, 44.898, 44.965, 45.027, 45.083, 45.134, 45.182,
45.227, 45.267, 45.304, 45.371
],
[
43.920, 44.081, 44.219, 44.339, 44.444, 44.538, 44.621, 44.694,
44.755, 44.811, 44.861, 44.907, 44.949, 44.987, 45.022, 45.054,
45.084, 45.112, 45.138, 45.183
]]])
#LPTMap.Np_data, LPTMap.alpha_data, LPTMap.PR_data = np.meshgrid(LPTMap.Nc_vals, LPTMap.alpha_vals, LPTMap.PR_vals, sparse=False)
LPTMap.Npts = LPTMap.NpMap.size
LPTMap.units = {}
LPTMap.units['NpMap'] = 'rpm'
LPTMap.units['WpMap'] = 'lbm/s'
# format for new regular grid interpolator:
LPTMap.param_data = []
LPTMap.output_data = []
LPTMap.param_data.append({
'name': 'alphaMap',
'values': LPTMap.alphaMap,
'default': 1.0,
'units': None
})
LPTMap.param_data.append({
[1.2313, 1.2363, 1.2353, 1.2283, 1.2154, 1.1968, 1.1730, 1.1449, 1.1125, 1.0764, 1.0366],
[1.3042, 1.3232, 1.3306, 1.3261, 1.3100, 1.2825, 1.2450, 1.1992, 1.1457, 1.0854, 1.0193],
[1.3479, 1.3775, 1.3917, 1.3901, 1.3728, 1.3402, 1.2944, 1.2376, 1.1708, 1.0956, 1.0133],
[1.3935, 1.4356, 1.4581, 1.4602, 1.4420, 1.4040, 1.3488, 1.2797, 1.1982, 1.1065, 1.0068],
[1.4654, 1.5127, 1.5381, 1.5409, 1.5209, 1.4789, 1.4178, 1.3413, 1.2513, 1.1501, 1.0406],
[1.5830, 1.6221, 1.6408, 1.6386, 1.6155, 1.5723, 1.5116, 1.4367, 1.3492, 1.2512, 1.1448],
[1.7258, 1.7494, 1.7565, 1.7469, 1.7208, 1.6787, 1.6229, 1.5555, 1.4778, 1.3912, 1.2969],
[1.8381, 1.8472, 1.8432, 1.8260, 1.7960, 1.7537, 1.7006, 1.6386, 1.5684, 1.4910, 1.4073],
[1.9316, 1.9312, 1.9197, 1.8973, 1.8642, 1.8209, 1.7687, 1.7091, 1.6427, 1.5702, 1.4923],
[2.0237, 2.0149, 1.9970, 1.9704, 1.9352, 1.8918, 1.8414, 1.7850, 1.7231, 1.6562, 1.5848],
[2.1043, 2.0885, 2.0659, 2.0366, 2.0008, 1.9588, 1.9114, 1.8595, 1.8033, 1.7433, 1.6797]]])
#FanMap.Nc_data, FanMap.alpha_data, FanMap.Rline_data = np.meshgrid(FanMap.Nc_vals, FanMap.alpha_vals, FanMap.Rline_vals, sparse=False)
FanMap.Npts = FanMap.NcMap.size
FanMap.units = {}
FanMap.units['NcMap'] = 'rpm'
FanMap.units['WcMap'] = 'lbm/s'
# format for new regular grid interpolator:
FanMap.param_data = []
FanMap.output_data = []
FanMap.param_data.append({'name': 'alphaMap', 'values': FanMap.alphaMap,
'default': 0, 'units': None})
FanMap.param_data.append({'name': 'NcMap', 'values': FanMap.NcMap,
'default': 0.99, 'units': 'rpm'})
FanMap.param_data.append({'name': 'RlineMap', 'values': FanMap.RlineMap,
'default': 2.2, 'units': None})
[6.2935, 6.1874, 6.0632, 5.9218, 5.7642, 5.5914, 5.4014, 5.1930, 4.9678],
[6.4390, 6.3324, 6.2162, 6.0908, 5.9568, 5.8145, 5.6627, 5.5004,
5.3284]],
[[1.7567, 1.7543, 1.7358, 1.7016, 1.6524, 1.5890, 1.5103, 1.4157, 1.3078],
[2.4783, 2.4610, 2.4214, 2.3602, 2.2787, 2.1784, 2.0579, 1.9167, 1.7585],
[2.9457, 2.9160, 2.8640, 2.7907, 2.6972, 2.5851, 2.4528, 2.2995, 2.1288],
[3.4221, 3.4045, 3.3514, 3.2639, 3.1441, 2.9945, 2.8132, 2.5998, 2.3609],
[4.0790, 4.0961, 4.0527, 3.9500, 3.7915, 3.5820, 3.3193, 3.0047, 2.6517],
[4.8443, 5.0097, 5.0468, 4.9535, 4.7350, 4.4035, 3.9581, 3.4109, 2.8024],
[5.2659, 5.4656, 5.5180, 5.4201, 5.1776, 4.8046, 4.3011, 3.6821, 2.9960],
[5.9603, 5.8925, 5.7804, 5.6258, 5.4313, 5.2000, 4.9289, 4.6166, 4.2701],
[6.2935, 6.1874, 6.0632, 5.9218, 5.7642, 5.5914, 5.4014, 5.1930, 4.9678],
[6.4390, 6.3324, 6.2162, 6.0908, 5.9568, 5.8145, 5.6627, 5.5004,
5.3284]]])
AXI5.units = {}
AXI5.units['NcMap'] = 'rpm'
AXI5.units['WcMap'] = 'lbm/s'
AXI5.Npts = AXI5.NcMap.size
# format for new regular grid interpolator:
AXI5.param_data = []
AXI5.output_data = []
AXI5.param_data.append({
'name': 'alphaMap',
'values': AXI5.alphaMap,
'default': 0,
'units': None
],
[
29.856, 29.948, 30.013, 30.059, 30.091, 30.113, 30.128, 30.139,
30.145, 30.149, 30.150, 30.150, 30.150, 30.150, 30.150, 30.150,
30.150, 30.150, 30.150, 30.150
],
[
29.799, 29.870, 29.920, 29.955, 29.979, 29.997, 30.009, 30.017,
30.023, 30.026, 30.028, 30.029, 30.029, 30.029, 30.029, 30.029,
30.029, 30.029, 30.029, 30.029
]]])
#HPT1269.Np_data, HPT1269.alpha_data, HPT1269.PR_data = np.meshgrid(HPT1269.Nc_vals, HPT1269.alpha_vals, HPT1269.PR_vals, sparse=False)
HPT1269.Npts = HPT1269.NpMap.size
HPT1269.units = {}
HPT1269.units['NpMap'] = 'rpm'
HPT1269.units['WpMap'] = 'lbm/s'
# format for new regular grid interpolator:
HPT1269.param_data = []
HPT1269.output_data = []
HPT1269.param_data.append({
'name': 'alphaMap',
'values': HPT1269.alphaMap,
'default': 1.0,
'units': None
})
HPT1269.param_data.append({
149.899, 149.899, 149.899, 149.899, 149.899, 149.899
],
[
145.352, 145.680, 145.905, 146.061, 146.169, 146.244, 146.293,
146.324, 146.339, 146.344, 146.344, 146.344, 146.344, 146.344,
146.344, 146.344, 146.344, 146.344, 146.344, 146.344
],
[
140.863, 141.131, 141.310, 141.428, 141.503, 141.547, 141.567,
141.569, 141.569, 141.569, 141.569, 141.569, 141.569, 141.569,
141.569, 141.569, 141.569, 141.569, 141.569, 141.569
]]])
#LPT2269.Np_data, LPT2269.alpha_data, LPT2269.PR_data = np.meshgrid(LPT2269.Nc_vals, LPT2269.alpha_vals, LPT2269.PR_vals, sparse=False)
LPT2269.Npts = LPT2269.NpMap.size
LPT2269.units = {}
LPT2269.units['NpMap'] = 'rpm'
LPT2269.units['WpMap'] = 'lbm/s'
# format for new regular grid interpolator:
LPT2269.param_data = []
LPT2269.output_data = []
LPT2269.param_data.append({
'name': 'alphaMap',
'values': LPT2269.alphaMap,
'default': 1.0,
'units': None
})
LPT2269.param_data.append({
[3.5, 3.55, 3.65, 3.75, 3.9, 4.1, 4.25]]) # volts, Open Circuit Voltage
# Thevenin components for Polarization Losses
battery.tC_Th = np.array([[2000., 2500., 22200., 1200., 20000., 14000., 10000.],
[1200., 20000., 40000., 20000., 30000., 20000., 25000.],
[30000., 35000., 49000., 27500., 50000., 29000., 30000.]]) # farads
battery.tR_Th = np.array([[0.0110, 0.0070, 0.0050, 0.0040, 0.0040, 0.0040, 0.0038],
[0.0030, 0.0029, 0.0029, 0.0020, 0.0029, 0.0025, 0.0024],
[0.0017, 0.0017, 0.0017, 0.0017, 0.0018, 0.0017, 0.0017]]) # ohm Ohmic Losses
battery.tR_0 = np.array([[0.0118, 0.0110, 0.0115, 0.0109, 0.0109, 0.0115, 0.0120],
[0.0090, 0.0090, 0.0091, 0.0089, 0.0092, 0.0089, 0.0089],
[0.0085, 0.0086, 0.0082, 0.0083, 0.0086, 0.0085, 0.0086]]) # ohm
battery.units = {}
battery.units['T_bp'] = 'degC'
battery.units['tU_oc'] = 'V'
battery.units['tC_Th'] = 'F'
battery.units['tR_Th'] = 'ohm'
battery.units['tR_0'] = 'ohm'
# format for new regular grid interpolator:
battery.param_data = []
battery.output_data = []
battery.param_data.append({'name': 'T_bp', 'values': battery.T_bp,
'default': 20.0, 'units': 'degC'})
battery.param_data.append({'name': 'SOC_bp', 'values': battery.SOC_bp,
'default': 0.5, 'units': None})
[2.3966, 2.4656, 2.4877, 2.4620, 2.3896,
2.2735, 2.1155, 1.9199, 1.6959],
[2.5861, 2.7183, 2.7803, 2.7686, 2.6839,
2.5309, 2.3141, 2.0425, 1.7336],
[2.6763, 2.9474, 3.1087, 3.1461, 3.0564,
2.8474, 2.5297, 2.1261, 1.6760],
[3.0338, 3.2416, 3.3490, 3.3485, 3.2401,
3.0314, 2.7305, 2.3538, 1.9309],
[3.5470, 3.5566, 3.5254, 3.4543, 3.3449,
3.2000, 3.0207, 2.8093, 2.5717],
[3.7752, 3.7298, 3.6683, 3.5912, 3.4994,
3.3935, 3.2737, 3.1402, 2.9944],
[3.8744, 3.8211, 3.7601, 3.6915, 3.6159, 3.5334, 3.4440, 3.3476, 3.2447]]])
AXI3_2.units = {}
AXI3_2.units['NcMap'] = 'rpm'
AXI3_2.units['WcMap'] = 'lbm/s'
AXI3_2.Npts = AXI3_2.NcMap.size
# format for new regular grid interpolator:
AXI3_2.param_data = []
AXI3_2.output_data = []
AXI3_2.param_data.append({'name': 'alphaMap', 'values': AXI3_2.alphaMap,
'default': 0, 'units': None})
AXI3_2.param_data.append({'name': 'NcMap', 'values': AXI3_2.NcMap,
'default': 1.0, 'units': 'rpm'})
AXI3_2.param_data.append({'name': 'RlineMap', 'values': AXI3_2.RlineMap,
1.4738, 1.4439, 1.4106, 1.3739, 1.3341, 1.2914
],
[
1.6184, 1.6123, 1.6024, 1.5888, 1.5717, 1.5510,
1.5273, 1.5010, 1.4722, 1.4410, 1.4076, 1.3720
],
[
1.6748, 1.6655, 1.6538, 1.6399, 1.6238, 1.6054,
1.5852, 1.5633, 1.5399, 1.5150, 1.4886, 1.4608
],
[
1.7224, 1.7122, 1.7008, 1.6882, 1.6744, 1.6595,
1.6434, 1.6266, 1.6089, 1.5904, 1.5711, 1.5510
]]])
NCP01.units = {}
NCP01.units['NcMap'] = 'rpm'
NCP01.units['WcMap'] = 'lbm/s'
# Get number of points in map
NCP01.Npts = NCP01.NcMap.size
# format for new regular grid interpolator:
NCP01.param_data = []
NCP01.output_data = []
NCP01.param_data.append({
'name': 'alphaMap',
'values': NCP01.alphaMap,
'default': 0,