Hi
I am modelling the hydrodynamics in the Gulf of Guayaquil, a large estuary. The seaward boundary is defined by a TPXO model. The upstream boundary by two boundaries with imposed discharges. For Bathymetry I am using a combination of the GEBCO database and data from local institutions.
When running the simulations, the model shows clear signs of numerical unstability close to the seaward boundary (see video attached).
I looked around on the forum and so far learned that the constant viscosity turbulence model with a not-too-small coefficient. While this already improved the stability compared to using a k-epsilon model, it is still not very stable. Time step is 3 sec and the stabilities are mostly appearing at the area where the spatial resolutions is the largest so I don't expect any Courant number violations.
Do you have any tips on further improving this?
I added the input cas file and resulting output video. Please let me know if I should any other files as well.
video:
drive.google.com/file/d/1ruGsnIM19Ke4Ghj...uts/view?usp=sharing
Kind regards,
Ignace Pelckmans