Introduction to BASEMENT Basic Simulation Environment for Computation of Environmental Flow and Natural Hazard Simulation Numerical Hydraulics Autumn semester 2016 Prof. Dr. Markus Holzner Author: Pascal Corso Input: Lukas Schmidt Pascal Corso 08.11.2016 1
BASEMENT, what for? Description Numerical simulation software developed (and still under development) at the Laboratory of Hydraulics, Hydrology and Glaciology (VAW) of ETH Zurich Applications Hydrodynamics of rivers, oceans and sediment transport Flood prevention and safety River revitalization Water quality estimation Pollution treatment Pascal Corso 08.11.2016 2
Physical models Hydrodynamic Saint-Venant equations (1D) Shallow water equations (2D) Sub-surface flow - Richard s equations (3D) Closure relations Sediment transport Scalar transport equation (suspended sediments) Bedload sediment transport Lateral transport Gravity-induced transport Pascal Corso 08.11.2016 3
Shallow water equations (SWE) Assumptions made to get the SWE: Hydrostatic distribution of pressure Small slope of the channel bottom Steady-state resistance laws applicable for unsteady flows h: water depth u: depth-averaged velocity in x direction v: depth-averaged velocity in y direction 3 unknowns τ Bx, τ By : bed shear stress τ xx, τ xy, τ yx,τ yy : depth-averaged viscous and turbulent stresses D xx, D xy, D yx,d yy : momentum dispersion terms z B : bottom elevation Closure equations Pascal Corso 08.11.2016 4
Computational grid (pre-processing) BASEchain 1D BASEmesh 2D Pascal Corso 08.11.2016 5
Numerical methods for SWE Finite Volume (1/2) Finite Volume Method based on the integral form of the SWE set in conservative form: Temporal term Flux term Source terms (bed shear stress (friction), bed slope, viscous and turbulent stresses) Discretization of equations and solving method Predictor step (neglected source term): Riemann problem solver for capturing discontinuities (wave, flood propagation) Corrector step: Second-order explicit time-marching scheme Temporal and source terms constant on the control volume Special treatments for source terms Pascal Corso 08.11.2016 6
Numerical methods Finite Volume (2/2) Boundary conditions (BC) Number of BC dependent on the type of SWE If Fr 1, sub-critical flow elliptical system of equations 2 inflow and 1 ouflow conditions needed If Fr > 1, super-critical flow hyperbolic system of equations 3 inflow and 0 outflow conditions required Inflow: hydrograph Outlet: free water elevation (zero gradient), weir, gate, h-q relation, Z- hydrograph Closure relations for bed friction, viscous and turbulent stresses (Boussinesq eddy-viscosity) Pascal Corso 08.11.2016 7
Installation of BASEMENT v. 2.6 Download BASEMENT software from the Numerical Hydraulics webpage: Unzip the file and go to the BASEMENT_v2.6_linux_ubuntu14.04 directory: >> cd yourpath/basement_v2.6_linux_ubuntu14.04 Run the following command: >> sudo./setup Install the required packages: >> sudo apt-get install libvtk5.8 >> sudo apt-get f install Install dependencies: >> sudo dpkg i basement-2.6-2967.deb You can now start BASEMENT GUI by typing basement into the terminal Pascal Corso 08.11.2016 8
Installation of QGIS and BASEmesh plugin Download and install QGIS (Quantum Geographic Information System) by typing in a terminal: >> sudo apt-get install qgis Start QGIS and load the QGIS plugin manager by choosing Manage and Install Plugins in the menu Plugins Go to Settings tab, click on Add and enter: Name URL BASEmesh repository http://people.ee.ethz.ch/~basement/qgis_plugins/qgis_plugins.xml Install QTOpenGL bindings by typing in a terminal: >> sudo apt-get install python-qt4-gl Restart QGIS Go to New or Get More tab in the QGIS plugin manager and install BASEMesh plugin Pascal Corso 08.11.2016 9
Installation of ParaView ParaView is used for the post-processing of data The installation is straightforward, you just need to enter: >> sudo apt-get install paraview ParaView reads VTK format files Pascal Corso 08.11.2016 10
Example: Hydrodynamics of Flaz river (1/3) Pascal Corso 08.11.2016 11
Example: Hydrodynamics of Flaz river (2/3) Inlet and outlet discharge for the steady simulation Inflow hydrograph (BC) for the unsteady simulation Pascal Corso 08.11.2016 12
Example: Hydrodynamics of Flaz river (3/3) Pascal Corso 08.11.2016 13
Example: Morphological changes of Flaz river bed Pascal Corso 08.11.2016 14
Thanks for your attention! Questions?