open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi JMH,

I am still not getting a good mass conservation. Could you kindly take a quick look at my case file attached here, kindly let me know if some thing seems wrong.

Best regards,
Sumit

Hello,

Try this instead of your last section of keywords :

OPTION FOR THE TREATMENT OF TIDAL FLATS : 1
TREATMENT OF NEGATIVE DEPTH : 2 / not mandatory but ensures exact mass conservation and positive depths

TYPE OF ADVECTION = 1;5
SUPG OPTION = 0;0 / replaced by next line
COMPATIBILITY OF FREE SURFACE GRADIENT : 0.8
TREATMENT OF THE LINEAR SYSTEM = 2
SOLVER = 1
IMPLICITATION FOR DEPTH = 1. / faster and more stable
IMPLICITATION FOR VELOCITY = 1.
CONTINUITY CORRECTION = YES
MASS-LUMPING ON H = 1.
H CLIPPING = NO
PARALLEL PROCESSORS = 48


You can also have a velocity diffusivity larger than 1.E-3, in a real river 0.1 to 1 is OK.

The main problem for mass conservation was H clipping, when you clip a negative depth it adds water in your domain.

Best regards,

Jean-Michel Hervouet

Thanks a lot JMH, it fixed the mass conservation issue.