Skip to content
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

Only gravity driven Outlet flow #549

Merged
merged 27 commits into from
Sep 13, 2023
Merged

Only gravity driven Outlet flow #549

merged 27 commits into from
Sep 13, 2023

Conversation

SouthEndMusic
Copy link
Collaborator

@SouthEndMusic SouthEndMusic commented Aug 25, 2023

Fixes #548.

Currently the outlet switches on and off in a discrete way due to a switch in sign of the level difference. This can result in the outlet rapidly turning on and off:

flipperende_outlet

To solve this problem, a reduction factor should be applied to the outlet flow, based on the level difference. This reduction factor will have to be taken into account in the implementation of the analytical Jacobian (with an unsatisfying effort to improvement ratio). This I think will be a recurring problem until we are able to switch to an AD Jacobian.

@SouthEndMusic SouthEndMusic marked this pull request as draft August 25, 2023 12:21
@SouthEndMusic SouthEndMusic changed the title Implement hard stop of outlet flow Only gravity driven Outlet flow Aug 25, 2023
@SouthEndMusic
Copy link
Collaborator Author

SouthEndMusic commented Sep 1, 2023

I made this plot of the smooth reduction factor:

reduction_factor_plot

how can I add it to equations.qmd? It has been added now.

@SouthEndMusic SouthEndMusic marked this pull request as ready for review September 4, 2023 18:00
@SouthEndMusic
Copy link
Collaborator Author

I found something peculiar when working on the test model for the Outlet:

outlet_model

outlet_constraints

In the current state it behaves as expected, that is: there is no outlet flow until the level boundary reaches the minimum crest level, and the basin level converges to the level boundary level. However, when the basin starts completely empty instead of with a tiny amount of water, it doesn't fill:

outlet_constraints_noflow

I haven't been able to figure out why this happens.

@SouthEndMusic SouthEndMusic requested a review from visr September 4, 2023 18:14
I prefer to avoid yellow since it is the color of selected nodes in QGIS.
@visr visr mentioned this pull request Sep 13, 2023
@visr visr merged commit 82b7279 into main Sep 13, 2023
@visr visr deleted the physical_outlet_flow branch September 13, 2023 14:37
visr added a commit that referenced this pull request Sep 14, 2023
Fixes #548.

Currently the outlet switches on and off in a discrete way due to a
switch in sign of the level difference. This can result in the outlet
rapidly turning on and off:

![flipperende_outlet](https://github.com/Deltares/Ribasim/assets/74617371/d834ffa9-ac48-4fb1-99c5-732369f43c35)

To solve this problem, a reduction factor should be applied to the
outlet flow, based on the level difference. This reduction factor will
have to be taken into account in the implementation of the analytical
Jacobian (with an unsatisfying effort to improvement ratio). This I
think will be a recurring problem until we are able to switch to an AD
Jacobian.

---------

Co-authored-by: Martijn Visser <[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

Successfully merging this pull request may close these issues.

Let outlet only flow from higher to lower level
2 participants