Georg_Brantegger/Python-DT_Slot_3
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

preparation for demo video

Browse Source
This commit is contained in:
Brantegger Georg
2022-09-19 15:26:01 +02:00
parent 8917444461
commit cec300973a
1 changed files with 19 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files

37
Druckrohrleitung/Druckstoß Visualisierung.ipynb
View File

@@ -2,7 +2,7 @@
"cells": [ "cells": [
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 39, "execution_count": 1,
"metadata": {}, "metadata": {},
"outputs": [], "outputs": [],
"source": [ "source": [
@@ -22,7 +22,7 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 40, "execution_count": 2,
"metadata": {}, "metadata": {},
"outputs": [], "outputs": [],
"source": [ "source": [
@@ -77,7 +77,7 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 41, "execution_count": 3,
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
@@ -94,7 +94,7 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 42, "execution_count": 4,
"metadata": {}, "metadata": {},
"outputs": [], "outputs": [],
"source": [ "source": [
@@ -111,11 +111,12 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 43, "execution_count": 5,
"metadata": {}, "metadata": {},
"outputs": [], "outputs": [],
"source": [ "source": [
"# initialization for timeloop\n", "# initialization for timeloop\n",
"reservoir.set_influx = 0.\n",
"\n", "\n",
"level_vec = np.zeros_like(t_vec)\n", "level_vec = np.zeros_like(t_vec)\n",
"level_vec[0] = reservoir.get_current_level()\n", "level_vec[0] = reservoir.get_current_level()\n",
@@ -155,7 +156,7 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 44, "execution_count": 6,
"metadata": {}, "metadata": {},
"outputs": [], "outputs": [],
"source": [ "source": [
@@ -194,7 +195,7 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 45, "execution_count": 7,
"metadata": {}, "metadata": {},
"outputs": [], "outputs": [],
"source": [ "source": [
@@ -255,7 +256,7 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 46, "execution_count": 8,
"metadata": {}, "metadata": {},
"outputs": [], "outputs": [],
"source": [ "source": [
@@ -266,11 +267,11 @@
"axs2[0,0].set_ylabel(r'$p$ [mWS]')\n", "axs2[0,0].set_ylabel(r'$p$ [mWS]')\n",
"axs2[0,0].set_ylim([0.9*np.min(pressure_conversion(p_boundary_res,pUnit_calc,pUnit_conv)),1.1*np.max(pressure_conversion(p_boundary_res,pUnit_calc,pUnit_conv))])\n", "axs2[0,0].set_ylim([0.9*np.min(pressure_conversion(p_boundary_res,pUnit_calc,pUnit_conv)),1.1*np.max(pressure_conversion(p_boundary_res,pUnit_calc,pUnit_conv))])\n",
"\n", "\n",
"axs2[1,1].set_title('Velocity Reservoir')\n", "axs2[1,0].set_title('Velocity Reservoir')\n",
"axs2[1,1].plot(t_vec,v_boundary_res)\n", "axs2[1,0].plot(t_vec,v_boundary_res)\n",
"axs2[1,1].set_xlabel(r'$t$ [$\\mathrm{s}$]')\n", "axs2[1,0].set_xlabel(r'$t$ [$\\mathrm{s}$]')\n",
"axs2[1,1].set_ylabel(r'$v$ [$\\mathrm{m}/\\mathrm{s}$]')\n", "axs2[1,0].set_ylabel(r'$v$ [$\\mathrm{m}/\\mathrm{s}$]')\n",
"axs2[1,1].set_ylim([-1.1*np.max(v_boundary_res),1.1*np.max(v_boundary_res)])\n", "axs2[1,0].set_ylim([-1.1*np.max(v_boundary_res),1.1*np.max(v_boundary_res)])\n",
"\n", "\n",
"axs2[0,1].set_title('Pressure Turbine')\n", "axs2[0,1].set_title('Pressure Turbine')\n",
"axs2[0,1].plot(t_vec,pressure_conversion(p_boundary_tur,pUnit_calc,pUnit_conv))\n", "axs2[0,1].plot(t_vec,pressure_conversion(p_boundary_tur,pUnit_calc,pUnit_conv))\n",
@@ -278,11 +279,11 @@
"axs2[0,1].set_ylabel(r'$p$ [mWS]')\n", "axs2[0,1].set_ylabel(r'$p$ [mWS]')\n",
"axs2[0,1].set_ylim([0.9*np.min(pressure_conversion(p_boundary_tur,pUnit_calc,pUnit_conv)),1.1*np.max(pressure_conversion(p_boundary_tur,pUnit_calc,pUnit_conv))])\n", "axs2[0,1].set_ylim([0.9*np.min(pressure_conversion(p_boundary_tur,pUnit_calc,pUnit_conv)),1.1*np.max(pressure_conversion(p_boundary_tur,pUnit_calc,pUnit_conv))])\n",
"\n", "\n",
"axs2[1,0].set_title('Velocity Turbine')\n", "axs2[1,1].set_title('Velocity Turbine')\n",
"axs2[1,0].plot(t_vec,v_boundary_tur)\n", "axs2[1,1].plot(t_vec,v_boundary_tur)\n",
"axs2[1,0].set_xlabel(r'$t$ [$\\mathrm{s}$]')\n", "axs2[1,1].set_xlabel(r'$t$ [$\\mathrm{s}$]')\n",
"axs2[1,0].set_ylabel(r'$v$ [$\\mathrm{m}/\\mathrm{s}$]')\n", "axs2[1,1].set_ylabel(r'$v$ [$\\mathrm{m}/\\mathrm{s}$]')\n",
"axs2[1,0].set_ylim([-0.1,1.05*np.max(v_boundary_tur)])\n", "axs2[1,1].set_ylim([-0.1,1.05*np.max(v_boundary_tur)])\n",
"\n", "\n",
"fig2.tight_layout()\n", "fig2.tight_layout()\n",
"plt.show()" "plt.show()"