Back to Input Deck Cards
Back to CHEMISTRY
Back to SORPTION
ISOTHERM_REACTIONS¶
Specifies parameters for a sorption reaction defined by an isotherm (e.g. linear, Langmuir, Freundlich).
Required Cards:¶
- ISOTHERM_REACTIONS
Opens the block.
- <Species_Name>
Name of primary species that sorbs.
- DISTRIBUTION_COEFFICIENT <float>
The value of KD [kgwater / m3 bulk]. For KD values provided in e.g. mL / g, convert to kgw / m3 bulk by multiplying by rhowater and rhobulk , with proper units conversion. If KD_MINERAL_NAME is included, the units of KD become [mLwater / grock].
Optional Cards:¶
- TYPE <string>
Type of isotherm, where the options for <string> include:
LINEAR : Linear isotherm (KD) – Csorb = KD Caq
LANGMUIR : Langmuir isotherm – Csorb = K Caq b / ( 1 + K Caq)
FREUNDLICH : Freundlich isotherm – Csorb = KD (Caq) 1/n
- LANGMUIR_B <float>
b coefficient for Langmuir isotherm. Automatically sets the type to LANGMUIR.
- FREUNDLICH_N <float>
n exponent in Freundlich isotherm. Automatically sets the type to FREUNDLICH.
- KD_MINERAL_NAME <string>
Name of mineral, the volume fraction of which will be used in conjunction with the distribution coefficient to calculate the PFLOTRAN internal distribution coefficient of kg water per cubic meter bulk. The distribution coefficient then has units of Lwater / kgsoil. The conversion equation is KD (kgwater / m3bulk) = KD (mLwater / gsoil) * water density * (1 - porosity) * soil particle density * 10-3 * mineral volume fraction. Note that the mineral volume fraction has no units in this case; it is solely a scaling parameter for varying the KD in space.
Examples¶
ISOTHERM_REACTIONS Tracer DISTRIBUTION_COEFFICIENT 0.25d3 / Tracer2 DISTRIBUTION_COEFFICIENT 0.5d3 FREUNDLICH_N 1.67 ! 1/n = 0.6 / END ISOTHERM_REACTIONS Tracer ! Kd units = mL water/kg soil ! assuming: ! water_density = 1000 kg/m^3 ! soil particle density = 2500 kg/m^3 ! porosity = 0.25 DISTRIBUTION_COEFFICIENT 0.133333 KD_MINERAL_NAME Calcite / END