example Langmuir isotherm

[Charlie]

# Here is an example for an Langmuir isotherm using solution spread block.  Results can be compared for simulations using different parameters such as log_k..

TITLE Langmuir isotherm (inspired by Table 7.6 in Zhu and Anderson, 2002)
 
SELECTED_OUTPUT
      -file Langmuir.txt
      -reset false
 
USER_PUNCH
      -headings Total_Ca   Ca+2   Adsorbed_Ca
      -start
10 PUNCH TOT("Ca") MOL("Ca+2") MOL("LCa+2")
      -end
 
USER_GRAPH
      -chart_title "Langmuir Isotherm"
      -axis_titles "Ca solution (mol/l)"  "Ca sorbed(mol/l)"
      -start
10 PLOT_XY (MOL("Ca+2")),(MOL("LCa+2")
      -end
 
SURFACE_MASTER_SPECIES
      L     L
 
SURFACE_SPECIES
      L=L
      Log_k       0.0
      L + Ca+2 = LCa+2
      Log_k 2.0
 
SURFACE 1
      #-equilibrate with solution 1  #not needed for many solutions below
      -no_edl
      L     0.02                    # number of sites
END
 
SOLUTION_SPREAD
      -units mg/kgw
      -pH         8.0   
      # must be tab separated.. pH taken from above
Name  Ca    pH    Na    Sr    K     S(6)  Cl    Description
1     10          10700 15    300   2500  18000 charge      A
2     100         10700 15    300   2500  18000 charge      B
3     1000        10700 15    300   2500  18000 charge      C
4     3000        10700 15    300   2500  18000 charge      D
5     5000        10700 15    300   2500  18000 charge      E
6     8000        10700 15    300   2500  18000 charge      F
7     10000       10700 15    300   2500  18000 charge      G

END
 
USE solution 1; USE surface 1
END
 
USE Solution 2; USE Surface 1
END
 
USE Solution 3; USE Surface 1
END
 
USE Solution 4; USE Surface 1
END
 
USE Solution 5; USE Surface 1
END
 
USE Solution 6; USE Surface 1
END
 
USE Solution 7; USE Surface 1
END