Thermodynamic modeling methods

At this point in the inside-out method, the revised column profiles of temperature and phase compositions are used in the outer loop with the complex SRK thermodynamic models to compute updates of the approximate K and H constants. Then only one inner-loop iteration is required to obtain satisfactory convergence of the energy equations. The K and H constants are again updated in the outer loop. After one inner-loop iteration, the approximate K and H constants are found to be sufficiently close to the SRK values that overall convergence is achieved. Thus, a total of only 3 outer-loop iterations and 4 inner-loop iterations are required. 

The implicit LS, ML and Constrained LS (CLS) estimation methods are now used to synthesize a systematic approach for the parameter estimation problem when no prior knowledge regarding the adequacy of the thermodynamic model is available. Given the availability of methods to estimate the interaction parameters in equations of state there is a need to follow a systematic and computationally efficient approach to deal with all possible cases that could be encountered during the regression of binary VLE data. The following step by step systematic approach is proposed (Englezos et al. 1993)... 

Higuchi and Lachman [122] pioneered the approach of improving drug stability by complexation. They showed that aromatic esters can be stabilized in aqueous solutions in the presence of xanthines such as caffeine. Thus, the half-lives of benzocaine, procaine hydrochloride, and tetracaine are increased by approximately two- to fivefold in the presence of 2.5% caffeine. This increase in stability is attributed to the formation of a less reactive complex between caffeine and the aromatic ester. Professor K. A. Connors has written a comprehensive textbook that describes methods for the measurement of binding constants for complex formation in solution—along with discussions of pertinent thermodynamics, modeling statistics,... 

Considerable use has been made of the thermodynamic perturbation and thermodynamic integration methods in biochemical modelling, calculating the relative Gibbs energies of binding of inhibitors of biological macromolecules (e.g. proteins) with the aid of suitable thermodynamic cycles. Some applications to materials are described by Alfe et al. [11]. 

The experimental basis of sorption studies includes structural data (SANS, SAXS, USAXS), isopiestic vapor sorption isotherms,i and capillary isotherms, measured by the method of standard porosimetry. i 2-i44 Thermodynamic models for water uptake by vapor-equilibrated PEMs have been suggested by various groupThe models account for interfacial energies, elastic energies, and entropic contributions. They usually treat rate constants of interfacial water exchange and of bulk transport of water by diffusion and hydraulic permeation as empirical functions of temperature. 

Thermodynamic modelling of solution phases lies at the core of the CALPHAD method. Only rarely do calculations involve purely stoichiometric compounds. The calculation of a complex system which may have literally 100 different stoichiometric substances usually has a phase such as the gas which is a mixture of many components, and in a complex metallic system with 10 or 11 alloying elements it is not unusual for all of the phases to involve solubility of the various elements. Solution phases will be defined here as any phase in which there is solubility of more than one component and within this chapter are broken down to four types (1) random substitutional, (2) sublattice, (3) ionic and (4) aqueous. Others types of solution phase, such as exist in polymers or complex organic systems, can also be modelled, but these four represent the major types which are currently available in CALPHAD software programmes. 

Most of this book concerns the development and application of theoretical thermodynamic models, as these are the basis of the CALPHAD method. However, none of this would be possible without the existence of the computational methods and software which allow these models to be applied in practice. In essence, the issues involved in computational methods are less diverse and mainly revolve around Gibbs energy minimisation. In addition, there are optimiser codes which are used for the thermodynamic assessment of phase equilibria. The essential aim of these codes is to reduce the statistical error between calculated phase equilibria, thermodynamic properties and the equivalent experimentally measured quantities. 

Surfactant Activity in Micellar Systems. The activities or concentrations of individual surfactant monomers in equilibrium with mixed micelles are the most important quantities predicted by micellar thermodynamic models. These variables often dictate practical performance of surfactant solutions. The monomer concentrations in mixed micellar systems have been measured by ultraf i Itration (I.), dialysis (2), a combination of conductivity and specific ion electrode measurements (3), a method using surface tension of mixtures at and above the CMC <4), gel filtration (5), conductivity (6), specific ion electrode measurements (7), NMR <8), chromatograph c separation of surfactants with a hydrophilic substrate (9> and by application of the Bibbs-Duhem equation to CMC data (iO). Surfactant specific electrodes have been used to measure anionic surfactant activities in single surfactant systems (11.12) and might be useful in mixed systems. ... 

Special classes of apparatus are used for the determination of particular thermodynamic properties, such as activity coefficients at infinite dilution, Henry s constants, or partial molar volumes at infinite dilution [105,106]. These data, together with a thermodynamic model, can be used for the calculation of the compositions of the coexisting phases at equilibrium, and for that reason - in this context - these methods are considered as indirect methods of measurement. 

Method of Cao, Knudsen, Fredenslund, and Rasmussen The method of Cao et al. [23] is based on a statistical thermodynamic model for pure liquids and liquid mixtures. It requires the input of the compound properties VM and AHv and two... 

The group contribution method UNIFAC [18] is based on the UNIQUAC thermodynamic model as well. It thus suffers from the same thermodynamic approximations as UNIQUAC, especially for strong interactions in the infinite dilution limit. 

Roberts refined the fragment constant method for estimating log Kow for LAS (Roberts, 1989), as well as for CnEOms (Roberts, 1991) by introducing factors which correct for the effective size of the hydrophobic moiety of the surfactant. Kiewiet (1996) employed a thermodynamic model to estimate log Kow of CnEOms. However, the quality of the... 

Many polymer blends or block polymer melts separate microscopically into complex meso-scale structures. It is a challenge to predict the multiscale structure of polymer systems including phase diagram, morphology evolution of micro-phase separation, density and composition profiles, and molecular conformations in the interfacial region between different phases. The formation mechanism of micro-phase structures for polymer blends or block copolymers essentially roots in a delicate balance between entropic and enthalpic contributions to the Helmholtz energy. Therefore, it is the key to establish a molecular thermodynamic model of the Helmholtz energy considered for those complex meso-scale structures. In this paper, we introduced a theoretical method based on a lattice model developed in this laboratory to study the multi-scale structure of polymer systems. First, a molecular thermodynamic model for uniform polymer system is presented. This model can... 

To establish the molecular thermodynamic model for uniform systems based on concepts from statistical mechanics, an effective method by combining statistical mechanics and molecular simulation has been recommended (Hu and Liu, 2006). Here, the role of molecular simulation is not limited to be a standard to test the reliability of models. More directly, a few simulation results are used to determine the analytical form and the corresponding coefficients of the models. It retains the rigor of statistical mechanics, while mathematical difficulties are avoided by using simulation results. The method is characterized by two steps (1) based on a statistical-mechanical derivation, an analytical expression is obtained first. The expression may contain unknown functions or coefficients because of mathematical difficulty or sometimes because of the introduced simplifications. (2) The form of the unknown functions or unknown coefficients is then determined by simulation results. For the adsorption of polymers at interfaces, simulation was used to test the validity of the weighting function of the WDA in DFT. For the meso-structure of a diblock copolymer melt confined in curved surfaces, we found from MC simulation that some more complex structures exist. From the information provided by simulation, these complex structures were approximated as a combination of simple structures. Then, the Helmholtz energy of these complex structures can be calculated by summing those of the different simple structures. 

---