Skip to content

anuga sed crash with inlet operators #1

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
mikilterribile opened this issue Apr 11, 2024 · 3 comments
Open

anuga sed crash with inlet operators #1

mikilterribile opened this issue Apr 11, 2024 · 3 comments

Comments

@mikilterribile
Copy link

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
Copy link

r0hinik commented Sep 27, 2024

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

@mikilterribile
Copy link
Author

mikilterribile commented Sep 27, 2024

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
Copy link

r0hinik commented Oct 4, 2024

thank you so much for your quick response! can you share whatever you did on email [email protected]?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants