open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

in v6p0 coupled sediment transport calculations and the bed load formula of Van Rijn, we had deposition higher than the Mount Everest and erosion deeper than the Mariana Trench.
So the bed load formula of Van Rijn - bedload_vanrijn.f - in v6p0 seems not be correct.
In the calculation of the non-dimensional skin friction
TETAP = 0.5D0 * TOB%R(I) * MU%R(I)**2 / C1 with C1 = DENS * GRAV * DM
C1 should be: C1 = 1000.0D0 * DENS * GRAV * DM
Additionally in this context I don't understand why the factor 0.5D0 and why to square the skin friction correction factor in the formula for TETAP.

In v6p1 the bed load formula of Van Rijn seems to be correct!
TETAP is calculated elsewhere, but I could'nt find the location?


Best regards,

Clemens

Hello,

I am not able to check bed load formulas, but I think that in version 6.1 TETAP is built at the beginning of subroutine bedload_formula.

With best regards,

Jean-Michel Hervouet

Hello,

thank you for the hint.
So in v6p0 the non-dimensional skin friction in the bed load formula of Van Rijn should be calculated as in v6p1 and according to the Sisyphe manual:
TETAP = TOB%R(I) * MU%R(I) / C1 with C1 = XMVE * DENS * GRAV * DM

Best regards,
Clemens