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 649 records · Page 36

Ab Initio Studies of Fluorine and Chlorine Oxide and Nitrogen Oxide Species of Interest in Stratospheric Chemistry

The quality of fundamental vibrational frequencies determined using the CCSD(T) method (singles and doubles coupled-cluster theory plus a perturbational estimate of the effects of connected triple excitations) is shown to be very good, usually predicting band centers to within plus or minus 8 per centimeter. This approach is applied to several molecules of interest in atmospheric chemistry, such as HNO, cis-FONO, cis-ClONO, and ClOOH. The HNO molecule displays a large and unusual anharmonicity in the H-N stretch. For the calculation of ultraviolet (UV) spectra, the linear response CCSD (LRCCSD) approach (which is equivalent to EOM-CCSD) has been shown to yield vertical excitation energies that are accurate to approximately equal to 0.1 eV for singly excited electronic states. This method together with more approximate methods is used to examine the UV spectra of several molecules important in stratospheric chemistry, including HOCl, Cl2O, ClONO2, HONO2, ClOOCl, ClOOH, and HOOH.

Lee, Timothy J.↗

Inflow velocity perturbations due to fuselage effects in the presence of a fully interactive wake

The velocity field at locations above the plane of the rotor disk of a generic helicopter configuration has been measured experimentally and modeled using two computational models. The first computational model represents only the portion of the velocity field due to the fuselage, using only the potential due to freestream onset flow. The second model represents all elements of the problem in an incompressible field: the rotor, the rotor wake, and the non-lifting fuselage. The first model is used to assess the fidelity requirement for fuselage paneling. Comparisons are made between a much simplified body shape, an accurate body shape, a body shape with a crude hub/shaft model and wind-tunnel measurements. The second model is used to assess the relative magnitude of the velocity perturbations at the rotor disk due to the effects of the fuselage. A fully interactive time-stepping wake method is used to predict the velocities at the rotor disk. These velocities are computed with and without the fuselage in the model and compared with measured wind-tunnel inflow velocities. The difference between the velocities computed with and without the fuselage is then compared with the perturbation velocity field computed by the first method.

Berry, John D.↗

Stability and noise spectra of relative Loran-C frequency comparisons

Relative comparisons of Loran-C frequency transmissions between the master station of Catanzaro (Simeri Crichi) and the X, Z slave stations of Estartit (Spain) and Lampedusa (Italy) are carrying out by the GG LORSTA monitor station of the Mediterranean Sea Loran-C chain. These comparisons are able to emphasize the relative and, under certain conditions, the absolute rate of the emitting standard frequencies of the slave stations and some relevant statistical properties of the Loran-C Method for frequency transmission and time synchronization. The stability of each Loran-C frequency standard transmission is subject to perturbations, more or less known, due to the propagation medium and other causes. Following the Allan (1966) method for data processing, the performance of the relative rate of frequency of the transmissions of the X, Z slave stations are described calculating the standard deviation of a set of N frequency measurements from its mean averaged during sampling times. This standard deviation is designated as the measure of the stability of the Loran-C frequency transmission.

Proverbio, E.↗

Analysis of crustal structure of Venus utilizing residual Line-of-Sight (LOS) gravity acceleration and surface topography data. A trial of global modeling of Venus gravity field using harmonic spline method

To construct Venus' gravity disturbance field (or gravity anomaly) with the spacecraft-observer line of site (LOS) acceleration perturbation data, both a global and a local approach can be used. The global approach, e.g., spherical harmonic coefficients, and the local approach, e.g., the integral operator method, based on geodetic techniques are generally not the same, so that they must be used separately for mapping long wavelength features and short wavelength features. Harmonic spline, as an interpolation and extrapolation technique, is intrinsically flexible to both global and local mapping of a potential field. Theoretically, it preserves the information of the potential field up to the bound by sampling theorem regardless of whether it is global or local mapping, and is never bothered with truncation errors. The improvement of harmonic spline methodology for global mapping is reported. New basis functions, a singular value decomposition (SVD) based modification to Parker & Shure's numerical procedure, and preliminary results are presented.

Fang, Ming↗

High zonal harmonics of rapidly rotating planets

A new perturbation expansion is derived for the structure of rotating bodies in hydrostatic equilibrium. The method uses an expansion of the density on Legendre polynomial functions of angle, and can be developed analytically in a manner analogous to the standard level-surface perturbation theory. The new theory proceeds from a prescribed pressure-density relation rather than from a prescribed density distribution, and is both simpler and more physically transparent than the level-surface approach. High zonal harmonics are shown to arise via a transfer function involving derivatives of the interior sound velocity, and via mixing of multipole density components in the outer shell of the planet. Sample calculations for polytropic sequences are presented, as well as standard gravity models for Jupiter and Saturn. Mathematical subleties of the theory are discussed in an appendix.

Hubbard, W. B.↗

Rapid and Accurate Computation of Invariant Tori, Manifolds, and Connections Near Mean Motion Resonances in the Periodically Perturbed Planar Circular Restricted 3-Body Problem Models

When the planar circular restricted 3-body problem (RTBP) is periodically perturbed, most unstable resonant periodic orbits become invariant tori. In this study, we 1) develop a quasi-Newton method which simultaneously solves for the tori and their center, stable, and unstable directions; 2) implement continuation by both perturbation as well as rotation numbers; 3) compute Fourier-Taylor parameterizations of the stable and unstable manifolds; 4) globalize these manifolds; 5) compute homoclinic and heteroclinic connections. Our methodology improves on efficiency and accuracy compared to prior studies, and applies to a variety of periodic perturbations. We demonstrate the tools on the planar elliptic RTBP.

de la Llave, Rafael↗

Vibrations in small Mg clusters

The equilibrium geometries, binding energies, and harmonic frequencies of Mg3 and Mg4 have been determined using large atomic natural orbital basis sets in conjunction with high levels of electron correlation. The correlation treatments comprise multireference configuration interaction (MRCI), singles and doubles coupled-cluster (CCSD) theory and the CCSD(T) extension that includes a perturbational estimate of connected triple excitations. As with previous studies of small Be clusters, the CCSD(T) method is found to reproduce the MRCI results with a remarkable degree of accuracy and at a fraction of the computational cost. Using the CCSD(T) method, full quartic force fields for Mg3 and Mg4 have been determined and anharmonic analyses have been performed using second-order perturbation theory. Vibration-rotation interaction constants and centrifugal distortion constants have also been determined.

Lee, Timothy J.↗

Comparison of the quadratic configuration interaction and coupled cluster approaches to electron correlation including the effect of triple excitations

The recently proposed quadratic configuration interaction (QCI) method is compared with the more rigorous coupled cluster (CC) approach for a variety of chemical systems. Some of these systems are well represented by a single-determinant reference function and others are not. The finite order singles and doubles correlation energy, the perturbational triples correlation energy, and a recently devised diagnostic for estimating the importance of multireference effects are considered. The spectroscopic constants of CuH, the equilibrium structure of cis-(NO)2 and the binding energies of Be3, Be4, Mg3, and Mg4 were calculated using both approaches. The diagnostic for estimating multireference character clearly demonstrates that the QCI method becomes less satisfactory than the CC approach as non-dynamical correlation becomes more important, in agreement with a perturbational analysis of the two methods and the numerical estimates of the triple excitation energies they yield. The results for CuH show that the differences between the two methods become more apparent as the chemical systems under investigation becomes more multireference in nature and the QCI results consequently become less reliable. Nonetheless, when the system of interest is dominated by a single reference determinant both QCI and CC give very similar results.

Taylor, Peter R.↗

Numerically-simulated formation and propagation of interplanetary shocks

The present numerical method for simulating the formation and propagation of interplanetary shocks is based on the shock-capturing finite difference scheme of Lax (1950) and Lax and Wendroff (1960), as well as the recent method of NEAR characteristics of Nakagawa (1980, 1981). Attention is given to examples which strongly suggest that all the shocked solar wind plasma parameters due to given physical perturbations, such as flare-generated shocks, can be predicted through the use of this method; the method is, however, limited to the supersonic and super-Alfvenic flow.

Wu, S. T.↗

Computational aspects of sensitivity calculations in linear transient structural analysis

The calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, and transient response problems is studied. Several existing sensitivity calculation methods and two new methods are compared for three example problems. Approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite model. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models. This was found to result in poor convergence of stress sensitivities in several cases. Two semianalytical techniques are developed to overcome this poor convergence. Both new methods result in very good convergence of the stress sensitivities; the computational cost is much less than would result if the vibration modes were recalculated and then used in an overall finite difference method.

Greene, W. H.↗

On the computation of the transonic perturbation flow field around two- and three-dimensional oscillating wings

A finite difference method for solving the unsteady flow about harmonically oscillating wings is investigated. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. Solutions are obtained using relaxation procedures. It is determined that there is a limit on reduced frequency, which is a function of Mach number and size of mesh region, above which the relaxation procedures will not converge. It is found that row line relaxation is more efficient than column relaxation and results are presented for a rectangular wing in harmonic pitch.

Weatherill, W. H.↗

An asymptotic induced numerical method for the convection-diffusion-reaction equation

A parallel algorithm for the efficient solution of a time dependent reaction convection diffusion equation with small parameter on the diffusion term is presented. The method is based on a domain decomposition that is dictated by singular perturbation analysis. The analysis is used to determine regions where certain reduced equations may be solved in place of the full equation. Parallelism is evident at two levels. Domain decomposition provides parallelism at the highest level, and within each domain there is ample opportunity to exploit parallelism. Run time results demonstrate the viability of the method.

Scroggs, Jeffrey S.↗

Orbit determination and prediction, and computer programs

Orbit determination and prediction programs are needed to generate ephemerides for the satellite. Orbit determination is from tracking data consisting of angles only, and is based on a modified version of a method by R. E. Briggs and J. W. Slowey of the Smithsonian Institution. Trends in the data due to perturbations from a Keplerian orbit are removed before this process, and estimates of the orbital elements from individual passes are combined statistically to produce refined estimates. Ephemeris calculation is by a semi-analytic method in which deviations from a Keplerian orbit are obtained by integrating the perturbing forces. The programs to implement these procedures have been written for both the IBM 7090 and the IBM 1620 computers.

EPHEMERIS↗

The Vibrational Frequencies of CaO2, ScO2, and TiO2: A Comparison of Theoretical Methods

The vibrational frequencies of several states of CaO2, ScO2, and TiO2 are computed at using density functional theory (DFT), the Hatree-Fock approach, second order Moller-Plesset perturbation theory (MP2), and the complete-active-space self-consistent-field theory. Three different functionals are used in the DFT calculations, including two hybrid functionals. The coupled cluster singles and doubles approach including the effect of unlinked triples, determined using perturbation theory, is applied to selected states. The Becke-Perdew 86 functional appears to be the cost effective method of choice, although even this functional does not perform well for one state of CaO2. The MP2 approach is significantly inferior to the DFT approaches.

Rosi, Marzio↗

Approximation of expectation values.

Approximation of expectation values for nonenergy properties via perturbation theory, obtaining values for molecular polarizability

POLARIZATION CHARACTERISTICS↗

On the accuracy of Whitham's method

The steady flow of an ideal gas past a conical body is studied by the method of matched asymptotic expansions and by Whitham's method in order to assess the accuracy of the latter. It is found that while Whitham's method does not yield a correct asymptotic representation of the perturbation field to second order in regions where the flow ahead of the Mach cone of the apex is disturbed, it does correctly predict the changes of the second-order perturbation quantities across a shock (the first-order shock strength). The results of the analysis are illustrated by a special case of a flat, rectangular plate at incidence.

Zahalak, G. I.↗