SEARCH · Search NASA

Results for “polynomial filtering”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A DATA EFFICIENT SPARSE MODELING FRAMEWORK FOR POWER ESTIMATION IN WATER TREATMENT SENSING OPERATIONS

With increasing freshwater scarcity, advanced process design mechanisms such as Closed-Circuit Reverse Osmosis (CCRO) and Digital/Physical Twin systems are gaining traction in water treatment and reuse operations. While digital and physical twin models enable improved system insight and control, their development is often expensive and computationally intensive, requiring large volumes of synthetic or experimental data to characterize underlying process dynamics. This work introduces a sparse surrogate modeling framework to estimate power consumption from measured flow and pressure variables, along with their nonlinear polynomial and interaction expansions. To ensure model reliability and reduce overfitting, a two-stage pipeline is proposed. First, a dynamic data filtering algorithm is employed to remove uninformative observations and transient operational states. Second, a sparse penalized regression technique is applied to select a minimal set of parsimonious features. The proposed model achieves high sparsity, retaining only 7 out of 34 candidate features (≈79.41% sparsity) while delivering a root mean square error (RMSE) of 0.072 on the test dataset.

Mukherjee, Subrata [ORNL] (ORCID:0000000309930338)↗

NASA NTRS1987

SIR-B ocean-wave enhancement with fast Fourier transform techniques

Shuttle Imaging Radar (SIR-B) imagery is Fourier filtered to remove the estimated system-transfer function, reduce speckle noise, and produce ocean scenes with a gray scale that is proportional to wave height. The SIR-B system response to speckled scenes of uniform surfaces yields an estimate of the stationary wavenumber response of the imaging radar, modeled by the 15 even terms of an eighth-order two-dimensional polynomial. Speckle can also be used to estimate the dynamic wavenumber response of the system due to surface motion during the aperture synthesis period, modeled with a single adaptive parameter describing an exponential correlation along track. A Fourier filter can then be devised to correct for the wavenumber response of the remote sensor and scene correlation, with subsequent subtraction of an estimate of the speckle noise component. A linearized velocity bunching model, combined with a surface tilt and hydrodynamic model, is incorporated in the Fourier filter to derive estimates of wave height from the radar intensities corresponding to individual picture elements.

Tilley, David G.↗

NASA NTRS1985

Solution of the two-dimensional spectral factorization problem

An approximation theorem is proven which solves a classic problem in two-dimensional (2-D) filter theory. The theorem shows that any continuous two-dimensional spectrum can be uniformly approximated by the squared modulus of a recursively stable finite trigonometric polynomial supported on a nonsymmetric half-plane.

Lawton, W. M.↗

NASA NTRS1991

How to observe the gyre to global-scale variability in satellite altimetry - Signal attenuation by orbit error removal

Formulas analogous to the frequency response functions for commonly used filters in orbit error removal are analytically derived to devise observational strategies for the large-scale oceanic variability and to decipher the signal contents of previous results. These include the polynomial orbit error approximations, i.e., the linear, bias-only and quadratic corrections, and the sinusoidal orbit error approximations (the purely sinusoidal correction, and the sinusoid-and-bias correction). It is shown that the frequency response function for a polynomial correction is a function of the ratio of wavelength/track length and to retain 90 percent or more of the signal at a certain wavelength, the ratio must be less than 0.65 (for the quadratic case), 0.90 (linear), and 1.54 (bias-only).

Tai, Chang-Kou↗

NASA NTRSDate not supplied

Field-Programmable Gate Array Implementation of a Single Photon-Counting Receive Modem

We present a field-programmable gate array (FPGA) implementation of a single photon-counting receive modem for a pulse position modulated signal. The modem is compliant with the Consultative Committee for Space Data Systems (CCSDS) High Photon Efficiency (HPE) Optical Communications Coding and Synchronization standard and is capable of a maximum data rate of 267 Mbps. The system is designed on a commercial off-the-shelf FPGA platform and utilizes superconducting nanowire single photon counting detectors, analog to digital converters (ADC s) to sample the detectors, and two FPGAs. Symbol timing recovery, photon counting, convolutional deinterleaving, and codeword synchronization areis performed in the first FPGA. The second FPGA performs iterative decoding on each codeword of the serially concatenated pulse position modulated (SCPPM) signal. A digital filter is included to compensate for timing jitter of the detector, and the decoder throughput can be modified adjusted through reconfigurable parallelization. The decoder also implements a resource-efficient, algorithmic polynomial interleaver and deinterleaver. Both FPGAs can be reconfigured to switch between pulse position modulation (PPM)-16 and PPM-32 with code rates 1/3, 1/2, and 2/3. In this paper, we describe the receiver architecture and FPGA implementation of the timing recovery loop and SCPPM decoder, FPGA utilization for the different modes, and receive modem characterization test results.

Field-programmable Gate Array↗

DOE OSTI2022

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

47 Tuc in Rubin Data Preview 1. Exploring Early LSST Data and Science Potential

We present analyses of the early data from Rubin Observatory’s Data Preview 1 (DP1) for the field of the globular cluster 47 Tuc. The DP1 data set for 47 Tuc includes four nights of observations from the Rubin Commissioning Camera (LSSTComCam), covering multiple bands (ugriy). We address challenges of crowding in the inner region of the cluster and toward the SMC in DP1, and demonstrate improved star–galaxy separation by fitting fifth-degree polynomials to the stellar loci in color–color diagrams and applying multidimensional sigma clipping. We compile a catalog of 3576 probable 47 Tuc member stars selected via a combination of isochrone, Gaia proper-motion, and color–color space matched filtering. We explore the sources of photometric scatter in the 47 Tuc color–color sequence, evaluating contributions from various potential sources, including differential extinction within the cluster. Finally, of the 72 well-characterized variables in the field, we recover three known variable stars, including two RR Lyrae and one eclipsing binary, in the coadd-based object catalog, and identify 62 in the difference image-based object catalog. Although the DP1 lightcurves have sparse temporal sampling, they appear to follow the patterns of densely sampled literature lightcurves well. Despite some data limitations for crowded-field stellar analysis, DP1 demonstrates the promising scientific potential for future LSST data releases.

Choi, Yumi [NSF National Optical-Infrared Astronom↗

NASA NTRS1991

A new technique for enlargement and reconstruction of digital sensor imagery

A new image enlargement and reconstruction routine, the digital image enlarging balanced reconstruction algorithm (DIEBRA), is described. Using a highly modified form of 'balanced' two-dimensional polynomial interpolation, this program enlarges digital imagery, creating synthetic high spatial resolution images. Statistical analysis shows the DIEBRA-generated imagery to be significantly closer to true high spatial resolution imagery at all frequencies than imagery generated by a cubic convolutional reconstruction filter.

Rhea, W. Joseph↗

NASA NTRS1978

Algorithms for isolating worst case systematic data errors

Two separate algorithms are derived for testing filter sensitivity to systematic data errors. One algorithm provides the absolute minimum Euclidean norm data error for a given estimate component error. The second algorithm can be used to find the minimum norm data error which can be generated by restricted degree Legendre polynomials. A specific very long baseline interferometry (VLBI) baseline estimation is analyzed with the algorithm. It is found that the local vertical is the most sensitive component to error in the data space. The efficiency of a data error sequence linear in elevation angle is within 7% that of the absolute worst case sequence. Elevation angle dependent errors are explored and the special case of a mismodeled troposphere is treated.

Curkendall, D. W.↗

NASA NTRS1984

The importance of nonlinearity on the higher harmonic control of helicopter vibration

The effect of nonlinearity on the higher harmonic control (HHC) of helicopter vibration is investigated using a nonlinear aeroelastic simulation. A nonlinear solution is proposed which relates the HHC inputs to vibration outputs on the basis of a Volterra functional series. The Volterra series solution is shown to reduce to a vector polynomial equation relating HHC inputs to vibration outputs at any harmonic frequency. The nonlinear transfer relationship is identified from a nonlinear vibration analysis computer program, and the identification model is examined in detail. Improvements to current HHC algorithms are presented, and several Kalman filter divergence problems are quantified.

Molusis, J. A.↗

NASA NTRS1982

Radar image processing for rock-type discrimination

Image processing and enhancement techniques for improving the geologic utility of digital satellite radar images are reviewed. Preprocessing techniques such as mean and variance correction on a range or azimuth line by line basis to provide uniformly illuminated swaths, median value filtering for four-look imagery to eliminate speckle, and geometric rectification using a priori elevation data. Examples are presented of application of preprocessing methods to Seasat and Landsat data, and Seasat SAR imagery was coregistered with Landsat imagery to form composite scenes. A polynomial was developed to distort the radar picture to fit the Landsat image of a 90 x 90 km sq grid, using Landsat color ratios with Seasat intensities. Subsequent linear discrimination analysis was employed to discriminate rock types from known areas. Seasat additions to the Landsat data improved rock identification by 7%.

Blom, R. G.↗

NASA NTRS1977

Undulation and anomaly estimation using Geos-3 altimeter data without precise satellite orbits

Geos-3 altimeter data was used to obtain an approximate geoid undulation which was contaminated by long wavelength errors caused primarily by altimeter bias and orbit error. This error was reduced by fitting, with a low degree polynomial, the raw undulation data to the undulations implied by the GEM 7 potential coefficients in an adjustment process that included conditions on tracks that cross. These adjusted undulations were used to construct a geoid map in the Geos-3 calibration area using a least squares filter to remove the remaining noise in the undulations.

Rummel, R.↗

NASA NTRS1990

Inflight estimation of gyro noise

A method is described and demonstrated for estimating single-axis gyro noise levels in terms of the Farrenkopf model parameters. This is accomplished for the Cosmic Background Explorer (COBE) by comparing gyro-propagated attitudes with less accurate single-frame solutions and fitting the squared differences to a third-order polynomial in time. Initial results are consistent with the gyro specifications, and these results are used to determine limits on the duration of batches used to determine attitude. Sources of error are discussed, and guidelines for a more elegant implementation, as part of a batch estimator or filter, are included for future work.

Filla, O. H.↗

NASA NTRS1981

Optimal spatial filtering and transfer function for SAR ocean wave spectra

The impulse response of the SAR system is not a delta function and the spectra represent the product of the underlying image spectrum with the transform of the impulse response which must be removed. A digitally computed spectrum of SEASAT imagery of the Atlantic Ocean east of Cape Hatteras was smoothed with a 5 x 5 convolution filter and the trend was sampled in a direction normal to the predominant wave direction. This yielded a transform of a noise-like process. The smoothed value of this trend is the transform of the impulse response. This trend is fit with either a second- or fourth-order polynomial which is then used to correct the entire spectrum. A 16 x 16 smoothing of the spectrum shows the presence of two distinct swells. Correction of the effects of speckle is effected by the subtraction of a bias from the spectrum.

Beal, R. C.↗

NASA NTRS1991

A fast algorithm for control and estimation using a polynomial state-space structure

One of the major problems associated with the control of flexible structures is the estimation of system states. Since the parameters of the structures are not constant under varying loads and conditions, conventional fixed parameter state estimators can not be used to effectively estimate the states of the system. One alternative is to use a state estimator which adapts to the condition of the system. One such estimator is the Kalman filter. This filter is a time varying recursive digital filter which is based upon a model of the system being measured. This filter adapts the model according to the output of the system. Previously, the Kalman filter has only been used in an off-line capacity due to the computational time required for implementation. With recent advances in computer technology, it is becoming a viable tool for use in the on-line environment. A distributed Kalman filter implementation is described for fast estimation of the state of a flexible arm. A key issue, is the sensor structure and initial work on a distributed sensor that could be used with the Kalman filter is presented.

Shults, James R.↗

DOE OSTI2024

Testing the LSST Difference Image Analysis Pipeline Using Synthetic Source Injection Analysis

Abstract We evaluate the performance of the Legacy Survey of Space and Time Science Pipelines Difference Image Analysis (DIA) on simulated images. By adding synthetic sources to galaxies on images, we trace the recovery of injected synthetic sources to evaluate the pipeline on images from the Dark Energy Science Collaboration Data Challenge 2. The pipeline performs well, with efficiency and flux accuracy consistent with the signal-to-noise ratio of the input images. We explore different spatial degrees of freedom for the Alard–Lupton polynomial-Gaussian image subtraction kernel and analyze for trade-offs in efficiency versus artifact rate. Increasing the kernel spatial degrees of freedom reduces the artifact rate without loss of efficiency. The flux measurements with different kernel spatial degrees of freedom are consistent. We also here provide a set of DIA flags that substantially filter out artifacts from the DIA source table. We explore the morphology and possible origins of the observed remaining subtraction artifacts and suggest that given the complexity of these artifact origins, a convolution kernel with a set of flexible bases with spatial variation may be needed to yield further improvements.

Liu, S. (ORCID:0000000244612143)↗

NASA NTRS1994

R(sub nu)-dependent optical and near-ultraviolet extinction

We have derived extinctions A(lambda)/A(V) at the wavelengths of the uvby filters for 22 stars, with a range of values of R(sub nu), from the sample of Cardelli, Clayton, & Mathis (1989, hereafter CCM). We have fit these extinctions, and also UBVRIJHKL, IUE and ANS extinction measurements, with linear relations A(lambda)/A/(V) = a+b/R(sub nu) and fit a and b as a function of x(=1/lambda) with polynomials to obtain an R(sub nu)-dependent mean extinction law (A(x)/A(V) = a(x) + b(x)/R(sub nu))in the optical and near-ultraviolet (1.1/micrometer less than or equal to 3.3/micrometer). This law is virtually identical to the CCM extinction law for large values of R(sub nu)(R(sub nu) approximately 5) but is slightly lower in the near-ultraviolet for smaller R(sub nu) (R(sub nu) approximately 3). The extinction law presented here agrees much better with a high-resolution extinction curve for the diffuse interstellar medium (R(sub nu) approximately 3.1), presented by Bastiaansen (1992), than CCM. The deviations of individual extinction curves from the mean are dominated by observational errors. The wavelength resolution of this work is not high enough to show evidence for or against the existence of very broad structure in optical extinction curves.

O'Donnell, James E.↗

NASA NTRS2000

Global Design Optimization for Aerodynamics and Rocket Propulsion Components

Modern computational and experimental tools for aerodynamics and propulsion applications have matured to a stage where they can provide substantial insight into engineering processes involving fluid flows, and can be fruitfully utilized to help improve the design of practical devices. In particular, rapid and continuous development in aerospace engineering demands that new design concepts be regularly proposed to meet goals for increased performance, robustness and safety while concurrently decreasing cost. To date, the majority of the effort in design optimization of fluid dynamics has relied on gradient-based search algorithms. Global optimization methods can utilize the information collected from various sources and by different tools. These methods offer multi-criterion optimization, handle the existence of multiple design points and trade-offs via insight into the entire design space, can easily perform tasks in parallel, and are often effective in filtering the noise intrinsic to numerical and experimental data. However, a successful application of the global optimization method needs to address issues related to data requirements with an increase in the number of design variables, and methods for predicting the model performance. In this article, we review recent progress made in establishing suitable global optimization techniques employing neural network and polynomial-based response surface methodologies. Issues addressed include techniques for construction of the response surface, design of experiment techniques for supplying information in an economical manner, optimization procedures and multi-level techniques, and assessment of relative performance between polynomials and neural networks. Examples drawn from wing aerodynamics, turbulent diffuser flows, gas-gas injectors, and supersonic turbines are employed to help demonstrate the issues involved in an engineering design context. Both the usefulness of the existing knowledge to aid current design practices and the need for future research are identified.

Shyy, Wei↗