Jacques Franc

Added validation example for Problem 3.x of Class (2009)

1

This PR is following #1723 in adding a validation case on Johansen mesh as proposed in the third problem (Problem 3.x) of the benchmark paper: ["A benchmark study on problems...

Companion PR with https://github.com/GEOSX/thirdPartyLibs/pull/193 adding a target for Sherlock Stanford cluster configuration: (CPU-version) - centos 7.9.2009 - gcc-10.1.0 - openmpi 4.1.2 - openblas 0.3.10 (GPU-version) - cuda11.5.0 [Debug only -...

[Bug] Failed compilation w/ Cuda 11.5 in `benchmarkOuterProduct` from LvArray on Sherlock

3

**Describe the bug** Compilation error while using compiling cuda-on version of GEOSX in *Sherlock* docker replica. related to #2010 and https://github.com/GEOSX/thirdPartyLibs/pull/193 **To Reproduce** 1. Pull docker image https://hub.docker.com/r/geosx/sherlock-gcc10.1.0-openmpi4.1.2-cuda11.5.0-openblas0.3.10-zlib1.2.11 2. Try...

Making relative permeabilities tensorial

1

## Aim This PR aims at making relative permeabilities (aka relperm) a diagonal tensor, so one can vary their characteristics (exponents, max values for Brooks and Corey(BC) or Van Genuchten(VG))...

Refactoring flow `StackVariables` to use `DerivOffset`

1

As it is used in `isothermalCompositionalMultiphaseFVMKernel` or `MultifluidBase`, derivative offset helps readability and would made extension easier. The refactoring though left out for now - `dTrans_dDisJump` as this derivative is...

Refactoring diffusion and tpfa weights calculation to enable phase indexed arrays (4d)

5

Like what was discuss [here](https://github.com/GEOS-DEV/GEOS/pull/2678#discussion_r1329317111) about allowing phase-indexed tensor in the tpfa transmissivities/weights calculation to be able to keep consistency with usual transmissivity computed on intrinsinc permeability. This is also...

Adding dispersion model for SPE11

1

Continuing the work initiated in #2678, the `LinearIsotropicDispersion` constitutive class init the `BaseDispersion`'s `m_dispersivity` cell-wise diagonal tensor field to the sole longitudinal dispersion coefficient. In order to get a cell-centered...

This would allow extension to the function triggered logic to `HaltEvent`. Practical use is _e.g._ to stop simulation and trigger vtk output as the CO2 front reach a certain cell...

Producing residual maps

4

This PR intends at implementing a mechanism to produce residual maps in addition to other fields produce by the `CompositionalMultiphaseFVM` solver. The desirable features are: - Storing residual systematically at...