Input file mcp keyword $ISOTHERM

Description

This keyword specifies, if the component sorbs in the phase in which it is transported simple retardation. Linear, Freundlich and Langmuir sorption are accounted for. No mass is transferred to the immobile phase. Aqueous and sorbed concentrations are termed C and CS, respectively.

Parameters

$ISOTHERM 0

$ISOTHERM
1 K_D

$ISOTHERM
2 K_fr e

$ISOTHERM
3 K L

Isotherm Method

• -1: No isotherm
• 0: get KD from curve derivative. Give curve_number. curve derivative is taken at the element_concentration.
• 1: Henry Isotherm (Linear isotherm). Give method linear_distribution_coefficient.

  \[ C_S = K_D C \]

• 2: Freundlich Isotherm. Give method Freundlich_coefficient freundlich_exponent.

  \[ C_S = K_{fr} C^e \]

• 3: Langmuir Isotherm. Give param1 param2.

  \[ C_S = \frac{KC}{1+KC} \]

• 4: Linear Isotherm, for fracture surface CMCD
• 5: two-rate model

Documentation from doxygen

None.

Used in benchmarks

Anisotropy/moleculardiffusion/soil_layer.mcp

 23:   $ISOTHERM
 24-    2  1e-3 0.9 ;
   ...

Anisotropy/permeability/soil_layer.mcp

 23:   $ISOTHERM
 24-    2  1e-3 0.9 ;
   ...

C/1d_analyt/1d_1.mcp

 12:  ; $ISOTHERM
 13-  ;  1  1.0e-03 ;
   ...
 22:   $ISOTHERM
 23-    2  1e-3 1.0 ;  0.9
   ...
 34:   $ISOTHERM
 35-    3  1e-3 0.5 ;
   ...

C/diffusion/Diff_HTO_test.mcp

 29:    $ISOTHERM
 30-    1 0.003  ;lineare Isotherme
   ...

C/diffusion/diff_aniso.mcp

 23:   $ISOTHERM
 24-    2  1e-3 0.9 ;
   ...

C/sorption/Freundlich/HC_sorp_Freundl_1D.mcp

 10:   $ISOTHERM
 11-    2  6.8e-6 0.99
   ...

C/sorption/Henry/HC_sorp_henry_1D.mcp

 10:   $ISOTHERM
 11-    1  6.8e-4
   ...

C/sorption/Langmuir/HC_sorp_langmuir_1D.mcp

 10:   $ISOTHERM
 11-    3 6.8e-5 1
   ...

C/sorption_decay/HC_decay_sorp_henry_1Du.mcp

 10:   $ISOTHERM
 11-    1  6.8e-5
   ...

FLUID_MOMENTUM/1d_hex/1d_hex.mcp

 14:   $ISOTHERM
 15-    1  1e-3 ;
   ...
 26:   $ISOTHERM
 27-    2  1e-3 0.9 ;
   ...
 38:   $ISOTHERM
 39-    3  1e-3 0.5 ;
   ...
 50:   $ISOTHERM
 51-    1 0.0 ;
   ...

FLUID_MOMENTUM/1d_line/1d_line.mcp

 14:   $ISOTHERM
 15-    1  1e-3 ;
   ...
 24:   $ISOTHERM
 25-    2  1e-3 0.9 ;
   ...
 36:   $ISOTHERM
 37-    3  1e-3 0.5 ;
   ...

FLUID_MOMENTUM/1d_pri/1d_pri.mcp

 14:   $ISOTHERM
 15-    1  1e-3 ;
   ...
 24:   $ISOTHERM
 25-    2  1e-3 0.9 ;
   ...
 36:   $ISOTHERM
 37-    3  1e-3 0.5 ;
   ...

FLUID_MOMENTUM/1d_quad/1d_quad.mcp

 14:   $ISOTHERM
 15-    1  1e-3 ;
   ...
 24:   $ISOTHERM
 25-    2  1e-3 0.9 ;
   ...
 36:   $ISOTHERM
 37-    3  1e-3 0.5 ;
   ...

FLUID_MOMENTUM/1d_tet/1d_tet.mcp

 14:   $ISOTHERM
 15-    1  1e-3 ;
   ...
 24:   $ISOTHERM
 25-    2  1e-3 0.9 ;
   ...
 36:   $ISOTHERM
 37-    3  1e-3 0.5 ;
   ...

FLUID_MOMENTUM/1d_tri/1d_tri.mcp

 15:   $ISOTHERM
 16-    1  1e-3 ;
   ...
 25:   $ISOTHERM
 26-    2  1e-3 0.9 ;
   ...
 37:   $ISOTHERM
 38-    3  1e-3 0.5 ;
   ...
 49:   $ISOTHERM
 50-    -1 ;
   ...

H/sat_2D/H_sat_flow_K_ortho.mcp

 23:   $ISOTHERM
 24-    2  1e-3 0.9 ;
   ...

H_us/Dual/dual_vl.mcp

 25:   $ISOTHERM
 26-    2  1e-3 0.9 ;
   ...

NUMERICS/DISCRETE_FEATURES/FractureNetwork/TH_3D_fracture_matrix.mcp

  8:   $ISOTHERM
  9-    1  4e-4
   ...

NUMERICS/DISCRETE_FEATURES/FractureNetwork/tracer_2.5D_kluefte.mcp

  8:   $ISOTHERM
  9-    1  4e-4
   ...

RWPT/2DGrains/2d_grains.mcp

 11:   $ISOTHERM
 12-    5  0.01  0.1  0.001
   ...

RWPT/3DGrain/3d_grain.mcp

 11:   $ISOTHERM
 12-    5  0.01  0.1  0.001
   ...

RWPT/Forchheimer/forchheimer_rwpt.mcp

 11:   $ISOTHERM
 12-    5  0.9  0.1  0.001
   ...

RWPT/Harter/colloid_t.mcp

 11:   $ISOTHERM
 12-    5  0.9  0.1  0.001
   ...

RWPT/HomoCube/3DRWPTCubTet.mcp

 11:   $ISOTHERM
 12-    5  0.01  0.1  0.001
   ...

RWPT/Veri1000/quad_homo.mcp

 11:   $ISOTHERM
 12-    5  0.01  0.1  0.001
   ...

THM/deco_2015_b2s1/deco_b2s1.mcp

 12:   $ISOTHERM
 13-    1  1e-3 ; sorption
   ...

Source code reading this keyword from Input files

FEM/rfmat_cp.cpp line 505

 505: 			in.clear();
 506: 		}                         // subkeyword found
 507: 		                          // subkeyword found
 508: 		if(line_string.find("$ISOTHERM") != std::string::npos)
 509: 		{
 510: 			in.str(GetLineFromFile1(rfd_file));
 511: 			in >> isotherm_model;

Last modified: 2016-01-21 13:48:59.637653

This page has been generated automatically. Do not edit it manually! It will be overwritten in regular intervals.

If you want to edit the handwritten part of this page, please do so here.