Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Second-order”

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 1,441 records · Page 80

Some Numerical Experiments on Detonation Wave Propagation

In this paper we present the results of a series of numerical experiments done on the propagation and initiation of a detonation wave. The calculations are performed in one-dimension, with considerable grid resolution. Of particular interest are the following questions: (1) the nature of periodic and chaotic instabilities generated by the wave; (2) the influence of the grid resolution on these instabilities; (3) the influence of the 'quality' of the numerical scheme; and (4) the influence of 'noise'. In the calculations, we use a second-order Total Variation Diminishing (TVD) scheme as the basic numerical method, with grid spacings as low as a fraction of a micron. Detonations waves are generated at the closed end of a tube, and allowed to propagate for approximately 20 cm. The required energy for successful initiation of the detonation will be measured for different cases of grid resolution and numerical schemes. A modified version of the TVD scheme has also been devised, which allows for much lower numerical diffusion of the radical species in the exponentially growing region behind the shock. The effect of this modification will be demonstrated. Oscillations in peak pressure and induction length are seen to develop in some cases: the oscillations can go through a sequence of modes, from a regular, high frequency mode to a low frequency mode with period doubling. A chaotic regime can also be obtained. General conclusions on the quality of algorithms will be presented. We will also discuss the performance of a version of the code developed on the IBM SP2 parallel computer.

Cambier, Jean-Luc↗

Capturing the Value: Earth Applications of Space Human Factors Research

This paper details how the Space Human Factors/Life Sciences program at Ames Research Center (ARC) has provided, and continues to provide, a variety of Earth-based benefits. These benefits will be considered under five categories: aeronautics, space-like environments, general applications, human/automation interaction, and methodology. The human factors work at ARC includes a range of activities whose products serve the aerospace community. Some areas of research focus specifically on aeronautical requirements; others are driven by space needs. However, the symbiosis between these two domains allows a sharing of resources, and the insights and experimental results gathered in one domain can often be applied in the other. Aeronautics is an industry whose survival is generally viewed as critical to American competitiveness, and where benefits can result in a very high payoff. The ability to apply space-initiated research to aeronautical requirements represents one example of bringing space benefits down to Earth. The second-order value of space human factors research goes well beyond the aerospace community. Spaceflight shares with a number of other activities certain environmental characteristics that drive human factors engineering design and procedural specification. Spaceflight is an isolated activity, conducted under severely confined conditions, with a high level of risk, and where provisions are restricted and opportunities for outside help are limited. A number of Earth-based activities including submarines and other naval vessels, oil rigs, remote weather stations, and scientific and polar expeditions, share many of these characteristics. These activities serve as testbeds for space-related research and, in turn, space-related research provides beneficial insight to the conduct of these activities.

Connors, Mary M.↗

Numerical Experiments on Homogeneous Strained Turbulence Subjected to Coriolis Force

Homogeneous turbulent flows with various combinations of strain-rate, rotation rate and coriolis force capture some important aspects of more complex flows with streamline curvature and rotation. Presently, a situation is considered in which as a box of turbulence rotates, strain axes rotate with it. This is to be contrasted with the elliptic streamline flow in which the strain axes are fixed in an inertial frame. The elliptic flow is known to exhibit (inviscid) growth of turbulent energy and one might expect even more rapid growth with the strain-axes following the box. Instead, it is found that the sign of the Reynolds shear stress is reversed leading to a negative production term for turbulent energy. Partial understanding of the phenomenon is obtained from a consideration of the rotation of inertial waves relative to the strain axes as well as the "pressure-less" RDT argument put forward by Cambon etal. [J. Fluid Mech, 278, 175]. Some comparisons with the predictions of second-order closure models will be presented.

Shariff, K.↗

A Comparison of the Accuracy of Different Functionals

The errors in the computed geometries, zero-point energies, and atomization energies of molecules containing only first and second row atoms are compared for several levels of theory, including Hartree-Fock, second-order Moller-Plesset perturbation theory (MP2), and density functional theory (DFT) using five different functionals, including two hybrid functionals. The 6-31G* and 6-311+G(3df,2p) basis set are used. Overall, the B3LYP hybrid functional yields the best results.

Bauschlicher, Charles W., Jr.↗

The Calculation of Accurate Harmonic Frequencies of Large Molecules: The Polycyclic Aromatic Hydrocarbons, a Case Study

The vibrational frequencies and infrared intensities of naphthalene neutral and cation are studied at the self-consistent-field (SCF), second-order Moller-Plesset (MP2), and density functional theory (DFT) levels using a variety of one-particle basis sets. Very accurate frequencies can be obtained at the DFT level in conjunction with large basis sets if they are scaled with two factors, one for the C-H stretches and a second for all other modes. We also find remarkably good agreement at the B3LYP/4-31G level using only one scale factor. Unlike the neutral PAHs where all methods do reasonably well for the intensities, only the DFT results are accurate for the PAH cations. The failure of the SCF and MP2 methods is caused by symmetry breaking and an inability to describe charge delocalization. We present several interesting cases of symmetry breaking in this study. An assessment is made as to whether an ensemble of PAH neutrals or cations could account for the unidentified infrared bands observed in many astronomical sources.

Bauschlicher, Charles W., Jr.↗

Positive Anharmonicities: The Oxonide Anion as an Example

An accurate ab initio quartic force field for the ozonide anion has been determined at the singles and doubles coupled-cluster level of theory that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), using the augmented valence triple-zeta correlation consistent one-particle basis set. Convergence of the harmonic frequencies with respect to the one-particle basis set has been demonstrated by computing quadratic force fields at the CCSD(T) level using augmented valence double-zeta and quadruple-zeta basis sets. Fundamental vibrational frequencies have been determined via second-order ro-vibrational perturbation theory and also using exact variational methods. Agreement is very good, and they both show that the antisymmetric O-O stretch, v 3, possesses a positive anharmonic correction (that is, the fundamental frequency is predicted to be higher in energy than the harmonic frequency). Comparison of the O_3 and O3 quartic force fields shows that the positive anharmonic correction is the result of a particularly large F3333 symmetry internal coordinate force constant. The reasonableness of this force constant has been tested by computing both F33 and F3333 at the CCSD, CCSD(T), and CASPT2 levels of theory. A discussion of known positive anharmonicities for stretching modes is presented.

Lee, Timothy J.↗

The General Equations of Equilibrium of Rotationally Symmetric Membranes and some Static Solutions for Uniform Centrifugal Loading

The general equations of equilibrium of a rotationally symmetrical body undergoing large arbitrary deflections but small strains are derived by means of the principle of virtual work. All deflections are referred to a reference surface which, in general, may be taken to coincide with any intermediate state of the membrane. Static solutions of the general equations for steadily rotating membranes are considered for the cases in which the final or the initial configuration of the membrane is specified. For problems in which the final configuration is specified, conditions are given under which such a configuration is obtainable from a stress-free initial state. Such conditions being satisfied, it is shown that the stresses and radial displacements are determinate while axial displacements may be obtained by quadrature. For problems in which the initial configuration is given, the determination of the stresses and deflections is reduced to the solution of a second-order nonlinear differential equation for a stress function. An approximate method of solution based on boundary-layer considerations is presented for the case of a spinning sphere of constant thickness. Specific examples are given for both classes of problems.

Simmonds, James G.↗

An ab initio benchmark study of the H + CO --> HCO reaction

The H + CO --> HCO reaction has been characterized with correlation consistent basis sets at five levels of theory in order to benchmark the sensitivities of the barrier height and reaction ergicity to the one-electron and n-electron expansions of the electronic wave function. Single and multireference methods are compared and contrasted. The coupled cluster method RCCSD(T) was found to be in very good agreement with Davidson-corrected internally-contracted multireference configuration interaction (MRCI+Q). Second-order Moller-Plesset perturbation theory (MP2) was also employed. The estimated complete basis set (CBS) limits for the barrier height (in kcal/mol) for the five methods, including harmonic zero-point energy corrections, are MP2, 4.66; RCCSD, 4.78; RCCSD(T), 4.15; MRCI, 5.10; and MRCI+Q, 4.07. Similarly, the estimated CBS limits for the ergicity of the reaction are: MP2, -17.99; RCCSD, -13.34; RCCSD(T), -13.79; MRCI, -11.46; and MRCI+Q, -13.70. Additional basis set explorations for the RCCSD(T) method demonstrate that aug-cc-pVTZ sets, even with some functions removed, are sufficient to reproduce the CBS limits to within 0.1-0.3 kcal/mol.

NASA Discipline Exobiology↗

Air oxidation of hydrazine. 1. Reaction kinetics on natural kaolinites, halloysites, and model substituent layers with varying iron and titanium oxide and O- center contents

Air oxidation of hydrazine was studied by using a group of kaolinites, halloysites, and substituent oxides as models for the tetrahedral and octahedral sheets. The rate was found to be linear with oxygen. The stoichiometry showed that oxygen was the primary oxidant and that dinitrogen was the only important nitrogen-containing product. The rates on kaolinites were strongly inhibited by water. Those on three-dimensional silica and gibbsite appeared not to be. That on a supposedly layered silica formed from a natural kaolinite by acid leaching showed transitional behavior--slowed relative to that expected from a second-order reaction relative to that on the gibbsite and silica but faster than those on the kaolinites. The most striking result of the reaction was the marked increase in the rate of reaction of a constant amount of hydrazine as the amount of clay was increased. The increase was apparent (in spite of the water inhibition at high conversions) over a 2 order of magnitude variation of the clay weight. The weight dependence was taken to indicate that the role of the clay is very important, that the number of reactive centers is very small, or that they may be deactivated over the course of the reaction. In contrast to the strong dependence on overall amount of clay, the variation of amounts of putative oxidizing centers, such as structural Fe(III), admixed TiO2 or Fe2O3, or O- centers, did not result in alteration of the rate commensurate with the degree of variation of the entity in question. Surface iron does play some role, however, as samples that were pretreated with a reducing agent were less active as catalysts than the parent material. These results were taken to indicate either that the various centers interact to such a degree that they cannot be considered independently or that the reaction might proceed by way of surface complexation, rather than single electron transfers.

NASA Discipline Number 52-20↗

Space proton transport in one dimension

An approximate evaluation procedure is derived for a second-order theory of coupled nucleon transport in one dimension. An analytical solution with a simplified interaction model is used to determine quadrature parameters to minimize truncation error. Effects of the improved method on transport solutions with the BRYNTRN data base are evaluated. Comparisons with Monte Carlo benchmarks are given. Using different shield materials, the computational procedure is used to study the physics of space protons. A transition effect occurs in tissue near the shield interface and is most important in shields of high atomic number.

NASA Discipline Number 45-10↗

Characteristics of the Shuttle Orbiter Leeside Flow During A Reentry Condition

A study of the leeside flow characteristics of the Shuttle Orbiter is presented for a reentry flight condition. The flow is computed using a point-implicit, finite-volume scheme known as the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA). LAURA is a second-order accurate, laminar Navier-Stokes solver, incorporating finite-rate chemistry with a radiative equilibrium wall temperature distribution and finite-rate wall catalysis. The resulting computational solution is analyzed in terms of salient flow features and the surface quantities are compared with flight data.

Kleb, William L.↗

A dynamical model for reflex activated head movements in the horizontal plane

We present a controls systems model of horizontal-plane head movements during perturbations of the trunk, which for the first time interfaces a model of the human head with neural feedback controllers representing the vestibulocollic (VCR) and the cervicocollic (CCR) reflexes. This model is homeomorphic such that model structure and parameters are drawn directly from anthropomorphic, biomechanical and physiological studies. Using control theory we analyzed the system model in the time and frequency domains, simulating neck movement responses to input perturbations of the trunk. Without reflex control, the head and neck system produced a second-order underdamped response with a 5.2 dB resonant peak at 2.1 Hz. Adding the CCR component to the system dampened the response by approximately 7%. Adding the VCR component dampened head oscillations by 75%. The VCR also improved low-frequency compensation by increasing the gain and phase lag, creating a phase minimum at 0.1 Hz and a phase peak at 1.1 Hz. Combining all three components (mechanics, VCR and CCR) linearly in the head and neck system reduced the amplitude of the resonant peak to 1.1 dB and increased the resonant frequency to 2.9 Hz. The closed loop results closely fit human data, and explain quantitatively the characteristic phase peak often observed.

NASA Discipline Neuroscience↗

The Parametrized Post-Newtonian Gravitational Redshift

A derivation of the gravitational redshift effect to order c ^(-4) is presented. The calculation isperformed within the framework of the parametrized post-Newtonian formalism for analyzing metrictheories of gravity, which includes corrections to second- order in the Newtonian potential,gravitomagnetic contributions, and preferred-frame terms. We briefly discuss how to generalize ourresults to include possible violations of local Lorentz invariance or local position invariance which canarise in nonmetric theories. Our results are useful for analyzing possible new redshift experimentswhich may be sensitive to second-order effects, such as a close solar flyby mission.

Krisher, T. P.↗

The Effects of Noise on Polarimetric SAR Data

Polarimetric SAR data can provide a great deal of information about the scattering behavior of the surface under observation. Polarimetric SAR systems often measure the scattering matrices of the areas under observation in linear polarizations (H and V). From the scattering matrix commonly used forms such as the covariance matrix and the Stokes matrix can be easily derived. Other measures derived from polarimetric SAR data include correlation coefficients between scattering matrix terms and the mode and variance of phase differences between scattering matrix terms. The effects of additive system noise on these measurements is not often considered in the literature on this subject. In this paper, the effects of additive system noise on measurements derived from polarimetric SAR data will be examined. It will be shown how first-order noise effects can be removed and how second-order noise effects can be reduced for some measurements...

Freeman, A.↗

A Parallel Incompressible Navier-Stokes Solver with Multigrid Iterations

We developed a parallel, numerically accurate and stable, and computationally efficient finate-difference incompressible Navier-Stokes (N-S) fluid flow solver. The solver runs on both sequential and massively parallel computers. The numerical method used here is a second-order projection method on a staggered grid. The code is highly modular and it can be used either as a stand-alone flow solver and or a template code which can be adapted or expanded to a specific application. Numerical results and parallel performances of our code on Intel Delta and Paragon are reported.

Navier-Stokes solver↗

Resonance Cases and Small Divisors in a Third Integral of Motion

In this paper a general discussion of the resonance cases in an axially-symmetric potential field is presented, when the unperturbed frequencies in the radial and z direction have a rational ratio. The general form of the third integral is not valid in these cases because of the appearance of divisors of the form (m(exp 2)P-n(exp 2)Q), which become zero in the resonance cases. However, a new isolating integral of the unperturbed case is available, and this can be used to construct a third integral in the form of a power series and eliminate all secular terms. Three cases are distinguished, (alpha) m+n>4, (beta) m+n=4, and (gamma) m+n<4. In the first case the orbits are rather similar to those of the general irrational case. In the third case th e orbits show a quite peculiar character, which, however, can be explained rather accurately by a first-order theory of the third integral. Numerical integrations were made for the cases P=16Q, 4P=9Q, and P=4Q. The third integral, given in first- or second-order approximation, is rather well conserved. Case beta and the cases of small divisors, when m(exp 2)P-n(exp 2)Q is near zero but not equal to zero, are discussed in Paper II.

Contopoulos, George↗

An analysis of digital phase-locked loops

This report focuses on second-order digital phase-locked loops (DPLLs) with uniformly sampled input, an amplitude-insenstive phase extractor and a conventional loop filter.

Thomas, J. B.↗

Higher-Order Ionospheric Effects on the GPS Observables and Means of Modeling Them

Based on realistic modeling of the electron density of the ionosphere and using a dipole moment approximation for the earth magnetic field, we are able to estimate the effect of the ionosphere on the GPS signal for a ground user. Our study implies reduction of the second-order ionospheric error to less than 2 mm on the average and, therefore, potential improvement of orbit determination and baseline solutions.

ionosphere↗