Computational procedures condensation

Because of its diatomic nature and permanent quadrupole moment, the physisorp-tion of nitrogen at 77 K presents special problems. The application of DFT is facilitated if the molecules are assumed to be spherical, which was the approach originally adopted by Seaton et al. (1989) and also by Lastoskie et al. (1993). The analytical procedures already outlined in Chapter 7 (Section 7.6) do not depend on the meniscus curvature and are in principle applicable to both capillary condensation and micropore filling. The non-local version of the mean field theory (NLDFT), which was used by Lastoskie, gave excellent agreement with computer simulation when applied to the carbon slit pore model. However, as pointed out earlier, these computational procedures are not entirely independent since they involve the same model parameters. 

For multicomponent mixtures, graphical representations of properties, as presented in Chapter 3, cannot be used to determine equilibrium-stage requirements. Analytical computational procedures must be applied with thermodynamic properties represented preferably by algebraic equations. Because mixture properties depend on temperature, pressure, and phase composition(s), these equations tend to be complex. Nevertheless the equations presented in this chapter are widely used for computing phase equilibrium ratios (K-values and distribution coefficients), enthalpies, and densities of mixtures over wide ranges of conditions. These equations require various pure species constants. These are tabulated for 176 compounds in Appendix I. By necessity, the thermodynamic treatment presented here is condensed. The reader can refer to Perry and Chilton as well as to other indicated sources for fundamental classical thermodynamic background not included here. 

A computer solution was obtained as follows. The only initial assumptions are a condenser outlet temperature of 65 F and a bottoms-prodiict temperature of 165 F, The bubble-point temperature of the feed is computed as 123,5 F, In the initiahzation procedure, the constants A and B in (13-106) for inner-loop calcu-... 

Samples are injected into the vaporizer by a metering pump or manually with septum injection the manual injection procedure is intended for method development. The sample gas mixture then passes through the chromatographic column where the sample compounds separate. Fractions pass through the thermal conductivity detector and then to a condenser collection manifold where up to five fractions can be collected. Complete control of the system is achieved via a mini-computer. 

The need to reliably describe liquid systems for practical purposes as condensed matter with high mobility at a given finite temperature initiated attempts, therefore, to make use of statistical mechanical procedures in combination with molecular models taking into account structure and reactivity of all species present in a liquid and a solution, respectively. The two approaches to such a description, namely Monte Carlo (MC) simulations and molecular dynamics (MD), are still the basis for all common theoretical methods to deal with liquid systems. While MC simulations can provide mainly structural and thermodynamical data, MD simulations give also access to time-dependent processes, such as reaction dynamics and vibrational spectra, thus supplying — connected with a higher computational effort — much more insight into the properties of liquids and solutions. 

In order to estimate the kinetic parameters for the addition and condensation reactions, the procedure proposed in [11, 14] has been used, where the rate constant kc of each reaction at a fixed temperature of 80°C is computed by referring it to the rate constant k° at 80°C of a reference reaction, experimentally obtained. The ratio kc/k°, assumed to be temperature independent, can be computed by applying suitable correction coefficients, which take into account the different reactivity of the -ortho and -para positions of the phenol ring, the different reactivity due to the presence or absence of methylol groups and a frequency factor. In detail, the values in [11] for the resin RT84, obtained in the presence of an alkaline catalyst and with an initial molar ratio phenol/formaldehyde of 1 1.8, have been adopted. Once the rate constants at 80°C and the activation energies are known, it is possible to compute the preexponential factors ko of each reaction using the Arrhenius law (2.2). 

---