Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient filter”

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 73 records · Page 4

Systems and Methods for Peak-Seeking Control

A computerized system and method for peak-seeking-control that uses a unique Kalman filter design to optimize a control loop, in real time, to either maximize or minimize a performance function of a physical object ("plant"). The system and method achieves more accurate and efficient peak-seeking-control by using a time-varying Kalman filter to estimate both the performance function gradient (slope) and Hessian (curvature) based on direct position measurements of the plant, and does not rely upon modeling the plant response to persistent excitation. The system and method can be naturally applied in various applications in which plant performance functions have multiple independent parameters, and it does not depend upon frequency separation to distinguish between system dimensions.

Ryan, John J↗

Coming of Age: Polarization as a Probe of Plant Canopy Water Status

We tested the hypothesis that the relative water content (RWC) of the sunlit leaves in a plant canopy may be estimated from polarized canopy imagery. Recently (IGARSS, July 27-31, 2015, Milan, Italy), we reported the results of laboratory polarization measurements of single detached leaves during dry down. We found that RWC was linearly related to the ratio of the reflectance of the interior of the leaf and the leaf transmittance. Here we report application of the laboratory results to estimate RWC for sunlit leaves in a plant canopy. Using a commercial-off-the-shelf (COTS) Nikon 810 camera with Nikkor 300 mm lens and Polaroid type HN-22 linear polarizer, we photographed in the principle plane a plant canopy displaying a gradient of water stress and collected, at each of multiple points along the gradient, two images, one with the polarization filter oriented for maximum scene response and a second with the filter oriented for minimum scene response. We converted the digital values in the two images to reflectance factor with reference to images of a white, flat, horizontal Spectralon surface. We classified the polarization imagery, identifying reflecting leaves, transmitting leaves, other sunlit vegetation and shadows. For each image pair we normalized the leaf internal reflectance by dividing by the cosine of the angle of incidence of the sunlight on the leaf, selected the leaf maximum transmittance in the scene and divided to obtain the ratio reflectance/transmittance, which we compared with leaf RWC. We determined the leaf relative water content by harvesting a section of leaf and immediately placing it in a sealed container in an ice chest. Later in the laboratory the leaf sample was weighed, rehydrated, weighed, dried and again weighed. RWC was determined using the standard formula. Our experimental results support our hypothesis, suggesting that the RWC of sunlit leaves in a plant canopy may be estimated from analysis of polarization imagery collected by a COTS camera system. Unlike remotely sensed estimates of canopy equivalent water thickness, our estimates of the RWC of sunlit canopy leaves provide leaf physiological information. We propose RWC estimates based upon sunlit leaves are more relevant to assessing the water status of a plant canopy than would be RWC estimates based upon large FOV canopy measurements.

estimates↗

Absolute acceleration measurements on STS-50 from the Orbital Acceleration Research Experiment (OARE)

Orbital Acceleration Research Experiment (OARE) data on Space Transportation System (STS)-50 have been examined in detail during a 2-day time period. Absolute acceleration levels have been derived at the OARE location, the orbiter center-of-gravity, and at the STS-50 spacelab Crystal Growth Facility. During the interval, the tri-axial OARE raw telemetered acceleration measurements have been filtered using a sliding trimmed mean filter in order to remove large acceleration spikes (e.g., thrusters) and reduce the noise. Twelve OARE measured biases in each acceleration channel during the 2-day interval have been analyzed and applied to the filtered data. Similarly, the in situ measured x-axis scale factors in the sensor's most sensitive range were also analyzed and applied to the data. Due to equipment problem(s) on this flight, both y- and z-axis sensitive range scale factors were determined in a separate process using orbiter maneuvers and subsequently applied to the data. All known significant low-frequency corrections at the OARE location (i.e., both vertical and horizontal gravity-gradient, and rotational effects) were removed from the filtered data in order to produce the acceleration components at the orbiter center-of-gravity, which are the aerodynamic signals along each body axis. Results indicate that there is a force being applied to the Orbiter in addition to the aerodynamic forces. The OARE instrument and all known gravitational and electromagnetic forces have been reexamined, but none produces the observed effect. Thus, it is tentatively concluded that the orbiter is creating the environment observed. At least part of this force is thought to be due to the Flash Evaporator System.

Blanchard, Robert C.↗

Preliminary OARE absolute acceleration measurements on STS-50

On-orbit Orbital Acceleration Research Experiment (OARE) data on STS-50 was examined in detail during a 2-day time period. Absolute acceleration levels were derived at the OARE location, the orbiter center-of-gravity, and at the STS-50 spacelab Crystal Growth Facility. The tri-axial OARE raw acceleration measurements (i.e., telemetered data) during the interval were filtered using a sliding trimmed mean filter in order to remove large acceleration spikes (e.g., thrusters) and reduce the noise. Twelve OARE measured biases in each acceleration channel during the 2-day interval were analyzed and applied to the filtered data. Similarly, the in situ measured x-axis scale factors in the sensor's most sensitive range were also analyzed and applied to the data. Due to equipment problem(s) on this flight, both y- and z- axis sensitive range scale factors were determined in a separate process (using the OARE maneuver data) and subsequently applied to the data. All known significant low-frequency corrections at the OARE location (i.e., both vertical and horizontal gravity-gradient, and rotational effects) were removed from the filtered data in order to produce the acceleration components at the orbiter's center-of-gravity, which are the aerodynamic signals along each body axes. Results indicate that there is a force of unknown origin being applied to the Orbiter in addition to the aerodynamic forces. The OARE instrument and all known gravitational and electromagnetic forces were reexamined, but none produce the observed effect. Thus, it is tentatively concluded that the Orbiter is creating the environment observed.

Blanchard, Robert C.↗

Thermal Effects on Far-Field Distributed Acoustic Strain-Rate Sensors

Summary Fiber-optic cables cemented outside of the casing of an unconventional well measure crosswell strain changes during fracturing of neighboring wells with low-frequency distributed acoustic sensing (LF-DAS). As a hydraulic fracture intersects an observation well instrumented with fiber-optic cables, the fracture fluid injected at ambient temperatures can cool a section of the sensing fiber. Often, LF-DAS and distributed temperature sensing (DTS) cables are run in tandem, enabling the detection of such cooling events. The increasing use of LF-DAS for characterizing unconventional hydraulic fracture completions demands an investigation of the effects of temperature on the measured strain response by LF-DAS. Researchers have demonstrated that LF-DAS can be used to extract the temporal derivative of temperature for use as a differential-temperature-gradient sensor. However, differential-temperature-gradient sensing is predicated on the ability to filter strain components out of the optical signal. In this work, beginning with an equation for optical phase shift of LF-DAS signals, a model relating strain, temperature, and optical phase shift is explicitly developed. The formula provides insights into the relative strength of strain and temperature effects on the phase shift. The uncertainty in the strain-rate measurements due to thermal effects is estimated. The relationship can also be used to quantify uncertainties in differential-temperature-gradient sensors due to strain perturbations. Additionally, a workflow is presented to simulate the LF-DAS response accounting for both strain and temperature effects. Hydraulic fracture geometries are generated with a 3D fracture simulator for a multistage unconventional completion. The fracture width distributions are imported by a displacement discontinuity method (DDM) program to compute the strain rates along an observation well. An analytic model is used to approximate the temperature in the fracture. Using the derived formulae for optical phase shift, the model outputs are then used to compute the LF-DAS response at a fiber-optic cable, enabling the generation of waterfall plots including both strain and thermal effects. The model results suggest that before, during, and immediately following a fracture intersecting a well instrumented with fiber, the strain on the fiber drives the LF-DAS signal. However, at later times, as completion fluid cools the observation well, the temperature component of the LF-DAS signal can be equal to or exceed the strain component. The modeled results are compared to a published field case in an attempt to enhance the interpretation of LF-DAS waterfall plots. Finally, we propose a sensing configuration to identify the events when “wet fractures” (fractures with fluids) intersect the observation well.

Engineering↗

A Residuals Approach to Filtering, Smoothing and Identification for Static Distributed Systems

An approach for state estimation and identification of spatially distributed parameters embedded in static distributed (elliptic) system models is advanced. The method of maximum likelihood is used to find parameter values that maximize a likelihood functional for the system model, or equivalently, that minimize the negative logarithm of this functional. To find the minimum, a Newton-Raphson search is conducted that from an initial estimate generates a convergent sequence of parameter estimates. For simplicity, a Gauss-Markov approach is used to approximate the Hessian in terms of products of first derivatives. The gradient and approximate Hessian are computed by first arranging the negative log likelihood functional into a form based on the square root factorization of the predicted covariance of the measurement process. The resulting data processing approach, referred to here by the new term of predicted data covariance square root filtering, makes the gradient and approximate Hessian calculations very simple. A closely related set of state estimates is also produced by the maximum likelihood method: smoothed estimates that are optimal in a conditional mean sense and filtered estimates that emerge from the predicted data covariance square root filter.

Rodriguez, G.↗

Field-level reconstruction from foreground-contaminated 21-cm maps

Current and upcoming 21-cm experiments will soon be able to map 21-cm spatial fluctuations in three dimensions for a wide range of redshifts. However, bright foreground contamination and the nature of radio interferometry create significant challenges, making it difficult to access rich cosmological information from the Fourier modes that lie within the “foreground wedge”. Here, in this work, we introduce two approaches aiming to reconstruct the full 21-cm density field, including the missing modes in the wedge: (a) a field-level inference under an effective field theory (EFT) framework; (b) a diffusion-based deep generative model trained on simulations. Under the EFT framework, we implement a fully differentiable forward model that maps the initial conditions of matter fluctuations to the observed, foreground-filtered 21-cm maps. This enables a gradient-based sampler to simultaneously sample the initial conditions and bias parameters, allowing a physically motivated mode reconstruction. Alternatively, we apply a variational diffusion model to perform 21-cm density reconstruction at the map level. Our model is trained on semi-numerical simulations over a wide range of astrophysical parameters. Our results from both approaches should provide improved cosmological constraints from the field level and also enable cross-correlation between experiments that have little or no overlapping modes.

cosmological perturbation theory↗

Observations of cosmic rays in the outer heliosphere

Outward propagating cosmic ray modulation features seen in the outer heliosphere include Forbush decreases, the minimum of the 11 year cycle, and intervals of set recovery rate. The best theoretical model may not be the same for all these features. Some modulation features have been modified or have disappeared between Pioneer 11 and Pioneer 10 as if by the action of an interpalnetary low-pass filter. Azimuthal symmetry still applies, and the radial gradient for nucleons is greater than approximately 500 MeV nucleon is approximately 1 percent AU.

Fillius, Walker↗

Resolving experimental biases in the interpretation of diffusion experiments with a user-friendly numerical reactive transport approach

The reactive transport code CrunchClay was used to derive effective diffusion coefficients (D e ), clay porosities (ε), and adsorption distribution coefficients (K D ) from through-diffusion data while considering accurately the influence of unavoidable experimental biases on the estimation of these diffusion parameters. These effects include the presence of filters holding the solid sample in place, the variations in concentration gradients across the diffusion cell due to sampling events, the impact of tubing/dead volumes on the estimation of diffusive fluxes and sample porosity, and the effects of O-ring-filter setups on the delivery of solutions to the clay packing. Doing so, the direct modeling of the measurements of (radio)tracer concentrations in reservoirs is more accurate than that of data converted directly into diffusive fluxes. While the above-mentioned effects have already been described individually in the literature, a consistent modeling approach addressing all these issues at the same time has never been described nor made easily available to the community. A graphical user interface, CrunchEase, was created, which supports the user by automating the creation of input files, the running of simulations, and the extraction and comparison of data and simulation results. While a classical model considering an effective diffusion coefficient, a porosity and a solid/solution distribution coefficient (D e –ε–K D ) may be implemented in any reactive transport code, the development of CrunchEase makes it easy to apply by experimentalists without a background in reactive transport modeling. CrunchEase makes it also possible to transition more easily from a D e –ε–K D modeling approach to a state-of-the-art process-based understanding modeling approach using the full capabilities of CrunchClay, which include surface complexation modeling and a multi-porosity description of the clay packing with charged diffuse layers.

58 GEOSCIENCES↗

Understanding the role of segmentation on process-structure–property predictions made via machine learning

Here, the present study investigated the effect of porosity surface determination methods on performance of machine learning models used to predict the tensile properties of AlSi10Mg processed by laser powder bed fusion from micro-computed tomography data. Machine learning models applied in this work include support vector machines, neural networks, decision trees, and Bayesian classifiers. The effects of isosurface thresholding and local gradient approaches for porosity segmentation, as well as image filtering schemes, on model precision were evaluated for samples produced under differing levels of global energy density.

36 MATERIALS SCIENCE↗

Interaction of feel system and flight control system dynamics on lateral flying qualities

An investigation of feel system and flight control system dynamics on lateral flying qualities was conducted using the variable stability USAF NT-33 aircraft. Experimental variations in feel system natural frequency, force-deflection gradient, control system command architecture type, flight control system filter frequency, and control system delay were made. The experiment data include pilot ratings using the Cooper-Harper (1969) rating scale, pilot comments, and tracking performance statistic. Three test pilots served as evaluators. The data indicate that as the feel system natural frequency is reduced lateral flying qualities degrade. At the slowest feel system frequency, the closed-loop response becomes nonlinear with a 'bobweight' effect apparent in the feel system. Feel system influences were essentially independent of the control system architecture. The flying qualities influence due to the feel system was different than when the identical dynamic systenm was used as a flight control system element.

Bailey, Randall E.↗

Three-axis attitude determination via Kalman filtering of magnetometer data

A three-axis Magnetometer/Kalman Filter attitude determination system for a spacecraft in low-altitude Earth orbit is developed, analyzed, and simulation tested. The motivation for developing this system is to achieve light weight and low cost for an attitude determination system. The extended Kalman filter estimates the attitude, attitude rates, and constant disturbance torques. Accuracy near that of the International Geomagnetic Reference Field model is achieved. Covariance computation and simulation testing demonstrate the filter's accuracy. One test case, a gravity-gradient stabilized spacecraft with a pitch momentum wheel and a magnetically-anchored damper, is a real satellite on which this attitude determination system will be used. The application to a nadir pointing satellite and the estimation of disturbance torques represent the significant extensions contributed by this paper. Beyond its usefulness purely for attitude determination, this system could be used as part of a low-cost three-axis attitude stabilization system.

Martel, Francois↗

Infrared Rugates by Molecular Beam Epitaxy

Rugates are optical structures that have a sinusoidal index of refraction (harmonic gradient-index field). As their discrete high/ low index filter counterparts, they can be used as narrow rejection band filters. However, since rugates do not have abrupt interfaces, they tend to have a smaller absorption, hence deliver a higher in band reflectivity. The absence of sharp interfaces makes rugates even more desirable for high-energy narrow band reflectors. In this application, the lack of a sharp interface at the maximum internal standing wave electric field results in higher breakdown strengths. Our method involves fabricating rugates, with molecular beam epitaxy, on GaAs wafers as an Al(x)Ga(1-x)As single-crystal film.

Rona, M.↗

Clutter filter design considerations for Airborne Doppler radar detection of windshear

The problem of clutter rejection when processing down-looking Doppler radar returns from a low altitude airborne platform is a paramount problem. With radar as a remote sensor for detecting and predicting windshear in the vicinity of an urban airport, dynamic range requirements can exceed 50 dB because of high clutter to signal ratios. This presentation describes signal processing considerations in the presence of distributed and/or discrete clutter interference. Previous analyses have considered conventional range cell processing of radar returns from a rigidly mounted radar platform using either the Fourier or the pulse-pair method to estimate average windspeed and windspeed variation within a cell. Clutter rejection has been based largely upon analyzing a particular environment in the vicinity of the radar and employing a variety of techniques to reduce interference effects including notch filtering, Fourier domain line editing, and use of clutter maps. For the airborne environment the clutter characteristics may be somewhat different. Conventional clutter rejection methods may have to be changed and new methods will probably be required to provide useful signal to noise ratios. Various considerations are described. A major thrust has been to evaluate the effect of clutter rejection filtering upon the ability to derive useful information from the post-filter radar data. This analysis software is briefly described. Finally, some ideas for future analysis are considered including the use of adaptive filtering for clutter rejection and the estimation of windspeed spatial gradient directly from radar returns as a means of reducing the effects of clutter on the determination of a windshear hazard.

Baxa, Ernest G., Jr.↗

A Priori Analyses of Three Subgrid-Scale Models for One-Parameter Families of Filters

The decay of isotropic turbulence a compressible flow is examined by direct numerical simulation (DNS). A priori analyses of the DNS data are then performed to evaluate three subgrid-scale (SGS) models for large-eddy simulation (LES): a generalized Smagorinsky model (M1), a stress-similarity model (M2), and a gradient model (M3). The models exploit one-parameter second- or fourth-order filters of Pade type, which permit the cutoff wavenumber k(sub c) to be tuned independently of the grid increment (delta)x. The modeled (M) and exact (E) SGS-stresses are compared component-wise by correlation coefficients of the form C(E,M) computed over the entire three-dimensional fields. In general, M1 correlates poorly against exact stresses (C < 0.2), M3 correlates moderately well (C approx. 0.6), and M2 correlates remarkably well (0.8 < C < 1.0). Specifically, correlations C(E, M2) are high provided the grid and test filters are of the same order. Moreover, the highest correlations (C approx.= 1.0) result whenever the grid and test filters are identical (in both order and cutoff). Finally, present results reveal the exact SGS stresses obtained by grid filters of differing orders to be only moderately well correlated. Thus, in LES the model should not be specified independently of the filter.

Pruett, C. David↗

Full Gradient Solution to Adaptive Hybrid Control

This paper focuses on the adaptation mechanisms in adaptive hybrid controllers. Most adaptive hybrid controllers update two filters individually according to the filtered-reference least mean squares (FxLMS) algorithm. Because this algorithm was derived for feedforward control, it does not take into account the presence of a feedback loop in the gradient calculation. This paper provides a derivation of the proper weight vector gradient for hybrid (or feedback) controllers that takes into account the presence of feedback. In this formulation, a single weight vector is updated rather than two individually. An internal model structure is assumed for the feedback part of the controller. The full gradient is equivalent to that used in the standard FxLMS algorithm with the addition of a recursive term that is a function of the modeling error. Some simulations are provided to highlight the advantages of using the full gradient in the weight vector update rather than the approximation.

Bean, Jacob↗

Full Gradient Solution to Adaptive Hybrid Control

This paper focuses on the adaptation mechanisms in adaptive hybrid controllers. Most adaptive hybrid controllers update two filters individually according to the filtered reference least mean squares (FxLMS) algorithm. Because this algorithm was derived for feedforward control, it does not take into account the presence of a feedback loop in the gradient calculation. This paper provides a derivation of the proper weight vector gradient for hybrid (or feedback) controllers that takes into account the presence of feedback. In this formulation, a single weight vector is updated rather than two individually. An internal model structure is assumed for the feedback part of the controller. The full gradient is equivalent to that used in the standard FxLMS algorithm with the addition of a recursive term that is a function of the modeling error. Some simulations are provided to highlight the advantages of using the full gradient in the weight vector update rather than the approximation.

Bean, Jacob↗

Atmospheric model development in support of SEASAT. Volume 3: Prediction models

As part of the SEASAT program of NASA, a set of four hemispheric, atmospheric models are developed. The four models use a polar stereographic grid in the horizontal and a sigma coordinate in the vertical. Conservation forms of the difference equations are integrated using either a fifteen or four minute time step on a 381 km or 127 km grid (at 60 deg N) for the 63 x 63 or 187 x 187 models respectively. A nonlinear pressure smoother is used along with momentum and temperature diffusion to help control computational noise. The horizontal boundary conditions are insulated slippery walls. Centered time differencing with time averaging of the pressure gradient force term in the momentum equations is used. Robert time filtering on of the temperature and moisture solutions is used for computational stability. The moisture and heat source/sink terms are modeled. Terms representing evaporation and large scale condensation, sensible heat exchange, parameterized cumulus convection and precipitation, and solar and terrestrial radiation are included. Initialization of the models is based on a pattern conservation technique to obtain objective analysis of the state parameter structure from the surface to 50 mb.

Wellck, R. E.↗