Ka 10^-3.49*1[M] "Acid equilibrium constant, H2CO3<->H + HCO3"
KD 715 "CO2 hydration equilibrium constant"
pKa_H2CO3_1 -log10((Ka/1[M])/KD) "pKa, first dissociation step"
pKa_H2CO3_2 10.3 "pKa, second dissociation step"
pKa_H3PO4_1 -log10(7.5e-3) "pKa, first dissociation step"
pKa_H3PO4_2 -log10(6.2e-8) "pKa, second dissociation step"
K_H3PO4_1 10^-pKa_H3PO4_1*1[M] "Acid equilibrium constant, first dissociation step"
K_H3PO4_2 10^-pKa_H3PO4_2*1[M] "Acid equilibrium constant, second dissociation step"
K_H2CO3_1 10^-pKa_H2CO3_1*1[M] "Acid equilibrium constant, first dissociation step"
K_H2CO3_2 10^-pKa_H2CO3_2*1[M] "Acid equilibrium constant, second dissociation step"
pH0 7.4 "Inlet pH"
cH0 (10^-pH0)*1[M] "Inlet H concentration"
cMg0 1.5[mM] "Inlet Mg concentration"
pKw 14.94-0.04209*(-273.15+T)/1[K]+1.718E-4*(-273.15+T)^2/1[K^2] pKw
Kw (10^-pKw)*(1[M])^2 "Water autoprotolysis equilibrium constant"
cOH0 Kw/cH0 "Inlet OH concentration"
cHCO30 25[mM] "Inlet HCO3 concentration"
cCO20 cH0*cHCO30/K_H2CO3_1 "Inlet CO2(aq) concentration"
cCO30 cHCO30*K_H2CO3_2/cH0 "Inlet CO3 concentration"
cCO2tot cHCO30+cCO20+cCO30 "Total concentration of CO2 species"
cHPO40 1[mM] "Inlet PO4 concentration"
cH2PO40 cHPO40*cH0/K_H3PO4_2 "Inlet H2PO4 concentration"
cH3PO40 cH2PO40*cH0/K_H3PO4_1 "Inlet H3PO4 concentration"
cH3PO4tot cHPO40+cH2PO40+cH3PO40 "Total concentration of H3PO4 species"
cAux0 -(2*cMg0+cH0-cH2PO40-2*cHPO40-cHCO30-2*cCO30-cOH0) "Inlet concentration of remaining ions"
T 37[degC] Temperature