Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative methods”

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.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO reversible reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HNO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X double bond H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate (Delta)H(sub f, sup o) value for HNO using state-of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +/- 0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗
NASA NTRS1994

A Coupled-Cluster Study of XON (X=H, F, Cl), and the XON (left and right arrow) XNO Transition States

The XON molecules (X=H, F, and Cl) have been studied using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, CCSD(T), in conjunction with a double polarized triple-zeta basis set. The equilibrium geometries, dipole moments, harmonic vibrational frequencies and infrared intensities have been predicted. The X-O bond distance is shown to be abnormally long for X=F and Cl, and this is attributed to the degree of ionic bonding and the stability of NO(+). Based on Mulliken population analyses, it is shown that there is a significant degree of X(-), ON(+) ionic bonding character for FON and ClON, whereas for HON the ionic character is reduced and best described as H(+), NO(-). The stability of the XON molecule relative to the XNO isomer is shown to increase in the order HON<FON<ClON, although even ClON is 27.0 plus or minus 1.0 kcal/mol (0 K) less stable than ClNO. The XON (left and right arrow) XNO transition states are also investigated. FON and ClON possess the lowest barrier heights to isomerization (7.2 plus or minus 1.0 kcal/mol at 0 K), but these are sufficiently large to suggest that isomerization should not be rapid at low temperatures.

Lee, Timothy↗
NASA NTRS1996

On the Reaction of FNO2 with CH3, t-butyl, and C13H21

Theoretical studies are reported for the reaction of FN02 with the radicals CH3, t-butyl, and C13H21, which are templates for the radical site of a hydrogenated diamond surface. All structures axe fully optimized using density functional theory (DFT) based on the B3LYP functional. Calibration calculations axe performed for CH3 + FNO2 using the coupled cluster approach, the the internally contracted multireference configuration. interaction method, and second order perturbation theory based upon complete active space SCF (CASSCF) reference wave function. These calibration calculations support the B3LYP approach for the calculation of bond energies, but show the B3LYP barrier is too low. Combining the calibration calculations with the larger clusters yields our best estimate of a barrier of about 10 kcal/mol. for the reaction of FNO2 with a radical site on hydrogenated diamond.

Thuemmel, H. T.↗
NASA NTRS2002

Gas-Phase Reaction Pathways and Rate Coefficients for the Dichlorosilane-Hydrogen and Trichlorosilane-Hydrogen Systems

As part of NASA Ames Research Center's Integrated Process Team on Device/Process Modeling and Nanotechnology our goal is to create/contribute to a gas-phase chemical database for use in modeling microelectronics devices. In particular, we use ab initio methods to determine chemical reaction pathways and to evaluate reaction rate coefficients. Our initial studies concern reactions involved in the dichlorosilane-hydrogen (SiCl2H2--H2) and trichlorosilane-hydrogen (SiCl2H-H2) systems. Reactant, saddle point (transition state), and product geometries and their vibrational harmonic frequencies are determined using the complete-active-space self-consistent-field (CASSCF) electronic structure method with the correlation consistent polarized valence double-zeta basis set (cc-pVDZ). Reaction pathways are constructed by following the imaginary frequency mode of the saddle point to both the reactant and product. Accurate energetics are determined using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations (CCSD(T)) extrapolated to the complete basis set limit. Using the data from the electronic structure calculations, reaction rate coefficients are obtained using conventional and variational transition state and RRKM theories.

Dateo, Christopher E.↗
NASA NTRS2002

Independent Confirmation of the Pioneer 10 Anomalous Acceleration

I perform an independent analysis of radio Doppler tracking data from the Pioneer 10 spacecraft for the time period 1987-1994. All of the tracking data were taken from public archive sources, and the analysis tools were developed independently by myself. I confirm that an apparent anomalous acceleration is acting on the Pioneer 10 spacecraft, which is not accounted for by present physical models of spacecraft navigation. My best fit value for the acceleration, including corrections for systematic biases and uncertainties, is (8.60 plus or minus 1.34) x 10(exp -8) centimeters per second, directed towards the Sun. This value compares favorably to previous results. I examine the robustness of my result to various perturbations of the analysis method, and find agreement to within plus or minus 5%. The anomalous acceleration is reasonably constant with time, with a characteristic variation time scale of greater than 70 yr. Such a variation timescale is still too short to rule out on-board thermal radiation effects, based on this particular Pioneer 10 data set.

Markwardt, Craig B.↗
NASA NTRS2013

Good Practices in Free-energy Calculations

As access to computational resources continues to increase, free-energy calculations have emerged as a powerful tool that can play a predictive role in drug design. Yet, in a number of instances, the reliability of these calculations can be improved significantly if a number of precepts, or good practices are followed. For the most part, the theory upon which these good practices rely has been known for many years, but often overlooked, or simply ignored. In other cases, the theoretical developments are too recent for their potential to be fully grasped and merged into popular platforms for the computation of free-energy differences. The current best practices for carrying out free-energy calculations will be reviewed demonstrating that, at little to no additional cost, free-energy estimates could be markedly improved and bounded by meaningful error estimates. In energy perturbation and nonequilibrium work methods, monitoring the probability distributions that underlie the transformation between the states of interest, performing the calculation bidirectionally, stratifying the reaction pathway and choosing the most appropriate paradigms and algorithms for transforming between states offer significant gains in both accuracy and precision. In thermodynamic integration and probability distribution (histogramming) methods, properly designed adaptive techniques yield nearly uniform sampling of the relevant degrees of freedom and, by doing so, could markedly improve efficiency and accuracy of free energy calculations without incurring any additional computational expense.

Pohorille, Andrew↗
NASA NTRS2014

Vibrational Frequencies and Spectroscopic Constants for 1(sup 3)A' HNC and 1(sup 3)A' HOC+ from High-Accuracy Quartic Force Fields

The spectroscopic constants and vibrational frequencies for the 1(sup 3)A' states of HNC, DNC, HOC+, and DOC+ are computed and discussed in this work. The reliable CcCR quartic force field based on high-level coupled cluster ab initio quantum chemical computations is exclusively utilized to provide the anharmonic potential. Then, second order vibrational perturbation theory and vibrational configuration interaction methods are employed to treat the nuclear Schroedinger equation. Second-order perturbation theory is also employed to provide spectroscopic data for all molecules examined. The relationship between these molecules and the corresponding 1(sup 3)A' HCN and HCO+ isomers is further developed here. These data are applicable to laboratory studies involving formation of HNC and HOC+ as well as astronomical observations of chemically active astrophysical environments.

Astrochemistry↗
NASA NTRS2017

Role of Forcing Uncertainty and Background Model Error Characterization in Snow Data Assimilation

Accurate specification of the model error covariances in data assimilation systems is a challenging issue. Ensemble land data assimilation methods rely on stochastic perturbations of input forcing and model prognostic fields for developing representations of input model error covariances. This article examines the limitations of using a single forcing dataset for specifying forcing uncertainty inputs for assimilating snow depth retrievals. Using an idealized data assimilation experiment, the article demonstrates that the use of hybrid forcing input strategies (either through the use of an ensemble of forcing products or through the added use of the forcing climatology) provide a better characterization of the background model error, which leads to improved data assimilation results, especially during the snow accumulation and melt-time periods. The use of hybrid forcing ensembles is then employed for assimilating snow depth retrievals from the AMSR2 (Advanced Microwave Scanning Radiometer 2) instrument over two domains in the continental USA with different snow evolution characteristics. Over a region near the Great Lakes, where the snow evolution tends to be ephemeral, the use of hybrid forcing ensembles provides significant improvements relative to the use of a single forcing dataset. Over the Colorado headwaters characterized by large snow accumulation, the impact of using the forcing ensemble is less prominent and is largely limited to the snow transition time periods. The results of the article demonstrate that improving the background model error through the use of a forcing ensemble enables the assimilation system to better incorporate the observational information.

assimilation↗
NASA NTRS2017

Satellite Constellation Orbit Design to Enable a Space-Based Radio Interferometer

Two different design methods for a networked constellation of N small satellites are presented. In the first method, the (linear) Clohessy-Wilthsire equations are used as an initial design tool, followed by conversion to a two-body model. Discrepancies between the linear and nonlinear solutions are minimized in the conversion process. The second method utilizes invariant manifold theory, by perturbing a reference trajectory in different directions along the center eigenvectors. Both methods require 5N parameters to fully define a constellation. In a relative, rotating frame of the reference path, the spacecraft appear as periodic ellipses of varying sizes. Deployment, reconfiguration using propulsive maneuvers, and station keeping costs for an example mission scenario are addressed.

Chien, Steve A.↗
NASA NTRS2019

Highly Accurate Quartic Force Field and Rovibrational Spectroscopic Constants for the Azirinyl Cation (c-C2NH2+) and Its Isomers

Azirinyl cation is an aromatic cyclic molecule that is isoelectronic with cyclopropenylidene, c-C3H2, and c-C3H3+. Cyclopropenylidene has been shown to be ubiquitous, existing in many different astrophysical environments. Given the similar chemistry between C and N, and the relative abundances between C and N in astrophysical environments, it is expected that there should be aromatic ringed molecules that incorporate N in the ring, but as yet, no such molecule has been identified. To address this issue, the present study uses high levels of electronic structure theory to compute a highly accurate quartic force field (QFF) for azirinyl cation and its two lowest lying isomers, the cyanomethyl and isocyanomethyl cations. The theoretical approach uses the singles and doubles coupled-cluster method that includes a perturbative correction for connected triple excitations, CCSD(T), together with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core-correlation. The QFF is then used in a second order vibrational perturbation theory analysis (VPT2) to compute the fundamental vibrational frequencies and rovibrational spectroscopic constants for all three C2NH2+ isomers. The reliability of the VPT2 vibrational frequencies is tested by comparison to vibrational configuration interaction (VCI) calculations and excellent agreement is found between the two approaches. Fundamental vibrational frequencies and rovibrational spectroscopic constants for all singly substituted 13C, 15N, and D isotopologues are also reported. It is expected that the highly accurate spectroscopic data reported herein will be useful in the identification of these cations in high resolution experimental or astronomical observations.

Partha P Bera↗
NASA NTRS1991

Resonant frequency calculations using a hybrid perturbation-Galerkin technique

A two-step hybrid perturbation Galerkin technique is applied to the problem of determining the resonant frequencies of one or several degree of freedom nonlinear systems involving a parameter. In one step, the Lindstedt-Poincare method is used to determine perturbation solutions which are formally valid about one or more special values of the parameter (e.g., for large or small values of the parameter). In step two, a subset of the perturbation coordinate functions determined in step one is used in Galerkin type approximation. The technique is illustrated for several one degree of freedom systems, including the Duffing and van der Pol oscillators, as well as for the compound pendulum. For all of the examples considered, it is shown that the frequencies obtained by the hybrid technique using only a few terms from the perturbation solutions are significantly more accurate than the perturbation results on which they are based, and they compare very well with frequencies obtained by purely numerical methods.

Geer, James F.↗
NASA NTRS1991

Resonant frequency calculations using a hybrid perturbation-Galerkin technique

A two-step hybrid perturbation Galerkin technique is applied to the problem of determining the resonant frequencies of one or several degrees of freedom nonlinear systems involving a parameter. In one step, the Lindstedt-Poincare method is used to determine perturbation solutions which are formally valid about one or more special values of the parameter (e.g., for large or small values of the parameter). In step two, a subset of the perturbation coordinate functions determined in step one is used in Galerkin type approximation. The technique is illustrated for several one degree of freedom systems, including the Duffing and van der Pol oscillators, as well as for the compound pendulum. For all of the examples considered, it is shown that the frequencies obtained by the hybrid technique using only a few terms from the perturbation solutions are significantly more accurate than the perturbation results on which they are based, and they compare very well with frequencies obtained by purely numerical methods.

Geer, James F.↗
NASA NTRS1998

Computational Investigation of Experimental Interaction Impedance Obtained by Perturbation for Helical Traveling Wave Tube Structures

Conventional methods used to measure the cold- test interaction impedance of helical slow-wave structures involve perturbing a helical circuit with a cylindrical dielectric rod placed on the central axis of the circuit It has been shown that the difference in resonant frequency or axial phase shift between the perturbed and unperturbed circuits can be related to the interaction impedance. However, because of the complex configuration of the helical circuit, deriving this relationship involves several approximations. With the advent of accurate three-dimensional (3-D) helical circuit models, these standard approximations can be fully Investigated. This paper addresses the most prominent approximations made in the analysis for measured interaction impedance by Lagerstrom and investigates their accuracy using the 3-D simulation code MAFIA. It is shown that a more accurate value of interaction impedance can be obtained by using 3-D computational methods rather than performing costly and time consuming experimental cold-test measurements.

Kory, Carol L.↗
NASA NTRS1973

Effects of resonant tesseral gravity coefficients on Viking-type orbits.

A special perturbation technique has been developed for the long-period and secular motion of a Viking-type orbiter. The technique uses a method of singly averaging the perturbations over the mean anomaly to determine changes in orbital elements due to atmospheric drag, solar radiation pressure, solar gravity, and the asymmetrical Mars gravitational field. The technique has been specialized to treat resonant tesseral gravity effects resulting from the commensurability of the satellite orbit period with the rotational period of the primary. Results are presented which indicate that orbiter lifetimes can be significantly reduced by resonant tesserals for orbits near the critical inclination. Conversely, it is shown that lifetimes can be extended via capture in resonant orbits which are stable against small drag forces.

Compton, H. R.↗
NASA NTRS1980

Evolution of solar magnetic fields - A new approach to MHD initial-boundary value problems by the method of nearcharacteristics

A method of analysis for the MHD initial-boundary problem is presented in which the model's formulation is based on the method of nearcharacteristics developed by Werner (1968) and modified by Shin and Kot (1978). With this method, the physical causality relationship can be traced from the perturbation to the response as in the method of characteristics, while achieving the advantage of a considerable reduction in mathematical procedures. The method offers the advantage of examining not only the evolution of nonforce free fields, but also the changes of physical conditions in the atmosphere accompanying the evolution of magnetic fields. The physical validity of the method is demonstrated with examples, and their significance in interpreting observations is discussed.

Nakagawa, Y.↗