open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everyone,

a colleague ask me to post the following information about cohesive sediment in 3D.

Regards,

Constantin

---

"Just in case anyone wishes to use cohesive sediment in 3D and might have crushed planes, in flused subroutine (that calculates fluxes at the bottom)
there are the following lines, that do not give the correct results:

00213 IF(IPBOT%I(I).NE.NPLAN-1) THEN
00214 ! DEPOSITION ON THE FIRST FREE PLANE WITH LOCAL VELOCITY
00215 I3D=I+IPBOT%I(I)*NPOIN2
00216 FLUDPT(I) = WC(I3D)*MAX(1.D0-TOB(I)/MAX(TOCD,1.D-6),0.D0)
00217 ELSE

This is not as intended, the brakes around TOB(I)/MAX(TOCD,1.D-6) are missing it should be as the same formula as in line 00224 for non-crushed planes (which I guess most people will have), just with I3D as argument for setting velocity :
FLUDPT(I) = WC(I3D)*MAX(1.D0-(TOB(I)/MAX(TOCD,1.D-6)),0.D0)

Hello Constantin,

Thanks for this information.

Does this missing brakes really change something ?
In fact with and withtout brakes the calculation is :
(See Calcul_T3D.PNG)



For (1) and (2) the max is still applied for 1.D0 – TOB(I)/MAX(TOCD,1.D-6).
(3) presents a problem, but this is not the case.

In case I'm wrong could you explain me ?

Best regards,
P.N