added file to simulate KW Hammer and added visualization for pressure surge presentation

This commit is contained in:
Brantegger Georg
2022-09-08 09:15:27 +02:00
parent da4c7beb8f
commit b9634e378b
4 changed files with 788 additions and 12849 deletions

View File

@@ -124,7 +124,11 @@
"v_boundary_res[0] = v_old[0]\n",
"v_boundary_tur[0] = v_old[-1] \n",
"p_boundary_res[0] = p_old[0]\n",
"p_boundary_tur[0] = p_old[-1]\n"
"p_boundary_tur[0] = p_old[-1]\n",
"\n",
"v_boundary_tur[:np.argmin(np.abs(t_vec-100))] = v_old[-1] \n",
"v_boundary_tur[np.argmin(np.abs(t_vec-100)):] = 0\n",
"\n"
]
},
{
@@ -171,7 +175,7 @@
" \n",
" # set boundary conditions for the next timestep of the characteristic method\n",
" p_boundary_res[it_pipe] = reservoir.get_current_pressure()\n",
" v_boundary_tur[it_pipe] = flux_init/Pip_area\n",
" # v_boundary_tur[it_pipe] = flux_init/Pip_area\n",
"\n",
" # the the boundary conditions in the pipe.object and thereby calculate boundary pressure at turbine\n",
" pipe.set_boundary_conditions_next_timestep(p_boundary_res[it_pipe],v_boundary_tur[it_pipe])\n",