User Forum of Software BASEMENT
| BASEMENTBasic Simulation Environment for computation of environmental flow and natural hazard simulationLaboratory of Hydraulics, Hydrology and Glaciology (VAW)ETH Zurich |  |
You are not logged in.
- Topics: Active | Unanswered
#1 2021-07-27 16:28:38
- valtre
- User
- Registered: 2021-02-15
- Posts: 2
Lateral coupling: between sequential 1D and 2D domains
Good afternoon.
I have 5 sequential 1D domains in my simulation and now I want to lateral connect the third one with a 2D mesh which should be the floodplain. In the error file there's no error, but the simulation doesn't work. When I run the .bcm there's open a grey window which report: "Generic warning: In Y:\buildtrees\vtk\src\ad8774e85e-85de89847a.clean\Rendering\Core\vtkTexture.cxx. line 32
Error: no override found for 'vtkTexture' ".
This is my command file:
PROJECT {
title = BS
author = V3
date = 10/02/2021
}
DOMAIN {
multiregion = 2
PHYSICAL_PROPERTIES {
gravity = 9.81
viscosity = 1e-06
rho_fluid = 1000
}
BASECHAIN_1D {
region_name = 81-196
GEOMETRY {
type = basement
file = from81to196.bmg
cross_section_order = (24.64295 24.44249 24.24159 24.04814 23.84542 23.65059 23.51563 23.50159 23.45892 23.44 23.43 23.41953 23.41404 23.38951 23.30776 23.20108 23.09145 22.99182 22.91054 22.80883 22.70789 22.58664 22.48145 22.37506 22.26598 22.15415 22.0177 21.91621 21.82842 21.81569 21.814 21.813 21.80114 21.79484 21.79037 21.77704 21.67754 21.56865 21.45028 21.36933 21.26467 21.15538 21.0616 20.97295 20.89675 20.79082 20.69837 20.58941 20.49769 20.38773 20.27613 20.19788 20.182 20.17495 20.15862 20.06162 19.95897 19.81858 19.80827 19.806 19.805 19.79577 19.78648 19.78297 19.68753 19.57957 19.48376 19.38126 19.27612 19.16282 19.03916 18.94501 18.86108 18.75953 18.66854 18.57661 18.55117 18.51137 18.48532 18.38357 18.31354 18.20717 18.10244 18.09418 18.08806 18.07615 17.99019 17.90219 17.77451 17.67825 17.58305 17.5089 17.41275 17.34136 17.23343 17.11245 16.99371 16.82817 16.69246 16.67474 16.65428 16.56937 16.45595 16.33298 16.23289 16.112 15.99967 15.85467 15.49769 15.4108 15.27551 15.15603 15.06668 14.92322 14.80526 14.72582 14.633 14.50618 14.39381 14.30453 14.20251 14.17477 14.14395 14.0538 13.95107 13.79677 13.7436 13.63166 13.53061 13.46252 13.41837 13.39004 13.35757 13.27108 13.265 13.260 13.24943 13.1644 13.05518 12.94875 12.82582 12.71601 12.61265 12.52555 12.42797 12.32036 12.21457 12.09161 11.99878 11.97195 11.94603 11.91966 11.8983 11.89194 11.890)
}
HYDRAULICS {
PARAMETER {
minimum_water_depth = 0.01
riemann_solver = roe
simulation_scheme = explicit
SECTION_COMPUTATION {
type = iteration
precision = 0.01
}
}
FRICTION {
type = manning
default_friction = 0.02857
}
BOUNDARY {
type = hydrograph
file = S_Torbole_hydrograph.dat
string = upstream
slope = 10
}
BOUNDARY {
string = downstream
type = coupling_hqrelation
name = downstream_1
file = SD.txt
}
INITIAL {
type = continue
file = 81-196_initial.dat
q_out = 400
WSE_out = 115.53
}
}
TIMESTEP {
CFL = 0.95
maximum_time_step = 60
initial_time_step = 1.05
total_run_time = 1000000
}
OUTPUT {
output_time_step = 1000
console_time_step = 100
restart_time_step = 1E32
}
}
BASECHAIN_1D {
region_name = 254-308
GEOMETRY {
type = basement
file = from254to308.bmg
cross_section_order = (6.29 6.28431 6.21136 6.08023 6.006020 5.886000 5.87928 5.77783 5.677410 5.573320 5.478000 5.37783 5.27954 5.167840 5.07673 4.97323 4.87716 4.83544 4.82 4.81 4.79496 4.78511 4.73661 4.70169 4.58662 4.47228 4.37106 4.27166 4.1652 4.07808 3.95945 3.87263 3.76812 3.67188 3.56586 3.47766 3.38283 3.26565 3.17135 3.06308 2.96543 2.85307 2.77395 2.674000 2.66411 2.65648 2.64497 2.622800 2.609710 2.60631 2.592270 2.58913 2.57601 2.471410 2.35068 2.259980 2.140420 1.946230 1.838510 1.47554 1.374370 1.232130 1.1338 1.02248 0.946040 0.821910 0.722470 0.642300 0.537950 0.424740 0.411180 0.405 )
}
HYDRAULICS {
PARAMETER {
minimum_water_depth = 0.01
riemann_solver = roe
simulation_scheme = explicit
SECTION_COMPUTATION {
type = iteration
precision = 0.01
}
}
FRICTION {
type = manning
default_friction = 0.02857
}
BOUNDARY {
string = downstream
type = coupling_hqrelation
slope = 1.5
name = downstream_4
file = AT.txt
}
INITIAL {
type = continue
file = 254-308_initial.dat
}
BOUNDARY {
type = coupling_hydrograph
string = upstream
slope = 1.5
name = upstream_4
precision = 0.01
}
}
TIMESTEP {
CFL = 0.95
maximum_time_step = 60
initial_time_step = 1
total_run_time = 1000000
}
OUTPUT {
output_time_step = 1000
console_time_step = 1000
restart_time_step = 1E32
}
}
COUPLINGS {
COUPLING {
type = sequential
upstream_subdomain = 81-196
downstream_subdomain = 196-214
HYDRAULICS {
downstream_interface = upstream_2
upstream_interface = downstream_1
}
}
COUPLING {
type = sequential
upstream_subdomain = 196-214
downstream_subdomain = 214-254
HYDRAULICS {
upstream_interface = downstream_2
downstream_interface = upstream_3
}
}
COUPLING {
type = sequential
upstream_subdomain = 214-254
downstream_subdomain = 254-308
HYDRAULICS {
upstream_interface = downstream_3
downstream_interface = upstream_4
}
}
COUPLING {
type = sequential
upstream_subdomain = 254-308
downstream_subdomain = 308-312
HYDRAULICS {
upstream_interface = downstream_4
downstream_interface = upstream_5
}
}
COUPLING_LATERAL {
basechain = 254-308
baseplane = allagabile_1
connections_automatic = no
connections_filename = lateral_connections_created.dat
side_weir_factor = 0.95
my_weir = 0.58
weir_from_dyke1D = yes
damping_factor = 0.5
connections_index = (1 2)
connections_distance = 30
}
PARAMETER {
max_time_level = 1
}
}
BASECHAIN_1D {
region_name = 196-214
GEOMETRY {
type = basement
file = from196to214.bmg
cross_section_order = (11.889 11.88244 11.88022 11.87523 11.86849 11.77979 11.68846 11.59023 11.46978 11.33601 11.21215 11.11063 11.00521 10.88468 10.76433 10.68062 10.5744 10.47446 10.3832 10.27866 10.18075 10.10884 10.08961 10.085)
}
HYDRAULICS {
PARAMETER {
minimum_water_depth = 0.01
riemann_solver = roe
simulation_scheme = explicit
SECTION_COMPUTATION {
type = iteration
precision = 0.01
}
}
FRICTION {
type = manning
default_friction = 0.02857
}
BOUNDARY {
type = coupling_hydrograph
name = upstream_2
string = upstream
precision = 0.01
slope = 22
}
INITIAL {
type = continue
file = 196-214_initial.dat
}
BOUNDARY {
name = downstream_2
string = downstream
type = coupling_hqrelation
file = DC.txt
}
}
TIMESTEP {
CFL = 0.95
maximum_time_step = 60
initial_time_step = 1
total_run_time = 1000000
}
OUTPUT {
output_time_step = 1000
console_time_step = 1000
restart_time_step = 1E32
}
}
BASECHAIN_1D {
region_name = 214-254
GEOMETRY {
type = basement
file = from214to254.bmg
cross_section_order = (10.084 10.07787 10.06108 9.98199 9.89099 9.79053 9.70642 9.67674 9.597 9.49667 9.40281 9.29833 9.19959 9.1039 9.03407 9.0062 8.99669 8.995340 8.99106 8.9888 8.971590 8.95678 8.84878 8.73571 8.64633 8.55389 8.4554 8.35236 8.27991 8.26862 8.25828 8.23946 8.15917 8.07224 8.05455 8.04833 8.03569 7.95164 7.85796 7.7571 7.64838 7.54343 7.45372 7.31707 7.22784 7.13687 6.99511 6.99277 6.87361 6.76701 6.65712 6.653 6.652 6.64513 6.5681 6.470090 6.39983 6.30281 6.30)
}
HYDRAULICS {
PARAMETER {
minimum_water_depth = 0.01
riemann_solver = roe
simulation_scheme = explicit
SECTION_COMPUTATION {
type = iteration
precision = 0.01
}
}
FRICTION {
type = manning
default_friction = 0.02857
}
BOUNDARY {
name = upstream_3
string = upstream
type = coupling_hydrograph
precision = 0.01
slope = 23
}
INITIAL {
type = continue
file = 214-254_initial.dat
}
BOUNDARY {
name = downstream_3
string = downstream
type = coupling_hqrelation
file = CA.txt
}
}
TIMESTEP {
CFL = 0.95
maximum_time_step = 60
initial_time_step = 1
total_run_time = 1000000
}
OUTPUT {
console_time_step = 1000
output_time_step = 1000
restart_time_step = 1E32
}
}
BASECHAIN_1D {
region_name = 308-312
GEOMETRY {
type = basement
file = from308to312.bmg
cross_section_order = (0.400 0.39526 0.380279 0.294150 0.188330 0.080990 0.018330 0.017 0.016 0)
}
HYDRAULICS {
PARAMETER {
minimum_water_depth = 0.01
riemann_solver = roe
simulation_scheme = explicit
SECTION_COMPUTATION {
type = iteration
precision = 0.01
}
}
FRICTION {
type = manning
default_friction = 0.02857
}
BOUNDARY {
name = upstream_5
string = upstream
type = coupling_hydrograph
slope = 1.5
}
BOUNDARY {
name = downstream_5
string = downstream
type = hqrelation
slope = 1.5
}
INITIAL {
type = continue
file = 308-312_initial.dat
}
}
TIMESTEP {
CFL = 0.95
initial_time_step = 1
maximum_time_step = 60
total_run_time = 1000000
}
OUTPUT {
output_time_step = 1000
console_time_step = 1000
restart_time_step = 1E32
}
}
BASEPLANE_2D {
region_name = allagabile_1
GEOMETRY {
type = 2dm
file = allagabile_arco_v1.2dm
}
HYDRAULICS {
PARAMETER {
minimum_water_depth = 0.01
}
FRICTION {
type = manning
index = (1)
friction = (0.0333)
input_type = index_table
default_friction = 0.0333
wall_friction = off
}
INITIAL {
type = dry
}
}
TIMESTEP {
CFL = 0.95
total_run_time = 400000
}
OUTPUT {
console_time_step = 100
SPECIAL_OUTPUT {
type = element_history
output_time_step = 100
element_values = (depth abs_velocity wse)
element_ids = (82 83 84)
}
SPECIAL_OUTPUT {
type = BASEviz
output_time_step = 10
}
}
}
}Could anyone help me?
Furthermore I need to start the 2d domain with initial type=dry
but the 1D domains need initial condition continue, so the time start simulation are different and this will be an error. How can I do?
Kind regards
Offline
#2 2021-11-23 15:49:30
- Matthias Bürgler
- Developer
- From: Zürich
- Registered: 2019-04-04
- Posts: 150
Re: Lateral coupling: between sequential 1D and 2D domains
Hi
The error message is not related to the coupling-setup of the model, but to the SPECIAL_OUTPUT of type BASEviz. BASEviz seems to be affected by a bug in version 2.8.1. This will be fixed in the next release of version 2.
Best
Matthias
Offline