Sunday, March 19, 2017

Bridge collapse (Peru) & stream bank erosion

Last days a video of a bridge collapse in Lima (Peru) became viral. The video shows the moment when the scour of a bridge gets scoured and the bridge collapses. Discussions about the collapse of this bridge became a trending topic in several social networks, especially in Peru.

Some discussions tried to find someone or something to blame and others focussed on whether the bridge was properly designed or not. In this post we will not discuss such topics. Whether there is someone to blame or not; whether it was properly designed or not, the fact is that the video a good example of stream bed erosion.

Due to space limitation, this post will be an introduction to stream bank erosion. In a future post we may simulate in detail the erosion process of the video and analyse whether a properly designed rip rap would have prevented the erosion (and collapse) or not.
Fig 1. Bridge collapse due to stream bank erosion in Peru (Source youtube)

Rivers are dynamic systems that change over time. One of the processes defining such change is stream bank erosion. All rivers have stream bank erosion. Even the so-called stable rivers have eroding banks. Certain events such a flooding usually trigger stream bank erosion

Stream bank erosion becomes even more complicated when the steam has bends inducing secondary flows. Even small bends create a vertical velocities profile and a helicoidally flow known as centre region cell. Besides this helicoidally flow, the outer bank of the river also experiences a so called outer bank cell flow. This outer bank cell flow is a small helicoidally flow with a direction that opposes the main helicoidally flow. Although the outer bank cell is small and weak, it plays an important role in stream bank erosion.
Fig 2. Secondary flows (Source: Blanckaert and Vriend 2004)

The full solution of the stream bank erosion is quite complicated. Thus, simplified approaches were developed based on the main mechanisms controlling stream bank erosion. The mechanisms controlling stream bank erosion may be divided in two groups: Scour and Mass failure
  • Mass failure is the process when large chunks of bank material collapse into the river. Sometimes mass failure is a consequence of local scour removal.
  • Scour failure is the direct removal of bank material

There are two main modelling approaches for analysing bank erosion due to the mentioned processes:
  • Bank failure. Bank failure is a geotechnical based model that evaluates the bank stability and its critical failure plane. Once the soil gets saturated due to the flood, several theoretical failure planes are analysed considering its resisting forces and the driving forces, in order to define the one with the lowest safety factor. If the critical safety factor is lower than one, then all the soil above the plane will be eroded by mass failure.
Fig 3. Mass failure (Source: Carey B 2014)
  • Toe scour. Toe scour considers the bank material removed by the flow. The shear stress between the flow and the bank toe is evaluated considering the critical share stress and the erodibility of the bank. If the shear stress is higher than the critical shear stress, a portion of the bank will be eroded. After the toe gets eroded, a second failure mechanism may occur. Due to the toe erosion, the soil will look like a cantilever. Depending on the size of toe scour and the volume of soil above the scour whole, a cantilever shear failure may occur.
Fig 4. Toe scour (Source: Carey B 2014)

If we get more information about the physical characteristics of the bridge and the river, in a further post we may simulate the scour process that collapsed this bridge, and analyze in more detail this case.

Monday, March 6, 2017

Automating HEC-RAS tool

Some years ago, (early 2010) I found a way to break the RAS code and to automate some features of HEC-RAS. Thus, I was able to perform sequential modelling HEC-HMSàHEC-RASàSobek1D2D and to perform uncertainty analysis within cloud (Moya Quiroga et al., 2013). At that time there was no literature about the RAS code, so I had to try several options in the Controller class and find out how they work. In October 2015, “breaking the RAS code” was published (Goodwell, C. 2015) and I hurried into getting my copy (I was among the first ones, so it included a 4 colour pen).


HEC-HMS, HEC-RAS, Sobek2D
Fig 1. HEC-RAS in sequential modelling (Source: Moya Quiroga et al., 2013)

Documentation in the book gave me important tips and I was able to develop several scripts for automating many pre-processing and post-processing tasks; not only updating parameters but also modifying geometry and post-processing results. 

HEC-RAS used to be a great tool, but now it became wonderful. However, all my scripts work under the console (no graphical user interface GUI). I developed them and I know how to use them. Thus, I decided to include a user interface in order to make more universal tools. Few years ago I released the first version with a graphic user interface GUI (AHYDRA), a tool for automating hydraulic analysis. Debugging and coding takes time, so I released and didn’t improve it. Few months ago I decided to update it and to release new improved version was released. This version includes new features and some bugs were fixed.

New automatization options for HEC-RAS

The first application that I got for this tool was to improve the boundary conditions in HEC-RAS. This  new version automatically updates data and simulates HEC-RAS considering:
·   Downstream boundary conditions (energy slope). Important for analysing the extent of the uncertainties due to BC and to select the best BC for our case.
·   Manning roughness. Important to analyse the range of potential manning roughness. Besides, it allows analysing a spatial variable Manning. For instance, the Manning range is between 0.03-0.04; we can simulate one cross section with 0.03, the next cross section with 0.034 and so on. Thus, we will be able to perform an uncertainty analysis of manning.
·     Upstream discharge boundary condition. Inflow is usually assumed as a fixed value. However, most inflow discharge data are based on steady state stage-discharge curves prone to errors. Thus, we can perform a sensitivity or uncertainty analysis of the inflow.

Breaking the HEC-RAS code
Fig 2. User interface of AHYDRA

Single and multiple simulations (Monte Carlo).

  •  The single simulations option (at the left of the GUI) allows to automate one simulation by changing one value of the mentioned possibilities.
  • Besides, this version includes the possibility to perform a Mont Carlo analysis by running several simulations and randomly updating the desired parameter within a user specified range. This option is located on the right side of the GUI.
river uncertainty

Rip rap design


The single simulation option also calculates the rip rap size (d50 in SI units) for each cross section, by considering the simulated hydraulic conditions and the USGS rip rap design criterion.

Installaing AHYDRA to automate HEC-RAS

This version includes an installer *.exe. Just double click it and follow the install wizard. After installation AHYDRA will be accessible via desktop and start up menu.

download for HEC-RAS
Scan the QR code with your tablet/smartphone to fill the free AHYDRA request form.