Jump to content
Geochemist's Workbench Support Forum

Tom Meuzelaar

Members
  • Posts

    486
  • Joined

  • Last visited

  • Days Won

    1

Posts posted by Tom Meuzelaar

  1. Hi Brian:

     

    I was unable to remove the sodium basis species from the cation exchange template Db (IonEx_revised.sdat) using TEdit. After deleting the default basis species (which also removed the two basis species), I was unable to add a new basis species (calcium)- when I select 'Ca' and click apply, nothing happens. The same behavior for any other species.

     

    I am using 11.0.5.

     

    I was able to manually edit the DB, so it's fine, but I thought I should report it.

     

    Tom

  2. Thanks Brian- I think we'll wait for the patch. We've got neutral species for numerous metals; I tried a quick global search and replace in the Minteq thermo dataset- X1T does indeed accept the revised thermo database and the surface dataset with neutral species, but renaming the revised minteq thermo database causes errors with how the alkalinity constraint is applied. I could leave the database name untouched, but then we have problems that will cascade through a number of our other models.

     

    Tom

  3. Hi Brian:


    I'm modeling simple linear selenium sorption in X1T using v. 11.0.4, and the Minteq database. I've created a surface dataset (attached 'Kd_modified...' with Kd values for several selenium species. It works fine until I add a neutral species with '(aq)' suffix. Then I get the attached error- this occurs in both React and X1T.


    I haven't tried replicating with thermo.dat, but I'm wondering if the extra space between species and '(aq)' in Minteq is causing the problem. Perhaps it's something else I've overlooked.


    Thanks in advance,


    Tom

    post-8-0-40181300-1466195274_thumb.png

    Kd_modified_Se_w_neutral_species.sdat

  4. Hello:

     

    This isn't necessarily an issue with charge balance. It could be that you've defined a system that is far from equilibrium, and thus has trouble converging. If the pH value is important, I recommend starting with 3.5, so the system initially converges, then sliding the pH (in the Reactants pane) down to 3.1, so that you have a stable chemistry.

     

    Hope that helps,

     

    Tom Meuzelaar

    RockWare, Inc.

  5. Hi:

     

    The time function in React is primarily used for kinetic mineral dissolution and precipitation. Without specifying a kinetic mineral rate law, it doesn't matter what time length you specify, you'll always get the same reaction results.

     

    Hope that helps,

     

    Tom Meuzelaar

    RockWare, Inc.

  6. Hi Olivia:

     

    We are wondering if a firewall or anti-virus program are temporarily disabling the software that runs your hardware lock. Is it possible to disable your firewall and/or antivirus software to see if that makes a difference?

     

    If it does not, can you go to the following address on your browser: http://localhost:1947, and click on the HASP keys menu on the left hand side? Please send the resulting screenshot to my attention via email.

     

    Regards,

     

    Tom Meuzelaar

    RockWare, Inc.

  7. Hi Ken:

     

    I'd start by taking a look at polythermal reaction path models in the Reaction Modeling User's Guide. You should be able to view the pH vs. temperature trends in Gtplot. Dealing with variable pressures is going to be a bit trickier, given that the database is compiled for a single confining pressure (1 atm). For most aqueous species and minerals, the equilibrium constant data does not change significantly, but for gases it is a different story altogether, and you may need to make some database modifications.

     

    Hope that helps get you started,

     

    Tom Meuzelaar

    RockWare, Inc.

  8. Hello:

     

    Yes, I've used the RC_helper.h and RC_helper.cpp source files without too much trouble using Visual Studio. I had to adjust a few settings in the default setup for VS, including:

     

    1. In project settings, remove explicit use of the Unicode character set.

     

    2. Turn off use of precompiled headers

     

    However, these settings are specific to Visual Studio, so they may or may not help you. I can tell you that there is nothing wrong with the RC_helper.h and RC_helper.cpp code, and that your problem is probably somewhere in the configuration of the application you are using for code development.

     

    Unfortunately, my response time will be much slower after this, as I am about to step out for a 4-day holiday weekend.

     

    Regards,

     

    Tom Meuzelaar

    RockWare, Inc.

  9. Hello:

     

    I recommend checking your syntax- as the error messages suggest,this is likely cause by missing a semicolon somewhere. If this occurs in one line of code, all of the other lines of code can become uninterpretable.

     

    I hope that helps,

     

    Tom Meuzelaar

    RockWare,Inc.

  10. Hi Ruth:

     

    This should work fine- just be sure to:

     

    • constrain the oxygen concentration in the system
    • specify a pH
    • convert your total As concentration to the arsenic basis species

     

    GWB will speciate all of your arsenic among all dissolved arsenic species in the thermo database according to their equilibrium constants.

     

    Let me know how I can be of further help.

     

    Regards,

     

    Tom Meuzelaar

    RockWare, Inc.

  11. Dear Christoph:

     

    The reason for this is that the temperature dependence of the rate constant is defined as a function of the activation energy (Ea) and pre-exponential factor (A) as follows:

     

    kinetic_temperature_dependence.jpg

    If you have a kinetic rate law that includes both a rate constant and activation energy, you might consider using the Custom Rate Law features in GWB.

     

    Hope that helps,

     

    Tom Meuzelaar

    RockWare, Inc.

  12. Hi Emma:

     

    React can work with waters that have negative pH (see attached script), but I don't think it will take a negative pH as input. Instead, use the equivalent proton activity to constrain your model.

     

    I hope that helps.

     

    Tom Meuzelaar

    RockWare, Inc.

  13. Hi Ross:

     

    Yes, your rate law evaluates to zero in the output file. I can't spot the issue right off- troubleshooting rate law scripts can tedious and time-consuming. What I recommend is that you simplify the script and first test your neutral, acid and base rates individually. Then put them together into your combined rate law.

     

    Hope that helps,

     

    Tom Meuzelaar

    RockWare, Inc.

  14. Hi Adam:

     

    You can do this using React, a simple mixing model, and the Pickup command. First add your acid groundwater to the React Basis. Speciate this water, and use the Config - Pickup - Reactants command to load the speciated water into the Reactants tab. Next load your surface water into the Basis tab. You now have a model simulating gradual titration of acidic groundwater into surface water.

     

    You'll run into a problem due to the fact that your surface water does not have all the components of your acid groundwater (Br-, Al, Fe, Mn). Remedy this by either excluding these components from your acid groundwater analysis, or by adding acceptably low concentrations for these components to your surface water in the Basis.

     

    I hope that helps,

     

    Tom Meuzelaar

    RockWare, Inc.

  15. Hi Ross:

     

    It turns our that the rate law script, as written in the v8 Reaction Modeling user's guide (page 56), reflects command line syntax. The actual code for the direct script editor looks exactly like the code you would put in a .bas file- you only need the code in between the brackets- not the 'kinetic Albite rate_law = { }' part.

     

    I hope that makes sense...

     

    We'll update the User's Guide to be a little more clear.

     

    Regards,

     

    Tom Meuzelaar

    RockWare, Inc.

×
×
  • Create New...