Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kernel”

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,351 records · Page 75

Data-enabled structure–property mappings for lanthanide-activated inorganic scintillators

The use of machine learning (ML) models toward development of validated structure–property relationships for two fundamental properties of activated inorganic scintillators for high energy radiation detection, namely the light yield (LY) and the decay time constant, is explored. The ML models are built on easily accessible proxies of materials—interchangeably referred to as features, descriptors or fingerprints—that are carefully selected on the basis of a physical understanding of the scintillation mechanism. Our study indicates that the developed physics-based ML models employing kernel ridge regression (KRR) and AdaBoost algorithm applied on top of a decision tree-based regression are able to “learn” the underlying design rules in a multi-dimensional feature space and thereby enable reasonably accurate predictions of the two target properties on unseen compounds (i.e., on a held-out test set). For instance, within a set of twenty-five cerium- or europium-doped scintillator materials, our analysis reveals a strong correlation between the average ionic part of the dielectric constant and the LY, irrespective of the specific chemistry of the compounds, indicating that the average ionic part of the dielectric constant is a particularly relevant descriptor toward prediction of the LY. Our results also demonstrate that, despite the use of small training datasets, the developed models are able to quickly distinguish high performing chemistries from those with relatively poor performance and therefore can play a crucial role in screening of new compounds with an attractive combination of targeted properties. Finally, the present study provides necessary motivation for future efforts involving ML models with relatively large training datasets, vast feature space explorations, and experimental design in search of promising novel scintillator chemistries.

36 MATERIALS SCIENCE↗

Method of Green’s Potentials for Elliptic PDEs in Domains with Random Apertures

Problems with topological uncertainties appear in many fields ranging from nano-device engineering to the design of bridges. In many of such problems, a part of the domains boundaries is subjected to random perturbations making inefficient conventional schemes that rely on discretization of the whole domain. Here, we study elliptic PDEs in domains with boundaries comprised of a deterministic part and random apertures, and apply the method of modified potentials with Green’s kernels defined on the deterministic part of the domain. This approach allows to reduce the dimension of the original differential problem by reformulating it as a boundary integral equation posed on the random apertures only. The multilevel Monte Carlo method is then applied to this modified integral equation and its optimal ϵ -2 asymptotical complexity is shown. Finally, we provide the qualitative analysis of the proposed technique and support it with numerical results.

97 MATHEMATICS AND COMPUTING↗

Improved Architectures and Training Algorithms for Deep Operator Networks

Operator learning techniques have recently emerged as a powerful tool for learning maps between infinite-dimensional Banach spaces. Trained under appropriate constraints, they can also be effective in learning the solution operator of partial differential equations (PDEs) in an entirely self-supervised manner. In this work we analyze the training dynamics of deep operator networks (DeepONets) through the lens of Neural Tangent Kernel theory, and reveal a bias that favors the approximation of functions with larger magnitudes. To correct this bias we propose to adaptively re-weight the importance of each training example, and demonstrate how this procedure can effectively balance the magnitude of back-propagated gradients during training via gradient descent. We also propose a novel network architecture that is more resilient to vanishing gradient pathologies. Taken together, our developments provide new insights into the training of DeepONets and consistently improve their predictive accuracy by a factor of 10-50x, demonstrated in the challenging setting of learning PDE solution operators in the absence of paired input-output observations.

97 MATHEMATICS AND COMPUTING↗

Persistent Sampling: Enhancing the Efficiency of Sequential Monte Carlo

Sequential Monte Carlo (SMC) samplers are powerful tools for Bayesian inference but suffer from high computational costs due to their reliance on large particle ensembles for accurate estimates. We introduce persistent sampling (PS), an extension of SMC that systematically retains and reuses particles from all prior iterations to construct a growing, weighted ensemble. By leveraging multiple importance sampling and resampling from a mixture of historical distributions, PS mitigates the need for excessively large particle counts, directly addressing key limitations of SMC such as particle impoverishment and mode collapse. Crucially, PS achieves this without additional likelihood evaluations-weights for persistent particles are computed using cached likelihood values. This framework not only yields more accurate posterior approximations but also produces marginal likelihood estimates with significantly lower variance, enhancing reliability in model comparison. Furthermore, the persistent ensemble enables efficient adaptation of transition kernels by leveraging a larger, decorrelated particle pool. Experiments on high-dimensional Gaussian mixtures, hierarchical models, and non-convex targets demonstrate that PS consistently outperforms standard SMC and related variants, including recycled and waste-free SMC, achieving substantial reductions in mean squared error for posterior expectations and evidence estimates, all at reduced computational cost. PS thus establishes itself as a robust, scalable, and efficient alternative for complex Bayesian inference tasks.

Karamanis, Minas↗

Nekbone performance on GPUs with OpenACC and CUDA Fortran implementations

We present a hybrid GPU implementation and performance analysis of Nekbone, which represents one of the core kernels of the incompressible Navier–Stokes solver Nek5000. The implementation is based on OpenACC and CUDA Fortran for local parallelization of the compute-intensive matrix–matrix multiplication part, which significantly minimizes the modification of the existing CPU code while extending the simulation capability of the code to GPU architectures. Our discussion includes the GPU results of OpenACC interoperating with CUDA Fortran and the gather–scatter operations with GPUDirect communication. We demonstrate performance of up to 552 Tflops on 16, 384 GPUs of the OLCF Cray XK7 Titan.

Computer Science↗

Quantitative Studies of Machining-Induced Microstructure Alteration and Plastic Deformation in AISI 316 Stainless Steel Using EBSD

Surface machining can result in a heavily plastic deformed layer in the subsurface. The microstructure alterations and the extent of plastic deformation of a milled 316 austenitic stainless steel were investigated using electron backscatter diffraction (EBSD) technique in this study. The microstructure alteration of the deformed layer was characterized by the generation of grain and subgrain boundaries. The plastic deformation was evaluated by kernel average misorientation (KAM) and grain average misorientation (GAM). The local deformation within grains was assessed by the misorientation from adjacent point (MFAP). The EBSD results indicated gradients of microstructure and plastic deformation along the depth direction. High-angle grain boundaries were generated within the depth of 40 μm, and subgrain boundaries were produced up to the depth of 150 μm. Estimated by KAM and GAM, the plastic strain was about 0.7 at the depth of 25 μm and decreased to 0.1 at the depth of 115 μm The MFAP analysis revealed that the local deformation varied within individual grains. Additionally, the change of micro-hardness along the depth direction showed a strong correlation with KAM and GAM values.

36 MATERIALS SCIENCE↗

Influence of Thermomechanically Controlled Processing on Microstructure and Hydrogen Induced Cracking Susceptibility of API 5L X70 Pipeline Steel

The effect of different thermomechanical controlled processing routes on susceptibility of X70 pipeline steel to hydrogen induced cracking (HIC) have been studied. Two X70 pipeline steel specimens labelled WE and WD were investigated. These specimens have the same chemical composition, but they were processed with seperate thermomechanical treatments parameters. Microstructural examinations showed that WE consists of mainly acicular ferrite and polygonal ferrite, while WD consists of acicular ferrite and bainitic ferrite. After subjecting both specimens to hydrogen charging for 12 and 16 h in 0.2 M sulfuric acid and 3 g/L ammonium thiocyanate, early onset of HIC was observed in specimen WD. Post-hydrogen charging microstructural evaluation showed the nucleation of discontinuous cracks in WD after 12 h of charging. However, extended charging for up to 16 h resulted in HIC along the mid-thickness region of both specimens. Hydrogen diffusion across specimen WE was better than that of specimen WD. Therefore, hydrogen trapping at grain boundaries, banded deformed grains, inclusions and secondary phases such as martensite and cementite aided initiation and propagation of HIC in specimens. Nevertheless, the adverse effect of these features on HIC risks was more prominent in specimen WD compared to specimen WE. The Vickers microhardness values measured in WD (349.6 HV) and WE (307.4 HV) suggest that WD is harder than WE; and higher kernel average misorientation of 0.66° in WD than in WE (0.58°) shows higher dislocation density in WD. The results from slow strain rate tensile test confirmed that specimen WD was stronger and more susceptible to HIC than specimen WE. It was concluded that microstructural phases developed during thermomechanical processing improved strength in WD at the expense of its crack resistance, while WE with lower strength showed more ductility and higher resistance to HIC.

36 MATERIALS SCIENCE↗

Electron Beam Welding of Nimonic 80A Superalloy: Microstructure Evolution and EBSD Study After Aging Heat Treatment

The purpose of this research was to evaluate integrity and microstructure changes of the electron beam-welded Nimonic 80A and the base metal after aging heat treatment process, i.e., 700 °C/72 h/air cooling. Here, microstructure, microhardness and tensile were investigated by various tools, i.e., scanning electron microscope with energy-dispersive spectroscopy and electron backscatter diffraction, x-ray diffraction analysis and Vickers hardness. Evaluation of microstructure shows that the weld zone mainly consists of columnar dendrites. The effect of aging heat treatment on the microstructure morphology and phase structure is very limited. Also, it was observed that there are obvious segregation and precipitates in intracrystalline and along the grain boundaries after aging heat treatment. Electron backscatter diffraction results show that the weld zone and base metal have a clear difference in the grain boundary characterization and kernel average misorientation at room temperature and after aging heat treatment, which leads to a difference in microhardness of the weld zone. In addition, the relationship between microhardness and tensile strength is discussed.

36 MATERIALS SCIENCE↗

Treatment of near-incompressibility in meshfree and immersed-particle methods

We propose here new projection methods for treating near-incompressibility in small and large deformation elasticity and plasticity within the framework of particle and meshfree methods. Using the $\bar{B}$ and $\bar{F}$ techniques as our point of departure, we develop projection methods for the conforming reproducing kernel method and the immersed-particle or material point-like methods. The methods are based on the projection of the dilatational part of the appropriate measure of deformation onto lower-dimensional approximation spaces, according to the traditional $\bar{B}$ and $\bar{F}$ approaches, but tailored to meshfree and particle methods. The presented numerical examples exhibit reduced stress oscillations and are free of volumetric locking and hourglassing phenomena.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A MUSCL-SCNI approach for meshfree modeling of shock waves in fluids

Here, a stable and nodally integrated meshfree formulation for modeling shock waves in fluids is developed. The reproducing kernel approximation is employed to discretize the conservation equations for compressible flow, and a flux vector splitting approach is applied to allow proper numerical treatments for the advection and pressure parts, respectively, based on the characteristics of each flux term. To capture the essential shock physics in fluids, including the Rankine–Hugoniot jump conditions and the entropy condition, local Riemann enrichment is introduced under the stabilized conforming nodal integration (SCNI) framework. Meanwhile, numerical instabilities associated with the advection flux are eliminated by adopting a modified upwind scheme. To further enhance accuracy, a MUSCL-type method is introduced in conjunction with an oscillation limiter to avoid Gibbs phenomenon and ensure monotonic piecewise linear reconstruction in the smooth region. The present meshfree formulation is free from tunable artificial parameters and is capable of capturing shock and rarefaction waves without over/undershoots. Finally, several numerical examples are analyzed to demonstrate the effectiveness of the proposed MUSCL-SCNI approach in meshfree modeling of complex shock phenomena, including shock diffraction, shock–vortex interaction, and high energy explosion processes.

42 ENGINEERING↗

Exploring temporal community evolution: algorithmic approaches and parallel optimization for dynamic community detection

Abstract Dynamic (temporal) graphs are a convenient mathematical abstraction for many practical complex systems including social contacts, business transactions, and computer communications. Community discovery is an extensively used graph analysis kernel with rich literature for static graphs. However, community discovery in a dynamic setting is challenging for two specific reasons. Firstly, the notion of temporal community lacks a widely accepted formalization, and only limited work exists on understanding how communities emerge over time. Secondly, the added temporal dimension along with the sheer size of modern graph data necessitates new scalable algorithms. In this paper, we investigate how communities evolve over time based on several graph metrics under a temporal formalization. We compare six different algorithmic approaches for dynamic community detection for their quality and runtime. We identify that a vertex-centric (local) optimization method works as efficiently as the classical modularity-based methods. To its advantage, such local computation allows for the efficient design of parallel algorithms without incurring a significant parallel overhead. Based on this insight, we design a shared-memory parallel algorithm DyComPar , which demonstrates between 4 and 18 fold speed-up on a multi-core machine with 20 threads, for several real-world and synthetic graphs from different domains.

97 MATHEMATICS AND COMPUTING↗

Effects of heat treatment and build orientation on the evolution of $ϵ$ and $α'$ martensite and strength during compressive loading of additively manufactured 304L stainless steel

In this paper, the effect of heat-treatment and build orientation on martensitic phase transformation in additively manufactured (AM) 304L stainless steel is studied and compared with conventionally produced wrought material. The relationships between observed martensitic transformations and material microstructures and their effects on mechanical strength are established through experimental observations. In situ high-energy X-ray powder diffraction measurements were performed to monitor the evolution of $ϵ$ and $α'$ martensite during compressive loading of stainless steel. Electron backscatter diffraction (EBSD) was used to provide insight on initial grain morphology, crystallographic misorientation within grains, and crystallographic texture. Heat treatment alters the microstructure of AM samples creating different initial conditions. This difference in starting microstructure resulted in variability in martensitic transformation during compressive deformation. The rate of martensitic transformation decreased for AM samples treated with temperatures up to 1100°C, after which the AM microstructures recrystallized, resulting in increased rate of martensitic transformation for those samples treated at higher temperatures. It was also observed that aligning the axis of compression with the AM build direction resulted in a lower rate of strain-induced martensite formation as opposed to aligning the compression axis perpendicular to it. More favorable distribution of crystal orientations in the latter loading orientation promoted martensitic transformation. These and additional experimental observations from EBSD in terms of kernel average misorientation, mean grain orientation spread, and mean crystallite size reveal strong microstructural effects on strength of additively manufactured metallic materials.

304L stainless steel↗

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↗

Emerging hierarchical dislocation structures: Insights from scanning electron microscopy-electron backscatter diffraction in situ tensile testing and multifractal analysis

Understanding the evolution of dislocation structures during plastic deformation is critical for predicting the mechanical performance of metallic materials. In this work, we applied in situ scanning electron microscopy/electron backscatter diffraction tensile testing combined with multifractal (MF) analysis to assess deformation-induced dislocation structure evolution in solution-annealed 304 L stainless steel, both in its as-received and neutron-irradiated states (5.4 displacements per atom). The analysis of kernel average misorientation patterns revealed the formation of hierarchical dislocation arrangements that exhibit clear MF scaling behavior. Despite pronounced visual differences between nonirradiated and irradiated specimens—most notably, the appearance of dislocation channels after irradiation—the singularity spectra suggest that both conditions give rise to similar underlying hierarchical structures. MF analysis provides a quantitative measure of the spatial complexity and self-organization of dislocation patterns, highlighting the accelerated emergence and evolution of the dislocation structures in irradiated polycrystalline materials, as well as the limitation of their spatial extent. The findings indicate that irradiation not only modifies microstructure but also alters correlation-driven dislocation organization. More generally, they demonstrate that MF analysis is a powerful tool for probing mesoscale deformation mechanisms.

Dislocation structures↗

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗

Thermal scattering law of C 2 H 4 n : Integrating experimental data with DFT calculations

Improvements in determination of the thermal scattering law of moderator materials (measuring, calculating and validating) are important for accurate prediction of neutron thermalization in nuclear systems. In this paper a methodology for producing thermal scattering libraries from the experimental data for polyethylene C 2 H 4 n is discussed. Double differential scattering cross section (DDSCS) experiments were performed at the Spallation Neutron Source of Oak Ridge National Laboratory (SNS ORNL). New scattering kernel evaluations, based on phonon spectrum for C 2 H 4 n , are created using the NJOY2016 code. Two different methods were used: direct and indirect geometry neutron scattering at ARCS and SEQUOIA, and VISION instruments, respectively, where the phonon spectrum was derived from the dynamical structure factor S(Q, ω ) obtained from the measured DDSCS. In order to compare and validate the newly created library, the experimental setup was simulated using MCNP6.1. Compared with the current ENDF/B-VII.1, the resulting RPI C 2 H 4 n libraries improved both double differential scattering and total scattering cross sections. A set of criticality benchmarks containing C 2 H 4 n from HEU-MET-THERM resulted in an overall improved calculation of K eff , although the libraries should be tested against benchmarks more sensitive to C 2 H 4 n . The DFT + oClimax method is used and is shown to be most comprehensive method for analysis of moderator materials. Finally, the importance of DFT + oClimax method lies in the fact that it can be validated against all data measured at VISION, ARCS and SEQUOIA, and experimental total scattering cross section measurements.

74 ATOMIC AND MOLECULAR PHYSICS↗

Experimental validation of the mechanistic scale-up methodology of gas–solid spouted beds using radioactive particle tracking (RPT)

The very high-temperature reactors (VHTRs) are highly ranked among candidates of Generation IV of nuclear reactors due to their high efficiency, safety, the resistance to proliferation, and reliability. The VHTRs are preferentially fueled by Tristructural-isometric (TRISO) coated fuel particles which has fuel kernels of fissile material coated by four coating layers: a porous buffer pyrolysis carbon layer (buffer PyC), an inner dense pyrolysis carbon layer (IPyC), a silicon carbide layer (SiC) and an outer dense pyrolysis carbon layer (OPyC). The heart of the operation and safety of the VHTRs significantly depends on the reliability of the coating layers of TRISO particles to retain metallic and gaseous fission products within the particles. The technique used for coating TRISO particles are gas-solids spouted beds via chemical vapor deposition (CVD). Fabrication of high-quality low-defect TRISO fuel particles fuel at larger scale spouted beds is required to support the commercialization of the VHTRs. In this work, our new developed mechanistic scale-up methodology of gas-solids spouted beds based on matching the radial profile of gas-holdup has been demonstrated and validated using radioactive particle tracking (RPT). Two spouted beds of small and larger scales were used in the study. Three sets of conditions were carried out which include the conditions of the reference case in the large scale, conditions that provide similar gas holdup radial profile to that of the reference case and conditions that provided dissimilar gas holdup radial in the small-scale spouted beds. The results confirm the validation of the scale-up methodology in terms of the dimensionless values of the spout diameter, cumulative probability distribution of the solids particles penetration into the spout, fraction of cycle time in each region of the bed, the radial profiles of the dimensionless values of the root-mean-square particle velocities and solids eddy diffusivity. Finally, the results further advance the knowledge and understanding of the gas-solids spouted beds provide deeper insight into their solids dynamics and presenting important benchmarking data for validating computational fluid dynamics codes and models. At last, procedures are established for the implementation of the new scale-up methodology.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗