OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-26 03:00:17 -06:00
Code Issues Packages Projects Releases Wiki Activity
3,885 Commits 24 Branches 101 Tags 97 MiB
60fd4add69
Commit Graph

533 Commits

Author SHA1 Message Date
Liu Ming
c2cdc7ec17 add dead poro volume functionality. 2013-12-12 22:46:29 +08:00
Liu Ming
2eaf24decf fix errors for debugging. 2013-12-12 21:58:25 +08:00
Liu Ming
25160d019a add adsorption term for polymer equation. 
rewrite some function for simplify.
 2013-12-12 21:24:47 +08:00
Liu Ming
0558184439 modify the function for the polymer term. 2013-12-12 20:29:40 +08:00
Liu Ming
720e7ebbef fix derivative of effective relperm respect to concentration. 2013-12-12 17:44:00 +08:00
Liu Ming
a7db4f3cf0 add adsorption and effective relperm functions. 2013-12-11 22:52:11 +08:00
Liu Ming
11ea882ed8 remove some output commits. 
get values from command line for exapmes.
 2013-12-11 18:27:05 +08:00
Liu Ming
a358da7afa fix the polymer source bug, 
warnning: water initial saturation should not bt zero when
running this simulator.
 2013-12-10 20:36:30 +08:00
Liu Ming
6fc24236df fix FullyImplicitTwophasePolymersolver constructor problem. 
add some input commits for debugging.
 2013-12-09 22:57:44 +08:00
Liu Ming
fb12565ddf make PolymerState return numPhases. 2013-12-09 20:34:23 +08:00
Liu Ming
e226b5d95e Simulator for fullyimplicit two phase flow with polymer 2013-12-09 20:33:05 +08:00
Liu Ming
080116c66b Add fully implicit solver for incomp two phase with polymer 
and the polymer properties based on AD.
 2013-12-06 23:35:13 +08:00
Liu Ming
eded8f735c chang include header from opm/autodiff to opm/polymer/fullyimplicit 2013-12-06 22:25:22 +08:00
Liu Ming
f6acccd396 incompressible two phase fully implicit simulator. 2013-12-06 22:13:29 +08:00
Andreas Lauser
44e0848127 make headers autonomously includable and fix a few harmless clang warnings 2013-09-10 19:41:49 +02:00
Andreas Lauser
4859e77c79 explicitly include <iostream> in files which use std::cout or std::cerr 2013-09-05 12:39:11 +02:00
Andreas Lauser
05af2ab6ff Replace MESSAGE by OPM_MESSAGE 2013-09-03 15:43:58 +02:00
Andreas Lauser
7e8d941a6f Replace THROW by OPM_THROW 2013-09-03 15:43:58 +02:00
Andreas Lauser
6a6d0c9b28 Replace the ASSERT and ASSERT2 macros by assert 2013-09-03 15:43:51 +02:00
Atgeirr Flø Rasmussen
d065a908d3 Fix issue 27, bug in call to viscosity(). 2013-06-10 10:31:18 +02:00
Atgeirr Flø Rasmussen
7d5d314fc0 Update to match opm-core changes. 2013-03-18 13:10:32 +01:00
Atgeirr Flø Rasmussen
c16cb6987b Renames classes in line with opm-core changes. 2013-03-15 18:33:41 +01:00
Atgeirr Flø Rasmussen
04abba4d5e Adapt to changed opm-core header name. 2013-03-15 17:55:24 +01:00
Atgeirr Flø Rasmussen
864cce5f5e Adapt to moved opm-core headers. 2013-03-14 13:59:31 +01:00
Atgeirr Flø Rasmussen
23cef788fc Adapt to changed include paths for some opm-core headers. 2013-03-12 11:08:33 +01:00
Atgeirr Flø Rasmussen
b4d5e03c61 Add eclipse output to simulator. 
Note: so far only for incompressible fluid simulator.
 2012-11-26 11:01:23 +01:00
Atgeirr Flø Rasmussen
b6b59918e1 Suppress output of multi-cell-block sizes and #iterations. 2012-11-26 11:00:20 +01:00
Atgeirr Flø Rasmussen
bcfcfdc930 Use minimum 5 digits for step numbers in output file names. 
Was 3 digits, created trouble for certain scripts reading
the files when more than 1000 steps were taken.
 2012-11-05 15:37:14 +01:00
Atgeirr Flø Rasmussen
2b292108ad Added well management to SimulatorCompressiblePolymer. 
Another polymer-specific overload of computeFractionalFlow()
has been added in support of this (this overload deals with
compressible properties.

Also added pressure normalization for situations with arbitrary
absolute pressure.
 2012-10-15 14:23:13 +02:00
Atgeirr Flø Rasmussen
f5861f0cc8 Added well management to SimulatorPolymer. 
A polymer-specific overload of computeFractionalFlow()
has been added in support of this.

Also added a minor feature: pressure normalization for
situations with arbitrary absolute pressure.
 2012-10-10 23:19:35 +02:00
Xavier Raynaud
9cf746b0b0 Removed TAB spaces. 2012-10-08 15:00:19 +02:00
Xavier Raynaud
20d7cf80ea Reintroduced Newton Simple (in s, c variables). 2012-10-08 13:07:22 +02:00
Xavier Raynaud
ef053a1883 Cleaned up single cell Newton solver. 2012-10-08 10:29:45 +02:00
Atgeirr Flø Rasmussen
2f14ed70c4 Modified compressible transport src term. 
Using surface volume for injection rather than reservoir volume.
Chase change in class TransportModelCompressiblePolymer and
function computeInjectedProduced().

Also changed interface to computeInjectedProduced() to take state
rather than individual state variables.
 2012-10-04 22:32:55 +02:00
Atgeirr Flø Rasmussen
f4abfd6cbb Fix mass balance report for multi-epoch case. 
Compute initial polymer amount instead of setting it to zero.
Error obvious with multiple epochs.
 2012-10-04 20:02:34 +02:00
Atgeirr Flø Rasmussen
e13e77a7bb Simulators handle WPOLYMER keyword. 
They do this by using the class PolymerInflowFromDeck
instead of PolymerInflowBasic if there is a WPOLYMER
keyword somewhere in the deck epochs. If there is no
WPOLYMER, the parameters 'poly_start_days' etc. will
be used to construct an instance of PolymerInflowBasic
instead.
 2012-10-04 16:00:33 +02:00
Atgeirr Flø Rasmussen
f2e7c0bece Made getInflowValues() a const method. 2012-10-04 15:59:15 +02:00
Atgeirr Flø Rasmussen
0428c4abb4 Moved PolymerInflow* classes to separate file. 
Also, implemented PolymerInflowFromDeck (but not yet used).
 2012-10-04 15:15:32 +02:00
Atgeirr Flø Rasmussen
6033e5ca69 Make polymer inflow specification by cell. 
Before, it was a single number for the whole domain. It could vary
with time, though.

Simulator behaviour should be identical before and after this commit.

This touches the classes TransportModel*Polymer, Simulator*Polymer,
and the computeInjectedProduced() functions. The PolymerInflow class
and its usage has been replaced with PolymerInflowInterface, and
a subclass PolymerInflowBasic has been created which provides the
old behaviour (using parameters poly_start_days, poly_end_days and
poly_amount).
 2012-10-04 13:38:22 +02:00
Atgeirr Flø Rasmussen
b32c212e56 Fix mass balance report. 
Now reports surface volumes for water and oil, instead of useless
(for mass balance) reservoir volumes. Also makes the case with
multiple transport steps more accurate.
 2012-10-04 11:25:42 +02:00
Atgeirr Flø Rasmussen
c6b76e715d Fix for case with incompressible rock. 
Made rock comp argument optional in computePolymerAdsorbed().
Check inside func for active rock comp.
 2012-10-03 16:17:54 +02:00
Atgeirr Flø Rasmussen
dbd4c9e7b5 Reinstate 10% ad-hoc c-interval overestimation. 
Similar to in TransportModelPolymer.
 2012-10-03 16:16:08 +02:00
Atgeirr Flø Rasmussen
25e528df0d Make output work the same as in other sims. 
I.e. putting each field in its own directory, adding the
output_vtk parameter etc.
 2012-10-03 16:02:27 +02:00
Atgeirr Flø Rasmussen
dc3b0dd6f4 Improve function comment (similar to in SimulatorPolymer). 2012-10-03 16:01:54 +02:00
Atgeirr Flø Rasmussen
753a0c6949 Reinstate 10% ad-hoc c-interval overestimation. 2012-10-03 16:00:51 +02:00
Atgeirr Flø Rasmussen
2990b0a3e8 Improve function comment. 2012-10-03 16:00:13 +02:00
Atgeirr Flø Rasmussen
349c81fb1e Stop writing ascii output twice. 2012-10-03 15:59:18 +02:00
Atgeirr Flø Rasmussen
678611166d Improve reporting of volume balance. 
Should now give more accurate results when taking multiple substeps.
 2012-10-03 14:25:20 +02:00
Atgeirr Flø Rasmussen
076164d3f3 Add "output_vtk" parameter, split output function in two. 
Vtk and Matlab output now happens in two different functions.
 2012-09-17 07:54:50 +02:00
Xavier Raynaud
4db4d715de Last saturation is stored in a mutable variable. 2012-09-10 11:03:59 +02:00
Halvor M. Nilsen
594eca3633 Added lines which is need to compile with PROFILING. Did not check if code brakes. 2012-09-08 11:11:23 +02:00
Xavier Raynaud
55486d3c03 Fixed source term (measured at reservoir condition). 2012-09-05 13:38:30 +02:00
Xavier Raynaud
938fc0f0d4 Fixed comment line. 2012-09-05 12:34:17 +02:00
Xavier Raynaud
7721defbd3 Fixed/Added documentation comments. 2012-09-05 11:45:54 +02:00
Xavier Raynaud
5f27612b05 Merged from upstream. 2012-09-05 11:11:06 +02:00
Xavier Raynaud
7fb671ffbf Changed computation of polymer production. Added functionalities for compressible case. 2012-09-05 11:09:45 +02:00
Bård Skaflestad
36e017d48d Honour initialisation order requirement. 
Introduced in change-set 85cde6f721a4.
 2012-09-05 10:27:37 +02:00
Xavier Raynaud
8e089dd80d Use porevolume instead of porosity, as argument for the transport solver. 2012-09-05 09:11:05 +02:00
Xavier Raynaud
fa84128ddc Fixed bug. 2012-09-05 09:07:47 +02:00
Xavier Raynaud
612a1ffc63 Fixed bug. 2012-09-03 16:46:02 +02:00
Xavier Raynaud
27dea09311 Added initializations in constructor. 2012-09-03 16:45:50 +02:00
Xavier Raynaud
8752f14015 Added executable example for compressible with polymer. 2012-09-03 13:52:59 +02:00
Xavier Raynaud
8bf900f1bf Update surface volume in transport solver. 2012-09-03 12:52:25 +02:00
Xavier Raynaud
708c56bb94 Simplified computation of adsorbed mass of polymer. 2012-09-03 12:49:04 +02:00
Xavier Raynaud
a335e1b149 Added Simulator for compressible polymer. Not tested. 2012-08-31 10:38:50 +02:00
Xavier Raynaud
8b40a4a68f Fixed comments. 2012-08-30 15:58:14 +02:00
Xavier Raynaud
99fe7bf46e Fixed preprocessor flag names in some header files. 2012-08-30 15:31:28 +02:00
Xavier Raynaud
92def1d514 Cleaned up ownership (private vs public) in class definition. 2012-08-30 09:09:23 +02:00
Xavier Raynaud
eb52ee58c1 Added gravity transport solver functionality. Not tested. 2012-08-29 13:49:02 +02:00
Xavier Raynaud
928513b561 Added compressible polymer transport solver. Not yet done. 2012-08-28 16:42:26 +02:00
Xavier Raynaud
9a271a6038 Added (fluid) compressible pressure solver for polymer. 2012-08-27 13:16:10 +02:00
Xavier Raynaud
72658ac91c Simplified code: Only one occurence of residual evaluation. 2012-08-23 17:18:19 +02:00
Xavier Raynaud
51e0512989 Removed non-used members in class ResidualEquation. 2012-08-23 16:51:51 +02:00
Xavier Raynaud
52a19192cb Added PROFILING preprocessor flag. 2012-08-23 14:58:26 +02:00
Halvor M. Nilsen
497822b893 Added posibility to use eclipse fluid for 1D case. Corrected bug in guess for newton solver 2012-07-03 11:05:51 +02:00
Halvor M. Nilsen
0d97945aba Removed safetyfactor with membervariable which for now is put to 1.0. Made jacobian calculations more transparent and fixed pointer error. 2012-06-26 08:37:08 +02:00
Halvor M. Nilsen
aba26a0583 Tried to make derivative genneral 2012-06-23 15:04:59 +02:00
Halvor M. Nilsen
b13416be86 Made new single cell solvers. 2012-06-15 14:40:07 +02:00
Atgeirr Flø Rasmussen
92bc55385d Update porosities before calling transport solver. 2012-06-14 21:53:56 +02:00
Xavier Raynaud
168c4c9e3d Optimized newton solver. 2012-06-14 16:52:03 +02:00
Xavier Raynaud
95c3756443 merge. 2012-06-14 14:35:06 +02:00
Atgeirr Flø Rasmussen
cd9be1030a Adapt to moved SimulatorTimer class, and use common SimulatorReport class. 2012-06-14 14:26:29 +02:00
Xavier Raynaud
6e7980ba47 fixed bug in relperm computation. 2012-06-14 14:01:53 +02:00
Xavier Raynaud
5f4b68fe33 merge profiling branch. 2012-06-14 13:39:12 +02:00
Xavier Raynaud
f3eb14e528 Added sc formulation for gradient method. 2012-06-14 13:37:08 +02:00
Atgeirr Flø Rasmussen
eb845fd45a Added reporting of time usage to interface, similar to SimulatorTwophase. 2012-06-14 13:24:38 +02:00
Atgeirr Flø Rasmussen
d9ea18dc9d Bugfix: fix dkr/ds matrix in IncompPropertiesDefaultPolymer::relperm(). 
Note: this class seems unused now. Can it be removed?
 2012-06-14 13:16:24 +02:00
Xavier Raynaud
9df87db47f fixed bug in computation of derivatives for polymer. 2012-06-14 09:35:45 +02:00
Xavier Raynaud
6620a1b5cd In profiling branch: merged last default branch. 2012-06-14 08:38:01 +02:00
Atgeirr Flø Rasmussen
502215832d Fix erroneous includes again. 2012-06-13 21:59:11 +02:00
Atgeirr Flø Rasmussen
7a40153e51 Added IncompTpfaPolymer class. 2012-06-13 19:56:44 +02:00
Xavier Raynaud
60f2fd2795 updated profiling branch. 2012-06-13 17:55:33 +02:00
Xavier Raynaud
a7e91a0bbc cleaned up TransportModelPolymer (default branch) 2012-06-13 17:52:07 +02:00
Xavier Raynaud
56ef0105af fixed wrong include. 2012-06-13 17:51:47 +02:00
Xavier Raynaud
f04f80e696 cleaned up TransportModelPolymer (profiling branch) 2012-06-13 17:36:36 +02:00
Xavier Raynaud
ab1591df40 fixed computation of derivatives in PolymerProperties.cpp (default branch) 2012-06-13 17:07:08 +02:00
Xavier Raynaud
4732ca3bbc Fixed computation of derivatives in PolymerProperties (profiling branch). 2012-06-13 17:01:23 +02:00
Xavier Raynaud
564ab4bc44 Removed some standard output messages. 2012-06-13 16:13:26 +02:00
Xavier Raynaud
20b87da8d5 Removed some standard output messages (profiling branch) 2012-06-13 16:12:32 +02:00
Xavier Raynaud
acaaada2d0 merge. 2012-06-13 15:56:44 +02:00
Xavier Raynaud
4c7c439dd7 merge. 2012-06-13 15:56:09 +02:00
Atgeirr Flø Rasmussen
005cda508c Improved docs. 2012-06-13 15:55:23 +02:00
Atgeirr Flø Rasmussen
876335908d Added class SimulatorPolymer. 2012-06-13 15:46:15 +02:00
Atgeirr Flø Rasmussen
42e862e7ec Added TransportModelPolymer::setPreferredMethod(). 2012-06-13 15:44:46 +02:00
Atgeirr Flø Rasmussen
b6001c8dac Added helper class PolymerInflow (moved from polymer_reorder.cpp). 2012-06-13 15:44:12 +02:00
Xavier Raynaud
9637906a35 Make compatible with incomptpfa (default branch). 2012-06-13 15:33:31 +02:00
Xavier Raynaud
a5de5341c5 Make compatible with incomptpfa interface. 2012-06-13 15:31:49 +02:00
Xavier Raynaud
6033527a82 Added pure Newton single cell method. Correction of Jacobi. sc formulation. 2012-06-13 15:11:30 +02:00
Atgeirr Flø Rasmussen
1e4db53e7f Bugfix: use water saturation vector where it is expected, not full saturation. 2012-06-13 14:35:18 +02:00
Atgeirr Flø Rasmussen
ffd33be55d Changed interface of TransportModelPolymer to expect full saturation, not just water. 2012-06-13 14:18:32 +02:00
Atgeirr Flø Rasmussen
c4b81d8da5 Removed two unnecessary arguments from TransportModelPolymer constructor. 2012-06-13 12:46:57 +02:00
Xavier Raynaud
022b1ccb11 Small changes in Newton/Gradient method. 2012-06-11 15:59:33 +02:00
Xavier Raynaud
02c5935865 Added profiling branch. Compute number of residual evaluations. 2012-06-11 14:44:21 +02:00
Xavier Raynaud
9eb9ba372d merge 2012-06-11 14:40:55 +02:00
Xavier Raynaud
6a9dc03085 Change copy assignment to reference. 2012-06-07 13:45:09 +02:00
Xavier Raynaud
ad1734265f Added documentation. 2012-06-06 17:03:19 +02:00
Atgeirr Flø Rasmussen
3660ec8e8c Follow headers that were moved in opm-core. 2012-06-05 15:45:41 +02:00
Atgeirr Flø Rasmussen
c3dac20d65 Use new interface for regula falsi, switch to non-throwing error policy. 2012-05-24 10:27:27 +02:00
Xavier Raynaud
91fe1c5d63 Removed useless commented regions. 2012-05-23 16:24:08 +02:00
Xavier Raynaud
4b805c191e Added bound checks to Newton column solver. 2012-05-23 16:22:48 +02:00
Xavier Raynaud
263a41b84b treat case when searching direction is small in Newton reorder solver. 2012-05-22 17:43:18 +02:00
Xavier Raynaud
afca287757 In Newton reordering solver, set smin to zero. 2012-05-22 17:22:31 +02:00
Xavier Raynaud
6a12fe0208 Added default concentration if guess_old_solution is false. 2012-05-22 17:21:12 +02:00
Xavier Raynaud
a9a93d795e Include concentration in error estimation for gauss seidel. 2012-05-22 10:45:43 +02:00
Xavier Raynaud
8e0936fadd Fixed bug in gauss seidel gravity solver. 2012-05-22 09:51:32 +02:00
Xavier Raynaud
48adf98b5a Corrected some mistakes in Gauss Seidel gravitation transport solver. 2012-05-21 16:41:29 +02:00
Xavier Raynaud
12bb246c06 corrected bug in SinglePointUpwindTwoPhasePolymer. 2012-05-21 16:35:39 +02:00
Atgeirr Flø Rasmussen
c115ba2c28 Adapt to renaming and other small modifications of state and initialization. 2012-05-16 15:59:37 +02:00
Xavier Raynaud
360ec5f86d merge. 2012-05-16 15:09:36 +02:00
Xavier Raynaud
76e82a718e fixed mistake in computation of one of the residuals. 2012-05-16 15:09:07 +02:00
Xavier Raynaud
60599b9dbf Implemented Gauss-Seidel solver for polymer. Not tested! 2012-05-16 15:06:14 +02:00
Xavier Raynaud
c0368fa61d Fixed problem induced by s=0 in gravitation transport solver. 2012-05-15 16:09:07 +02:00
Atgeirr Flø Rasmussen
088c61ce00 Added check for single-cell columns. 2012-05-14 11:04:00 +02:00
Atgeirr Flø Rasmussen
67a6674495 Update for simplified column data structure. 2012-05-11 09:41:01 +02:00
Xavier Raynaud
1eb75314e6 Make polymer_reorder.cpp compliant with well management. Not tested! 2012-05-10 14:42:08 +02:00
Xavier Raynaud
0a0ca3bbd3 Added compressibility to polymer reorder transport solver. 2012-05-10 09:42:35 +02:00
Xavier Raynaud
9010f2da42 Fixed bug in Newton solver for polymer. 2012-04-26 14:16:40 +02:00
Xavier Raynaud
6748a6e5e9 some code optimization (instead of mu_w_eff, use its inverse). 2012-04-26 09:25:46 +02:00
Xavier Raynaud
af79b9555c Fixed bug in computation of effectiveTotalMobility. 2012-04-25 09:23:59 +02:00
Xavier Raynaud
3ecb7bf88c Removed compiler warnings. 2012-04-25 08:51:32 +02:00
Xavier Raynaud
ab961dda94 more efficient polymer code (replaced std::vector by double*). 2012-04-25 08:39:39 +02:00
Xavier Raynaud
6eb68c9652 Fixed bug in SinglePointUpwindTwoPhasePolymer. 2012-04-24 11:20:06 +02:00
Xavier Raynaud
e867c53722 polymerUtilities uses now PolymerProperties to compute effective mobilities. 2012-04-23 11:49:05 +02:00
Xavier Raynaud
4794348740 added computeMc function in PolymerProperties. 2012-04-23 10:27:27 +02:00
Xavier Raynaud
88941610ed Corrected a bug in computation of effective mobilities. 2012-04-23 10:07:51 +02:00
Xavier Raynaud
a0794117c9 Cleaned code for computation of residual, many changes, not tested! 2012-04-19 17:21:08 +02:00
Atgeirr Flø Rasmussen
eb502faae7 Moved AdHocProps to new Inc.Props.DefaultPolymer class. 
Also moved some functions in polymer_reorder.cpp to be more similar to spu_2p.cpp.
 2012-04-12 14:10:47 +02:00
Atgeirr Flø Rasmussen
c99e396ed1 Silence warnings related to ad-hoc polymer fluid wrapper class. 2012-04-12 09:07:15 +02:00
Atgeirr Flø Rasmussen
6e26837dc7 Silence warnings resulting from unneeded redeclaration of dps etc. 2012-04-12 09:01:09 +02:00
Atgeirr Flø Rasmussen
b197c6b86e Added rock compressibility terms to residuals (bracketing method only). 2012-04-11 16:04:04 +02:00
Atgeirr Flø Rasmussen
d71742871b Expand concentration interval in which we look for solution. 
This is an ad-hoc fix. We should investigate if (or how much)
the concentration may be above from the injected value.
 2012-04-10 14:08:32 +02:00
Atgeirr Flø Rasmussen
dce3068343 Modified some assertions. 
- In a gravity segregation setting, total upwind mobility may be zero.
 - Make assertion more readable by removing negation.
 2012-04-10 14:05:39 +02:00
Atgeirr Flø Rasmussen
8b8d10d58f Comment unused variables to silence warnings. 2012-04-10 10:25:29 +02:00
Atgeirr Flø Rasmussen
adf40b19c1 Fixed initialization bug in PolymerProperties. 2012-04-10 10:25:07 +02:00
Atgeirr Flø Rasmussen
9fadfc316e Added clarification. 2012-04-02 15:45:34 +02:00
Atgeirr Flø Rasmussen
0e178199ce Added simulator state class for polymer case. 2012-04-02 15:44:32 +02:00
Xavier Raynaud
515bcb485e In transport-reorder Newton, limit interval check for s, c (reset if needed). 2012-03-28 13:10:53 +02:00
Xavier Raynaud
18e45caf4f Added enum to choose gradient computation method. By default, analytic. 2012-03-28 10:14:24 +02:00
Atgeirr Flø Rasmussen
724645b3f8 Introduced enum for adsorption behaviour to make code more self-documenting. 2012-03-27 15:33:48 +02:00
Xavier Raynaud
252fada2b4 Adapted to change in gravitation solver interface. 2012-03-27 13:05:53 +02:00
Xavier Raynaud
b846b15f1f Corrected typo: adsorbtion -> adsorption. 2012-03-27 11:26:55 +02:00
Xavier Raynaud
acedd4b596 Added adsorbtion in residual computation. 2012-03-27 11:20:20 +02:00
Xavier Raynaud
6d9f0ea3ec Added support for adsorbtion index (desorbtion allowed or not). 2012-03-26 16:37:39 +02:00
Atgeirr Flø Rasmussen
294c8b03be Bugfix in calculation of adsorbed polymer for balance reports. 2012-03-26 15:33:38 +02:00
Xavier Raynaud
023cbb228e Fixed bug for computeMc. 2012-03-26 13:45:14 +02:00
Xavier Raynaud
4e3f486d13 Simplified computation of coefficient mc. 2012-03-26 11:31:12 +02:00
Xavier Raynaud
01aab0270b Fixed bug in lapack call for band matrices. 2012-03-26 11:29:19 +02:00
Xavier Raynaud
1836247dab Added actual resistance factor. Not tested. 2012-03-23 17:40:03 +01:00
Xavier Raynaud
4c77354d56 Changed the definition of dead pore space. 2012-03-23 16:08:41 +01:00
Xavier Raynaud
62a1c1ad98 Corrected various bugs in gravitation solver. 2012-03-20 17:43:25 +01:00
Xavier Raynaud
49b9f340ac Added functions to compute mobilities and coefficient mc (with derivatives). 2012-03-20 11:58:09 +01:00
Xavier Raynaud
542ed94b84 removed bugs. can compile. Not tested. 2012-03-19 16:33:32 +01:00
Xavier Raynaud
b6b720f338 Removed useless comment. 2012-03-19 11:30:00 +01:00
Xavier Raynaud
2daf397268 removed dmcds because the coefficient mc does not depend on s. 2012-03-19 11:28:55 +01:00
Xavier Raynaud
b513b129ce changed syntax for BandMatrixCoef. 2012-03-19 09:56:20 +01:00
Xavier Raynaud
dc0771ee3e merge 2012-03-15 18:10:10 +01:00
Xavier Raynaud
efe0408b71 Gravity column solver with polymer. 2012-03-15 18:09:29 +01:00
Xavier Raynaud
54844f60c1 Computes s and c-residual and its Jacobian at a cell interface. 2012-03-15 18:08:30 +01:00
Atgeirr Flø Rasmussen
b120aa6c90 Bugfix: corrected computation of polymer contained in reservoir. 2012-03-14 21:12:54 +01:00
Atgeirr Flø Rasmussen
1554bfff11 Still working on adsorption computations. 2012-03-14 16:13:59 +01:00
Atgeirr Flø Rasmussen
efc03f1596 Remove unneeded dps argument from computePolymerAdsorbed(). 2012-03-14 15:59:56 +01:00
Atgeirr Flø Rasmussen
2f06b301bf Change function name to computePolymerAdsorbed() [with a d, not b]. 2012-03-14 15:57:28 +01:00
Atgeirr Flø Rasmussen
ad470a7796 Bugfix: dead pore space in polymer mass calc. Added computePolymerAbsorbed(). 2012-03-14 15:54:36 +01:00
Atgeirr Flø Rasmussen
bf44f9f72e Added computePolymerMass() and polymer-aware overload of computeInjectedProduced(). 2012-03-14 14:03:07 +01:00
Xavier Raynaud
039c05d9f5 Added comments. 2012-03-07 09:59:04 +01:00
Bård Skaflestad
d12d444411 Include <cmath> for std::pow(). 2012-03-06 21:59:02 +01:00
Atgeirr Flø Rasmussen
ab324ed02e Merged. 2012-03-06 16:09:27 +01:00
Atgeirr Flø Rasmussen
f40aab431b Added polymerUtilitites.hpp and .cpp, implementing polymer-dependent total mobilities etc. 2012-03-06 16:07:27 +01:00
Atgeirr Flø Rasmussen
1919284446 Added effectiveInvVisc() method. 2012-03-06 16:06:00 +01:00
Xavier Raynaud
c3e89e6b5c Added "not so successfull" case in Newton method. 2012-03-06 15:25:37 +01:00
Xavier Raynaud
196c29522d Introduced ResidualEquation class to compute value of residual and derivatives. 2012-03-06 14:58:19 +01:00
Xavier Raynaud
e11cea7432 Set tolerance Newton solver according to tolerance in TransportPolymer. 2012-03-06 09:31:11 +01:00
Xavier Raynaud
801f975b5b Implemented case det(jacobian)=0 in Newton method. 2012-03-06 09:09:37 +01:00
Xavier Raynaud
02296478b1 Minor code optimization. 2012-03-06 09:00:05 +01:00
Xavier Raynaud
93503662ce Fixed bug for Newton method. Added initial residual check for Bracketing method. 2012-03-05 15:32:22 +01:00
Xavier Raynaud
a91b2d991d Use tolerance consistantly. 2012-03-05 12:58:09 +01:00
Xavier Raynaud
f783f9eebe Removed unused code. 2012-03-05 12:50:13 +01:00
Xavier Raynaud
c3d6cc429a Ensure min saturation is max(dead pore space, connate water saturation). 2012-03-05 12:44:31 +01:00
Xavier Raynaud
ceef4bbdcd Parameter single_cell_method is now enum. Some optimization. 
Redundant computation of fractional flow and mc removed.
 2012-03-05 11:56:23 +01:00
Atgeirr Flø Rasmussen
8877fa6ed0 Untabify. 2012-03-02 15:14:29 +01:00
Atgeirr Flø Rasmussen
ea1bd3ec2d Cleaned up PolymerProperties class a little. 2012-03-02 15:14:14 +01:00
Xavier Raynaud
f78730cb94 Cleaned up code for Splitting method to find zero of s-c residuals. 2012-02-29 15:56:11 +01:00
Xavier Raynaud
a33e3d9db0 rock density is obtained from PLYROCK in PolymerProperties::getFromDeck. 2012-02-29 14:58:51 +01:00
Xavier Raynaud
07f136befb Transport solver reads polymer data from PolymerProperties. 2012-02-29 13:27:20 +01:00
Xavier Raynaud
5ef45c59e3 Cleaned up PolymerProperties. 2012-02-29 13:25:51 +01:00
Xavier Raynaud
123762caa9 Added class PolymerProperties with a constructor which takes argument from an eclipse parser. 2012-02-29 11:27:59 +01:00
Xavier Raynaud
b0fdc4db7d Added Newton step as first step in Splitting s-c residual solver. 2012-02-28 17:36:29 +01:00
Xavier Raynaud
ee52e354b1 Bugfixed in computation of the gradient of the residual (from formulae). 2012-02-27 15:49:04 +01:00
Xavier Raynaud
5c78dc8960 Cleaned up code for Splitting method for residual zero finding. 2012-02-27 10:39:21 +01:00
Xavier Raynaud
a26656c6e5 Cleaned up code for Splitting method in zero finding residual. 2012-02-27 10:10:32 +01:00
Atgeirr Flø Rasmussen
753ebedc0d Removed "hello polymer" test program. 2012-02-26 20:20:37 +01:00
Atgeirr Flø Rasmussen
40e2ccd1f0 Removed legacy reorder solver code. 2012-02-26 20:17:22 +01:00
Xavier Raynaud
c63d817332 Added varying bounded box for allowable values of c and c in the splitting residual solver. 2012-02-24 17:35:47 +01:00
Xavier Raynaud
f0fc7bf3c0 Merge. 2012-02-23 17:00:43 +01:00
Xavier Raynaud
266b451715 Fix splitting method residual solver. Add piecewise linear curve to handle the boundaries of the acceptable domain for s and c. 2012-02-23 16:59:17 +01:00
Atgeirr Flø Rasmussen
95618aecac Implemented solveMultiCell() by the nonlinear Gauss-Seidel method. 2012-02-22 14:55:34 +01:00
Atgeirr Flø Rasmussen
21f7022afd Tolerance and max iters are constructor arguments. Uses satRange() properly. 2012-02-20 17:07:42 +01:00
Xavier Raynaud
1763e8afd7 Merged. 2012-02-20 09:37:54 +01:00
Xavier Raynaud
dd324478de Added s-c splitting solver for single cell problem. 2012-02-20 09:27:22 +01:00
Atgeirr Flø Rasmussen
d041c98f21 Added solveMultiCell() method to match change in interface. Currently throwing. 2012-02-15 22:52:32 +01:00
Atgeirr Flø Rasmussen
a3dfdad3a6 Use PolymerData struct as moved to namespace Opm and defined in TransportModelPolymer.hpp. 2012-02-12 23:45:43 +01:00
Atgeirr Flø Rasmussen
8cb2af77e7 Added TransportModelPolymer class. 2012-02-12 23:30:57 +01:00
Atgeirr Flø Rasmussen
094ede0855 Set minimum s to 0.2 for s-iterations (temporary measure). 2012-02-09 09:05:49 +01:00
Atgeirr Flø Rasmussen
240d69da57 Separated nonlinear controls for c and s residuals, set bracket fields. 2012-02-08 14:02:46 +01:00
Atgeirr Flø Rasmussen
46f29380b3 Added more (inactive) debug output. 2012-02-08 13:43:22 +01:00
Atgeirr Flø Rasmussen
37ce9d0db3 Improved treatment of inflow (source) polymer concentration: 
- solver now takes parameter for inflow_c
 - simulator calls polymerInflowAtTime(double time) for its value
 - currently starts polymer injection after 4 days (0.4 pv)
Also added some (inactive) debug printouts.
Still not working when injecting polymer.
 2012-02-07 09:09:34 +01:00
Atgeirr Flø Rasmussen
d17cfd1372 First compiling version of full polymer residual. Not yet working. 2012-02-06 14:01:56 +01:00
Atgeirr Flø Rasmussen
0ad4fab860 Minor bugfix, now reproducing non-polymer results. 2012-02-06 10:45:54 +01:00
Atgeirr Flø Rasmussen
f46587bf3c Structure for c-solve in place, still not implemented residual_c(). 2012-02-06 08:54:38 +01:00
Atgeirr Flø Rasmussen
143861745a Initial work on supporting polymer transport. Work in progress. 2012-02-03 16:39:55 +01:00
Atgeirr Flø Rasmussen
e1df60c04d Copied and renamed files and functions from opm-core as base for polymer reorder solver. 2012-02-02 16:35:55 +01:00
Bård Skaflestad
08671a1d6d Coalesce constructors using a default value. 2012-02-01 16:40:15 +01:00
Bård Skaflestad
a785d9488e Create minimum buildable Polymer module for OPM. 2012-02-01 16:16:30 +01:00