Hello,
For a thermal plume or dissolved chemical species in an estuarine environment we consider that it is always a matter of months, to get the data, set up the mesh and the model, deal with the problem of open boundaries at the offshore boundary. For example, the accurate calibration of a model in an estuary sometimes requires that the friction coefficient depends on time, to take into account modifications of the bottom roughness due to the tide and the ripples formation.
As for Telemac-3D, it has been intensively used for such studies, e.g. in the Gironde Estuary, Loire Estuary, and many others in the UK, Germany, etc., for what I know. For thermal plumes and salinity or any dissolved matter that modifies the density, we find that only non hydrostatic Navier-Stokes equations can give the correct concentration. This is clearly shown by our test case called "lock exchange", which we also tried with other software.
Thanks to domain decomposition, the computation time can be highly reduced with a parallel computer.
With best regards,
Jean-Michel Hervouet