Commit Graph

333 Commits

Author SHA1 Message Date
Bård Skaflestad
afeda5a95a Use pure htrans strategy for computing fluxes/if-pressures.
Still wrong.  Wrong: R-O-N-G.  Wrong.
2010-11-11 14:47:17 +01:00
Bård Skaflestad
112c4503f3 Merge from upstream. 2010-11-11 09:41:16 +01:00
Atgeirr Flø Rasmussen
777f3129e0 Updated after nf parameter was removed. 2010-11-11 09:16:49 +01:00
Bård Skaflestad
1ffc7449f0 Compute real Darcy flux in cfs_tpfa_press_flux().
The previous edition made the grave mistake of summing the mass fluxes
per face.  This does not make sense.

Pointed out by: Jostein.R.Natvig@sintef.no
2010-11-10 14:31:30 +01:00
Bård Skaflestad
5600f8d505 Add flux and face-pressure computations.
Specifically, introduce utilities compute_fpress() and
compute_flux().  The former is needed to implement the latter across
external boundary faces.  Moreover, interface pressure values are
needed to evaluate fluid properties on faces (specifically R/B).

Add small gateway routine, cfs_tpfa_fpress(), to allow callers to
recover interface pressure values.  Re-implement cfs_tpfa_press_flux()
in terms of compute_fpress() and compute_flux().

Also, add fields 'fpress' and 'fpaccum' to struct cfs_tpfa_impl.
2010-11-09 12:52:31 +01:00
Bård Skaflestad
7d668267cf Remove unused variable. 2010-11-09 10:03:42 +01:00
Bård Skaflestad
d7a80b1082 Merge from jrn. 2010-11-08 14:13:09 +01:00
Bård Skaflestad
fbed044c94 Don't enforce p[0]=0 for non-Neumann problems. 2010-11-06 19:16:23 +01:00
Jostein R. Natvig
43ec21e382 Remove unnecessary argument in mim_ip_simple_all; fix affected code;
as suggested by bska.
2010-11-05 22:05:55 +01:00
Bård Skaflestad
223681b8c4 Use a 'P' sign that is somehow related to the equation we're solving
The former definition was inherited from the hybrid solver in MRST and
not really applicable to compressible TPFA.
2010-11-05 16:58:19 +01:00
Bård Skaflestad
8780df0a55 Tentatively add boundary condition support.
We currently handle only Dirichlet (pressure) BC's, and the feature is
largely untested.  More work remaining.
2010-11-04 19:05:41 +01:00
Bård Skaflestad
7b692205ae Merge from upstream sources. 2010-11-03 18:50:45 +01:00
Bård Skaflestad
08ae534fa5 Correct silly mistake in defining hface phase contribs.
We were only advancing the phase contrib once per cell, but need a
different contribution for each face for each cell...
2010-11-03 18:08:02 +01:00
Bård Skaflestad
1640925cd1 Move compressible support utilities to cfs_tpfa module.
We may wish to generalise this part into real utility functions but
for now, leave the functionality where it is most directly needed.
Rename the functions to highlight relationship to compr. tpfa.
2010-11-03 14:55:23 +01:00
Bård Skaflestad
a0208b30dc Add compressible TPFA solver to build. 2010-11-03 14:50:26 +01:00
Atgeirr Flø Rasmussen
d6fd317c66 Added code to check the grid adapter, normals in C now area-weighted. 2010-11-02 14:25:29 +01:00
Atgeirr Flø Rasmussen
d8d0693f62 Added (inactive) test code to check the behaviour of the grid adapter. 2010-11-02 14:24:29 +01:00
Bård Skaflestad
83c9d73252 Add files for implementing compressible TPFA assembly. 2010-11-01 18:45:06 +01:00
Bård Skaflestad
601d5c9099 Add files for implementing compressible TPFA assembly. 2010-11-01 18:44:40 +01:00
Bård Skaflestad
220d3f36ab Request Fortran compiler support to accomodate LAPACK checks. 2010-11-01 14:39:48 +01:00
Bård Skaflestad
8cd3579ab1 Request "all" Automake warnings. 2010-11-01 14:39:03 +01:00
Bård Skaflestad
dbb2111f30 Distribute and install TPFAPressureSolver.hpp . 2010-11-01 14:26:18 +01:00
Bård Skaflestad
aa97ec444b Don't pretend to expand SVN keywords. 2010-11-01 13:33:54 +01:00
Bård Skaflestad
3f8ddf9699 Merge from Atgeirr. 2010-11-01 13:30:49 +01:00
Bård Skaflestad
a70a97f003 Don't use uninitialised variable. 2010-11-01 10:07:49 +01:00
Bård Skaflestad
a13c9cd0d9 Move source files to sub-dir 'src'. 2010-10-29 15:08:09 +02:00
Atgeirr Flø Rasmussen
11c67f8c0c Removed unneeded member variables. 2010-10-29 14:57:11 +02:00
Bård Skaflestad
288062511e Implement tools for computing accumulation term. 2010-10-29 14:44:10 +02:00
Bård Skaflestad
7f6646677e Build a second, no-inst, edition of lib. for inclusion purposes. 2010-10-29 14:26:57 +02:00
Bård Skaflestad
cfe50e70ad Move factoring matrix closer to solving systems.
Factoring the matrix and building the right-hand sides are
order-independent.  Thus, there is no cost in moving the factoring
closer to its use in actually solving the systems.
2010-10-29 11:17:30 +02:00
Bård Skaflestad
9205fc318e Re-purpose tpfa_compr_htran_mult to computing compr. htrans.
In other words, assume that transmissibility effects are already
included in the 'xf' parameter.  Rename tpfa_compr_htran_mult() to
tpfa_compr_htrans() to reflect new purpose.
2010-10-29 11:07:37 +02:00
Bård Skaflestad
7c7cc5cffa Ipiv is output from DGETRF. No need to initialise. 2010-10-28 14:22:55 +02:00
Bård Skaflestad
e6958e5ac0 Fix another fallout from renaming invAc -> luAc. 2010-10-28 14:03:53 +02:00
Bård Skaflestad
560e281738 Merge from Atgeirr. 2010-10-28 14:01:17 +02:00
Atgeirr Flø Rasmussen
050312f1b2 Modified comments. 2010-10-28 11:56:07 +02:00
Bård Skaflestad
ee05efd2cb Fix spelling error left over in cset df6e16aa6dd5.
Pointy hat: bska
2010-10-28 11:56:07 +02:00
Bård Skaflestad
677329acf6 Add low-level utilities for computing compr. htran multipliers. 2010-10-28 11:53:40 +02:00
Atgeirr Flø Rasmussen
95808ecee7 Added C++ interface to incompressible tpfa solver. 2010-10-28 11:39:56 +02:00
Atgeirr Flø Rasmussen
6dacf10b60 Added extern "C" clauses. 2010-10-28 11:39:23 +02:00
Atgeirr Flø Rasmussen
2d912ff0a9 Set the 'n' member of the matrix (number of columns). 2010-10-28 11:38:41 +02:00
Bård Skaflestad
38c8de2f4d Declare DGETRF and DGETRS for compressible support. 2010-10-28 10:51:59 +02:00
Bård Skaflestad
3c2f060b04 Add simple utility function for computing many small matvecs.
This is simply the obvious loop around DGEMV.  We may wish to move it
elsewhere if need arises, but for now the function finds utility only
in computing compressible transmissibility multipliers.
2010-10-28 10:07:19 +02:00
Bård Skaflestad
667eefad0b M-x delete-trailing-whitespace . 2010-10-26 22:30:54 +02:00
Bård Skaflestad
6f003d091c Kill EOF blank lines. 2010-10-26 22:26:25 +02:00
Bård Skaflestad
dcb548c07f Update ifsh_assemble() comment to reflect effective terms. 2010-10-26 18:39:06 +02:00
Atgeirr Flø Rasmussen
56822daebb Added implementation header to SOURCES variable as recommended by automake. 2010-10-26 14:37:47 +02:00
Atgeirr Flø Rasmussen
1add223064 Moved and removed files. 2010-10-26 14:19:07 +02:00
Atgeirr Flø Rasmussen
411dafca63 Fixed for interface changes. 2010-10-26 13:58:38 +02:00
Atgeirr Flø Rasmussen
474c2ef855 Transformed fsh->cfsh when meaning compressible flow solver hybrid. 2010-10-26 13:54:49 +02:00
Bård Skaflestad
e073bfa081 Merge from Atgeirr. 2010-10-26 13:16:27 +02:00