end of day 29.06.2022
This commit is contained in:
@@ -2,7 +2,7 @@
|
||||
"cells": [
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": null,
|
||||
"execution_count": 1,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -14,7 +14,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": null,
|
||||
"execution_count": 2,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -98,7 +98,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": null,
|
||||
"execution_count": 3,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -107,6 +107,15 @@
|
||||
"axs1[0,1].plot(t_vec,v_1)\n",
|
||||
"axs1[1,0].plot(t_vec,p_np1)\n",
|
||||
"axs1[1,1].plot(t_vec,v_np1)\n",
|
||||
"axs1[0,0].set_xlabel(r'$t$ [$\\mathrm{s}$]')\n",
|
||||
"axs1[0,0].set_ylabel(r'$p$ [Pa]')\n",
|
||||
"axs1[0,1].set_xlabel(r'$t$ [$\\mathrm{s}$]')\n",
|
||||
"axs1[0,1].set_ylabel(r'$v$ [$\\mathrm{m}/\\mathrm{s}$]')\n",
|
||||
"axs1[1,0].set_xlabel(r'$t$ [$\\mathrm{s}$]')\n",
|
||||
"axs1[1,0].set_ylabel(r'$p$ [Pa]')\n",
|
||||
"axs1[1,1].set_xlabel(r'$t$ [$\\mathrm{s}$]')\n",
|
||||
"axs1[1,1].set_ylabel(r'$v$ [$\\mathrm{m}/\\mathrm{s}$]')\n",
|
||||
"\n",
|
||||
"axs1[0,0].set_title('Pressure Reservoir')\n",
|
||||
"axs1[0,1].set_title('Velocity Reservoir')\n",
|
||||
"axs1[1,0].set_title('Pressure Turbine')\n",
|
||||
@@ -117,7 +126,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": null,
|
||||
"execution_count": 4,
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"source": [
|
||||
@@ -141,7 +150,10 @@
|
||||
"fig2,axs2 = plt.subplots(2,1)\n",
|
||||
"axs2[0].set_title('Pressure distribution in pipeline')\n",
|
||||
"axs2[1].set_title('Velocity distribution in pipeline')\n",
|
||||
"\n",
|
||||
"axs2[0].set_xlabel(r'$x$ [$\\mathrm{m}$]')\n",
|
||||
"axs2[0].set_ylabel(r'$p$ [Pa]')\n",
|
||||
"axs2[1].set_xlabel(r'$x$ [$\\mathrm{m}$]')\n",
|
||||
"axs2[1].set_ylabel(r'$p$ [Pa]')\n",
|
||||
"lo_00, = axs2[0].plot(pl_vec,pipe.p_old,marker='.')\n",
|
||||
"lo_01, = axs2[1].plot(pl_vec,pipe.v_old,marker='.')\n",
|
||||
"axs2[0].set_ylim([-20*p0,20*p0])\n",
|
||||
@@ -183,6 +195,14 @@
|
||||
"axs3[0,1].set_title('Velocity Reservoir')\n",
|
||||
"axs3[1,0].set_title('Pressure Turbine')\n",
|
||||
"axs3[1,1].set_title('Velocity Turbine')\n",
|
||||
"axs3[0,0].set_xlabel(r'$t$ [$\\mathrm{s}$]')\n",
|
||||
"axs3[0,0].set_ylabel(r'$p$ [Pa]')\n",
|
||||
"axs3[0,1].set_xlabel(r'$t$ [$\\mathrm{s}$]')\n",
|
||||
"axs3[0,1].set_ylabel(r'$v$ [$\\mathrm{m}/\\mathrm{s}$]')\n",
|
||||
"axs3[1,0].set_xlabel(r'$t$ [$\\mathrm{s}$]')\n",
|
||||
"axs3[1,0].set_ylabel(r'$p$ [Pa]')\n",
|
||||
"axs3[1,1].set_xlabel(r'$t$ [$\\mathrm{s}$]')\n",
|
||||
"axs3[1,1].set_ylabel(r'$v$ [$\\mathrm{m}/\\mathrm{s}$]')\n",
|
||||
"fig3.tight_layout()\n",
|
||||
"plt.show()"
|
||||
]
|
||||
|
||||
Reference in New Issue
Block a user