Thermodynamic equilibrium, computational

The predictions of CO are based on a thermodynamic equilibrium computer program in which the chemical equilibrium compositions for assigned thermodynamic states are calculated. Most often temperature and pressure, calculated from the main CFD code, are used to specify the thermodynamic state. The method for evaluating the equilibrium compositions is based on the minimization of Gibbs energy as mentioned earlier (Wang et al., 2006). 

As can be seen, conditions in a flowing reactor, even the simplest such as a tube, may be far from the thermodynamic equilibrium conditions predicted by the equilibrium computer programs. However, in the diffusion controlled range, it is possible to use as the driving force for diffusion, the difference between an assumed equi-... 

In this work, we developed the safeguard active-set method by modifying the active-set method for thermodynamic equilibrium in order to include the classical nucleation theory. At tn, assume that the partition ( (r ), M(t ), N(t ) and the crystallization time tciyst(t ) forM(t ) are known. For a new feed vector and RH at Vu compute W(tn+i), M(t i), N(t + )) and tciyst(t +i) as follows ... 

When the simulation of deep-well temperatures, pressures, and salinities is imposed as a condition, the number of codes that may be of value is reduced to a much smaller number. Nordstrom and Ball121 recommend six references as covering virtually all the mathematical, thermodynamic, and computational aspects of chemical-equilibrium formulations (see references 123-128). Recent references on modeling include references 45, 63, 70, 129, and 130. 

McFarland et al. recently [1] published the results of studies carried out on 22 crystalline compounds. Their water solubilities were determined using pSOL [21], an automated instrument employing the pH-metric method described by Avdeef and coworkers [22]. This technique assures that it is the thermodynamic equilibrium solubility that is measured. While only ionizable compounds can be determined by this method, their solubilities are expressed as the molarity of the unionized molecular species, the intrinsic solubility, SQ. This avoids confusion about a compound s overall solubility dependence on pH. Thus, S0, is analogous to P, the octanol/water partition coefficient in both situations, the ionized species are implicitly factored out. In order to use pSOL, one must have knowledge of the various pKas involved therefore, in principle, one can compute the total solubility of a compound over an entire pH range. However, the intrinsic solubility will be our focus here. There was one zwitterionic compound in this dataset. To obtain best results, this compound was formulated as the zwitterion rather than the neutral form in the HYBOT [23] calculations. 

We then use a Feautrier scheme [4] to perform spectral line formation calculations in local thermodynamic equilibrium approximation (LTE) for the species indicated in table 1. At this stage we consider only rays in the vertical direction and a single snapshot per 3D simulation. Abundance corrections are computed differentially by comparing the predictions from 3D models with the ones from ID MARCS model stellar atmospheres ([2]) generated for the same stellar parameters (a microturbulence = 2.0 km s-1 is applied to calculations with ID models). 

Such information has been stored within a database in the ECES system. Then, using user input in the form of equation (1), ECES writes the expression for computing the thermodynamic equilibrium constant as a function of temperature to a file where it will eventually become part of a program to solve the many equilibria that might describe a complex system. 

Nitric oxide has a very low ionization potential and could ionize at flame temperatures. For a normal composite solid propellant containing C—H—O—N—Cl—Al, many more products would have to be considered. In fact if one lists all the possible number of products for this system, the solution to the problem becomes more difficult, requiring the use of advanced computers and codes for exact results. However, knowledge of thermodynamic equilibrium constants and kinetics allows one to eliminate many possible product species. Although the computer codes listed in Appendix I essentially make it unnecessary to eliminate any product species, the following discussion gives one the opportunity to estimate which products can be important without running any computer code. 

The last few years have seen a minor revolution in determining solar and stellar abundances (Asplund, 2005). Much of the previous work assumed that the spectral lines originate in local thermodynamic equilibrium (LTE), and the stellar atmosphere has been modeled in a single dimension. Since 2000, improved computing power has permitted three-dimensional modeling of the Sun s atmosphere and non-LTE treatment of line formation. The result has been significant shifts in inferred solar abundances. 

Gas detonation at reduced initial pressures were studied by Vasil ev et al (Ref 8). They point out the errors in glibly comparing ideal lossless onedimensional computations with measurements made in 3-dimensicnal systems. We quote In an ideal lossless detonation wave, the Chapman-Jouguet plane is identified with the plane of complete chemical and thermodynamic equilibrium. As a rule, in a real detonation wave the Chapman-Jouguet state is assumed to be the gas state behind the front, where the measurable parameters are constant, within the experimental errors. It is assumed that, in the one-dimensional model of the detonation wave in the absence of loss, the conditions in the transient rarefaction wave accompanying the Chapman-Jouguet plane vary very slowly if the... 

A puzzling problem was posed by Levinthal many years ago.329 We usually assume that the peptide chain folds into one of the most stable conformations possible. However, proteins fold very rapidly. Even today, no computer would be able, in our lifetime, to find by systematic examination the thermodynamically most stable conformation.328 It would likewise be impossible for a folding protein to "try out" more than a tiny fraction of all possible conformations. Yet folded and unfolded proteins often appear to be in a thermodynamic equilibrium Experimental results indicate that denatured proteins are frequently in equilibrium with a compact denatured state or "molten globule" in which hydrophobic groups have become clustered and some secondary structures exists.330-336 From this state the polypeptide may rearrange more slowly through other folding intermediates to the final "native conformation."3363 3361 ... 

ORM assumes that the atmosphere is in local thermodynamic equilibrium this means that the temperature of the Boltzmann distribution is equal to the kinetic temperature and that the source function in Eq. (4) is equal to the Planck function at the local kinetic temperature. This LTE model is expected to be valid at the lower altitudes where kinetic collisions are frequent. In the stratosphere and mesosphere excitation mechanisms such as photochemical processes and solar pumping, combined with the lower collision relaxation rates make possible that many of the vibrational levels of atmospheric constituents responsible for infrared emissions have excitation temperatures which differ from the local kinetic temperature. It has been found [18] that many C02 bands are strongly affected by non-LTE. However, since the handling of Non-LTE would severely increase the retrieval computing time, it was decided to select only microwindows that are in thermodynamic equilibrium to avoid Non-LTE calculations in the forward model. 

Figure 1. Thermodynamic equilibrium in atmospheres of varying elemental proportions. The ternary diagram shows all compositions of systems containing carbon, hydrogen, and oxygen (each point represents 100% of the three components). Lower curves indicate the potential formation of solid carbon if equilibrium could be attained. Dashed curve holds at 500°K., the continuous one at 700°K. The upper lines indicate the asphalt threshold, the dashed one at 500° K., and the continuous one at 700° K. Above this threshold, thermodynamic equilibrium favors the formation of large proportions of polycyclic aromatic compounds ( asphalt ) ana a lesser increase of most of the other families of compounds. The dots through points A to C indicate the points used in the computations for Figure 2 (6).

In this section we collect some computed results when simulating industrial steam reformers. We compare the actual plant outputs with those obtained by simulation using the three models that we have developed earlier. We investigate both the close to thermodynamic equilibrium case and the far from thermodynamic equilibrium case. 

The input of the problem requires total analytically measured concentrations of the selected components. Total concentrations of elements (components) from chemical analysis such as ICP and atomic absorption are preferable to methods that only measure some fraction of the total such as selective colorimetric or electrochemical methods. The user defines how the activity coefficients are to be computed (Davis equation or the extended Debye-Huckel), the temperature of the system and whether pH, Eh and ionic strength are to be imposed or calculated. Once the total concentrations of the selected components are defined, all possible soluble complexes are automatically selected from the database. At this stage the thermodynamic equilibrium constants supplied with the model may be edited or certain species excluded from the calculation (e.g. species that have slow reaction kinetics). In addition, it is possible for the user to supply constants for specific reactions not included in the database, but care must be taken to make sure the formation equation for the newly defined species is written in such a way as to be compatible with the chemical components used by the rest of the program, e.g. if the species A1H2PC>4+ were to be added using the following reaction ... 

Since equilibrium constants are defined by ion activities, which are defined by their concentrations and coefficients (see equation 4.15), they do not include ion pairing or complexation effects. In a multi-ion and multiligand solution, where ion pairing is common, it is necessary to use thermodynamic equilibrium constants to convert the ion-pair concentrations to concentrations of free ions. This equilibrium constant (Kc) is defined by concentrations, making it useful to compute ion speciation. The thermodynamic equilibrium constant (K q) used in calculating Kc is based on the following conditions I = 0 m, 25°C and 1 atm. Thus, Kc is defined by the following equation ... 

The susceptibility and polarizability differ in that the latter deals with single >i taolecules, and the former with an entire medium. Further, die polarizability can jjjbe defined for the system in a particular quantum state, whereas the susceptibility j ienerally refers to the bulk system, often in thermodynamic equilibrium. Typically, ltsn [212] the susceptibility includes the number of particles per unit volume p and pjl. average over populated system energy levels. Below we compute the polariz-oility and susceptibility for neither of these cases. Rather, we extend the standard rition of the polarizability to include a molecule initially in a superposition state. 

The strategy of design, illustrated in Figure 8.1, consists of an evolutionary search of the feasible design space by means of a systematic combination of thermodynamic analysis, computer simulation and only limited experiments. The approach is generic for developing a RD process, at least for similar systems. The first element of similarity is the existence of an equilibrium reaction with water as product This raises the problem of possible aqueous-phase segregation. The second element is the similarity of thermodynamics properties over a class of substrates. However, while the fatty acids and fatty esters manifest a certain... 

---