The Equilibrium Flux Method for the calculation of flows with non-equilbrium chemical reactions

The Equilibrium Flux Method [1] is a kinetic theory based finite volume method for calculating the flow of a compressible ideal gas. It is shown here that, in effect, the method solves the Euler equations with added pseudo-dissipative terms and that it is a natural upwinding scheme. The method can be easily modified so that the flow of a chemically reacting gas mixture can be calculated. Results from the method for a one-dimensional non-equilibrium reacting flow are shown to agree well with a conventional continuum solution. Results are also presented for the calculation of a plane two-dimensional flow, at hypersonic speed, of a dissociating gas around a blunt-nosed body.

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

Analysis and application of an equilibrium model for in vitro bioassay systems with three components: Receptor, hormone and hormone-binding-protein

A simple theoretical framework is presented for bioassay studies using three component in vitro systems. An equilibrium model is used to derive equations useful for predicting changes in biological response after addition of hormone-binding-protein or as a consequence of increased hormone affinity. Sets of possible solutions for receptor occupancy and binding protein occupancy are found for typical values of receptor and binding protein affinity constants. Unique equilibrium solutions are dictated by the initial condition of total hormone concentration. According to the occupancy theory of drug action, increasing the affinity of a hormone for its receptor will result in a proportional increase in biological potency. However, the three component model predicts that the magnitude of increase in biological potency will be a small fraction of the proportional increase in affinity. With typical initial conditions a two-fold increase in hormone affinity for its receptor is predicted to result in only a 33% increase in biological response. Under the same conditions an Ii-fold increase in hormone affinity for receptor would be needed to produce a two-fold increase in biological potency. Some currently used bioassay systems may be unrecognized three component systems and gross errors in biopotency estimates will result if the effect of binding protein is not calculated. An algorithm derived from the three component model is used to predict changes in biological response after addition of binding protein to in vitro systems. The algorithm is tested by application to a published data set from an experimental study in an in vitro system (Lim et al., 1990, Endocrinology 127, 1287-1291). Predicted changes show good agreement (within 8%) with experimental observations. (C) 1998 Academic Press Limited.

Experimental study of phase equilibria in the "FeO"-ZnO-(CaO+SiO2) system with CaO/SiO2 weight ratios of 0.33, 0.93, and 1.2 in equilibrium with metallic iron

The pseudoternary sections FeO-ZnO-(CaO + SiO2) with CaO/SiO2 weight ratios of 0.33, 0.93, and 1.2 in equilibrium with metallic iron have been experimentally investigated in the temperature range from 1000 degreesC to 1300 degreesC (1273 to 1573 K). The liquidus surfaces in these pseudoternary sections have been experimentally determined in the composition range from 0 to 33 wt pct ZnO and 30 to 70 wt pct (CaO + SiO2). The sections contain primary-phase fields of wustite (FexZn1-xO1+y), zincite (ZnzFe1-zO), fayalite (Fu(w)Zn(2-w)SiO(4)), melilite (Ca2ZnuFe1-uSi2O7), willemite (ZnvFe2-vSiO4), dicalcium silicate (Ca2SiO4), pseudowollastonite and wollastonite (CaSiO3), and tridymite (SiO2). The phase equilibria involving the liquid phase and the solid solutions-have also been measured.

Non-equilibrium energy and momentum accommodation coefficients of Ar atoms scattered from Ni(001) in the thermal regime: A molecular dynamics study

Molecular dynamics simulations are used to study energy and momentum transfer of low-energy Ar atoms scattered from the Ni(001) surface. The investigation concentrates on the dependence of these processes on incident energy, angles of incidence and surface temperature. Energy transfer exhibits a strong dependence on the surface temperature, at incident energies below 500 meV, and incident angles close to specular incidence. Above 500 meV, the surface temperature dependence vanishes, and a limiting value in the amount of energy transferred to the surface is attained. Momentum exchange is investigated in terms of tangential and normal components. Both components exhibit a weak surface temperature dependence, but they have opposite behaviours at all incidence angles. In each component, momentum can be lost or gained following the interaction with the surface. (C) 1997 Elsevier Science B.V.

Studies of solute self-association by sedimentation equilibrium: allowance for effects of thermodynamic non-ideality beyond the consequences of nearest-neighbor interactions

A sedimentation equilibrium study of a-chymotrypsin self-association in acetate-chloride buffer, pH 4.1 I 0.05, has been used to illustrate determination of a dimerization constant under conditions where thermodynamic non-ideality is manifested beyond the consequences of nearest-neighbor interactions. Because the expressions for the experimentally determinable interaction parameters comprise a mixture of equilibrium constant and excluded volume terms, the assignment of reasonable magnitudes to the relevant virial coefficients describing non-associative cluster formation is essential for the evaluation of a reliable estimate of the dimerization constant. Determination of these excluded volume parameters by numerical integration over the potential-of-mean-force is shown to be preferable to their calculation by approximate analytical solutions of the integral for this relatively small enzyme monomer with high net charge (+ 10) under conditions of low ionic strength (0.05 M). (C) 2001 Elsevier Science B.V. All rights reserved.

The equilibrium penetration of monolayers. Is equilibrium really established?

In studying the penetration of water-soluble surfactants into water-insoluble monolayers the main theoretical problem is to find a relationship that would enable the amount of surfactant that has entered the monolayer to be calculated from a set of equilibrium surface pressure-area isotherms. Despite many attempts, no current theory gives satisfactory results when applied to experimental data (Langmuir 14 (1998) 2148). One possible reason is that equilibrium had not been established when the surface pressure-area curves were measured. The three experiments reported here suggest that equilibrium is extremely difficult to establish in such systems when the area is low or the surface pressure is high. The essence of these experiments is to try to reach the same final condition by two different routes. In the first route, the one nearly always used in equilibrium penetration measurements, the surfactant is injected under the expanded monolayer, which is then slowly compressed in steps, with time allowed at each step for a steady surface pressure to be attained. In the second procedure, the monolayer is first compressed to a high surface pressure and the surfactant then injected. A stepped expansion isotherm may then be observed. Surface pressure-area per monolayer molecule isotherms, reflection spectra, and slow neutron reflectivity data all show the same pattern: if the surfactant was allowed to penetrate while the monolayer was in an expanded state, it was not completely removed when the monolayer was compressed; but if the monolayer was in a highly compressed state when exposed to the surfactant little penetration took place until the film was expanded. There thus appear to be very large energy barriers to the ejection of surfactant from a compressed monolayer and to the penetration of surfactant into a compressed monolayer. Although these experiments have some limitations, it now seems likely that at least some of the penetration data used in evaluating the various thermodynamic treatments of equilibrium penetration were not equilibrium data. (C) 2001 Elsevier Science B.V. All rights reserved.

H-1 NMR study of the states of water in equilibrium Poly(HEMA-co-THFMA) hydrogels

The spin-spin relaxation times, T-2, of hydrated samples of poly(hydroxymethyl methacrylate), PHEMA, poly(tetrahydrofurfuryl methacrylate),PTHFMA, and the,corresponding HEMA-THFMA copolymers have been examined to probe the states of,the imbibed water in these polymers. The decay in the transverse magnetization of water. in fully hydrated samples of PHEMA, PTHFMA, and copolymers of HEMA and THFMA was described by a multiexponential function. The short component of T-2 was interpreted as water molecules that were strongly interacting with the polymer chains. The intermediate component of T-2 was assigned to water residing in the porous structure of the samples. The long component of T-2 was believed to arise from water residing in the remnants of cracks formed in the polymer network during water sorption.

PFG-NMR measurements of the self-diffusion coefficients of water in equilibrium poly(HEMA-co-THFMA) hydrogels

The self-diffusion coefficients for water in a series of copolymers of 2-hydroxyethyl methacrylate, HEMA, and tetrahydrofurfuryl methacrylate, THFMA, swollen with water to their equilibrium states have been studied at 310 K using PFG-NMR. The self-diffusion coefficients calculated from the Stejskal-Tanner equation, D-obs, for all of the hydrated polymers were found to be dependent on the NMR storage time, as a result of spin exchange between the proton reservoirs of the water and the polymers, reaching an equilibrium plateau value at long storage times. The true values of the diffusion coefficients were calculated from the values of D-obs, in the plateau regions by applying a correction for the fraction of water protons present, obtained from the equilibrium water contents of the gels. The true self-diffusion coefficient for water in polyHEMA obtained at 310 K by this method was 5.5 x 10(-10) m(2) s(-1). For the copolymers containing 20% HEMA or more a single value of the self-diffusion coefficient was found, which was somewhat larger than the corresponding values obtained for the macroscopic diffusion coefficient from sorption measurements. For polyTHFMA and copolymers containing less than 20% HEMA, the PFG-NMR stimulated echo attenuation decay curves and the log-attenuation plots were characteristic of the presence of two diffusing water species. The self-diffusion coefficients of water in the equilibrium-hydrated copolymers were found to be dependent on the copolymer composition, decreasing with increasing THFMA content.

Exact theory of sedimentation equilibrium made useful

The exact description of the thermodynamics of solutions has been used to describe, without approximation, the distribution of all the components of an incompressible solution in a centrifuge cell at sedimentation equilibrium. Thermodynamic parameters describing the interactions between solute components of known molar mass can be obtained by direct analysis of the experimental data. Interpretation of the measured thermodynamic parameters in terms of molecular interactions requires that an arbitrary distinction be made between nonassociative forces, like hard-sphere volume-exclusion and mean-field electrostatic repulsion or attraction, and specific short-range forces of association that give rise to the formation of molecular aggregates. Provided the former can be accounted for adequately, the effects of the latter can be elucidated in the form of good estimates of the equilibrium constants for the reactions of aggregation.