Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “filter 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 721 records · Page 40

Generating Control Commands From Gestures Sensed by EMG

An effort is under way to develop noninvasive neuro-electric interfaces through which human operators could control systems as diverse as simple mechanical devices, computers, aircraft, and even spacecraft. The basic idea is to use electrodes on the surface of the skin to acquire electromyographic (EMG) signals associated with gestures, digitize and process the EMG signals to recognize the gestures, and generate digital commands to perform the actions signified by the gestures. In an experimental prototype of such an interface, the EMG signals associated with hand gestures are acquired by use of several pairs of electrodes mounted in sleeves on a subject s forearm (see figure). The EMG signals are sampled and digitized. The resulting time-series data are fed as input to pattern-recognition software that has been trained to distinguish gestures from a given gesture set. The software implements, among other things, hidden Markov models, which are used to recognize the gestures as they are being performed in real time. Thus far, two experiments have been performed on the prototype interface to demonstrate feasibility: an experiment in synthesizing the output of a joystick and an experiment in synthesizing the output of a computer or typewriter keyboard. In the joystick experiment, the EMG signals were processed into joystick commands for a realistic flight simulator for an airplane. The acting pilot reached out into the air, grabbed an imaginary joystick, and pretended to manipulate the joystick to achieve left and right banks and up and down pitches of the simulated airplane. In the keyboard experiment, the subject pretended to type on a numerical keypad, and the EMG signals were processed into keystrokes. The results of the experiments demonstrate the basic feasibility of this method while indicating the need for further research to reduce the incidence of errors (including confusion among gestures). Topics that must be addressed include the numbers and arrangements of electrodes needed to acquire sufficient data; refinements in the acquisition, filtering, and digitization of EMG signals; and methods of training the pattern- recognition software. The joystick and keyboard simulations were chosen for the initial experiments because they are familiar to many computer users. It is anticipated that, ultimately, interfaces would utilize EMG signals associated with movements more nearly natural than those associated with joysticks or keyboards. Future versions of the pattern-recognition software are planned to be capable of adapting to the preferences and day-today variations in EMG outputs of individual users; this capability for adaptation would also make it possible to select gestures that, to a given user, feel the most nearly natural for generating control signals for a given task (provided that there are enough properly positioned electrodes to acquire the EMG signals from the muscles involved in the gestures).

Wheeler, Kevin R.↗

Image sharpening for mixed spatial and spectral resolution satellite systems

Two methods of image sharpening (reconstruction) are compared. The first, a spatial filtering technique, extrapolates edge information from a high spatial resolution panchromatic band at 10 meters and adds it to the low spatial resolution narrow spectral bands. The second method, a color normalizing technique, is based on the ability to separate image hue and brightness components in spectral data. Using both techniques, multispectral images are sharpened from 30, 50, 70, and 90 meter resolutions. Error rates are calculated for the two methods and all sharpened resolutions. The results indicate that the color normalizing method is superior to the spatial filtering technique.

Hallada, W. A.↗

Plastic work constrained elastoplastic topology optimization

An elastoplastic topology optimization framework for limiting plastic work generation while maximizing stiffness is presented. The kinematics and constitutive model are based on finite strain linear isotropic hardening plasticity, and the balance laws are solved using a total Lagrangian finite element formulation. Aggregation of the specific plastic work combined with an adaptive normalization scheme efficiently constrains the maximum specific plastic work. The optimization problem is regularized using an augmented partial differential equation filter, and is solved by the method of moving asymptotes where path-dependent sensitivities are derived using the adjoint method. The numerical examples show a clear dependence on the optimized maximum stiffness structures for different levels of constrained specific plastic work. It is also shown that due to the history dependency of the plasticity, the load path significantly influences the structural performance and optimized topology.

42 ENGINEERING↗

Development of a hardened THz energy meter for use on the kilojoule-scale, short-pulse OMEGA EP laser

A highly adaptable and robust terahertz (THz) energy meter is designed and implemented to detect energetic THz pulses from high-intensity (>10 18 W/cm 2 ) laser–plasma interactions on the OMEGA EP. THz radiation from the laser driven target is detected by a shielded pyrometer. A second identical pyrometer is used for background subtraction. The detector can be configured to detect THz pulses in the 1 mm to 30 μm (0.3- to 10-THz) range and pulse energies from joules to microjoules via changes in filtration, aperture size, and position. Additional polarization selective filtration can also be used to determine the THz pulse polarization. The design incorporates significant radiation and electromagnetic pulse shielding to survive and operate within the OMEGA EP radiation environment. We describe the design, operational principle, calibration, and testing of the THz energy meter. The pyrometers were calibrated using a benchtop laser and show linear sensitivity to up to 1000 nJ of absorbed energy. Furthermore, the initial results from four OMEGA EP THz experiments detected up to ~15 μJ at the detector, which can correspond to hundreds of mJ depending on THz emission and reflection models.

47 OTHER INSTRUMENTATION↗

Characterization of core neutrals using passive measurements of the D α spectrum near the X -point in the DIII-D tokamak

New spectroscopic measurements of deuterium Balmer-α emission are used to infer the spatial distribution of neutral particles near the X-point in a diverted high confinement mode plasma. The charge exchange neutral spectroscopy (CENS) diagnostic—recently installed on the DIII-D tokamak—uses 15 lines of sight extending from the edge of the confined region (ρ ≈ 0.7) to the X-point in lower single-null plasmas. Each CENS line of sight is spectrally resolved around the deuterium Balmer-α wavelength (6561 Å) to measure the Doppler shift, line broadening, and intensity of Dα emission. Thermal line broadening is used to identify emission from neutrals undergoing charge-exchange (CX) with high temperature ions in the confined plasma. This spectral information allows a more accurate determination of the neutral density deeper inside the confined plasma compared to traditional filter-based diagnostics. We present two methods of analyzing CENS measurements: (1) fitting the spectrum from each view-chord separately, and (2) a tomographic method for inverting the neutral density over a 2D region of space using all CENS views collectively. The neutral density profile is found to decay exponentially in the radial direction with two scale-lengths. In the pedestal region the neutral density decays at a rate approximately equal to the local mean-free-path for CX collisions, $L_{n_\mathrm{D0}}\approx \lambda_\text{CX}$. Further inside the plasma the neutral density decays at a rate equal to the mean-free-path for ionization, $L_{n_\mathrm{D0}}\approx \lambda_\text{inz.}$. The separatrix $n_\mathrm{D0}$ value is found to be approximately $2\times10^{15}$ m−3. Based on results from the 2D inversion, the density of neutrals is found to increase along the separatrix approaching the X-point by a factor of $\approx 5\times$ over the region covered by the CENS diagnostic.

X-point neutrals↗

Precipitation mapping with an airborne synthetic aperture imaging radar

Use of an airborne synthetic-aperture radar (SAR) of 3.1-cm wavelength for the detection and mapping of precipitation echoes is described. Difficulties in mapping related to the incoherence of the echoes are discussed; methods of increasing resolution by Doppler filtering are suggested. Observations made with a vertical side-looking SAR system flying above a stratiform storm were found to resemble vertical cross sections obtained with conventional radar showing snow streamers, the bright band and rain. However, the actual section sampled by the SAR after filtering for zero Doppler is shown to depend upon both the relative motion of the particles and the relative winds below the aircraft. Knowledge of the relative winds, deduced from Doppler spectrum variations, and of the particle fall speeds permits the generation of horizontal and vertical cross sections displaced from the aircraft's position.

Atlas, D.↗

Analysis of time dependent phenomena observed with the LPSP OSO-8 instrument

Data obtained by the Laboratoire de Physique Stellaire et Planetaire's ultraviolet spectrometer onboard the OSO-8 spacecraft were analyzed in an effort to dynamically model the solar chromosphere as an aid in enhancing knowledge of the dynamical processes themselves and of spectral line formation in the dynamic chromosphere. Repeated spectral scans of strong, optically thick resonance lines formed in the solar chromosphere were examined for indications of oscillatory velocities and intensities among other indications of velocity which were studied, the blue peak is reasonably well defined, and the position of a parabolic filter fitted by the least squares method was used to define it. Observed chromospheric oscillation periods are discussed as well as the variations in altitude of the emitting region which result primarily from the motion up and down during the oscillation.

Leibacher, J. W.↗

Fast algorithm for calculating chemical kinetics in turbulent reacting flow

This paper addresses the need for a fast batch chemistry solver to perform the kinetics part of a split operator formulation of turbulent reacting flows, with special attention focused on the solution of the ordinary differential equations governing a homogeneous gas-phase chemical reaction. For this purpose, a two-part predictor-corrector algorithm which incorporates an exponentially fitted trapezoidal method was developed. The algorithm performs filtering of ill-posed initial conditions, automatic step-size selection, and automatic selection of Jacobi-Newton or Newton-Raphson iteration for convergence to achieve maximum computational efficiency while observing a prescribed error tolerance. The new algorithm, termed CREK1D (combustion reaction kinetics, one-dimensional), compared favorably with the code LSODE when tested on two representative problems drawn from combustion kinetics, and is faster than LSODE.

Radhakrishnan, K.↗

ENSO Bred Vectors in Coupled Ocean-Atmosphere General Circulation Models

The breeding method has been implemented in the NASA Seasonal-to-Interannual Prediction Project (NSIPP) Coupled General Circulation Model (CGCM) with the goal of improving operational seasonal to interannual climate predictions through ensemble forecasting and data assimilation. The coupled instability as cap'tured by the breeding method is the first attempt to isolate the evolving ENSO instability and its corresponding global atmospheric response in a fully coupled ocean-atmosphere GCM. Our results show that the growth rate of the coupled bred vectors (BV) peaks at about 3 months before a background ENSO event. The dominant growing BV modes are reminiscent of the background ENSO anomalies and show a strong tropical response with wind/SST/thermocline interrelated in a manner similar to the background ENSO mode. They exhibit larger amplitudes in the eastern tropical Pacific, reflecting the natural dynamical sensitivity associated with the presence of the shallow thermocline. Moreover, the extratropical perturbations associated with these coupled BV modes reveal the variations related to the atmospheric teleconnection patterns associated with background ENSO variability, e.g. over the North Pacific and North America. A similar experiment was carried out with the NCEP/CFS03 CGCM. Comparisons between bred vectors from the NSIPP CGCM and NCEP/CFS03 CGCM demonstrate the robustness of the results. Our results strongly suggest that the breeding method can serve as a natural filter to identify the slowly varying, coupled instabilities in a coupled GCM, which can be used to construct ensemble perturbations for ensemble forecasts and to estimate the coupled background error covariance for coupled data assimilation.

Yang, S. C.↗

Orion MPCV Touchdown Detection Threshold Development and Testing

A robust method of detecting Orion Multi-Purpose Crew Vehicle (MPCV) splashdown is necessary to ensure crew and hardware safety during descent and after touchdown. The proposed method uses a triple redundant system to inhibit Reaction Control System (RCS) thruster firings, detach parachute risers from the vehicle, and transition to the post-landing segment of the Flight Software (FSW). An in-depth trade study was completed to determine optimal characteristics of the touchdown detection method resulting in an algorithm monitoring filtered, lever-arm corrected, 200 Hz Inertial Measurement Unit (IMU) vehicle acceleration magnitude data against a tunable threshold using persistence counter logic. Following the design of the algorithm, high fidelity environment and vehicle simulations, coupled with the actual vehicle FSW, were used to tune the acceleration threshold and persistence counter value to result in adequate performance in detecting touchdown and sufficient safety margin against early detection while descending under parachutes. An analytical approach including Kriging and adaptive sampling allowed for a sufficient number of finite element analysis (FEA) impact simulations to be completed using minimal computation time. The combination of a persistence counter of 10 and an acceleration threshold of approximately 57.3 ft/s2 resulted in an impact performance factor of safety (FOS) of 1.0 and a safety FOS of approximately 2.6 for touchdown declaration. An RCS termination acceleration threshold of approximately 53.1 ft/s(exp)2 with a persistence counter of 10 resulted in an increased impact performance FOS of 1.2 at the expense of a lowered under-parachutes safety factor of 2.2. The resulting tuned algorithm was then tested on data from eight Capsule Parachute Assembly System (CPAS) flight tests, showing an experimental minimum safety FOS of 6.1. The formulated touchdown detection algorithm will be flown on the Orion MPCV FSW during the Exploration Flight Test 1 (EFT-1) mission in the second half of 2014.

Daum, Jared↗

Direct Simulation of Flow over A Wall-Mounted Hump

Simulating separated flows at high Reynolds numbers using Reynolds Averaged Navier-Stokes (RANS) modeled equations remains a challenge in aeronautics. The main hindrance to progress stems from the lack of extended de- tailed data pertinent to the root cause of failure of RANS models, as well as the little progress in RANS modeling innovations in the past several decades. The goal of the current effort is to generate data for separated flow at a Reynolds numbers where conventional models are challenged. We use Direct Numerical Simulations to model turbulent flow over the wall-mounted hump configuration to investigate the physics of flow separation and boundary layer recovery, as well as provide data relevant to the modeling community. A chord-based Reynolds number of Re(sub c) = 47,500 is considered with a turbulent inflow profile of Re(sub θ) = 1,400 (θ /c = 3%). We use FDL3DI, a code that solves the compressible Navier-Stokes equations using high-order compact-difference scheme and filter, with the standard re- cycling/rescaling method of generating turbulent boundary layers as inflow to the computational domain. Two differ- ent configurations of the upper-wall are analyzed for two sets of boundary conditions (slip and no-slip). The results are compared with the available higher Re(sub c) (= 936,000, Re(sub θ) = 7,200, θ /c = 0.77%) experiment for major flow features. The simulated lower Rec allows for DNS-like mesh resolutions, and adequately wide spans. The results from these simulations show earlier separation and delayed reattachment compared to Re(sub c) = 936,000, and significantly higher skin friction in the forebody of the hump. We also find that the upper-wall shape and boundary condition influence pressure distribution over the hump, whereas skin friction is only influenced by the boundary condition.

Wall-Mounted↗

Biomass Estimation from Simulated GEDI, ICESat-2 and NISAR Across Environmental Gradients in Sonoma County, California

Estimates of the magnitude and distribution of aboveground carbon in Earth’s forests remain uncertain, yet knowledge of forest carbon content at a global scale is critical for forest management in support of climate mitigation. In light of this knowledge gap, several upcoming spaceborne missions aim to map forest aboveground biomass, and many new biomass products are expected from these datasets. As these new missions host different technologies, each with relative strengths and weaknesses for biomass retrieval, as well as different spatial resolutions, consistently comparing or combining biomass estimates from these new datasets will be challenging. This paper presents a demonstration of an inter-comparison of biomass estimates from simulations of three NASA missions (GEDI, ICESat-2 and NISAR) over Sonoma county in California, USA. We use a high resolution, locally calibrated airborne lidar map as our reference dataset, and emphasize the importance of considering uncertainties in both reference maps and spaceborne estimates when conducting biomass product validation. GEDI and ICESat-2 were simulated from airborne lidar point clouds, while UAVSAR’s L-band backscatter was used as a proxy for NISAR. To estimate biomass for the lidar missions we used GEDI’s footprint-level biomass algorithms, and also adapted these for application to ICESat-2. For UAVSAR, we developed a locally trained biomass model, calibrated against the ALS reference map. Each mission simulation was evaluated in comparison to the local reference map at its native product resolution (25 m, 100m transect, and 1 ha) yielding RMSEs of 57%, 75%, and 89% for GEDI, NISAR, and ICESat-2 respectively. RMSE values increased for GEDI’s power beam during simulated daytime conditions (64%), coverage beam during nighttime conditions (72%), and coverage beam daytime conditions (87%). We also test the application of GEDI’s biomass modeling framework for estimation of biomass from ICESat-2, and fine that ICESat-2 yields reasonable biomass estimates, particularly in relatively short, open canopies. Results suggest that while all three missions will produce datasets useful for biomass mapping, tall, dense canopies such as those found in Sonoma County present the greatest challenges for all three missions, while steep slopes also prove challenging for single-date SAR based biomass retrieval. Our methods provide guidance for the inter-comparison and validation of spaceborne biomass estimates through the use of airborne lidar reference maps, and could be repeated with on-orbit estimates in any area with high quality field plot and ALS data. These methods allow for regional interpretations and filtering of multi-mission biomass estimates toward improved wall-to-wall biomass maps through data fusion.

GEDI↗

Collective Analysis of Alpha Particle Losses Due to Self-Absorption by Mass Loading on Radioactive Particulate Glass Fiber Filters

In this study, we derived a relationship between filter mass loading and the percent loss during analysis using the mass loading data collected from six previous studies of self-absorption. Components of mass loading include particulate dust, radioactive particulates, and filter material. In a research report published in 1984, Higby calculated a minimum burial depth for an alpha particle to be lost due to absorption (100% loss) of about 3.7 mg/cm 2 based on calculations for the range of 239-Pu alpha particles in glass fiber filters. From there, Higby concluded that a correction factor of 0.85 assumes approximately 15% losses in the count rate of both alpha and beta particles. In 2000, Luetzelschwab et al. recommended assuming a 40% loss at a loading of 3.3 mg/cm 2 and a 28% loss for a loading of 2.3 mg/cm 2 which included the frontal face mass of the filter. More recently, the 100% losses due to absorption were reported to be in the 10 mg/cm 2 range. Presented here is a trinomial relationship method of relating percent loss due to self-absorption to filter mass loading, based on data reported by Higby, Luetzelschwab et al., Huang et al., Barnett et al., Smith et al., and Hogue et al. Under normal operating conditions at the stacks monitored by Effluent Management, the mass loading of sample filters averages 0.09 ± 0.12 (2s) mg/cm 2 (excluding negative values and outliers) and ranges from 0 mg/cm 2 to 0.24 mg/cm 2 . Based on current mass loading results for Effluent Management stack sample filters, the forced-zero trinomial relationship method estimated self-absorption losses of less than 5%. Because American National Standards Institute/Health Physics Society N13.1-2011 guidelines indicate a correction factor should be used when the penetration of radioactive material into the collection media or self-absorption of radiation by the material collected would reduce the count rate by more than 5%, it is possible continued application of a correction factor to the Effluent Management stack samples is no longer necessary. Nevertheless, continuing to assign a correction factor at the 5% threshold (i.e., 0.95) would be a conservative approach.

36 MATERIALS SCIENCE↗

Design-Filter Selection for H2 Control of Microgravity Isolation Systems: A Single-Degree-of-Freedom Case Study

Many microgravity space-science experiments require active vibration isolation, to attain suitably low levels of background acceleration for useful experimental results. The design of state-space controllers by optimal control methods requires judicious choices of frequency-weighting design filters. Kinematic coupling among states greatly clouds designer intuition in the choices of these filters, and the masking effects of the state observations cloud the process further. Recent research into the practical application of H2 synthesis methods to such problems, indicates that certain steps can lead to state frequency-weighting design-filter choices with substantially improved promise of usefulness, even in the face of these difficulties. In choosing these filters on the states, one considers their relationships to corresponding design filters on appropriate pseudo-sensitivity- and pseudo-complementary-sensitivity functions. This paper investigates the application of these considerations to a single-degree-of-freedom microgravity vibration-isolation test case. Significant observations that were noted during the design process are presented. along with explanations based on the existent theory for such problems.

Hampton, R. David↗