Return to Kaj Thomsen

Publication list of Kaj Thomsen

Contributions to books

Chemical Thermodynamics for Industry, (RSC books, 2004), Edited by T.M. Letcher
Chapter on "Thermodynamics of Electrolyte Systems of Industry" by
Kaj Thomsen


Kemiske Enhedsoperationer 5. udgave (2004, Chemical Unit Operations, Language: Danish)
Authors: Karsten H. Clement, Peder Fangel, Anker Jensen,
Kaj Thomsen

Rettelsesliste


Electrolyte Solutions: Thermodynamics, Crystallization, Separation methods
Teaching material for course on Aqueous Electrolyte Thermodynamics.
Author: Kaj Thomsen


Papers, reports, and thesis

Leila Faramarzi, Georgios M. Kontogeorgis, Kaj Thomsen, Erling H. Stenby, Extended UNIQUAC model for thermodynamic modeling of CO2 absorption in aqueous alkanolamine solutions, Fluid Phase Equilibria, 282(2009)121–132

ABSTRACT: The extended UNIQUAC model [Thomsen, Rasmussen, Chem. Eng. Sci. 54 (1999) 1787-1802] was applied to the thermodynamic representation of carbon dioxide absorption in aqueous monoethanolamine (MEA), methyldiethanolamine (MDEA) and varied strength mixtures of the two alkanolamines (MEA-MDEA). For these systems, altogether 13 interaction model parameters are adjusted. Out of these parameters, 11 are temperature dependent. All the essential parameters of the model are simultaneously regressed to a collective set of data on the single MEA and MDEA systems. Different types of data are used for modeling and they cover a very wide range of conditions. Vapor-liquid equilibrium (VLE) data for the aqueous alkanolamine systems containing CO2 in the pressure range of 3-13000 kPa and temperatures of 25-200°C are used. The model is also regressed with the VLE and freezing point depression data of the binary aqueous alkanolamine systems (MEA-water and MDEA-water). The two just mentioned types of data cover the full concentration range of alkanolamines from extremely dilute to almost pure. The experimental freezing point depression data down to the temperature of -20°C are used. Experimental excess enthalpy (HE) data of the binary MEA-water and MDEA-water systems at 25, 40, 65 and 69°C are used as well. In order to enhance the calculation of the infinite dilution activity coefficients of MEA and MDEA, the pure alkanolamines vapor pressure data in a relevant temperature range (up to almost 230°C) are included in the parameter estimation process. The previously unavailable standard state properties of the alkanolamine ions appearing in this work i.e. MEA protonate, MEA carbamate and MDEA protonate are determined. The concentration of the species in both MEA and MDEA solutions containing CO2 are predicted by the model and in the case of MEA compared to NMR spectroscopic data. Using only one set of parameters for correlation of different thermodynamic properties, the model has represented the experimental data with good precision.



Philip L. Fosbøl; Kaj Thomsen; Erling H. Stenby, Modeling of the Mixed Solvent Electrolyte System CO2-Na2CO3-NaHCO3-Monoethylene Glycol-Water, Industrial & Engineering Chemistry Research, 48(2009)4565-4578


ABSTRACT: The Extended UNIQUAC electrolyte activity coefficient model has been correlated to 751 experimental solid-liquid equilibrium (SLE), vapor-liquid equilibrium (VLE), and excess enthalpy data for the mixed solvent CO2-NaHCO3-Na2CO3-Mono ethylene glycol(MEG)-H2O electrolyte system. The model was validated by predicting the excess heat capacity. The model is consistent and one set of binary interaction parameters are used for calculating all the properties between 50 and 90 °C. The model is compared to experimental data of infinite dilution activity coefficient measurements of MEG and may be used for calculating activities, solubility, salt hydrate properties, pH, and CO2 solubility in the NaHCO3-Na2CO3-MEG-H2O system. A simple density model of NaHCO3-Na2CO3-NaCl-MEG-water is given.



Philip Loldrup Fosbøl, Kaj Thomsen; Erling Halfdan Stenby, Solubility Measurements in the Mixed Solvent Electrolyte System Na2CO3-NaHCO3-Monoethylene Glycol-Water, Industrial & Engineering Chemistry Research, vol: 48(4), p. 2218-2228 (2009).


ABSTRACT: 212 new data points for solubility and density have been measured in the mixed solvent electrolytic CO2-NaHCO3-Na2CO3-Mono ethylene glycol (MEG)-water system. Measurements were conducted at 2 to 60 °C at atmospheric pressure. An overview of methods available for determining the carbon dioxide and sodium content is given and the reverse Schreinemakers method was chosen. The method gives the amounts of all the species in the liquid. The solutions have a tendency to form meta-stable mixtures from which sodium carbonate salt hydrates may precipitate. 30 meta-stable data points are reported.


Philip L. Fosbøl, Kaj Thomsen, and Erling H. Stenby, Reverse Schreinemakers method for experimental analysis of mixed-solvent electrolyte systems, J. Solution Chem. 38(2009)1–14.


ABSTRACT: A method based on Schreinemakers tie-line theory of 1893 is derived for determining the composition and phase amounts in solubility experiment for multi-solvent electrolyte systems. The method uses the lever rule in reverse compared to Schreinemakers wet residue method and is therefore called the reverse Schreinemakers(RS) method. The method is based on simple mass balance principles similar to the wet residues method. It allows for accurate determination of mixed-solvent phase composition even though solvent may precipitate as complexes of solvent and salt. Discrepancy between determining composition of salt mixtures by pH titration is discussed and the derived method significantly improves the obtained result from titration. Furthermore the method reduces the required experimental work for analysis of phase composition. The method is applicable to multi-solvent systems and may be used for determination of solid composition, similar to Schreinemakers original rest method. An example calculation in the Na2CO3-NaHCO3-MEG-H2O system is presented.


Lars Jensen, Kaj Thomsen, and Nicolas von Solms, “Propane Hydrate Nucleation: Experimental Investigation and Correlation” (Chemical Engineering Science 63(2008)3069-3080)

ABSTRACT: In this work the nucleation kinetics of propane gas hydrate has been investigated experimentally using a stirred batch reactor. The experiments have been performed isothermally recording the pressure as a function of time. Experiments were conducted at different stirring rates, but in the same supersaturation region. The experiments showed that the gas dissolution rate rather than the induction time of propane hydrate is influenced by a change in agitation. This was especially valid at high stirring rates when the water surface was severely disturbed. Addition of polyvinylpyrrolidone to the aqueous phase was found to reduce the gas dissolution rate slightly. However the induction times were prolonged quite substantially upon addition of polyvinylpyrrolidone. The induction time data were correlated using a newly developed induction time model based on crystallization theory also capable of taking into account the presence of additives. In most cases reasonable agreement between the data and the model could be obtained. The results revealed that especially the effective surface energy between propane hydrate and water is likely to change when the stirring rate varies from very high to low. The prolongation of induction times according to the model is likely to be due to a change in the nuclei-substrate contact angle.


Philip L. Fosbøl, Kaj Thomsen, and Erling H. Stenby , A Review and Recommended Thermodynamic Properties of FeCO3, (Accepted for publication in Corrosion Engineering Science and Technology, 2008)

ABSTRACT: An extensive review of entropy, enthalpy of formation and Gibbs energy of formation, heat capacity, aqueous solubility and solubility constant of FeCO3 is given. A consistent set of thermodynamic properties for FeCO3 and relevant aqeous species is selected and recommended for use. Speciation schemes for aqeous FeCO3 are reviewed and evaluated. Issues related to supersaturation of FeCO3 are discussed. Works on the thermal decomposition of FeCO3 are presented and an overview of measured solubility and synthesis of FeCO3 is given.


Zheng Guo, Bena-Marie Lue, Kaj Thomsen, Anne Boye Strunge Meyer, and Xuebing Xu, Predictions of Flavonoid Solubility in Ionic Liquids by COSMO-RS: Experimental Verification, Structural Elucidation, and Solvation Characterization,  (Green Chemistry,  9(2007)1362-1373))

ABSTRACT: Predictions of the solubility of flavonoids in a large variety of ionic liquids (ILs) with over 1800 available structures were examined based on COSMO-RS computation. The results shows that the solubilities of flavonoids are strongly anion-dependent; and the ILs could be classified into 3 groups according to the values of esculin solubility. The predictions were experimentally verified by the measurement of the solubilities of esculin and rutin in 12 ILs with varying anion and cation parts. It is shown that predicted and experimental results generally have a good agreement, indicating the correctness of the physics of COSMO-RS as an experimentally independent approach, and the application potential to preselect favourable structures from a large pool of available ILs. Importantly, this work first systemically demonstrated that COSMO-RS derived parameters, electrostatic misfit, H-bonding, and van der Waals interaction energy, are capable of and effectively characterize the complicated multiple interactions in IL system. Force field analysis shows that H-bonding interaction is the most dominant interaction for ILs (followed by electrostatic misfit and van der Waals interactions) to determine the solubility of flavonoids, and anionic part has greater effect on the overall H-bonding capability of the IL. Based on solvation interaction energy with flavonoids and hydrogen bond basicity derived from COSMO-RS, ILs were scaled and categorized qualitatively and quantitatively, which may be of general value to better understand the solvation behaviors of ILs. Conversely, the improved understanding and knowledge of the correct association of specific interactions with defined groups of ILs may help to mix and match these moieties into an optimal structure for a particular application. Therefore, the data obtained in this study may be important for the tailoring of the desired structures of ILs used as the media for efficient enzymatic esterification of flavonoids.


Yi Lin, Kaj Thomsen and Jean-Charles de Hemptinne, Multi Component Equations of State for Electrolytes”,  (AIChE Journal, 53(4)(2007)989-1005)

ABSTRACT: Four equations of state have been implemented and evaluated for multi-component electrolyte solutions at different temperatures and pressures. The equations contain terms accounting for short-range and long-range interactions in electrolyte solutions. Short range interactions are described by one of the three equations of state, Peng-Robinson, Soave-Redlich-Kwong, or Cubic-Plus-Association (CPA). Long range interactions are described by either the simplified mean spherical approximation (MSA) solution of the Ornstein–Zernicke equation or the simplified Debye-Hückel term. An optional Born term is added to these electrostatic terms. The resulting electrolyte equations of state were tested by determining the optimal model parameters for the multi-component test system consisting of H2O, Na+, H+, Ca2+, Cl-, OH-, SO42-. In order to describe the thermodynamics of this multi-component system, ion specific parameters were determined. The parameters in the equations of state were fitted to experimental data consisting of apparent molar volumes, osmotic coefficients, mean ionic activity coefficients, and solid-liquid equilibrium data. The results of the parameter fitting are presented. The ability of the equations of state to reproduce the experimental data is demonstrated. The performance of the equations of state for multi-component systems is compared and analyzed in view of the various short range and long range terms employed.


Kaj Thomsen, Jørgen Peter Jensen, Peter Simonsen, Bo Sander, Reuse of Alkali from Fly Ash from Biomass Combustion,  Report written in Danish language on research project sponsored by PSO (Danish Power Plants).

ABSTRACT: Experiments to leach potassium salts from fly ash from straw combustion were performed. The fly ash was produced in the bio mass boiler at the Avedøre power plant in south western Copenhagen, Denmark. The fly ash contained approximately 90 mass percent water soluble material. When the fly ash was dissolved at low pH, a slightly higher solubility was found. 100 gram fly ash consisted typically of 9 gram insoluble material, 9 gram calcium phosphate, 29 gram potassium sulfate and 53 gram of potassium chloride. In addition, 100 gram of fly ash contained approximately 1 mg of cadmium, corresponding to a concentration of cadmium of 10 ppm in the fly ash. Fly ash from the bio mass boiler at the Avedøre power plant apparently has a significantly larger content of potassium salts than fly ash from other boilers.
The Extended UNIQUAC thermodynamic model was used for calculating relevant phase diagrams and calculations of the necessary amount of water required for dissolving all the KCl and all KCl + K2SO4 of the fly ash. This theoretical minimum amount of water was calculated at a range of temperatures between 10 and 100°C. The amount of water required at 100°C was less than half of that required at 10°C. Experiments were performed in order to find a feasible method for separating the potassium salts of the fly ash from the ash residue and especially from the soluble cadmium salts found in the fly ash.
Experiments with counter current leaching of fly ash in a fluid bed gave unsatisfactory results. Apparently there was a lack of contact between the wash water and the ash. In addition, sedimentation was very slow resulting in an incomplete separation of wash water and ash residue. Experiments with ion exchange by adding CaCl2 to the wash water and successive precipitation of gypsum or anhydrite gave unsatisfactory results. Process simulation had shown that by this method the necessary amount of washing water could be decreased. This is due to the higher solubility of potassium chloride compared to potassium sulfate. By using this reduced amount of washing water with CaCl2 a viscous, muddy slurry was obtained. It was not possible to separate this slurry into a liquid phase and a solid phase.
It was found that a good separation was obtained by intense stirring followed by centrifuging or by filtering the slurry in a filter press. The lab experiments were performed as a counter current leaching process in four stages. Depending on the applied amount of water relative to the theoretical minimum amount of water a complete separation could be obtained in two or more stages. More stages are required if less water is used. If the leaching is performed in less stages, a larger amount of water is required.
If the washing process was performed at pH 1, the phosphate fraction could successively be obtained as a fine, white powder of almost pure calcium phosphate by regulating pH to 12. If the washing process was conducted at a pH between 1 and 5 and pH was successively adjusted to 12, cadmium could not be detected in the brine by current analysis methods. All cadmium was precipitated as phosphate.


Kaj Thomsen, Duc Thoung Vu, Mette Stenby, Jørgen Peter Jensen, Peter Simonsen and Bo Sander, Leaching of Nutrient Salts from Fly Ash from Biomass Combustion, (Proceedings from 14th European Conference and Technology Exhibition on Biomass for Energy, Industry and Climate Protection, October 2005, p. 1273-1276)

ABSTRACT: Methods to selectively leach nutrient salts from fly ash, while leaving cadmium un-dissolved were studied. Temperature, pH, water to fly ash ratio are all expected to influence the kinetics and the equilibrium boundaries for this process. Three different leaching methods were investigated. The first method was a counter current moving bed process in four stages. The ash was kept in filter bags and leached with water that was introduced into the bags at 40-50°C. In the second method, fly ash and water was brought into contact in a partially fluidized bed. The third method was a counter current moving bed process with agitation/centrifugation. It was found that a satisfactory leaching of the nutrient salts could be achieved with the third method using only two or three stages, depending on the water to fly ash ratio. It is an advantage to perform the process at temperatures above 50°C as the amount of water used for leaching at this temperature can be reduced compared to the amount of water needed at 25°C.


Ada Villafáfila García, Kaj Thomsen, and Erling H. Stenby, Prediction of Mineral Scale Formation in Geothermal and Oilfield operations using the Extended UNIQUAC Model. Part II: Carbonate Scaling Minerals,  (Geothermics 35(2006)239-284)

ABSTRACT: Two additional parameters to account for the pressure dependency of solubility are added to the Extended UNIQUAC model presented by Thomsen and Rasmussen (1999). The improved model has been used for correlation and prediction of vapor-liquid-solid equilibrium for different carbonate systems (CaCO3, MgCO3, BaCO3 and SrCO3) causing scale problems. The solubility of NaCl and CO2 in pure water, and the solubility of CO2 in solutions of different salts (NaCl and Na2SO4) have also been correlated. The temperature and pressure range covered is from 0 to 250°C and from 1 to 1000 bar, respectively. The results show that the Extended UNIQUAC model, with the proposed pressure parameters, is able to represent binary (NaCl-H2O, MCO3-H2O and CO2-H2O), ternary (MCO3-CO2-H2O, CO2-NaCl-H2O and CO2-Na2SO4-H2O) and quaternary (CO2-NaCl-Na2SO4-H2O) solubility data within the experimental accuracy in the temperature range from 0 to 250°C, and the pressure range from 1 to 1000 bar. M stands for Ca2+, Mg2+, Ba2+ and Sr2+


Søren Gregers Christensen and Kaj Thomsen, Representation of volumetric data of electrolyte solutions at varying concentrations and temperatures,  (Internal report)

ABSTRACT: A modification of the Masson equation combined with Young’s rule based on ion specific parameters has been applied to volumetric data for mixtures of (H+, Na+, K+, NH4+, Ca++, Mg++) (Cl-,NO3-,SO4--). The parameters have been regressed from data in the temperature range 0 – 100°C and concentration range 0 – 11.8 mol/kg but are shown to be valid up to saturation. The model only requires 5 parameters per ion in the entire range of concentration and temperature. The model is easily applied to multi component mixtures, and it is shown that the relative errors of the predicted results in both ternary and quaternary systems are well within the experimental accuracy of the data.


Kaj Thomsen, Modeling Electrolyte Solutions with the extended universal quasichemical (UNIQUAC) Model, (Presented at the 11th ISSP in Aveiro, Portugal, 2004) (Journal of Pure and Applied Chemistry, 77(2005)531-542, issue 3 )

ABSTRACT: The extended universal quasichemical (UNIQUAC) model is a thermodynamic model for solutions containing electrolytes and non-electrolytes. The model is a gibbs excess function consisting of a Debye-Hückel term and a standard UNIQUAC term. The model only requires binary, ion specific interaction parameters. A unique choice of standard states makes the model able to reproduce solid-liquid, vapor-liquid, and liquid-liquid phase equilibria as well as thermal properties of electrolyte solutions using one set of parameters.


Søren Gregers Christensen and Kaj Thomsen, Experimental measurement and modeling of the distribution of solvent and ions between an aqueous phase and an ion exchange resin, (Fluid Phase Equilibria, 228-229(2005)247-260).

ABSTRACT: The distribution of solutes and solvent between an aqueous solution of salt and an ion exchange resin has been measured at ambient temperature. The experiments have been performed for aqueous solutions of KNO3, KCl, Ca(NO3)2 and CaCl2 in the concentration range of 0-3N. The absorption has been measured for 3 gel type and 3 macroreticular resins with a degree of crosslinking varying from 10.5 to 18.5%. The experimental results have been modeled with the Extended UNIQUAC model combined with an elastic term taking the elastic properties of the resin structure into account. The model shows very good predictions with varying degree of crosslinking, and the deviations between model results and experimental data are all within the experimental error.


Ada Villafáfila García, Kaj Thomsen, Prediction of Mineral Scale Formation in Geothermal and Oilfield Operations using the Extended UNIQUAC Model. Part I: Sulphate Scaling Minerals, and Erling H. Stenby (Geothermics, 34(2005)61-97)

Pressure parameters are added to the Extended UNIQUAC model presented by Thomsen and Rasmussen (1999). The improved model has been used for correlation and prediction of solid-liquid equilibrium (SLE) of scaling minerals (CaSO4, CaSO4·2H2O, BaSO4 and SrSO4) at temperatures up to 300°C and pressures up to 1000 bar. The results show that the Extended UNIQUAC model, with the proposed pressure parameters is able to represent binary, ternary and quaternary solubility data within the experimental accuracy in the temperature range from -20 to 300°C, and the pressure range from 1 to 1000 bar.


Kaj Thomsen, Maria Iliuta, and Peter Rasmussen "Extended UNIQUAC model for correlation and prediction of vapor-liquid-liquid-solid equilibria in aqueous salt systems containing non-electrolytes. Part B. Alcohol (Ethanol, Propanols, Butanols) - water - salt systems". (Chemical Engineering Science 59(2004)3631-3647, issue 17)

ABSTRACT:  The Extended UNIQUAC model is an electrolyte model formed by combining the original UNIQUAC model, the Debye-Hückel law, and the Soave-Redlich-Kwong equation of state. The model only requires binary, temperature dependent interaction parameters. It has previously been used to describe the excess Gibbs energy for aqueous electrolyte mixtures and aqueous electrolyte systems containing methanol. It has been found to be an adequate model for representing solid-liquid-vapor equilibrium and thermal property data for strongly non-ideal systems. In this work the model is extended to aqueous salt systems containing higher alcohols. The calculations are based on an extensive database consisting of salt solubility data, vapor liquid equilibrium data, and liquid-liquid equilibrium data for solvent mixtures and for mixed solvent-electrolyte systems.
The application of this model to represent the vapor-liquid-liquid-solid equilibria in aqueous systems containing various non-electrolytes (ethanol, 1-propanol, 2-propanol, 1-butanol, 2-butanol, 2-methyl 1-propanol, 2-methyl 2-propanol) and various ions (Na+, K+, NH4+, Cl-, NO3-, SO42-, SO32-, HSO3-, CO32-,and HCO3-) shows the capability of the model to accurately represent the phase behavior of these kinds of systems.


Kaj Thomsen,  "Thermodynamics of Electrolyte Systems of Industry", Published in Chemical Thermodynamics for Industry, (RSC books, 2004), Edited by T.M. Letcher


Søren Gregers Christensen and Kaj Thomsen, “Modeling of Vapor-Liquid-Solid Equilibria in Acidic Aqueous Solutions”  (Ind. & Eng. Chem. Res. 42(2003)4260-4268, issue 18)

ABSTRACT: The phase behavior (vapor - liquid equilibria (VLE) and solid – liquid equilibria (SLE)) and thermal properties of aqueous solutions of ions like (K+, Na+, NH4+, Ca2+, Cl-) in the presence of phosphoric acid (H3PO4, H2PO4-, HPO42- ) and nitric acid (HNO3, NO3-) are described by means of the Extended UNIQUAC model. Model parameters are evaluated on the basis of more than 2000 experimental data points. There is good agreement between calculated and experimental data points. The model parameters are valid in the temperature range from -18 - +122°C and in the concentration range up to 12 molal for the acids HNO3 and H3PO4.


Søren Gregers Christensen and Kaj Thomsen, "Production of fertilizer salts"  (Dansk Kemi, 83(2) (2002)18-19)


Maria C. Iliuta, Kaj Thomsen and Peter Rasmussen "Modeling of heavy metal salt solubility using the Extended UNIQUAC model" (AIChE Journal, 48(11)(2002)2664-2689)

ABSTRACT: Solid-liquid equilibria in complex systems involving a heavy metal cation (Mn2+, Fe2+, Co2+, Ni2+, Cu2+, or Zn2+) and one or more ions for which Extended UNIQUAC parameters have been published previously are modeled using the Extended UNIQUAC model. Model parameters are determined on the basis of a data bank with more than 4000 experimental data points for binary and ternary systems. The parameters are generally valid in the temperature range from the cryohydratic point to the boiling point of the respective solutions.


Raphaël Huyghe, Peter Rasmussen, and Kaj Thomsen, "Solid-Liquid Equilibria for the Binary Mixtures 1,4-Xylene + Ethylbenzene and 1,4-Xylene + Toluene". (Chemical Engineering Communications 191(8)(2004)1017-1023)

ABSTRACT: Solid-liquid equilibrium (SLE) data for the binary mixtures 1,4-xylene + ethylbenzene, and 1,4-xylene + toluene have been measured using differential scanning calorimetry (DSC) in the temperature range from 133.15 K to 293.15 K.


Morten Mejlholm, Kaj Thomsen, Peter Rasmussen, Jørgen Vergod, Freddy Knudsen, Hugo Høyer: "SODIUM CHLORIDE DIHYDRATE - A POTENTIAL CAUSE OF SLIPPERY ACCIDENTS" (Presented at the XIth PIARC International Winter Road Congress, Sapporo, Japan, January 28-31, 2002. Proceedings of the XIth PIARC International Winter Road Congress, Sapporo, Japan (2002))

ABSTRACT: From a thermodynamic point of view, it can be expected that sodium chloride dihydrate (hydrohalite, NaCl·2H2O) will form on winter roads under certain conditions at temperatures below 0.1°C. In order to elucidate whether or not the formation of hydrohalite on the pavement can explain the phenomenon of ice appearing to be resistant to road salt, a comparative study has been made on a number of different surfaces measuring the friction index. The friction measurements were performed with a Portable Skid-Resistance Tester. Discontinuous surfaces consisting of small islands of hydrohalite was classified as potentially slippery surfaces. It is therefore possible that the formation of hydrohalite contributes to accidents on slippery roads.


Selva Pereda, Kaj Thomsen and Peter Rasmussen, "Vapor-Liquid-Solid Equilibria of Sulfur Dioxide in Aqueous Electrolyte Solutions" Chemical Engineering Science 55(2000)2663-2671.

ABSTRACT: The Extended UNIQUAC model for electrolyte systems, combined with the Soave-Redlich-Kwong equation of state is used to describe the complex vapor-liquid-solid equilibria of sulfur dioxide in electrolyte solutions. Model parameters based on 1500 experimental data points are presented. The parameters are applicable in the temperature range 0 to 110 °C, concentrations up to saturation and pressures up to 30 bar. This validity range corresponds to the experimental data used for the evaluation of parameters.


Maria Iliuta, Kaj Thomsen and Peter Rasmussen, "Extended UNIQUAC model for correlation and prediction of vapour-liquid-solid equilibria in aqueous salt systems containing non-electrolytes . Part A. Methanol - water - salt systems " Chemical Engineering Science, 55(2000)2673-2686

ABSTRACT: The Extended UNIQUAC model has previously been used to describe the excess Gibbs energy for aqueous electrolyte mixtures. It is an electrolyte model formed by combining the original UNIQUAC model, the Debye-Hückel law and the Soave-Redlich-Kwong equation of state. In this work the model is extended to aqueous salt systems containing non-electrolytes in order to demonstrate its ability in representing solid-liquid-vapour (SLV) equilibrium and thermal property data for these strongly non-ideal systems. The model requires only pure component and binary temperature dependent interaction parameters. The calculations are based on an extensive database consisting of salt solubility data in pure and mixed solvents, VLE data for solvent mixtures and mixed solvent - electrolyte systems and thermal properties for mixed solvent solutions. Application of the model to the methanol - water system in the presence of several ions (Na+, K+, NH4+, Cl-, NO3-, SO42-, CO32-, and HCO3-) shows that the Extended UNIQUAC model is able to give an accurate description of VLE and SLE in ternary and quaternary mixtures, using the same set of binary interaction parameters. The capability of the model to predict accurately the phase behaviour of methanol - water - three salts systems is illustrated.


K. Thomsen and P. Rasmussen, "Thermodynamic Model for the Ammonia-Water System" ("Steam Water and Hydrothermal Systems: Physics and Chemistry Meeting the Needs of Industry" Proceedings of the 13th International Conference on the Properties of Water and Steam, Editors: P.G. Hill, P. Tremaine, D. Irish, and P.V. Balakrishnan, NRC Press, Ottawa, 2000, p. 118-125)

ABSTRACT: The ammonia-water system is described by the Extended UNIQUAC model, which is an electrolyte model, formed by combining the original UNIQUAC model, the Debye-Hückel law and the Soave-Redlich-Kwong equation of state. The model is limited to temperatures below the critical temperature of ammonia. Vapor-liquid equilibria are described within the experimental accuracy. The accuracy of enthalpy calculations is better than ± 100 J mol-1, and heat capacity calculations deviate less than ± 1.0% from experimental data. The accurate description of the thermal properties is achieved by taking speciation equilibria into account. Model parameters valid in the temperature range 0 - 130°C are given.


Kaj Thomsen and Peter Rasmussen: "Modeling of vapor - liquid - solid Equilibria in gas - aqueous electrolyte systems" Chemical Engineering Science Vol. 54(1999)1787-1802

Abstract: A thermodynamic model for the description of vapor-liquid-solid equilibria is introduced. This model is a combination of the extended UNIQUAC model for electrolytes and the Soave-Redlich-Kwong cubic equation of state. The model has been applied to aqueous systems containing ammonia and/or carbon dioxide along with various salts. Model parameters valid in the temperature range 0 - 110 °C, the pressure range from 0 - 100 bar, and the concentration range up to approximately 80 molal ammonia are given. The model parameters were evaluated on the basis of more than 7000 experimental data points.


Kaj Thomsen, Peter Rasmussen, and Rafiqul Gani: "Simulation and optimization of fractional crystallization processes". Chemical Engineering Science, Vol. 53(1998)1551-1564.

Abstract: A general method for the calculation of various types of phase diagrams for aqueous electrolyte mixtures is outlined. It is shown how the thermodynamic equilibrium precipitation process can be used to satisfy the operational needs of industrial crystallizer/centrifuge units. Examples of simula tion and optimization of fractional crystallization processes are shown. In one of these examples a process with multiple steady states is analyzed. The thermodynamic model applied for describing the highly non-ideal aqueous electrolyte systems is the Extended UNIQUAC model.


Kaj Thomsen: "Aqueous electrolytes: model parameters and process simulation". Ph.D. thesis, (1997). Department of Chemical Engineering, Technical University of Denmark.  Download thesis in pdf format

Abstract: This thesis deals with aqueous electrolyte mixtures. The Extended UNIQUAC model is being used to describe the excess Gibbs energy of such solutions. Extended UNIQUAC parameters for the twelve ions Na+, K+, NH4+, H+, Cl-, NO3-, SO42-, HSO4-, OH-, CO32-, HCO3-, and S2O82- are estimated. A computer program including a steady state process simulator for the design, simulation, and optimization of fractional crystallization processes is presented.


Kaj Thomsen, Peter Rasmussen, and Rafiqul Gani: "Correlation and Prediction of Thermal Properties and Phase Behaviour for a Class of Aqueous Electrolyte Systems". Chemical Engineering Science, Vol. 51(1996)3675-3683.

Abstract: An extended UNIQUAC model is used to describe phase behaviour (VLE, SLE) and thermal properties (heat of mixing, heat capacity) for aqueous solutions containing ions like (Na+, K+, H+) (Cl-, NO3-, SO42-, OH-, CO32-, HCO3-). A linear temperature dependence of the binary interaction parameters allows good agreement with experimental data in the temperature range 0-110°C.


Kaj Thomsen, Rafiqul Gani, and Peter Rasmussen: "Synthesis and analysis of processes with electrolyte mixtures". Computers and Chemical Engineering 19S(1995)S27-S32

Abstract: A computer aided system for synthesis, design and simulation of crystallization and fractional crystallization processes with electrolyte mixtures is presented. The synthesis methodology is based on the use of computed solubility diagrams for the corresponding electrolyte systems. For a specified crystallizer and product(s), the process flowsheet along with conditions of operation is determined while for a specified feed mixture and product(s), the process flowsheet together with the type of crystallizer is determined. In order to verify the "determined" flowsheet an option to perform steady state simulation is also provided. Examples highlighting the various features of the computer aided system are presented.