Help with instabilities in Celeris WebGPU
Hello, I am trying to run a simulation using Celeris WebGPU but I am observing large instabilities when refining the mesh. I was wondering if anybody had seen something similar and if maybe you were able to overcome this problem?
The simulation is around the San Antonio port and Maipo river in Chile. I am forcing a 14s period, 31.4cm amplitude wave from the west boundary of this bathymetry as a prototype. This is how the topo/bathy looks like:
You can see that in the north west portion of the domain there is a deep water region (dark-colored) that corresponds to the geological trench.
When I run it with resolutions of dx=10m or higher everything goes smoothly. But if I change to dx=4 for example, even with Courant numbers as small as 0.01 I get the result below:
I am guessing it may have to do with the narrow portion of the trench, where the bathymetry is deeper and with greater slopes, so I tried smoothing it with a gaussian filter but got a similar result.
Or maybe I am missing something. Could you help me inspect this error please?
Here are the files for the coarse mesh dx = 10m case:
https://1drv.ms/f/c/629ccbf9e520fe1c/EpdvRpL_3w9FuyHFNo9JF34B7P0l2fRiBij3H-knNLKH0A?e=DNQrXw
And these are for the fine mesh dx = 4m:
https://1drv.ms/f/c/629ccbf9e520fe1c/EsSll3mo5qxLvub0och0m9wBS3Lm5j5pp96F1VSwnbq_cg?e=DBhfIa
Thanks in advance, and please do not hesitate to ask for more details.
Thanks for your explanation, I see what you mean and it is very interesting to know that discrete BT models have this challenge.
I will keep refining with h / dx <~50 then, thank you very much.