[c1] Add default debris-flow friction parameters to config file

debrisframe/c1Tif/c1TifCfg.ini

@@ -48,85 +48,81 @@ timeDependentRelease = True
 # will include the points located on the polygon line)
 thresholdPointInRel = 0.01
 
+#++++++++++++Time stepping parameters
+# End time [s]
+tEnd = 600
+
 #++++++++++++Friction model
 # add the friction using an explicit formulation (1)
 # or use an implicit method (0)
 explicitFriction = 0
 # friction type. Available options are:
-# samosAT - standard setup
-# samosATSmall - setup for release volumes < 25k m³
-# samosATMedium - setup for release volumes < 60k m³
-# samosATAuto - chooses samosAT, samosATSmall or samosATMedium based on release volume and volClass
 # Coulomb
 # Voellmy
-# VoellmyMinShear
 # spatialVoellmy
-# CoulombMinShear
-# wetsnow
 ## Rheological models for debris-flow simulations
 # OBrienAndJulien
 # HerschelAndBulkley
-# Ostwald
 # Bingham
-# Newton
 # Please note that each type has their own/separate parameters:
 # https://docs.avaframe.org/en/latest/theoryCom1DFA.html#friction-model
-frictModel = Bingham
+frictModel = Voellmy
 #+++++++++++++Voellmy friction model
 muvoellmy = 0.155
-xsivoellmy = 4000.
+xsivoellmy = 200
 #+++++++++++++Coulomb friction model
-mucoulomb = 0.155
+mucoulomb = 0.47
 #+++++++++DEBRIS FLOW properties
 # density of sediment [kg/m³]
 rhoSediment = 2650
 # mean sediment size [m]
 ## TODO  adjusting
-sizeSediment = 0.002
+sizeSediment = 0.001
 # maximum concentration of sediment particles [m³/m³]
 cvMaxSediment = 0.615
 # volumetric sediment concentration [m³/m³]
-## TODO adjusting
-cvSediment = 0.5
+cvSediment = 0.4
 #++++++++++++++Parameter Rheological models
+# Basically, there are two options for defining the yield shear stress, tauy,
+# and the dynamic viscosity, eta:
+# - set the values directly (default) or
+# - compute them using the alpha and beta parameters (see documentation).
+# Set the eta and tauy values to 0 if you want to use the option
+# that computes eta and tauy using the alpha and beta parameters.
 #++++++++++++O´Brien and Julien
 # viscosity parameter
-## TODO adjusting, default value O´Brien and Julien (1985)
-alpha1EtaObrienAndJulien = 0.650
-## TODO adjusting, default value O´Brien and Julien (1985)
-beta1EtaObrienAndJulien = 16.81
+etaObrienAndJulien = 25
+# set etaObrienAndJulien to 0 to compute eta from alpha1EtaObrienAndJulien and beta1EtaObrienAndJulien
+alpha1EtaObrienAndJulien = 0.14
+beta1EtaObrienAndJulien = 12.96
 # yield shear stress parameter
-## TODO adjusting, default value O´Brien and Julien (1985)
-alpha2TauyObrienAndJulien = 0.00886
- ## TODO adjusting, default value O´Brien and Julien (1985)
-beta2TauyObrienAndJulien = 13.11
+tauyObrienAndJulien = 200
+# set tauyObrienAndJulien to 0 to compute tauy from alpha2TauyObrienAndJulien and beta2TauyObrienAndJulien
+alpha2TauyObrienAndJulien = 0.2511
+beta2TauyObrienAndJulien = 16.7
 # parameter, default value Takahashi (1978)
 alphaObrienAndJulien = 0.01
 #++++++++++++Herschel and Bulkley
-# viscosity parameter
-## TODO adjusting, default value O´Brien and Julien (1985)
-alpha1EtaHerschelAndBulkley = 0.650
- ## TODO adjusting, default value O´Brien and Julien (1985)
-beta1EtaHerschelAndBulkley = 16.81
 # yield shear stress parameter
- ## TODO adjusting, default value O´Brien and Julien (1985)
-alpha2TauyHerschelAndBulkley = 0.00886
- ## TODO adjusting, default value O´Brien and Julien (1985)
-beta2TauyHerschelAndBulkley = 13.11
+tauyHerschelAndBulkley = 200
+# set tauyHerschelAndBulkley to 0 to compute tauy from alpha2TauyHerschelAndBulkley and beta2TauyHerschelAndBulkley
+alpha2TauyHerschelAndBulkley = 0.2511
+beta2TauyHerschelAndBulkley = 16.7
+# consistency factor K
+kHerschelAndBulkley = 25
 # flow exponent
-## TODO adjusting
-nHerschelAndBulkley = 1.5
+nHerschelAndBulkley = 2
 #++++++++++++Bingham
 # viscosity parameter
-## TODO adjusting, default value O´Brien and Julien (1985)
-alpha1EtaBingham = 0.650
- ## TODO adjusting, default value O´Brien and Julien (1985)
-beta1EtaBingham = 16.81
+etaBingham = 200
+# set etaBingham to 0 to compute eta from alpha1EtaBingham and beta1EtaBingham
+alpha1EtaBingham = 0.2511
+beta1EtaBingham = 16.7
 # yield shear stress parameter
- ## TODO adjusting, default value O´Brien and Julien (1985)
-alpha2TauyBingham = 0.00886
- ## TODO adjusting, default value O´Brien and Julien (1985)
-beta2TauyBingham = 13.11
+tauyBingham = 2000
+# set tauyBingham to 0 to compute tauy from alpha2TauyBingham and beta2TauyBingham
+alpha2TauyBingham = 1.613
+beta2TauyBingham = 17.81
 
 #++++++++++++ Entrainment Erosion Energy
 # Used to determine speed loss via energy loss due to entrained mass

debrisframe/data/debrisParabolaChanneled/Inputs/REL/releaseParabolaChanneled.prj

@@ -0,0 +1 @@
+PROJCS["ETRS_1989_Austria_Lambert",GEOGCS["GCS_ETRS_1989",DATUM["D_ETRS_1989",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",400000.0],PARAMETER["False_Northing",400000.0],PARAMETER["Central_Meridian",13.3333333333333],PARAMETER["Standard_Parallel_1",49.0],PARAMETER["Standard_Parallel_2",46.0],PARAMETER["Latitude_Of_Origin",47.5],UNIT["Meter",1.0]]

debrisframe/data/debrisParabolaChanneled/Inputs/REL/releaseParabolaChanneledTimeDep.csv

@@ -0,0 +1,13 @@
+timestep,thickness,velocity
+0,0.965,3.7
+20,7.33,7.4
+40,10.328,8.3
+60,9.295,8
+80,8.64,7.7
+100,7.29,7.3
+120,5.979,6.9
+140,5.148,6.5
+160,3.805,6
+180,2.718,5.3
+200,1.819,4.5
+220,0.544,3.1

debrisframe/debrisframeCfg.ini

@@ -1,11 +1,11 @@
-### Config File - This file contains the main settings for Avaframe
-# This file is part of Avaframe.
+### Config File - This file contains the main settings for Debrisframe
+# This file is part of Debrisframe.
 
 
 [MAIN]
 # Path to avalanche directory
 #TODO: change avalancheDir to debrisDir
-avalancheDir = ../../AvaFrame/avaframe/data/avaParabolaTimeDep
+avalancheDir = data/debrisGenericTopoTimeDep
 
 # number of CPU cores to use for the computation of com1DFA
 # possible values are: