ogs5 -- Dual-porosity -- DUAL_RICHARDS -- Multiple Materials

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.).

  1. 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

  1. 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

  1. 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

  1. 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

Hi Johannes,

  the materials have to be present in the mesh, i.e. you have to assign the materials to the corresponding parts of the mesh (material properties are assigned to the elements). Please have a look how the mesh file looks in the H_us/dual_van benchmark.

CU, Georg

···

Am 06.02.2018 um 10:50 schrieb ‘Johannes Boog’ via ogs-users:

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.).
      1. 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

        2. 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

  1. 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

  1. 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

  You received this message because you are subscribed to the Google Groups "ogs-users" group.

  To unsubscribe from this group and stop receiving emails from it, send an email to ogs-users+unsubscribe@googlegroups.com.

  For more options, visit [https://groups.google.com/d/optout](https://groups.google.com/d/optout).

Hello Georg,

material ID’s in the dual_van.msh are all the same. In case of DUAL_RICHARDS, you do not add material IDs in the mesh file for matrix and fracture this is just defined in the *.mmp file.

In my single RICHARDS_FLOW model, I assigned all three materials in the mesh file, of course. Now it is a question how to define matrix and fracture part for all three materials.

Best,

···

Am Dienstag, 6. Februar 2018, 11:24:02 MEZ hat Georg Kosakowski georg.kosakowski@psi.ch Folgendes geschrieben:

Hi Johannes,

  the materials have to be present in the mesh, i.e. you have to

assign the materials to the corresponding parts of the mesh
(material properties are assigned to the elements). Please have a
look how the mesh file looks in the H_us/dual_van benchmark.

CU, Georg

  Am 06.02.2018 um 10:50 schrieb 'Johannes Boog' via ogs-users:

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.).

      1. 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

        2. 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

  1. 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

  1. 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

  You received this message because you are subscribed to the Google

Groups “ogs-users” group.

  To unsubscribe from this group and stop receiving emails from it,

send an email to ogs-users+unsubscribe@googlegroups.com.

  For more options, visit [https://groups.google.com/d/optout](https://groups.google.com/d/optout).

You received this message because you are subscribed to the Google Groups “ogs-users” group.

To unsubscribe from this group and stop receiving emails from it, send an email to ogs-users+unsubscribe@googlegroups.com.

For more options, visit https://groups.google.com/d/optout.