open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear All,

Hi.
I am trying to simulate the scouring phenomena around a underwater structure.
I have attached all the input files.
I am keep get "exceeding maximum iterations" errors and I don't know how to solve the problem. Could someone please take a look at my files and let me know any problems in the files?
Thank you.

Hello,

Even without coupling with sisyphe, your case doesn't run.
I advise you to first set-up your telemac2d model alone.
You're trying to prescribe velocities of 20cm/s for 15cm of water depth over a domain of 1m long (element size of about 1cm). Is it really what you want? If yes, I think you should decrease your time step for the solver to converge.

You're using 2 sisyphe subroutines as user fortran, but you declare it as a "FRICTION DATA FILE". I think it's rather FORTRAN FILE keyword, otherwise you must provide a valid friction data file.

"PRECONDITIONING FOR DIFFUSION OF TRACERS" is useless as you don't use any tracer.

Your sisyphe steering file contains both BEDLOAD and SUSPENSION processes. For bedload computation, you're using Soulsby/Van Rijn formula which already accounts for suspension effects, so you're considering it twice by turning on SUSPENSION too. I advise you to turn off the suspension, or to choice another bedload formula that does not account for suspension effects.

Regards,
Laurent

Hi Laurent,

There is no problem with using Soulsby van Rijn and SUSPENSION = yes in the recent versions of TELEMAC-MASCARET. SISYPHE checks if suspension is calculated. If that is the case, the suspended load part of the total load formulations is ignored. See the last bit of code in bedload_formula.f

Michiel

Yes, thank you Michiel for pointing it out to me!

Hi,

as Laurent points out, the time step for the computations needs to be reduced. With the velocities and grid size, it will have to be around 0.01s or so.

In addition to his comments:
As the unit of the grain diameter is meter, you are specifying a grain size of 1mm. Is that correct?

Also, with the depths and velocities you are using and no wave forcing, Soulsby van Rijn will predict there is no sediment transport at all (very likely not even around your structure). You might be better off using van Rijn or another transport predictor.