Thermodynamics framework

The bubble and dew-point temperature calculations have been implemented by the FORTRAN IV subroutine BUDET and the pressure calculations by subroutine BUDEP, which are described and listed in Appendix F. These subroutines calculate the unknown temperature or pressure, given feed composition and the fixed pressure or temperature. They provide for input of initial estimates of the temperature or pressure sought, but converge quickly from any estimates within the range of validity of the thermodynamic framework. Standard initial estimates are provided by the subroutines. 

Flash calculations for these mixtures usually require four to eight iterations. Cases 5 and 6 in Table 1 have feeds of this type, including noncondensable components in Case 6. Within the limits of the thermodynamic framework used here, no case has been encountered where FLASH has required more than 12 iterations for satisfactory convergence. 

The treatments that are concerned in more detail with the nature of the adsorbed layer make use of the general thermodynamic framework of the derivation of the Gibbs equation (Section III-5B) but differ in the handling of the electrochemical potential and the surface excess of the ionic species [114-117]. The derivation given here is after that of Grahame and Whitney [117]. Equation III-76 gives the combined first- and second-law statements for the surface excess quantities... 

We present a molecular theory of hydration that now makes possible a unification of these diverse views of the role of water in protein stabilization. The central element in our development is the potential distribution theorem. We discuss both its physical basis and statistical thermodynamic framework with applications to protein solution thermodynamics and protein folding in mind. To this end, we also derive an extension of the potential distribution theorem, the quasi-chemical theory, and propose its implementation to the hydration of folded and unfolded proteins. Our perspective and current optimism are justified by the understanding we have gained from successful applications of the potential distribution theorem to the hydration of simple solutes. A few examples are given to illustrate this point. 

For obvious reasons, we need to introduce surface contributions in the thermodynamic framework. Typically, in interface thermodynamics, the area in the system, e.g. the area of an air-water interface, is a state variable that can be adjusted by the observer while keeping the intensive variables (such as the temperature, pressure and chemical potentials) fixed. The unique feature in selfassembling systems is that the observer cannot adjust the area of a membrane in the same way, unless the membrane is put in a frame. Systems that have self-assembly characteristics are conveniently handled in a setting of thermodynamics of small systems, developed by Hill [12], and applied to surfactant self-assembly by Hall and Pethica [13]. In this approach, it is not necessary to make assumptions about the structure of the aggregates in order to define exactly the equilibrium conditions. However, for the present purpose, it is convenient to take the bilayer as an example. 

Needless to say, a simplified model leads to corresponding thermodynamic quantities, i.e. not all correlations are included. However, the thermodynamic framework itself is fully internally consistent. This is an important observation, because such a model can for this reason be of use to establish the thermodynamic feasibility of what-if questions. Full control over the absolute deviations from the true thermodynamic behaviour is unfortunately not possible. The approach ignores important (cooperative) fluctuations, and it is expected that especially near phase transitions the approach may give only qualitative results. In particular, comparison of SCF results with experiments or with simulation data can lead to insights into how rigorous the method is. 

The first rigorous method for weak electrolyte solutions was that of Edwards et al. ( 5). Because comparisons with the models of other workers will be made, the thermodynamic framework will be outlined and the assumptions that were made stated. For a single solute which dissociates only in the aqueous solution, the model is based on four principles ... 

Recently, there have been a number of significant developments in the modeling of electrolyte systems. Bromley (1), Meissner and Tester (2), Meissner and Kusik (2), Pitzer and co-workers (4, ,j5), and" Cruz and Renon (7j, presented models for calculating the mean ionic activity coefficients of many types of aqueous electrolytes. In addition, Edwards, et al. (8) proposed a thermodynamic framework to calculate equilibrium vapor-liquid compositions for aqueous solutions of one or more volatile weak electrolytes which involved activity coefficients of ionic species. Most recently, Beutier and Renon (9) and Edwards, et al.(10) used simplified forms of the Pitzer equation to represent ionic activity coefficients. 

Table 1. Thermodynamic Framework of Representation of Vapor-Liquid Equilibria of Weak Electrolytes...

In this paper I will discuss the tools we have developed and are continuing to develop and enhance. In particular, I will describe the basic thermodynamic framework which we have synthesized - nd how this framework is implemented into a computer software system that ... 

In developing the thermodynamic framework for ECES, we attempted to synthesize computer software that would correctly predict the vapor-liquid-solid equilibria over a wide range of conditions for multicomponent systems. To do this we needed a good basis which would make evident to the user the chemical and ionic equilibria present in aqueous systems. We chose as our cornerstone the law of mass action which simply stated says "The product of the activities of the reaction products, each raised to the power indicated by its numerical coefficient, divided by the product of the activities of the reactants, each raised to a corresponding power, is a constant at a given temperature. ... 

I have outlined the basic thermodynamic framework we adopted and in many cases where we have extended the work of several of the researchers present at this svmposium. The particular structure was chosen for several reasons including the following ... 

Which ever path is followed, the thermodynamic framework which I earlier described is used as the basis along with data obtained from a large data-base created by using the Data Regression Program Block or Data Estimation Block, to finally describe the system. 

To use the ECES Program Builder, the user creates computer input as shown in Figure 1. The input consists of three parts The names of the components in the inlet stream to the flash unit, the names of the species present in all the outlet streams - vapor, aqueous, or solids, and the chemical and/or ionic equilibria of interest. From this input the Program Builder Block using the imbedded thermodynamic framework described earlier writes the model description as given in Figure 2. 

These expressions were developed automatically from the user input and the embedded thermodynamic framework. In addition a separate routine was written to compute the equilibrium constants as a function of emperature and the B,C, and D coefficients for each salt pair. If in the associated ECES database, certain species were not already present, a message such as... 

These sort of problems make it difficult to obtain reliable high temperature data on the aqueous chemistry of transition metal ions. Unfortunately the necessary timescales for even the simpler experimental studies are frequently too long for a Ph.D. student to make reasonable progress in 3 years from scratch or for industrial researchers to make much reportable progress before the patience of those supporting the work is exhausted. Results can be reported far more rapidly from, for example, corrosion experiments and since corrosion theories are in general of so little predictive value, each relevant alloy/electrolyte combination needs its own study. In such circumstances it is hardly surprising that thermodynamic studies have been (with a few notable exceptions) relatively poorly supported, while corrosion data continue to be amassed without any reliable thermodynamic framework within which to understand them. 

Once this common thermodynamic framework is established for the solubility of U02 under nominally reducing conditions, we have to ascertain the most probable pathway for the oxidative alteration of U02 spent fuel in geological repository conditions. There is a large body of evidence on the processes involved in the oxidative alteration of natural uraninites and unirradiated U02. Long-term unsaturated tests performed by Wronckiewicz et al. (1992) on groundwater from Yucca Mountain (the so-called J-13 groundwater), indicated that the formation of schoepite, as described by process (20) and (21), occurs, but is a transient event and that the alteration proceeds towards the precipitation of... 

Application of the collected thermodynamic data to model the oxidative alteration pathway of U02 under repositoiy conditions by using the PHREEQC code (Parkhurst Appelo 1999) is given in Fig. 1 la and b. Once the thermodynamic framework is set for the geochemical evolution of the repositoiy system, we have to take into consideration that for many of the processes involved, there will be some kinetic constraints. This is illustrated by Table 2, where a comparison of the expected lifetime for some of the phases expected in the repositoiy system is made. 

In this study, a thermodynamic framework has been presented for the calculation of vapor-liquid equilibria for binary solvents containing nonvolatile salts. From an appropriate definition of a pseudobinary system, infinite dilution activity coefficients for the salt-containing system may be estimated from a knowledge of vapor pressure lowering, salt-free infinite dilution activity coefficients, and a single system-dependent constant. Parameters for the Wilson equation may be determined from the infinite dilution activity coefficients. 

In addition, several alternative formulations of thermodynamic geometry have been presented, starting from entropy-based (or other) fundamental equations (see Sections 5.4 and 5.5). From the equilibrium thermodynamics viewpoint, these alternative formulations are completely equivalent, and each could be considered a special case of the general transformations outlined in Section 11.4. Nevertheless, each alternative may suggest distinct statistical-mechanical origins, Riemannian paths, or other connotations that make it preferable for applications outside the equilibrium thermodynamics framework. 

In the last ten years, much work in this field has been performed in the groups of Radosz [78-86] and McHugh [87-92]. It is not the purpose of this chapter to review all the results obtained by the quoted researchers. It is important to note that all these results can be successfully placed in the thermodynamic framework described above. 

The constitutive equations use a thermodynamic framework, that in fact embodies not only purely mechanical aspects, but also transfers of masses between the phases and diffusion of matter through the extrafibrillar phase. Since focus is on the chemo-mechanical couplings, we use experimental data that display different salinities. The structure of the constitutive functions and the state variables on which they depend are briefly motivated. Calibration of material parameters is defined and simulations of confined compression tests and of tree swelling tests with a varying chemistry are described and compared with available data in [3], The evolution of internal entities entering the model, e.g. the masses and molar fractions of water and ions, during some of these tests is also documented to highlight the main microstructural features of the model. 

Abstract In a thermodynamic framework which exploits the entropy inequality to obtain constitutive equations, it is common practice to assume charge neutrality and enforce this restriction using Lagrange multipliers. In this paper we show that the Lagrange multiplier used to enforce charge neutrality does not correspond to any known physical parameter, raising the question of whether charge neutrality can really be enforced. 

Fang, X., Pan, T., and Sosnick, T. R. (1999). A thermodynamic framework and coopera-tivity in the tertiary folding of a Mg2+-dependent ribozyme. Biochemistry 38,... 

In Chapters 6, 7, and 8, the thermodynamic framework is successively apphed to phase transformations of single-component systems, chemical reactions, and ideal solutions. Included are discussions of the thermodynamics of open systems, the phase rule, and colligative properties. Chapter 9 gives the framework for discussing nonideal multicomponent systems and describes a... 

Equations (l)-(4) provide the basic statistical thermodynamic framework necessary to deal with the protein folding problem. Several years ago, Freire and Biltonen (1978a) showed that scanning calorimetry data could be used to evaluate the protein folding/unfolding partition function experimentally by a double integration procedure ... 

Allen [240], Costanzo et al. [241], and Krajcinovic [242-244] that thermodynamics of appropriate internal state variables for damage in composites under mechanical loads were addressed. Similar to the metal plasticity theoreticians, these researchers employed Coleman and Gurtin s [11] thermodynamic framework to determine the kinetic equations. However, unlike the metal internal state variable community that could quantify the evolution equations for dislocations and damage, these polymer-based composites theoreticians did not propose evolutionary rate equations, but just damage state equations. 

The key result in terms of a sensor is that specific interactions, as sought for biochemical sensors (6), may be sufficiently strong that a coordination-type model applies. Note that this does not contradict the activity arguments of the previous section, but is a special case within the general thermodynamic framework. Under these special circumstances, the polymer will be "saturated with the target species, and film composition will not depend on solution concentration, except at a very low level. 

---