Click here to donate to keep PhreeqcUsers open

Quote from: Seven on July 26, 2017, 03:20:54 AMmy initial solution can be represented by:SOLUTION 1 Trail Creek - Irrigation temp 5.47 pH 7.35 pe 4 redox pe units ppm density 1 Alkalinity 149 Ca 44 Cl 0.73 K 0.38 Mg 11 Na 0.74 S(6) 8.8 -water 1 # kgand the final solution would beSOLUTION 2 DSP sampler temp 7.8 pH 6.7 pe 4 redox pe units ppm density 1 Alkalinity 379 Ca 130 Cl 1.2 K 0.65 Mg 2 Na 27 S(6) 21.8 -water 1 # kgsolution species would be:NH4+ + 3H2O = NO3- + 10H+ + 8e-NH4+ = NH3 + H+the hope is to plot the resulting mole concentrations similar to the Drew_2006 example.Reference to examples would be appreciated

my initial solution can be represented by:SOLUTION 1 Trail Creek - Irrigation temp 5.47 pH 7.35 pe 4 redox pe units ppm density 1 Alkalinity 149 Ca 44 Cl 0.73 K 0.38 Mg 11 Na 0.74 S(6) 8.8 -water 1 # kgand the final solution would beSOLUTION 2 DSP sampler temp 7.8 pH 6.7 pe 4 redox pe units ppm density 1 Alkalinity 379 Ca 130 Cl 1.2 K 0.65 Mg 2 Na 27 S(6) 21.8 -water 1 # kgsolution species would be:NH4+ + 3H2O = NO3- + 10H+ + 8e-NH4+ = NH3 + H+the hope is to plot the resulting mole concentrations similar to the Drew_2006 example.Reference to examples would be appreciated

SOLUTION_SPECIES# eliminate N2(aq)2 NO3- + 12 H+ + 10 e- = N2 + 6 H2O #-log_k 207.08 -log_k 0RATESNitrification 10 o2 = TOT("O(0)")20 amm = TOT("N(-3)")30 k = 1/8640040 rate = k*o2/(1e-5 + o2)*Amm/(1e-5 + amm)*amm50 SAVE rate*TIMESOLUTION-units mmol/kgwpH 7 chargeN(-3) 1C 1ENDUSE solution 1REACTIONO2 12.5e-3 in 10USER_GRAPH 1 -headings time O(0) Amm N(5) -axis_titles "Moles of O2 reacted" "Moles per kilogram water" "" -initial_solutions false -connect_simulations true -plot_concentration_vs x -start10 GRAPH_X RXN20 GRAPH_Y TOT("O(0)"), TOT("N(-3)"), TOT("N(5)") -end -active trueEND