Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gaussian process”

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 1,315 records · Page 73

Development of Thermal Tomography Imaging System for In-Service Nondestructive Evaluation of Additively Manufactured Components

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of nuclear reactor parts. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process, which can introduce internal material flaws, such as pores and anisotropy. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures due to exposure to high temperature, radiation and corrosive environment in a nuclear reactor. Thermal tomography (TT) provides a capability for non-destructive evaluation of sub-surface defects in arbitrary size structures. We investigate TT of AM stainless steel 316L specimen with imprinted internal porosity defects using relatively low-cost, small form factor infrared (IR) camera based on uncooled micro bolometer detector. Sparse coding related K-means singular value decomposition (SVD) machine learning, image processing algorithms are developed to improve quality of TT images through removal of Additive white Gaussian noise without blurring the images. Following initial qualification of an AM component for deployment in a nuclear reactor, a compact TT can also be used for in-service nondestructive evaluation (NDE). With capability to perform in-service NDE of the AM component lifecycle, TT data can be used for development of a component digital twin. MIT research reactor (MITR) has been identified as a representative environment to perform investigation of feasibility of in-service NDE during shutdown. Preliminary study of radiation level in the MITR viewport, with geometric parameters potentially compatible with compact TT system, has been conducted.

36 MATERIALS SCIENCE↗

Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials (Second Annual Progress Report)

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of nuclear reactor parts. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process, which can introduce internal material flaws, such as pores and anisotropy. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures due to exposure to high temperature, radiation and corrosive environment in a nuclear reactor. Thermal tomography (TT) provides a capability for non-destructive evaluation of sub-surface defects in arbitrary size structures. Thermal tomography is a computational method for heat diffusion-based imaging of solids, which provides 3D visualization of data from flash thermography measurements. We investigate thermal tomography imaging and nondestructive evaluation of stainless steel and nickel super alloy metallic structures produced with laser powder bed fusion (LPBF) additive manufacturing (AM) process. Metallic structures produced with LPBF contain defects, and there are limited capabilities to evaluate these structures non-destructively. Thermal tomography reconstruction of 3D apparent spatial effusivity provides information about AM structure geometry and internal material flaws. We study performance of thermal tomography in imaging of metallic structures through COMSOL computer simulations of transient heat transfer, and through reconstruction of data obtained from experimental measurements. Reconstruction of internal defects is investigated using a stainless steel 316L specimen with flat bottom hole (FBH) indentations, and Inconel 718 plate produced with laser powder bed fusion (LPBF) method, which contains imprinted hemispherical shape low density regions containing non-sintered metallic powder. The FBH’s have the same sizes as the imprinted defects in the LPBF specimens, but offer better imaging contrast. Thermal tomography reconstructions provide visualizations of internal defects, and allow for estimation of their sizes and locations. Detection sensitivity of TT is limited by noises. We investigate separation of signal from noise in thermography images using several machine learning (ML) methods, including new spatio-temporal blind source separation (STBSS) and spatio-temporal sparse dictionary learning (STSDL) methods. Performance of the ML methods is benchmarked using thermography data obtained from imaging stainless steel 316L and Inconel 718 specimens produced LPBF method with imprinted calibrated porosity defects. The ML methods are ranked by F-score and execution runtime. Finally, we investigate TT of AM stainless steel 316L specimen with imprinted internal porosity defects using relatively low-cost, small form factor infrared (IR) camera based on uncooled micro bolometer detector. Sparse coding related K-means singular value decomposition (SVD) machine learning, image processing algorithms are developed to improve quality of TT images through removal of Additive white Gaussian noise without blurring the images. Following initial qualification of an AM component for deployment in a nuclear reactor, a compact TT can also be used for in-service nondestructive evaluation (NDE). With capability to perform in-service NDE of the AM component lifecycle, TT data can be used for development of a component digital twin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Digital simulation of continuous error models with application to an instrument landing system error

A digital simulation of the continuous error of the localized beam of a conventional instrument landing system is discussed. The digital simulation was developed during the analysis of space shuttle navigation capabilities. A discrete mathematical model for use on a digital computer is described. The model generates an output random sequence which is equivalent, for simulation purposes, to the desired random process. The model is a system of difference equations driven by a zero-mean Gaussian random sequence.

Merrick, R. B.↗

A co-designed equalization, modulation, and coding scheme

The commercial impact and technical success of Trellis Coded Modulation seems to illustrate that, if Shannon's capacity is going to be neared, the modulation and coding of an analogue signal ought to be viewed as an integrated process. More recent work has focused on going beyond the gains obtained for Average White Gaussian Noise and has tried to combine the coding/modulation with adaptive equalization. The motive is to gain similar advances on less perfect or idealized channels.

Peile, Robert E.↗

X-ray light curves of active galactic nuclei are phase incoherent

We compute the Fourier phase spectra for the light curves of five low-luminosity active galactic nuclei observed by EXOSAT. There is no statistically significant phase coherence in any of them. This statement is equivalent, subject to a technical caveat, to a demonstration that their fluctuation statistics are Gaussian. Models in which the X-ray output is controlled wholly by a unitary process undergoing a nonlinear limit cycle are therefore ruled out, while models with either a large number of randomly excited independent oscillation modes or nonlinearly interacting spatially dependent oscillations are favored. We also demonstrate how the degree of phase coherence in light curve fluctuations influences the application of causality bounds on internal length scales.

Krolik, Julian↗

Galaxy Correlation Function and Local Density from Photometric Redshifts Using the Stochastic Order Redshift Technique (SORT)

The stochastic order redshift technique (SORT) is a simple, efficient, and robust method to improve cosmological redshift measurements. The method relies upon having a small (∼10 per cent) reference sample of high-quality redshifts. Within pencil- beam-like sub-volumes surrounding each galaxy, we use the precise dN/dz distribution of the reference sample to recover new redshifts and assign them one-to-one to galaxies such that the original rank order of redshifts is preserved. Preserving the rank order is motivated by the fact that random variables drawn from Gaussian probability density functions with different means but equal standard deviations satisfy stochastic ordering. This process is repeated for sub-volumes surrounding each galaxy in the survey. This results in every galaxy being assigned multiple ‘recovered’ redshifts from which a new redshift estimate is determined. An earlier paper applied SORT to a mock Sloan Digital Sky Survey at z 0.2 and accurately recovered the two-point correlation function (2PCF) on scales > 4 h−1Mpc. In this paper, we test the performance of SORT in surveys spanning the redshift range 0.75 < z < 2.25. We used two mock surveys extracted from the Small MultiDark–Planck and Bolshoi–Planck N-body simulations with dark matter haloes that were populated by the Santa Cruz semi-analytic model. We find that SORT overall improves redshift estimates, accurately recovers the redshift-space 2PCF ξ (s) on scales > 2.5 h−1Mpc, and provides improved local density estimates in regions of average or higher density, which may allow for improved understanding of how galaxy properties relate to their environments.

James Kakos↗

A new stationary gridline artifact suppression method based on the 2D discrete wavelet transform

Purpose: In digital x-ray radiography, an antiscatter grid is inserted between the patient and the image receptor to reduce scattered radiation. If the antiscatter grid is used in a stationary way, gridline artifacts will appear in the final image. In most of the gridline removal image processing methods, the useful information with spatial frequencies close to that of the gridline is usually lost or degraded. In this study, a new stationary gridline suppression method is designed to preserve more of the useful information. Methods: The method is as follows. The input image is first recursively decomposed into several smaller subimages using a multiscale 2D discrete wavelet transform. The decomposition process stops when the gridline signal is found to be greater than a threshold in one or several of these subimages using a gridline detection module. An automatic Gaussian band-stop filter is then applied to the detected subimages to remove the gridline signal. Finally, the restored image is achieved using the corresponding 2D inverse discrete wavelet transform. Results: The processed images show that the proposed method can remove the gridline signal efficiently while maintaining the image details. The spectra of a 1D Fourier transform of the processed images demonstrate that, compared with some existing gridline removal methods, the proposed method has better information preservation after the removal of the gridline artifacts. Additionally, the performance speed is relatively high. Conclusions: The experimental results demonstrate the efficiency of the proposed method. Compared with some existing gridline removal methods, the proposed method can preserve more information within an acceptable execution time.

60 APPLIED LIFE SCIENCES↗

Coupled wave-packets for non-adiabatic molecular dynamics: a generalization of Gaussian wave-packet dynamics to multiple potential energy surfaces

Accurate simulation of the non-adiabatic dynamics of molecules in excited electronic states is key to understanding molecular photo-physical processes. Here we present a novel method, based on a semiclassical approximation, that is as efficient as the commonly used mean field Ehrenfest or ad hoc surface hopping methods and properly accounts for interference and decoherence effects. This novel method is an extension of Heller's thawed Gaussian wave-packet dynamics that includes coupling between potential energy surfaces. By studying several standard test problems we demonstrate that the accuracy of the method can be systematically improved while maintaining high efficiency. The method is suitable for investigating the role of quantum coherence in the non-adiabatic dynamics of many-atom molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determination of Cross-Sectional Area of Focused Picosecond Gaussian Laser Beam

Measurement of the waist diameter of a focused Gaussian-beam at the 1/e(sup 2) intensity, also referred to as spot size, is key to determining the fluence in laser processing experiments. Spot size measurements are also helpful to calculate the threshold energy and threshold fluence of a given material. This work reports an application of a conventional method, by analyzing single laser ablated spots for different laser pulse energies, to determine the cross-sectional area of a focused Gaussian-beam, which has a nominal pulse width of approx. 10 ps. Polished tungsten was used as the target material, due to its low surface roughness and low ablation threshold, to measure the beam waist diameter. From the ablative spot measurements, the ablation threshold fluence of the tungsten substrate was also calculated.

Ledesma, Rodolfo↗

SUPERPOSITION OF POLYTROPES IN THE INNER HELIOSHEATH

This paper presents a possible generalization of the equation of state and Bernoulli's integral when a superposition of polytropic processes applies in space and astrophysical plasmas. The theory of polytropic thermodynamic processes for a fixed polytropic index is extended for a superposition of polytropic indices. In general, the superposition may be described by any distribution of polytropic indices, but emphasis is placed on a Gaussian distribution. The polytropic density–temperature relation has been used in numerous analyses of space plasma data. This linear relation on a log–log scale is now generalized to a concave-downward parabola that is able to describe the observations better. The model of the Gaussian superposition of polytropes is successfully applied in the proton plasma of the inner heliosheath. The estimated mean polytropic index is near zero, indicating the dominance of isobaric thermodynamic processes in the sheath, similar to other previously published analyses. By computing Bernoulli's integral and applying its conservation along the equator of the inner heliosheath, the magnetic field in the inner heliosheath is estimated, B ∼ 2.29 ± 0.16 μG. The constructed normalized histogram of the values of the magnetic field is similar to that derived from a different method that uses the concept of large-scale quantization, bringing incredible insights to this novel theory.

79 ASTRONOMY AND ASTROPHYSICS↗

Exploring the energy landscape of RBMs: reciprocal space insights into bosons, hierarchical learning and symmetry breaking

Deep generative models have become ubiquitous due to their ability to learn and sample from complex distributions. Despite the proliferation of various frameworks, the relationships among these models remain largely unexplored, a gap that hinders the development of a unified theory of AI learning. In this work, we address two central challenges: clarifying the connections between different deep generative models and deepening our understanding of their learning mechanisms. We focus on Restricted Boltzmann Machines (RBMs), a class of generative models known for their universal approximation capabilities for discrete distributions. By introducing a reciprocal space formulation for RBMs, we reveal a connection between these models, diffusion processes, and systems of coupled bosons. Our analysis shows that at initialization, the RBM operates at a saddle point, where the local curvature is determined by the singular values of the weight matrix, whose distribution follows the Marc̆enko-Pastur law and exhibits rotational symmetry. During training, this rotational symmetry is broken due to hierarchical learning, where different degrees of freedom progressively capture features at multiple levels of abstraction. This leads to a symmetry breaking in the energy landscape, reminiscent of Landau’s theory. This symmetry breaking in the energy landscape is characterized by the singular values and the weight matrix eigenvector matrix. We derive the corresponding free energy in a mean-field approximation. We show that in the limit of infinite size RBM, the reciprocal variables are Gaussian distributed. Our findings indicate that in this regime, there will be some modes for which the diffusion process will not converge to the Boltzmann distribution. To illustrate our results, we trained replicas of RBMs with different hidden layer sizes using the MNIST dataset. Our findings not only bridge the gap between disparate generative frameworks but also shed light on the fundamental processes underpinning learning in deep generative models.

97 MATHEMATICS AND COMPUTING↗

Quantitative Mapping of Nanoscale Chemical Dynamics in Sub-sampled Operando (S)TEM Images Using Spatio-Temporal Analytics

One of the main limitations in the use of operando scanning transmission electron microscopes to study dynamic chemical processes is the effect of the electron beam on the kinetics of the reaction being observed. Here we demonstrate that a flexible Gaussian mixture model can be used to extract quantitative information directly from sub-sampled images, i.e.~images where not all the pixels in the image are illuminated with the beam. The use of this method is demonstrated for the charge/discharge cycle of a Li-battery, where the deposition and dissolution of Li at the anode can be accurately quantified for sub-sampling levels down to 1%. Performing operando imaging using a small fraction of the pixels means that the observations can significantly reduce the effect of the beam, automatically increase the imaging speed and decrease the total data transfer rate required. Such new software capabilities offer the potential to significantly widen the application of operando hardware approaches to study nanoscale dynamics in materials.

Stanfill, Bryan A.↗

Femtosecond laser-induced periodic surface structure on the Ti-based nanolayered thin films

Laser-induced periodic surface structures (LIPSSs) and chemical composition changes of Ti-based nanolayered thin films (Al/Ti, Ni/Ti) after femtosecond (fs) laser pulses action were studied. Irradiation is performed using linearly polarized Ti:Sapphire fs laser pulses of 40 fs pulse duration and 800 nm wavelength. The low spatial frequency LIPSS (LSFL), oriented perpendicular to the laser polarization with periods slightly lower than the irradiation wavelength, was typically formed at elevated laser fluences. On the contrary, high spatial frequency LIPSS (HSFL) with uniform period of 155 nm, parallel to the laser light polarization, appeared at low laser fluences, as well as in the wings of the Gaussian laser beam distribution for higher used fluence. LSFL formation was associated with the material ablation process and accompanied by the intense formation of nanoparticles, especially in the Ni/Ti system. The composition changes at the surface of both multilayer systems in the LSFL area indicated the intermixing between layers and the substrate. Concentration and distribution of all constitutive elements in the irradiated area with formed HSFLs were almost unchanged.

36 MATERIALS SCIENCE↗

Residual Optically Stimulated Luminescent (OSL) Signals for Al 2 O 3 :C and a Readout System with Reproducible Partial Signal Clearance

Optically stimulated luminescent dosimeters (OSLDs) are devices that, when stimulated with light, emit light in proportion to the integrated ionizing radiation dose. The stimulation of OSL material results in the loss of a small fraction of signal stored within the dosimetric traps. Previous studies have investigated the signal loss due to readout stimulation, and the optical annealing of OSLDs. This study builds on former research by examining the behavior of OSL signal after annealing, exploring the functionality of a previously developed signal loss model, and comparing uncertainties for dosimeters reused with or without annealing. For a completely annealed dosimeter, the minimum signal level was 56 ± 8 counts, and readings followed a Gaussian distribution. For dosimeters above this signal level, the fractional signal loss due to the reading process has a linear relationship with the calculated signal. At low signal levels (below 20,000 counts), in this OSLD system,calculated signal percent errors increase significantly but otherwise are on average 0.72 ± 0.27%, 0.40 ± 0.19%, 0.33±0.12%, and 0.24±0.07% for 30, 75, 150, and 300 readings. Theoretical calculations of uncertainties showed that annealing before reusing dosimeters allows for dose errors below 1% with as few as 30 readings. Reusing dosimeters multiple times increases the dose errors especially with low numbers of readouts, so theoretically around 300 readings would be necessary to achieve errors around 1% or below in most scenarios. Note that these dose errors do not include the error associated with the signal-to-dose conversion factor.

61 RADIATION PROTECTION AND DOSIMETRY↗

Characterization of groove density variation of VLS gratings with ALS XROL LTP-II in different operation modes

The long trace profiler, LTP-II, available at the Advanced Light Source (ALS) X-ray Optics Laboratory (XROL), was recently upgraded by replacing a multimode diode laser light source with a single-mode, wavelength-stabilized, fibercoupled diode laser system. The upgrade enables us to reliably characterize the lateral variation of groove density of variable-line-spacing (VLS) x-ray diffraction gratings. Here, we discuss the LTP-II performance with an example of measurements with a VLS grating with the groove density at the grating center of 300 lines/mm. For the measurements, we use the LTP-II in two different operation arrangements, the single Gaussian beam and the pencil beam interferometer arrangements. For each operation arrangement, we apply two data processing algorithms: with calculating the centroid position and with determining the position of a characteristic features of the detected beam intensity distributions. Additionally, we discuss the observed strong correlation between the LTP-II modes of operation and the resulted (extracted) groove density variations. We also speculate on possible origin of the correlation.

42 ENGINEERING↗

Score-Based Physics-Informed Neural Networks for High-Dimensional Fokker–Planck Equations

The Fokker-Planck (FP) equation is a foundational partial differential equation (PDE) in stochastic processes involving Brownian motions. However, the curse of dimensionality (CoD) poses a formidable challenge when dealing with high-dimensional FP equations. Although Monte Carlo simulation and (vanilla) Physics-Informed Neural Networks (PINNs) have shown the potential to tackle CoD, both methods exhibit significant numerical errors in high dimensions when dealing with the probability density function (PDF) associated with Brownian motion. The point-wise PDF values tend to decrease exponentially as dimensionality increases, surpassing the precision of numerical simulations and resulting in substantial errors. In addition, due to its massive sampling, Monte Carlo fails to offer fast sampling. Modeling the logarithm likelihood (LL) via vanilla PINNs transforms the FP equation into a notoriously difficult Hamilton-Jacobi-Bellman (HJB) equation, which is impractical for PINN learning, whose error grows rapidly with dimension. To this end, we propose a novel approach utilizing a score-based solver to fit the score function in stochastic differential equations (SDEs). The score function, defined as the gradient of the LL, plays a fundamental role in inferring LL and PDF and enables fast SDE sampling, offering an effective means to overcome the CoD. Three fitting methods, Score Matching (SM), Sliced Score Matching (SSM), and Score-PINN, are introduced, each contributing unique advantages in computational complexity, accuracy, and generality. The proposed score-based SDE solver operates in two stages: first, employing score matching or Score-PINN to acquire the score function; and second, solving the LL via an ordinary differential equation (ODE) using the obtained score function. Comparative evaluations across these methods showcase varying trade-offs. The proposed methodology is evaluated across diverse SDEs, including anisotropic Ornstein-Uhlenbeck processes, geometric Brownian motion, and Brownian motion with varying eigenspace. We also test various distributions, including Gaussian, Log-normal, Laplace, and Cauchy distributions. The numerical results demonstrate the score-based SDE solver’s stability, speed, and performance across different experimental settings, solidifying its potential as a solution to CoD for high-dimensional FP equations.

97 MATHEMATICS AND COMPUTING↗

A General Method for Estimating Bulk 2D Projections of Ice Particle Shape: Theory and Applications

The orientation of falling ice particles directly influences estimates of microphysical and radiative bulk quantities as well as in situ retrievals of size, shape, and mass. However, retrieval efforts and bulk calculations often incorporate very basic orientations or ignore these effects altogether. To address this deficiency, this study develops a general method for projecting bulk distributions of particle shape for arbitrary orientations. The Amoroso distribution provides the most general bulk aspect ratio distribution for gamma-distributed particle axis lengths. The parameters that govern the behavior of this aspect ratio distribution depend on the assumed relationship between mass, maximum dimension, and aspect ratio. Individual spheroidal geometry allows for eccentricity quantities to linearly map onto ellipse analogs, whereas aspect ratio quantities map nonlinearly. For particles viewed from their side, this analytic distinction leads to substantially larger errors in projected aspect ratio than for projected eccentricity. Distribution transformations using these mapping equations and numerical integration of projection kernels show that both truncation of size distributions and changes in Gaussian dispersion can alter the modality and shape of projection distributions. As a result, the projection process can more than triple the relative entropy between the spheroidal and projection distributions for commonly assumed model and orientation parameters. This shape uncertainty is maximized for distributions of highly eccentric particles and for particles like aggregates that are thought to fall with large canting-angle deviations. As a result, the methods used to report projected aspect ratios and the corresponding values should be questioned.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Preliminary results from the COBE differential microwave radiometers - Large angular scale isotropy of the cosmic microwave background

Preliminary but precise micowave maps are presented of the sky, and thus of the early universe, derived as the first results from the Differential Microwave Radiometers experiment aboard COBE. The dipole anisotropy attributed to the motion of the solar system with respect to the CMB reference frame shows strongly in all six sky maps and is consistent with a Doppler-shifted thermal spectrum. The best-fitted dipole has amplitude 3.3 + or - 0.2 mK in the direction (alpha, delta) = 11.2 h + or - 0.2 h, -7 deg + or - 2 deg (J2000) or (l,b) = 265 deg + or - 2 deg, 48 deg + or - 2 deg. There is no clear evidence in the maps for any other large angular-scale feature. Limits on Delta T/T0 of 3 x 10 to the -5th (T0 = 2.735 K), 4 x 10 to the -5th, and 4 x 10 to the -5th are found for the rms quadrupole amplitude, monochromatic fluctuations, and Gaussian fluctuations, respectively. These measurements place the most severe constraints to date on many potential physical processes in the early universe.

Smoot, G. F.↗