A Computer Program to Solve Aqueous Inorganic Equilibrium Systems Using Pitzer's Method to Determine Activity Coefficients

Total Page:16

File Type:pdf, Size:1020Kb

A Computer Program to Solve Aqueous Inorganic Equilibrium Systems Using Pitzer's Method to Determine Activity Coefficients AN ABSTRACT OF THE THESIS OF Barry Kelly for the degree of Master of Science in Chemical Engineering/Forest Products presentedon August 12, 1988. Title: A Computer Program to Solve Aqueous Inorganic Equilibrium Systems Using Pitzer's Method to Determine Activity Coefficients. Abstract approved: Redacted for Privacy W. J. FredericK/ Redacted for Privacy r 6/ R. V. Mrazek Chemical equilibrium is a major factor inmany natural and industrial systems. The ability to predict ionic equilibrium is extremely valuable in industry in solvingproblems without the costs involved in full scale experiments. An example of an application in the pulp and paper industry is the predictionof non-process element solubility in process streams. A computer program, ISIS, was developed to estimate the solubil- ity of inorganic salts in aqueous, inorganic solutions. The model incorporates a two step Gibbs free energy minimization algorithmand Pitzer's method for ionic activity coefficient prediction. The program, ISIS, was tested on threecases. The test cases were: KC1 solubility in a NaC1 and water solution, NaC1 solubility in a MgC12 and water solution, and CaSO4 solubility ina NaC1 and water solution. The equilibrium predictions in the test caseswere very good, with the mean of the absolute value of the relativeerror ranging from 4.3% for the first case, 17% for the secondcase, to 3.8% in the third case. ISIS accurately predicted the solid phase amount and chemical composition in each of the test cases. A parametric study was conducted on the three testcases to examine the effects of the activity coefficient predictor, and of uncertainty in chemical potentials, and third virial coefficientsand mixing terms in the Pitzer activity coefficient prediction model. The effects of the activity coefficient predictorwere determined by comparing predictions made with the Pitzer modelversus an extended Debye-Huckel equation and an ideal solution assumption. The dif- ferences were great, with the other activity coefficientpredictors resulting in errors greater than five times theerror with the Pitzer activity coefficient model. The effect of the chemical potentials was large, especially in the trace species case, CaSO4 ina NaC1 and water solution. A relative change of less than a tenth ofa percent in the solid species chemical potential resulted inan increased error of ten times the original error. The effect of assuming the third virial coefficients and the mixingterms to be zero in the Pitzer activity coefficient predictor can be large. Errors up to 90% in the mean activity coefficient were found. It was concluded that the computer program ISIS couldpredict accurately the solubility of inorganic salts inaqueous, inorganic solutions. The accuracy of the prediction would be greatly affected by the accuracy of the chemical potentials and the availabilityof the third virial coefficient and mixing terms in the Pitzeractivity coefficient predictor. Future work is recommended in collecting accurate chemical potentials and Pitzer interaction terms to increase the database for ISIS and similar programs. The inclusion of ISIS into a steady state chemical process simulator is also recommended. A Computer Program to Solve Aqueous Inorganic Equilibrium Systems Using Pitzer's Method to Determine Activity Coefficients by Barry Kelly A THESIS Submitted to Oregon State University in partial fulfillment of the requirements for the degree of Master of Science Completed August 12, 1988 Commencement June 1989 APPROVED: Redacted for Privacy 13r/ofesdoVICh-e-mical Engineeringin charge of major Redacted for Privacy ProfeSZ;rof Forest Productsj4charge of major Redacted for Privacy HeadofOh,lical Engineering Department Redacted for Privacy Head of Forest/Products Department Redacted for Privacy e School1' Date thesis is presented: August 12. 1988 ACKNOWLEDGMENT The completion of this work represents the support and education I received from many individuals. I would like to acknowledge just a few of the many. Specifically, I wish to thank Dr. William J. Frederick and Dr. Robert V. Mrazek for their guidance and assistance throughout the research and writing of project. A special thanks goes to my family. The combination of support from my mother and father, constant harassment frommy brothers and sisters, and continual help from my fiance, JoAnne, createda driving force towards the completion of the project of incrediblemagnitude. I would also like to thank Dr. Biermann, Dr. Humphrey, Dr. Levien, and Dr. Nelson for having served on my committee, and Sadie Airth for her help with this manuscript. TABLE OF CONTENTS Ent INTRODUCTION 1 THEORY 5 A. Gibbs Free Energy Minimization 5 B. Activity Coefficient Calculation 9 RESULTS 21 A. Test Cases 22 B. Parametric Study 28 CONCLUSION 48 FUTURE WORK 50 BIBLIOGRAPHY 51 APPENDICES Appendix A: Calculation of the Electrostatic Mixing Terms 54 Appendix B: Average ISIS Run Times 56 Appendix C: ISIS User Manual 57 Appendix D: ISIS Operations Manual 86 Appendix E: ISIS Source Code 124 Appendix F: ISIS Database Files 199 Appendix G: PITZER Source Code 211 LIST OF FIGURES Figure Ent 1. Solubility of a slightly soluble salt, CaSO4, in NaC1 and water solutions (Silcock, 1979; Block and Waters, 1968; Marshel and Slusher, 1966) and calculated solubilities from several computer programs 3 2. Solubility of KC1 in NaC1 and water solutions at 25°C (Silcock, 1979) and calculated solubility from ISIS computer program 23 3. Solubility of NaCl in MgC12 and water solutions at 25°C (Silcock, 1979) and calculated solubility from ISIS computer program 24 4. Solubility of Nadi in MgC12 and water solutions at 25°C (Silcock, 1979) calculated solubility from ISIS computer program. Ionic strength, and absolute value of the relative deviations between the ISIS prediction and the experimental points versus MgC12 concentration 25 5. Solubility of CaSO4 in NaCl and water solutions at 25°C (Block and Waters, 1968; Marshel and Slusher, 1966; Silcock, 1979) and calculated solubility from ISIS computer program 27 6. Solubility of KC1 in NaCl and water solutions 15 25°C (Silcock, 1979) and calculated solubility from ISIS computer program using Pitzer, extended Dedye-Huckel, and ideal solution activity coefficient prediction models 29 7. Solubility of Nadi in MgC12 and water solutions at 25°C (Silcock, 1979) and calculated solubility from ISIS computer program using Pitzer, extended Dedye-Huckel, and ideal solution activity coefficient production models 30 8. Solubility of CaSO4 in NaC1 and water solutions at 25°C (Block and Waters, 1968; Marshel and Slusher, 1966; Silcock, 1979) and calculated solubility from ISIS computer program using Pitzer, extended Debye-Huckel, and ideal solution activity coefficient prediction models 32 9. Solubility of KC1 in NaC1 and water solutions at 25°C (Silcock, 1979) and calculated solubility from ISIS computer program using the Pitzer activity coefficient predictor and varying chemical potentials 34 Figure Page 10. Solubility of NaC1 in MgC12 and water solutions at 25°C (Silcock, 1979) and calculated solubility from ISIS computer program using the Pitzer activity coefficient predictor and varying chemical potentials 36 11. Solubility of CaSO4 in NaC1 and water solutions at 25°C (Block and Waters, 1968; Marshel and Slusher, 1966; Silcock, 1979) and calculated solubility from ISIS computer program using the Pitzer activity coefficient predictor and varying chemical potentials 39 12. Percent error in mean activity coefficient of NaC1 versus NaC1 ionic fraction for Case 1 45 13. Percent error in mean activity coefficient of KC1 versus KC1 ionic fraction for Case 1 45 14. Percent error in mean activity coefficient of NaC1 versus NaC1 ionic fraction for Case 2 46 15. Percent error in mean activity coefficient of MgC12 versus MgC12 ionic fraction for Case 2 46 16. Percent error in mean activity coefficient of NaC1 versus NaC1 ionic fraction for Case 3 47 17. Percent error in mean activity coefficient of CaSO4 versus CaSO4 ionic fraction for Case 3 47 LIST OF TABLES Table Page 1. Empirical Constants in the Pitzer Activity Coefficient Equations 16 2. Comparison of Chemical Potentials Used inCase1 35 3. Comparison of Chemical Potentials Used inCase2 37 4. Comparison of Chemical Potentials Used inCase3 37 NOMENCLATURE aw activity of water Actual actual solubility data from literature A, B, ai, bi constants in the extended Debye-Huckel equation (17) ai activity of species i Aki number of moles of element k in species i AS Pitzer-Debye-Huckel limiting slope of osmotic coefficient, Equation (31) b empirical parameter in Equation (30) B' MX derivative of BMXwith respect to ionic strength BOmx second virial coefficient term, Equation (22) MX second virial coefficient term, Equation (23) COmx third virial coefficient constant for salt MX dw density of water at 25 °C D dielectric constant of water at 25 °C e charge of electron Ek total moles of element k in system f(I) Debye-Huckel term, Equation (19) F sum of Debye-Huckel term and derivatives of the second virial coefficient terms with respect to ionic strength, Equation (30) G Gibbs free energy CEx excess Gibbs free energy AGj Gibbs free energy change of reaction j I ionic strength, Equation (18) Ifmx ionic fraction of salt MX, Equation (56) total number of all species ISIS value calculated by program ISIS Boltzmann constant kt total number of elements K search parameter m molality of salt mi molality of species i n number of data points compared ni number of moles of species i nw number of moles of water No Avogadro
Recommended publications
  • A Thermodynamic Model of Aqueous Electrolyte Solution Behavior And
    A thermodynamic model of aqueous electrolyte solution behavior and solid-liquid equilibrium in the Li-H-Na-K-Cl-OH-H2O system to very high concentrations (40 molal) and from 0 to 250 C Arnault Lassin, Christomir Christov, Laurent André, Mohamed Azaroual To cite this version: Arnault Lassin, Christomir Christov, Laurent André, Mohamed Azaroual. A thermodynamic model of aqueous electrolyte solution behavior and solid-liquid equilibrium in the Li-H-Na-K-Cl-OH-H2O system to very high concentrations (40 molal) and from 0 to 250 C. American journal of science, American Journal of Science, 2015, 815 (8), pp.204-256. 10.2475/03.2015.02. hal-01331967 HAL Id: hal-01331967 https://hal-brgm.archives-ouvertes.fr/hal-01331967 Submitted on 15 Jun 2016 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. A thermodynamic model of aqueous electrolyte solution behavior and solid- liquid equilibrium in the Li-H-Na-K-Cl-OH-H2O system to very high concentrations (40 molal) and from 0 to 250 °C Arnault Lassin1*, Christomir Christov2,3, Laurent André1, Mohamed Azaroual1 1 BRGM – 3 avenue C. Guillemin – Orléans, France ([email protected], [email protected], [email protected]) 2 GeoEco Consulting 2010, 13735 Paseo Cevera, San Diego, CA 92129, USA ([email protected]) 3 Institute of General and Inorganic Chemistry, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria ([email protected] ) * Corresponding author.
    [Show full text]
  • Downloaded Using the “Download” Table
    applied sciences Article A Review of Geochemical Modeling for the Performance Assessment of Radioactive Waste Disposal in a Subsurface System Suu-Yan Liang 1, Wen-Sheng Lin 2,*, Chan-Po Chen 3, Chen-Wuing Liu 4 and Chihhao Fan 5 1 Department of Civil Engineering, National Taiwan University, Taipei 10617, Taiwan; [email protected] 2 Hydrotech Research Institute, National Taiwan University, Taipei 10617, Taiwan 3 Institute of Environmental Engineering, National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan; [email protected] 4 Department of Water Resources, Taoyuan City Government, Taoyuan 33001, Taiwan; [email protected] 5 Department of Bioenvironmental Systems Engineering, National Taiwan University, Taipei 10617, Taiwan; [email protected] * Correspondence: [email protected] Abstract: Radionuclides are inorganic substances, and the solubility of inorganic substances is a major factor affecting the disposal of radioactive waste and the release of concentrations of radionuclides. The degree of solubility determines whether a nuclide source migrates to the far field of a radioactive waste disposal site. Therefore, the most effective method for retarding radionuclide migration is to reduce the radionuclide solubility in the aqueous geochemical environment of subsurface systems. In order to assess the performance of disposal facilities, thermodynamic data regarding nuclides in water–rock systems and minerals in geochemical environments are required; the results obtained from the analysis of these data can provide a strong scientific basis for maintaining safety performance Citation: Liang, S.-Y.; Lin, W.-S.; to support nuclear waste management. The pH, Eh and time ranges in the environments of disposal Chen, C.-P.; Liu, C.-W.; Fan, C.
    [Show full text]
  • Section 1: Identification
    RPP-RPT-50703 Rev.01A 8/23/2016 - 10:10 AM 1 of 77 Release Stamp DOCUMENT RELEASE AND CHANGE FORM Prepared For the U.S. Department of Energy, Assistant Secretary for Environmental Management By Washington River Protection Solutions, LLC., PO Box 850, Richland, WA 99352 Contractor For U.S. Department of Energy, Office of River Protection, under Contract DE-AC27-08RV14800 DATE: TRADEMARK DISCLAIMER: Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or any agency thereof or its contractors or subcontractors. Printed in the United States of America. Aug 23, 2016 1. Doc No: RPP-RPT-50703 Rev. 01A 2. Title: Development of a Thermodynamic Model for the Hanford Tank Waste Simulator (HTWOS) 3. Project Number: ☒ N/A 4. Design Verification Required: ☐ Yes ☒ No 5. USQ Number: ☒ N/A 6. PrHA Number Rev. ☒ N/A Clearance Review Restriction Type: public 7. Approvals Title Name Signature Date Checker CREE, LAURA H CREE, LAURA H 08/16/2016 Clearance Review RAYMER, JULIA R RAYMER, JULIA R 08/23/2016 Document Control Approval MANOR, TAMI MANOR, TAMI 08/23/2016 Originator BRITTON, MICHAEL D BRITTON, MICHAEL D 08/16/2016 Quality Assurance DELEON, SOSTEN O DELEON, SOSTEN O 08/16/2016 Responsible Manager CREE, LAURA H CREE, LAURA H 08/16/2016 8. Description of Change and Justification Updated the reduced chemical potential coefficient vlaues for Na2SO4 and Na2SO4·10H2O in Table A.1, as the original values were incorrect.
    [Show full text]
  • 3. Activity Coefficients of Aqueous Species 3.1. Introduction
    3. Activity Coefficients of Aqueous Species 3.1. Introduction The thermodynamic activities (ai) of aqueous solute species are usually defined on the basis of molalities. Thus, they can be described by the product of their molal concentrations (mi) and their γ molal activity coefficients ( i): γ ai = mi i (77) The thermodynamic activity of the water (aw) is always defined on a mole fraction basis. Thus, it can be described analogously by product of the mole fraction of water (xw) and its mole fraction λ activity coefficient ( w): λ aw = xw w (78) It is also possible to describe the thermodynamic activities of aqueous solutes on a mole fraction (x) basis. However, such mole fraction-based activities (ai ) are not the same as the more familiar (m) molality-based activities (ai ), as they are defined with respect to different choices of standard λ states. Mole fraction based activities and activity coefficients ( i), are occasionally applied to aqueous nonelectrolyte species, such as ethanol in water. In geochemistry, the aqueous solutions of interest almost always contain electrolytes, so mole-fraction based activities and activity co- efficients of solute species are little more than theoretical curiosities. In EQ3/6, only molality- (m) based activities and activity coefficients are used for such species, so ai always implies ai . Be- cause of the nature of molality, it is not possible to define the activity and activity coefficient of (x) water on a molal basis; thus, aw always means aw . Solution thermodynamics is a construct designed to approximate reality in terms of deviations from some defined ideal behavior.
    [Show full text]
  • Review of Uranyl Mineral Solubility Measurements
    Available online at www.sciencedirect.com J. Chem. Thermodynamics 40 (2008) 335–352 www.elsevier.com/locate/jct Review of uranyl mineral solubility measurements Drew Gorman-Lewis *, Peter C. Burns, Jeremy B. Fein University of Notre Dame, Department of Civil Engineering and Geological Sciences, 156 Fitzpatrick Hall, Notre Dame, IN 46556, United States Received 19 September 2007; received in revised form 10 December 2007; accepted 12 December 2007 Available online 23 December 2007 Abstract The solubility of uranyl minerals controls the transport and distribution of uranium in many oxidizing environments. Uranyl minerals form as secondary phases within uranium deposits, and they also represent important sinks for uranium and other radionuclides in nuclear waste repository settings and at sites of uranium groundwater contamination. Standard state Gibbs free energies of formation can be used to describe the solubility of uranyl minerals; therefore, models of the distribution and mobility of uranium in the environ- ment require accurate determination of the Gibbs free energies of formation for a wide range of relevant uranyl minerals. Despite dec- ades of study, the thermodynamic properties for many environmentally-important uranyl minerals are still not well constrained. In this review, we describe the necessary elements for rigorous solubility experiments that can be used to define Gibbs free energies of formation; we summarize published solubility data, point out difficulties in conducting uranyl mineral solubility experiments, and identify areas of future research necessary to construct an internally-consistent thermodynamic database for uranyl minerals. Ó 2007 Elsevier Ltd. All rights reserved. Keywords: Uranium; Uranyl minerals; Solubility; Gibbs free energy; Thermodynamics Contents 1.
    [Show full text]
  • Estimates of the Solubilities of Waste Element Radionuclides in Waste Isolation Pilot Plant Brines: a Report by the Expert Panel on the Source Term
    SEP f7 1996 SANDIA REPORT -'• <j . i, • i' •) V t=, LrJ SAND96-0098 • UC-721 Unlimited Release Printed May 1996 Estimates of the Solubilities of Waste Element Radionuclides in Waste Isolation Pilot Plant Brines: A Report by the Expert Panel on the Source Term David E. Hobart, Carol J. Bruton, Frank J. Millero, I-Ming Chou, Kathleen M. Trauth, D. Richard Anderson Prepared by Sandia National Laboratories Albuquerque, New Mexico 87185 and Livermore, California 94550 for the United States Department of Energy under Contract DE-AC04-94AL85000 , Approved for public release; distribution is unlimited. '• vi. >;'. •• HI :•.:!:!(':'•; ';' -6" 11 ii ''I 'illllii'i'lliilll'il <i H.lilKHl'Ill'H ,i| '' ',,.i, ';; '1| ilitii Hilili'ij'/'l iVlii li !•", il I. i.ll (iiiiiiliili'iiliim.iill iiidliiu../ C2R33 II SF2900Q(8-81) OF THIS DOCUMENT !§- Issued by Sandia National Laboratories, operated for the United States Department of Energy by Sandia Corporation. NOTICE: This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Govern- ment nor any agency thereof, nor any of their employees, nor any of their contractors, subcontractors, or their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, prod- uct, or process disclosed, or represents that its use would not infringe pri- vately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government, any agency thereof or any of their contractors or subcontractors.
    [Show full text]
  • Luginin's Therttgchanlcal Tabor-Atom Op Moscow St*Ta Unfverhtg International Synpaslun an Calortnetra and Ctieniual Thornadyimitii»
    </•/.//•".' DedicAted to the centenary of the Luginin's tHerttGcHanlcal tabor-atom oP Moscow St*ta UnfverHtg International Synpaslun an Calortnetra and Ctieniual ThornadyiMitii» June 33-38, 1991, noseow, USSR ABSTRACTS Dedicated to tke centenary ef the Luginiii's tHernoohanlcftl l«bar*tor* of ПОБВОИ Stata *nlvar»tt«i International SyMpasiuM MI C*lorinetra. «ml Cbaninal ИмгатАцмиИав June 83 - 3Bf 1991t Некоей• USSR ABSTRACTS Печатаетоя по решение оргкомитета международного симпо­ зиума по термохимии и термодинамике Оформление литературно-иэдатехьокое агенто!о Pi Злинина 1 CHEMICAL THERMODYNAMICS IN THE FUTURE DEVELOPMENT OF CHEMISTRY INCLUDING ENVIRONMENTAL PROBLEMS Leo Brewer Department of Chemistry, University of California, and Materials and Chemical Sciences Division, Lawrence Berkeley Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA. High-temperature science will be facing unusual challenges in the decade ahead. New technological advances are requiring new materials with unusual properties that will either be prepared by high-temperature techniques or will need to have long-term stability at high temperatures in various environments. One of the major driving forces for new materials arises from the increasing public concern about environmental pollution. Equipment using volatile fluids that can survive up to the stratosphere and destroy the ozone will have to be replaced. Processes that emit sulfur oxides win have to be modified to reduce sulfur emission to very low values. The efficiency of solar energy devices wilt have to be improved and nuclear power plants win have to be designed to make serious accidents extremely unlikely so that energy production by combustion to carbon dioxide is greatly reduced. Many other examples can be given of the need for new materials.
    [Show full text]
  • Estimation of the Pitzer Parameters for 1–1, 2–1, 3–1, 4–1, and 2–2 Single Electrolytes at 25 °C
    This is a repository copy of Estimation of the Pitzer Parameters for 1–1, 2–1, 3–1, 4–1, and 2–2 Single Electrolytes at 25 °C. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/102328/ Version: Accepted Version Article: Simoes, M.C., Hughes, K.J., Ingham, D.B. et al. (2 more authors) (2016) Estimation of the Pitzer Parameters for 1–1, 2–1, 3–1, 4–1, and 2–2 Single Electrolytes at 25 °C. Journal of Chemical & Engineering Data, 61 (7). pp. 2536-2554. ISSN 0021-9568 https://doi.org/10.1021/acs.jced.6b00236 This document is the Accepted Manuscript version of a Published Work that appeared in final form in Journal of Chemical Engineering & Data, copyright © American Chemical Society after peer review and technical editing by the publisher. To access the final edited and published work see http://dx.doi.org/10.1021/acs.jced.6b00236 Reuse Unless indicated otherwise, fulltext items are protected by copyright with all rights reserved. The copyright exception in section 29 of the Copyright, Designs and Patents Act 1988 allows the making of a single copy solely for the purpose of non-commercial research or private study within the limits of fair dealing. The publisher or other rights-holder may allow further reproduction and re-use of this version - refer to the White Rose Research Online record for this item. Where records identify the publisher as the copyright holder, users can verify any specific terms of use on the publisher’s website. Takedown If you consider content in White Rose Research Online to be in breach of UK law, please notify us by emailing [email protected] including the URL of the record and the reason for the withdrawal request.
    [Show full text]
  • 36. Aqua Module
    HSC 8 - Aqua November 19, 2014 Research Center, Pori / Petri Kobylin, Justin 14021-ORC-J 1 (49) Salminen, Matias Hultgren 36. Aqua Module 36.1. Introduction ............................................................................................................................. 3 36.2. Aqua module interface in HSC 8 ............................................................................................. 4 36.3. Input sheets ............................................................................................................................ 5 36.3.1. Binary system sheet........................................................................................................ 7 36.4. Output sheets ......................................................................................................................... 8 36.5. Multi-component examples ..................................................................................................... 9 36.6. Aqueous Equilibria ................................................................................................................ 11 36.7. General calculation principles for water systems .................................................................. 12 36.8. Excess enthalpy and heat capacities .................................................................................... 13 36.9. Summary of the equations used in HSC Aqua ...................................................................... 14 36.10. Conclusions .........................................................................................................................
    [Show full text]
  • Empirical Equations for Activity and Osmotic Coefficients N
    Eastern Illinois University The Keep Masters Theses Student Theses & Publications 1986 Empirical Equations for Activity and Osmotic Coefficients N. Ragunathan Eastern Illinois University This research is a product of the graduate program in Chemistry at Eastern Illinois University. Find out more about the program. Recommended Citation Ragunathan, N., "Empirical Equations for Activity and Osmotic Coefficients" (1986). Masters Theses. 2641. https://thekeep.eiu.edu/theses/2641 This is brought to you for free and open access by the Student Theses & Publications at The Keep. It has been accepted for inclusion in Masters Theses by an authorized administrator of The Keep. For more information, please contact [email protected]. THESIS REPRODUCTION CERTIFICATE TO: Graduate Degree Candidates who have written formal theses. SUBJECT: Permission to reproduc;e theses. The University Library is receiving a number of requests from other institutions asking permission to reproduce dissertations for inclusion in their library holdings. Although no copyright laws are involved, we feel that professional courtesy demands that permission be obtained from the author before we allow theses to be copied. Piease sign one of the following statements: Booth Library of Eastern Illin�is University has my permission to lend my thesis to a reputable college or university for the purpose of copying it for inclusion in that institution's library qr research holdings. Date Author I respectfully request Booth Librar_y of Eastern Illinois University not allow my thesis be reproduced because �- Date Author m Empirical Equations For Activity And Osmotic Coefficients. (TITLE) BY N.Ragunathan THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF Master cif Science, Department of Chemistry.
    [Show full text]
  • A Computer Program Incorporating Pitzer's Equations for Calculation of Geochemical Reactions in Brines
    A COMPUTER PROGRAM INCORPORATING PITZER'S EQUATIONS FOR CALCULATION OF GEOCHEMICAL REACTIONS IN BRINES By L.N. Plummer, D.L. Parkhurst, G.W. Fleming, and S.A. Dunkle U.S. GEOLOGICAL SURVEY Water-Resources Investigations Report 88-4153 Reston, Virginia 1988 DEPARTMENT OF THE INTERIOR DONALD PAUL MODEL, Secretary U.S. GEOLOGICAL SURVEY Dallas L. Peck, Director For additional information write to: Copies of this report can be purchased from: Chief, Branch of Regional Research, U.S. Geological Survey Northeastern Region Books and Open-File Reports Section U.S. Geological Survey Federal Center, Bldg. 810 432 National Center A\O. Box 25425 12201 Sunrise Valley Drive Denver, Colorado 80225 Reston, Virginia 22092 Additional copies of the software described in this report are available from: WATSTORE PROGRAM OFFICE U.S. Geological Survey 437 National Center Reston, VA 22092 -11- CONTENTS page Abstract 1 Introduction - ---- - - - -- - - --- -- - ---- -- 1 Pitzer equations ------- - - - - - - - - - -- -- -- - 2 Data base 5 Extensions of the data base 6 Internal consistency of the data base -- - - - - -- - 6 Literature sources of Pitzer interaction parameters - - - - - -- --- 11 Scale convention of activity coefficients - - - -- --- 11 Precautions and limitations ------- ----- -- - - --- - - -- - -- 14 Numerical method ------------------------ --- -- -- - ---- 15 Basic equations ---- -- - 15 Restrictions on the Newton-Raphson approach - - --- -- - - 17 Scaling the Newton-Raphson results ------- - - -- ___-- __ ___ 17 Description of input -- -- --
    [Show full text]
  • Activity Coefficients of a Simplified Seawater Electrolyte at Varying Salinity (5–40) and Temperature (0 and 25 °C) Using Monte Carlo Simulations
    Gothenburg University Publications Activity coefficients of a simplified seawater electrolyte at varying salinity (5–40) and temperature (0 and 25 °C) using Monte Carlo simulations This is an author produced version of a paper published in: Marine Chemistry (ISSN: 0304-4203) Citation for the published paper: Ulfsbo, A. ; Abbas, Z. ; Turner, D. (2015) "Activity coefficients of a simplified seawater electrolyte at varying salinity (5–40) and temperature (0 and 25 °C) using Monte Carlo simulations". Marine Chemistry, vol. 171 pp. 78-86. http://dx.doi.org/10.1016/j.marchem.2015.02.006 Downloaded from: http://gup.ub.gu.se/publication/214853 Notice: This paper has been peer reviewed but does not include the final publisher proof- corrections or pagination. When citing this work, please refer to the original publication. (article starts on next page) Activity coefficients of a simplified seawater electrolyte at varying salinity (5–40) and temperature (0 and 25°C) using Monte Carlo simulations Adam Ulfsbo*, Zareen Abbas, David R. Turner Department of Chemistry and Molecular Biology, University of Gothenburg, SE-412 96 Gothenburg, Sweden *Corresponding author: [email protected] (Adam Ulfsbo) Phone: +46 31 786 9053 1 Abstract + 2+ 2+ - 2- Mean salt activity coefficients of a simplified seawater electrolyte (Na , Mg , Ca , Cl , SO4 , - 2- HCO3 , CO3 ) at varying salinity (5-40) and temperature (0-25 °C) were estimated by Monte Carlo (MC) simulations, and compared with Pitzer calculations. The MC simulations used experimentally determined dielectric constants of water at different temperatures, and optimal agreement with the experimental data and the Pitzer calculations was achieved by adjusting the ionic radii.
    [Show full text]