{import:NavigationMenu}
h1. Introduction
Deltares Systems commitment to quality control and quality assurance has leaded them to develop a formal and extensive procedure to verify the correct working of all of their geotechnical engineering tools. An extensive range of benchmark checks have been developed to check the correct functioning of each tool. During product development these checks are run on a regular basis to verify the improved product. These benchmark checks are provided in the following sections, to allow the users to overview the checking procedure and verify for themselves the correct functioning of *D-FLOW SLIDE*.
The benchmarks for Deltares Systems are subdivided into five separate groups as described below:
* *Group 1 -- Benchmarks from literature (exact solution)*
Simple benchmarks for which an exact analytical result is available from literature.
* *Group 2 -- Benchmarks from literature (approximate solution)*
More complex benchmarks described in literature for which an approximate solution is known.
* *Group 3 -- Benchmarks from spreadsheets*
Benchmarks which test program features using Excel spreadsheets.
* *Group 4 -- Benchmarks generated by the program itself*
Benchmarks for which the reference results are generated using *D-FLOW SLIDE*.
* *Group 5 -- Benchmarks compared with other programs*
Benchmarks for which the results of *D-FLOW SLIDE* are compared with the results of other programs.
As much as software developers would wish they could, it is impossible to prove the correctness of any non-trivial program. Re-calculating all the benchmarks and making sure the results are as they should be will prove to some degree that the program works as it should. Nevertheless there will always be combinations of input values that will cause the program to crash or produce wrong results. Hopefully by using the verification procedure the number of times this occurs will be limited.
The benchmarks will all be described to such detail that reproduction is possible at any time. In some cases, when the geometry is too complex to describe, the input file of the benchmark is needed. The results are presented in text format with each benchmark description.
The input files belonging to the benchmarks can be downloaded from those pages.
h1. Overview of the benchmarks
Legend:
(/) = Results of *D-Flow Slide* and results of the Benchmark are identical.
(x) = Results of *D-Flow Slide* and results of the Benchmark differ.
|| Group || File name || Input file \\
D-Flow Slide (*.fsx) || Input file \\
SLIQ2D 1.1.3.1 (*.slq) || Input file \\
SLIQ2D-Dos (*.sli) || Title || Global (VTV) || Simple (CUR-113) || Detailed (TR) || Advanced (SLIQ2D) ||
| 1 | bm1-1 | | | | Study Case described in "Technisch Rapport Voorland Zettingsvloeiing" | (/) | | (/) | |
| 2 | bm2-1 | | | | Spui dike - hmp 63.9 (location Nieuw Beijerland) | (/) | | (/) | |
| | bm2-2 | | | | Spui dike - hmp 65.0 (between locations Oud Beijerland and Nieuw Beijerland) | (/) | | (/) | |
| | bm2-3 | | | | Spui dike - hmp 67.8 (location Oud Beijerland) | (/) | | (x) ^a^ | |
| 3 | bm3-1 | | | | Global check with traject: step 1 = no, step 3 = yes | | | | |
| | bm3-2 | | | | Global check with traject: step 1 = no, step 3 = no | | | | |
| | bm3-3 | | | | Global check with traject: | | | | |
| 4 | bm4-1 | | | | Test on the level indicator | | | | |
| | bm4-2 | | | | Test on the units | | | | |
| 5 | bm5-1a | | | | Comparison with SLIQ2D - Case LGZM1 | | | | (/) |
| | bm5-1b | | | | Comparison with SLIQ2D - Case LGZM2 | | | | (/) |
| | bm5-1c | | | | Comparison with SLIQ2D - Case LGZM3 | | | | (/) |
| | bm5-1d | | | | Comparison with SLIQ2D - Case LGZM4 | | | | (/) |
| | bm5-1e | | | | Comparison with SLIQ2D - Case SIMPLETA | | | | (/) |
| | bm5-1f | | | | Comparison with SLIQ2D - Case LG1D5N5H | | | | (/) |
| | bm5-1g | | | | Comparison with SLIQ2D - Case HBPZBUI3 | | | | (/) |
| | bm5-2 | | | | Comparison with SLIQ2D-Windows - 2 layers partially saturated with fixed slope angle | | | (/) | |
| | bm5-3 | | | | Comparison with DZettingsVloeiing - Study Case described in "Technisch Rapport Voorland Zettingsvloeiing" | | | (/) | |
^a^ The detailed check with *D-Flow Slide* not succeeds whereas it should succeed acc. to the benchmark.
h1. Group 3: Benchmarks from spreadsheets
This section contains tests that are missing in the other groups, for example for the CUR-113 method.
h1. Group 4: Benchmarks generated by the program itself
This section contains testtests for which the reference results are generated with *D-FLOW SLIDE*.
h2. 4.1 Test on the level indicator
h5. Description
To check that the level indicator is correctly inputted in *D-FLOW SLIDE*, the results of the two following calculations should lead to the same results:
* Benchmark 4-1a uses the same input as benchmark 1-1a, refer to paragraph 1.1 above.
* Benchmark 4-1b uses the same input as benchmark 4-1a except that the Y coordinates of the geometry are shifted by \-20 m so that all points have a negative Y coordinate.
h2. 4.2 Test on the units
h5. Description
In the Units tab of *D-FLOW SLIDE*, it is possible to specify the unit of different parameters:
* Fraction: \-, %, ^o^/~oo~ or ppm
* Length: m, mm, cm, inch, ft or km
* Tiny length: xxxm
* Angle: deg, rad, grad, tan or cot
* Weigth: kN/m{^}3^, N/m{^}3^, lb/in{^}3^, MN/m{^}3^, lb/ft{^}3^
* Pressure: kN/m{^}2^, N/m{^}2^, Pa, kPa, MPa, kN/cm{^}2^, psi
* Permeability: m/s, m/min, m/hr, m/day
To test the correctness of the convertion , 7 benchmarks are created using the following set of units:
|| Benchmark name || Fraction unit || Length unit || Tiny length unit || Angle unit || Weigth unit || Pressure unit || Permeability unit ||
| bm4-2a | - | m | xxxm | deg | kN/m{^}3^ | kN/m{^}2^ | m/s |
| bm4-2b | % | mm | xxxm | rad | N/m{^}3^ | N/m{^}2^ | m/min |
| bm4-2c | ^o^/~oo~ | cm | xxxm | grad | lb/in{^}3^ | Pa | m/hr |
| bm4-2d | ppm | inch | xxxm | tan | MN/m{^}3^ | kPa | m/day |
| bm4-2e | - | ft | xxxm | cot | lb/ft{^}3^ | MPa | m/s |
| bm4-2f | - | km | xxxm | deg | kN/m{^}3^ | kN/cm{^}2^ | m/s |
| bm4-2g | - | m | xxxm | deg | kN/m{^}3^ | psi | m/s |
h5. *D-FLOW SLIDE* inputs
h5. *D-FLOW SLIDE* results
h1. Group 5: Benchmarks compared with other programs
h2. 5.1 Comparison with SLIQ2D (DOS and Windows) - One fully saturated layer with variable slope angle
h5. 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*, using both DOS and Windows versions of this program.
Because of the limitations in the DOS version of *SLIQ2D*, only one fully saturated layer is inputted. The geometry and the material properties for each cases are given in the table below.
The original slope angle is 1:1.25 and is set to "variable" so that the program will search (for each point) for the most unfavorable slope.
|| Benchmark name || Original file name (SLIQ2D-DOS) || Soil type || Slope height (m) || Slope angle || n (%) || n{~}min~ (%) || n ~max~ (%) || Eps ~voldm0~ || m || r || s ~2~ || s ~max~ || v || k ~so~ (kN/m ^2^) || u || Unit weight (kN/m ^3^) || D ~r~ (%) ||
| bm5-1a | LGZM1 | 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 | LGZM2 | 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 | LGZM3 | 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 | LGZM4 | 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 | SIMPLETA | 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 | LG1D5N5H | 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.333 |
| bm5-1g | HBPZBUI3 | 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 |
h5. Benchmark 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. Nota also that the automatic generation of the variable values of the slope angle used by *SLIQ2D-DOS* is not exactly the same as in *SLIQ2D-Windows* or *D-FLOW SLIDE* (which use 2 extra slopes compared to *SLIQ2D-DOS*: 1:1.4 and 1:1.3). The figures below show the critical slope angles, for each calculated points.
Results of *SLIQ2D-Windows* \- Overview of the unstable points:
| bm5-1a: \\ !Embankment SLIQ2D - bm5-1a.bmp|width=400! | bm5-1b: \\ !Embankment SLIQ2D - bm5-1b.bmp|width=400! | bm5-1c: \\ !Embankment SLIQ2D - bm5-1c.bmp|width=400! |
| bm5-1d: \\ !Embankment SLIQ2D - bm5-1d.bmp|width=400! | bm5-1e: \\ !Embankment SLIQ2D - bm5-1e.bmp|width=400! | bm5-1f: \\ !Embankment SLIQ2D - bm5-1f.bmp|width=400! |
| bm5-1g: \\ !Embankment SLIQ2D - bm5-1g.bmp|width=400! |
Results of *SLIQ2D-Windows* \- Values of the critical slope angle for each calculated point:
| bm5-1a: \\ !Results SLIQ2D - bm5-1a.bmp|width=800! | bm5-1b: \\ !Results SLIQ2D - bm5-1b.bmp|width=800! |
| bm5-1c: \\ !Results SLIQ2D - bm5-1c.bmp|width=800! | bm5-1d: \\ !Results SLIQ2D - bm5-1d.bmp|width=800! |
| bm5-1e: \\ !Results SLIQ2D - bm5-1e.bmp|width=800! | bm5-1f: \\ !Results SLIQ2D - bm5-1f.bmp|width=800! |
| bm5-1g: \\ !Results SLIQ2D - bm5-1g.bmp|width=800! | |
h5. *D-FLOW SLIDE* results
Maximum relative error between *SLIQ2D-DOS* and *D-FLOW SLIDE*:
|| Case || Slope angle at which instability occurs (Tan Alpha) || Tan Alpha_p || Ko_p || Normal stress p || Deviatoric stress q ||
| bm5-1a | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| bm5-1b | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| bm5-1c | 0.00 | 0.54 | 0.13 | 0.00 | 0.00 |
| bm5-1d | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| bm5-1e | 0.00 | 0.00 | 0.15 | 0.00 | 0.02 |
| bm5-1f | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| bm5-1g | 0.00 | 0.00 | 0.12 | 0.00 | 0.02 |
Results between *SLIQ2D-Windows* and *D-FLOW SLIDE* are exactly the same for the tangent angle, as show in the table below.
Maximum relative error between *SLIQ2D-Windows* and *D-FLOW SLIDE* for the tangent angle:
|| Case || Tan Alpha ||
| bm5-1a | 0.00 % |
| bm5-1b | 0.00 % |
| bm5-1c | 0.00 % |
| bm5-1d | 0.00 % |
| bm5-1e | 0.00 % |
| bm5-1f | 0.00 % |
| bm5-1g | 0.00 % |
h2. 5.2 Comparison with SLIQ2D-Windows - 2 layers partially saturated with fixed slope angle
h5. 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 Alpha = 0.3).
|| Layer || Height (m) || Soil type || n (%) || n ~min~ (%) || n ~max~ (%) || Eps ~voldm0~ || m || r || s ~2~ || s ~max~ || v || k ~so~ (kN/m ^2^) || u || Unit weight (kN/m ^3^) || D ~r~ (%) ||
| 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 | 5 | 1.23 | 1.35 | 1.25 | 39460 | 1.33 | 6.312 | 20 |
h5. Benchmark results
Results of *SLIQ2D-Windows* for bm5-2 - Overview of the unstable points:
!Embankment SLIQ2D - bm5-2.bmp|width=800!
h5. *D-FLOW SLIDE* results
!Results D-FLOW SLIDE - bm5-2.bmp|width=800!
Results between *SLIQ2D-Windows* and *D-FLOW SLIDE* are exactly the same for all calculated parameters, as shown in the table below.
Maximum relative error between *SLIQ2D-Windows* and *D-FLOW SLIDE*:
|| Parameter || Max. relative error (%) ||
| Tan Alpha_p | 0.00 |
| Ko_p | 0.00 |
| Normal stress p | 0.00 |
| Deviatoric stress q | 0.00 |
h2. 5.3 Comparison with DZettingsVloeiing - Study Case described in "Technisch Rapport Voorland Zettingsvloeiing"
h5. Description
This benchmark is intended to verify the detailed method by comparing *D-FLOW SLIDE* results with those from the program *DZettingsVloeiing*. The same input as in benchmark 1-1 is used.
h5. *DZettingsVloeiing* results
Input and Results of *DZettingsVloeiing* for bm5-3 (= bm1-1):
!Voorland zettingvloeing - Benchmark 1.bmp|width=800!
h5. *D-FLOW SLIDE* results
Results between *DZettingsVloeiing* and *D-FLOW SLIDE* are exactly the same for all calculated parameters, as shown in the table below.
xxx |