Jump to content
Geochemist's Workbench Support Forum

Rob

Members
  • Posts

    2
  • Joined

  • Last visited

Everything posted by Rob

  1. Rob

    Surface complexation

    Hi Tom, I am still having problems with my model. When I model the system in React, with sorbing surfaces included, it crashes as one iron solid phase is depleted and the second one approaches its maximum. When the model is run without sorbing surfaces, it runs to completion, but at the point where all of the ferrihydrite is gone there is a change in pH and Eh that the moeld with sorbing surfaces doesn't get through, i.e. it crashes at that point. I have attached the sorbing basis species and surface species I am trying to model as a .txt file. The general model is simple, I am trying to model adsorption of As on two iron surfaces under decreasing redox potential. The log Ks are from the literature and are consistent with the thermo.com.v8.r6+ database. When I model this in React, I start with 5 free grams of Fe(OH)3 and titrate 50 mmol acetic acid to drive the Eh to reducing conditions. My React script is include at the bottom of the attached .txt file. To get Fe(OH)3 as the metastable phase I have to suppress the other iron oxides. Do you have any suggestions that can help the model get through the point where it crashes? rroot_REACT_061509.txt
  2. I am using React to model the sorption of arsenic on ferrihydrite and magnetite as conditions become reducing using constants from Dixit and Hering (2003). I have split the iron surface constants as ferrihydite >FeOH etc. and magnetite as >MOH etc. Can React handle 2 different mineral surface sites? If SCM is omiited, the model runs as predicted, i.e. ferrihydrite reductively dissolves and magnetite becomes the stable iron solid. However if I include SCM the model crashes as ferrihydrite dissolves and magnetite becomes the stable mineral. This happens even if I remove all possible sorbing ions.
×
×
  • Create New...