DATABASE c:\phreeqc\database\T_H.DAT # from www.hydrochemistry.eu # Cd sorption on OC in Christensen's xpt PRINT; -reset false; -user_print true SURFACE 1 # 1 g loamy sand = 0.007 g HA, multiply site conc's for 1 g HA with 7e-3 # Specific surface = 46.5e3 m2/g for 1 mM CaCl2 (I = 0.003 mol/L) H_a 4.97e-6 46.5e3 0.007 H_b 4.97e-6; H_c 4.97e-6; H_d 4.97e-6 H_e 2.483e-6; H_f 2.483e-6; H_g 2.483e-6; H_h 2.483e-6 H_ab 8.278e-7; H_ad 8.278e-7; H_af 8.278e-7; H_ah 8.278e-7 H_bc 8.278e-7; H_be 8.278e-7; H_bg 8.278e-7; H_cd 8.278e-7 H_cf 8.278e-7; H_ch 8.278e-7; H_de 8.278e-7; H_dg 8.278e-7 -equil 1 -diffuse_layer SOLUTION 1 pH 6.0 Ca 1; Cd 1e-6; Cl 2 charge -water 0.05 USER_PRINT -start 10 H_Cd = mol("H_aCd+") + mol("H_bCd+") + mol("H_cCd+") + mol("H_dCd+") 20 H_Cd = H_Cd + mol("H_eCd+") + mol("H_fCd+") + mol("H_gCd+") + mol("H_hCd+") 30 H_Cd = H_Cd + mol("H_abCd") + mol("H_adCd") + mol("H_afCd") + mol("H_ahCd") 40 H_Cd = H_Cd + mol("H_bcCd") + mol("H_beCd") + mol("H_bgCd") + mol("H_cdCd") 50 H_Cd = H_Cd + mol("H_cfCd") + mol("H_chCd") + mol("H_deCd") + mol("H_dgCd") 60 print " ug Cd/g =", H_Cd * 112.4e6 * 0.05,\ " ug Cd/L =", tot("Cd") * 112.4e6,\ " Kd (L/kg) =", H_Cd / tot("Cd") * 50 70 print " ug Cd/g in DL =", edl("Cd") * 112.4e6, chr$(13) 80 dl_Ca = edl("Ca")*2 - edl("water")*tot("Ca") * 2 90 dl_Cl = edl("Cl") - edl("water")*tot("Cl") 100 print "Excess meq Ca in DL =", dl_Ca 110 print "Excess meq Cl in DL =", dl_Cl, "Ca - Cl charge in DL =", dl_Ca-dl_Cl 120 print "Surface charge =", edl("Charge") -end END