Jump to content
Geochemist's Workbench Support Forum

Vincent

Members
  • Content count

    5
  • Joined

  • Last visited

Community Reputation

0 Neutral

About Vincent

  • Rank
    Newbie
  • Birthday 06/17/1978

Contact Methods

  • Website URL
    http://comp.uark.edu/~vchevrie/

Profile Information

  • Gender
    Male
  • Location
    Fayetteville, Arkansas, USA
  1. Hello, In my laboratory we have developed a dynamic model using Matlab. Our code requires calculating some thermodynamic parameters like the concentration of some elements during evaporation and we thought about using GWB. However, I don't know how to create the pipes mentioned in the Reference Manual of GWB to be able to control GWB from our code. Does anyone has some experience in creating programs or routines that allow controlling GWB from Matlab and would be willing to share their experience? Also does someone knows if the helper functions provided for C++ or TCL exist for Matlab? Thank you very much Best regards Vincent
  2. "Admittedly reading the data in from database file and converting it to the right format might be more tricky (and could pose a compatibility problem with old database files, but perhaps one can just append the new coefficients after the existing ones and just read treat missing coefficients as zeroes)." The GWB thermodynamic databases based on the Pitzer model (Phrqpitz) can be modified to add new Pitzer coefficients. I have been doing this for a while and it work pretty well. "If this is programmed in a modular way with derivatives of the original T function propagating into activity coefficient (and osmotic coefficient) expressions, replacing the original T dependence and its derivatives by a new one should be no big deal actually." Sorry but I am not sure I understand what you mean. The program is defined not to read various values of the coefficients at various temperatures (the way the solubility constants are defined) but to read the fit coefficients X0, c1, c2, c3 and c4 and recalculate the parameters from these and the equation. This is why I have difficulties defining the parameters if the temperature dependency does not follow the predefined fit equation. Vincent
  3. Thanks for both answers. I agree it may look a bit complicated to redesign the program, I was just wondering if there was a simple solution. I also agree that 8 different terms can look a bit like over-complicating the situation. I realized I forgot to precise that not all of them are used. It is a general equation used by some author to fit the data. I think the highest number of parameters I have seen is 4. However, in some cases the T3 or the T4 or the 1/T2 were used. Therefore, even if I did not necessarily need the complete equation, I needed a way to define them in GWB. The whole idea is that I use compilations of parameters and their temperature dependency is based of this general equation. I could track down the original data, and do the fit with the GWB equation, but I preferred to save a lot of time and use these parameters. In any case it would be useful to have custom defined virial coefficients, just in case, and especially because their temperature dependency is always based on empirical equations. Best regards Vincent
  4. Hello, I am currently creating a new database, based on the phrqpitz database. In the database, the virial coefficients (beta_0, beta_1, beta_2, C_Phi, theta, lambda and psi) temperature dependency is defined by a polynomial equation of the form: X = X0 + c1*(T - T0) + c2*(1/T - 1/T0) + c3*ln(T/T0) + c4*(T2-T02) where T0 = 298.15K, T is the temperature in K and X0 is the value of the parameter at 25oC. For the coefficients I am using I need more parameters since the coefficients are fitted by a general equation of the form: X = X0 + c1*(T - T0) + c2*(1/T - 1/T0) + c3*ln(T/T0) + c4*(T2-T02) + c5*(T3-T03) + c6*(T4-T04) + C7*(1/T2-1/T02) Is there any way I can modify the equation used in GWB to adapt it to my database, rather than recalculating all the parameters to find a temperature dependency fitted by the standard GWB equation? Although this is feasible, it is very time consuming and often I can't find a standard equation fitting as well my general equation. Thank you very much in advance Vincent
  5. Hello, I'm relatively new in the use of GWB (version 6.0.3). I'm especially interested in React to trace the evolutions of evaporating solutions and as function of temperature. I know that for example it is possible to determine the precipitating parageneses with increasing temperature. However, my models aim at simulating surface processes on Mars. Therfore, I'd like to simulate the processes occuring under 0 C. I know that saline solutions at low temperatire have a lower eutectic temperature. In some cases some eutectic solutions can remain liquid down to -70 C (Ferric sulfate for example). I previously modelled such behaviors on simple salts using the Pitzer model. I know there is a Pitzer database working very well for evaporation but I would like to be able to use it in models at low temperature. Is there a way I could model/program this in the command section? Also two quick questions: - How can I calculate the water-to-rock ratio during evaporation other than manually by adding all the precipitated phases and dividing by the amount of water - How can I extract the water activity (for example as a text file) calculated using the Pitzer model after running a model? Thank you very much. Vincent
×