Dear ogs users,
I want to model water flow and mass transport in a dual-porosity 2d domain using ogs5 and DUAL_RICHARDS flow.
I set up a simple model similar to the benchmark “H_us/dual_van” from the 2012 benchmark book “THMC processes in fractured porous media”/chapter 6.2. with one material property bloc for the matrix , one for the fracture (*.mmp file, 1. example see below). It works.
However, my real model domain contains three materials which you usually define in the *.mmp file when using single RICHARDS_FLOW (2. example see below).
My Question: How do you define the matrix and the fracture part for all three materials? Is that possible at all? Just putting all material blocs twice (alternating matrix/fracture and/or all matrix then all fracture) did not work (see example 3./4.).
- MMP Example: simple DUAL_RICHARDS model with 1 matrix and 1 fracture:
// Matrix
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.3
$PERMEABILITY_TENSOR
ISOTROPIC 2.86e-11
$PERMEABILITY_SATURATION
4 0.0 0.93 0.44
$CAPILLARY_PRESSURE
4 4.5
$MASS_DISPERSION
1 0.8 0.2
$DENSITY
1 1800.0
$TORTUOSITY
1 1.0
$TRANSFER_COEFFICIENT
1.e+2
// Fracture
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.38
$PERMEABILITY_TENSOR
ISOTROPIC 1.0e-8
$PERMEABILITY_SATURATION
4 0.1 0.93 0.44
$CAPILLARY_PRESSURE
4 810
$MASS_DISPERSION
1 0.8 0.2
$DENSITY
1 1800.0
$TORTUOSITY
1 1.0
#STOP
- MMP Example: single RICHARDS_FLOW with 3 materials in domain
// inlet
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.38
$PERMEABILITY_TENSOR
ISOTROPIC 1.0e-07
$PERMEABILITY_SATURATION
4 0.01 0.99 0.78
$CAPILLARY_PRESSURE
4 1780
$DENSITY
1 1800.0
// main
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.38
$PERMEABILITY_TENSOR
ISOTROPIC 1.0e-8
$PERMEABILITY_SATURATION
4 0.1 0.93 0.44
$CAPILLARY_PRESSURE
4 810
$DENSITY
1 1800.0
// outlet
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.38
$PERMEABILITY_TENSOR
ISOTROPIC 1.0e-07
$PERMEABILITY_SATURATION
4 0.01 0.99 0.78
$CAPILLARY_PRESSURE
4 1780
$DENSITY
1 1800.0
#STOP
- Example DUAL_RICHARDS for three materials
// DUAL_RICHARDS_3materials
// Inlet -matrix
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.3
$PERMEABILITY_TENSOR
ISOTROPIC 2.86e-11
$PERMEABILITY_SATURATION
4 0.0 0.93 0.44
$CAPILLARY_PRESSURE
4 4.5
$DENSITY
1 1800.0
$TRANSFER_COEFFICIENT
1.e+2
// Inlet -fracture
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.38
$PERMEABILITY_TENSOR
ISOTROPIC 1.0e-07
$PERMEABILITY_SATURATION
4 0.01 0.99 0.78
$CAPILLARY_PRESSURE
4 1780
$DENSITY
1 1800.0
// main - matrix
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.3
$PERMEABILITY_TENSOR
ISOTROPIC 2.86e-11
$PERMEABILITY_SATURATION
4 0.0 0.93 0.44
$CAPILLARY_PRESSURE
4 4.5
$DENSITY
1 1800.0
$TRANSFER_COEFFICIENT
1.e+2
// main - fracture
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.38
$PERMEABILITY_TENSOR
ISOTROPIC 1.0e-8
$PERMEABILITY_SATURATION
4 0.1 0.93 0.44
$CAPILLARY_PRESSURE
4 810
$DENSITY
1 1800.0
// outlet - matix
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.3
$PERMEABILITY_TENSOR
ISOTROPIC 2.86e-11
$PERMEABILITY_SATURATION
4 0.0 0.93 0.44
$CAPILLARY_PRESSURE
4 4.5
$DENSITY
1 1800.0
$TRANSFER_COEFFICIENT
1.e+2
// outlet - fracture
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.38
$PERMEABILITY_TENSOR
ISOTROPIC 1.0e-07
$PERMEABILITY_SATURATION
4 0.01 0.99 0.78
$CAPILLARY_PRESSURE
4 1780
$DENSITY
1 1800.0
#STOP
- Example
// DUAL_RICHARDS_3materials
// Inlet -matrix
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.3
$PERMEABILITY_TENSOR
ISOTROPIC 2.86e-11
$PERMEABILITY_SATURATION
4 0.0 0.93 0.44
$CAPILLARY_PRESSURE
4 4.5
$DENSITY
1 1800.0
$TRANSFER_COEFFICIENT
1.e+2
// main - matrix
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.3
$PERMEABILITY_TENSOR
ISOTROPIC 2.86e-11
$PERMEABILITY_SATURATION
4 0.0 0.93 0.44
$CAPILLARY_PRESSURE
4 4.5
$DENSITY
1 1800.0
$TRANSFER_COEFFICIENT
1.e+2
// outlet - matix
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.3
$PERMEABILITY_TENSOR
ISOTROPIC 2.86e-11
$PERMEABILITY_SATURATION
4 0.0 0.93 0.44
$CAPILLARY_PRESSURE
4 4.5
$DENSITY
1 1800.0
$TRANSFER_COEFFICIENT
1.e+2
// Inlet -fracture
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.38
$PERMEABILITY_TENSOR
ISOTROPIC 1.0e-07
$PERMEABILITY_SATURATION
4 0.01 0.99 0.78
$CAPILLARY_PRESSURE
4 1780
$DENSITY
1 1800.0
// main - fracture
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.38
$PERMEABILITY_TENSOR
ISOTROPIC 1.0e-8
$PERMEABILITY_SATURATION
4 0.1 0.93 0.44
$CAPILLARY_PRESSURE
4 810
$DENSITY
1 1800.0
// outlet - fracutre
#MEDIUM_PROPERTIES
$GEOMETRY_DIMENSION
2
$GEOMETRY_AREA
1.2
$POROSITY
1 0.38
$PERMEABILITY_TENSOR
ISOTROPIC 1.0e-07
$PERMEABILITY_SATURATION
4 0.01 0.99 0.78
$CAPILLARY_PRESSURE
4 1780
$DENSITY
1 1800.0
#STOP