Extended Corresponding States Methods

Huber and Ely expanded the extended corresponding states models reported by Leach and Ely to predict the thermodynamic properties of mixtures assuming the mixture behaves as a hypothetical equivalent pure substance. To determine mixture properties, the states or properties of the mixture, identified by the subscript x, and those of a reference fluid, designated by subscript ref, must be in correspondence  

Corresponding states are found with shape factors, and 6, that relate the reduced properties of the mixture to those of the reference fluid. The thermodynamic properties for the mixture are reduced by the pseudo-critical parameters pc,x and defined by 

---

For higher gas densities, the Lee-Kesler method described below provides excellent predictions for nonpolar and slightly polar fluids. Extended four-parameter corresponding-states methods are available for polar and slightly associating compounds. 

Pai-Panandiker, R. S., Nieto de Castro, C. A., Marrucho, I. M., and Ely, J. E, Development of an extended corresponding states principle method for volumetric property predictions based on Lee-Kesler reference fluid, Int. J. Themwphys., 23, 771-785 (2002). 

The multi-fluid approach can always be used with corresponding states methods for well-defined mixtures. In the one-fluid approach, however, a mixing rule must be proposed for each of the input parameters. For the Petersen et al. [22] corresponding states model discussed earlier, the following relations [31] are used to extend the model to mixtures ... 

The quantities 9 and cp which appear in these equations are referred to as shape factors and we refer to this method as the extended corresponding-states theory (ECST). Several review papers have been published that focus on this approach, the most extensive of which are those of Leland and Chappelear, Rowlinson and Watson and Mentzer et al ... 

The extended corresponding-states theory incorporating shape factors can predict thermodynamic properties of mixtures precisely, especially if the exact or saturation- boundary methods are used. Due to the higher level of complexity of this method, recent applications have been limited to studies where very low uncertainty is of importance. Advances in the development of precise equations of state for a wide variety of substances (for example, heavier hydrocarbons, refrigprants and polar compounds such as alcohols, water and ammonia) have enabled researchers to extend the extended corresponding states methodology to a wide variety of systems and this trend will continue in the future. 

The extended principle of corresponding state methods are being further developed. Wider experience with this approach might produce earlier adoption of new methods of prediction because most practicing engineers have a better intuitive understanding of these methods. 

Lee and Kesler (1975) extended the Benidict- Webb-Rubin equation to a wider variety of substances, using the principle of corresponding states. The method was modified further by Plocker et al. (1978). 

For transition elements like Pt and Au the linear orbital extension to the LAPW method [43] has been used. We have employed the procedure proposed in [20], in which the 5p-states for Au and Pt are included in the core for total energy calculations, but corresponding local orbitals are also included in the basis for the valence states in order to allow the basis functions for the actual valence electrons to orthogonalize to the extended core states. 

At low and moderate pressures, the viscosity of a gas is nearly independent of pressure and can be correlated for engineering purposes as a function of temperatnre only. Eqnations have been proposed based on kinetic theory and on corresponding-states principles these are reviewed in The Properties of Gases and Liquids [15], which also inclndes methods for extending the calculations to higher pressures. Most methods contain molecular parameters that may be fitted to data where available. If data are not available, the parameters can be estimated from better-known quantities such as the critical parameters, acentric factor, and dipole moment. The predictive accuracy for gas viscosities is typically within 5%, at least for the sorts of small- and medinm-sized, mostly organic, molecules used to develop the correlations. 

In 1973 Andersen and Woolley [1.25] extended the LCMTO method to molecular calculations. At the end of their paper they introduced that choice of MTO tail, i.e. proportional to J = 9i/j/9E, which in a natural fashion ensured orthogonality to the core states and at the same time led to an accurate and elegant formulation of linear methods. The resulting, technique was immediately developed in a paper by Andersen [1.26] which, in a condensed form, contains most of what one need know about the simple concepts of linear band theory. Thus, we find here the KKR equation within the atomic-sphere approximation at this stage is called ASM the LCMTO secular matrix, latter called the LMTO matrix the energy-independent structure constants and the canonical bands and the Laurent expansion of the logarithmic-derivative function and the corresponding potential parameters. 

For nonassociating molecules, the methods of Tsonopoulos [14] and Hayden and O Connell [15] yield pretty good results for the estimation of the second virial coefficient. The first method works with the three-parameter corresponding-states principle (Section 2.4.4) the second method takes into account several molecular effects. In case of vapor phase association (Section 13.2), the model of Hayden-O Connell is extended by a chemical equilibrium term, whose parameters have to be fitted to experimental data. 

A major difference between ODE s and SDE s is the inclusion of a stochastic term in the latter, which allows uncertainty to be accomodated, and which, if the constitutive equations of interest are reformulated as additional state equations, allows estimates of the corresponding state variables to be computed from the experimental data. The specific approach used for this purpose involves parameter estimation and subsequent state estimation by means of methods based on the extended Kalman filter (EKF). 

The initial conditions for the quantum dynamics corresponding to those in the preceding section are as follows. Suppose that a total wavefunction totai(r, i , t) is expanded as in Eq. (6.104) and we are interested only in the nuclear wavepackets xi R,t) and x2 R,t). We propagate them with the extended split operator method. [77] The initial nuclear wavepackets xi R,t) are chosen to be a coherent-type Gaussian function only on the adiabatic ground state as... 

The Lee-Kesler method for correlating gas phase properties (reference equation and corresponding states approach) has been extended to vapor-liquid equilibria calculations. The calculation of accurate phase equilibria results is more difficult than for the properties of single phase fluids. 

The thermal conductivity of a multicomponent mixture of monatomic species therefore requires a knowledge of the diermal conductivity of the pure components and of three quantities characteristic of the unlike interaction. The final three quantities may be obtained by direct calculation from intermolecular potentials, whereas the interaction thermal conductivity, Xgg, can also be obtained by means of an analysis of viscosity and/or diffusion measurements through equations (4.112) and (4.125) or by the application of equation (4.122) to an analysis of the thermal conductivity data for all possible binary mixtures, or by a combination of both. If experimental data are used in the prediction it may be necessary to estimate both and This is readily done using a realistic model potential or the correlations of the extended law of corresponding states (Maitland et al. 1987). Generally, either of these procedures can be expected to yield thermal conductivity predictions with an accuracy of a few percent for monatomic systems. Naturally, all of the methods of evaluating the properties of the pure components and the quantities characteristic of binary interactions that were discussed in the case of viscosity are available for use here too. 

Several empirical methods to extend the critical enhancement term to mixtures have been explored (see also Chapters 6 and 15). The parameters in the mode-coupling approach mentioned above have been made composition-dependent and the results are reasonable. Luettmer-Strathmann (1994) has recently reported a new mode-coupling solution which describes the critical enhancements to the transport properties of fluid mixtures. Corresponding states algorithms, based on the mode-coupling solutions, have also been used to describe the thermal conductivity of mixtures (Huber et al 1992). 

The same methods have been extended to the dense fluid region by corresponding-states approaches. Errors of 5 to 10% are typical for estimations in the dense fluid region. 

We will state in this chapter the mathematical task of parameter identification and discuss the corresponding numerical methods. Techniques from various branches of numerical mathematics are required, e.g. numerical solution of differential equations, numerically solving nonlinear problems especially large-scale constrained nonlinear least squares problem. Thus, some of the methods discussed in the previous chapters will reappear here. We will see how parameter identification problems can be treated efficiently by boundary value problem (BVP) methods and extend the discussion of solution techniques for initial value problems (IVPs) to those for BVPs. 

As the load varies, it will be assumed that the contact patch will pass through a one-parameter family of states, as shown schematically in Fig. 3.2. This assumption will be justified later on the basis that it enables the problem to be solved. Furthermore, it will be shown that the one-parameter family of states is in fact the family of possible elastic states. The fact that C t) is a one-parameter family means that the explicit formalism developed for repetitive expansion and contraction in Sects. 2.6 and 3.10 may be used, as opposed to the more general method summarized in Sect. 2.6 in the context of the Extended Correspondence Principle, which is applicable to any situation where the boundary regions are expanding and contracting in time. 

---