5.1 Comparison with SLIQ2D-Windows - Fully saturated slope with 1 layer
Description
The benchmarks in this paragraph are intended to verify the advanced method by comparing D-FLOW SLIDE results with those from the older program SLIQ2D Windows.
In this benchmark, a fully saturated layer is inputted. The geometry and the material properties for each cases are given in the table below.
Case | Soil type | Slope height | Slope angle | n | nmin | nmax | εvoldm0 | m | r | s2 | smax | v | kso | u | γ | Dr |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
bm5-1a | Sand clay | 21.9 | 1:1.25 | 45.5 | 50 | 35 | 0.0025 | 1.7 | 7 | 1.28 | 1.7 | 1 | 50000 | 1 | 8.856 | 30 |
bm5-1b | Sand clay | 21.9 | 1:1.25 | 45.5 | 50 | 35 | 0.0035 | 1.7 | 7 | 1.28 | 1.7 | 1 | 50000 | 1 | 8.856 | 30 |
bm5-1c | Sand clay | 21.9 | 1:1.25 | 45.5 | 50 | 35 | 0.005 | 1.7 | 7 | 1.28 | 1.7 | 1 | 50000 | 1 | 8.856 | 30 |
bm5-1d | Sand clay | 21.9 | 1:1.25 | 45.5 | 50 | 35 | 0.006 | 1.7 | 7 | 1.28 | 1.7 | 1 | 50000 | 1 | 8.856 | 30 |
bm5-1e | Sea sand | 10 | 1:1.25 | 47.4 | 50 | 37 | 0.0092 | 3 | 7 | 1.23 | 1.35 | 1 | 39460 | 1.33 | 8.547 | 20 |
bm5-1f | Sand clay | 20.9 | 1:1.25 | 40.2 | 50 | 35 | 0.0022 | 1.7 | 7 | 1.18 | 1.4 | 1 | 85000 | 1 | 9.717 | 65.33 |
bm5-1g | Sea sand | 22 | 1:1.25 | 45.5 | 50 | 35 | 0.0054 | 2 | 7 | 1.25 | 1.4 | 1.25 | 50000 | 1 | 8.856 | 30 |
SLIQ 2D-Windows results
In SLIQ2D-Windows only the value of the tangent angle is available in the output. The other parameters (such as Ko, the normal stress p and the deviatoric stress q) are not saved and can only be red on the screen, making the comparison difficult. The figures below show the position of the unstable points and the points considered as instable (i.e. stable points surrounded with instable points) as option Fill Holes is active.
Results of SLIQ2D-Windows - Overview of the unstable points for bm5-1a:
Results of SLIQ2D-Windows - Overview of the unstable points for bm5-1b:
Results of SLIQ2D-Windows - Overview of the unstable points for bm5-1c:
Results of SLIQ2D-Windows - Overview of the unstable points for bm5-1d:
Results of SLIQ2D-Windows - Overview of the unstable points for bm5-1e:
Results of SLIQ2D-Windows - Overview of the unstable points for bm5-1f:
Results of SLIQ2D-Windows - Overview of the unstable points for bm5-1g:
D-FLOW SLIDE results
Results between SLIQ2D-Windows and D-FLOW SLIDE are exactly the same for the position of the imaginary surface and for the number of instable points, as show in the tables below.
Comparison between SLIQ2D-Windows and D-FLOW SLIDE for benchmark 5-1a:
| SLIQ2D-Windows | D-Flow Slide | Error |
---|---|---|---|
Imaginary geometry: | |||
Corrected tangent angle [-] | 0.8 | 0.8 | 0.00 % |
Corrected bottom length [m] | 27.375 | 27.375 | 0.00 % |
Corrected height [m] | 21.9 | 21.9 | 0.00 % |
Corrected embankment length [m] | 35.057 | 35.057 | 0.00 % |
Offset top [m] | 0 | 0 | 0.00 % |
Offset bottom [m] | 0 | 0 | 0.00 % |
Imaginary elevation of water [m] | 0 | 0 | 0.00 % |
Points: | |||
Number of instable points | 125 | 125 | 0.00 % |
Number of stable points surrounded with instable points | 4 | 4 | 0.00 % |
Comparison between SLIQ2D-Windows and D-FLOW SLIDE for benchmark 5-1b:
| SLIQ2D-Windows | D-Flow Slide | Error |
---|---|---|---|
Imaginary geometry: | |||
Corrected tangent angle [-] | 0.8 | 0.8 | 0.00 % |
Corrected bottom length [m] | 27.375 | 27.375 | 0.00 % |
Corrected height [m] | 21.9 | 21.9 | 0.00 % |
Corrected embankment length [m] | 35.057 | 35.057 | 0.00 % |
Offset top [m] | 0 | 0 | 0.00 % |
Offset bottom [m] | 0 | 0 | 0.00 % |
Imaginary elevation of water [m] | 0 | 0 | 0.00 % |
Points: | |||
Number of instable points | 177 | 177 | 0.00 % |
Number of stable points surrounded with instable points | 1 | 1 | 0.00 % |
Comparison between SLIQ2D-Windows and D-FLOW SLIDE for benchmark 5-1c:
| SLIQ2D-Windows | D-Flow Slide | Error |
---|---|---|---|
Imaginary geometry: | |||
Corrected tangent angle [-] | 0.8 | 0.8 | 0.00 % |
Corrected bottom length [m] | 27.375 | 27.375 | 0.00 % |
Corrected height [m] | 21.9 | 21.9 | 0.00 % |
Corrected embankment length [m] | 35.057 | 35.057 | 0.00 % |
Offset top [m] | 0 | 0 | 0.00 % |
Offset bottom [m] | 0 | 0 | 0.00 % |
Imaginary elevation of water [m] | 0 | 0 | 0.00 % |
Points: | |||
Number of instable points | 232 | 232 | 0.00 % |
Number of stable points surrounded with instable points | 10 | 10 | 0.00 % |
Comparison between SLIQ2D-Windows and D-FLOW SLIDE for benchmark 5-1d:
| SLIQ2D-Windows | D-Flow Slide | Error |
---|---|---|---|
Imaginary geometry: | |||
Corrected tangent angle [-] | 0.8 | 0.8 | 0.00 % |
Corrected bottom length [m] | 27.375 | 27.375 | 0.00 % |
Corrected height [m] | 21.9 | 21.9 | 0.00 % |
Corrected embankment length [m] | 35.057 | 35.057 | 0.00 % |
Offset top [m] | 0 | 0 | 0.00 % |
Offset bottom [m] | 0 | 0 | 0.00 % |
Imaginary elevation of water [m] | 0 | 0 | 0.00 % |
Points: | |||
Number of instable points | 196 | 196 | 0.00 % |
Number of stable points surrounded with instable points | 80 | 80 | 0.00 % |
Comparison between SLIQ2D-Windows and D-FLOW SLIDE for benchmark 5-1e:
| SLIQ2D-Windows | D-Flow Slide | Error |
---|---|---|---|
Imaginary geometry: | |||
Corrected tangent angle [-] | 0.8 | 0.8 | 0.00 % |
Corrected bottom length [m] | 12.5 | 12.5 | 0.00 % |
Corrected height [m] | 10 | 10 | 0.00 % |
Corrected embankment length [m] | 16.008 | 16.008 | 0.00 % |
Offset top [m] | 0 | 0 | 0.00 % |
Offset bottom [m] | 0 | 0 | 0.00 % |
Imaginary elevation of water [m] | 0 | 0 | 0.00 % |
Points: | |||
Number of instable points | 195 | 195 | 0.00 % |
Number of stable points surrounded with instable points | 18 | 18 | 0.00 % |
Comparison between SLIQ2D-Windows and D-FLOW SLIDE for benchmark 5-1f:
| SLIQ2D-Windows | D-Flow Slide | Error |
---|---|---|---|
Imaginary geometry: | |||
Corrected tangent angle [-] | 0.8 | 0.8 | 0.00 % |
Corrected bottom length [m] | 26.125 | 26.125 | 0.00 % |
Corrected height [m] | 20.9 | 20.9 | 0.00 % |
Corrected embankment length [m] | 33.456 | 33.456 | 0.00 % |
Offset top [m] | 0 | 0 | 0.00 % |
Offset bottom [m] | 0 | 0 | 0.00 % |
Imaginary elevation of water [m] | 0 | 0 | 0.00 % |
Points: | |||
Number of instable points | 188 | 188 | 0.00 % |
Number of stable points surrounded with instable points | 1 | 1 | 0.00 % |
Comparison between SLIQ2D-Windows and D-FLOW SLIDE for benchmark 5-1g:
| SLIQ2D-Windows | D-Flow Slide | Error |
---|---|---|---|
Imaginary geometry: | |||
Corrected tangent angle [-] | 0.8 | 0.8 | 0.00 % |
Corrected bottom length [m] | 27.5 | 27.5 | 0.00 % |
Corrected height [m] | 22 | 22 | 0.00 % |
Corrected embankment length [m] | 35.217 | 35.217 | 0.00 % |
Offset top [m] | 0 | 0 | 0.00 % |
Offset bottom [m] | 0 | 0 | 0.00 % |
Imaginary elevation of water [m] | 0 | 0 | 0.00 % |
Points: | |||
Number of instable points | 208 | 208 | 0.00 % |
Number of stable points surrounded with instable points | 1 | 1 | 0.00 % |
5.2 Comparison with SLIQ2D-Windows - Partially saturated slope with 2 layers
Description
This benchmark is intended to verify the advanced method by comparing D-FLOW SLIDE results with those from the older program SLIQ2D-Windows, for a 21 m height channel composed of 2 layers where the top layer is partially saturated. The slope of the channel is fixed to 1:3.333 (i.e. Tan α = 0.3).
Layer | Height | Soil type | n | nmin | nmax | εvoldm0 | m | r | s2 | smax | v | kso | u | γ | Dr |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Top sand | 12 | Sand clay | 45.5 | 50 | 35 | 0.006 | 1.7 | 7 | 1.28 | 1.7 | 1 | 50000 | 1 | 8.856 | 30 |
Bottom sand | 9 | Sand clay | 47.4 | 50 | 37 | 0.0092 | 3 | 7 | 1.23 | 1.35 | 1.25 | 39460 | 1.33 | 6.312 | 20 |
Benchmark results
Results of SLIQ2D-Windows for bm5-2 - Overview of the unstable points:
D-FLOW SLIDE results
Results between SLIQ2D-Windows and D-FLOW SLIDE are exactly the same for the position of the imaginary surface and for the number of instable points, as show in the tables below.
Comparison between SLIQ2D-Windows and D-FLOW SLIDE for benchmark 5-2:
| SLIQ2D-Windows | D-Flow Slide | Error |
---|---|---|---|
Imaginary geometry: | |||
Corrected tangent angle [-] | 0.301 | 0.301 | 0.00 % |
Corrected bottom length [m] | 81 | 81 | 0.00 % |
Corrected height [m] | 24.387 | 24.387 | 0.00 % |
Corrected embankment length [m] | 84.592 | 84.592 | 0.00 % |
Offset top [m] | 11 | 11 | 0.00 % |
Offset bottom [m] | 0 | 0 | 0.00 % |
Imaginary elevation of water [m] | 3.387 | 3.387 | 0.00 % |
Points: | |||
Number of instable points | 96 | 96 | 0.00 % |
Number of stable points surrounded with instable points | 4 | 4 | 0.00 % |
5.3 Comparison with HMBreach - One sand layer, variable slope with channel depth
Description
This benchmark is intended to verify the advanced method for breaching by comparing D-FLOW SLIDE results with those from the program HMBreach.
Only one material is used with the following properties:
- Friction angle φ = 32°
- Diameter D50 = 200 μm
- Ration D50/D15 = 1.75 which leads to D15 = 115 μm
- Porosity n = 40%
- Density of sand grains ρsand = 2650 kg/m3 which leads to γsand = ρsand × g / 1000 = 26 kN/m3 (with g = 9.81236)
- Density of water ρwater = 1000 kg/m3 which leads to γwater = ρwater × g / 1000 = 9.81 kN/m3 (with g = 9.81236)
- Temperature 15°
The soil profile along the channel is composed of 10 layers:
Layer | Thickness [m] | Slope angle [°] |
---|---|---|
1 | 0.5 | 26 |
2 | 0.5 | 26 |
3 | 1 | 15 |
4 | 2 | 15 |
5 | 2 | 15 |
6 | 2 | 10 |
7 | 2 | 10 |
8 | 2 | 5 |
9 | 5 | 5 |
10 | 5 | 2 |
HMBreach results
In HMBreach, the initialization height is varied manually from 0.10 m until 1.00 m with an interval of 0.10 m. Note: a more accurate computation is obtained with interval 0.01 m.
The value of the sand transport at the top (SZ0) and the toe (SZtoe) of the channel are given in the table below and compared with the values calculated by D-Flow Slide.
The critical initialization height (corresponding to the initialization height for which the ratio SZtoe / SZ0 is higher than10) calculated by HMBreach is expected to be equal to hcrit = 0.60 m (see table below).
D-FLOW SLIDE results
Results of HMBreach and D-FLOW SLIDE are very close when the initialization height (h0) is less than 0.6 m, as shown in the table below. Above 0.6 m, the sand transport at channel toe differs. However, the critical initialization height corresponding to the initialization height for which the ratio SZtoe / SZ0 is higher than 10, is correct.
Two cases are considered with D-FLOW SLIDE :
Case A: the allowable critical height (hallowable) is 0.50 m. As the critical initialization height is higher than the allowable critical height, the slope is safe. So the Advanced check for breaching passes.
Case B: the allowable critical height (hallowable) is 0.60 m. As the critical initialization height is equal to the allowable critical height, it means that the probability that a breach flow slide will occur is high. So the Advanced check for breaching fails.
Benchmark 5-3: Sand transports results for HMBreach and D-FLOW SLIDE:
Initial. height | HMBreach | D-FLOW SLIDE | Relative error | ||||||
---|---|---|---|---|---|---|---|---|---|
h0 [m] | SZtoe [kg/s] | SZ0 [kg/s] | SZtoe / SZ0 [-] | SZtoe [kg/s] | SZ0 [kg/s] | SZtoe / SZ0 [-] | SZtoe [%] | SZ0 [%] | SZtoe / SZ0 [%] |
0.1 | 1.256 | 1.002 | 1.254 | 1.255 | 1.004 | 1.250 | 0.08 | 0.20 | 0.32 |
0.2 | 1.361 | 2.003 | 0.679 | 1.360 | 2.009 | 0.677 | 0.07 | 0.30 | 0.30 |
0.3 | 1.441 | 3.005 | 0.480 | 1.444 | 3.013 | 0.479 | 0.21 | 0.27 | 0.21 |
0.4 | 1.609 | 4.007 | 0.402 | 1.602 | 4.018 | 0.399 | 0.44 | 0.27 | 0.75 |
0.5 | 1.793 | 5.009 | 0.358 | 1.788 | 5.022 | 0.356 | 0.28 | 0.26 | 0.56 |
0.6 | 263.352 | 6.010 | 43.817 | 224.543 | 6.027 | 37.258 | 17.28 | 0.28 | 17.60 |
0.7 | 314.820 | 7.012 | 44.898 | 275.447 | 7.031 | 39.175 | 14.29 | 0.27 | 14.61 |
0.8 | 343.467 | 8.014 | 42.860 | 299.110 | 8.036 | 37.223 | 14.83 | 0.27 | 15.14 |
0.9 | 364.280 | 9.015 | 40.407 | 317.751 | 9.040 | 35.149 | 14.64 | 0.28 | 14.96 |
1 | 378.509 | 10.017 | 37.787 | 331.657 | 10.045 | 33.019 | 14.13 | 0.28 | 14.44 |
Benchmark 5-3: Criterium results for HMBreach and D-FLOW SLIDE:
Case | HMBreach | D-FLOW SLIDE | Error [%] | |
---|---|---|---|---|
Case A: hallowable = 0.50 m | Critical initialization height [m] | 0.60 | 0.60 | 0.00 |
Succeeded | Pass | Pass | OK | |
Case B: hallowable = 0.60 m | Critical initialization height [m] | 0.60 | 0.60 | 0.00 |
Succeeded | Fail | Fail | OK |
With HMbreach can be computed that the critical initiation height in this case is 0.56 m. The graph below shows the slope geometry (red line) and the sand transport along the sand slope (blue line).
The increasing sand transport rate is an indication that the erosion is sufficient to maintain a self-accelerating turbidity current. The initial breach can retrogress and grow in height and finally erode the foreshore completely (Breach flow slide).
In D Flow slide the critical height is 0.6 m