Skip to content

WIP add WCD, GCDI, GCDM, WIPG, WIPL summary for F and R#4477

Merged
bska merged 1 commit into
OPM:masterfrom
totto82:summary_co2store
Feb 24, 2023
Merged

WIP add WCD, GCDI, GCDM, WIPG, WIPL summary for F and R#4477
bska merged 1 commit into
OPM:masterfrom
totto82:summary_co2store

Conversation

@totto82

@totto82 totto82 commented Feb 22, 2023

Copy link
Copy Markdown
Member

These are summary output that is used for CO2STORE simulations.

Depends on OPM/opm-common#3408

Comment thread ebos/ecloutputblackoilmodule.hh Outdated
this->fip_[Inplace::Phase::GAS][globalDofIdx] += gasInPlaceWater;
}

Scalar criticalGas = 0.1; //TODO get this from the tables

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@bska What is the best way to get the critical gas saturation here? Make it a member in this class and initialize from eclState via tables and satprop code?

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

What is the best way to get the critical gas saturation here

Just to be sure: By "critical gas saturation" you mean the saturation at which the gas phase becomes mobile—i.e., exceeds the critical relative permeability tolerance ("TOLCRIT")? If so, you typically need to care about end-point scaling and/or hysteresis which varies by cell and/or system state. In that case, I think the most general approach is to go through

simulator_.problem().materialLawManager()

and effectively do the same that we do when computing the relative permeability. Storing this as a data member will just duplicate what we already have.

@totto82 totto82 requested a review from bska February 22, 2023 11:28
@totto82 totto82 marked this pull request as ready for review February 23, 2023 08:19

@bska bska left a comment

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This looks good to me. With the caveat that I don't know how/if this works for the imbibition case I'm prepared to merge this. On a tangential note, I didn't notice that the immobile CO2 in the gas phase was tagged as CO2InGasPhaseInMob using the name InMob. At some point we might consider renaming that to CO2InGasPhaseImmob or similar instead.

Comment on lines +1112 to +1113
const auto& scaledDrainageInfo
= simulator_.problem().materialLawManager()->oilWaterScaledEpsInfoDrainage(globalDofIdx);

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Are we intentionally deferring imbibition/hysteresis?

@bska

bska commented Feb 24, 2023

Copy link
Copy Markdown
Member

jenkins build this serial please

@bska

bska commented Feb 24, 2023

Copy link
Copy Markdown
Member

I need to merge this now to restore the build following OPM/opm-common#3408 being merged. I'm still a little concerned about what to do with the imbibition curve, but we'll work on that as a refinement. For now we'll use the current version.

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.

2 participants