Jump to content
Geochemist's Workbench Support Forum

Brian Farrell

  • Posts

  • Joined

  • Last visited

  • Days Won


Brian Farrell last won the day on August 16 2019

Brian Farrell had the most liked content!

Profile Information

  • Gender

Recent Profile Visitors

The recent visitors block is disabled and is not being shown to other users.

Brian Farrell's Achievements


Rookie (2/14)

  • First Post Rare
  • Collaborator Rare
  • Posting Machine Rare
  • Dedicated Rare
  • Week One Done Rare

Recent Badges



  1. Hi Polly, If I follow the benchmark while using thermo_phreeqc.tdat, I get almost the same results as PHREEQC-OGS. The only difference is a slightly wider zone of dolomite stability, and consistent with that, a few shifted characteristic points in the Mg++, Ca++, C, and pH curves. If you can find the PSI-NAGRA database used in the benchmark, your results might be even closer. You might only need to alter the log Ks of calcite and dolomite, or you could convert the entire PSI-NAGRA dataset from phreeqc to GWB format in TEdit if desired. Here are a few comments on your current input file. Calcite should be swapped into the basis on the Initial pane, as you’ve done, since it is in equilibrium with the pore water at the start of the simulation, but its abundance should be set to 2e-4 free mol/l to match the original simulation. Calcite should not additionally be set as a simple reactant. I don’t believe a non-default value for the diffusion coefficient is warranted, although its effect is negligible in this simulation. There is no reason to set an inert volume here. There is no reason to change the convergence criterion “epsilon” variable, as far as I can tell. You can set Nx to 100 to match the original simulation and to speed up program execution. I took care of these issues and didn’t have any convergence issues using thermo_phreeqc.tdat, but if you continue to have difficulty you can attach a revised version. If you want to overlay different results (mineral abundances, fluid composition, pH) in a single plot, you’ll have to do so in Excel, PowerPoint, Adobe Illustrator, or some other graphics program, as you’ve done. Finally, regarding your previously posted custom thermo dataset, you had 3 basis species SiO2(aq), Si(OH)4, and H4SiO4, along with an aqueous species H4SiO4(aq). Those species all represent the same thing, so there should only be one of them. I’d be very careful to clean up those entries and any dependent reactions so that you have a complete, consistent, and non-redundant set of reactions. Hope this helps, Brian Farrell Aqueous Solutions
  2. Hi Demetra, The ionic strength in most modern programs is calculated from and consistent with the distribution of species in solution. There is no way to override it in SpecE8. Please note, ionic strength here is not calculated from the component (bulk) concentrations. It accounts for the presence of ion pairs or complexes like MgCl+ and MgSO4 that cause the solution to have a lower “concentration of charge” than it would if the Mg++, Cl-, and SO4-- ions were fully dissociated. In other words, in ion association theory the activity reflects the activity coefficient as well as the distribution of mass. Hope this helps, Brian Farrell Aqueous Solutions LLC
  3. I have no specific insight about trevorite (perhaps someone else on the forum will), but in calculating equilibrium models it's important to remember that the most thermodynamically stable minerals might not be likely to form, especially within relatively short timeframes. Using your experience to suppress minerals that fall in this category is common in geochemical modeling. Hope this helps, Brian Farrell Aqueous Solutions LLC
  4. I believe that the dissolved Al concentration vs. Kaolinite equilibrium example is intended to demonstrate how alternative conceptual models might be set up in SpecE8, not to indicate a general rule to follow in speciation modeling. For more info, please see section 6.2 Amazon River water in Craig Bethke's Geochemical and Biogeochemical Reaction Modeling text. It's probably a good idea to first launch SpecE8 for a few of your samples and look at the entire set of calculation results, but if you have a large group of samples that are pretty similar, you're probably better off calculating a few analytes of interest from within GSS. So, instead of going to Analysis > Launch..., try going to Analysis > Calculate with Spece8..., then selecting your analytes of interest. In other words, you might look at the complete distribution of species in solution and the saturation indices for a large set of minerals for a few samples at first, but then when you know what you're looking for (e.g. S.I. of calcite) you can pick and choose to calculate just that analyte for all of your samples, and add it directly to your GSS file. Either way, you can specify trailer commands (for example, "swap Kaolinite for Al+++", "1 free g Kaolinite") in GSS that will be used for all samples when you launch SpecE8 or calculate values with SpecE8. Please see 3.3.7 Calculated values and 3.5 Launching SpecE8 and React in the GWB Essentials Guide as well as the SpecE8 chapter in the GWB Command Reference. Regarding convergence, I'd check to make sure your units and constraints are correct, first of all. I'd also check whether settings for the extent of redox equilibrium/disequilibrium to use in the calculation are reasonable (see 7.3 Redox disequilibrium in the GWB Essentials Guide). It can sometimes be helpful to swap in species that are expected to be dominant under the conditions of interest, to help the program converge. Al+++, the default basis species for Al, is likely to be dominant at low pH, but not at high pH values, in which case it might be better to swap in Al(OH)4-, for example. This last option is straightforward in SpecE8, but not GSS, so let us know if you'd like to pursue it further. Hope this helps, Brian Farrell Aqueous Solutions LLC
  5. Hi Erik, A user can specify under File > Preferences… in most apps (or the Settings tab of the GWB dashboard) the default thermo dataset to be something other than thermo.tdat. One can also specify a default dataset using an environment variable (see 2.3 Thermodynamic datasets and 2.8 Environment variables in the GWB Essentials Guide). Yes, a single .tdat file can be loaded at a time. If you incorporate surface chemistry in a calculation, however, you load one or more surface (.sdat) datasets in addition to the thermo dataset. Yes, that's correct. The input at a basic level consists of a description of the temperature and chemical composition of the system. Beyond constraining the original basis entries, you might use basis swaps to prescribe equilibrium with gases of known partial pressure or minerals in equilibrium with the fluid (e.g. swap O2(g) for O2(aq); P O2(g) = .2 bars). The pH setting is just a special constraint for the H+ basis entry. Additional options can be specified as well (constraints for isotopic calculations, suppressing certain species, controlling the extent of redox disequilibrium, and so on). You can find the complete set of input commands in the GWB Command Reference. Each app (e.g. SpecE8, React, ...) has its own chapter. You might also look at the example files included within the Script folder of the software. Each is described in one of the User’s Guides (e.g. SpecE8 in the GWB Essentials Guide, React in the GWB Reaction Modeling Guide). GWB won't be able to do anything with components it doesn't know about. So you could not tell GWB about them, or get an error if you do, I suppose. A user can edit the thermo dataset before a run if additional reactions are needed. The species and components and so on included in the output file will all be defined in the thermo dataset in use. The output also includes a few specially calculated parameters (e.g. total dissolved solids, carbonate alkalinity, fluid density, and so on) that are not defined in the thermo dataset. Hope this helps, Brian Farrell Aqueous Solutions
  6. Hi Geoff, I'm guessing you launched React at some point prior to mixing. GSS passes the "precip off" command to React in this case by default, so that calculations more closely match those you might make within the spreadsheet, which doesn't account for precipitation (SpecE8 does the work within GSS). The command is passed via a "trailer command", which you can see in the Analysis > Launch... or Analysis > Options... dialogs. The trailer commands (you might choose a specific ion for charge balancing, for example) in general are saved for other calculations within GSS, like calculating analytes and mixing, but unfortunately this command is valid only for React. You can go to Analysis > Options... before you begin mixing to remove the "precip off" command. Hope this helps, Brian Farrell Aqueous Solutions
  7. Hello Jerome, I realize this post is quite old, but I happened to come across it and realized some updated information might be useful. In your post from March 1, 2019, you show formation reactions for O2 and H2 and their log Ks, as taken from Thermoddem's phreeqc-format dataset (the first two of three reactions you provide). You claimed these are for formation of the gaseous species, but they are undoubtedly for aqueous species. I couldn't comment on the contents of the phreeqc data at the time, but I am more familiar with phreeqc datasets now. The first two reactions you supplied are within the SOLUTION_SPECIES block and contain the "CO2_llnl_gamma" keyword, which indicates the method to use for calculating activity coefficients of aqueous species. Reactions for gas dissolution appear within a separate PHASES block. Hopefully that clears up the matter. FYI, the TEdit app in GWB 2021 can now read in phreeqc datasets and convert them to GWB format, for use with the rest of the GWB programs. You might wish to compare reactions balanced in Rxn or pe-pH diagrams constructed in Act2 with the native GWB Thermoddem dataset and that converted from phreeqc. Regards, Brian
  8. Dear GWB users, A problem was detected reading in thermo dataset entries that were edited by hand in a text editor when tabs are used in place of spaces in an entry, such as between a mineral's name and its type: Calcite type= carbonate formula= CaCO3 mole vol.= 36.9 cc mole wt.= 100.0892 g The user would receive an error like "Could not suppress Calcite" or "Could not identify reactant species Calcite". The issue has been fixed and a release candidate installer is available to get affected users up and running right away: GWB 15.0.1rc1 installer The fix will eventually be available in an official maintenance release. At that time, users with the release candidate installed will be automatically notified of the update. If you encounter any more issues, please notify us here. We appreciate your patience and hope you enjoy using The Geochemist's Workbench. Sincerely, Brian Farrell Aqueous Solutions LLC
  9. Hi Karen, We've made a fix to the problem you've encountered. Please download and install this "release candidate" and try using your customized thermo dataset. 15.0.1rc1 installer I don't know why X1t would crash for you, but if it happens again you might try reinstalling and restarting your computer. And you can also try recompiling your rate law function after reinstalling 15.0.1rc1. Please let us know if you have any more issues and we'll work to fix them as soon as possible. Regards, Brian
  10. Hi Karen, We put a lot of work into adding flexibility to thermo datasets and improving TEdit, but unfortunately we introduced a problem handling minerals that include tabs, rather than spaces, in their entry. Assuming you're using the same thermo dataset that you had trouble modifying previously, I've gone to your Na2EX2 and CaEX2 additions and simply replaced the tab characters between the name and type fields with spaces. If you have made more modifications, you can make a copy and replace the tabs with spaces by hand, or do it automatically with an editor like Notepad++ (go to Edit > Blank Operations > TAB to Space). That should take care of your problem while we work on a software fix. In the future, we'll be able to diagnose issues like this much more quickly if you actually provide example datasets and input files that demonstrate the problem. And we'll appreciate it if you can limit your comments to a single post, and understand that it takes time to see, read, work through, and respond to queries. We apologize for the inconvenience and appreciate your patience as we work to rectify this issue. Regards, Brian Brian Farrell Aqueous Solutions thermo_ionx_TabsToSpaces.tdat
  11. Dear GWB12 users, We are pleased to announce our latest maintenance release for GWB12, GWB 12.0.7. The 12.0.7 update provides a fix to calculated analytes in GSS when temperature is given in Fahrenheit; fixes a bug saving Pcrit, Tcrit, and omega in TEdit; and provides fixes for all known issues. Update from 12.0.0-12.0.6 at no charge to ensure you have all the newest features and bug fixes. Existing installations should automatically update to this release, unless auto-update is disabled. In that case, users should update their installations from the Help menu of any GWB app. Regards, Brian Farrell Aqueous Solutions
  12. Dear GWB users, We have good news to start the year: GWB 2021 is here! GWB 2021 brings SIT activity coefficients, an all new thermo data editor, the ability to import thermo datasets from PhreeqC, and much more! Your GWB subscription will upgrade automatically within three days, or simply click “Check for updates” on the GWB dashboard. Not subscribing yet? Take advantage of one of our flexible plans or upgrade that old GWB license to a modern subscription. Contact us today for a quote, or go shopping at our online store and you’ll be up and running in minutes. Sincerely, Brian Farrell Aqueous Solutions LLC
  13. Hi DI Xie, Can you please provide a screenshot of the error message that you receive when you try to deactivate? Thanks, Brian Farrell Aqueous Solutions
  14. Dear GWB users, We completed some updates for the forum software earlier this week. If you were not able to post while it was undergoing maintenance, please try again now. Regards, Brian Farrell Aqueous Solutions LLC
  15. Hi, Act2 only displays samples whose temperature is close to that of the diagram you've calculated (+/- 10 C). The samples at 76 C plot because the diagram is drawn at 76 C, but the samples that are 30-40 C will not plot. Regards, Brian Farrell Aqueous Solutions LLC
  • Create New...