Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “frequency estimation”

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,783 records · Page 99

Characterizing Oscillatory Bursts in Single-Trial EEG Data

Oscillatory bursts in numerous bands ranging from low (theta) to high frequencies (e.g., gamma) undoubtedly play an important role in cortical dynamics. Largely because of the inadequacy of existing analytic techniques. however, oscillatory bursts and their role in cortical processing remains poorly understood. To study oscillatory bursts effectively one must be able to isolate them and characterize them in the single trial. We describe a series of straightforward analysis techniques that produce useful indices of burst characteristics. First, stimulus-evoked responses are estimated using Differentially Variable Component Analysis (dVCA), and are subtracted from the single-trial. The single-trial characteristics of the evoked responses are stored to identify possible correlations with burst activity. Time-frequency (T-F), or wavelet, analyses are then applied to the single trial residuals. While T-F plots have been used in recent studies to identify and isolate bursts, we go further by fitting each burst in the T-F plot with a two-dimensional Gaussian. This provides a set of burst characteristics, such as, center time. burst duration, center frequency. frequency dispersion. and amplitude, all of which contribute to the accurate characterization of the individual burst. The burst phase can also be estimated. Burst characteristics can be quantified with several standard techniques (e.g.. histogramming and clustering), as well as Bayesian techniques (e.g., blocking) to allow a more parametric description analysis of the characteristics of oscillatory bursts, and the relationships of specific parameters to cortical excitability and stimulus integration.

Knuth, K. H.↗

Experimental Feedback Control of Flow Induced Cavity Tones

Discrete-time, linear quadratic methods were used to design feedback controllers for reducing tones generated by flow over a cavity. The dynamics of a synthetic jet actuator mounted at the leading edge of the cavity as observed by two microphones in the cavity were modeled over a broad frequency range using state space models computed from experimental data. Variations in closed loop performance as a function of model order, control order, control bandwidth, and state estimator design were studied using a cavity in the Probe Calibration Tunnel at NASA Langley. The controller successfully reduced the levels of multiple cavity tones at the tested flow speeds of Mach 0.275, 0.35, and 0.45. In some cases, the closed loop results were limited by excitation of sidebands of the cavity tones, or the creation of new tones at frequencies away from the cavity tones. Nonetheless, the results validate the combination of optimal control and experimentally-generated state space models, and suggest this approach may be useful for other flow control problems. The models were not able to account for non-linear dynamics, such as interactions between tones at different frequencies.

Cabell, Randolph H.↗

On-Orbit Lunar Modulation Transfer Function Measurements for the Moderate Resolution Imaging Spectroradiometer

Spatial quality of an imaging sensor can be estimated by evaluating its modulation transfer function (MTF) from many different sources such as a sharp edge, a pulse target, or bar patterns with different spatial frequencies. These well-defined targets are frequently used for prelaunch laboratory tests, providing very reliable and accurate MTF measurements. A laboratory-quality edge input source was included in the spatial-mode operation of the Spectroradiometric Calibration Assembly (SRCA), which is one of the onboard calibrators of the Moderate Resolution Imaging Spectroradiometer (MODIS). Since not all imaging satellites have such an instrument, SRCA MTF estimations can be used as a reference for an on-orbit lunar MTF algorithm and results. In this paper, the prelaunch spatial quality characterization process from the Integrated Alignment Collimator and SRCA is briefly discussed. Based on prelaunch MTF calibration using the SRCA, a lunar MTF algorithm is developed and applied to the lifetime on-orbit Terra and Aqua MODIS lunar collections. In each lunar collection, multiple scan-directionMoon-to-background transition profiles are aligned by the subpixel edge locations from a parametric Fermi function fit. Corresponding accumulated edge profiles are filtered and interpolated to obtain the edge spread function (ESF). The MTF is calculated by applying a Fourier transformation on the line spread function through a simple differentiation of the ESF. The lifetime lunar MTF results are analyzed and filtered by a relationship with the Sun-Earth-MODIS angle. Finally, the filtered lunarMTF values are compared to the SRCA MTF results. This comparison provides the level of accuracy for on-orbit MTF estimations validated through prelaunch SRCA measurements. The lunar MTF values had larger uncertainty than the SRCA MTF results; however, the ratio mean of lunarMTF fit and SRCA MTF values is within 2% in the 250- and 500-m bands. Based on the MTF measurement uncertainty range, the suggested lunar MTF algorithm can be applied to any on-orbit imaging sensor with lunar calibration capability.

Choi, Taeyong↗

Estimation of the longitudinal and lateral-directional aerodynamic parameters from flight data for the NASA F/A-18 HARV

This progress report presents the results of an investigation focused on parameter identification for the NASA F/A-18 HARV. This aircraft was used in the high alpha research program at the NASA Dryden Flight Research Center. In this study the longitudinal and lateral-directional stability derivatives are estimated from flight data using the Maximum Likelihood method coupled with a Newton-Raphson minimization technique. The objective is to estimate an aerodynamic model describing the aircraft dynamics over a range of angle of attack from 5 deg to 60 deg. The mathematical model is built using the traditional static and dynamic derivative buildup. Flight data used in this analysis were from a variety of maneuvers. The longitudinal maneuvers included large amplitude multiple doublets, optimal inputs, frequency sweeps, and pilot pitch stick inputs. The lateral-directional maneuvers consisted of large amplitude multiple doublets, optimal inputs and pilot stick and rudder inputs. The parameter estimation code pEst, developed at NASA Dryden, was used in this investigation. Results of the estimation process from alpha = 5 deg to alpha = 60 deg are presented and discussed.

Napolitano, Marcello R.↗

Flight Investigation of Prescribed Simultaneous Independent Surface Excitations for Real-Time Parameter Identification

Near real-time stability and control derivative extraction is required to support flight demonstration of Intelligent Flight Control System (IFCS) concepts being developed by NASA, academia, and industry. Traditionally, flight maneuvers would be designed and flown to obtain stability and control derivative estimates using a postflight analysis technique. The goal of the IFCS concept is to be able to modify the control laws in real time for an aircraft that has been damaged in flight. In some IFCS implementations, real-time parameter identification (PID) of the stability and control derivatives of the damaged aircraft is necessary for successfully reconfiguring the control system. This report investigates the usefulness of Prescribed Simultaneous Independent Surface Excitations (PreSISE) to provide data for rapidly obtaining estimates of the stability and control derivatives. Flight test data were analyzed using both equation-error and output-error PID techniques. The equation-error PID technique is known as Fourier Transform Regression (FTR) and is a frequency-domain real-time implementation. Selected results were compared with a time-domain output-error technique. The real-time equation-error technique combined with the PreSISE maneuvers provided excellent derivative estimation in the longitudinal axis. However, the PreSISE maneuvers as presently defined were not adequate for accurate estimation of the lateral-directional derivatives.

Moes, Timothy R.↗

Ion Number Density Quantification Utilizing Pulsing Frequency in Negative Differential Resistance (NDR) Regime of Microplasma Operation

A mathematical model is presented for quantifying the ion number density over a range of discharge current in microplasma discharges. The ion number density is determined from measured discharge voltage, current, pulsing frequency, and trace impurity that is injected into the system deliberately. Experiments are conducted in pure helium with nitrogen being introduced as the trace species. The ion number density is estimated over a range of discharge current and is compared with multidimensional simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Estimates of Long Term Surface Soil Moisture in the Midwestern U.S. Derived from Satellite Microwave Observations

Soil moisture is a key component of the water and energy balances of the Earth's surface, and has been identified as a parameter of significant potential for improving the accuracy of large-scale land surface-atmosphere interaction models. However, soil moisture is often somewhat difficult to measure accurately in both space and time, especially at large spatial scales. Soil moisture is highly variable, and while point measurements are typically quite accurate, subsequent areal averaging of these measurements often leads to large errors. Since remotely sensed land surface observations are already a spatially averaged or areally integrated value, they are a logical input parameter to regional or larger scale land process models. A database of long-term soil moisture was compared to satellite microwave observations over test sites in the Midwestern United States. Ground measurements of average volumetric surface soil moisture in the top ten cm were made bimonthly at 19 locations throughout the state of Illinois. Nighttime microwave brightness temperatures were observed at a frequency of 6.6 GHz, by the Scanning Multichannel Microwave Radiometer (SMMR), onboard the Nimbus 7 satellite. The life of the SMMR instrument spanned from Nov. 1978 to Aug. 1987. At 6.6 GHz, the instrument provided a spatial resolution of approximately 150 km, and a temporal frequency over the test area of about 3 nighttime orbits per week. Vegetation radiative transfer characteristics, such as the canopy transmissivity, were estimated from vegetation indices such as the Normalized Difference Vegetation Index (NDVI) and the 37 GHz Microwave Polarization Difference Index (MPDI). Because the time of satellite coverage does not always coincide with the ground measurements of soil moisture, the existing ground data were used to calibrate a water balance for the top IO cm surface layer in order to interpolate daily surface moisture values. Such a climate-based approach is often more appropriate for estimating large-area average soil moisture because meteorological data are generally more spatially representative than isolated point measurements of soil moisture, Passive microwave remote sensing presents the greatest potential for providing regular spatially representative estimates of surface soil moisture at global scales. Real time estimates should improve weather and climate modelling efforts, while the development of historical data sets will provide necessary information for simulation and validation of long-term climate and global change studies.

Owe, M.↗

Solar Dynamics Observatory On-Orbit Jitter Testing, Analysis, and Mitigation Plans

The Solar Dynamics Observatory (SDO) was designed to understand the Sun and the Sun s influence on Earth. SDO was launched on February 11, 2010 carrying three scientific instruments: the Atmospheric Imaging Assembly (AIA), the Helioseismic and Magnetic Imager (HMI), and the Extreme Ultraviolet Variability Experiment (EVE). Both AIA and HMI are sensitive to high frequency pointing perturbations and have sub-arcsecond level line-of-sight (LOS) jitter requirements. Extensive modeling and analysis efforts were directed in estimating the amount of jitter disturbing the science instruments. To verify the disturbance models and to validate the jitter performance prior to launch, many jitter-critical components and subassemblies were tested either by the mechanism vendors or at the NASA Goddard Space Flight Center (GSFC). Although detailed analysis and assembly level tests were performed to obtain good jitter predictions, there were still several sources of uncertainties in the system. The structural finite element model did not have all the modes correlated to test data at high frequencies (greater than 50 Hz). The performance of the instrument stabilization system was not known exactly but was expected to be close to the analytical model. A true disturbance-to-LOS observatory level test was not available due to the tight schedule of the flight spacecraft, the cost in time and manpower, difficulties in creating gravity negation systems, and risks of damaging flight hardware. To protect the observatory jitter performance against model uncertainties, the SDO jitter team devised several on-orbit jitter reduction plans in addition to reserve margins on analysis results. Since some of these plans severely restricted the capabilities of several spacecraft components (e.g. wheels and High Gain Antennas), the SDO team performed on-orbit jitter tests to determine which jitter reduction plans, if any, were necessary to satisfy science LOS jitter requirements. The SDO on-orbit jitter tests were designed to satisfy the following four objectives: 1. Determine the acceptable reaction wheel operational speed range during Science Mode. 2. Determine HGA algorithm jitter parameters (number of stagger steps and enable/disable no-steprequests). 3. Determine acceptable EVE instrument filter wheels spin rates. 4. Determine if AIA instrument filter wheels excite the first AIA telescope structural mode. This paper provides detailed information on the SDO wheel jitter test plan, shows on-orbit jitter measurements and how ground predictions compare to those measurements, and describes the final jitter mitigation plan executed on SDO.

Liu, Kuo-Chia (Alice)↗

Excitations for Rapidly Estimating Flight-Control Parameters

A flight test on an F-15 airplane was performed to evaluate the utility of prescribed simultaneous independent surface excitations (PreSISE) for real-time estimation of flight-control parameters, including stability and control derivatives. The ability to extract these derivatives in nearly real time is needed to support flight demonstration of intelligent flight-control system (IFCS) concepts under development at NASA, in academia, and in industry. Traditionally, flight maneuvers have been designed and executed to obtain estimates of stability and control derivatives by use of a post-flight analysis technique. For an IFCS, it is required to be able to modify control laws in real time for an aircraft that has been damaged in flight (because of combat, weather, or a system failure). The flight test included PreSISE maneuvers, during which all desired control surfaces are excited simultaneously, but at different frequencies, resulting in aircraft motions about all coordinate axes. The objectives of the test were to obtain data for post-flight analysis and to perform the analysis to determine: 1) The accuracy of derivatives estimated by use of PreSISE, 2) The required durations of PreSISE inputs, and 3) The minimum required magnitudes of PreSISE inputs. The PreSISE inputs in the flight test consisted of stacked sine-wave excitations at various frequencies, including symmetric and differential excitations of canard and stabilator control surfaces and excitations of aileron and rudder control surfaces of a highly modified F-15 airplane. Small, medium, and large excitations were tested in 15-second maneuvers at subsonic, transonic, and supersonic speeds. Typical excitations are shown in Figure 1. Flight-test data were analyzed by use of pEst, which is an industry-standard output-error technique developed by Dryden Flight Research Center. Data were also analyzed by use of Fourier-transform regression (FTR), which was developed for onboard, real-time estimation of the derivatives.

Moes, Tim↗

Using Dense Array Waveform Correlations to Build a Velocity Model with Stochastic Variability

In an ongoing effort to improve 3D seismic-wave propagation modeling for frequencies up to 10 Hz, we used cross correlations between vertical-component waveforms from an underground chemical explosion to estimate the statistical properties of small-scale velocity heterogeneities. The waveforms were recorded by a dense 2D seismic array deployed during the Source Physics Experiments for event number 5 (SPE-5) in a series of six underground chemical explosions, conducted at the Nevada National Security Site. The array consisted of 996 geophones with a 50–100 m grid spacing, deployed at the SPE site at the north end of the Yucca Flat basin. The SPE were conducted to investigate the generation and propagation of seismic and acoustic waves from underground explosions. Comparisons of decay rates of waveform cross correlations as function of interstation distance, computed for observed and synthetic seismograms from the SPE-5 chemical explosion, were used to constrain statistical properties of correlated stochastic velocity perturbations representing small-scale heterogeneities added to a geology-based velocity model of the Yucca Flat basin. Using comparisons between recorded and simulated waveform cross correlations, we were able to recover sets of statistical properties of small-scale velocity perturbations in the velocity model that produce the best-fit between the recorded and simulated ground motion. The stochastic velocity fluctuations in the velocity model that produced the smallest misfits have a horizontal correlation distance of between 400 and 800 m, a vertical correlation distance between 100 and 200 m, and a standard deviation of 10% from the nominal model velocity in the alluvium basin layers. They also have a horizontal correlation distance of 1000 m, a vertical correlation distance of 250 m, and a standard deviation of 6% in the underlying and consolidated sedimentary layers, up to a depth of 4 km. Comparisons between observed and simulated wavefields were used to assess the proposed small-scale heterogeneity enhancements to the Yucca Flat basin model. We found that adding a depth-resolved stochastic variability to the geology-based velocity model improves the overall performance of ground-motion simulations of an SPE-5 explosion in the modeled frequency range up to 10 Hz. The results may be applicable to other similar basins.

58 GEOSCIENCES↗

The role of photospheric magnetic fields in the variation of solar oscillation eigenfrequencies

The frequencies of solar oscillations have been found to be higher at the peak of the solar cycle, their variation being consistent with that expected from a change in the solar structure confined to the Sun's outer layers. In this work, a technique is introduced for estimating the effect of the highly inhomogeneous photospheric magnetic field on solar oscillation eigenfrequencies, and its results are compared with the solar data. It is found that the observed steep rise in frequency shift at low frequency and the leveling of or decline at high frequency are well reproduced by the model. The correct order of magnitude for the frequency changes is also reproduced, at least at low frequency. The model implies a particular functional form for the variation of the frequency changes with order and degree, and it is suggested that future advances in observation and theory will therefore allow this model to be confirmed or refuted.

Rosenthal, C. S.↗

Estimates of the wavenumber wavelet power spectrum of magnetic fluctuations during magnetic reconnection

Fluctuation analyses of experimental observations generally lack high temporal resolution and are in frequency-space f, contrary to theoretical efforts in wavenumber-space k. This is due to the inherent limits of the Fourier transform, though it is prominent due to the ease of diagnostic implementation. Advances in wavelet-based analysis have provided relief due to its temporal resolution, but in its common use, is still hard to compare to theoretical models. By using the two-point correlation technique in conjunction with large data sets, a wavelet power spectrum in wavenumber-space can be created. Dubbed the wavenumber wavelet power spectrum, this spectrum relates wavenumber to power in time. Further, this analysis technique more closely connects characterizations of experimentally observed fluctuations with other system parameters and theoretical predictions. In this article, we develop the wavenumber wavelet power spectrum using magnetic fluctuations caused by tearing instability driven magnetic reconnection in reproducible, high temperature laboratory plasmas. These dynamic magnetic fluctuations generated in reversed field pinch plasmas are broadband, ranging from the low frequency, 10's of kHz, up to the ion gyroradii frequencies, 100's of kHz. The dominant fluctuations have poloidal and toroidal mode numbers (m,n)=(1,6−10) and can grow to 2%–3% of the mean magnetic field. During these reconnection events, ions, and electrons are energized, magnetic fluctuation amplitudes increase, plasma flow is halted, and the toroidal magnetic flux increases, all on a semi-periodic basis. The newly developed spectrum provides better temporal resolution of spectrum characteristics to correlate with these particle energization phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Autonomous frequency domain identification: Theory and experiment

The analysis, design, and on-orbit tuning of robust controllers require more information about the plant than simply a nominal estimate of the plant transfer function. Information is also required concerning the uncertainty in the nominal estimate, or more generally, the identification of a model set within which the true plant is known to lie. The identification methodology that was developed and experimentally demonstrated makes use of a simple but useful characterization of the model uncertainty based on the output error. This is a characterization of the additive uncertainty in the plant model, which has found considerable use in many robust control analysis and synthesis techniques. The identification process is initiated by a stochastic input u which is applied to the plant p giving rise to the output. Spectral estimation (h = P sub uy/P sub uu) is used as an estimate of p and the model order is estimated using the produce moment matrix (PMM) method. A parametric model unit direction vector p is then determined by curve fitting the spectral estimate to a rational transfer function. The additive uncertainty delta sub m = p - unit direction vector p is then estimated by the cross spectral estimate delta = P sub ue/P sub uu where e = y - unit direction vectory y is the output error, and unit direction vector y = unit direction vector pu is the computed output of the parametric model subjected to the actual input u. The experimental results demonstrate the curve fitting algorithm produces the reduced-order plant model which minimizes the additive uncertainty. The nominal transfer function estimate unit direction vector p and the estimate delta of the additive uncertainty delta sub m are subsequently available to be used for optimization of robust controller performance and stability.

Yam, Yeung↗

Estimating Canopy Water Content of Chaparral Shrubs Using Optical Methods

California chaparral ecosystems are exceptionally fire adapted and typically are subject to wildfire at decadal to century frequencies. The hot dry Mediterranean climate summers and the chaparral communities of the Santa Monica Mountains make wildfire one of the most serious economic and life-threatening natural disasters faced by the region. Additionally, the steep fire-burned hillsides are subject to erosion, slumpage, and mud slides during the winter rains. The Santa Monica Mountain Zone (SMMZ) is a 104,000 ha eastwest trending range with 607 m of vertical relief and located in the center of the greater Los Angeles region. A series of fires in the fall of 1993 burned from Simi Valley to Santa Monica within a few hours. Developing techniques to monitor fire hazard and predict the spread of fire is of major concern to the region. One key factor in the susceptibility to fire is the water content of the vegetation canopy. The development of imaging spectrometry and remote sensing techniques may constitute a tool to provide this information.

Ustin, Susan L.↗

Estimating Canopy Water Content of Chaparral Shrubs Using Optical Methods

California chaparral ecosystems are exceptionally fire adapted and typically are subject to wildfire at decadal to century frequencies. The hot dry Mediterranean climate summers and the chaparral communities of the Santa Monica Mountains make wildfire one of the most serious economic and life-threatening natural disasters faced by the region.

ecosystems↗

Nanoflare Heating Frequency of an X-ray Bright Point Observed by MaGIXS

Nanoflares have been considered to be one of the most likely candidates for heating the solar corona to multi-million kelvin temperatures. Individual nanoflares are difficult to detect with today's instruments, but their presence may be established by comparing simulated nanoflare-heated plasma emissions to observed emissions. We present a simulation of emission from an X-ray Bright Point (XBP) detected by the MaGIXS, as well as simultaneous observations from SDO/AIA and Hinode/XRT. To simulate the XBP loops, we utilize the HYDRAD code. The length and magnetic field strength of these loops are determined using potential field extrapolation of SDO/HMI's observed photospheric magnetogram. Each loop is considered to be heated by random nanoflares, the amplitude and frequency of which are governed by the length of the loop and the strength of the magnetic field. The simulated outputs are used to estimate the intensity of spectrally pure maps of Fe-18, Fe-17, Ne-9, O-8, O-9, Ne-9, and so on, which is then compared to the intensity determined from MaGIXS observations. In addition, we derived the intensity maps obtained by AIA and XRT and compared them to the observed data. The composite distribution of the delay time of the nanoflares for which the simulated loops morphology and intensities match with observation shows a peak at 200s-500s, indicating that most of the nanoflares have a high/intermediate frequency.

coronal heating↗

Evolution of global snow drought characteristics from 1850 to 2100

Abstract Seasonal snow is an integral part of the global water supply and storage system. Snow droughts impact ecological, agricultural, and urban systems by altering the amount and timing of meltwater delivery. These droughts are characterized by a lack of on-the-ground snow (snow water equivalent, SWE) that can be caused by low total precipitation (dry drought) or low proportion of precipitation falling as snowfall (warm drought), often combined with an early melt. The standardized SWE index (SWEI) ranks the current status of SWE for a given location compared to a baseline condition and identifies the existence, but not the cause, of snow drought. In this work, we use estimates of SWE, temperature, and precipitation from nine coupled model intercomparison project phase 6 (CMIP6) models to quantify the frequency, severity, and type of snow droughts globally for historical and future scenarios. Compared to a historical baseline (1850–1900) total snow drought frequency more than doubles under socio-economic pathway (SSP)2-4.5 and SSP5-8.5; all of the increase in snow drought frequency comes from an increase in warm droughts. The probability distribution of future SWEI in major snowy basins around the world are likely to be centered on more negative values, which corresponds to more severe drought and, with only moderate changes in distribution spread, more frequent drought. CMIP6 simulations pinpoint snow drought as an emerging global threat to water resources and highlight the need to explore higher resolution future models that better capture complex mountain topography, wildland fires, and snow-forest interactions.

CMIP6↗