Hi,
I hope all is well with you.
When I use an arbitrary DEM with steep slopes (in some areas exceeding 70 degrees), I notice a kind of leakage in the model. Even after increasing the number of boundary particles to 12, the leakage persists. I suspect this issue might be related to the way the boundary particles are defined, as they are not perpendicular to the DEM.
Is there a way to resolve this issue?
Thank you for your help!
Best regards,
Alireza
Hi @AlirezaZarei,
indeed, the DEM geometry is not appropriate for steep slopes. The approach used to compute the distance between particle and DEM is built on the assumption that the projection of the particle position to the DEM surface is approximately vertical. I would say that anything steeper than 45 degs is likely to encounter issues (of course, the steeper the slope, the worse the issues).
Presently in GPUSPH we don’t have logic to fill (or unfill) «arbitrary» geometries, which is something that affects also STL meshes with multi-layer boundary models, so I’m afraid there isn’t a simple solution to your problem, unless you can isolate the steeper section of your DEM to realize them using a different geometric entity (e.g. rectangles or something like that).
If you have an example DEM that is giving you the issue, can you send me by email a test case and the DEM so I can give it a quick look?