TELEMAC-2D + SISYPHE

Finally, riverbed change calculation.

Coupling SISYPHE with TELEMAC-2D (v8p1r1).

Prepare 2 'cas' and 2 'cli' files for sis and t2d respectively. If you set a boundary condition as unsteady sediment discharge, it will be recommended to prepare 2 'liq' files for better interpretability.

The setting that bothered me till the end was 'cli'.

＊I've only tried bed-lord, so the following is limited to that. 

**********************************

cli (boundary condition file )

**********************************

First, set the boundary conditions. The user manual (v8p1) contains incorrect expressions . You must enter a non-zero value when setting the time series sediment inflow.

Time-series of sediment discharge

The SISYPHE’s boundary condition file must contain the flags as shown below:

4 5 5 0.0 0.0 0.0 0.0 4 0.0 0.0 0.0 565 1

4 5 5 1.0 0.0 0.0 0.0 4 0.0 0.0 0.0 565 1

Regardless of 1.0 and so on, the time series data in the liq file will take precedence. 

The keyword PRESCRIBED SOLID DISCHARGES must be also included in the steering file,with an arbitrary (but not 0.0) value.

If you set solid discharge =1.0 m3/s (without porosity)  in the 'liq' file and NON COHESIVE BED POROSITY = 0.4, printouts information about solid is 1.66666 m3/s (with porosity) at the time step. In manual,

When the keyword PRESCRIBED SOLID DISCHARGES is used, the mass balance provided in the listing printouts information accounts for the pores = Qb/(1-λ), with λ the porosity.

 

*********************************

liquid file（from the forum）

*********************************

liquid Input file for bed load or suspended load

http://www.opentelemac.org/index.php/community-old/forum4/17-sisyphe/12336-liquid-boundaries-file-syntax-for-multible-suspended-load-con

T QG(1) QG(2) CG(3) CG(1,1) CG(1,2) CG(2,1) CG(2,2) CG(3,1) CG(3,2)

(1) is the number of the boundary

How to IMPOSE a time-series of BED-LOAD discharge ?

http://www.openmascaret.org/index.php/assistance/forum5/17-sisyphe/12417-how-to-impose-a-time-series-of-bed-load-discharge

I managed to run a simulation with Liquid Boundary File for Sisyphe (version V8P0)

As you did, it's necessary to put QG instead of Q2BOR. Remenber that the unit is m3/s (without porosity).

 

*********************************

rigid bed

*********************************

In manual, 

Non-erodible bed everywhere
Replace in noerod.f the line of code: 

CALL OV(’X=Y+C ’,ZR,ZF,ZF,-100.D0,NPOIN)

with

\texttt{CALL OV(’X=Y+C ’,ZR,ZF,ZF,0.D0,NPOIN)}

 But source shows below

CALL OV('X=Y+C   ', X=ZR, Y=ZF, C=-100.D0, DIM1=NPOIN)

So, I fixed it like ↓

CALL OV('X=Y+C   ', X=ZR, Y=ZF, C=0.D0, DIM1=NPOIN)

I'm checking it now.