**Background Information**

**Questions:**

**What friction models are used in RAMMS?****When and why does the Voellmy model work well?****Why use a hydrograph for debris flow modelling?****What numerical solution technique does RAMMS use?****How long does it take to perform a simulation?**

For** **more information about important terms please consult the glossary.** **

** **

RAMMS employs a Voellmy-fluid friction model. This model divides the frictional resistance into two parts: a dry-Coulomb type friction (coefficient μ ) that scales with the normal stress and a velocity squared drag (coefficient ξ) . The frictional resistance S (Pa) is then

where ρ is the flow density, g gravitational acceleration, φ the slope angle, H the flow height and U the flow velocity. This model has found wide application in the simulation of mass movements, especially snow avalanches. The Voellmy model has been in use in

**Cohesion**

Since Version 1.6.20 the basic Voellmy equation has been modified to include cohesion:

where *C* is the cohesion of the flowing material. Unlike a standard Mohr-Coulomb type relation this formula ensures that *S*→0 when both *N*→0 and *U*→0. It increases the shear stress and therefore causes the avalanche or debris flow to stop earlier, depending on the value of *C*.

This formula has been established using chute experiments with flowing snow (Platzer et al., 2007a and Platzer et al., 2007b) and real scale experiments with debris flows in Illgraben (VS). Snow has different cohesive properties depending on snow temperature. Wet snow avalanches have higher cohesion values; dry snow avalanches have lower cohesion values.

Cohesion can help reduce spurious numerical diffusion in runout zones, providing a clearer delineation of the deposition zone.

Cohesion values (unit Pascal) may be entered in the Mu/Xi tab of the Run Simulation window. Recommended values may be found in the following:

- Avalanche, dry snow: 0 - 100 Pa
- Avalanche, wet snow: 100 - 300 Pa
- Debris Flow: 0 - 2000 Pa

Please use cohesion values with care!

**Curvature**

Since Version 1.6.20, the normal force *N* now includes centrifugal forces arising from the terrain curvature. We use the method proposed by Fischer et al. (2012) which was specifically developed for RAMMS. The centrifugal acceleration ƒ is both a function of the avalanche velocity and terrain curvature. The acceleration is calculated according to

where μ is the vector μ = (u,v), consisting of the avalanche velocity in the x- and y-directions. The matrix **K** describes the track curvature in all directions, including the track “twist”. The centrifugal force is then

which is added to the normal force *N*. Typically this increases the friction, causing the avalanche to slow down in tortuous and twisted flow paths. It can change the location of the deposition once the flow leaves the gully. Curvature may be activated/deactivated via the menu ‘Help → Advanced… → Curvature’.

**Avalanches:**

The Voellmy model – coupled with the calibrated parameters – can be used to (1) predict the runout distance and (2) predict the maximum flow velocity of extreme, large snow avalanches. This is one of the important research results from the Vallée de la Sionne test site.

The Voellmy parameters that we recommend describe the front of a dry-snow avalanche. Because the front defines the runout distance and maximum velocity the Voellmy model will work.

However, the Voellmy model will not describe the avalanche flow behind the front, at the tail of the avalanche. Here, measurements show an increase in the friction (a rapid decrease in speed). This effect causes avalanches to elongate and eventually deposit mass. Therefore, the Voellmy model will not predict the deposition behaviour. The Voellmy model has difficulties to predict the runout of small avalanches, which sometimes begin immediately to deposit or “to starve”. Of course, small avalanches can be modelled using higher μ and ξ values, but this is a very ad-hoc approach.

**Debris flow:**

In all RAMMS versions up to Version 1.5.01, an ENO (Essentially Non-Oscillatory) scheme was used to numerically solve the governing differential equations (Christen et al., 2010). However, the numerical solution was implemented on strictly orthogonal grids. This improves computational speed, but introduces numerical instabilities especially in steep terrain. The new Version 1.6.20 uses the same second order ENO scheme, but now on general quadrilateral grid. This new scheme improves numerical stability, but slows the computational speed somewhat. The introduction of this stable ENO scheme allows us to use lower H_cutoff values minimizing mass loss during calculations. The standard value of H_cutoff is 0.000001 m..

The time required to simulate an avalanche or a debris flow is a function of the finite volume grid resolution and the size of the calculation domain. Typically we use 5m resolutions and the simulations require around 10 minutes. We usually perform the initial simulations at 10 m resolution and therefore we have results in 1 or 2 minutes. When we have a solution that we like we might take a look at the problem at 2m resolution.