Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kernel function”

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 757 records · Page 42

Improved Cross Section Generation Capability of Griffin in FY22

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of deterministic steady-state transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication as well as transient solvers for spatial dynamics with the improved quasi-static method. Griffin uses cross section data in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. In recent years, the MC 2 -3 modules have been added to Griffin for fast reactor cross section generation, and the self-shielding application programming interface (SSAPI) was implemented in the ISOXML module for thermal reactor cross section generation. The on-the-fly slowing down method and double-heterogeneity treatment have been implemented to SSAPI and verified against particulate fuel-bearing graphite-moderated thermal reactor problems with high accuracy. This year, work has been focused on improving the cross section generation capability of ISOXML and streamlining the cross section generation procedures. In addition, the form function data were added to ISOXML in order to support the pin power reconstruction capability that was newly implemented in Griffin in this fiscal year. To facilitate the cross section generation using MC 2 -3 and SSAPI in Griffin, the cross section generation workflows have been set up for both fast and thermal spectrum reactors. The MOOSE action system tool was devised for fast spectrum problems, and the MOOSE stochastic tool was adopted to the branch calculation procedure for thermal spectrum problems. Meanwhile, to ensure the accuracy of group-constants, the thermal up-scattering kernel calculator accounting for resonance scattering was implemented in ISOXML, demonstrating the accurate computation of a Doppler-broadened scattering kernel of any Legendre order within a reasonable timescale. Other aspects of ISOXML, such as deletion solver and data, documentation, ISOXML file management, and the interface for Mixture, were improved as well.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Nonlinear Matrix Approximation with Radial Basis Function Components

We introduce and investigate matrix approximation by decomposition into a sum of radial basis function (RBF) components. An RBF component is a generalization of the outer product between a pair of vectors, where an RBF function replaces the scalar multiplication between individual vector elements. Even though the RBF functions are positive definite, the summation across components is not restricted to convex combinations and allows us to compute the decomposition for any real matrix that is not necessarily symmetric or positive definite. We formulate the problem of seeking such a decomposition as an optimization problem with a nonlinear and non-convex loss function. Several modern versions of the gradient descent method, including their scalable stochastic counterparts, are used to solve this problem. We provide extensive empirical evidence of the effectiveness of the RBF decomposition and that of the gradient-based fitting algorithm. While being conceptually motivated by singular value decomposition (SVD), our proposed nonlinear counterpart outperforms SVD by drastically reducing the memory required to approximate a data matrix with the same L2 error for a wide range of matrix types. For example, it leads to 2 to 6 times memory save for Gaussian noise, graph adjacency matrices, and kernel matrices. Moreover, this proximity-based decomposition can offer additional interpretability in applications that involve, e.g., capturing the inner low-dimensional structure of the data, retaining graph connectivity structure, and preserving the acutance of images.

Rebrova, Elizaveta↗

Hierarchical Gaussian process-based Bayesian optimization for materials discovery in high entropy alloy spaces

Bayesian optimization (BO) is a powerful and data-efficient method for iterative materials discovery and design, particularly valuable when prior knowledge is limited, underlying functional relationships are complex or unknown, and the cost of querying the materials space is significant. Traditional BO methodologies typically utilize conventional Gaussian Processes (cGPs) to model the relationships between material inputs and properties, as well as correlations within the input space. However, cGP-BO approaches often fall short in multi-objective optimization scenarios, where they are unable to fully exploit correlations between distinct material properties. Leveraging these correlations can significantly enhance the discovery process, as information about one property can inform and improve predictions about others. Here, this study addresses this limitation by employing advanced kernel structures to capture and model multi-dimensional property correlations through multi-task (MTGPs) or deep Gaussian Processes (DGPs), thus accelerating the discovery process. We demonstrate the effectiveness of MTGP-BO and DGP-BO in rapidly and robustly solving complex materials design challenges that occur within the context of complex multi-objective optimization over FCC FeCrNiCoCu high entropy alloy (HEA) spaces, where traditional cGP-BO approaches fail. Furthermore, we highlight how the differential costs associated with querying various material properties can be strategically leveraged to make the materials discovery process more cost-efficient.

36 MATERIALS SCIENCE↗

Stress-intensity factors for a thick-walled cylinder containing an annular imbedded or external or internal surface crack

The elastostatic axisymmetric problem for a long thick-walled cylinder containing a ring-shaped internal or edge crack is considered. Using the standard transform technique the problem is formulated in terms of an integral equation which has a simple Cauchy kernel for the internal crack and a generalized Cauchy kernel for the edge crack as the dominant part. As examples the uniform axial load and the steady-state thermal stress problems have been solved and the related stress intensity factors have been calculated. Among other findings the results show that in the cylinder under uniform axial stress containing an internal crack the stress intensity factor at the inner tip is always greater than that at the outer tip for equal net ligament thicknesses and in the cylinder with an edge crack which is under a state of thermal stress the stress intensity factor is a decreasing function of the crack depth, tending to zero as the crack depth approaches the wall thickness.

Erdol, R.↗

Runtime Thread-Block Optimization for Custom Multistream CUDA Kernels for the Glenn Research Center Communication Analysis Suite

In preparation of the return of humans to the Moon with the coming Artemis missions, NASA scientists must evaluate proposed landing site locations for terrain and communications viability. The Glenn Research Center Communication Analysis Suite (GCAS) combines sophisticated communication network models with accurate lunar terrain to access sites across the Moon’s south pole. Given the importance of proper site selection to crew safety and mission success, many locations need to be analyzed resulting in a large computational load needing to be performed. To meet the growing project demands, development has begun to improve the runtime efficiency of GCAS with GPU parallelization by way of multi-stream CUDA kernels. One of the most prominent factors in kernel optimization is the proper selection of thread-block dimensions in order to maximize the concurrent operation on the device. Typically, thread-block dimensions are optimized by hand requiring many stages of benchmarking and iteration. Additionally, given the main conditions to optimization are the physical GPU architecture and problem size, these optimal dimensions are non-portable and fragile in their scope. As such, a novel optimization routine was developed to generate the optimal thread-block dimensions during runtime with considerations to hardware specifications and problem size resolving the issues of portability and enabling the function of more dynamic routines.

Aden Bergstresser↗

Radiation effects induced in pin photodiodes by 40- and 85-MeV protons

PIN photodiodes were bombarded with 40- and 85-MeV protons to a fluence of 1.5 x 10 to the 11th power p/sq cm, and the resulting change in spectral response in the near infrared was determined. The photocurrent, dark current and pulse amplitude were measured as a function of proton fluence. Changes in these three measured properties are discussed in terms of changes in the diode's spectral response, minority carrier diffusion length and depletion width. A simple model of induced radiation effects is presented which is in good agreement with the experimental results. The model assumes that incident protons produce charged defects within the depletion region simulating donor type impurities.

Becher, J.↗

A thick-walled cylinder with an axisymmetric internal or edge crack

The paper considers the elastostatic axisymmetric problem for a long thick-walled cylinder containing a ring-shaped internal or edge crack. Using the standard transform technique the problem is formulated in terms of an integral equation which has a simple Cauchy kernel for the internal crack and a generalized Cauchy kernel for the edge crack as the dominant part. As examples the uniform axial load and the steady-state thermal stress problems have been solved and the related stress-intensity factors have been calculated. Among other findings the results show that in the cylinder under uniform axial stress containing an internal crack the stress-intensity factor at the inner tip is always greater than that at the outer tip for equal net ligament thicknesses and in the cylinder with an edge crack which is under a state of thermal stress the stress-intensity factor is a decreasing function of the crack depth, tending to zero as the crack depth approaches the wall thickness.

Erdol, R.↗

Control Architecture for Robotic Agent Command and Sensing

Control Architecture for Robotic Agent Command and Sensing (CARACaS) is a recent product of a continuing effort to develop architectures for controlling either a single autonomous robotic vehicle or multiple cooperating but otherwise autonomous robotic vehicles. CARACaS is potentially applicable to diverse robotic systems that could include aircraft, spacecraft, ground vehicles, surface water vessels, and/or underwater vessels. CARACaS incudes an integral combination of three coupled agents: a dynamic planning engine, a behavior engine, and a perception engine. The perception and dynamic planning en - gines are also coupled with a memory in the form of a world model. CARACaS is intended to satisfy the need for two major capabilities essential for proper functioning of an autonomous robotic system: a capability for deterministic reaction to unanticipated occurrences and a capability for re-planning in the face of changing goals, conditions, or resources. The behavior engine incorporates the multi-agent control architecture, called CAMPOUT, described in An Architecture for Controlling Multiple Robots (NPO-30345), NASA Tech Briefs, Vol. 28, No. 11 (November 2004), page 65. CAMPOUT is used to develop behavior-composition and -coordination mechanisms. Real-time process algebra operators are used to compose a behavior network for any given mission scenario. These operators afford a capability for producing a formally correct kernel of behaviors that guarantee predictable performance. By use of a method based on multi-objective decision theory (MODT), recommendations from multiple behaviors are combined to form a set of control actions that represents their consensus. In this approach, all behaviors contribute simultaneously to the control of the robotic system in a cooperative rather than a competitive manner. This approach guarantees a solution that is good enough with respect to resolution of complex, possibly conflicting goals within the constraints of the mission to be accomplished by the vehicle(s).

Huntsberger, Terrance↗

Fission Product and Actinide Distribution in the SiC Layer of Fissioned Triso Fuel Particles

Tri-structural Isotropic (TRISO) encapsulated fuel particles are used in next generation high-temperature gas-cooled nuclear reactors. The integrity of the fuel over a wide range of temperatures is integral to the safety of an operating nuclear reactor. It has been observed that fission products are released at high temperatures that simulate hypothesized accident scenarios. The particles under study were from the AGR-1 fuel development and qualification program. They consist of an enriched UCO fuel kernel surrounded by a dense carbon buffer layer (100 {mu}m), an inner pyrolytic carbon (IPyC) layer (40 {mu}m), a silicon carbide (SiC) layer (35 {mu}m), and an outer pyrolytic carbon (OPyC) layer (40 {mu}m). The SiC layer acts as the main diffusion barrier to fission products, but has been demonstrated to fail due to release of Ag through intact particles and local corrosion from metallic fission products such as Pd. In this work, we use transmission electron microscopy and diffraction techniques to study the structure and dispersion of release products within the SiC of TRISO particles as a function of annealing temperature used to simulate margin temperatures for in-reactor accident scenarios. In order to improve our understanding of fission product retention as a function of temperature, we have the unique opportunity to examine irradiated, fissioned TRISO particles to study post-irradiation effects on these materials. This data will contribute to current models used to predict the behavior of fission products in the SiC layer of TRISO fuels.

07 ISOTOPE AND RADIATION SOURCES↗

ENDF/B-VIII.1: Updated Nuclear Reaction Data Library for Science and Applications

The ENDF/B-VIII.1 library is the newest recommended evaluated nuclear data file by the Cross Section Evaluation Working Group (CSEWG) for use in nuclear science and technology applications, and incorporates advances made in the six years since the release of ENDF/B-VIII.0. Among key advances made are that the 239 Pu file was reevaluated by a joint international effort and that updated 16,18 O, 19 F, 28–30 Si, 50–54 Cr, 55 Mn, 54,56,57 Fe, 63,65 Cu, 139 La, 233,235,238 U, and 240,241 Pu neutron nuclear data from the IAEA coordinated INDEN collaboration were adopted. Over 60 neutron dosimetry cross sections were adopted from the IAEA's IRDFF-II library. In addition, the new library includes significant changes for 3 He, 6 Li, 9 Be, 51 V, 88 Sr, 103 Rh, 140,142 Ce, Dy, 181 Ta, Pt, 206–208 Pb, and 234,236 U neutron data, and new nuclear data for the photonuclear, charged-particle and atomic sublibraries. Numerous thermal neutron scattering kernels were reevaluated or provided for the very first time. On the covariance side, work was undertaken to introduce better uncertainty quantification standards and testing for nuclear data covariances. The significant effort to reevaluate important nuclides has reduced bias in the simulations of many integral experiments with particular progress noted for fluorine, copper, and stainless steel containing benchmarks. Data issues hindered the successful deployment of the previous ENDF/B-VIII.0 for commercial nuclear power applications in high burnup situations. These issues were addressed by improving the 238 U and 239,240,241 Pu evaluated data in the resonance region. The new library performance as a function of burnup is similar to the reference ENDF/B-VII.1 library.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

PIN photodiodes irradiated with 40- and 85-MeV protons

PIN photodiodes were bombarded with 40- and 85-MeV protons to a fluence of 1.5 x 10 to the 11th p/sq cm, and the resulting change in spectral response in the near infrared was determined. The photocurrent, dark current, and pulse amplitude were measured as a function of proton fluence. Changes in these three measured properties are discussed in terms of changes in the diode's spectral response, minority carrier diffusion length, and depletion width. A simple model of induced radiation effects is presented which is in good agreement with the experimental results. The model assumes that incident protons produce charged defects within the depletion region simulating donor type impurities.

Reft, C. S.↗

The impact of anisotropic redshift distributions on angular clustering

A leading way to constrain physical theories from cosmological observations is to test their predictions for the angular clustering statistics of matter tracers, a technique that is set to become ever more central with the next generation of large imaging surveys. Interpretation of this clustering requires knowledge of the projection kernel, or the redshift distribution of the sources, and the typical assumption is an isotropic redshift distribution for the objects. However, variations in the kernel are expected across the survey footprint due to photometric variations and residual observational systematic effects. Here, we develop the formalism for anisotropic projection and present several limiting cases that elucidate the key aspects. We quantify the impact of anisotropies in the redshift distribution on a general class of angular two-point statistics. In particular, we identify a mode-coupling effect that can add power to auto-correlations, including galaxy clustering and cosmic shear, and remove it from certain cross-correlations. If the projection anisotropy is primarily at large scales, the mode-coupling depends upon its variance as a function of redshift; furthermore, it is often of similar shape to the signal. In contrast, the cross-correlation of a field whose selection function is anisotropic with another one featuring no such variations — such as CMB lensing — is immune to these effects. We discuss explicitly several special cases of the general formalism including galaxy clustering, galaxy-galaxy lensing, cosmic shear and cross-correlations with CMB lensing, and publicly release a code to compute the biases.

79 ASTRONOMY AND ASTROPHYSICS↗

Notos - a galaxy tool to analyze CpN observed expected ratios for inferring DNA methylation types

Background: DNA methylation patterns store epigenetic information in the vast majority of eukaryotic species. The relatively high costs and technical challenges associated with the detection of DNA methylation however have created a bias in the number of methylation studies towards model organisms. Consequently, it remains challenging to infer kingdom-wide general rules about the functions and evolutionary conservation of DNA methylation. Methylated cytosine is often found in specific CpN dinucleotides, and the frequency distributions of, for instance, CpG observed/expected (CpG o/e) ratios have been used to infer DNA methylation types based on higher mutability of methylated CpG. Results: Predominantly model-based approaches essentially founded on mixtures of Gaussian distributions are currently used to investigate questions related to the number and position of modes of CpG o/e ratios. These approaches require the selection of an appropriate criterion for determining the best model and will fail if empirical distributions are complex or even merely moderately skewed. We use a kernel density estimation (KDE) based technique for robust and precise characterization of complex CpN o/e distributions without a priori assumptions about the underlying distributions. Conclusions: We show that KDE delivers robust descriptions of CpN o/e distributions. For straightforward processing, we have developed a Galaxy tool, called Notos and available at the ToolShed, that calculates these ratios of input FASTA files and fits a density to their empirical distribution. Based on the estimated density the number and shape of modes of the distribution is determined, providing a rational for the prediction of the number and the types of different methylation classes. Notos is written in R and Perl.

59 BASIC BIOLOGICAL SCIENCES↗

SU-G-BRC-02: A Novel Multi-Criteria Optimization Approach to Generate Deliverable Intensity-Modulated Radiation Therapy (IMRT) Treatment Plans

Purpose: To present a novel multi-criteria optimization (MCO) solution approach that generates treatment plans with deliverable apertures using column generation. Methods: We demonstrate our method with 10 locally advanced head-and-neck cancer cases retrospectively. In our MCO formulation, we defined an objective function for each structure in the treatment volume. This resulted in 9 objective functions, including 3 distinct objectives for primary target volume, high-risk and low-risk target volumes, 5 objectives for each of the organs-at-risk (OARs) (two parotid glands, spinal cord, brain stem and oral cavity), and one for the non-target non-OAR normal tissue. Conditional value-at-risk (CVaR) constraints were utilized to ensure at least certain fraction of the target volumes receiving the prescription doses. To directly generate deliverable plans, column generation algorithm was embedded within our MCO approach for aperture shape generation. Final dose distributions for all plans were generated using a Monte Carlo kernel-superposition dose calculation. We compared the MCO plans with the clinical plans, which were created by clinicians. Results: At least 95% target coverage was achieved by both MCO plans and clinical plans. However, the average conformity indices of clinical plans and the MCO plans were 1.95 and 1.35, respectively (31% reduction, p<0.01). Compared to the conventional clinical plan, the proposed MCO method achieved average reductions in left parotid mean dose of 5% (p=0.06), right parotid mean dose of 18% (p<0.01), oral cavity mean dose of 21% (p=0.03), spinal cord maximum dose of 20% (p<0.01), brain stem maximum dose of 61% (p<0.01), and normal tissue maximum dose of 5% (p<0.01), respectively. Conclusion: We demonstrated that the proposed MCO method was able to obtain deliverable IMRT treatment plans while achieving significant improvements in dosimetric plan quality.

60 APPLIED LIFE SCIENCES↗

Automatic tools for system testing

As spacecraft control and other space-related ground systems become increasingly complex, the effort required in testing and validation also increases. Implementation of a spacecraft control system normally involves a number of incremental deliveries. In addition kernel or general purpose software may also be involved, which must itself be considered in the integration and testing program. Tools can be used to assist this testing. These can reduce the effort required or alternatively they can ensure that for a given level of effort, a better job is done. Great benefit could be derived by automating certain types of testing (interactive software) which up to now has been performed manually at a terminal. This paper reports on an on-going study. The study examines means of automating spacecraft control system testing, evaluates relevant commercial tools and aims to prototype basic automatic testing functions.

Peccia, N. M.↗

Deep learning the properties of inorganic perovskites

We report the ability to accurately and quickly predict the stability of materials and their structural and electronic properties remains a grand challenge in materials science. Density functional theory is widely used as a means of predicting these material properties, but is known to be computationally expensive and scales as the cube of the number of electrons in the material’s unit cell. In this article, for a previously published dataset of inorganic perovskites, we show that a single neural network model using only the elemental properties of the compounds’ constituents can predict lattice constants to within 0.1 Å, heat of formation to within 0.2 eV, and band gaps to within 0.7 eV RMSE. We also compare the performance of the trained network to two widely used regression techniques, namely random forest and Kernel ridge regression, and find that the neural network’s predictions are more accurate for each of the properties. The simultaneous accurate prediction of multiple key properties of technologically relevant materials is promising for rational design and optimization in known and novel chemical spaces.

36 MATERIALS SCIENCE↗

The temperature and density structures of an X-ray flare during the decay phase

The X-ray flare of 9 August 1973 was characterized by a spatially small kernel structure which persisted throughout its duration. The decay phase of this flare was observed in the objective grating mode of the X-ray telescope aboard the Skylab. Data analysis was carried out by scanning the images with a microdensitometer, converting the density arrays to energy using laboratory film calibration data and taking cross sections of the energy images. The 9 August flare shows two distinct periods in its decay phase, involving both cooling and material loss. The objective grating observations reveal that the two phenomena are separated in time. During the earlier phase of the flare decay, the distribution of emission measure as a function of temperature is changing, the high temperature component of the distribution being depleted relative to the cooler body of plasma. As the decay continues, the emission measure distribution stabilizes and the flux diminishes as the amount of material at X-ray emitting temperatures decreases.

Silk, J. K.↗

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗