Commit Graph

470 Commits

Author SHA1 Message Date
Bård Skaflestad
552bfd5f00 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
27ddb568ba Remove unused variable. 2010-11-09 10:03:42 +01:00
Bård Skaflestad
80f771d674 Merge from jrn. 2010-11-08 14:13:09 +01:00
Bård Skaflestad
43ddeaeca8 Don't enforce p[0]=0 for non-Neumann problems. 2010-11-06 19:16:23 +01:00
Jostein R. Natvig
ee6fe14e48 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
dba90f41bd 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
650de8a127 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
99383e8f02 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
Jostein R. Natvig
7459c5988e Add bare-bones explicit mobility-weighted upwind solver for two-phase
incompressible flow.
2010-11-02 09:10:36 +01:00
Bård Skaflestad
a24685af0a Add files for implementing compressible TPFA assembly. 2010-11-01 18:45:06 +01:00
Bård Skaflestad
5c6f149485 Add files for implementing compressible TPFA assembly. 2010-11-01 18:44:40 +01:00
Bård Skaflestad
ffaa7c4585 Move source files to sub-dir 'src'. 2010-10-29 15:08:09 +02:00
Atgeirr Flø Rasmussen
5274eafcc1 Added extern "C" clauses. 2010-10-28 11:39:23 +02:00
Atgeirr Flø Rasmussen
496bbb7faf Set the 'n' member of the matrix (number of columns). 2010-10-28 11:38:41 +02:00
Bård Skaflestad
a05de2d311 M-x delete-trailing-whitespace . 2010-10-26 22:30:54 +02:00
Bård Skaflestad
a13e6a63d5 Complete cset a68849062624 (Update comment...).
Pointy hat: bska
2010-10-26 12:35:35 +02:00
Bård Skaflestad
0633a430db Update comment to reflect reality. 2010-10-26 12:33:07 +02:00
Bård Skaflestad
c918c34e1f Don't update in place. 2010-10-26 12:28:31 +02:00
Bård Skaflestad
8d528e13d1 Implement simple gravity handling.
Produces expected results on standard verification test.
2010-10-26 11:57:00 +02:00
Bård Skaflestad
bdd538c764 Remove zero eigenval. Prepare for adding gravity. 2010-10-25 20:28:34 +02:00
Bård Skaflestad
7ee6eb7758 The ifs_tpfa_press_flux() function does not need 'src'. 2010-10-24 21:22:29 +02:00
Bård Skaflestad
d130aef655 Add (incomp) pressure system assembler based on TPFA. 2010-10-24 21:20:15 +02:00
Bård Skaflestad
4274f20c27 Add utilities for computing derived fluid quantities.
These are (currently) entirely geared towards pressure solvers.
2010-10-20 19:14:07 +02:00
Bård Skaflestad
664f1ac2fe Add facilities for computing/updating gpress/Binv.
We were already computing the inverse IP, but now centralise the
gpress as well.  Moreover, the mim_ip_*_update() functions will assist
in the pressure solvers accepting effective inner products and gravity
pressures only.
2010-10-14 14:20:11 +02:00
Bård Skaflestad
edd470dd7d Move implementation of vector_zero() to central location.
There is no need for each file to contain a separate, though
trivial, implementation of this feature.
2010-10-14 13:58:47 +02:00
Bård Skaflestad
ef4d884d12 Use canonical include guards.
Suggested by atgeirr.  Template:

   OPM_<FILENAME>_HEADER_INCLUDED
2010-10-13 18:35:15 +02:00
Atgeirr Flø Rasmussen
799dc936cd Made all C headers includeable from C++. 2010-10-12 07:44:02 +00:00
Atgeirr Flø Rasmussen
6530a16785 Added copyright block to all source code files. 2010-10-12 07:25:46 +00:00
Bård Skaflestad
35c0ef8ee9 Remove delusional assertion. While its predecessor asserted basic
sanity, the current check would never be satisfied in a reasonable
  grid.
2010-10-06 21:08:56 +00:00
Bård Skaflestad
4357d2564c Correct another off-by-one in the CSR start pointer implementation. 2010-10-04 17:43:45 +00:00
Atgeirr Flø Rasmussen
14f682a0b4 Documented mim_ip_simple_all. 2010-10-01 12:21:24 +00:00
Bård Skaflestad
284683b00d Add a routine for allocating a CSR matrix with known number of rows
and non-zeros.  Note that this routine is for allocation only.
  Caller must create and manage sparsity structure if this is being
  used in a global assembly process.

  This routine finds utility in MsMFEM BF construction.
2010-09-30 21:37:14 +00:00
Atgeirr Flø Rasmussen
1e2e18bda2 Added subrepo for libmimetic. It is currently on internal server. 2010-09-27 07:48:31 +00:00
Atgeirr Flø Rasmussen
767b9f3b8b Added extern "C" statements for C++ compatibility. 2010-09-24 09:07:50 +00:00
Bård Skaflestad
f3e5ad1723 Use M's notion of the MAT_SIZE_T when building for M. 2010-09-17 17:43:07 +00:00
Bård Skaflestad
fca8a08b29 Separate well connections (typically static) from well controls
(often dynamically changing throughout simulation).
2010-09-15 17:27:59 +00:00
Bård Skaflestad
222d6da348 Re-factor creation and destruction of cell->well mapping out to
separate module.  Update callers.
2010-09-15 12:13:53 +00:00
Bård Skaflestad
31c5ff58b5 Don't state that we reverse-insert (i.e., push-forward) the bin
elements when we in fact use jrn's much more straight-forward
  push-back mode.
2010-09-14 12:21:30 +00:00
Bård Skaflestad
293e51f94b Insert "H1 lines" to summarise purpose. Document a few parameters
and return values as well.

  While here, convert CSR build to push-back mode.
2010-09-14 11:19:08 +00:00
Bård Skaflestad
1540fb3949 Remove special MATLAB compilation support. This is provided by blas_lapack.h. 2010-09-13 16:11:49 +00:00
Bård Skaflestad
5e3dc1daea Add simple type for describing boundary conditions for the flow
equation.
2010-09-13 15:46:06 +00:00
Bård Skaflestad
371fd66ced Add initial cut at well definition. 2010-09-10 17:14:00 +00:00
Bård Skaflestad
6e330e6e98 Add simple support routines for building CSR matrices. 2010-09-09 16:48:03 +00:00
Bård Skaflestad
2bced196f6 Move solvers/common solvers/euler and solvers/mimetic to dune-porsol 2010-08-23 09:23:09 +00:00
Bård Skaflestad
3e76680f4d Silence compiler by adding prototype for dfs() (i.e., #include
"dfs.h").  Also, don't require build(er|system) to always #define
  the preprocessing symbol "TEST" (to a zero or non-zero value).
2010-08-20 21:44:49 +00:00
Bård Skaflestad
a32d350270 Assert copyright for 2010. 2010-08-20 19:16:59 +00:00
Bård Skaflestad
7bcdf7363e M-x delete-trailing-whitespace . 2010-08-20 19:05:27 +00:00
Bård Skaflestad
7976cffaf3 Correct a supremely embarrassing typo: When counting the number of
cell-to-cell connections (i.e., when calling count_block_conns()),
  inspect the pointer array (pc2c) rather than the cell array itself
  (c2c).  Clue-by-four supplied by gdb(1).

  Move self connection handling to partition_create_c2c().

  Adjust dfs() 'work' array size to comply with new requirements of
  -r5088.
2010-08-20 17:24:26 +00:00
Jostein R. Natvig
d31614d608 Fix embarrasing error. 2010-08-20 14:18:22 +00:00
Bård Skaflestad
d2c603b980 Add files missing in -c5083.
Pointy hat: bska
2010-08-19 22:30:54 +00:00