Benchmark Problem #6: Three-dimensional subaerial granular slide

This benchmark problem is based on the 3D laboratory experiments of Mohammed and Fritz (2012)

1. Experimental Setup

The landslide tsunami experiments were conducted in the tsunami wave basin at Oregon State University in Corvallis. The landslides are deployed off a plane slope built on one end of the wave basin as shown in Figure 1. The landslide material is deployed in a box measuring , with a volume of  and weight of approximately .

Figure 1. Landslide Tsunami Generator Setup

The plane slope has a slope angle of  corresponding to a slope of 1 vertical to 2 horizontal. Two coordinate systems are defined to characterize the slide motion and wave propagation independently. Both the coordinate systems can be unified through the location of the shoreline, which is dependent on the water depth. The wave gauge locations in the basin with reference to the toe of the slope is shown in Figure 2.

Figure 2. Wave gauge location in the wave basin relative to the wave basin.

The slope coordinate system  and the wave gauge coordinate system  is shown in Figure 3  for a reference water depth of . The locations of the wave gauges are with reference to the toe of the slope shown in Figure 2. The toe and the shoreline can be related to each other knowing the slope angle of the hill slope as shown in Figure 3. Similarly, the shoreline location with reference to the slope coordinate system can be related to the water depth as shown in Figure 3.

Figure 3. Two coordinate systems  and  shown for the slide and wave measurements. The -direction is considered to be perpendicular outside the image.

The wave gauges numbered by indices 1-25 in the wave basin are located at  locations in the tsunami wave basin with reference to the hill slope toe and the wave basin coordinate system as shown in Table 1

Table 1. Wave gauge location with reference to the toe.

WG#

1

24.1

0

2

24.1

2.108

3

14.0

0.0

4

14.0

1.225

5

14.0

3.232

6

14.0

5.374

7

8.5

0.0

8

8.5

0.744

9

8.5

1.962

10

8.5

3.236

11

8.5

4.907

12

5.12

0.0

13

5.12

0.448

14

5.12

1.182

15

5.12

1.965

16

5.12

2.956

17

5.12

5.12

18

5.12

8.868

19

3.9

2.252

20

3.9

3.9

21

3.9

6.755

Runup Wave Gauges

22

0

2.0

23

0

2.6

24

0

3.8

25

0

5.6

 

The location of the wave gauges with reference to the shoreline varies with the water depth. Considering the slope planar hill (1:2) the following increments shown in Table 2 are added to the “X” coordinate system for the water depths.

Table 2. Increment added to the wave gauge map for appropriate water depth.

0.3

0.6

0.6

1.2

1.2

2.4

 

A total of twelve experimental trials are presented along with this document. It should be noted that the slide deposits information are not available for a water depth of . The PX case corresponds to a pressure in the pneumatic pistons of the landslide tsunami generator  with and . For each water depth and P combination, the location of the slide impact with the shoreline in the slope coordinate system changes. The corresponding shoreline location in the slope coordinate system at impact for all the cases is shown in Table 3.  The corresponding slide front velocity can be extracted from the data on slide front velocity as a function of distance  as shown in Figure 4. The time of impact of the slide front with the shoreline is also shown in Table 3. Since the time coordinate of the wave gauges are initialized to the slide motion initiation, the  time allows to fix the moment of slide impact with the water body.

Table 3. LTG Experimental Cases

P145

30

2.663

0.44

60

3.356

0.37

120

2.014

0.3

 

 

P116

30

2.663

0.47

60

3.356

0.39

120

2.014

0.33

 

 

P87

30

2.663

0.49

60

3.356

0.4

120

2.014

0.37

 

 

P58

30

2.663

0.52

60

3.356

0.44

120

2.014

0.4

 

2.   Summary of Coordinate Systems

Wave Gauges

Spatial: Distances of the wave gauges are with reference to the toe of the hill slope. Incremental distances,  shown in Table 2, are added to the “X” coordinate system for the appropriate water depth to obtain the wave gauge locations with reference to the shoreline.

Temporal: The time coordinate in the wave gauge history is with the reference to the initiation of the slide motion. In order to adjust the time coordinate to the moment of slide impact, the time of impact is specified in Table 3.

Slide Velocity, Slide Thickness and Slide Width

Spatial: The slide velocity, thickness and slide width are with reference to the slope coordinate system. In this system, the origin is at the initial location of the slide box along the slide center. The location of impact with the shoreline changes based on the slide location, motion and case as shown in Table 3.

Temporal: The time history of the slide thickness is with reference to the moment of impact of the slide front with the shoreline location.  correspond to the moment of the slide impact with the shoreline.

Slide Deposit

Spatial: The slide deposit data is with reference to the slope toe and the wave basin bottom. The location  correspond to the location of the toe while  correspond to the bottom of the wave basin and the hill slope.

3.   Data Files

All data files referred to below my be found in a zip file LTG_Experiment_Data.zip. The data is stored as Matlab .mat files.

Wave Gauges

File Name: “wg_h_30.mat”, “wg_h_60.mat”, “wg_h_120.mat”

File Description: Time history of water surface elevation at an array of wave gauges in the wave basin.

File Variables: “wg

Variable Description: The wave gauge data is in [cm] units. The wave gauge data are with reference to the wave basin [spatial] and slide initiation [temporal]. Wave gauges 22-25 are measured along the slope since the wave gauges are inclined along the slope.

Variable Dimension:

Dimension Values:

 

 

 

LTG cases corresponding to the 4 P cases shown in Table 3

4

Time Array with

5000

Index of wave gauges as shown in Figure 2.

25

 

The  locations of the wave gauges relative to the toe of the slope in the wave basin is shown in Table 1. The actual distances of the wave gauges  from the shoreline can be calculated based on the water depth and hence the shoreline location from the toe. For ex, for a water depth of 60cm, the shoreline is located at 1.2m from the toe. Hence, adding 1.2m to  in Table 1 gives the actual locations of the offshore wave gauges.

Slide Velocity

File Name: “slide_velocity.mat

File Description: Evolution of slide front velocity with slide front location along the slope in the slope coordinate syste,

File Variables: “v58”,”v58_h”,”v87”,”v87_h”,”v116”,”v116_h”,”v145”,”v145_h”

Variable Description: The velocity data is presented as slide velocity [m/s] versus the slide front [m]. The slide velocity is presented in the slope following coordinate system [spatial]. There is no time reference here. The velocity of the LTG slide box and the released granular landslide is included in this data.

Variable Dimension:

Dimension Values:

 

 

 

Rows indicate

1: slide front location

2: slide front velocity

2

# points in the evolution data

Variable

 

vXXX” contains the slide front velocity along the slide front location for the case corresponding to “XXX = P145,P116,P87,P58”. “vXXX” contains the velocity information corresponding to slide volume of 0.756 m3. “vXXX_h” contains the velocity information corresponding to slide volume of 0.378 m3. The information in the variables vXXX and vXXX_h is shown in Figure 4.

Figure 4. Slide front velocity along the slope with the slide front location.

Slide Thickness

File Name: “Slide_Thick_F.mat

File Description: Maximum slide thickness along the slide center along the slide propagation on the slope for slide volume of 0.756.

File Variables: “XFM_P145”,”YFM_P145”, “XFM_P116”,”YFM_P116”, “XFM_P87”,”YFM_P87”, “XFM_P58”,”YFM_P58”

Variables Description:”X” represents the location on the slope [m] and “Y” represents the maximum slide thickness at the “X” location [m]. The “X” is in the slope coordinate system [spatial]. There is no time reference.

 

File Name: “Slide_Thick_h.mat

File Description: Maximum slide thickness along the slide center along the slide propagation on the slope for slide volume of 0.378.

File Variables: “XHM_P145”,”YHM_P145”, “XHM_P116”,”YHM_P116”, “XHM_P87”,”YHM_P87”, “XHM_P58”,”YHM_P58”

Variables Description:”X” represents the location on the slope [m] and “Y” represents the maximum slide thickness [m] at the “X” location. The “X” is in the slope coordinate system [spatial]. There is no time reference.

The maximum slide thickness along the slide propagation on the slope in the slope following coordinate is shown in Figure 5.

Figure 5. Maximum slide thickness along the slide slope location for slide volumes of  and .

File Name: “Thick_PX_FM.mat” [ 4 files with varying values of X]

File Description: Time history of slide thickness along the slide center for the 4 cases corresponding to “X=P145, P116, P87, P58” for slide volume of 0.786 .

File Variables: “x”,”t”,”S

Variables Description:

                  “x” : slide propagation along the slope [m]

                  “t” : time of measurement adjusted to time  of impact (t=0) [s]

                  “S” : slide thickness with dimensions  (i.e length of “x” times length of “t”)[m]

The “x” coordinate is in the slope coordinate system [spatial]. The time is with reference to the moment of impact ( is the moment of slide impact with the shoreline) [temporal]

 

File Name: “Thick_PX_HM.mat” [ 4 files with varying values of X]

File Description: Time history of slide thickness along the slide center for the 4 cases corresponding to “X=P145, P116, P87, P58” for slide volume of 0.378.

File Variables: “x”,”t”,”S

Variables Description:

                  “x” : slide propagation along the slope [m]

                  “t” : time of measurement adjusted to time of impact (t=0) [s]

                  “S” : slide thickness with dimensions  (i.e length of “x” times length of “t”) [m]

The “x” coordinate is in the slope coordinate system [spatial]. The time is with reference to the moment of impact ( is the moment of slide impact with the shoreline) [temporal]

Slide Width

File Name: “Slide_Width.mat” [ 1 file]

File Description: Edge of maximum slide spread along the slope in the lateral direction. The maximum spread is the same for all the “P” cases and volumes.

File Variables: “X”,”Y

Variables Description:

                  “X” : slide propagation along the slope in slope coordinate system [m]

                  “Y” : Lateral extent of the slide on the slope [m]

The slide motion is with reference to the slope coordinate system [spatial].

The slide width in the file is shown in Figure 6.

Figure 6. Extent of the maximum slide spread on the slope in the lateral direction

Underwater Slide Deposit

File Name: “Deps_PX_h60_FM.mat” [ 4 files with varying values of X]

File Description: Measured underwater slide deposit for the 4 cases corresponding to “X=P145, P116, P87, P58” for slide volume of  at a water depth of

File Variables: “X”,”Y”,”Z”

Variables Description:

                  “X” : horizontal coordinate with reference to the toe [0 value is the toe of the slope] [m]

                  “Y” : lateral coordinate along the slope [0 value is the center of the slide][m]

                  “Z” : two-dimensional array with values of the slide deposits above the basin and slope bottom[m]

The data is with reference to the toe of the slope [spatial]. There is no temporal reference.

 

File Name: “Deps_PX_h120_FM.mat” [ 4 files with varying values of X]

File Description: Measured underwater slide deposit for the 4 cases corresponding to “X=P145, P116, P87, P58” for slide volume of at a water depth of .

File Variables: “X”,”Y”,”Z”

Variables Description:

                  “X” : horizontal coordinate with reference to the toe [0 value is the toe of the slope] [m]

                  “Y” : lateral coordinate along the slope [0 value is the center of the slide] [m]

                  “Z” : two-dimensional array with values of the slide deposits above the basin and slope bottom [m]

The data is with reference to the toe of the slope [spatial]. There is no temporal reference.

 

4. Benchmark Problem:

The benchmark here will consist in simulating the time series of free surface elevations at the 21 labeled wave gauge and 4 runup gauge locations as indicated in Figure 2, along with details of slide geometry and motion. The benchmark is optional, and a model-data comparison table will be constructed for all 12 tests. In order to build up some model-model comparisons prior to the workshop, it would be helpful to concentrate initially on some of the "low pressure" (P58) tests, as these are likely to be the easiest to run.