Search NASA⌕ Search

SEARCH · Search NASA

Results for “Residual point distribution”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A comprehensive study of non-adaptive and residual-based adaptive sampling for physics-informed neural networks

Physics-informed neural networks (PINNs) have shown to be effective tools for solving both forward and inverse problems of partial differential equations (PDEs). PINNs embed the PDEs into the loss of the neural network using automatic differentiation, and this PDE loss is evaluated at a set of scattered spatio-temporal points (called residual points). The location and distribution of these residual points are highly important to the performance of PINNs. However, in the existing studies on PINNs, only a few simple residual point sampling methods have mainly been used. Here, we present a comprehensive study of two categories of sampling for PINNs: non-adaptive uniform sampling and adaptive nonuniform sampling. We consider six uniform sampling methods, including (1) equispaced uniform grid, (2) uniformly random sampling, (3) Latin hypercube sampling, (4) Halton sequence, (5) Hammersley sequence, and (6) Sobol sequence. We also consider a resampling strategy for uniform sampling. To improve the sampling efficiency and the accuracy of PINNs, we propose two new residual-based adaptive sampling methods: residual-based adaptive distribution (RAD) and residual-based adaptive refinement with distribution (RAR-D), which dynamically improve the distribution of residual points based on the PDE residuals during training. Hence, we have considered a total of 10 different sampling methods, including six non-adaptive uniform sampling, uniform sampling with resampling, two proposed adaptive sampling, and an existing adaptive sampling. We extensively tested the performance of these sampling methods for four forward problems and two inverse problems in many setups. Our numerical results presented in this study are summarized from more than 6000 simulations of PINNs. Here, we show that the proposed adaptive sampling methods of RAD and RAR-D significantly improve the accuracy of PINNs with fewer residual points for both forward and inverse problems. Furthermore, the results obtained in this study can also be used as a practical guideline in choosing sampling methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Gradient-enhanced physics-informed neural networks for forward and inverse PDE problems

Deep learning has been shown to be an effective tool in solving partial differential equations (PDEs) through physics-informed neural networks (PINNs). PINNs embed the PDE residual into the loss function of the neural network, and have been successfully employed to solve diverse forward and inverse PDE problems. However, one disadvantage of the first generation of PINNs is that they usually have limited accuracy even with many training points. Here, we propose a new method, gradient-enhanced physics-informed neural networks (gPINNs), for improving the accuracy of PINNs. gPINNs leverage gradient information of the PDE residual and embed the gradient into the loss function. We tested gPINNs extensively and demonstrated the effectiveness of gPINNs in both forward and inverse PDE problems. Our numerical results show that gPINN performs better than PINN with fewer training points. Additionally, we combined gPINN with the method of residual-based adaptive refinement (RAR), a method for improving the distribution of training points adaptively during training, to further improve the performance of gPINN, especially in PDEs with solutions that have steep gradients.

42 ENGINEERING↗
NASA NTRS2007

Accuracy Analysis on Large Blocks of High Resolution Images

Although high altitude frequencies effects are removed at the time of basic image generation, low altitude (Yaw) effects are still present in form of affinity/angular affinity. They are effectively removed by additional parameters. Bundle block adjustment based on properly weighted ephemeris/altitude quaternions (BBABEQ) are not enough to remove the systematic effect. Moreover, due to the narrow FOV of the HRSI, position and altitude are highly correlated making it almost impossible to separate and remove their systematic effects without extending the geometric model (Self-Calib.) The systematic effects gets evident on the increase of accuracy (in terms of RMSE at GCPs) for looser and relaxed ground control at the expense of large and strong block deformation with large residuals at check points. Systematic errors are most freely distributed and their effects propagated all over the block.

Passini, Richardo M.↗
NASA NTRS1981

The MSS control point location error filter for LANDSAT-D

The filter produces the maximum likelihood estimates for average values of the spacecraft position and attitude errors during a single scene. The quality of the filter performance is characterized by the maximal cross and along track residual errors for which probability distributions can be calculated analytically for a given pattern of control points. The filter with an automatic selection of the best set of estimates provides geodetic correction at 90% of pixels with residual errors less than 40 m for four or more control points and the mean squared measurement errors of the order of 20-25 m. The same accuracy can be preserved for eight or more control points and measurement errors of 30-35 m.

Levine, I.↗
DOE OSTI2022

KaRMMa – kappa reconstruction for mass mapping

ABSTRACT We present KaRMMa, a novel method for performing mass map reconstruction from weak-lensing surveys. We employ a fully Bayesian approach with a physically motivated lognormal prior to sample from the posterior distribution of convergence maps. We test KaRMMa on a suite of dark matter N-body simulations with simulated DES Y1-like shear observations. We show that KaRMMa outperforms the basic Kaiser–Squires mass map reconstruction in two key ways: (1) our best map point estimate has lower residuals compared to Kaiser–Squires; and (2) unlike the Kaiser–Squires reconstruction, the posterior distribution of KaRMMa maps is nearly unbiased in all summary statistics we considered, namely: one-point and two-point functions, and peak/void counts. In particular, KaRMMa successfully captures the non-Gaussian nature of the distribution of κ values in the simulated maps. We further demonstrate that the KaRMMa posteriors correctly characterize the uncertainty in all summary statistics we considered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Spatial variation of the low-temperature micro-photoluminescence of THM-grown CdZnTeSe

The quaternary compound CdZnTeSe (CZTS) has emerged as a next-generation detector material. Unlike CdZnTe (CZT), CZTS has the distinct advantages of little or no sub-grain boundary networks and a much lower concentration of Te inclusions, plus better compositional homogeneity. Thus, the material is expected to offer better spatial charge-transport homogeneity compared to CZT. However, an inhomogeneous distribution of point defects in both CZT and CZTS, such as Cd vacancies and residual impurities, can impose additional spatial inhomogeneity of the resistivity, carrier mobilities, and carrier lifetimes of the material. To investigate the spatial distribution of such defects in CZTS, low-temperature photoluminescence (PL) spectroscopic studies were performed at different positions along a single crystalline CZTS sample surface, which was grown by the traveling heater method (THM). In conclusion, the intensity variation of the PL emission of excitons bound to a neutral acceptor defect (A 0 , X), which is likely dominated by an acceptor like Cu-related level, with respect to the PL emission line from a neutral donor-bound exciton (D 0 , X) exhibited significant spatial variation, while the peak energy positions were approximately uniform due to the high compositional and bandgap homogeneity of the material.

II-VI semiconductor↗
DOE OSTI2024

Uncertainty Quantification in CO2 Trapping Mechanisms: A Case Study of PUNQ-S3 Reservoir Model Using Representative Geological Realizations and Unsupervised Machine Learning

Evaluating uncertainty in CO2 injection projections often requires numerous high-resolution geological realizations (GRs) which, although effective, are computationally demanding. This study proposes the use of representative geological realizations (RGRs) as an efficient approach to capture the uncertainty range of the full set while reducing computational costs. A predetermined number of RGRs is selected using an integrated unsupervised machine learning (UML) framework, which includes Euclidean distance measurement, multidimensional scaling (MDS), and a deterministic K-means (DK-means) clustering algorithm. In the context of the intricate 3D aquifer CO2 storage model, PUNQ-S3, these algorithms are utilized. The UML methodology selects five RGRs from a pool of 25 possibilities (20% of the total), taking into account the reservoir quality index (RQI) as a static parameter of the reservoir. To determine the credibility of these RGRs, their simulation results are scrutinized through the application of the Kolmogorov–Smirnov (KS) test, which analyzes the distribution of the output. In this assessment, 40 CO2 injection wells cover the entire reservoir alongside the full set. The end-point simulation results indicate that the CO2 structural, residual, and solubility trapping within the RGRs and full set follow the same distribution. Simulating five RGRs alongside the full set of 25 GRs over 200 years, involving 10 years of CO2 injection, reveals consistently similar trapping distribution patterns, with an average value of Dmax of 0.21 remaining lower than Dcritical (0.66). Using this methodology, computational expenses related to scenario testing and development planning for CO2 storage reservoirs in the presence of geological uncertainties can be substantially reduced.

Mahjour, Seyed Kourosh↗
NASA NTRS1984

Comparative assessment of LANDSAT-D MSS and TM data quality for mapping applications in the Southeast

Rectifications of multispectral scanner and thematic mapper data sets for full and subscene areas, analyses of planimetric errors, assessments of the number and distribution of ground control points required to minimize errors, and factors contributing to error residual are examined. Other investigations include the generation of three dimensional terrain models and the effects of spatial resolution on digital classification accuracies.

Source record↗
NASA NTRS1985

The X-ray halo of Cygnus X-1

Four Einstein HRI images of Cygnus X-1 were examined for the presence of a halo due to scattering of X-rays by interstellar grains. The analysis technique exploits the intrinsic aperiodic variability of the source to map the point response function of the optics. A residual, nonvariable component to the surface brightness distribution (comprising approximately more than 12 percent of the source flux) is interpreted as a scattered halo. The halo flux does not reflect the short term time variability of the central source as it is smoothed by differential time delays of order days. The Cygnus X-1 halo is consistent with those of other sources derived in previous studies using different techniques. Comparison is made with a scattering model, and the sensitivity of the halo flux to maximal grain size is demonstrated.

Bode, M. F.↗
NASA NTRS2008

On Combining Thermal-Infrared and Radio-Occultation Data of Saturn's Atmosphere

Radio-occultation and thermal-infrared measurements are complementary investigations for sounding planetary atmospheres. The vertical resolution afforded by radio occultations is typically approximately 1 km or better, whereas that from infrared sounding is often comparable to a scale height. On the other hand, an instrument like CIRS can easily generate global maps of temperature and composition, whereas occultation soundings are usually distributed more sparsely. The starting point for radio-occultation inversions is determining the residual Doppler-shifted frequency, that is the shift in frequency from what it would be in the absence of the atmosphere. Hence the positions and relative velocities of the spacecraft, target atmosphere, and DSN receiving station must be known to high accuracy. It is not surprising that the inversions can be susceptible to sources of systematic errors. Stratospheric temperature profiles on Titan retrieved from Cassini radio occultations were found to be very susceptible to errors in the reconstructed spacecraft velocities (approximately equal to 1 mm/s). Here the ability to adjust the spacecraft ephemeris so that the profiles matched those retrieved from CIRS limb sounding proved to be critical in mitigating this error. A similar procedure can be used for Saturn, although the sensitivity of its retrieved profiles to this type of error seems to be smaller. One issue that has appeared in inverting the Cassini occultations by Saturn is the uncertainty in its equatorial bulge, that is, the shape in its iso-density surfaces at low latitudes. Typically one approximates that surface as a geopotential surface by assuming a barotropic atmosphere. However, the recent controversy in the equatorial winds, i.e., whether they changed between the Voyager (1981) era and later (after 1996) epochs of Cassini and some Hubble observations, has made it difficult to know the exact shape of the surface, and it leads to uncertainties in the retrieved temperature profiles of one to a few kelvins. This propagates into errors in the retrieved helium abundance, which makes use of thermal-infrared spectra and synthetic spectra computed with retrieved radio-occultation temperature profiles. The highest abundances are retrieved with the faster Voyager-era winds, but even these abundances are somewhat smaller than those retrieved from the thermal-infrared data alone (albeit with larger formal errors). The helium abundance determination is most sensitive to temperatures in the upper troposphere. Further progress may include matching the radio-occultation profiles with those from CIRS limb sounding in the upper stratosphere.

Flasar, F. M.↗
DOE OSTI2024

Point-of-use filtration units as drinking water distribution system sentinels

Abstract Municipal drinking water distribution systems (DWDSs) and associated premise plumbing (PP) systems are vulnerable to proliferation of opportunistic pathogens, even when chemical disinfection residuals are present, thus presenting a public health risk. Monitoring the structure of microbial communities of drinking water is challenging because of limited continuous access to faucets, pipes, and storage tanks. We propose a scalable household sampling method, which uses spent activated carbon and reverse osmosis (RO) membrane point-of-use (POU) filters to evaluate mid- to long-term occurrence of microorganisms in PP systems that are relevant to consumer exposure. As a proof of concept, POU filter microbiomes were collected from four different locations and analyzed with 16S rRNA gene amplicon sequencing. The analyses revealed distinct microbial communities, with occasional detection of potential pathogens. The findings highlight the importance of local, and if possible, continuous monitoring within and across distribution systems. The continuous operation of POU filters offers an advantage in capturing species that may be missed by instantaneous sampling methods. We suggest that water utilities, public institutions, and regulatory agencies take advantage of end-of-life POU filters for microbial monitoring. This approach can be easily implemented to ensure drinking water safety, especially from microbes of emerging concerns; e.g., pathogenic Legionella and Mycobacterium species.

42 ENGINEERING↗
DOE OSTI2025

Risk Assessment in a Chemical Laboratory Following an Explosive Incident Involving a Novel Diazonium Compound: Retrospective Analysis and Lessons Learned

Diazonium compounds are synthetically useful in the production of dyes and textiles, however they are highly explosive under dry conditions. Explosion prevention becomes more difficult when new diazonium compounds are synthesized, because while some syntheses include a counterion to increase their stability, this is not always a reliable method to prevent an explosive incident. Due to the uncertainty surrounding the explosiveness of different diazonium compounds, it is important to understand how to safely clean up after an incident and how to determine when it is safe to return a laboratory to typical operational use, particularly when the incident involves a novel compound where a standard does not exist for instrument calibration. Here, an explosive event is discussed involving the synthesis of 4-bromo-benzenediazonium-2-carboxylate. Following the explosive incident and 3-step cleanup, which involved a precautionary neutralization step, samples were collected from the fume hood where the incident occurred. Because the incident involved an unstable, novel compound that is not commercially available and was deemed unsafe to resynthesize for instrument calibration, we assessed the risk of further explosion by analyzing for the stable decomposition products. Mass spectrometry analysis confirmed that the residue in the fume hood contained 5-bromosalicylic acid, a decomposition product of 4-bromo-benzenediazonium-2-carboxylate. Samples were taken from multiple points in the fume hood and analyzed to estimate the spatial distribution of the decomposition product. Based on this analysis, we inferred that the primary decomposition product was far more abundant than residual energetic, indicating the energetic had been consumed or neutralized to a trace quantity where the risk of further explosion was low. Furthermore, the steps presented here─specifically, initial neutralization and then analyzing the spatial distribution of expected decomposition products to assess risk when a novel explosive material is detonated in a confined space─were our approach to assess further risk following an explosion due to a novel diazonium compound without the need for any further handling or resynthesis of the energetic. Here, we present our approach and critically analyze these steps by discussing retrospective lessons learned and alternative analytical approaches.

Computer simulations↗
NASA NTRS1978

Ion tracking in an electrostatic potential distribution

An ion tracking system utilizing a single detector was designed to resolve two dimensional potential distributions in a vacuum. The ion detector consisted of a continuous dynode electron multiplier tube mounted behind a series of aperture plates. The detector was located on a track below a wedge shaped structure used to produce the desired potential distribution. A 150 volt electron beam aimed at selected regions, ionized the residual gas in the chamber. The ions were ejected from their source points by electrostatic forces causing them to drift downward toward the detector. Measurements of particle energy, detector position, and entry angle of the incoming particle flux were used to initialize a computer trajectory-tracing program. The program assumes a known theoretically calculated potential distribution near the wedge for computation of particle trajectories, providing a basis for evaluating the data. Special attention was given to characterization of the ion detector.

Ross, D. P.↗
NASA NTRS1982

Residual stress in high modulus carbon fibers

The modulus and residual strain in carbon fibers are measured by successively electrochemically milling away the fiber surface. Electrochemical etching is found to remove the carbon fiber surface very uniformly, in contrast to air and wet oxidation. The precision of fiber diameter measurements is improved by using a laser diffraction technique instead of optical methods. More precise diameter measurements reveal that past correlations of diameter and fiber modulus are largely measurement artifact. The moduli of most carbon fibers decrease after the outer layers of the fibers are removed. Owing to experimental difficulties, the moduli and strengths of the fibers at their centers are not determined, and moduli are estimated on the basis of microstructure. The calculated residual stresses are found to be insensitive to these moduli estimates as well as the exact form of regression equation used to describe the moduli and residual strain distributions. Axial compressive residual stresses are found to be very high for some higher modulus carbon fibers. It is pointed out that the compressive stress makes the fibers insensitive to surface flaws when loaded in tension but it may initiate failure by buckling when loaded in compression.

Chen, K. J.↗
NASA NTRS1974

Lunar satellite techniques applicable to earth satellite geodesy

The techniques considered are related to gravity profiling from Doppler residuals, a solution for a surface mass distribution, and a dynamical approach for the estimation of random forces at each data point. The dynamical approach makes use of a new filtering method which uses sequential estimation theory to determine unknown and somewhat correlated forces which randomly perturb the orbit. It is concluded that lunar data reduction methods would be suitable for applications in connection with the GEOS-C/ATS-F mission.

Sjogren, W. L.↗
NASA NTRS2024

Deciphering the Spectra of Flowers to Map Landscape-scale Blooming Dynamics

Like leaves, floral coloration is driven by inherent optical properties, which are determined by pigments, scattering structure, and thickness. However, establishing the relative contribution of these factors to canopy spectral signals is usually limited to in-situ observations. Modeling flowering dynamics (e.g., blooming duration, spatial distribution) at the landscape scale may reveal insights into ecological processes and phenological adaptations to environmental changes. Multitemporal visible to shortwave infrared (VSWIR) imaging spectroscopy observations are especially suited for such efforts. Reflectance in this spectral range is sensitive to major flower pigments, flowering phenology traces, and biophysical differences between flowers and other plant parts. We explored how flowers contribute to spectral signals using a time series of imagery from the Airborne Visible InfraRed Imaging Spectrometer - Next Generation (AVIRIS-NG) collected as part of the SBG High-Frequency Time Series (SHIFT) campaign as a case study. Airborne data were collected weekly during the spring of 2022 across two natural reserves in California. Field spectra were gathered from blooming plots at leaf, flower, and canopy levels at two time points during the campaign. The processed data was used to investigate flowering species' spectro-temporal variation and spatial distribution using Spectral Mixture Residual, Gaussian clustering techniques, and a proposed narrow-band flowering index. Linear spectral unmixing allowed the computation of the weighted contribution of four major high-variance endmembers (leaves, flowers, soil, dark) and low-variance residual signal that comprises subtle spectral features used to track biophysical processes. The reflectance residual was projected on a low principal component basis to characterize flowering clusters' variation and spatial distribution based on the Gaussian mixture model, providing an uncertainty metric to assess the results. Mapping flowering events from modeling spectro-temporal dynamics throughout the season, from pre-blooming to post-flowering stages, allowed us to identify gradient variations in spectral features within the VSWIR spectral range linked to flowering pigments. Time series of the Mixture Residual Blooming Index and the Red-Edge Normalized Difference Vegetation Index revealed specific flowering and greenness phenophases across the two main species (Coreopsis gigantea, Artemisia californica) in the flowering areas. Overall, our approach opens opportunities for future satellite monitoring of floral cycles at broader scales.

Yoseline Angel↗
NASA NTRS1983

Stress studies in EFG

Stress distributions were calculated for a creep law to predict a rate of plastic deformation. The expected reduction in stresses is obtained. Improved schemes for calculating growth system temperature distributions were evaluated. Temperature field modeling examined the possibility of using horizontal temperature gradients to influence stress distribution in ribbon. The defect structure of 10 cm wide ribbon grown in the cartridge system was examined. A new feature is identified from an examination of cross sectional micrographs. It consists of high density dislocation bands extending through the ribbon thickness. A four point bending apparatus was constructed for high temperature study of the creep response of silicon, to be used to generate defects for comparison with as grown defects in ribbon. The feasibility of laser interferometric techniques for sheet residual stress distribution measurement is examined. The mathematical formalism for calculating residual stress from changes in surface topology caused by an applied stress in a rectangular specimen was developed, and the system for laser interferometric measurement to obtain surface topology data was tested on CZ silicon.

Source record↗
NASA NTRS1991

Diabatic heating rate estimates from European Centre for Medium-Range Weather Forecasts analyses

Vertically integrated diabatic heating rate estimates (H) calculated from 32 months of European Center for Medium-Range Weather Forecasts daily analyses (May 1985-December 1987) are determined as residuals of the thermodynamic equation in pressure coordinates. Values for global, hemispheric, zonal, and grid point H are given as they vary over the time period examined. The distribution of H is compared with previous results and with outgoing longwave radiation (OLR) measurements. The most significant negative correlations between H and OLR occur for (1) tropical and Northern-Hemisphere mid-latitude oceanic areas and (2) zonal and hemispheric mean values for periods less than 90 days. Largest positive correlations are seen in periods greater than 90 days for the Northern Hemispheric mean and continental areas of North Africa, North America, northern Asia, and Antarctica. The physical basis for these relationships is discussed. An interyear comparison between 1986 and 1987 reveals the ENSO signal.

Christy, John R.↗