-
Notifications
You must be signed in to change notification settings - Fork 91
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
Fix for black oil #2937
Fix for black oil #2937
Conversation
phaseCompFraction.value[ipOil][icOil] = zo; | ||
phaseCompFraction.value[ipOil][icGas] = zg; | ||
phaseCompFraction.value[ipOil][icOil] = zo / ( zo + zg ); | ||
phaseCompFraction.value[ipOil][icGas] = zg / ( zo + zg ); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I would have written this as for example zo / (1 - zw)
which should have no impact here but would have an impact on the derivatives. Have you looked at the Jacobian in this alternative formulation. I find these composition derivatives a bit confusing when honouring that they sum to 1.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
just checked - Jacobian is the same, number of linear iterations is the same, I will update with your version
Codecov ReportAll modified and coverable lines are covered by tests ✅
Additional details and impacted files@@ Coverage Diff @@
## develop #2937 +/- ##
========================================
Coverage 51.33% 51.33%
========================================
Files 968 968
Lines 86817 86819 +2
========================================
+ Hits 44570 44572 +2
Misses 42247 42247 ☔ View full report in Codecov by Sentry. |
* fix for black oil --------- Co-authored-by: Pavel Tomin <“paveltomin@users.noreply.github.com”> Co-authored-by: Matteo Cusini <49037133+CusiniM@users.noreply.github.com>
Fix small bug in black oil computations
One test case requires rebaseline, rebase PR here https://github.com/GEOS-DEV/integratedTests/pull/84