search

mperignon / anugaSed

Add-on package to ANUGA with sediment transport and vegetation-drag operators
MIT License
6 stars 2 forks source link

anuga sed crash with inlet operators #1

Open mikilterribile opened 7 months ago

mikilterribile commented 7 months ago

Hy everyone, I'm trying to use the anuga-sed solid transport operator in order to simulate a flooding event in a river. But it crashes, in each scenario it gives me the following error with the inlet operator: "assert current_volume >= 0.0 , 'Volume of watrer in inlet negative!' AssertionError: Volume of watrer in inlet negative!" Thanks in advance for your attention, Best regards, Michele Zucchelli.

r0hinik commented 1 month ago

hi mikilterribile, I am getting same error! have you resolved this issue ? If yes it would be great if you help me!

mikilterribile commented 1 month ago

Hi,

I haven't solved the problem yet. Instead, I found another developed operator for computing solid transport (I can't remember its name, but there was a past discussion on GitHub with Stephen Roberts where we talked about it). However, I was still confused about some aspects of that operator, particularly the manual balancing of erosion and deposition within the domain. This is especially problematic also because near the domain edges Anuga's divergence function doesn't work correctly, leading to significant false erosion.

So, starting from that operator, I developed my own custom operator that uses the Meyer-Peter & Müller formula to compute solid transport and addresses the issues mentioned earlier. This operator seems to work well, but it hasn't been officially tested against other models or field data yet (though it will be in the future). If you're interested, I'd be happy to share it with you.

Best regards, Michele Zucchelli

r0hinik commented 1 month ago

thank you so much for your quick response! can you share whatever you did on email rohinik@cdac.in?