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.
High energy electron beam studies dealing with nonlinear analysis of beam-plasma interactions, cyclotron harmonic instabilities, and frequency multiplication
Computation of a single geoidal height from gravity acceleration data formally requires that the latter be known everywhere on the earth. A computational procedure based on linear inverse theory for estimating geoidal heights from incomplete sets of data is presented. The same scheme can be used to estimate gravity accelerations from altimetry-derived geoids. The systematic error owing to lack of data and the choice of a particular inverse operator is described by using resolution functions and their spherical harmonic expansions. An rms value of this error is also estimated by assuming a spectrum for the unknown geoid. The influence of the size of the data region, the spacing between data, the filtering applied to the data, and the model weighting function chosen are all quantified in a spherical geometry. The examples presented show that when low degree spherical harmonic coefficients are available - from satellite orbit analysis - a band-passed version of the geoid can be constructed from local gravity data, even with a relatively restricted data set.
Geometric nonlinearities for large amplitude free and forced vibrations of circular plates are investigated. Inplane displacement and inertia are included in the formulation. Finite element method is used. Harmonic force matrix for nonlinear forced vibration analysis is introduced and derived. Various out-of-plane and inplane boundary conditions are considered. The relations of amplitude - frequency ratio for different boundary conditions and various loads conditions are presented.
The paper presents computational assessment of advanced numerical schemes for nonlinear acoustic problems related to combustion instabilities in liquid rocket engines. Several time-accurate, shock capturing schemes have been evaluated on a benchmark, closed-end resonant pipe flow problem. It involves the numerical solution of inviscid, compressible gas dynamics equations to predict acoustic wave propagation, wave steepening, formation of shocks, acoustic energy dissipation and wave-wall reflection for several hundred wave cycles. It was demonstrated that high accuracy TVD type schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time. The selected scheme was then applied to analyze the acoustic responses of resonant pipe-resonator, radial acoustic modes and hub-baffle configurations. Interesting observations of wave shape and damping characteristics have been drawn from presented computational studies.
The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.
Previous estimates on the degradations due to a subcarrier loop assume a square-wave subcarrier. This article provides a closed-form expression for the degradations due to the subcarrier loop when a finite number of harmonics are used to demodulate the subcarrier, as in the case of the buffered telemetry demodulator. We compared the degradations using a square wave and using finite harmonics in the subcarrier demodulation and found that, for a low loop signal-to-noise ratio, using finite harmonics leads to a lower degradation. The analysis is under the assumption that the phase noise in the subcarrier (SC) loop has a Tikhonov distribution. This assumption is valid for first-order loops.
A description of the joint test program objectives and scope is given together with the performance capability of the NASA Langley B-737 instrumented aircraft. The B-737 test run matrix conducted during the first 8 months of this 5-year program is discussed with a description of the different runway conditions evaluated. Some preliminary test results are discussed concerning the Electronic Recording Decelerometer (ERD) readings and a comparison of B-737 aircraft braking performance for different winter runway conditions. Detailed aircraft parameter time history records, analysis of ground vehicle friction measurements and harmonization with aircraft braking performance, assessment of induced aircraft contaminant drag, and evaluation of the effects of other factors on aircraft/ground vehicle friction performance will be documented in a NASA Technical Report which is being prepared for publication next year.
Surface pressure measurements on Mars have revealed a wide variety of atmospheric phenomena. The Mars Science Laboratory Rover Environmental Monitoring Station pressure sensor dataset is now the longest duration record of surface pressure on Mars. We use the first 2580 martian sols, nearly 4 Mars years,of measurements to identify atmospheric pressure waves with periods of tens of minutes to hours using wavelet analysis on residual pressure after the tidal harmonics are removed. We find these waves have a clear diurnal cycle with strongest activity in the early morning and late evening and a seasonal cycle with the strongest waves in the second half of the martian year (Ls= 180-360°). The strongest such waves of the entire mission occurred during the Mars Year 34 global dust storm. Comparable atmospheric waves are identified using atmospheric modeling with the MarsWRF general circulation model in a “nested” high spatial resolution mode. With the support of the modeling, we find these waves best fit the expected properties of inertia-gravity waves with horizontal wavelengths of O(100s) of km. PLAIN LANGUAGE SUMMARY Measuring air pressure from the surface of Mars has revealed a wide variety of atmospheric phenomena. The Curiosity rover's record or surface air pressure is now the longest yet made on Mars.We use the first ~8 years of Curiosity's pressure observations to look for atmospheric waves with periods of tens of minutes to hours. We find these waves have a clear pattern in their daily behavior with the strongest activity in the early morning and late evening and a seasonal cycle with the strongest waves in the second half of the martian year (Northern hemisphere fall and winter). The strongest such waves occurred in 2018 during a global dust storm. We find comparable waves in atmospheric modeling. With the support of modeling, we find these waves best fit the expected properties of buoyancy waves forced by airflow over topography with horizontal wavelengths of 100-1000 km.
A terrestrial upstream wave event which demonstrates multiple, ion cyclotron harmonic resonances between the interplanetary wave population and an observed proton beam is analyzed. The techniques and parameters employed in the data analysis are discussed, including the use of differential and band-pass filters. An upstream wave event demonstrating multiple harmonic waves is examined, and the instability analysis relevant to the ion beam observations thought to be responsible for that event is discussed. It is shown that an observed bi-Maxwellian ion beam is capable of generating right and left-hand polarized waves through ion cyclotron harmonic resonance.
Sikorsky Aircraft has performed analytical studies, design analyses, and risk reduction tests have been performed for Higher Harmonic Control (HHC) on the S-76. The S-76 is an 8 to 10,000 lb helicopter which cruises at 145 kts. Flight test hardware has been assembled, main servo frequency response tested and upgraded, aircraft control system shake tested and verified, open loop controllers designed and fabricated, closed loop controllers defined and evaluated, and rotors turning ground and flight tests planned for the near future. Open loop analysis shows that about 2 deg of higher harmonic feathering at the blade 75% radius will be required to eliminate 4P vibration in the cockpit.
The linear theory of electromagnetic ion beam instabilities for arbitrary angles of propagation is studied, with an emphasis on the conditions necessary to generate unstable modes at low harmonics of the ion cyclotron resonance condition. The present results extend the analysis of Smith et al. (1985). That paper considered only the plasma parameters at a time during which harmonic wave modes were observed in the earth's foreshock. The parameters of that paper are used as the basis of parametric variations here to establish the range of beam properties which may give rise to observable harmonic spectra. It is shown that the growth rates of both left-hand and right-hand cyclotron harmonic instabilities are enhanced by an increase in the beam temperature anisotropy and/or the beam speed. Decreases in the beam density and/or the core-ion beta reduce the overall growth of the cyclotron harmonic instabilities but favor the growth of these modes over the growth of the nonresonant instability and thereby enhance the observability of the harmonics.
Wavelet analysis is a relatively new technique that is an important addition to standard signal analysis methods. Unlike Fourier analysis that yields an average amplitude and phase for each harmonic in a dataset, the wavelet transform produces an instantaneous estimate or local value for the amplitude and phase of each harmonic. This allows detailed study of nonstationary spatial or time-dependent signal characteristics. The wavelet transform is discussed, examples are given, and some methods for preprocessing data for wavelet analysis are compared. By studying the dispersion of Yanai waves in a reduced gravity equatorial model, the usefulness of the transform is demonstrated. The group velocity is measured directly over a finite range of wavenumbers by examining the time evolution of the transform. The results agree well with linear theory at higher wavenumber but the measured group velocity is reduced at lower wavenumbers, possibly due to interaction with the basin boundaries.
The paper presents data on the accumulated V-band orbital phase curves of polarization in Cygnus X-1 which disagree with the analysis of Simmons, Aspin, and Brown (1980). No curves which conform to the canonical model were observed, and a noncanonical structure is indicated by F tests on third, fourth, and fifth harmonics, which have a probability of 0.85 as a group. The fifth harmonic has a probability of 0.99; in contrast, the fundamental harmonic has a probability of only 0.6, which is below the detection threshold. Therefore, major conclusions of Simmons et al based on a nonsymmetrical canonical model assuming that fundamental harmonics are present, but harmonics with n greater than 2 are absent, are unacceptable. Analysis of the author's data by Simmons et al is criticized on the grounds that he did not take into account the cumulative aspects of the data and did not recognize the need to average out the random noise generated by Cygnus X-1 over many orbital cycles.
The Space Shuttle Vehicle is assembled in the Vertical Assembly Building (VAB) at Kennedy Space Flight Center in Florida. The Vehicle is stacked on a Mobile Launch Platform (MLP) that weighs eight million pounds. A Crawler Transporter (CT) then carries the MLP and the stacked vehicle (12 million pounds total weight) to the launch complex located 5 miles away. This operation is performed at 0.9 mph resulting in a 4.5-hour transport. A recent test was performed to monitor the dynamic environment that was produced during rollout. It was found that the rollout is a harmonic-rich dynamic environment that was previously not understood. This paper will describe work that has been performed to estimate the forcing function that is produced in the transportation process. The rollout analysis team has determined that there are two families of harmonics of the drive train, which excite the system as a function of CT speed. There are also excitation sources, which are random or narrow-band in frequency and are not a function of CT speed. This presentation will discuss the application of the Sum of Weighted Accelerations Technique (SWAT) to further refine this understanding by estimating the forces and moments at the center-of-mass.
Induction motors having more than three phases were investigated to determine their suitability for electric vehicle applications. The objective was to have a motor with a current rating lower than that of a three-phase motor. The name chosen for these is high phase order (HPO) motors. Motors having six phases and nine phases were given the most attention. It was found that HPO motors are quite suitable for electric vehicles, and for many other applications as well. They have characteristics which are as good as or better than three-phase motors for practically all applications where polyphase induction motors are appropriate. Some of the analysis methods are presented, and several of the equivalent circuits which facilitate the determination of harmonic currents and losses, or currents with unbalanced sources, are included. The sometimes large stator currents due to harmonics in the source voltages are pointed out. Filters which can limit these currents were developed. An analysis and description of these filters is included. Experimental results which confirm and illustrate much of the theory are also included. These include locked rotor test results and full-load performance with an open phase. Also shown are oscillograms which display the reduction in harmonic currents when a filter is used with the experimental motor supplied by a non-sinusoidal source.
The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the nonlocal, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, such as those in the vicinity of roughness elements, steps, gaps, or corners. To compute the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 10^(7) to 10^(8) even in relatively simple flows. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear systems. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE. Additionally, inflow-resolvent analysis based on the HLNSE for discrete-roughness-induced disturbances in the nose tip of a blunt cone at Mach 6 demonstrates the importance of including the disturbance amplification along the near vicinity of the roughness element and separation region.
An efficient state-of-the-art method was developed to determine the unteady vibratory airloads produced by the interaction of the main rotor wake with a helicopter empennage. This method was incorporated into a computer program, Rotor Induced Empennage Vibration Analysis (RIEVA). The program requires the main rotor wake position and the strength of the vortices locates near the empennage surfaces. A nonlinear lifting surface analysis is utilized to predict the aerodynamic loads on the empennage surfaces in the presence of these concentrated vortices. The analysis was formulated to include all pertinent effects such as suction of the interacting vortices and the shed vorticity behind the empennage surfaces. The analysis employs a time domain solution. The output of the program consists of chordwise and spanwise airload distributions on the empennage surfaces. The airload distributions are harmonically analyzed and formulated for input into the Coupled Rotor/Airframe Vibration Analysis.
The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the non-local, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, as roughness elements, steps, gaps, or corners. To solve the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 107 to 108. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear system. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE.