A set of Direct Simulation Monte Carlo (DSMC) chemical-reaction models recently proposed by Bird and based solely on the collision energy and the vibrational energy levels of the species involved is applied to calculate nonequilibrium chemical-reaction rates for atmospheric reactions in hypersonic flows. The DSMC non-equilibrium model predictions are in good agreement with theoretical models and experimental measurements. The observed agreement provides strong evidence that modeling chemical reactions using only the collision energy and the vibrational energy levels provides an accurate method for predicting non-equilibrium chemical-reaction rates.
This report summarizes the work completed under the Laboratory Directed Research and Development (LDRD) project 09-1351, 'Computational Investigation of Thermal Gas Separation for CO{sub 2} Capture'. Thermal gas separation for a binary mixture of carbon dioxide and nitrogen is investigated using the Direct Simulation Monte Carlo (DSMC) method of molecular gas dynamics. Molecular models for nitrogen and carbon dioxide are developed, implemented, compared to theoretical results, and compared to several experimental thermophysical properties. The molecular models include three translational modes, two fully excited rotational modes, and vibrational modes, whose degree of excitation depends on the temperature. Nitrogen has one vibrational mode, and carbon dioxide has four vibrational modes (two of which are degenerate). These models are used to perform a parameter study for mixtures of carbon dioxide and nitrogen confined between parallel walls over realistic ranges of gas temperatures and nominal concentrations of carbon dioxide. The degree of thermal separation predicted by DSMC is slightly higher than experimental values and is sensitive to the details of the molecular models.
This report summarizes the work completed during FY2009 for the LDRD project 09-1332 'Molecule-Based Approach for Computing Chemical-Reaction Rates in Upper-Atmosphere Hypersonic Flows'. The goal of this project was to apply a recently proposed approach for the Direct Simulation Monte Carlo (DSMC) method to calculate chemical-reaction rates for high-temperature atmospheric species. The new DSMC model reproduces measured equilibrium reaction rates without using any macroscopic reaction-rate information. Since it uses only molecular properties, the new model is inherently able to predict reaction rates for arbitrary nonequilibrium conditions. DSMC non-equilibrium reaction rates are compared to Park's phenomenological non-equilibrium reaction-rate model, the predominant model for hypersonic-flow-field calculations. For near-equilibrium conditions, Park's model is in good agreement with the DSMC-calculated reaction rates. For far-from-equilibrium conditions, corresponding to a typical shock layer, the difference between the two models can exceed 10 orders of magnitude. The DSMC predictions are also found to be in very good agreement with measured and calculated non-equilibrium reaction rates. Extensions of the model to reactions typically found in combustion flows and ionizing reactions are also found to be in very good agreement with available measurements, offering strong evidence that this is a viable and reliable technique to predict chemical reaction rates.
Particle contamination is analyzed for a reticle in the inner pod of a carrier with particular emphasis on the effect of raising the cover of the inner pod before removing the reticle from the carrier at atmospheric pressure (not low pressure). Two mechanisms for particle transport into the gap between the base plate and the reticle are considered: injection and advection-diffusion. It is shown that injection is not an important mechanism but that advection-diffusion transport can carry particles deeply into the gap, where they can deposit on the reticle surface. Closed-form expressions are presented for the transmission probability that particles at the reticle edge are transported inward past the exclusion zone around the reticle perimeter. The gas flow in the gap that occurs during cover-raising is found by numerical simulation, and the closed-form expressions are applied to determine the probability of contamination for different cover-raising scenarios.
The convergence rate of a new direct simulation Monte Carlo (DSMC) method, termed 'sophisticated DSMC', is investigated for one-dimensional Fourier flow. An argon-like hard-sphere gas at 273.15K and 266.644Pa is confined between two parallel, fully accommodating walls 1mm apart that have unequal temperatures. The simulations are performed using a one-dimensional implementation of the sophisticated DSMC algorithm. In harmony with previous work, the primary convergence metric studied is the ratio of the DSMC-calculated thermal conductivity to its corresponding infinite-approximation Chapman-Enskog theoretical value. As discretization errors are reduced, the sophisticated DSMC algorithm is shown to approach the theoretical values to high precision. The convergence behavior of sophisticated DSMC is compared to that of original DSMC. The convergence of the new algorithm in a three-dimensional implementation is also characterized. Implementations using transient adaptive sub-cells and virtual sub-cells are compared. The new algorithm is shown to significantly reduce the computational resources required for a DSMC simulation to achieve a particular level of accuracy, thus improving the efficiency of the method by a factor of 2.
This report documents technical work performed to complete the ASC Level 2 Milestone 2841: validation of thermal models for a prototypical MEMS thermal actuator. This effort requires completion of the following task: the comparison between calculated and measured temperature profiles of a heated stationary microbeam in air. Such heated microbeams are prototypical structures in virtually all electrically driven microscale thermal actuators. This task is divided into four major subtasks. (1) Perform validation experiments on prototypical heated stationary microbeams in which material properties such as thermal conductivity and electrical resistivity are measured if not known and temperature profiles along the beams are measured as a function of electrical power and gas pressure. (2) Develop a noncontinuum gas-phase heat-transfer model for typical MEMS situations including effects such as temperature discontinuities at gas-solid interfaces across which heat is flowing, and incorporate this model into the ASC FEM heat-conduction code Calore to enable it to simulate these effects with good accuracy. (3) Develop a noncontinuum solid-phase heat transfer model for typical MEMS situations including an effective thermal conductivity that depends on device geometry and grain size, and incorporate this model into the FEM heat-conduction code Calore to enable it to simulate these effects with good accuracy. (4) Perform combined gas-solid heat-transfer simulations using Calore with these models for the experimentally investigated devices, and compare simulation and experimental temperature profiles to assess model accuracy. These subtasks have been completed successfully, thereby completing the milestone task. Model and experimental temperature profiles are found to be in reasonable agreement for all cases examined. Modest systematic differences appear to be related to uncertainties in the geometric dimensions of the test structures and in the thermal conductivity of the polycrystalline silicon test structures, as well as uncontrolled nonuniform changes in this quantity over time and during operation.