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

for pictures

Browse Source
This commit is contained in:
Brantegger Georg
2022-08-02 11:31:26 +02:00
parent 1fa71d4939
commit 1aee8debc4
5 changed files with 99 additions and 81 deletions
Download Patch File Download Diff File Expand all files Collapse all files

50
Untertweng.ipynb
View File

@@ -40,9 +40,9 @@
"alpha = np.arcsin(h_pipe/L) # Höhenwinkel der Druckrohrleitung \n",
"n = 50 # number of pipe segments in discretization # initial flow velocity [m/s]\n",
"f_D = 0.014 # Darcy friction factor\n",
"c = 500. # propagation velocity of the pressure wave [m/s]\n",
"c = 400. # propagation velocity of the pressure wave [m/s]\n",
"# consider prescribing a total simulation time and deducting the number of timesteps from that\n",
"nt = 2500 # number of time steps after initial conditions\n",
"nt = 9000 # number of time steps after initial conditions\n",
"\n",
"# derivatives of the pipeline constants\n",
"dx = L/n # length of each pipe segment\n",
@@ -126,7 +126,7 @@
"p_boundary_tur[0] = p_old[-1]\n",
"\n",
"LA_soll_vec = np.full_like(t_vec,T1.LA)\n",
"LA_soll_vec[500:]= 0\n",
"# LA_soll_vec[500:]= 0\n",
"LA_ist_vec = np.full_like(t_vec,T1.LA)\n",
"\n",
"\n"
@@ -169,8 +169,8 @@
"# loop through time steps of the pipeline\n",
"for it_pipe in range(1,pipe.nt+1):\n",
"\n",
" if it_pipe == 250:\n",
" V.set_influx(0.)\n",
" # if it_pipe == 250:\n",
" # V.set_influx(0.)\n",
"\n",
"# for each pipeline timestep, execute nt_eRK4 timesteps of the reservoir code\n",
" # set initial conditions for the reservoir time evolution calculted with e-RK4\n",
@@ -185,7 +185,7 @@
" T1.update_LA(LA_soll_vec[it_pipe])\n",
" T1.set_pressure(p_old[-1])\n",
"\n",
" LA_ist_vec[it_pipe] = T1.LA\n",
" LA_ist_vec[it_pipe] = T1.get_current_LA()\n",
"\n",
" # set boundary conditions for the next timestep of the characteristic method\n",
" p_boundary_res[it_pipe] = V.get_current_pressure()\n",
@@ -231,7 +231,7 @@
"name": "stdout",
"output_type": "stream",
"text": [
"0.0\n"
"0.7727272727272726\n"
]
}
],
@@ -296,28 +296,28 @@
"----------------------------- \n",
"Base area = 74.0 m² \n",
"Outflux area = 0.636 m² \n",
"Current level = 7.875725956447418 m\n",
"Current level = 8.020201323491174 m\n",
"Critical level low = 0.0 m \n",
"Critical level high = inf m \n",
"Volume in reservoir = -- m³ \n",
"Current influx = 0.0 m³/s \n",
"Current outflux = -0.1415386124341686 m³/s \n",
"Current outflux vel = -0.222 m/s \n",
"Current pipe pressure = 0.772 bar \n",
"Simulation timestep = 0.0004052 s \n",
"Current influx = 0.7727272727272726 m³/s \n",
"Current outflux = -3.529572535198196 m³/s \n",
"Current outflux vel = -5.548 m/s \n",
"Current pipe pressure = 0.479 bar \n",
"Simulation timestep = 0.0005065 s \n",
"Density of liquid = 1000 kg/m³ \n",
"----------------------------- \n",
"\n",
"9.22707730779877\n",
"10.57842865915012\n",
"11.92978001050147\n",
"13.281131361852822\n",
"14.632482713204173\n",
"15.983834064555523\n",
"17.335185415906874\n",
"18.686536767258225\n",
"20.037888118609576\n",
"21.389239469960927\n"
"9.47597527926513\n",
"10.931749235039087\n",
"12.387523190813043\n",
"13.843297146587\n",
"15.299071102360957\n",
"16.75484505813491\n",
"18.210619013908868\n",
"19.666392969682825\n",
"21.12216692545678\n",
"22.577940881230738\n"
]
}
],
@@ -332,7 +332,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"
},
@@ -351,7 +351,7 @@
"orig_nbformat": 4,
"vscode": {
"interpreter": {
"hash": "4a28055eb8a3160fa4c7e4fca69770c4e0a1add985300856aa3fcf4ce32a2c48"
"hash": "84fb123bdc47ab647d3782661abcbe80fbb79236dd2f8adf4cef30e8755eb2cd"
}
}
},