"Exception" problem with sediment transport

nicolagroff · 2022-03-29 16:08:22

Hi!

I'm running a simulation with sediment transport with a gate and a weir condition.
I got this error:

EXCEPTION -> layer top and layer bottom are identical (down_hyd: 1202.63; top_hyd:
1202.63)...this is not allowed., in file: ..\src\BASEplane\BPSedInput.cpp on
line 467, Version: 2.8 R5771, compiled at: May 31 2018, 15:28:00

What's wrong??

This is the morfology block:

MORPHOLOGY {
PARAMETER {
porosity = 35
density = 2650
control_volume_type = constant
control_volume_thickness = (0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1)
control_volume_thickness_index = (1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 25 26 30 99)
}
BEDMATERIAL {
GRAIN_CLASS {
diameters = (61)
}
MIXTURE {
name = M1
volume_fraction = (100)
}
SOIL_DEF {
name = S0
LAYER {
mixture = M1
bottom_elevation = -1
}
}
SOIL_DEF {
name = S1
LAYER {
mixture = M1
bottom_elevation = -1
}
}
SOIL_DEF {
name = S2
LAYER {
mixture = M1
bottom_elevation = -1
}
}
SOIL_DEF {
name = S3
LAYER {
mixture = M1
bottom_elevation = -0.5
}
}
SOIL_DEF {
name = S4
LAYER {
mixture = M1
bottom_elevation = -1
}
}
SOIL_DEF {
name = S5
LAYER {
mixture = M1
bottom_elevation = -1
}
}
SOIL_DEF {
name = S6
LAYER {
mixture = M1
bottom_elevation = -0.5
}
}
SOIL_DEF {
name = S7
}
SOIL_DEF {
name = S8
LAYER {
mixture = M1
bottom_elevation = -0.5
}
}
SOIL_DEF {
name = S9
LAYER {
mixture = M1
bottom_elevation = 0
}
}
SOIL_DEF {
name = S10
LAYER {
mixture = M1
bottom_elevation = -1
}
}
SOIL_DEF {
name = S11
}
SOIL_DEF {
name = S12
}
SOIL_DEF {
name = S13
LAYER {
mixture = M1
bottom_elevation = 0
}
}
SOIL_DEF {
name = S14
LAYER {
mixture = M1
bottom_elevation = 0
}
}
SOIL_DEF {
name = S15
}
SOIL_DEF {
name = S16
LAYER {
mixture = M1
bottom_elevation = -1
}
}
SOIL_DEF {
name = S25
LAYER {
mixture = M1
bottom_elevation = 0
}
}
SOIL_DEF {
name = S26
LAYER {
mixture = M1
bottom_elevation = 0
}
}
SOIL_DEF {
name = S30
LAYER {
mixture = M1
bottom_elevation = 0
}
}
SOIL_DEF {
name = S99
}
SOIL_ASSIGNMENT {
index = (0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 25 26 30 99)
soil = (S0 S1 S2 S3 S4 S5 S6 S7 S8 S9 S10 S11 S12 S13 S14 S15 S16 S25 S26 S30 S99)
type = index_table
}
}
INITIAL {
type = initial_mesh
}
BEDLOAD {
FORMULA {
bedload_formula = smartjaeggi
d30 = 22
d90 = 293
bedload_factor = 0.25
theta_critical = (0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05)
theta_critical_index = (1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 25 26 30 99)
}
BOUNDARY {
name = SEZ_IN
type = transport_capacity
string_name = SEZ_IN
fraction_boundary = 1
mixture = M1
factor = 0.25
}
BOUNDARY {
name = SEZ_OUT
type = IODown
string_name = SEZ_OUT
fraction_boundary = 0.9
}
PARAMETER {
}
DIRECTION {
lateral_transport_type = lateral_bed_slope
lateral_index = (5 6 7 9)
}
INNER_BOUNDARY {
type = open
string_name1 = PRESA_M
string_name2 = PRESA_V
}
}
GRAVITATIONAL_TRANSPORT {
index = (0 1 2 3 4 5 6 7 8 12 13 14 15 16 24 25 26)
angle_failure_dry = (40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40)
angle_failure_wetted = (30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30)
angle_failure_deposited = (25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25)
gravity_transport_on_cells = partially wetted
}
}
}
}

Thank you!

Nicola

Lukas Vonwiller · 2022-03-29 17:20:05

Hi Nicola

Does it work if you remove the Layers where you have a bottom_elevation = 0?

e.g. instead of

SOIL_DEF {
name = S9
LAYER {
mixture = M1
bottom_elevation = 0
}
}

try

SOIL_DEF {
name = S9
}

This results in a fixed bed for soil S9.
Hope this helps!
cheers, Lukas

nicolagroff · 2022-03-30 09:27:30

Hi Lukas,

It works!

Thank you very much!

Have a nice day

Nicola

User Forum of Software BASEMENT

#1 2022-03-29 16:08:22

"Exception" problem with sediment transport

#2 2022-03-29 17:20:05

Re: "Exception" problem with sediment transport

#3 2022-03-30 09:27:30

Re: "Exception" problem with sediment transport

Board footer