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 91 records · Page 5

Remote sensing of ice phenomena from orbit by signal correlation of multiple receiver responses

The method of signal correlation of microwave responses as applied to the measurement of Earth-surface ice temperatures from orbit is explained and summarized. Ice temperatures are estimated by a correlation function that is derived from the processes of a forward stepwise correlator. Subsets of the post-detected outputs of microwave receiving channels are combined in a multivariate cross-correlation function which operates as a spatial filter and serves to improve the spatial resolution of the thermal gradients in ice structures. The correlator is designed to selectively identify the correlative components among the microwave responses and to strongly suppress or cancel the non-correlative components appearing in the post-detected outputs.

Stacey, J. M.↗

AIRS Observations of DomeC in Antarctica and Comparison with Automated Weather Stations (AWS)

We compare the surface temperatures at Dome Concordia (DomeC) deduced from AIRS data and two Automatic Weather Stations at Concordia Station: AWS8989 , which has been in operation since December 1996, and AWS.it, for which data are available between January and November 2005. The AWS8989 readings are on average 3 K warmer than the AWS.it readings, with a warmer bias in the Antarctic summer than in the winter season. Although AIRS measures the skin brightness temperature, while the AWS reports the temperature of the air at 3 meter above the surface, the AIRS measurements agree well with the AWS.it readings for all data and separately for the summer and winter seasons, if data taken in the presence of strong surface inversions are filtered out. This can be done by deducing the vertical temperature gradient above the surface directly from the AIRS temperature sounding channels or indirectly by noting that extreme vertical gradients near the surface are unlikely if the wind speed is more than a few meters per second. Since the AIRS measurements are very well calibrated, the agreement with AWS.it is very encouraging. The warmer readings of AWS8989 are likely due to thermal contamination of the AWS8989 site by the increasing activity at Concordia Station. Data from an AWS.it quality station could be used for the evaluation of radiometric accuracy and stability of polar orbiting sounders at low temperatures. Unfortunately, data from AWS.it was available only for a limited time. The thermal contamination of the AWS8989 data makes long-term trends deduced from AWS8989 and possibly results about the rapid Antarctic warming deduced from other research stations on Antarctica suspect. AIRS is the first hyperspectral infrared sounder designed in support of weather forecasting and climate research. It was launched in May 2002 on the EOS Aqua spacecraft into a 704 km altitude polar sun-synchronous orbit. The lifetime of AIRS, estimated before launch to be at least 5 years is, based on the latest evaluation, limited by the amount of attitude control gas on the EOS Aqua spacecraft, which is expected to last through 2015.

atmospheric↗

Optimization of OT-MACH Filter Generation for Target Recognition

An automatic Optimum Trade-off Maximum Average Correlation Height (OT-MACH) filter generator for use in a gray-scale optical correlator (GOC) has been developed for improved target detection at JPL. While the OT-MACH filter has been shown to be an optimal filter for target detection, actually solving for the optimum is too computationally intensive for multiple targets. Instead, an adaptive step gradient descent method was tested to iteratively optimize the three OT-MACH parameters, alpha, beta, and gamma. The feedback for the gradient descent method was a composite of the performance measures, correlation peak height and peak to side lobe ratio. The automated method generated and tested multiple filters in order to approach the optimal filter quicker and more reliably than the current manual method. Initial usage and testing has shown preliminary success at finding an approximation of the optimal filter, in terms of alpha, beta, gamma values. This corresponded to a substantial improvement in detection performance where the true positive rate increased for the same average false positives per image.

false alarm reduction↗

Improved detection of low-frequency within-host variants from deep sequencing: A case study with human papillomavirus

Abstract High-coverage sequencing allows the study of variants occurring at low frequencies within samples, but is susceptible to false-positives caused by sequencing error. Ion Torrent has a very low single nucleotide variant (SNV) error rate and has been employed for the majority of human papillomavirus (HPV) whole genome sequences. However, benchmarking of intrahost SNVs (iSNVs) has been challenging, partly due to limitations imposed by the HPV life cycle. We address this problem by deep sequencing three replicates for each of 31 samples of HPV type 18 (HPV18). Errors, defined as iSNVs observed in only one of three replicates, are dominated by C→T (G→A) changes, independently of trinucleotide context. True iSNVs, defined as those observed in all three replicates, instead show a more diverse SNV type distribution, with particularly elevated C→T rates in CCG context (CCG→CTG; CGG→CAG) and C→A rates in ACG context (ACG→AAG; CGT→CTT). Characterization of true iSNVs allowed us to develop two methods for detecting true variants: (1) VCFgenie, a dynamic binomial filtering tool which uses each variant’s allele count and coverage instead of fixed frequency cut-offs; and (2) a machine learning binary classifier which trains eXtreme Gradient Boosting models on variant features such as quality and trinucleotide context. Each approach outperforms fixed-cut-off filtering of iSNVs, and performance is enhanced when both are used together. Our results provide improved methods for identifying true iSNVs in within-host applications across sequencing platforms, specifically using HPV18 as a case study.

Virology↗

High-gradient magnetic separation of colloidal uranium oxide particles from soil components in aqueous suspensions

The separation of uranium oxide (UO 2 ) particles from soil-surrogate particles in aqueous suspensions was achieved using filtration enhanced by a magnetic field. Enhanced attraction of paramagnetic UO 2 colloids to a ferromagnetic stainless-steel filter placed in a strong magnetic field arises because of the positive magnetic susceptibility of the particles and the high-gradient field generated near ferromagnetic fibers. Enhanced uptake of smaller particles over larger ones occurs through Brownian motion that promotes the collision of particles with the ferromagnetic fibers of the filter. Hence, this work focused on UO 2 particles in the colloidal size range. Experiments used a water-cooled electromagnet and an array of permanent magnets. Chemical analysis showed that the magnetic field increased the capture efficiency of uranium particles from a range of 27–53% with the magnet off up to 98% with the magnet on after a single pass of the suspension through the filter. Further, the recovery of the UO 2 particles from the filter, however, was more difficult to achieve. Small amounts of UO 2 , together with significant amounts of background SiO 2 particles, were removed from the filter during a first flush with the magnetic field on. A much larger recovery of UO 2 was not observed until a second out-of-field flush was performed, which also released some SiO 2 . The degree to which particle separation was enhanced through the use of multi-stage filtration compared to single pass-through filtration was also examined. A design was suggested that could be used to optimize the separation efficiency for a continuous process.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Outputs of paired Gabor filters summed across the background frame of reference predict the direction of movement

A cortical neural network that computes the visibility of shifts in the direction of movement is proposed. The network computes: (1) the magnitude of the position difference between the test and background patterns, (2) localized contrast differences at different spatial scales analyzed by computing temporal gradients of the difference and sum of the outputs of paired even- and odd-symmetric bandpass filters convolved with the input pattern, and (3) using global processes that pool the output from paired even- and odd-symmetric simple and complex cells across the spatial extent of the background frame of reference the direction a test pattern moved relative to a textured background. Evidence that magnocellular pathways are used to discriminate the direction of movement is presented. Since magnocellular pathways are used to discriminate the direction of movement, this task is not affected by small pattern changes such as jitter, short presentations, blurring, and different background contrasts that result when the veiling illumination in a scene changes.

Lawton, Teri B.↗

Atmospheric Modeling And Sensor Simulation (AMASS) study

A 4800 band synchronous communications link was established between the Perkin-Elmer (P-E) 3250 Atmospheric Modeling and Sensor Simulation (AMASS) system and the Cyber 205 located at the Goddard Space Flight Center. An extension study of off-the-shelf array processors offering standard interface to the Perkin-Elmer was conducted to determine which would meet computational requirements of the division. A Floating Point Systems AP-120B was borrowed from another Marshall Space Flight Center laboratory for evaluation. It was determined that available array processors did not offer significantly more capabilities than the borrowed unit, although at least three other vendors indicated that standard Perkin-Elmer interfaces would be marketed in the future. Therefore, the recommendation was made to continue to utilize the 120B ad to keep monitoring the AP market. Hardware necessary to support requirements of the ASD as well as to enhance system performance was specified and procured. Filters were implemented on the Harris/McIDAS system including two-dimensional lowpass, gradient, Laplacian, and bicubic interpolation routines.

Parker, K. G.↗

Part 2 of a Computational Study of a Drop-Laden Mixing Layer

This second of three reports on a computational study of a mixing layer laden with evaporating liquid drops presents the evaluation of Large Eddy Simulation (LES) models. The LES models were evaluated on an existing database that had been generated using Direct Numerical Simulation (DNS). The DNS method and the database are described in the first report of this series, Part 1 of a Computational Study of a Drop-Laden Mixing Layer (NPO-30719), NASA Tech Briefs, Vol. 28, No.7 (July 2004), page 59. The LES equations, which are derived by applying a spatial filter to the DNS set, govern the evolution of the larger scales of the flow and can therefore be solved on a coarser grid. Consistent with the reduction in grid points, the DNS drops would be represented by fewer drops, called computational drops in the LES context. The LES equations contain terms that cannot be directly computed on the coarser grid and that must instead be modeled. Two types of models are necessary: (1) those for the filtered source terms representing the effects of drops on the filtered flow field and (2) those for the sub-grid scale (SGS) fluxes arising from filtering the convective terms in the DNS equations. All of the filtered-sourceterm models that were developed were found to overestimate the filtered source terms. For modeling the SGS fluxes, constant-coefficient Smagorinsky, gradient, and scale-similarity models were assessed and calibrated on the DNS database. The Smagorinsky model correlated poorly with the SGS fluxes, whereas the gradient and scale-similarity models were well correlated with the SGS quantities that they represented.

Okongo, Nora↗

Treating convection in sequential solvers

The treatment of the convection terms in the sequential solver, a standard procedure found in virtually all pressure based algorithms, to compute the flow problems with sharp gradients and source terms is investigated. Both scalar model problems and one-dimensional gas dynamics equations have been used to study the various issues involved. Different approaches including the use of nonlinear filtering techniques and adoption of TVD type schemes have been investigated. Special treatments of the source terms such as pressure gradients and heat release have also been devised, yielding insight and improved accuracy of the numerical procedure adopted.

Shyy, Wei↗

Physical processes in comets

Post-Halley comets are known to be irregular objects with most nucleus activity very localized and with the dust coma capable of fragmentation and apparently being a source of gas. Older, 1-D strategies which assume steady isotropic outflow of material can give poor time-and-space averaged results, at best. With 2-D data, images through interference filters, it is hoped that: dust structures can be seen that give evidence of the proper geometry for data reduction; gradients along the axes of symmetry and evidence of fragmentation can be studied; and that evidence for gas abundance gradients associated with the dust can be found. High quality data from brighter comets can then be used to suggest improved data reduction procedures for fainter ones. To obtain such data, large image-quality interference filters were procured for use with a charge coupled device (CCD) camera at Lick Observatory, where the scale of the 1 m Nickel reflector is ideal for brighter comets. Whenever possible, data is taken simultaneously with other telescopes and equipment, especially spectroscopy at the Lick 3 m or infrared photometry at the Infrared Telescope Facility (IRTF) on Mauna Kea.

Newburn, Ray L., Jr.↗

A Priori Subgrid Scale Modeling for a Droplet Laden Temporal Mixing Layer

Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using a direct numerical simulation (DNS) database. The DNS is for a Reynolds number (based on initial vorticity thickness) of 600, with droplet mass loading of 0.2. The gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. Since Large Eddy Simulation (LES) of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be given by the filtered variables plus a correction based on the filtered standard deviation, which can be computed from the sub-grid scale (SGS) standard deviation. This model predicts unfiltered variables at droplet locations better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: Smagorinsky, gradient and scale-similarity. When properly calibrated, the gradient and scale-similarity methods give results in excellent agreement with the DNS.

Okongo, Nora↗

Procedures for Tuning a Multiresonator Photonic Filter

Two procedures have been devised for tuning a photonic filter that comprises multiple whispering-gallery mode (WGM) disk resonators. As used here, tuning signifies locking the filter to a specific laser frequency and configuring the filter to obtain a desired high-order transfer function. The main problem in tuning such a filter is how to select the correct relative loading of the resonators to realize a prescribed filter function. The first of the two procedures solves this problem. As temperature gradients develop during operation, the spectra of individual resonators tend to drift, primarily because of the thermorefractive effect. Thus, there arises the additional problem of how to adjust the tuning during operation to maintain the desired transfer function. The second of the two procedures solves this problem. To implement the procedures, it is necessary to incorporate the resonators into an apparatus like that of Figure 1. In this apparatus, the spectrum of each resonator can be adjusted individually, via the electro-optical effect, by adjusting a bias voltage applied to that resonator. In addition, the positions of the coupling prisms and resonators can be adjusted to increase or reduce the gaps between them, thereby reducing or increasing, respectively, the optical coupling between them. The optical power (Pi) in resonator i is monitored by use of a tracking photodiode. Another tracking diode monitors the power reflected from the input terminal (Pr), and still others monitor the input power (Pin) and output power (Po). The readings of these photodiodes are used to guide the tuning adjustments described in this paper.

Matsko, Andrey↗

Embedded training of neural-network subgrid-scale turbulence models

We report that the weights of a deep neural-network model are optimized in conjunction with the governing flow equations to provide a model for subgrid-scale stresses in a temporally developing plane turbulent jet at Reynolds number Re 0 = 6000 . The objective function for training is first based on the instantaneous filtered velocity fields from a corresponding direct numerical simulation, and the training is by a stochastic gradient descent method, which uses the adjoint Navier-Stokes equations to provide the end-to-end sensitivities of the model weights to the velocity fields. In-sample and out-of-sample testing on multiple dual-jet configurations show that its required mesh density in each coordinate direction for prediction of mean flow, Reynolds stresses, and spectra is half that needed by the dynamic Smagorinsky model for comparable accuracy. The same neural-network model trained directly to match filtered subgrid-scale stresses, without the constraint of being embedded within the flow equations during the training, fails to provide a qualitatively correct prediction. The coupled formulation is generalized to train based only on mean-flow and Reynolds stresses, which are more readily available in experiments. The mean-flow training provides a robust model, which is important, though a somewhat less accurate prediction for the same coarse meshes, as might be anticipated due to the reduced information available for training in this case. The anticipated advantage of the formulation is that the inclusion of resolved physics in the training increases its capacity to extrapolate. This is assessed for the case of passive scalar transport, for which it outperforms established models due to improved mixing predictions.

42 ENGINEERING↗

Low Dissipative High Order Shock-Capturing Methods using Characteristic-Based Filters

An approach which closely maintains the non-dissipative nature of classical fourth or higher- order spatial differencing away from shock waves and steep gradient regions while being capable of accurately capturing discontinuities, steep gradient and fine scale turbulent structures in a stable and efficient manner is described. The approach is a generalization of the method of Gustafsson and Olsson and the artificial compression method (ACM) of Harten. Spatially non-dissipative fourth or higher-order compact and non-compact spatial differencings are used as the base schemes. Instead of applying a scalar filter as in Gustafsson and Olsson, an ACM like term is used to signal the appropriate amount of second or third-order TVD or ENO types of characteristic based numerical dissipation. This term acts as a characteristic filter to minimize numerical dissipation for the overall scheme. For time-accurate computations, time discretizations with low dissipation are used. Numerical experiments on 2-D vortical flows, vortex-shock interactions and compressible spatially and temporally evolving mixing layers showed that the proposed schemes have the desired property with only a 10% increase in operations count over standard second-order TVD schemes. Aside from the ability to accurately capture shock-turbulence interaction flows, this approach is also capable of accurately preserving vortex convection. Higher accuracy is achieved with fewer grid points when compared to that of standard second-order TVD or ENO schemes. To demonstrate the applicability of these schemes in sustaining turbulence where shock waves are absent, a simulation of 3-D compressible turbulent channel flow in a small domain is conducted.

Yee, H. C.↗

Low Dissipative High Order Shock-Capturing Methods Using Characteristic-Based Filters

An approach which closely maintains the non-dissipative nature of classical fourth or higher- order spatial differencing away from shock waves and steep gradient regions while being capable of accurately capturing discontinuities, steep gradient and fine scale turbulent structures in a stable and efficient manner is described. The approach is a generalization of the method of Gustafsson and Oisson and the artificial compression method (ACM) of Harten. Spatially non-dissipative fourth or higher-order compact and non-compact spatial differencings are used as the base schemes. Instead of applying a scalar filter as in Gustafsson and Olsson, an ACM like term is used to signal the appropriate amount of second or third-order TVD or ENO types of characteristic based numerical dissipation. This term acts as a characteristic filter to minimize numerical dissipation for the overall scheme. For time-accurate computations, time discretizations with low dissipation are used. Numerical experiments on 2-D vortical flows, vortex-shock interactions and compressible spatially and temporally evolving mixing layers showed that the proposed schemes have the desired property with only a 10% increase in operations count over standard second-order TVD schemes. Aside from the ability to accurately capture shock-turbulence interaction flows, this approach is also capable of accurately preserving vortex convection. Higher accuracy is achieved with fewer grid points when compared to that of standard second-order TVD or ENO schemes. To demonstrate the applicability of these schemes in sustaining turbulence where shock waves are absent, a simulation of 3-D compressible turbulent channel flow in a small domain is conducted.

Yee, H. C.↗

A regularized-interface method as a unified formulation for simulations of high-pressure multiphase flows

The injection of multi-species fluids into high-pressure and high-temperature environments beyond the species' critical points is commonly found in engineering applications. At these conditions, for immiscible species, both subcritical interfacial dynamics and supercritical mixing can coexist due to variations in temperature around the mixture critical point. The modeling of these complex transcritical phenomena for large-scale configurations is so far not possible. To address this issue, we propose the Regularized-Interface Method (RIM) as a unified formulation that can describe both sub- and supercritical processes as well as the transition between them. The proposed method is derived via filtering of the nanoscale interface-resolving formulation based on van der Waals' linear gradient theory. Thus, this approach allows for the consistent modeling of interfacial dynamics that vanishes at supercritical conditions, while significantly reducing the temporal and spatial resolution constraints of the original nanoscale formulation. The resulting RIM formulation is examined in interface-capturing simulations of sub-, trans-, and supercritical fuel injection processes, involving droplets and jets. Furthermore, these results highlight the importance of resolving spatio-temporal transitions from subcritical interfacial dynamics to supercritical mixing in high-pressure multiphase simulations, in contrast to commonly employed diffused-interface methods, where interfacial dynamics are often neglected.

Interface capturing↗

Accelerating Noisy VQE Optimization with Gaussian Processes

Hybrid variational quantum algorithms, which combine a classical optimizer with evaluations on a quantum chip, are the most promising candidates to show quantum advantage on current noisy, intermediate-scale quantum (NISQ) devices. The classical optimizer is required to perform well in the presence of noise in the objective function evaluations, or else it becomes the weakest link in the algorithm. We introduce the use of Gaussian Processes (GP) as surrogate models to reduce the impact of noise and to provide high quality seeds to escape local minima, whether real or noise-induced. We build this as a framework on top of local optimizations, for which we choose Implicit Filtering (ImFil) in this study. ImFil is a state-of-the-art, gradient-free method, which in comparative studies has been shown to outperform on noisy VQE problems. The result is a new method: "GP+ImFil". We show that when noise is present, the GP+ImFil approach finds results closer to the true global minimum in fewer evaluations than standalone ImFil, and that it works particularly well for larger dimensional problems. Using GP to seed local searches in a multi-modal landscape shows mixed results: although it is capable of improving on ImFil standalone, it does not do so consistently and would only be preferred over other, more exhaustive, multistart methods if resources are constrained.

Muller, Juliane↗