First-principles binary diffusion coefficients for H, H2 and four normal alkanes + N2
- Sandia National Lab. (SNL-CA), Livermore, CA (United States)
- Argonne National Lab. (ANL), Argonne, IL (United States)
Collision integrals related to binary (dilute gas) diffusion are calculated classically for six species colliding with N2. The most detailed calculations make no assumptions regarding the complexity of the potential energy surface, and the resulting classical collision integrals are in excellent agreement with previous semiclassical results for H + N2 and H2 + N2 and with recent experimental results for C n H2n+2 + N2, n = 2–4. The detailed classical results are used to test the accuracy of three simplifying assumptions typically made when calculating collision integrals: (1) approximating the intermolecular potential as isotropic, (2) neglecting the internal structure of the colliders (i.e., neglecting inelasticity), and (3) employing unphysical R–12 repulsive interactions. The effect of anisotropy is found to be negligible for H + N2 and H2 + N2 (in agreement with previous quantum mechanical and semiclassical results for systems involving atomic and diatomic species) but is more significant for larger species at low temperatures. For example, the neglect of anisotropy decreases the diffusion coefficient for butane + N2 by 15% at 300 K. The neglect of inelasticity, in contrast, introduces only very small errors. Approximating the repulsive wall as an unphysical R–12 interaction is a significant source of error at all temperatures for the weakly interacting systems H + N2 and H2 + N2, with errors as large as 40%. For the normal alkanes in N2, which feature stronger interactions, the 12/6 Lennard–Jones approximation is found to be accurate, particularly at temperatures above –700 K where it predicts the full-dimensional result to within 5% (although with somewhat different temperature dependence). Overall, the typical practical approach of assuming isotropic 12/6 Lennard–Jones interactions is confirmed to be suitable for combustion applications except for weakly interacting systems, such as H + N2. For these systems, anisotropy and inelasticity can safely be neglected but a more detailed description of the repulsive wall is required for quantitative predictions. Moreover, a straightforward approach for calculating effective isotropic potentials with realistic repulsive walls is described. An analytic expression for the calculated diffusion coefficient for H + N2 is presented and is estimated to have a 2-sigma error bar of only 0.7%.
- Research Organization:
- Sandia National Lab. (SNL-CA), Livermore, CA (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA)
- Grant/Contract Number:
- AC04-94AL85000
- OSTI ID:
- 1237364
- Report Number(s):
- SAND-2015-0160J; JCPSA6; 560415; TRN: US1600393
- Journal Information:
- Journal of Chemical Physics, Vol. 141, Issue 12; ISSN 0021-9606
- Publisher:
- American Institute of Physics (AIP)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Understanding the breakdown of classic two-phase theory and spray atomization at engine-relevant conditions
|
journal | April 2016 |
Similar Records
Investigations of perturbative curved path computational schemes for the calculation of H/sub 2/O--N/sub 2/ collision broadening. [Dipole-quadrupole interaction, resonance]
Molecule--molecule collisions: Semiclassical calculations of rotational excitation for the H/sub 2/--CO/sub 2/ system using the exponential approximation