preparation for KW Arriach calculations
This commit is contained in:
Brantegger Georg
@@ -2,7 +2,7 @@
|
||||
"cells": [
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 188,
|
||||
"execution_count": 1,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -28,7 +28,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 189,
|
||||
"execution_count": 2,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -43,7 +43,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 190,
|
||||
"execution_count": 3,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -62,7 +62,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 191,
|
||||
"execution_count": 4,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -84,7 +84,7 @@
|
||||
"cmpr_val_UL_M2_LA = validation_data_UT['UL_T2_LA'].to_numpy(copy=True)/100.\n",
|
||||
"cmpr_val_UL_Ausl = validation_data_UT['Ausl'].to_numpy(copy=True)/100.\n",
|
||||
"\n",
|
||||
"cmpr_val_UL_Ausl[cmpr_val_UL_Ausl<0.05] = 0.\n",
|
||||
"# cmpr_val_UL_Ausl[cmpr_val_UL_Ausl<0.007] = 0.\n",
|
||||
"# cmpr_val_UL_Ausl[0] = 0.\n",
|
||||
"# cmpr_val_UL_Ausl[-1]= 0.\n",
|
||||
"\n",
|
||||
@@ -97,7 +97,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 192,
|
||||
"execution_count": 5,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -110,21 +110,21 @@
|
||||
"pUnit_conv = 'mWS' # [string] for pressure conversion in print statements and plot labels\n",
|
||||
"\n",
|
||||
" # for KW OL \n",
|
||||
"OL_T1_Q_nenn = 0.85 # [m³/s] nominal flux of turbine \n",
|
||||
"OL_T1_p_nenn = pressure_conversion(10.6,'bar',pUnit_calc) # [Pa] nominal pressure of turbine \n",
|
||||
"OL_T1_Q_nenn = 0.75 # [m³/s] nominal flux of turbine \n",
|
||||
"OL_T1_p_nenn = pressure_conversion(6.04,'bar',pUnit_calc) # [Pa] nominal pressure of turbine \n",
|
||||
"OL_T1_closingTime = 10. # [s] closing time of turbine\n",
|
||||
"\n",
|
||||
"OL_T2_Q_nenn = 0.85/2 # [m³/s] nominal flux of turbine \n",
|
||||
"OL_T2_p_nenn = pressure_conversion(10.6,'bar',pUnit_calc) # [Pa] nominal pressure of turbine \n",
|
||||
"OL_T2_Q_nenn = 0.75 # [m³/s] nominal flux of turbine \n",
|
||||
"OL_T2_p_nenn = pressure_conversion(6.04,'bar',pUnit_calc) # [Pa] nominal pressure of turbine \n",
|
||||
"OL_T2_closingTime = 10. # [s] closing time of turbine\n",
|
||||
"\n",
|
||||
" # for KW UL\n",
|
||||
"UL_T1_Q_nenn = 0.85 # [m³/s] nominal flux of turbine \n",
|
||||
"UL_T1_p_nenn = pressure_conversion(10.6,'bar',pUnit_calc) # [Pa] nominal pressure of turbine \n",
|
||||
"UL_T1_Q_nenn = 0.54 # [m³/s] nominal flux of turbine \n",
|
||||
"UL_T1_p_nenn = pressure_conversion(10.72,'bar',pUnit_calc) # [Pa] nominal pressure of turbine \n",
|
||||
"UL_T1_closingTime = 10. # [s] closing time of turbine\n",
|
||||
"\n",
|
||||
"UL_T2_Q_nenn = 0.85/2 # [m³/s] nominal flux of turbine \n",
|
||||
"UL_T2_p_nenn = pressure_conversion(10.6,'bar',pUnit_calc) # [Pa] nominal pressure of turbine \n",
|
||||
"UL_T2_Q_nenn = 0.96 # [m³/s] nominal flux of turbine \n",
|
||||
"UL_T2_p_nenn = pressure_conversion(10.72,'bar',pUnit_calc) # [Pa] nominal pressure of turbine \n",
|
||||
"UL_T2_closingTime = 10. # [s] closing time of turbine\n",
|
||||
"\n",
|
||||
" # for PI controller\n",
|
||||
@@ -137,7 +137,7 @@
|
||||
"Pip_length = (535.+478.) # [m] length of pipeline\n",
|
||||
"Pip_dia = 0.9 # [m] diameter of pipeline\n",
|
||||
"Pip_area = Pip_dia**2/4*np.pi # [m²] crossectional area of pipeline\n",
|
||||
"Pip_head = 105. # [m] hydraulic head of pipeline without reservoir\n",
|
||||
"Pip_head = 115. # [m] hydraulic head of pipeline without reservoir\n",
|
||||
"Pip_angle = np.arcsin(Pip_head/Pip_length) # [rad] elevation angle of pipeline \n",
|
||||
"Pip_n_seg = 50 # [-] number of pipe segments in discretization\n",
|
||||
"Pip_f_D = 0.014 # [-] Darcy friction factor\n",
|
||||
@@ -168,7 +168,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 193,
|
||||
"execution_count": 6,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -189,29 +189,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 194,
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
{
|
||||
"data": {
|
||||
"text/plain": [
|
||||
"[<matplotlib.lines.Line2D at 0x1ea83745220>]"
|
||||
]
|
||||
},
|
||||
"execution_count": 194,
|
||||
"metadata": {},
|
||||
"output_type": "execute_result"
|
||||
}
|
||||
],
|
||||
"source": [
|
||||
"%matplotlib qt5\n",
|
||||
"plt.figure()\n",
|
||||
"plt.plot(t_vec,val_UL_Ausl)\n"
|
||||
]
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 195,
|
||||
"execution_count": 7,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -248,7 +226,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 196,
|
||||
"execution_count": 8,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -316,7 +294,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 197,
|
||||
"execution_count": 9,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -354,27 +332,18 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 198,
|
||||
"execution_count": 10,
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
{
|
||||
"name": "stdout",
|
||||
"output_type": "stream",
|
||||
"text": [
|
||||
"0.1703345698765209\n"
|
||||
]
|
||||
}
|
||||
],
|
||||
"outputs": [],
|
||||
"source": [
|
||||
"flux_corr = (526.-444.)/trap_int(val_UL_Ausl,Pip_dt)\n",
|
||||
"print(flux_corr)\n",
|
||||
"flux_corr = (635.-444.)/trap_int(val_UL_Ausl,Pip_dt)\n",
|
||||
"\n",
|
||||
"# flux_corr = 0."
|
||||
]
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 199,
|
||||
"execution_count": 11,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -460,7 +429,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 200,
|
||||
"execution_count": 12,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -526,7 +495,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 201,
|
||||
"execution_count": 13,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -574,16 +543,16 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 202,
|
||||
"execution_count": 14,
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
{
|
||||
"data": {
|
||||
"text/plain": [
|
||||
"<matplotlib.legend.Legend at 0x1ea934c9bb0>"
|
||||
"<matplotlib.legend.Legend at 0x22076d48ac0>"
|
||||
]
|
||||
},
|
||||
"execution_count": 202,
|
||||
"execution_count": 14,
|
||||
"metadata": {},
|
||||
"output_type": "execute_result"
|
||||
}
|
||||
@@ -597,16 +566,16 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 203,
|
||||
"execution_count": 15,
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
{
|
||||
"data": {
|
||||
"text/plain": [
|
||||
"<matplotlib.legend.Legend at 0x1ea8d7dd310>"
|
||||
"<matplotlib.legend.Legend at 0x2207de5f370>"
|
||||
]
|
||||
},
|
||||
"execution_count": 203,
|
||||
"execution_count": 15,
|
||||
"metadata": {},
|
||||
"output_type": "execute_result"
|
||||
}
|
||||
@@ -623,7 +592,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 204,
|
||||
"execution_count": 16,
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
{
|
||||
@@ -631,7 +600,7 @@
|
||||
"output_type": "stream",
|
||||
"text": [
|
||||
"444.0\n",
|
||||
"448.8083535188542\n"
|
||||
"443.13951899767665\n"
|
||||
]
|
||||
}
|
||||
],
|
||||
@@ -643,7 +612,7 @@
|
||||
],
|
||||
"metadata": {
|
||||
"kernelspec": {
|
||||
"display_name": "Python 3.8.13 ('DT_Slot_3')",
|
||||
"display_name": "Python 3.8.13 ('Georg_DT_Slot3')",
|
||||
"language": "python",
|
||||
"name": "python3"
|
||||
},
|
||||
@@ -662,7 +631,7 @@
|
||||
"orig_nbformat": 4,
|
||||
"vscode": {
|
||||
"interpreter": {
|
||||
"hash": "4a28055eb8a3160fa4c7e4fca69770c4e0a1add985300856aa3fcf4ce32a2c48"
|
||||
"hash": "84fb123bdc47ab647d3782661abcbe80fbb79236dd2f8adf4cef30e8755eb2cd"
|
||||
}
|
||||
}
|
||||
},
|
||||
|
||||
Reference in New Issue
Block a user