User Forum of Software BASEMENT

BASEMENT
Basic Simulation Environment for computation of environmental flow and natural hazard simulation
Laboratory of Hydraulics, Hydrology and Glaciology (VAW)
ETH Zurich		Basement_Logo

You are not logged in.

Pages: 1

#1 2024-03-06 10:58:26

Rebecca_Ag
User
Registered: 2024-02-21
Posts: 12

Suspended Load release from a section

Good morning, everyone,

I intend to release suspended sediment, considering fixed bed load, from a specific section (middle of the channel).

I've learned that to achieve this, I need to input sediment data via the "INITIAL" file.

Could you please provide the exact syntax for this file? I've tried the one I found in other forum posts, but it doesn't seem to be working.

Should I input this syntax into an Excel file?

File example:

time 0.0
         C [m3/m3]  Mixture Name
'Qp1a'   0.001      unique
'Qp1b'   0.001      unique
'Qp1c'   0.001      unique
'Qp1d'   0.001      unique
'Qp1e'   0.001      unique
Should the cross-section names be enclosed in quotation marks ""?

May I request the initial file used in the TEST CASE ST_BC_2?

Thank you,

Rebecca

Offline

#2 2024-03-07 00:24:01

Matthias Bürgler
Developer
From: Zürich
Registered: 2019-04-04
Posts: 175

Re: Suspended Load release from a section

Hi Rebecca,

If I understand correctly, you want to provide the initial conditions for a suspended load simulation (initial concentration) via a file.

In this case, the file specified under 'SUSPENDED_LOAD/INITIAL/file' should have the following format:

time: 0                                                                      
                          C		      Mixture
              Qp1	  5.44571E-05          mixture
              Qp1.1       6.97497E-05          mixture        
              Qp1.2       8.88912E-05          mixture        
              Qp1.3       0.000112721          mixture        
              Qp1.4       0.000142225          mixture        
              Qp1.5       0.000178558          mixture        
              Qp1.6       0.000223054          mixture      
              Qp1.7       0.000277248          mixture

The test cases ST_BC_2a or ST_BC_2b including any initial files can be downloaded from the Test Case Section of official BASEMENT website.

I hope this helps.

Best regards,
Matthias

Offline

#3 2024-03-12 11:44:58

Rebecca_Ag
User
Registered: 2024-02-21
Posts: 12

Re: Suspended Load release from a section

Hello Matthias,

Unfortunately it's still not working.
Even if i copy paste my information in the file you sent me or in the initial files of the test case.
I can't figure out what's the problem.

I have a very regular geometry. Rectangular channel exactly as in the test case.

Offline

#4 2024-03-12 23:10:21

Matthias Bürgler
Developer
From: Zürich
Registered: 2019-04-04
Posts: 175

Re: Suspended Load release from a section

Hi Rebecca

Is there an error message before the simulation crashes? And does the model run if you use only the hydrodynamic block?

Regards,
Matthias

Offline

#5 2024-03-13 09:43:34

Rebecca_Ag
User
Registered: 2024-02-21
Posts: 12

Re: Suspended Load release from a section

Dear Matthias,
Thank you for your quick answer.

The hydrodynamic block runs and it reaches the steady state.
Also the morphodynamic block runs but only if i put the input as a BOUNDARY CONDITION FILE (in this case i can set only the duration of the "spill").

Once i set the INITIAL file in the suspended load block it just crashes, no error messages.

p.s i tried to modify the geometry of the test cases with my data but it does not run as well.


Thank you,
Rebecca

Offline

#6 2024-03-13 16:00:39

Rebecca_Ag
User
Registered: 2024-02-21
Posts: 12

Re: Suspended Load release from a section

Sorry still to bother you.
I see in the Geometry file this Syntax.

It's maybe this command that changes since the test case runs, and my model (that i set exactly in the same way) does not?
What are FUNIN VALUES?


// FUNIN values: 0.1	0.05	30 30 , , 30	/
/* generated section -> put it into your bmc file!
   ATTENTION: bmg expects friction definition in the HYDRAULICS > FRICTION block
SECTION_COMPUTATION {
    type = table ST_BC_2_table /
    min_interval = 0.05
    max_interval = 0.1
    table_name   = unused
}
*/
// FUNIN 20: 	Qp1	0	, ,,, , 1 /
// FUNIN 21: 	0	13	 	/
// FUNIN 21: 	0	1	/
// FUNIN 21: 	20	1	/
// FUNIN 21: 	20	13	 	/
CROSS_SECTION {
    name = Qp1
    distance_coord = 0
    friction_coefficients = (30)
    friction_ranges       = (
        (0,20)
    )
    node_coords = (
        (0,13),
        (0,1),
        (20,1),
        (20,13)
    )
}

Thank you,

Rebecca

Offline

#7 2024-03-18 00:18:42

Matthias Bürgler
Developer
From: Zürich
Registered: 2019-04-04
Posts: 175

Re: Suspended Load release from a section

Hi Rebecca

The lines in the geometry file starting with "//" or any lines between "/*" and "*/" are commented out and do affect how the geometry file is read.

From your description I cannot figure out the problem. If you can provide the complete files of a test case that reproduces the problem, I can try to have a look.

Best Regards,
Matthias

Offline

#8 2024-03-18 14:11:26

Rebecca_Ag
User
Registered: 2024-02-21
Posts: 12

Re: Suspended Load release from a section

Good Morning Matthias,

I copy here my Geometry File and the Initial file of the suspended load block.

Beside these two files, the model is exactly the same of the test case.

GEOMETRY : 

//---------------------------------------------------
// BASEMENT file
// generated by BASEview
//---------------------------------------------------
CROSS_SECTION {
	name                     = 5000
	distance_coord           = 0.000000
	left_point_global_coords = (-1152.917297, 2320.247198, 27.000000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,27.000000),
(0.000000,25.000000),
(16.000000,25.000000),
(16.000000,27.000000)
 )
}
CROSS_SECTION {
	name                     = 4500.00
	distance_coord           = 0.500000
	left_point_global_coords = (-700.294937, 2532.692485, 24.500000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,24.500000),
(0.000000,22.500000),
(16.000000,22.500000),
(16.000000,24.500000)
 )
}
CROSS_SECTION {
	name                     = 4000.00
	distance_coord           = 1.000000
	left_point_global_coords = (-247.672578, 2745.137772, 22.000000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,22.000000),
(0.000000,20.000000),
(16.000000,20.000000),
(16.000000,22.000000)
 )
}
CROSS_SECTION {
	name                     = 3500.00
	distance_coord           = 1.500000
	left_point_global_coords = (204.949781, 2957.583059, 19.500000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,19.500000),
(0.000000,17.500000),
(16.000000,17.500000),
(16.000000,19.500000)
 )
}
CROSS_SECTION {
	name                     = 3000.00
	distance_coord           = 2.000000
	left_point_global_coords = (657.572141, 3170.028345, 17.000000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,17.000000),
(0.000000,15.000000),
(16.000000,15.000000),
(16.000000,17.000000)
 )
}
CROSS_SECTION {
	name                     = 2500.00
	distance_coord           = 2.500000
	left_point_global_coords = (1110.194500, 3382.473632, 14.500000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,14.500000),
(0.000000,12.500000),
(16.000000,12.500000),
(16.000000,14.500000)
 )
}
CROSS_SECTION {
	name                     = 2000.00
	distance_coord           = 3.000000
	left_point_global_coords = (1562.816859, 3594.918919, 12.000000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,12.000000),
(0.000000,10.000000),
(16.000000,10.000000),
(16.000000,12.000000)
 )
}
CROSS_SECTION {
	name                     = 1500.00
	distance_coord           = 3.500000
	left_point_global_coords = (2015.439218, 3807.364206, 9.500000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,9.500000),
(0.000000,7.500000),
(16.000000,7.500000),
(16.000000,9.500000)
 )
}
CROSS_SECTION {
	name                     = 1000.00
	distance_coord           = 4.000000
	left_point_global_coords = (2468.061578, 4019.809493, 7.000000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,7.000000),
(0.000000,5.000000),
(16.000000,5.000000),
(16.000000,7.000000)
 )
}
CROSS_SECTION {
	name                     = 500.00
	distance_coord           = 4.500000
	left_point_global_coords = (2920.683937, 4232.254779, 4.500000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,4.500000),
(0.000000,2.500000),
(16.000000,2.500000),
(16.000000,4.500000)
 )
}
CROSS_SECTION {
	name                     = 0
	distance_coord           = 5.000000
	left_point_global_coords = (3373.306296, 4444.700066, 2.000000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,2.000000),
(0.000000,0.000000),
(16.000000,0.000000),
(16.000000,2.000000)
 )
}
//---------------------------------------------------
// BASEMENT file
// generated by BASEview
//---------------------------------------------------
CROSS_SECTION {
	name                     = 5000
	distance_coord           = 0.000000
	left_point_global_coords = (-1152.917297, 2320.247198, 27.000000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,27.000000),
(0.000000,25.000000),
(16.000000,25.000000),
(16.000000,27.000000)
 )
}
CROSS_SECTION {
	name                     = 4500.00
	distance_coord           = 0.500000
	left_point_global_coords = (-700.294937, 2532.692485, 24.500000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,24.500000),
(0.000000,22.500000),
(16.000000,22.500000),
(16.000000,24.500000)
 )
}
CROSS_SECTION {
	name                     = 4000.00
	distance_coord           = 1.000000
	left_point_global_coords = (-247.672578, 2745.137772, 22.000000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,22.000000),
(0.000000,20.000000),
(16.000000,20.000000),
(16.000000,22.000000)
 )
}
CROSS_SECTION {
	name                     = 3500.00
	distance_coord           = 1.500000
	left_point_global_coords = (204.949781, 2957.583059, 19.500000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,19.500000),
(0.000000,17.500000),
(16.000000,17.500000),
(16.000000,19.500000)
 )
}
CROSS_SECTION {
	name                     = 3000.00
	distance_coord           = 2.000000
	left_point_global_coords = (657.572141, 3170.028345, 17.000000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,17.000000),
(0.000000,15.000000),
(16.000000,15.000000),
(16.000000,17.000000)
 )
}
CROSS_SECTION {
	name                     = 2500.00
	distance_coord           = 2.500000
	left_point_global_coords = (1110.194500, 3382.473632, 14.500000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,14.500000),
(0.000000,12.500000),
(16.000000,12.500000),
(16.000000,14.500000)
 )
}
CROSS_SECTION {
	name                     = 2000.00
	distance_coord           = 3.000000
	left_point_global_coords = (1562.816859, 3594.918919, 12.000000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,12.000000),
(0.000000,10.000000),
(16.000000,10.000000),
(16.000000,12.000000)
 )
}
CROSS_SECTION {
	name                     = 1500.00
	distance_coord           = 3.500000
	left_point_global_coords = (2015.439218, 3807.364206, 9.500000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,9.500000),
(0.000000,7.500000),
(16.000000,7.500000),
(16.000000,9.500000)
 )
}
CROSS_SECTION {
	name                     = 1000.00
	distance_coord           = 4.000000
	left_point_global_coords = (2468.061578, 4019.809493, 7.000000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,7.000000),
(0.000000,5.000000),
(16.000000,5.000000),
(16.000000,7.000000)
 )
}
CROSS_SECTION {
	name                     = 500.00
	distance_coord           = 4.500000
	left_point_global_coords = (2920.683937, 4232.254779, 4.500000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,4.500000),
(0.000000,2.500000),
(16.000000,2.500000),
(16.000000,4.500000)
 )
}
CROSS_SECTION {
	name                     = 0
	distance_coord           = 5.000000
	left_point_global_coords = (3373.306296, 4444.700066, 2.000000)
	orientation_angle        = 115.143737
	main_channel_range       = (0.000000, 16.000000)
	friction_ranges          = (
 (0.000000,16.000000) 
)
	friction_coefficients    = ( 40 )
	node_coords              = (
  (0.000000,2.000000),
(0.000000,0.000000),
(16.000000,0.000000),
(16.000000,2.000000)
 )
}


SUSPENDED LOAD INITIAL FILE

time: 0.0
C    Mixture
5000    0.05    uniform
4500.00   0    uniform
4000.00   0    uniform
3500.00   0    uniform
3000.00   0    uniform
2500.00   0    uniform
2000.00   0    uniform
1500.00   0    uniform
1000.00   0    uniform
500.00   0    uniform
0   0    uniform


Thank you for helping me

Offline

#9 2024-03-18 22:36:03

Matthias Bürgler
Developer
From: Zürich
Registered: 2019-04-04
Posts: 175

Re: Suspended Load release from a section

Hi Rebecca

Please provide a complete test case (all necessary files) in a zip folder. Either email them to me or provide a link for download.

Regards,
Matthias

Offline

#10 2024-03-22 09:55:12

Rebecca_Ag
User
Registered: 2024-02-21
Posts: 12

Re: Suspended Load release from a section

Hi Matthias,

can you provide me an e-mail so i can send the zip folder?

Thank you,
Best,

Rebecca

Offline

#11 2024-04-04 10:23:44

Matthias Bürgler
Developer
From: Zürich
Registered: 2019-04-04
Posts: 175

Re: Suspended Load release from a section

Hi Rebecca,

The problem is a missing "end" at the bottom of the suspended load initial file:

time: 0.0
C    Mixture
5000    0.05    uniform
4500.00   0    uniform
4000.00   0    uniform
3500.00   0    uniform
3000.00   0    uniform
2500.00   0    uniform
2000.00   0    uniform
1500.00   0    uniform
1000.00   0    uniform
500.00   0    uniform
0   0    uniform
end # <-- this was missing

I hope this solves the issue.

Best regards,
Matthias

Offline

Pages: 1

Board footer

Powered by FluxBB