User Forum of Software BASEMENT

Yang He

User

Registered: 2024-06-02

Posts: 8

Issue with Simulation Stopping Unexpectedly in Fixed Bed River Model

Hi!

I have set up a model to simulate a fixed bed river with one tributary inputting suspended load into the main river. The model operates under "continue" conditions and is based on a model that has achieved hydraulic steady state.

I have included the Morphology and Suspended_load components in the model. However, when I attempt to run the simulation, it stops abruptly without any error notifications. Below is the console output information:

"BASEMENT v4.1.0 (v4.1.0)
BASEMD omp binary
-> Domain: Input file has been successfully parsed and passed preliminary
   validation
-> *** Start Simulation ***
-> OpenMP: Computation is executed in parallel using 8 threads
-> Domain Initialized
-> initialize domain done
-> initialize networkcommunication done
-> Initialization of SubDomain 'Inn'
-> TopoBMG: Input file has 405 warnings and 0 errors
-> Grid: building grid using 412 cross sections
-> Grid: ...done.
-> Inner weir: using constant weir level
-> Inner weir: using constant weir width
-> Domain: Starting time loop
"
After this, the button below changes from "Abort" to "Run" automatically, indicating that the model stops the simulation.

When I delete the Morphology and Suspended_load components, the model runs successfully. Below are the settings for the Morphology component for your reference:

"MORPHOLOGY": {
            "BEDMATERIAL": {
              "GRAIN_CLASS": {
                "diameters": [
                  35.0
                ]
              },
              "MIXTURE": [
                {
                  "name": "unique",
                  "volume_fraction": [
                    100.0
                  ]
                }
              ],
              "SOIL_ASSIGNMENT": {
                "index": [
                  1
                ],
                "soil": [
                  "fixed"
                ],
                "type": "index_table"
              },
              "SOIL_DEF": [
                {
                  "name": "fixed",
                  "tau_erosion_start": -1.0,
                  "LAYER": [
                    {
                      "bottom_elevation": 0.0,
                      "mixture": "unique"
                    }
                  ]
                }
              ]
            },
            "INITIAL": {
              "type": "initial_mesh"
            },
            "PARAMETER": {
              "density": 2650.0,
              "porosity": 37.0,
              "control_volume_thickness": 0.02,
              "control_volume_type": "constant",
              "create_new_layers": "off",
              "max_dz_table": 0.1,
              "min_theta_critic": 0.02,
              "wetting_fraction": 0.1
            },
            "SUSPENDED_LOAD": {
              "BOUNDARY": [
                {
                  "string": "upstream",
                  "type": "out_down"
                },
                {
                  "string": "downstream",
                  "type": "out_down"
                }
              ],
              "INITIAL": {
                "concentration": 0.0,
                "type": "global_value",
                "mixture": "unique"
              },
              "PARAMETER": {
                "scheme": "mdpm",
                "sediment_exchange": "off",
                "diffusion_factor": 0.0,
                "diffusion_factor_dynamic": 1,
                "sink_velocity_type": "vanrijn",
                "local_slope": "on"
              },
              "SOURCE": {
                "EXTERNAL_SOURCE": [
                  {
                    "cross_section": "cs_68",
                    "type": "sediment_discharge",
                    "file": "ssc.txt",
                    "name": "Spoel_E_flooding",
                    "mixture": "unique"
                  }
                ]
              }
            }
          },

Could you please help me understand why the simulation stops when including these components and how I might resolve this issue?

Thank you for your assistance.

Best regards,
Yang