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

Saturation cannot be lower than residual #458

Open
rolk opened this issue Jan 17, 2014 · 1 comment
Open

Saturation cannot be lower than residual #458

rolk opened this issue Jan 17, 2014 · 1 comment
Labels

Comments

@rolk
Copy link
Member

rolk commented Jan 17, 2014

There is a problem with having saturations lower than the residual saturation; the transport solver will alter the mass balance to get into the range of the rel.perm. curve.

Clearly, this is unphysical; I can have an aquifer filled with 100% brine which has no CO2, i.e. S_co2 = 0 initially, but which has a residual saturation, i.e. an amount of any CO2 that enters the cell later will be trapped.

Consider the following test case to elicit the problem: It consists of 5 blocks of unit size in a regular grid. There is an injector and a producer, but I am injecting with a rate of 0, so there should be no saturation change.

I run with a command like this (using commit adf421a of opm-core):

bin/sim_2p_incomp 458_ressat.param > 458_ressat.out

using these input files:

  • 458_ressat.param
  • 458_ressat.data

and these are the output files obtained:

  • 458_ressat.out
  • output/saturation/000.txt
  • output/saturation/001.txt

Notice that the saturation of CO2 is increased to slightly above the residual
saturation, in all blocks. If I change the lower bound of the rel.perm. curve,
then there is a corresponding change in the saturation.

@atgeirr
Copy link
Member

atgeirr commented Oct 21, 2015

Still valid, not resolved, not working on it, will keep open as a reminder but remove from release milestone.

@atgeirr atgeirr removed this from the Release 2015.10 milestone Oct 21, 2015
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants