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
Kjetil Olsen Lye
d1e8525c4d
merge
2012-02-27 17:47:55 +01:00
Kjetil Olsen Lye
c06401ef5b
Fixed some warnings
2012-02-27 17:46:06 +01:00
Bård Skaflestad
a0598eae4a
Expose internal half-transmissibility through read-only reference.
...
Useful to initialise the Opm::ImplicitTransport<> solver with gravity
effects.
2012-02-27 17:29:09 +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
Bård Skaflestad
f4df7aa1bd
Don't assume that num_phases==2 when calculating total mobility.
2012-02-27 10:07:00 +01:00
Atgeirr Flø Rasmussen
b058678235
Moved simulator utility functions to opm/core/utility/miscUtilities.hpp.
...
We want to avoid having hundreds of files with a single function, therefore
these functions have been lumped together. With time, a more sophisticated
organization of such free functions may be in order.
2012-02-26 21:05:19 +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
Atgeirr Flø Rasmussen
85daef9b5b
Added experimental code guarded by EXPERIMENT_GAUSS_SEIDEL #define.
2012-02-26 00:30:36 +01:00
Atgeirr Flø Rasmussen
3312aad6f5
Adding data for experiments, left alone if unused.
2012-02-26 00:29:39 +01:00
Atgeirr Flø Rasmussen
9667d09fe0
Added more inactive experimental code.
2012-02-25 22:39:01 +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
8c16722b84
Removed copied text accidentally left in.
2012-02-21 22:02:16 +01:00
Atgeirr Flø Rasmussen
9ed9f8b288
Fixed class comments, LinearSolverUmfpack::solve() properly returns a report.
2012-02-21 21:54:46 +01:00
Atgeirr Flø Rasmussen
87723a08a0
IncompTpfa now takes a linear solver as a constructor argument.
2012-02-21 21:45:04 +01:00
Atgeirr Flø Rasmussen
e52cae02ac
Made the solve() methods const.
2012-02-21 21:36:09 +01:00
Atgeirr Flø Rasmussen
8f8baa62c8
Using LinearSolverUmfpack for the linear solver.
2012-02-21 21:32:11 +01:00
Atgeirr Flø Rasmussen
3cb3d69d90
Added LinearSolverInterface and two subclasses, using Umfpack and Istl.
2012-02-21 21:27:15 +01:00
Atgeirr Flø Rasmussen
21f7022afd
Tolerance and max iters are constructor arguments. Uses satRange() properly.
2012-02-20 17:07:42 +01:00
Atgeirr Flø Rasmussen
22bf83dda7
Bugfix: must always initialize gpress_omegaweighted_, even if no gravity.
2012-02-20 17:05:28 +01:00
Atgeirr Flø Rasmussen
2aac2255b5
Ensuring gravity feature is used consistently.
2012-02-20 13:39:40 +01:00
Atgeirr Flø Rasmussen
e5a7157bc3
Added class IncompTpfa.
2012-02-20 13:23:01 +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
f9fd79dbcb
Moved writeVtkData() functions to utility/, added a test program.
2012-02-19 21:24:59 +01:00
Atgeirr Flø Rasmussen
16378af014
Added tol and maxit arguments. Using Incomp*::satRange() for initial s-bracket.
2012-02-17 14:17:16 +01:00
Atgeirr Flø Rasmussen
9901d20721
Added satRange() method to Incomp* classes and Saturation* classes.
...
So far not yet added to BlackoilPropertiesInterface.
2012-02-17 14:15:25 +01:00
Atgeirr Flø Rasmussen
75731efa9a
Added lots of (inactive) experimental code.
...
The code attempts to improve #iterations in the Gauss-Seidel-like
multicell solver by improving ordering. In general, experiment failed
to improve #iterations, except for one: totally random order was the
best (for the 100x100 case tried)!
2012-02-17 09:39:45 +01:00
Atgeirr Flø Rasmussen
b5eb8543aa
Added compute_sequence_graph(), which also returns the upwind graph.
2012-02-17 09:29:56 +01:00
Xavier Raynaud
Atgeirr Flø Rasmussen
Kjetil Olsen Lye
Bård Skaflestad