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.