Materials and Media#

Section author: Norbert Grunwald

NB! Works for HEAT_CONDUCTION and TH2M (with Phase Transitions) only so far.

This example shows how to build OpenGeoSys media definitions directly from a YAML-based Material Library using MaterialManager and MediaSet.

The YAML schema defines materials with:

  • name - unique material identifier

  • properties - sets of constitutive relations (with type, parameters, scope)

  • scope - indicates whether a property applies at phase or medium level

Together with the built-in process schemas (e.g. TH2M_PT), these building blocks allow you to construct full Media definitions including phases and components, and import them into an OGS project XML via set_media().

from pathlib import Path
from tempfile import mkdtemp

import ogstools as ot
from ogstools import definitions as defs

model_dir = Path(mkdtemp())

Example materials#

Materials are provided as YAML files in the built-in Material Library. For example, here are the definitions of “opalinus_clay” (solid) and “water” (fluid):

print((Path(defs.MATERIALS_DIR) / "opalinus.yml").read_text(encoding="utf-8"))
print((Path(defs.MATERIALS_DIR) / "water.yml").read_text(encoding="utf-8"))
---
group: medium
name: opalinus_clay
properties:
  biot_coefficient:
    - type: Constant
      unit: 0.6
      value: 1
  bishops_effective_stress:
    - type: BishopsPowerLaw
      unit: 1
      exponent: 1
  density:
    - type: Constant
      unit: kg/m³
      value: 2600
  permeability:
    - type: Constant
      unit: m²
      value: 1.e-21
  poisson_ratio:
    - type: Constant
      unit: 1
      value: 0.23
  porosity:
    - type: Constant
      unit: 1
      value: 0.15
  saturation:
    - type: SaturationVanGenuchten
      exponent: 0.2
      p_b: 48000000
      residual_gas_saturation: 0.01
      residual_liquid_saturation: 0.01
  specific_heat_capacity:
    - type: Constant
      unit: J/(kg·K)
      value: 995
  thermal_conductivity:
    # Thermal conductivity of the individual phase (e.g. solid or fluid),
    # used in local energy balance equations of that phase.
    - type: Constant
      scope: phase
      unit: W/(m·K)
      value: 1.7
    # Effective thermal conductivity of the full medium,
    # typically calculated from the phase conductivities using
    # porosity-weighted mixing models.
    - type: SaturationWeightedThermalConductivity
      mean_type: arithmetic_squareroot
      dry_thermal_conductivity: 1.36
      wet_thermal_conductivity: 1.7
      scope: medium
      unit: W/(m·K)
    - type: EffectiveThermalConductivityPorosityMixing
      scope: medium
      unit: W/(m·K)
  thermal_expansivity:
    - type: Constant
      unit: 1/K
      value: 1.7e-5
  tortuosity:
    - type: Constant
      unit: 1
      value: 1
  youngs_modulus:
    - type: Constant
      unit: Pa
      value: 6500000000
  relative_permeability:
    - type: RelativePermeabilityVanGenuchten
      residual_liquid_saturation: 0.01
      residual_gas_saturation: 0.01
      exponent: 0.2
      minimum_relative_permeability_liquid: 1e-12
      unit: 1
  relative_permeability_nonwetting_phase:
    - type: RelativePermeabilityNonWettingPhaseVanGenuchtenMualem
      residual_liquid_saturation: 0.01
      residual_gas_saturation: 0.01
      exponent: 0.2
      min_relative_permeability: 1e-5
      unit: 1

name: water
group: liquid
type: phase
properties:
  molar_mass:
    - type: Constant
      value: 0.018016
      unit: kg/mol
      source: "CODATA 2018 recommended values"
  specific_heat_capacity:
    - type: Constant
      value: 4187
      unit: J/(kg·K)
      source: "CRC Handbook of Chemistry and Physics"
  thermal_conductivity:
    - type: Constant
      value: 0.6
      unit: W/(m·K)
      source: "CRC Handbook of Chemistry and Physics"
  density:
    - type: Constant
      value: 1000
      unit: kg/m³
      source: "at 4 °C, standard value"
    - type: Linear
      reference_value: 1000
      independent_variables:
        concentration:
          reference_condition: 0.0
          slope: 2.0e-6
        temperature:
          reference_condition: 293.15
          slope: 1.5e-5
      source: "CRC Handbook of Chemistry and Physics"
  viscosity:
    - type: Constant
      value: 1.0e-3
      unit: Pa·s
      source: "at 20 °C, standard value"
  vapour_pressure:
    - type: ClausiusClapeyron
      triple_temperature: 273.16
      triple_pressure: 611.66
      critical_temperature: 647.096
      critical_pressure: 22.064e6
      reference_temperature: 373.15
      reference_pressure: 101325
      source: "IAPWS, 1995 formulation"
  diffusion:
    - type: Constant
      value: 1.0e-9
      unit: m²/s
      source: "typical value for water in porous media"
  specific_latent_heat:
    - type: Constant
      value: 2268000.0
      unit: J/kg
      source: "at 100 °C, standard value"

Media creation#

We build a MaterialManager from the built-in library, filter it with a schema, and construct a MediaSet object.

  • subdomain: the subdomain name (string, one per entry)

  • material: must match the name in the YAML file

  • material_ids: list of integers corresponding to the MatIDs in the mesh (allows grouping several mesh regions under one subdomain name)

db = ot.MaterialManager()

subdomains = [
    {
        "subdomain": "host_rock",
        "material": "opalinus_clay",
        "material_ids": [0, 3, 4],  # multiple MatIDs grouped under one name
    },
    {
        "subdomain": "buffer",
        "material": "bentonite",
        "material_ids": [1, 2],
    },
]

fluids = {
    "AqueousLiquid": "water",
    "Gas": "carbon_dioxide",
}  # required by TH2M_PT schema

filtered = db.filter(process="TH2M_PT", subdomains=subdomains, fluids=fluids)

media = ot.MediaSet(filtered)

Export to OGS Project XML#

The MediaSet is imported into an OGS Project instance via Project.set_media().

prj = ot.Project()
prj.set_media(media)

xml_file = model_dir / "material_test.prj"
prj.write_input(xml_file)
print(xml_file.read_text())
<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
    <mesh/>
    <processes>
        <process/>
    </processes>
    <time_loop>
        <processes/>
    </time_loop>
    <parameters/>
    <process_variables/>
    <nonlinear_solvers/>
    <linear_solvers/>
    <media>
        <medium id="0">
            <properties>
                <property>
                    <name>biot_coefficient</name>
                    <type>Constant</type>
                    <value>1</value>
                </property>
                <property>
                    <name>bishops_effective_stress</name>
                    <type>BishopsPowerLaw</type>
                    <exponent>1</exponent>
                </property>
                <property>
                    <name>permeability</name>
                    <type>Constant</type>
                    <value>1e-21</value>
                </property>
                <property>
                    <name>porosity</name>
                    <type>Constant</type>
                    <value>0.15</value>
                </property>
                <property>
                    <name>saturation</name>
                    <type>SaturationVanGenuchten</type>
                    <exponent>0.2</exponent>
                    <p_b>48000000</p_b>
                    <residual_gas_saturation>0.01</residual_gas_saturation>
                    <residual_liquid_saturation>0.01</residual_liquid_saturation>
                </property>
                <property>
                    <name>thermal_conductivity</name>
                    <type>SaturationWeightedThermalConductivity</type>
                    <mean_type>arithmetic_squareroot</mean_type>
                    <dry_thermal_conductivity>1.36</dry_thermal_conductivity>
                    <wet_thermal_conductivity>1.7</wet_thermal_conductivity>
                </property>
                <property>
                    <name>thermal_conductivity</name>
                    <type>EffectiveThermalConductivityPorosityMixing</type>
                </property>
                <property>
                    <name>tortuosity</name>
                    <type>Constant</type>
                    <value>1</value>
                </property>
                <property>
                    <name>relative_permeability</name>
                    <type>RelativePermeabilityVanGenuchten</type>
                    <exponent>0.2</exponent>
                    <minimum_relative_permeability_liquid>1e-12</minimum_relative_permeability_liquid>
                    <residual_gas_saturation>0.01</residual_gas_saturation>
                    <residual_liquid_saturation>0.01</residual_liquid_saturation>
                </property>
                <property>
                    <name>relative_permeability_nonwetting_phase</name>
                    <type>RelativePermeabilityNonWettingPhaseVanGenuchtenMualem</type>
                    <exponent>0.2</exponent>
                    <min_relative_permeability>1e-5</min_relative_permeability>
                    <residual_gas_saturation>0.01</residual_gas_saturation>
                    <residual_liquid_saturation>0.01</residual_liquid_saturation>
                </property>
            </properties>
            <phases>
                <phase>
                    <type>Solid</type>
                    <properties>
                        <property>
                            <name>density</name>
                            <type>Constant</type>
                            <value>2600</value>
                        </property>
                        <property>
                            <name>specific_heat_capacity</name>
                            <type>Constant</type>
                            <value>995</value>
                        </property>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>1.7</value>
                        </property>
                        <property>
                            <name>thermal_expansivity</name>
                            <type>Constant</type>
                            <value>1.7e-05</value>
                        </property>
                    </properties>
                </phase>
                <phase>
                    <type>AqueousLiquid</type>
                    <properties>
                        <property>
                            <name>specific_heat_capacity</name>
                            <type>Constant</type>
                            <value>4187</value>
                        </property>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>0.6</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>Constant</type>
                            <value>1000</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>Linear</type>
                            <reference_value>1000</reference_value>
                            <independent_variable>
                                <variable_name>concentration</variable_name>
                                <reference_condition>0.0</reference_condition>
                                <slope>2e-06</slope>
                            </independent_variable>
                            <independent_variable>
                                <variable_name>temperature</variable_name>
                                <reference_condition>293.15</reference_condition>
                                <slope>1.5e-05</slope>
                            </independent_variable>
                        </property>
                        <property>
                            <name>viscosity</name>
                            <type>Constant</type>
                            <value>0.001</value>
                        </property>
                    </properties>
                    <components>
                        <component>
                            <name>carbon_dioxide</name>
                            <properties>
                                <property>
                                    <name>specific_latent_heat</name>
                                    <type>Constant</type>
                                    <value>100.0</value>
                                </property>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                                <property>
                                    <name>diffusion</name>
                                    <type>Constant</type>
                                    <value>0.0</value>
                                </property>
                                <property>
                                    <name>henry_coefficient</name>
                                    <type>Constant</type>
                                    <value>1.0</value>
                                </property>
                            </properties>
                        </component>
                        <component>
                            <name>water</name>
                            <properties>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                            </properties>
                        </component>
                    </components>
                </phase>
                <phase>
                    <type>Gas</type>
                    <properties>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>0.6</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>IdealGasLawBinaryMixture</type>
                        </property>
                        <property>
                            <name>viscosity</name>
                            <type>Constant</type>
                            <value>0.001</value>
                        </property>
                    </properties>
                    <components>
                        <component>
                            <name>carbon_dioxide</name>
                            <properties>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                                <property>
                                    <name>molar_mass</name>
                                    <type>Constant</type>
                                    <value>0.018016</value>
                                </property>
                            </properties>
                        </component>
                        <component>
                            <name>water</name>
                            <properties>
                                <property>
                                    <name>vapour_pressure</name>
                                    <type>ClausiusClapeyron</type>
                                    <critical_pressure>22.064e6</critical_pressure>
                                    <critical_temperature>647.096</critical_temperature>
                                    <reference_pressure>101325</reference_pressure>
                                    <reference_temperature>373.15</reference_temperature>
                                    <triple_pressure>611.66</triple_pressure>
                                    <triple_temperature>273.16</triple_temperature>
                                </property>
                                <property>
                                    <name>diffusion</name>
                                    <type>Constant</type>
                                    <value>1e-06</value>
                                </property>
                                <property>
                                    <name>molar_mass</name>
                                    <type>Constant</type>
                                    <value>0.018016</value>
                                </property>
                                <property>
                                    <name>specific_latent_heat</name>
                                    <type>Constant</type>
                                    <value>2268000.0</value>
                                </property>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                            </properties>
                        </component>
                    </components>
                </phase>
            </phases>
        </medium>
        <medium id="3">
            <properties>
                <property>
                    <name>biot_coefficient</name>
                    <type>Constant</type>
                    <value>1</value>
                </property>
                <property>
                    <name>bishops_effective_stress</name>
                    <type>BishopsPowerLaw</type>
                    <exponent>1</exponent>
                </property>
                <property>
                    <name>permeability</name>
                    <type>Constant</type>
                    <value>1e-21</value>
                </property>
                <property>
                    <name>porosity</name>
                    <type>Constant</type>
                    <value>0.15</value>
                </property>
                <property>
                    <name>saturation</name>
                    <type>SaturationVanGenuchten</type>
                    <exponent>0.2</exponent>
                    <p_b>48000000</p_b>
                    <residual_gas_saturation>0.01</residual_gas_saturation>
                    <residual_liquid_saturation>0.01</residual_liquid_saturation>
                </property>
                <property>
                    <name>thermal_conductivity</name>
                    <type>SaturationWeightedThermalConductivity</type>
                    <mean_type>arithmetic_squareroot</mean_type>
                    <dry_thermal_conductivity>1.36</dry_thermal_conductivity>
                    <wet_thermal_conductivity>1.7</wet_thermal_conductivity>
                </property>
                <property>
                    <name>thermal_conductivity</name>
                    <type>EffectiveThermalConductivityPorosityMixing</type>
                </property>
                <property>
                    <name>tortuosity</name>
                    <type>Constant</type>
                    <value>1</value>
                </property>
                <property>
                    <name>relative_permeability</name>
                    <type>RelativePermeabilityVanGenuchten</type>
                    <exponent>0.2</exponent>
                    <minimum_relative_permeability_liquid>1e-12</minimum_relative_permeability_liquid>
                    <residual_gas_saturation>0.01</residual_gas_saturation>
                    <residual_liquid_saturation>0.01</residual_liquid_saturation>
                </property>
                <property>
                    <name>relative_permeability_nonwetting_phase</name>
                    <type>RelativePermeabilityNonWettingPhaseVanGenuchtenMualem</type>
                    <exponent>0.2</exponent>
                    <min_relative_permeability>1e-5</min_relative_permeability>
                    <residual_gas_saturation>0.01</residual_gas_saturation>
                    <residual_liquid_saturation>0.01</residual_liquid_saturation>
                </property>
            </properties>
            <phases>
                <phase>
                    <type>Solid</type>
                    <properties>
                        <property>
                            <name>density</name>
                            <type>Constant</type>
                            <value>2600</value>
                        </property>
                        <property>
                            <name>specific_heat_capacity</name>
                            <type>Constant</type>
                            <value>995</value>
                        </property>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>1.7</value>
                        </property>
                        <property>
                            <name>thermal_expansivity</name>
                            <type>Constant</type>
                            <value>1.7e-05</value>
                        </property>
                    </properties>
                </phase>
                <phase>
                    <type>AqueousLiquid</type>
                    <properties>
                        <property>
                            <name>specific_heat_capacity</name>
                            <type>Constant</type>
                            <value>4187</value>
                        </property>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>0.6</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>Constant</type>
                            <value>1000</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>Linear</type>
                            <reference_value>1000</reference_value>
                            <independent_variable>
                                <variable_name>concentration</variable_name>
                                <reference_condition>0.0</reference_condition>
                                <slope>2e-06</slope>
                            </independent_variable>
                            <independent_variable>
                                <variable_name>temperature</variable_name>
                                <reference_condition>293.15</reference_condition>
                                <slope>1.5e-05</slope>
                            </independent_variable>
                        </property>
                        <property>
                            <name>viscosity</name>
                            <type>Constant</type>
                            <value>0.001</value>
                        </property>
                    </properties>
                    <components>
                        <component>
                            <name>carbon_dioxide</name>
                            <properties>
                                <property>
                                    <name>specific_latent_heat</name>
                                    <type>Constant</type>
                                    <value>100.0</value>
                                </property>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                                <property>
                                    <name>diffusion</name>
                                    <type>Constant</type>
                                    <value>0.0</value>
                                </property>
                                <property>
                                    <name>henry_coefficient</name>
                                    <type>Constant</type>
                                    <value>1.0</value>
                                </property>
                            </properties>
                        </component>
                        <component>
                            <name>water</name>
                            <properties>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                            </properties>
                        </component>
                    </components>
                </phase>
                <phase>
                    <type>Gas</type>
                    <properties>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>0.6</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>IdealGasLawBinaryMixture</type>
                        </property>
                        <property>
                            <name>viscosity</name>
                            <type>Constant</type>
                            <value>0.001</value>
                        </property>
                    </properties>
                    <components>
                        <component>
                            <name>carbon_dioxide</name>
                            <properties>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                                <property>
                                    <name>molar_mass</name>
                                    <type>Constant</type>
                                    <value>0.018016</value>
                                </property>
                            </properties>
                        </component>
                        <component>
                            <name>water</name>
                            <properties>
                                <property>
                                    <name>vapour_pressure</name>
                                    <type>ClausiusClapeyron</type>
                                    <critical_pressure>22.064e6</critical_pressure>
                                    <critical_temperature>647.096</critical_temperature>
                                    <reference_pressure>101325</reference_pressure>
                                    <reference_temperature>373.15</reference_temperature>
                                    <triple_pressure>611.66</triple_pressure>
                                    <triple_temperature>273.16</triple_temperature>
                                </property>
                                <property>
                                    <name>diffusion</name>
                                    <type>Constant</type>
                                    <value>1e-06</value>
                                </property>
                                <property>
                                    <name>molar_mass</name>
                                    <type>Constant</type>
                                    <value>0.018016</value>
                                </property>
                                <property>
                                    <name>specific_latent_heat</name>
                                    <type>Constant</type>
                                    <value>2268000.0</value>
                                </property>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                            </properties>
                        </component>
                    </components>
                </phase>
            </phases>
        </medium>
        <medium id="4">
            <properties>
                <property>
                    <name>biot_coefficient</name>
                    <type>Constant</type>
                    <value>1</value>
                </property>
                <property>
                    <name>bishops_effective_stress</name>
                    <type>BishopsPowerLaw</type>
                    <exponent>1</exponent>
                </property>
                <property>
                    <name>permeability</name>
                    <type>Constant</type>
                    <value>1e-21</value>
                </property>
                <property>
                    <name>porosity</name>
                    <type>Constant</type>
                    <value>0.15</value>
                </property>
                <property>
                    <name>saturation</name>
                    <type>SaturationVanGenuchten</type>
                    <exponent>0.2</exponent>
                    <p_b>48000000</p_b>
                    <residual_gas_saturation>0.01</residual_gas_saturation>
                    <residual_liquid_saturation>0.01</residual_liquid_saturation>
                </property>
                <property>
                    <name>thermal_conductivity</name>
                    <type>SaturationWeightedThermalConductivity</type>
                    <mean_type>arithmetic_squareroot</mean_type>
                    <dry_thermal_conductivity>1.36</dry_thermal_conductivity>
                    <wet_thermal_conductivity>1.7</wet_thermal_conductivity>
                </property>
                <property>
                    <name>thermal_conductivity</name>
                    <type>EffectiveThermalConductivityPorosityMixing</type>
                </property>
                <property>
                    <name>tortuosity</name>
                    <type>Constant</type>
                    <value>1</value>
                </property>
                <property>
                    <name>relative_permeability</name>
                    <type>RelativePermeabilityVanGenuchten</type>
                    <exponent>0.2</exponent>
                    <minimum_relative_permeability_liquid>1e-12</minimum_relative_permeability_liquid>
                    <residual_gas_saturation>0.01</residual_gas_saturation>
                    <residual_liquid_saturation>0.01</residual_liquid_saturation>
                </property>
                <property>
                    <name>relative_permeability_nonwetting_phase</name>
                    <type>RelativePermeabilityNonWettingPhaseVanGenuchtenMualem</type>
                    <exponent>0.2</exponent>
                    <min_relative_permeability>1e-5</min_relative_permeability>
                    <residual_gas_saturation>0.01</residual_gas_saturation>
                    <residual_liquid_saturation>0.01</residual_liquid_saturation>
                </property>
            </properties>
            <phases>
                <phase>
                    <type>Solid</type>
                    <properties>
                        <property>
                            <name>density</name>
                            <type>Constant</type>
                            <value>2600</value>
                        </property>
                        <property>
                            <name>specific_heat_capacity</name>
                            <type>Constant</type>
                            <value>995</value>
                        </property>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>1.7</value>
                        </property>
                        <property>
                            <name>thermal_expansivity</name>
                            <type>Constant</type>
                            <value>1.7e-05</value>
                        </property>
                    </properties>
                </phase>
                <phase>
                    <type>AqueousLiquid</type>
                    <properties>
                        <property>
                            <name>specific_heat_capacity</name>
                            <type>Constant</type>
                            <value>4187</value>
                        </property>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>0.6</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>Constant</type>
                            <value>1000</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>Linear</type>
                            <reference_value>1000</reference_value>
                            <independent_variable>
                                <variable_name>concentration</variable_name>
                                <reference_condition>0.0</reference_condition>
                                <slope>2e-06</slope>
                            </independent_variable>
                            <independent_variable>
                                <variable_name>temperature</variable_name>
                                <reference_condition>293.15</reference_condition>
                                <slope>1.5e-05</slope>
                            </independent_variable>
                        </property>
                        <property>
                            <name>viscosity</name>
                            <type>Constant</type>
                            <value>0.001</value>
                        </property>
                    </properties>
                    <components>
                        <component>
                            <name>carbon_dioxide</name>
                            <properties>
                                <property>
                                    <name>specific_latent_heat</name>
                                    <type>Constant</type>
                                    <value>100.0</value>
                                </property>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                                <property>
                                    <name>diffusion</name>
                                    <type>Constant</type>
                                    <value>0.0</value>
                                </property>
                                <property>
                                    <name>henry_coefficient</name>
                                    <type>Constant</type>
                                    <value>1.0</value>
                                </property>
                            </properties>
                        </component>
                        <component>
                            <name>water</name>
                            <properties>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                            </properties>
                        </component>
                    </components>
                </phase>
                <phase>
                    <type>Gas</type>
                    <properties>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>0.6</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>IdealGasLawBinaryMixture</type>
                        </property>
                        <property>
                            <name>viscosity</name>
                            <type>Constant</type>
                            <value>0.001</value>
                        </property>
                    </properties>
                    <components>
                        <component>
                            <name>carbon_dioxide</name>
                            <properties>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                                <property>
                                    <name>molar_mass</name>
                                    <type>Constant</type>
                                    <value>0.018016</value>
                                </property>
                            </properties>
                        </component>
                        <component>
                            <name>water</name>
                            <properties>
                                <property>
                                    <name>vapour_pressure</name>
                                    <type>ClausiusClapeyron</type>
                                    <critical_pressure>22.064e6</critical_pressure>
                                    <critical_temperature>647.096</critical_temperature>
                                    <reference_pressure>101325</reference_pressure>
                                    <reference_temperature>373.15</reference_temperature>
                                    <triple_pressure>611.66</triple_pressure>
                                    <triple_temperature>273.16</triple_temperature>
                                </property>
                                <property>
                                    <name>diffusion</name>
                                    <type>Constant</type>
                                    <value>1e-06</value>
                                </property>
                                <property>
                                    <name>molar_mass</name>
                                    <type>Constant</type>
                                    <value>0.018016</value>
                                </property>
                                <property>
                                    <name>specific_latent_heat</name>
                                    <type>Constant</type>
                                    <value>2268000.0</value>
                                </property>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                            </properties>
                        </component>
                    </components>
                </phase>
            </phases>
        </medium>
        <medium id="1">
            <properties>
                <property>
                    <name>thermal_conductivity</name>
                    <type>Constant</type>
                    <value>0.6</value>
                </property>
                <property>
                    <name>porosity</name>
                    <type>Constant</type>
                    <value>0.4</value>
                </property>
                <property>
                    <name>permeability</name>
                    <type>Constant</type>
                    <value>1e-20</value>
                </property>
                <property>
                    <name>saturation</name>
                    <type>SaturationVanGenuchten</type>
                    <exponent>0.6</exponent>
                    <p_b>2.0e6</p_b>
                    <residual_gas_saturation>1e-15</residual_gas_saturation>
                    <residual_liquid_saturation>0.05</residual_liquid_saturation>
                </property>
                <property>
                    <name>bishops_effective_stress</name>
                    <type>BishopsPowerLaw</type>
                    <exponent>1</exponent>
                </property>
                <property>
                    <name>biot_coefficient</name>
                    <type>Constant</type>
                    <value>1.0</value>
                </property>
                <property>
                    <name>tortuosity</name>
                    <type>Constant</type>
                    <value>1.0</value>
                </property>
                <property>
                    <name>relative_permeability_nonwetting_phase</name>
                    <type>RelativePermeabilityNonWettingPhaseVanGenuchtenMualem</type>
                    <exponent>0.3</exponent>
                    <min_relative_permeability>1e-5</min_relative_permeability>
                    <residual_gas_saturation>1e-15</residual_gas_saturation>
                    <residual_liquid_saturation>0.05</residual_liquid_saturation>
                </property>
                <property>
                    <name>relative_permeability</name>
                    <type>RelativePermeabilityVanGenuchten</type>
                    <exponent>0.3</exponent>
                    <minimum_relative_permeability_liquid>1e-12</minimum_relative_permeability_liquid>
                    <residual_gas_saturation>1e-15</residual_gas_saturation>
                    <residual_liquid_saturation>0.05</residual_liquid_saturation>
                </property>
            </properties>
            <phases>
                <phase>
                    <type>Solid</type>
                    <properties>
                        <property>
                            <name>density</name>
                            <type>Constant</type>
                            <value>2500</value>
                        </property>
                        <property>
                            <name>specific_heat_capacity</name>
                            <type>Constant</type>
                            <value>900</value>
                        </property>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>0.6</value>
                        </property>
                        <property>
                            <name>thermal_expansivity</name>
                            <type>Constant</type>
                            <value>2e-05</value>
                        </property>
                    </properties>
                </phase>
                <phase>
                    <type>AqueousLiquid</type>
                    <properties>
                        <property>
                            <name>specific_heat_capacity</name>
                            <type>Constant</type>
                            <value>4187</value>
                        </property>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>0.6</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>Constant</type>
                            <value>1000</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>Linear</type>
                            <reference_value>1000</reference_value>
                            <independent_variable>
                                <variable_name>concentration</variable_name>
                                <reference_condition>0.0</reference_condition>
                                <slope>2e-06</slope>
                            </independent_variable>
                            <independent_variable>
                                <variable_name>temperature</variable_name>
                                <reference_condition>293.15</reference_condition>
                                <slope>1.5e-05</slope>
                            </independent_variable>
                        </property>
                        <property>
                            <name>viscosity</name>
                            <type>Constant</type>
                            <value>0.001</value>
                        </property>
                    </properties>
                    <components>
                        <component>
                            <name>carbon_dioxide</name>
                            <properties>
                                <property>
                                    <name>specific_latent_heat</name>
                                    <type>Constant</type>
                                    <value>100.0</value>
                                </property>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                                <property>
                                    <name>diffusion</name>
                                    <type>Constant</type>
                                    <value>0.0</value>
                                </property>
                                <property>
                                    <name>henry_coefficient</name>
                                    <type>Constant</type>
                                    <value>1.0</value>
                                </property>
                            </properties>
                        </component>
                        <component>
                            <name>water</name>
                            <properties>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                            </properties>
                        </component>
                    </components>
                </phase>
                <phase>
                    <type>Gas</type>
                    <properties>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>0.6</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>IdealGasLawBinaryMixture</type>
                        </property>
                        <property>
                            <name>viscosity</name>
                            <type>Constant</type>
                            <value>0.001</value>
                        </property>
                    </properties>
                    <components>
                        <component>
                            <name>carbon_dioxide</name>
                            <properties>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                                <property>
                                    <name>molar_mass</name>
                                    <type>Constant</type>
                                    <value>0.018016</value>
                                </property>
                            </properties>
                        </component>
                        <component>
                            <name>water</name>
                            <properties>
                                <property>
                                    <name>vapour_pressure</name>
                                    <type>ClausiusClapeyron</type>
                                    <critical_pressure>22.064e6</critical_pressure>
                                    <critical_temperature>647.096</critical_temperature>
                                    <reference_pressure>101325</reference_pressure>
                                    <reference_temperature>373.15</reference_temperature>
                                    <triple_pressure>611.66</triple_pressure>
                                    <triple_temperature>273.16</triple_temperature>
                                </property>
                                <property>
                                    <name>diffusion</name>
                                    <type>Constant</type>
                                    <value>1e-06</value>
                                </property>
                                <property>
                                    <name>molar_mass</name>
                                    <type>Constant</type>
                                    <value>0.018016</value>
                                </property>
                                <property>
                                    <name>specific_latent_heat</name>
                                    <type>Constant</type>
                                    <value>2268000.0</value>
                                </property>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                            </properties>
                        </component>
                    </components>
                </phase>
            </phases>
        </medium>
        <medium id="2">
            <properties>
                <property>
                    <name>thermal_conductivity</name>
                    <type>Constant</type>
                    <value>0.6</value>
                </property>
                <property>
                    <name>porosity</name>
                    <type>Constant</type>
                    <value>0.4</value>
                </property>
                <property>
                    <name>permeability</name>
                    <type>Constant</type>
                    <value>1e-20</value>
                </property>
                <property>
                    <name>saturation</name>
                    <type>SaturationVanGenuchten</type>
                    <exponent>0.6</exponent>
                    <p_b>2.0e6</p_b>
                    <residual_gas_saturation>1e-15</residual_gas_saturation>
                    <residual_liquid_saturation>0.05</residual_liquid_saturation>
                </property>
                <property>
                    <name>bishops_effective_stress</name>
                    <type>BishopsPowerLaw</type>
                    <exponent>1</exponent>
                </property>
                <property>
                    <name>biot_coefficient</name>
                    <type>Constant</type>
                    <value>1.0</value>
                </property>
                <property>
                    <name>tortuosity</name>
                    <type>Constant</type>
                    <value>1.0</value>
                </property>
                <property>
                    <name>relative_permeability_nonwetting_phase</name>
                    <type>RelativePermeabilityNonWettingPhaseVanGenuchtenMualem</type>
                    <exponent>0.3</exponent>
                    <min_relative_permeability>1e-5</min_relative_permeability>
                    <residual_gas_saturation>1e-15</residual_gas_saturation>
                    <residual_liquid_saturation>0.05</residual_liquid_saturation>
                </property>
                <property>
                    <name>relative_permeability</name>
                    <type>RelativePermeabilityVanGenuchten</type>
                    <exponent>0.3</exponent>
                    <minimum_relative_permeability_liquid>1e-12</minimum_relative_permeability_liquid>
                    <residual_gas_saturation>1e-15</residual_gas_saturation>
                    <residual_liquid_saturation>0.05</residual_liquid_saturation>
                </property>
            </properties>
            <phases>
                <phase>
                    <type>Solid</type>
                    <properties>
                        <property>
                            <name>density</name>
                            <type>Constant</type>
                            <value>2500</value>
                        </property>
                        <property>
                            <name>specific_heat_capacity</name>
                            <type>Constant</type>
                            <value>900</value>
                        </property>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>0.6</value>
                        </property>
                        <property>
                            <name>thermal_expansivity</name>
                            <type>Constant</type>
                            <value>2e-05</value>
                        </property>
                    </properties>
                </phase>
                <phase>
                    <type>AqueousLiquid</type>
                    <properties>
                        <property>
                            <name>specific_heat_capacity</name>
                            <type>Constant</type>
                            <value>4187</value>
                        </property>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>0.6</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>Constant</type>
                            <value>1000</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>Linear</type>
                            <reference_value>1000</reference_value>
                            <independent_variable>
                                <variable_name>concentration</variable_name>
                                <reference_condition>0.0</reference_condition>
                                <slope>2e-06</slope>
                            </independent_variable>
                            <independent_variable>
                                <variable_name>temperature</variable_name>
                                <reference_condition>293.15</reference_condition>
                                <slope>1.5e-05</slope>
                            </independent_variable>
                        </property>
                        <property>
                            <name>viscosity</name>
                            <type>Constant</type>
                            <value>0.001</value>
                        </property>
                    </properties>
                    <components>
                        <component>
                            <name>carbon_dioxide</name>
                            <properties>
                                <property>
                                    <name>specific_latent_heat</name>
                                    <type>Constant</type>
                                    <value>100.0</value>
                                </property>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                                <property>
                                    <name>diffusion</name>
                                    <type>Constant</type>
                                    <value>0.0</value>
                                </property>
                                <property>
                                    <name>henry_coefficient</name>
                                    <type>Constant</type>
                                    <value>1.0</value>
                                </property>
                            </properties>
                        </component>
                        <component>
                            <name>water</name>
                            <properties>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                            </properties>
                        </component>
                    </components>
                </phase>
                <phase>
                    <type>Gas</type>
                    <properties>
                        <property>
                            <name>thermal_conductivity</name>
                            <type>Constant</type>
                            <value>0.6</value>
                        </property>
                        <property>
                            <name>density</name>
                            <type>IdealGasLawBinaryMixture</type>
                        </property>
                        <property>
                            <name>viscosity</name>
                            <type>Constant</type>
                            <value>0.001</value>
                        </property>
                    </properties>
                    <components>
                        <component>
                            <name>carbon_dioxide</name>
                            <properties>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                                <property>
                                    <name>molar_mass</name>
                                    <type>Constant</type>
                                    <value>0.018016</value>
                                </property>
                            </properties>
                        </component>
                        <component>
                            <name>water</name>
                            <properties>
                                <property>
                                    <name>vapour_pressure</name>
                                    <type>ClausiusClapeyron</type>
                                    <critical_pressure>22.064e6</critical_pressure>
                                    <critical_temperature>647.096</critical_temperature>
                                    <reference_pressure>101325</reference_pressure>
                                    <reference_temperature>373.15</reference_temperature>
                                    <triple_pressure>611.66</triple_pressure>
                                    <triple_temperature>273.16</triple_temperature>
                                </property>
                                <property>
                                    <name>diffusion</name>
                                    <type>Constant</type>
                                    <value>1e-06</value>
                                </property>
                                <property>
                                    <name>molar_mass</name>
                                    <type>Constant</type>
                                    <value>0.018016</value>
                                </property>
                                <property>
                                    <name>specific_latent_heat</name>
                                    <type>Constant</type>
                                    <value>2268000.0</value>
                                </property>
                                <property>
                                    <name>specific_heat_capacity</name>
                                    <type>Constant</type>
                                    <value>4187</value>
                                </property>
                            </properties>
                        </component>
                    </components>
                </phase>
            </phases>
        </medium>
    </media>
</OpenGeoSysProject>

Total running time of the script: (0 minutes 0.111 seconds)