open TELEMAC-MASCARET The mathematically superior suite of solvers

I am trying to run a simulation with telemac3d with salinity as a tracer. This simulation is in the Gironde estuary on the southwest coast of France. Each simulation that I have ran either resulted in the salinity being pushed offshore with little to no vertical stratification or the model crashes. I have tried different combinations of both the k-epsilon and mixing length vertical turbulence models with the NON-HYDROSTATIC parameter both ON and OFF and the DENSITY LAW either 2 or 0. However, each time I set DENSITY LAW=2, the model crashed. I have also tried various numbers of horizontal planes (from 4 to 10) and timesteps (10 sec to 5 sec). Has anyone else come across this issue? I have attached one of my steering files for reference. This is a case using the MIXING LENGTH vertical turbulence model that ran, but the salinity iwas all pushed offshore. Thank you in advance for your help.

Best,
Lauren

Hello,

Yes, only the non hydrostatic version is relevant in that case.

You can try : SCHEME FOR ADVECTION OF VELOCITIES = 1
which will be less diffusive, but I see no real problem in your steering file, though the HORIZONTAL TURBULENCE MODEL would be better at 1 (constant value).

At least, to be sure that nothing is wrong and that your offshore boundary is boundary 1, you should see the salinity entering the domain when the tide is coming up ?

Regards,

Jean-Michel Hervouet