diff --git a/Ausgleichsbecken/Main_Program.ipynb b/Ausgleichsbecken/Main_Program.ipynb index 5cf46d7..9ca27cd 100644 --- a/Ausgleichsbecken/Main_Program.ipynb +++ b/Ausgleichsbecken/Main_Program.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": 16, + "execution_count": 38, "metadata": {}, "outputs": [], "source": [ @@ -21,16 +21,16 @@ }, { "cell_type": "code", - "execution_count": 17, + "execution_count": 39, "metadata": {}, "outputs": [], "source": [ "# define constants\n", "initial_level = 5. # m\n", - "initial_influx = 0.5 # m³/s\n", + "initial_influx = 0. # m³/s\n", "initial_outflux = 0. # m³/s\n", - "initial_pipeline_pressure = 1\n", - "initial_pressure_unit = 'bar'\n", + "initial_pipeline_pressure = 5.\n", + "initial_pressure_unit = 'mWS'\n", "conversion_pressure_unit = 'mWS'\n", "\n", "area_base = 1. # m²\n", @@ -46,7 +46,7 @@ }, { "cell_type": "code", - "execution_count": 18, + "execution_count": 40, "metadata": {}, "outputs": [], "source": [ @@ -57,7 +57,7 @@ "V.set_influx(initial_influx)\n", "V.set_outflux(initial_outflux)\n", "\n", - "converted_pressure, V.pressure_unit = pressure_conversion(initial_pipeline_pressure,input_unit = initial_pressure_unit, target_unit = conversion_pressure_unit)\n", + "converted_pressure,_ = pressure_conversion(initial_pipeline_pressure,input_unit = initial_pressure_unit, target_unit = 'Pa')\n", "V.pressure = converted_pressure\n", "\n", "time_vec = np.arange(0,total_max_time,simulation_timestep)\n", @@ -66,8 +66,7 @@ "level_vec = np.empty_like(time_vec)\n", "level_vec[0] = initial_level\n", "\n", - "pressure_vec = np.full_like(time_vec,converted_pressure)*((np.sin(time_vec/5)+1)*np.exp(-time_vec/50))\n", - " \n", + "pressure_vec = np.full_like(time_vec,converted_pressure)*((np.sin(time_vec)+1)*np.exp(-time_vec/50))\n", " \n", "i_max = -1\n", "\n", @@ -82,7 +81,15 @@ " if V.level < total_min_level:\n", " i_max = i\n", " break\n", - "\n", + "\n" + ] + }, + { + "cell_type": "code", + "execution_count": 41, + "metadata": {}, + "outputs": [], + "source": [ "\n", "fig1, (ax1, ax2, ax3, ax4) = plt.subplots(4, 1)\n", "fig1.set_figheight(10)\n", @@ -98,8 +105,8 @@ "ax2.set_xlabel(r'$t$ ['+V.time_unit+']')\n", "ax2.legend()\n", "\n", - "ax3.plot(time_vec[:i_max],pressure_vec[:i_max], label='Pipeline pressure at reservoir')\n", - "ax3.set_ylabel(r'$p_{pipeline}$ ['+V.pressure_unit+']')\n", + "ax3.plot(time_vec[:i_max],pressure_conversion(pressure_vec[:i_max],'Pa',conversion_pressure_unit)[0], label='Pipeline pressure at reservoir')\n", + "ax3.set_ylabel(r'$p_{pipeline}$ ['+conversion_pressure_unit+']')\n", "ax3.set_xlabel(r'$t$ ['+V.time_unit+']')\n", "ax3.legend()\n", "\n", @@ -133,7 +140,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" }, @@ -152,7 +159,7 @@ "orig_nbformat": 4, "vscode": { "interpreter": { - "hash": "4a28055eb8a3160fa4c7e4fca69770c4e0a1add985300856aa3fcf4ce32a2c48" + "hash": "84fb123bdc47ab647d3782661abcbe80fbb79236dd2f8adf4cef30e8755eb2cd" } } },