SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The angular dependence of the interaction between two rigid linear rotors is considered, and one of the two expansions in current use is shown to be inappropriate. Interpolation of tabulated radial interactions is also discussed, and numerical tests are presented which demonstrate the accuracy of Lagrange and spline algorithms as a function of the spacing between points.
Theoretical potentials for rotational excitation of H2O by He were constructed via several methods, all of which start with a large basis set SCF interaction. The semiempirical Hartree-Fock with damped dispersion model adds a damped long-range attraction with parameters adjusted to fit experimental total differential cross sections. Purely ab initio potentials add correlation energies obtained via perturbation theory (MP2 and MP4) or a variational method (ICF1). Scattering calculations were performed on all surfaces to compare wih available beam scattering and pressure broadening data and to assess sensitivity of state-to-state rates to uncertainties in the potential. From comparison with the limited experimental data, the ICF1 surface appears to be marginally better than the MP4 surface. Thermal rates calculated from this surface should be accurate to better than 50 percent, at least for the larger, more important rates.
The results of SCF-CI calculations of the intermolecular potential surface for CO and He are compared with pressure broadening data. Coupled states calculations were performed for both the electron gas and the more accurate SCF-CI potentials at collision energies from 15 to 1200/cm, and Boltzmann averages over collision energy were taken to compare with experimental results at liquid nitrogen (77 K), dry ice (195 K) and room (295 K) temperatures. The cross sections from the SCF-CI potential are shown to be in slightly better agreement with experiment than those from the electron gas potential, although both are consistent with the data to within experimental uncertainties. Closer examination of the experimental results reveals that this essential equivalence is due not to any inherent insensitivity of the pressure broadening technique, but to the uncertainties in the experimental values. It is thus concluded that more accurate experimental values are required to place meaningful constraints on the intermolecular potential.
The paper describes a practical application of Remler's (1971) method by which one constructs a set of phase shifts from high resolution measurements of the differential elastic scattering of protons by rare-gas atoms. These JWKB phase shifts are then formally inverted to determine the corresponding intermolecular potentials. The validity of the method is demonstrated by comparing an intermolecular potential obtained by direct inversion of experimental data with a fairly accurate calculation by Wolniewicz (1965).
Reduced second virial gas coefficients calculations by modified Lennard-Jones intermolecular potential function
In this work, Complex Robert-Bonamy calculations of half-widths and line shifts were done for N2-broadening of water for 1639 transitions in the rotational band using two models for the trajectories. The first is a model correct to second order in time, the Robert-Bonamy parabolic approximation. The second is the solution of Hamilton's equations. Both models use the isotropic part of the atom-atom potential to determine the trajectories. The present calculations used an intermolecular potential expanded to 20th order to assure the convergence of the half-widths and line shifts. The aim of the study is to assess if the difference in the half-widths and line shifts determined from the two trajectory models is greater than the accuracy requirements of the spectroscopic and remote sensing communities. The results of the calculations are compared with measurements of the half-widths and line shifts. It is shown that the effects of the trajectory model greatly exceed the needs of current remote sensing measurements and that line shape parameters calculated using trajectories determined by solving Hamilton's equations agree better with measurement.
Transport collision integrals and thermal diffusion factors evaluated, obeying 9-6 and 28-7 intermolecular potentials
The interatomic V(r) found by potential-model calculations is compared to an ab initio calculation of the intermolecular potential for NeH+ due to Peyerimhoff (1965). The reason for a rather significant difference in the two results is discussed, and a new method due to Remler (1971) involving Regge poles is applied. The method can be used to calculate the differential cross section, where the starting point for such calculations is based on intuition in light of the classical deflection function.
Cross sections for rotational excitation and spectral pressure broadening of HD, HCl, CO, and HCN due to collisions with low energy He atoms have been computed within the 'decoupled l-dominant' (DLD) approximation and are compared with accurate close coupling results and also with two similar approximations, the effective potential of Rabitz and the coupled states of McGuire and Kouri. DLD predictions of state-to-state cross sections are rather good, being only slightly less accurate than coupled states results. DLD is far superior to either the coupled states or effective potential methods for pressure broadening calculations, although it may not be uniformly of the quantitative accuracy desirable for obtaining intermolecular potentials from experimental data.
Several aspects of the kinetic models used in the collision phase of Monte Carlo direct simulations have been studied. Accurate molecular velocity distribution function predictions require a significantly increased number of computational cells in one maximum slope shock thickness, compared to predictions of macroscopic properties. The shape of the highly repulsive portion of the interatomic potential for argon is not well modeled by conventional interatomic potentials; this portion of the potential controls high Mach number shock thickness predictions, indicating that the specification of the energetic repulsive portion of interatomic or intermolecular potentials must be chosen with care for correct modeling of nonequilibrium flows at high temperatures. It has been shown for inverse power potentials that the assumption of variable hard sphere scattering provides accurate predictions of the macroscopic properties in shock waves, by comparison with simulations in which differential scattering is employed in the collision phase. On the other hand, velocity distribution functions are not well predicted by the variable hard sphere scattering model for softer potentials at higher Mach numbers.
The formalism for describing rotational excitation in collisions between symmetric top rigid rotors and spherical atoms is presented both within the accurate quantum close coupling framework and also the coupled states approximation of McGuire and Kouri and the effective potential approximation of Rabitz. Calculations are reported for thermal energy NH3-He collisions, treating NH3 as a rigid rotor and employing a uniform electron gas (Gordon-Kim) approximation for the intermolecular potential. Coupled states are found to be in nearly quantitative agreement with close coupling results while the effective potential method is found to be at least qualitatively correct. Modifications necessary to treat the inversion motion in NH3 are discussed.
While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].
Recent studies of low-Z molecular materials including hydrogen to multimegabar pressures (less than 300 GPa) have uncovered a range of phenomena relevant to understanding the nature of the interiors of the outer planets and their satellites. Synchrotron x ray diffraction measurements (to 42 GPa) have been used to determine the crystal structure of the solid (hexagonal-close packed) and equation of state. Sound velocities in fluid and solid hydrogen (to 24 GPa) have been inverted to obtain elastic constants and aggregate bulk and shear moduli. In addition, an improved intermolecular potential has been determined which fits both static and shock-wave data. Use of the new potential for the molecular envelope of Jupiter suggests the need for major revisions of existing Jovian models or a reanalysis of reported free oscillations for the planet. Studies at higher pressures (greater than 100 GPa) reveal a sequence of pressure-induced symmetry-breaking transitions in molecular hydrogen, giving rise to three high-pressure phases (1, 2, and 3). Phase 1 is the rotationally disordered hcp phase which persists from low pressure to well above 100 GPa at high temperature (e.g., 300 K). Phase 2 is a low-temperature, high-pressure phase (transition at 100 GPa and 77 K in H2) with spectral features indicative of partial rotational ordering and crystallographic distortion. The transition to Phase 3 at 150 GPa is accompanied by a weakening of the molecular bond, gradual changes in orientational ordering, strong enhancement of the infrared intramolecular vibrational absorption, and strong intermolecular interactions similar to those of ambient-pressure network solids. Studies of the phase diagram reveal a triple point near 130 K and 160 GPa. Higher pressure measurements of vibrational spectra place a lower bound of approximately 250 GPa on the predicted transition pressure for dissociation of molecular hydrogen to form a monatomic metal.
Physics-based modeling of hypersonic flows is predicated on the availability of chemical reaction rate coefficients and cross sections for the collisional processes. This approach has been built around the use of quantum mechanical calculations to describe the interaction between the colliding particles. In this approach a potential energy surface (PES) is computed by solving the electronic Schrödinger equation and collision cross sections are determined for that PES using classical, semiclassical or quantum mechanical scattering methods. The rate coefficients are computed by integrating the thermally weighted cross sections. State-to-state rate coefficients are determined by only integrating over a thermal distribution of collisional energies. Finally, thermal rate coefficients are determined by summation of the state-to-state rate coefficients for reactions of molecules in all relevant ro-vibrational energy levels. If the flow is in thermal non-equilibrium, the translational, vibrational and rotational energy modes can be represented in different ways: three unique temperatures can be used to describe the distributions, the populations of individual ro-vibrational energy levels can be determined by solving the Master Equation, or through the use of direct simulation in particle-based Monte Carlo sampling. The PES-to-rate coefficient approach had been proposed and attempted in the early days of digital computing, but it is only in the last 15 years that computer hardware and software have been up to the task of calculating accurate interatomic and intermolecular potentials.
Intermolecular potential energy functions, heat conduction in chemically reacting gases, and heat conductivity of polyatomic gases
Physical adsorption isotherm based on triangular site model, including neighbor interactions via Lennard-Jones intermolecular potentials
High temperature polyatomic gas and gas mixture transport properties, discussing thermoconductivity, multiple intermolecular potentials and plasmas
Benzene quadrupole moment estimates from second virial coefficient data with polarizability in intermolecular potential function and from molecular susceptibility anisotropy