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 145 records · Page 8

Determination of Total Magnetic Anomalies and Their Vertical Gradients of Swarm-A Satellite over Central Europe and Pannonian Basin

Our paper discusses the determination of total magnetic field anomalies derived from the Swarm–A satellite data; one of two low orbiting satellites of the three Swarm formations. This procedure requires several modifications. The first step is the conversion of the measured CDF data to the ASCII format. This step is followed by the selection of the data with K(sub p) index ≤ 1(sub +). The anomalies are determined by the removal of the IGRF from the resulting satellite data. There are two Swarm–A data sets: descending (6000) orbits and ascending (5688) orbits. For our study the descending orbits were used. The next step of calculations is to determine the difference of the two-dimensional linear field fitted to the Swarm–A anomalies and the anomalies given. These anomalies are filtered by Gaussian low-pass filter. The last step of the corrections is the removal of the direct component, zero spatial frequency, from the descending anomalies. The anomalies and their vertical gradients are qualitatively interpreted over Central Europe and the Pannonian Basin.

Kis, K.↗

Real-time breath analysis towards a healthy human breath profile

Abstract The direct analysis of molecules contained within human breath has had significant implications for clinical and diagnostic applications in recent decades. However, attempts to compare one study to another or to reproduce previous work are hampered by: variability between sampling methodologies, human phenotypic variability, complex interactions between compounds within breath, and confounding signals from comorbidities. Towards this end, we have endeavored to create an averaged healthy human ‘profile’ against which follow-on studies might be compared. Through the use of direct secondary electrospray ionization combined with a high-resolution mass spectrometry and in-house bioinformatics pipeline, we seek to curate an average healthy human profile for breath and use this model to distinguish differences inter- and intra-day for human volunteers. Breath samples were significantly different in PERMANOVA analysis and ANOSIM analysis based on Time of Day, Participant ID, Date of Sample, Sex of Participant, and Age of Participant ( p < 0.001). Optimal binning analysis identify strong associations between specific features and variables. These include 227 breath features identified as unique identifiers for 28 of the 31 participants. Four signals were identified to be strongly associated with female participants and one with male participants. A total of 37 signals were identified to be strongly associated with the time-of-day samples were taken. Threshold indicator taxa analysis indicated a shift in significant breath features across the age gradient of participants with peak disruption of breath metabolites occurring at around age 32. Forty-eight features were identified after filtering from which a healthy human breath profile for all participants was created.

60 APPLIED LIFE SCIENCES↗

An updated estimate of the Mu2e experiment sensitivity

The Mu2e experiment at Fermilab will search for the conversion of a negative muon into an electron inside the field of a nucleus. This process does not conserve charged-lepton flavour and is heavily suppressed in the Standard Model (SM), with a branching ratio < 10-50. Any evidence of it would be an undeniable evidence of new physics beyond the SM. The project sets out to achieve a single event sensitivity of $\sim 3 × 10^{-17}$ on the ratio between the probability for a conversion of a negative muon into an electron and the one for a muon capture by the nucleus. Such a sensitivity would represent a 4 orders of magnitude improvement on the previous upper limit for the process, making possible to test predictions of different extensions of the SM. Mu2e uses three superconducting solenoids to produce and measure the muon conversions. In the first solenoid, the Production Solenoid, pions and kaons are produced, together with other secondary products, by 8 GeV kinetic energy proton interactions in a tungsten target. A gradient magnetic field is specifically designed to direct low momentum particles into the Transport Solenoid, an S-shaped magnet that filters out particles with unwanted charge and momentum. Muons, produced by pion and kaon decays, finally reach the aluminum Stopping Target in the Detector Solenoid, where they eventually stop and convert. The result of the conversion process is a monochromatic electron of ~105 MeV/c momentum. The Detector Solenoid also hosts the two main detectors: a straw tube tracker and two CsI calorimeter disks, both providing measurement of event kinematics. A germanium detector and a LaBr crystal are located downstream of the Detector Solenoid to measure the X and gamma rays produced by the muon captures in the Stopping Target. A veto system of scintillators covering the Detector Solenoid and half of the Transport Solenoid is used to identify and reject cosmic rays interactions. With respect to the initial project, the Mu2e running plan has evolved to a staged configuration with 2 years at reduced intensity before the 2025 accelerator shutdown for the neutrino beam upgrade and 2 or 3 years at full intensity after that. This, together with geometry changes and a better knowledge of detector performances obtained by the first slice tests, has required a full revision of the signal over background selection that is the subject of this thesis. In order to achieve a new estimate of Mu2e sensitivity, the simulation of the data corresponding to the first 2 years of data acquisition has been performed. This includes both conversion electrons (CE) and the main sources of background: cosmics, decay-in orbits (DIO), radiative pion captures (RPC) and antiprotons. The characteristics of the signal and of the main backgrounds have been studied to define the best selection variables for CE. A special effort has been devoted to the evaluation of the antiproton background. The lack of experimental data for antiproton production cross section makes the systematic uncertainty on this background significant. A new parameterization of the cross section has been developed to fit the existing data and to provide a more reliable estimate of the systematic uncertainty by comparing the results of the old and the new model. A special effort has been devoted to the optimization of the antiprotons Monte Carlo generator. This study has also revealed that the dominant component within this background is represented by antiprotons produced in the opposite direction with respect to the Transport Solenoid entrance and then redirected to it by back-scattering processes in the Production Target. This ultimately highlights the sensitivity of the background estimate to the G4 handling of antiproton interactions in the 1 to 3 GeV/c momentum range. Finally, the final experimental sensitivity has been studied. The momentum and time selection have been optimized to obtain the 5-sigma discovery reach or, in case of no signal, the upper limit on conversion probability. The results confirm that, in the firs t two years of data taking, Mu2e will be able to improve the current experimental sensitivity for muon-to-electron conversion in an atomic field by more than 3 order of magnitudes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Description and discussion of the NRL TGDCC

The NRL thermal gradient diffusion cloud chamber (TGDCC) consists of two plates 7.5 cm in diameter separated by 1.25 cm and covered with saturated filter paper. The cylindrical wall separating the plates is glass. The top plate is at room temperature and the bottom plate is cooled with a thermoelectric cooler. The temperature difference is measured with several sets of thermocouples. The CCN concentration was determined from the video recording. This procedure of averaging the maximum count obtained on several successive recordings at the same supersaturation results in concentrations which are somewhat higher than concentrations calculated from an average across the plateau.

Hoppel, W. A.↗

Toward a gradiometer analytic model

A new model is developed to model the data type, formulate the filter structure, and perform covariance for an orbiting gravity gradiometer, which is much more elaborate and realistic than the earlier gravity gradient model of Sonnabend and McEneaney (1988). The main new features of the new model are a general inertia tensor for the floated instrument, air drag and radiation pressure models using cubic power spectra, and more reasonable kinematics.

Sonnabend, Dave↗

Damping of thermal acoustic oscillations in hydrogen systems

Acoustic waves initiated by a large temperature gradient along a tube are defined as thermal acoustic oscillations (TAOs). These oscillations have been damped by introducing such sound absorbing techniques as acoustic filters, resonators, etc.. These devices serve as an acoustic sink that is used to absorb or dissipate the acoustic energy thereby eliminating or damping such oscillations. Several empirical damping techniques, such as attaching a resonator as a side branch or inserting a wire in the tube, have been developed in the past and have provided reasonable success. However, the effect of connecting tube radius, length, and resonator volume on the damping of thermal acoustic oscillations has not been evaluated quantitatively. Further, these methods have not been effective when the oscillating tube radius was relatively large. Detailed theoretical analyses of these techniques including a newly developed method for damping oscillations in a tube of relatively large radius are provided in this presentation.

Gu, Youfan↗

Thin Semiconductor/Metal Films For Infrared Devices

Spectral responses of absorbers and reflectors tailored. Thin cermet films composites of metals and semiconductors undergoing development for use as broadband infrared reflectors and absorbers. Development extends concepts of semiconductor and dielectric films used as interference filters for infrared light and visible light. Composite films offer advantages over semiconductor films. Addition of metal particles contributes additional thermal conductivity, reducing thermal gradients and associated thermal stresses, with resultant enhancements of thermal stability. Because values of n in composite films made large, same optical effects achieved with lesser thicknesses. By decreasing thicknesses of films, one not only decreases weights but also contributes further to reductions of thermal stresses.

Lamb, James L.↗

Precipitation-Based ENSO Indices

In this study gridded observed precipitation data sets are used to construct rainfall-based ENSO indices. The monthly El Nino and La Nina Indices (EI and LI) measure the steepest zonal gradient of precipitation anomalies between the equatorial Pacific and the Maritime Continent. This is accomplished by spatially averaging precipitation anomalies using a spatial boxcar filter, finding the maximum and minimum averages within a Pacific and Maritime Continent domain for each month, and taking differences. EI and LI can be examined separately or combined to produce one ENSO Precipitation Index (ESPI). ESPI is well correlated with traditional sea surface temperature and pressure indices, leading Nino 3.4. One advantage precipitation indices have over more conventional indices, is describing the strength and position of the Walker circulation. Examples are given of tracking the impact of ENSO events on the tropical precipitation fields.

Adler, Robert↗

Contributions of Spherical Harmonics to Gravitational Moment

A scalar gravitational potential function expressed as a series of spherical harmonics frequently serves as the basis for a model of an astronomical body's gravitational field. The contribution of a generic spherical harmonic to gravitational gradient is expressed as a dyadic, which is then used to obtain an analytical expression in vector-dyadic form for the contribution to the moment of gravitational forces about the mass center of a small body such as a spacecraft. The expression developed for a harmonic's contribution to gravitational gradient can be applied in areas beyond the scope of the paper; for example, gravitational gradient plays an important role in the state propagation matrix and the state transition matrix that are used in spacecraft trajectory targeting and Kalman filtering. Additionally, it can be employed in numerical simulations of orbit determination based on measurements obtained with a gradiometer in low-Earth orbit. Contributions of spherical harmonics to gravitational moment may be of interest in connection with attitude control of a spacecraft in the vicinity of a body with an irregular shape, such as an asteroid. Normalized spherical harmonic coefficients up to degree and order 10 are obtained for the asteroid 216 Kleopatra and used in numerical evaluations of contributions to gravitational moment.

Roithmayr, Carlos M.↗

Reinforcement-based Program Induction in a Neural Virtual Machine

We present a neural virtual machine that can be trained to perform algorithmic tasks. Rather than combining a neural controller with non-neural memory storage as has been done in the past, this architecture is purely neural and emulates tape-based memory via fast associative weights (onestep learning). Here we formally define the architecture, and then extend the system to learn programs using recurrent policy gradient reinforcement learning based on examples of program inputs labeled with corresponding output targets, which are compared against actual output to generate a sparse reward signal. We describe the policy gradient training procedure used, and report its empirical performance on a number of smallscale list processing tasks, such as finding the maximum list element, filtering out certain elements, and reversing the order of the elements. These results show that program induction via reinforcement learning is possible using sparse rewards and solely neural computations.

Katz, Garrett E.↗

Multi-channel, multi-template event reconstruction for SuperCDMS data using machine learning

SuperCDMS SNOLAB uses kilogram-scale germanium and silicon detectors to search for dark matter. Each detector has Transition Edge Sensors (TESs) patterned on the top and bottom faces of a large crystal substrate, with the TESs electrically grouped into six phonon readout channels per face. Noise correlations are expected among a detector's readout channels, in part because the channels and their readout electronics are located in close proximity to one another. Moreover, owing to the large size of the detectors, energy deposits can produce vastly different phonon propagation patterns depending on their location in the substrate, resulting in a strong position dependence in the readout-channel pulse shapes. Both of these effects can degrade the energy resolution and consequently diminish the dark matter search sensitivity of the experiment if not accounted for properly. We present a new algorithm for pulse reconstruction, mathematically formulated to take into account correlated noise and pulse shape variations. This new algorithm fits N readout channels with a superposition of M pulse templates simultaneously - hence termed the N$\times$M filter. We describe a method to derive the pulse templates using principal component analysis (PCA) and to extract energy and position information using a gradient boosted decision tree (GBDT). We show that these new N$\times$M and GBDT analysis tools can reduce the impact from correlated noise sources while improving the reconstructed energy resolution for simulated mono-energetic events by more than a factor of three and for the 71Ge K-shell electron-capture peak recoils measured in a previous version of SuperCDMS called CDMSlite to $<$ 50 eV from the previously published value of $\sim$100 eV. These results lay the groundwork for position reconstruction in SuperCDMS with the N$\times$M outputs.

Albakry, M. F. [British Columbia U.; TRIUMF]↗

Toward Better Modeling of Supercritical Turbulent Mixing

study was done as part of an effort to develop computational models representing turbulent mixing under thermodynamic supercritical (here, high pressure) conditions. The question was whether the large-eddy simulation (LES) approach, developed previously for atmospheric-pressure compressible-perfect-gas and incompressible flows, can be extended to real-gas non-ideal (including supercritical) fluid mixtures. [In LES, the governing equations are approximated such that the flow field is spatially filtered and subgrid-scale (SGS) phenomena are represented by models.] The study included analyses of results from direct numerical simulation (DNS) of several such mixing layers based on the Navier-Stokes, total-energy, and conservation- of-chemical-species governing equations. Comparison of LES and DNS results revealed the need to augment the atmospheric- pressure LES equations with additional SGS momentum and energy terms. These new terms are the direct result of high-density-gradient-magnitude regions found in the DNS and observed experimentally under fully turbulent flow conditions. A model has been derived for the new term in the momentum equation and was found to perform well at small filter size but to deteriorate with increasing filter size. Several alternative models were derived for the new SGS term in the energy equation that would need further investigations to determine if they are too computationally intensive in LES.

Selle, Laurent↗

A new method for measuring the 3D turbulent velocity dispersion of molecular clouds

ABSTRACT The structure and star formation activity of a molecular cloud are fundamentally linked to its internal turbulence. However, accurately measuring the turbulent velocity dispersion is challenging due to projection effects and observational limitations, such as telescope resolution, particularly for clouds that include non-turbulent motions, such as large-scale rotation. Here, we develop a new method to recover the 3D turbulent velocity dispersion (σv,3D) from position–position–velocity (PPV) data. We simulate a rotating, turbulent, collapsing molecular cloud, and compare its intrinsic σv,3D with three different measures of the velocity dispersion accessible in PPV space: (1) the spatial mean of the 2nd-moment map, σi, (2) the standard deviation of the gradient/rotation-corrected 1st-moment map, σ(c − grad), and (3) a combination of (1) and (2), called the ‘gradient-corrected parent velocity dispersion’, $\sigma _{\mathrm{(p}-\mathrm{grad)}}=(\sigma _{\mathrm{i}}^2+\sigma _{(\mathrm{c}-\mathrm{grad)}}^2)^{1/2}$. We show that the gradient correction is crucial in order to recover purely turbulent motions of the cloud, independent of the orientation of the cloud with respect to the line of sight. We find that with a suitable correction factor and appropriate filters applied to the moment maps, all three statistics can be used to recover σv,3D, with method 3 being the most robust and reliable. We determine the correction factor as a function of the telescope beam size for different levels of cloud rotation, and find that for a beam full width at half-maximum f and cloud radius R, the 3D turbulent velocity dispersion can best be recovered from the gradient-corrected parent velocity dispersion via $\sigma _{v,\mathrm{3D}}= \left[(-0.29\pm 0.26)\, f/R + 1.93 \pm 0.15\right] \sigma _{\mathrm{(p}-\mathrm{grad)}}$ for f/R < 1, independent of the level of cloud rotation or LOS orientation.

Stewart, Madeleine↗

Seismic Spatial Gradients and Machine Learning-Based Classifiers for Explosion Monitoring (LDRD 218327)

This final report summarizes the work completed under the Laboratory Directed Research and Development (LDRD) project “Seismic Spatial Gradients as a Machine Learning-Based Classifier for Explosion Monitoring.” The overarching goal of the project was to explore the efficacy of using machine learning-based classification algorithms where the input data are the spatial gradient of the seismic wavefield collected at a single point on the Earth’s surface. The methods that I describe here are in direct contrast to conventional methods of seismic discrimination which typically rely on a spatially extended network of instruments and physics-based wavefield attributes such as, for example, the ratio between $\textit{P}$ and $\textit{S}$ waves. Rather, we use the spatial gradient of the seismic wavefield observed at a single point on the Earth’s surface and data processing approaches inspired by the machine learning community. We tested two algorithms, a neural network and a modified version of principal component analysis termed Spectrally Filtered Principal Component Analysis (SFPCA). To test these algorithms, we first conducted a series of numerical tests using synthetic data and then conducted a small-scale controlled field experiment. The tests using synthetic data showed that both algorithms had high success rates on gradiometric data, even when simulated noise was added to the signal. Furthermore, we found that using seismic spatial gradients increased the performance of our discrimination algorithms when compared to using just the traditional translational motion seismic data. The tests with field data also showed a high degree of discriminative success.

58 GEOSCIENCES↗

Recursive inverse kinematics for robot arms via Kalman filtering and Bryson-Frazier smoothing

This paper applies linear filtering and smoothing theory to solve recursively the inverse kinematics problem for serial multilink manipulators. This problem is to find a set of joint angles that achieve a prescribed tip position and/or orientation. A widely applicable numerical search solution is presented. The approach finds the minimum of a generalized distance between the desired and the actual manipulator tip position and/or orientation. Both a first-order steepest-descent gradient search and a second-order Newton-Raphson search are developed. The optimal relaxation factor required for the steepest descent method is computed recursively using an outward/inward procedure similar to those used typically for recursive inverse dynamics calculations. The second-order search requires evaluation of a gradient and an approximate Hessian. A Gauss-Markov approach is used to approximate the Hessian matrix in terms of products of first-order derivatives. This matrix is inverted recursively using a two-stage process of inward Kalman filtering followed by outward smoothing. This two-stage process is analogous to that recently developed by the author to solve by means of spatial filtering and smoothing the forward dynamics problem for serial manipulators.

Rodriguez, G.↗

Maximum likelihood tuning of a vehicle motion filter

This paper describes the use of maximum likelihood parameter estimation unknown parameters appearing in a nonlinear vehicle motion filter. The filter uses the kinematic equations of motion of a rigid body in motion over a spherical earth. The nine states of the filter represent vehicle velocity, attitude, and position. The inputs to the filter are three components of translational acceleration and three components of angular rate. Measurements used to update states include air data, altitude, position, and attitude. Expressions are derived for the elements of filter matrices needed to use air data in a body-fixed frame with filter states expressed in a geographic frame. An expression for the likelihood functions of the data is given, along with accurate approximations for the function's gradient and Hessian with respect to unknown parameters. These are used by a numerical quasi-Newton algorithm for maximizing the likelihood function of the data in order to estimate the unknown parameters. The parameter estimation algorithm is useful for processing data from aircraft flight tests or for tuning inertial navigation systems.

Trankle, Thomas L.↗

Airborne lidar research

The planetary boundary layer (PBL) is the lowest layer of the atmosphere in contact with the Earth's surface. It plays an important role in atmospheric circulation and dynamics by influencing surface fluxes of moisture, heat, and momentum. The PBL is characterized by its turbulence structure. Organized convection within the PBL may generate gravity waves within the free troposphere. These waves can propagate vertically, modifying global circulation. Remote sensing of the PBL using a downward-looking lidar is providing new insight into dynamic processes of this important layer and its connection with wave generation. An airborne lidar is an ideal tool to study the PBL. The lidar is installed on the NASA Electra aircraft. It consists of a frequency doubled Nd:YAG laser aligned with a 40 cm diameter telescope. As the laser beam propagates downward from the aircraft toward the surface, it is scattered by the aerosols and molecules in the atmosphere. Generally, there is a sharp gradient of aerosol scattering associated with the PBL top, with high aerosol scattering from within the PBL, and generally very low scattering from the free atmosphere above. It is this gradient in aerosol scattering that is used to visualize the structure of the PBL. A small portion of the scattered radiation is collected by the telescope, filtered, detected by a photomultiplier, and digitized at 100 nsec rate. The horizontal resolution is typically around 13 m, using a 10 Hz laser and a nominal aircraft moves, the resolution of features in the atmosphere below the aircraft is about 15 x 13 m. Results from several studies are presented and discussed.

Melfi, S. Harvey↗

Evaluation of a vortex-based subgrid stress model using DNS databases

The performance of a SubGrid Stress (SGS) model for Large-Eddy Simulation (LES) developed by Misra k Pullin (1996) is studied for forced and decaying isotropic turbulence on a 32(exp 3) grid. The physical viability of the model assumptions are tested using DNS databases. The results from LES of forced turbulence at Taylor Reynolds number R(sub (lambda)) approximately equals 90 are compared with filtered DNS fields. Probability density functions (pdfs) of the subgrid energy transfer, total dissipation, and the stretch of the subgrid vorticity by the resolved velocity-gradient tensor show reasonable agreement with the DNS data. The model is also tested in LES of decaying isotropic turbulence where it correctly predicts the decay rate and energy spectra measured by Comte-Bellot & Corrsin (1971).

Misra, Ashish↗