User Forum of Software BASEMENT

AlessandroCasalicchio

User

Registered: 2020-04-17

Posts: 8

Reduce the duration of the simulation

Hi

I'm doing a 2D morphological-hydraulical simulation but I have a problem with the too long duration of the simulation (I have already increase the max area of the mesh). 

To reduce the duration, first of all I tried to not apply the morphological block to some areas far from the riverbed (i didn't assign the index in the soil assignment block) but the program don't let me do this, because ask me to assign a soil to each index.

So, I ask you how can reduce the duration, maybe some of this attempts:

1) maybe I can increase the thickness of the control volume (parameter block)

2) maybe I can reduce the bottom elevation of some areas

Maybe there are some other things I can change that can help the program to do the simulation

Thank you

AlessandroCasalicchio

User

Registered: 2020-04-17

Posts: 8

Re: Reduce the duration of the simulation

This is the code:

PROJECT {
	title  = ProvaBasement
	author = AC
	date   = 2020
}
DOMAIN {
	PARALLEL {
		number_threads = 8
	}
	BASEPLANE_2D {
		region_name = Po
		GEOMETRY {
			type = 2dm
			file = prog_29_maggio.2dm
			STRINGDEF {
				name               = inflow
				node_ids           = (5 5996 6080 1445 1778 906 1772 2126 1449 2124 1775 5850 725 2423 2049 2300 1713 2304 2064 5737 905 2120 1770 2127 1450 2119 1771 5739 4 1716 1406 5721 724 5732 904 5719 3 2048 1714 2066 1407 2069 1715 1998 723 2052 1709 2059 903 2058 1448 2302 2 5992 5989 5991 1718 6077 5993 6075 1408 5981 5974 5977 5971 5979 5976 5982 1)
				upstream_direction = right
			}
			STRINGDEF {
				name               = outflow
				node_ids           = (17 6758 6701 6753 732 6719 912 6711 18 6704 913 6689 733 6631 6547 6635 19 6579 734 6626 6540 6629 20 6581 2277 6583 1660 6553 6554 735 6551 1458 6592 914 6538 1457 6566 21 6655 6569 1661 6574 6658 736 6665 6575 915 6683 6577 6670 22)
				upstream_direction = right
			}
			STRINGDEF {
				name               = Inflow_sed
				node_ids           = (5 5996 6080 1445 1778 906 1772 2126 1449 2124 1775 5850 725 2423 2049 2300 1713 2304 2064 5737 905 2120 1770 2127 1450 2119 1771 5739 4 1716 1406 5721 724 5732 904 5719 3 2048 1714 2066 1407 2069 1715 1998 723 2052 1709 2059 903 2058 1448 2302 2 5992 5989 5991 1718 6077 5993 6075 1408 5981 5974 5977 5971 5979 5976 5982 1)
				upstream_direction = right
			}
		}
		HYDRAULICS {
			FRICTION {
				type             = manning
				default_friction = 0.028
				wall_friction    = off
				input_type       = index_table
				index            = (1 2 3 4 5 6 7 )
				friction         = (0.028 0.03 0.06 0.1 0.028 0.03 0.06)
			}
			PARAMETER {
				riemann_solver       = hllc
				dynamic_depth_solver = off
			}
			TURBULENCE_MODEL {
				type                = algebraic
				kinematic_viscosity = 0.000001307
			}
			INITIAL {
				type = dry
			}
			BOUNDARY {
				name        = inflow
				string_name = inflow
				type        = hydrograph
				file        = Inflow_T200.txt
				slope       = 0.1
			}
			BOUNDARY {
				name        = outflow
				string_name = outflow
				type        = zhydrograph
				file        = zhyd.txt
			}
		}
		TIMESTEP {
			total_run_time    = 151200
			minimum_time_step = 120
			CFL               = 1.0
		}
		MORPHOLOGY {
			PARAMETER {
				porosity                       = 40
				density                        = 2650
				control_volume_type            = constant
				control_volume_thickness       = (0.4 0.4 0.4)
				control_volume_thickness_index = (1 2 3 )
			}
			INITIAL {
				type = initial_mesh
			}
			BEDMATERIAL {
				GRAIN_CLASS {
					diameters = (5 25 100)
				}
				MIXTURE {
					name            = mixture_alveo
					volume_fraction = (18 25 57 )
				}
				SOIL_DEF {
					name = Alveo
					LAYER {
						bottom_elevation = -5
						mixture          = mixture_alveo
					}
				}
				SOIL_DEF {
					name = Traversa
					LAYER {
						bottom_elevation = -0.0000000001
						mixture          = mixture_alveo
					}
				}
				SOIL_ASSIGNMENT {
					type  = index_table
					index = (1 2 3 4 5 6 7 )
					soil  = (Alveo Alveo Alveo Traversa Traversa Traversa Traversa)
				}
			}
			BEDLOAD {
				FORMULA {
					bedload_formula      = mpmh
					bedload_factor       = 1
					theta_critical_index = (1 2 3 )
					theta_critical       = (0.047 0.047 0.047 )
				}
				BOUNDARY {
					type        = sediment_discharge
					string_name = Inflow_sed
					mixture     = mixture_alveo
					file        = B_C_sed.txt
				}
				BOUNDARY {
					type        = IODown
					string_name = outflow
				}
				DIRECTION {
					lateral_transport_type = lateral_bed_slope
					lateral_index          = (1 )
				}
				PARAMETER {
				}
			}
		}
		OUTPUT {
			console_time_step = 120
			restart_time_step = 1E32
			SPECIAL_OUTPUT {
				type             = node_centered
				output_time_step = 3600
				values           = (wse  depth  velocity  deltaz  z_node )
				format           = sms
				threshold_wse    = -1000000
			}
			SPECIAL_OUTPUT {
				type             = balance
				balance_values   = (sediment  timestep )
				output_time_step = 3600
			}
		}
	}
}

Lukas Vonwiller

Developer

Registered: 2014-09-04

Posts: 127

Re: Reduce the duration of the simulation

Hi Alessandro

Have a look at the 'morph_cycle' in the TIMESTEP block!

Hope this helps!

AlessandroCasalicchio

User

Registered: 2020-04-17

Posts: 8

Re: Reduce the duration of the simulation

Ok thank you. I'll try