Bård Skaflestad
fb2192a7ed
Implement total rate constraint for wells.
...
Untested.
2011-02-16 10:41:43 +01:00
Ove Saevareid
c7950ff341
Correcting (hopefully:)) compressibility computations.
2011-02-11 11:32:52 +01:00
Bjørn Spjelkavik
492ad8425e
Now handles PVCDO.
2011-02-09 11:28:00 +01:00
Ove Saevareid
405900ac39
Bugfix: Correct number of connections also for the last well.
2011-02-07 17:52:32 +01:00
Atgeirr Flø Rasmussen
c1e686a3d6
Removed check for wells, no longer necessary since we support them now.
2011-02-03 15:44:36 +01:00
Bjørn Spjelkavik
a27015d00c
No longer converts units.
2011-02-03 15:41:20 +01:00
Bjørn Spjelkavik
4c10c38965
No longer converts units.
2011-02-03 15:40:41 +01:00
Atgeirr Flø Rasmussen
a0c54cb6de
Added well handling to cfs_tpfa_expl_mass_transport(), simplified interface.
2011-02-03 12:51:53 +01:00
Atgeirr Flø Rasmussen
c6835bb9a7
Now the relative volume discrepancy takes std::fabs(), so max_element() will be ok.
2011-02-02 10:12:12 +01:00
Atgeirr Flø Rasmussen
4f78892dc4
Added (commented out) code to dump tables.
2011-02-01 12:40:32 +01:00
Atgeirr Flø Rasmussen
6f9ef5a32c
Added output operator for easy dumping of tables.
2011-02-01 12:40:05 +01:00
Atgeirr Flø Rasmussen
1e27b84dd8
Changed order of phases and components.
2011-01-26 16:01:04 +01:00
Atgeirr Flø Rasmussen
a779a061cc
Corrected size of gpot_storage.
2011-01-26 16:00:08 +01:00
Ove Saevareid
5c358c2de4
MiscibilityLiveOil: Additional undersaturated table entries and a few fixes.
2011-01-25 11:52:50 +01:00
Atgeirr Flø Rasmussen
6c5520a9ad
Merge.
2011-01-24 14:37:32 +01:00
Atgeirr Flø Rasmussen
fba93287c9
Fixed sign of gravity input to pressure solver.
2011-01-24 14:34:38 +01:00
Bård Skaflestad
5c3b6432b2
Explicitly state that ratio->psum is a scalar per total connection.
2011-01-24 11:21:25 +01:00
Bård Skaflestad
079d2704a6
Merge from upstream.
2011-01-20 16:07:36 +01:00
Atgeirr Flø Rasmussen
3ecfac750a
Updated to match changed C interface.
2011-01-20 16:01:58 +01:00
Bård Skaflestad
e23b9e5aff
Satisfy all post-conditions of compute_densrat_update()
...
The total (accumulated) phase contributions must be available in
ratio->psum for each completion when leaving the function.
2011-01-20 10:49:11 +01:00
Bård Skaflestad
8efa1b88e2
Add compressible well completion flows.
...
Specifically, allocate storage for 'masstrans' and 'gravtrans' terms
per well completion (perforation), as well as compressible
transmissibilities (scalar per completion). Calculate compressible
quantities by treating each completion as an interface.
Introduce a new structure, 'completion_data', to collect static and
dynamic discretisation data pertaining to each completion (e.g.,
productivity indices, gravity potentials and density ratio
operators). Pass this structure, rather than individual fields, into
affected CFS_TPFA entry points.
Compile tested only.
2011-01-19 20:20:15 +01:00
Bård Skaflestad
570ca120b9
Restructure calculation of compressible terms.
...
Specifically, rename the obtuse structure 'disc_data' to the more targeted
'densrat_util' and hoist the compressible terms 'ctrans' and 'P' into
the 'cfs_tpfa_impl' structure. Moreover, rename the remaining fields
into something that makes sense in (almost) isolation. Update
compute_densrat_update() and cfs_tpfa_construct() accordingly.
This is in preparation of adding compressible well terms.
2011-01-19 12:08:12 +01:00
Ove Saevareid
419faa324f
BugReport: MiscibilityLiveOil::miscible_oil, several issues here ...
2011-01-18 14:33:55 +01:00
Bård Skaflestad
f5eaab04fe
Reorder internal functions for call-sequence consistency.
2011-01-18 13:49:35 +01:00
Bård Skaflestad
0a42f65a52
Add gravity contributions to flux.
...
Only really verified for incompressible data. There are likely to be
programming errors.
2011-01-16 21:21:06 +01:00
Bård Skaflestad
2e3a5fff64
Improve gravity handling in linear system.
...
Specifically, only assemble gravity contributions on internal faces or
external Dirichlet faces. Moreover, pay attention to direction of
gravity flux (in/out of cell) during assembly.
2011-01-14 20:30:07 +01:00
Ove Saevareid
f7dc641f4f
Bugfix
2011-01-14 17:15:04 +01:00
Atgeirr Flø Rasmussen
827bd17d5d
Added transmissibilities to transport solver. Fixed some initialization issues.
2011-01-13 16:04:19 +01:00
Bård Skaflestad
e7214faf0f
Re-install zero-level enforcement, but only for incompressible problems.
2011-01-13 14:19:07 +01:00
Bård Skaflestad
4d6a286ab4
Remove enforcement of pressure zero level for pure Neumann problems.
...
This direct manipulation of the linear system is inappropriate when
solving compressible flows where the absolute pressure level matters
greatly.
2011-01-13 13:06:45 +01:00
Atgeirr Flø Rasmussen
fefaa301e7
Removed zero-gravity check, since we now (supposedly) handle gravity.
2011-01-13 09:23:01 +01:00
Atgeirr Flø Rasmussen
c9e9b07c5f
Made phasemobc public since it is needed by the transport solver.
2011-01-13 09:21:17 +01:00
Atgeirr Flø Rasmussen
15438e67fc
Added gravity support. Not yet tested.
2011-01-12 09:02:02 +01:00
Atgeirr Flø Rasmussen
4752cecffa
Added surfaceDensities().
2011-01-11 16:11:50 +01:00
Atgeirr Flø Rasmussen
4ce8b87c0f
Ensure proper behaviour if well object is empty.
2011-01-07 14:43:53 +01:00
Atgeirr Flø Rasmussen
da647c6658
Work in progress adding wells to transport part.
2011-01-06 11:23:58 +01:00
Bård Skaflestad
30a7644e9f
Update for sub-repo state change.
2011-01-06 10:12:02 +01:00
Atgeirr Flø Rasmussen
d478df23cf
Added missing well argument in function call.
2011-01-06 09:28:51 +01:00
Atgeirr Flø Rasmussen
2e40522940
Properly initialize with wells. New init() overload with wells.
2011-01-05 13:36:54 +01:00
Atgeirr Flø Rasmussen
16b2d16452
Fixed usage of wctrl object.
2011-01-05 11:15:51 +01:00
Atgeirr Flø Rasmussen
378ab369ed
Implemented initWells(). Added whitespace.
2011-01-05 09:48:06 +01:00
Atgeirr Flø Rasmussen
74f4f1d26a
(Last commmit message was wrong.) Working on C++ interface for wells.
2011-01-03 15:04:27 +01:00
Atgeirr Flø Rasmussen
ab31293578
Modified to account for changing interface in opmpressure.
2011-01-03 14:45:36 +01:00
Bård Skaflestad
96bcb03b91
Correct blunders concerning presence of prescribed pressure values.
...
Specifically, split detection of prescribed reservoir pressure values
from prescribed well (i.e., BHP) pressure values. Previously, we
would not even assemble any well contributions if there were any
prescribed reservoir pressure values. Moreover, the return value from
assemble_cell_contrib() was exactly opposite of its intended value...
2010-12-08 14:50:19 +01:00
Atgeirr Flø Rasmussen
1f6f70c39d
Larger number of samples, power-of-two intervals.
2010-12-08 11:15:11 +01:00
Bård Skaflestad
b990022c93
Fix build.
...
Pointy hat to self.
2010-12-07 19:35:31 +01:00
Bård Skaflestad
c52a248ea4
Merge from upstream.
2010-12-07 18:29:01 +01:00
Bård Skaflestad
68f80cecd7
Update calls to account for addition of well support in CFS_TPFA.
2010-12-07 17:55:54 +01:00
Bård Skaflestad
506546f2d2
Add retrieval of well BHP and perforation fluxes.
...
Callers may pass NULLs if any given model does not contain wells.
2010-12-07 17:30:13 +01:00
Bård Skaflestad
7830e54fd1
Add tentative Peaceman well support.
...
Callers may pass NULLs in absence of wells in any given model. This
implementation assembles an equation for each well, irrespective of
well control type (BHP or RATE), and assumes that productivity indices
and perforation pressure drops account for multiphase effects.
2010-12-07 17:12:09 +01:00