Commit Graph
11092 Commits
Author SHA1 Message Date
Bård Skaflestad 2ad508d58b Expand flux calculation interface to account for presence of wells. 2011-11-25 15:27:08 +01:00
Bård Skaflestad a11c59f8dd Really remove all purported boundary condition handling. 2011-11-25 14:19:12 +01:00
Bård Skaflestad b5ebb988da Don't pretend to support boundary conditions.
We don't.
2011-11-25 14:17:07 +01:00
Bård Skaflestad 634d168b3f Assemble well equation form contributions per completion.
We assemble one equation for each well, with the equation being
trivial for wells controlled by bottom-hole targets.  Moreover, at
present we only support BHP and RESV controls.  The infrastructure is
in place to support more general controls (i.e., surface rate
targets), but such wells are currently not part of the 'well.h'
module.
2011-11-25 14:02:15 +01:00
Bård Skaflestad f45f6ea62c Prepare to assemble RESV control target. 2011-11-24 23:21:41 +01:00
Bård Skaflestad 1d4d16feab Assemble completion-to-cell contributions. Reestablish build. 2011-11-24 22:39:17 +01:00
Bård Skaflestad fcbc068091 Check-point commit to mark half-way support for well assembly.
Does not build.
2011-11-24 18:57:38 +01:00
Bård Skaflestad ffd1531f2e Implement csrmatrix_zero() in terms of vector_zero(). 2011-11-23 23:06:16 +01:00
Bård Skaflestad 4a158b89a9 Compute component fluxes across completions if the model includes wells.
This necessitates a public interface change whereby the caller is
required to pass the current well (bottom-hole) pressure values into
the assembler.
2011-11-23 19:25:47 +01:00
Bård Skaflestad 6ab5800483 Allocate backing store for well completion data.
Also, include well connections in system matrix.
2011-11-23 15:25:48 +01:00
Bård Skaflestad a1df963169 Take initial steps towards including wells in assembly.
Specifically, accept (and currently ignore), a WellCompletions
structure into the constructor, and aggregate all driving forces
(source terms, boundary conditions and all well-related structures)
into an assembler-specific "force" structure.  Accept a pointer to
such a structure into the assemble() function.  Currently ignored
except for source terms.
2011-11-23 13:20:37 +01:00
Bård Skaflestad 9a9b47a49a Name the well_t and well_control_t structures. 2011-11-23 11:23:30 +01:00
Bård Skaflestad 51111176db Guard against 'bc' being NULL.
This is in preparation of introducing compressible boundary conditions
as represented by the 'compr_bc' module.
2011-10-25 19:34:32 +02:00
Bård Skaflestad d1026c82d5 Don't use uninitialised variables. 2011-10-21 15:30:44 +02:00
Bård Skaflestad afdb88678f Fix remaining bugs for incompressible flows driven by sources.
In particular:

  - Zero sufficient portion of the ->mat_row to hold all connections
    of a cell in addition to the accumulation term.

  - Don't write the residual into the accumulation term of the
    Jacobian matrix row.
2011-10-20 13:59:15 +02:00
Bård Skaflestad 9392b3700f Don't pretend that the time-step is always one. 2011-10-19 22:28:41 +02:00
Bård Skaflestad 98abee6a79 Account for flow-direction signs when accumulating fluxes. 2011-10-19 22:10:22 +02:00
Bård Skaflestad fb1316802c Implement compressible (volume flux) source term assembly. 2011-10-19 20:09:51 +02:00
Bård Skaflestad 79b5e94127 Accept a set of compressible source terms rather than a double*.
While here, restore an essential semicolon lost in cset 11881e6cd650.
2011-10-19 19:54:17 +02:00
Bård Skaflestad 5fea1f8556 Remove unused function. 2011-10-19 19:48:46 +02:00
Bård Skaflestad 14533690e4 Add a simple representation of a compressible (volume flux) source term. 2011-10-19 19:44:32 +02:00
Bård Skaflestad 7f40b2a4e7 Clear Jacobian row before assembling local contributions. 2011-10-19 15:14:11 +02:00
Bård Skaflestad 3cf841fc71 Simplify offset calculation. 2011-10-19 15:12:16 +02:00
Bård Skaflestad 4163ea49ee Correct various indexing errors. 2011-10-19 15:11:21 +02:00
Bård Skaflestad 501b04317d Re-install interface pressure calculation in residual formulation. 2011-10-19 13:27:24 +02:00
Bård Skaflestad 2a833768b2 Remove further traces of well support.
Wells will be introduced at a later time.
2011-10-19 13:11:44 +02:00
Bård Skaflestad 6acd22f7d7 Hook residual formulation up to build. 2011-10-19 13:01:06 +02:00
Bård Skaflestad c8d7ca7fbd Merge cfs_tpfa backout. This restores backwards compatibility. 2011-10-18 23:22:52 +02:00
Bård Skaflestad 151adc6e3f Backed out changeset d83710dd6839 2011-10-18 23:22:26 +02:00
Bård Skaflestad 7d2f0b2b71 Merge 'compr_quant' backout.
This is another attempt at restoring backwards compatibility.
2011-10-18 23:20:57 +02:00
Bård Skaflestad f55b9f5412 Backed out changeset 9faef9e37070 2011-10-18 23:20:34 +02:00
Bård Skaflestad 62e73a8a67 Grab copy of 'compr_quant' module. 2011-10-18 23:19:53 +02:00
Bård Skaflestad 4e77f6cd56 Merge backout. 2011-10-18 23:18:29 +02:00
Bård Skaflestad d909e0f6d5 Backed out changeset 3b7e5d602aab
We need a copy of 'compr_quant' before restoring compat.
2011-10-18 23:18:22 +02:00
Bård Skaflestad d64bdc8997 Merge backout. 2011-10-18 23:17:33 +02:00
Bård Skaflestad 2134f1f91e Backed out changeset 6b1a90716ea5 2011-10-18 23:17:02 +02:00
Bård Skaflestad 000f71e12c Back out cset 5e3d75476d64. 2011-10-18 23:14:15 +02:00
Bård Skaflestad 704d7eab5c Backed out changeset 5e3d75476d64
In preparation of restoring backwards compatibility in 'cfs_tpfa' module.
2011-10-18 23:13:41 +02:00
Bård Skaflestad e81a424995 Grab copy of 'cfs_tpfa' module in preparation of restoring backwards compat. 2011-10-18 23:12:54 +02:00
Bård Skaflestad dd042b75a8 Make first attempt at implementing the residual/Jacobian formulation.
This is a step in the direction of having a true IMPES pressure solver
and to remove the 'experimental_jacobian' option in dune-porsol.
2011-10-18 20:46:05 +02:00
Bård Skaflestad 7fb517e192 Delete trailing whitespace. 2011-10-18 20:40:19 +02:00
Bård Skaflestad 94fc6de78a Rework compressibility representation.
Switch to storing a complete fluid-matrix derivative in the
compr_quantities rather than the total compressibility.  Maintain the
"volume discrepancy" field.  Also, add traditional memory management
functions.
2011-10-18 20:38:28 +02:00
Bård Skaflestad c3036166a9 Factor pressure (increment) assignment out of _press_flux().
The linear solution h->x is the pressure increment, not the actual
pressure value, so we cannot compute fluxes based on h->x alone.
2011-10-17 11:05:04 +02:00
Bård Skaflestad be6ae76d88 Don't reference (non-existent) Dune directories. Remove disabled code. 2011-10-13 19:38:33 +02:00
Halvor M. Nilsen e7d5d4bb27 added EulerUpstreamImplicit, added implicit traits: does notthing now 2011-10-11 11:42:40 +02:00
Halvor M. Nilsen b42bd08154 Move 'common' directory into 'dune'. 2011-10-07 10:54:25 +02:00
Bård Skaflestad 5ee2c381c3 Import glue code for calling into UMFPACK using an OPMPressure CSRMatrix. 2011-10-05 20:03:07 +02:00
Bård Skaflestad 01eb7a2214 Add an overload suited to pointer arguments (primitive types). 2011-10-05 19:54:12 +02:00
Bård Skaflestad 4f488a08e7 Work around different 'const' rules in C. 2011-10-04 21:21:30 +02:00
Bård Skaflestad 55cbd9143b Don't reference unneeded interfaces. 2011-10-04 10:10:44 +02:00