Jump to content
Geochemist's Workbench Support Forum

Brian Farrell

Admin
  • Posts

    1,060
  • Joined

  • Last visited

  • Days Won

    11

Everything posted by Brian Farrell

  1. Hello, Thanks for providing your .dis and .bud files. X2t currently expects CONSTANT DELR, DELC, TOP, and BOTM. Your script has varying values for TOP, which is giving it trouble. It looks like the Reactive Transport Modeling Guide specifies that delta x and delta y in X2t need to be constant (taken from your DELC and DELR), but not the TOP and BOTM. We'll get the appendix updated. By the way, I also noticed your .dis file has multiple layers. X2t can only use a single layer, but as mentioned in the appendix, it will pick the uppermost layer. Hope this helps, Brian Farrell Aqueous Solutions LLC
  2. Dear GWB users, We are pleased to announce our latest maintenance release, GWB 12.0.3. The 12.0.3 update features support for Windows 10 2018 Creators Update; smaller installer size; better font and symbol support for east Asian locales; various improvements to plotting apps; correction of a glitch affecting some influx/efflux calculations in X1t and X2t; correction to an auditing issue when modeling retardation within stagnant zones; and fixes to a few issues in GSS in which copy/paste operations involve an “as” unit conversion. Update from 12.0.0 - 12.0.2 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
  3. I'm glad to hear you're all set. By the way, you should now be able to post on the main GWB forum. Regards, Brian Farrell Aqueous Solutions LLC
  4. Hi Ellie, There are two aspects to this question. First, in any run in which you have supersaturated minerals, React will report calculations for the metastable fluid (the first Xi = 0) block, as well as the "true" equilibrium state, in which supersaturated minerals have been allowed to precipitate (the second Xi = 0 block). Section 2.3 Initial system in the GWB Reaction Modeling Guide describes the differences between these two blocks of output, and suggests using input from a SpecE8 example (seawater.sp8) in React to compare the differences. In that example, the precipitation of dolomite and quartz from seawater alters the pH from its initial value. The scope of your modeling work will dictate which block of results is more appropriate. Second, if you trace an actual reaction path of some sort, for example, by titrating mass as a simple reactants, sliding the temperature, or sliding gas fugacity, React will report in its text output file results from various steps in the calculation (Xi = .1, Xi = .2, ..., Xi = 1). If you slide temperature from 0 to 100 C, for example, Xi = 0 corresponds to the initial system, Xi = .1 corresponds to the system heated to 10 C, and Xi = 1 corresponds to the end of the simulation, the system at 100 C. Understanding the entire reaction path is of course important, but often people care only about the results at the end of the calculation. Averaging values from different steps would not be useful. I assume these steps are what you're referring to, not the actual iterations within any single reaction step. Please let me know if you have any more questions. Regards, Brian
  5. Hi Odeta, I’m including some bibliographic information from the GWB’s thermo_minteq.tdat dataset below: This dataset is the thermodynamic database from Visual MINTEQ release 2.40, for use with the GWB programs. The dataset has been compiled from files thermo.mdb, analyt.mdb, and comp.mdb. The Visual MINTEQ database is maintained by Jon Petter Gustafsson, KTH (Royal Institute of Technology), Stockholm, Sweden. The conversion code was written by Daniel Saalfeld and Craig Bethke, University of Illinois, September 2003. This version was compiled in December 2005 by Jon Petter Gustafsson. It’s not entirely clear to me whether the results of this compilation were merely recommendations or were actually included in an update to Minteq. And as far as I’m aware, Visual Minteq diverged from the original Minteq some time ago. In any case, we are not specialists in thermodynamic data and we don’t change our thermodynamic datasets from release to release, apart from adding support for new features. Our datasets are fully editable by our users. You are free to modify any of the equilibrium constants in the datasets, add new reactions, adjust redox coupling, and so on. We built a graphical editor called TEdit for modifying the thermo datasets, but you can alternatively use any text editor. For more information on TEdit, please see Chapter 9 in the GWB Essentials Guide. Regards, Brian Farrell Aqueous Solutions
  6. Hi Ellie, No worries. That's what I was going to suggest. It sounds like you're on the right track. Regards, Brian
  7. Hi, Before you dive in too deep with your own data, I think it would be a good idea to practice some of the skills you’ll need using examples from the User’s Guide. I think the most important parts of the problem are: -defining the initial system (Sections 7.1 and 7.2 in the GWB Essentials Guide; Sections 2.1 and 2.3 in the GWB Reaction Modeling Guide) -redox disequilibrium (Section 7.3 in the GWB Essentials Guide) -sliding fugacity paths (Section 3.5 in the GWB Reaction Modeling Guide) -picking up chemical systems (Section 3.9 in the GWB Reaction Modeling Guide) -sliding temperature paths (Section 3.4 in the GWB Reaction Modeling Guide) As a first step, I would define a fluid containing Fe++, HS-, Ni++, H+, and so on, with a negligible amount of HCO3-, then verify that you precipitate your desired mineral. Next, you can fold in your CO2 addition, pick up the results, and then use those as the starting point for a sliding temperature path (the GWB does indeed account for the effects of temperature on pH). Hope this helps, Brian Farrell Aqueous Solutions LLC
  8. Dear Yanlu, Thanks for explaining your situation and for providing the error messages you’ve received. I just sent you a file that should reset your activation. If you need to clean your OS in the future, you can deactivate your license beforehand, then reactivate when you're done. That will be a little easier than running the reset file I sent. Regards, Brian Farrell Aqueous Solutions LLC
  9. Hi Ben, You’re very welcome. One idea would be to add a kinetic rate law for any mineral present but not in equilibrium with the initial system. Of course, you’d need meaningful kinetic parameters. If you prefer the “backwards” configuration, you can do that too. In any titration path (minerals added as simple reactants to a fluid, for example) the reactant mass is added incrementally throughout the calculation (100 steps by default). So if you react 1 kg of a mineral with ~1 kg of solution, you’ll first add 10 g, then another 10 g, and so on, until you’ve reacted the entire mass. In this way, you’re effectively starting at an infinite water:rock ratio and decreasing it to 1 at the end of the reaction path. And as before, you can plot the amount of mass reacted on your x or y axis. This is obviously different conceptually from the previous suggestion, but I think the term leaves a lot of room for interpretation. I'm not positive how EQ3/6 conceptualizes the problem, but after a quick search I found an excerpt from Low-Grade Metamorphism of Mafic Rocks edited by Peter Schiffman and Howard W. Day. The description on page 90 seems pretty similar to what I've described in this reply. Hope this helps, Brian
  10. Dear Avery, The thermo_ymp.R2.tdat dataset was first made available for the GWB12 release. As such, it follows the “jul17” dataset format introduced with GWB12. The format adds support for including formulae for aqueous species; includes factors for calculating fugacity coefficients for the gas species; and designates a default method for the calculation of fugacity coefficients and partial pressure. The dataset will not work with GWB releases before 12. For more information on thermo dataset formats, please see the “Thermo Datasets” chapter of the GWB Reference Manual. The “Legacy formats” section should be especially relevant. Regards, Brian Farrell Aqueous Solutions LLC
  11. Dear Benjamin, There are probably a variety of ways to conceptualize the problem, but combining React’s “flush” configuration with the “reactants times” feature is probably one of the more elegant ways to do so. A flush model is traced from the reference frame of the rock through which a fluid migrates. The migrating fluid displaces existing fluid from the system. Please note there is no “water-to-rock ratio” variable in the plot output. However, the “H2O reacted” (solvent) or “Mass reacted” (solvent + solutes) variables available under Reactants properties (or Chemical parameters for GWB releases before 12) are probably the best ways to describe the progress of reaction. If you export the plot as an enhanced metafile to PowerPoint, it’s a simple matter to adjust the label. For more information and an example of the flush configuration, please see section 3.3, Flush model, in the GWB Reaction Modeling Guide. More good examples can be found in section 3.4, Polythermal reaction paths, of the Reaction Modeling Guide, as well as section 19.4, Dolomitization of a limestone, in Craig Bethke’s Geochemical and Biogeochemical Reaction Modeling text. Hope this helps, Brian Farrell Aqueous Solutions LLC
  12. You are very welcome. I hope it goes well. Brian
  13. You should be able to plot user defined analytes on XY plots, but the special plots only accept true basis species (e.g. B(OH)3, Cl-, Li+) and in some cases chemical parameters (e.g. pH and TDS in the Durov diagram). Thanks for your suggestions to improve our software. We receive many requests from our users (some more feasible than others), so it's hard to keep everyone happy, but we can certainly add your request to the list to consider as we plan the next release. Regards, Brian
  14. Hello, Currently there is no way to add multipliers on a ternary diagram, apart from manually adjusting the concentration values in your GSS spreadsheet. Regards, Brian Farrell Aqueous Solutions LLC
  15. Hi Yongtae Kim, The Dzombak and Morel two-layer model is a very thorough compilation of cation and anion sorption reactions involving hydrous ferric oxide. It is in no way complete, however. The FeOH.sdat dataset included with the GWB contains reactions and log Ks that Dzombak and Morel derived from actual experiments, while the FeOH+.sdat includes some estimates for additional sorbing species. I believe there is a reaction for SiO2 sorption in the FeOH+ dataset, but not in the FeOH dataset, for example, because they didn’t have any experimental data available. You should keep in mind that the Dzombak and Morel datasets are examples. The mineral they studied (and its surface properties) might be close to the iron mineral in your experiments, but they are not necessarily the same. In your models, the SiO2(aq) and HCO3- components do not significantly decrease the efficiency of arsenic removal. As I mentioned above, the FeOH.sdat dataset contains no sorption reactions for SiO2, while neither FeOH.sdat nor FeOH+.sdat include any reactions for HCO3- species. They do, however, have reactions for As(III) species and PO4---, which is why you observe both of those components to sorb. You can view the datasets yourself using the TEdit app. TEdit can also be used to add reactions that you need. You’ll need to supply a reaction for each surface complex along with an equilibrium constant, expressed as a log K. There should be information for sorption of HCO3- to iron oxyhydroxides in the literature, and there may be additional information for SiO2 as well. Or, you can perform sorption experiments and derive log Ks yourself. For more information on TEdit, please see chapter 9 Using TEdit in the GWB Essentials Guide. If you have a redox equilibrium model enabled (the default), the arsenic will be present as both As(V) and As(III), but mostly As(V), reflecting the high DO concentration you’ve specified. If you believe that the arsenic should be present in only one valence, you should decouple the redox reactions involving arsenic and add only the appropriate basis entry. For more information, please see 2.4 Redox couples and 7.3 Redox disequilibrium in the GWB Essentials Guide. By the way, if your reactor remains in contact with the atmosphere throughout the experiment, you should fix the fugacity of O2(g) in your model (and CO2(g) as well, when you include the HCO3- component). If it is closed to the atmosphere, then you most likely should not fix the fugacity of these gases. Please see section 3.5 Buffered paths in the GWB Reaction Modeling Guide for more information. Hope this helps, Brian Farrell Aqueous Solutions LLC
  16. I'm glad to hear that you can draw the water stability limits now. I hope you enjoy using the software. Brian
  17. Hi Xinyu, Act2 needs reactions for O2(g) and H2(g) to draw the water stability limits. The thermo_minteq.tdat dataset does not contain a reaction for H2(g), however, so you’ll have to add it. Because the reaction for H2(g) is typically written in terms of H2(aq), there are actually two reactions you need to add: H2(aq) + .5 O2(aq) = H2O (in the redox couples section) H2(g) = H2(aq) (in the gases section) These are the reactions you’ll find in the GWB’s default dataset, thermo.tdat. Each of these reactions has equilibrium constants (expressed as log Ks) at multiple principal temperatures. The principal temperatures differ between the two datasets, so you’ll need to make sure to specify the log Ks for the correct temperatures. In older versions of the GWB, you had to edit the datasets by hand in a text editor like Notepad, so this was easy to overlook. Now you can simply copy and paste (or drag and drop) an entry from one dataset to another using the TEdit app. Best of all, TEdit will automatically find matching principal temperatures and insert the log Ks in the correct locations. For more on using TEdit to copy reactions, please see 9.2.6 Transferring dataset entries in the GWB Essentials Guide. For an older discussion of a similar topic, please see this thread. Hope this helps, Brian Farrell Aqueous Solutions LLC
  18. Hi Mauricio, The GWB does not currently support surface precipitation as developed in Dzombak and Morel’s generalized two-layer model. Regards, Brian Farrell Aqueous Solutions LLC
  19. Hi Alero, After thinking about this some more, you might be able to put together something like you described with a modified thermo dataset. I'd start by adding h+ as a redox couple, using the reaction h+ = H+ + .25 O2(aq) - .5 H2O. Set the charge of h+ to 1, its mass to that of an electron, and make it really unstable. In React, set up the chloride oxidation problem with a kinetic rate law. Decouple the h+ as well as the ClO4-, set very small initial amounts of h+ and ClO4-, then enter the reaction above for the kinetic redox reaction and assume a rate constant. If you add the h+ as a simple reactant, you should observe the Cl- oxidizing to ClO4- and pH decreasing with time, like in the reaction you described. Hope this helps, Brian
  20. Dear GWB users, We are pleased to announce our latest maintenance release, GWB 12.0.2. The 12.0.2 update features improvements to special plots; addresses issues with Xtplot contours and color maps; updates and fixes GWB Plugin feature and documentation, including support for Visual Studio 2017; adds fixes for ChemPlugin wrappers for Fortran, Java, and MatLab; and provides fixes for all known issues. Update from 12.0.1 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
  21. Hi Coralie, You can overlay contours of a single variable on a predominance map in P2plot. Currently you can plot a predominance map of one basis entry or element at a time, however. If you want to overlay the predominance bounds from other basis entries or elements you will need to do so in another program like PowerPoint or Illustrator. Regards, Brian
  22. Hi Alero, The software does not currently support photoelectrochemical reactions like that you've described. Regards, Brian Farrell Aqueous Solutions LLC
  23. I'm happy to hear it helped. Cheers, Brian
  24. Hi Coralie, There are various ways to interact with P2plot and the other plotting apps. You can double-click on a specific plot aspect, like the "Rxn progress (x)" axis variable, to open the "Axis Range and Variables" dialog. Or, you can right-click on "Rxn progress (x)" and pick "X Axis Range" to open the same dialog. Finally, you can use the menubar along the top of P2plot to access the dialogs. Under Format, choose "Axis Range and Variables...". From that dialog, you can pick your axis variables, set linear or log scales, tick increments, units, and so on. You can add contours from Format -> "Contour...". I recommend reading the Phase2 and P2plot chapters in the GWB Reaction Modeling Guide, which you can access from the Help menu of any GWB app. At the very least, you should check out the first couple examples (7.8 Example: Speciation diagram and 7.9 Example: Mineral solubility) in the guide to see the process of calculating and rendering the diagrams. The thermo dataset that you set in the Preferences dialog will become the default dataset every subsequent time you open a GWB app, but it won't load that dataset into an instance of the app that is already open. To do that, you would go to File -> Open -> Thermo Data... You can verify at any time the dataset that is currently loaded by going to File -> File Properties -> Thermo Data. Hope this helps, Brian
  25. Hi Coralie, It sounds like you're using Act2 and Tact for your diagram construction, but I think what you need is Phase2, a new app we developed for GWB12. From the GWB Reaction Modeling Guide: A Phase2 diagram differs from a simple diagram of the type constructed by Act2 and Tact in that each point in the diagram represents the complete solution to the equations describing the system’s distribution of mass. If you were to use Phase2 to calculate an Eh-pH diagram, for example, you would find that, unlike the result from Act2, the boundary lines are curved, rather than linear. Because of the calculation’s completeness, some Phase2 diagrams differ qualitatively from their Act2 counterparts. For this reason, the diagrams are sometimes referred to as “true” Eh-pH or “true” activity diagrams. Phase2 diagrams, furthermore, can be plotted over a wide range of variable choices, and in a variety of ways. You can plot the results of Phase2 calculations as maps of species predominance or mineral assemblages, and you can color map, mask, or contour any variable. You can additionally render cross-sections through the diagrams. Because a Phase2 calculation is conceptually similar to React, you can use any units to constrain your chemical system, and you can choose a variety of axis variables not possible in Act2 or Tact. For more information, please see our Phase2 webpage and the documentation on Phase2 and P2plot in the GWB Reaction Modeling Guide. Regards, Brian Farrell Aqueous Solutions LLC
×
×
  • Create New...