Skip to content

Instantly share code, notes, and snippets.

@tadeu

tadeu/ph.m

Created June 13, 2013 22:28
Show Gist options

  • Select an option

    Learn more about clone URLs
  • Save tadeu/5777936 to your computer and use it in GitHub Desktop.

Select an option

Learn more about clone URLs
Save tadeu/5777936 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
ph.m
solve for (C_co2_o, C_h2s_o, C_h, C_hco3, C_co3, C_h2s_w), m_co2 = V_g * P_co2 + V_o * C_co2_o + V_w * (C_co2_w + C_hco3 + C_co3), m_h2s = V_g * P_h2s + V_o * C_h2s_o + V_w * (C_h2s_w + C_hs), alk_tot = C_hco3 + 2 * C_co3 + C_hs + C_oh - C_h, K1 = C_h * C_hco3 / C_co2_w, K2 = C_h * C_co3 / C_hco3, K_h2s = C_hs * C_h / C_h2s_w
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment