Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed data structures”

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 289 records · Page 16

Chemical cartography with LAMOST and Gaia reveal azimuthal and spiral structure in the Galactic disc

Chemical Cartography, or mapping, of our Galaxy has the potential to fully transform our view of its structure and formation. In this work, we use chemical cartography to explore the metallicity distribution of OBAF-type disc stars from the LAMOST survey and a complementary sample of disc giant stars from Gaia DR3. We use these samples to constrain the radial and vertical metallicity gradients across the Galactic disc. We also explore whether there are detectable azimuthal variations in the metallicity distribution on top of the radial gradient. For the OBAF-type star sample from LAMOST, we find a radial metallicity gradient of Δ[Fe/H]/ΔR ~–0.078 ± 0.001 dex kpc –1 in the plane of the disc and a vertical metallicity gradient of Δ[Fe/H]/ΔZ ~–0.15 ± 0.01 dex kpc –1 in the solar neighbourhood. The radial gradient becomes shallower with increasing vertical height, while the vertical gradient becomes shallower with increasing Galactocentric radius, consistent with other studies. We also find detectable spatially dependent azimuthal variations on top of the radial metallicity gradient at the level of ~0.10 dex. Interestingly, the azimuthal variations appear be close to the Galactic spiral arms in one data set (Gaia DR3) but not the other (LAMOST). Furthermore, these results suggest that there is azimuthal structure in the Galactic metallicity distribution and that in some cases it is co-located with spiral arms.

79 ASTRONOMY AND ASTROPHYSICS↗

Combining In-situ Diagnostics and Data Analytics for Discovery of Process-Structure-Property Relationships in AM parts – A Step Toward Digital Twins

In-situ additive manufacturing (AM) diagnostic tools (e.g., optical/infrared imaging, acoustic, etc.) already exist to correlate process anomalies to printed part defects. This current work aimed to augment existing capabilities by: 1) Incorporating in-situ imaging w/ machine learning (ML) image processing software (ORNL- developed "Peregrine") for AM process anomaly detection 2) Synchronizing multiple in-situ sensors for simultaneous analysis of AM build events 3) Correlating in-situ AM process data, generated part defects and part mechanical properties The key R&D question investigated was to determine if these new combined hardware/software tools could be used to successfully quantify defect distributions for parts build via SNL laser powder bed fusion (LPBF) machines, aiming to better understand data-driven process-structure-property- performance relationships. High resolution optical cameras and acoustic microphones were successfully integrated in two LPBF machines and linked to the Peregrine ML software. The software was successfully calibrated on both machines and used to image hundreds of layers of multiple builds to train the ML software in identifying printed part vs powder. The software's validation accuracy to identify this aspect increased from 56% to 98.8% over three builds. Lighting conditions inside the chamber were found to significantly impact ML algorithm predictions from in-situ sensors, so these were tailored to each machine's internal framework. Finally, 3D part reconstructions were successfully generated for a build from the compressed stack of layer-wise images. Resolution differences nearest and furthest from the optical camera were discussed. Future work aims to improve optical resolution, increase process anomalies identified, and integrate more sensor modalities.

36 MATERIALS SCIENCE↗

Window Observables for Benchmarking Parton Distribution Functions

Global analysis of collider and fixed-target experimental data and calculations from lattice quantum chromodynamics (QCD) are used to gain complementary information on the structure of hadrons. We propose novel “window observables” that allow for higher precision cross-validation between the different approaches, a critical step for studies that wish to combine the datasets. Global analyses are limited by the kinematic regions accessible to experiment, particularly in a range of Bjorken-𝑥, and lattice QCD calculations also have limitations requiring extrapolations to obtain the parton distributions. We provide two different window observables that can be defined within a region of 𝑥 where extrapolations and interpolations in global analyses remain reliable and where lattice QCD results retain sensitivity and precision.

lattice QCD↗

Window observables for benchmarking parton distribution functions

Global analysis of collider and fixed-target experimental data and calculations from lattice quantum chromodynamics (QCD) are used to gain complementary information on the structure of hadrons. We propose novel ``window observables'' that allow for higher precision cross-validation between the different approaches, a critical step for studies that wish to combine the datasets. Global analyses are limited by the kinematic regions accessible to experiment, particularly in a range of Bjorken-x, and lattice QCD calculations also have limitations requiring extrapolations to obtain the parton distributions. We provide two different ``window observables'' that can be defined within a region of x where extrapolations and interpolations in global analyses remain reliable and where lattice QCD results retain sensitivity and precision.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NMR_DIFFSIM v. 1.0

The software generates and analyzes nuclear magnetic resonance spin diffusion simulations for the purpose of measuring domain sizes or spin diffusion constants in a heterogenous sample. Its primary capabilities are: 1. It solves the spin magnetization diffusion equation to simulate the results of a spin diffusion experiment using any structural input geometry, including structure models generated by molecular dynamics or other computational or experimental methods. 2. It compares the results of simulations on an input structure to experimental data, providing a data fitting algorithm to extract physical information from the experimental data via the simulation. 3. It can provide statistical measures to compare correlations between structural models, to determine domain clustering and connectedness, and to provide goodness-of-fit measures including chord length distributions, chi-squared values, Pearson coefficients and pair correlation functions. 4. It provides output reports, data files, and movies to visualize simulation results.

Sorte, Eric↗

Memory access optimization for particle operations in computational fluid dynamics-discrete element method simulations

Computational Fluid Dynamics - Discrete Element Method is used to model gas-solid systems in several applications in energy, pharmaceutical and petrochemical industries. Computational performance bottlenecks often limit the problem sizes that can be simulated at industrial scale. The data structures used to store several millions of particles in such large-scale simulations have a large memory footprint that does not fit into the processor cache hierarchies on current high-performance-computing platforms, leading to reduced computational performance. This paper specifically addresses this aspect of memory access bottlenecks in industrial scale simulations. The use of space-filling curves to improve memory access patterns is described and their impact on computational performance is quantified in both shared and distributed memory parallelization paradigms. The Morton space filling curve applied to uniform grids and k-dimensional tree partitions are used to reorder the particle data-structure thus improving spatial and temporal locality in memory. The performance impact of these techniques when applied to two benchmark problems, namely the homogeneous-cooling-system and a fluidized-bed, are presented. We report these optimization techniques lead to approximately two-fold performance improvement in particle focused operations such as neighbor-list creation and data-exchange, with ~ 1.5 times overall improvement in a fluidization simulation with 1.27 million particles.

97 MATHEMATICS AND COMPUTING↗

Systems and methods for generating a molecular dynamic graded lattice structure and their application to additive manufacturing

Systems and methods for generating molecular dynamic graded lattice structures that can be used as infill for additively manufactured articles. Molecular dynamically generated lattice infill is based on force balancing a node distribution instead of a circle packing. Field data can be utilized to adjust the spacing of the node distribution according to a force balance equilibrium model that accounts for the field expected to be experienced by the article being additively manufactured. The resultant non-uniform honeycomb structures from force-balancing robustly and efficiently address the connection issues with traditional non-uniform lattice structures.

Kim, Seokpum↗

Inverse design of two-dimensional graphene/h-BN hybrids by a regressional and conditional GAN

Design of materials with desired properties is currently laborious and heavily relies on intuition of researchers through a trial-and-error process. To tackle this challenge, in this work we propose a novel regressional and conditional generative adversarial network (RCGAN) for inverse design of representative two-dimensional materials, the graphene and boron-nitride (BN) hybrids. RCGAN incorporates a supervised regressor network, thus overcoming the common technical barrier in the traditional unsupervised GANs, which cannot generate data when fed with continuous and quantitative labels. RCGAN can autonomously generate graphene/BN hybrids given any target bandgap values. These structures are distinguished from the ones used for training and exhibit high diversity for a given bandgap. Moreover, they exhibit high fidelity, yielding bandgaps within ~10% MAE F of the desired bandgaps as validated by density functional theory (DFT) calculations. Analysis by the principle component analysis (PCA) and modified locally linear embedding (MLLE) reveals that the generator has successfully generated structures following the statistical distribution of the real structures. It implies the possibility of the RCGAN in recognizing physical rules hidden in the high-dimensional data. The novel strategy for designing regressional GAN architecture together with the successful application to inverse design of materials would inspire further exploration in research fields beyond materials.

36 MATERIALS SCIENCE↗

Analysis of TPA Pulsed-Laser-Induced Single-Event Latchup Sensitive-Area

Two-photon absorption (TPA) testing is employed to analyze the laser-induced latchup sensitive-volume (SV) of a specially designed test structure. This method takes into account the existence of an onset region in which the probability of triggering latchup transitions from zero to one as the laser pulse energy increases. This variability is attributed to pulse-to-pulse variability, uncertainty in measurement of the pulse energy, and variation in local carrier density and temperature. For each spatial position, the latchup probability associated with a given energy is calculated from multiple pulses. The latchup probability data are well-described by a Weibull distribution. The results show that the area between p-n-p-n cell structures is more sensitive than the p+ and n+ source areas, and locations far from the well contacts are more sensitive than those near the contact region. The transition from low probability of latchup to high probability is more abrupt near the source contacts than it is for the surrounding areas.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Local cation ordering in compositionally complex Ruddlesden–Popper n = 1 oxides

The Ruddlesden–Popper (RP) layered perovskite structure is of great interest due to its inherent tunability, and the emergence and growth of the compositionally complex oxide (CCO) concept endows the RP family with further possibilities. Here, a comprehensive assessment of thermodynamic stabilization, local order/disorder, and lattice distortion was performed in the first two reported examples of lanthanum-deficient La n+1 B n O 3n+1 (n = 1, B = Mg, Co, Ni, Cu, Zn) obtained via various processing conditions. Chemical short-range order (CSRO) at the B-site and the controllable excess interstitial oxygen (δ) in RP-CCOs are uncovered by neutron pair distribution function analysis. Reverse Monte Carlo analysis of the data, Metropolis Monte Carlo simulations, and extended x-ray absorption fine structure analysis implies a modest degree of magnetic element segregation on the local scale. Further, ab initio molecular dynamics simulations results obtained from special quasirandom structure disagree with experimentally observed CSRO but confirm Jahn–Teller distortion of CuO 6 octahedra. These findings highlight potential opportunities to control local order/disorder and excess interstitial oxygen in layered RP-CCOs and demonstrate a high degree of freedom for tailoring application-specific properties. They also suggest a need for expansion of theoretical and data modeling approaches in order to meet the innate challenges of CCO and related high-entropy phases.

36 MATERIALS SCIENCE↗

Near L-edge Single and Multiple Photoionization of Singly Charged Iron Ions

Absolute cross-sections for m-fold photoionization (m=1, …, 6) of Fe{sup +} by a single photon were measured employing the photon–ion merged-beams setup PIPE at the PETRA III synchrotron light source, operated by DESY in Hamburg, Germany. Photon energies were in the range 680–920 eV, which covers the photoionization resonances associated with 2p and 2s excitation to higher atomic shells as well as the thresholds for 2p and 2s ionization. The corresponding resonance positions were measured with an uncertainty of ±0.2 eV. The cross-section for Fe{sup +} photoabsorption is derived as the sum of the individually measured cross-sections for m-fold ionization. Calculations of the Fe{sup +} absorption cross-sections were carried out using two different theoretical approaches, Hartree–Fock including relativistic extensions and fully relativistic multiconfiguration Dirac–Fock. Apart from overall energy shifts of up to about 3 eV, the theoretical cross-sections are in good agreement with each other and with the experimental results. In addition, the complex de-excitation cascades after the creation of inner-shell holes in the Fe{sup +} ion were tracked on the atomic fine-structure level. The corresponding theoretical results for the product charge-state distributions are in much better agreement with the experimental data than previously published configuration-average results. The present experimental and theoretical results are valuable for opacity calculations and are expected to pave the way to a more accurate determination of the iron abundance in the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗

Understanding protein-complex assembly through grand canonical maximum entropy modeling

Inside a cell, heterotypic proteins assemble in inhomogeneous, crowded systems where the abundance of these proteins vary with cell types. While some protein complexes form putative structures that can be visualized with imaging, there are far more protein complexes that are yet to be solved because of their dynamic associations with one another. Nevertheless, it is possible to infer these protein complexes through a physical model. However, it is often not clear to physicists what kind of data from biology is necessary for such a modeling endeavor. Here, we aim to model these clusters of coarse-grained protein assemblies from multiple subunits through the constraints of interactions among the subunits and the chemical potential of each subunit. We obtained the constraints on the interactions among subunits from the known protein structures. We inferred the chemical potential that dictates the particle number distribution of each protein subunit from the knowledge of protein abundance from experimental data. Guided by the maximum entropy principle, we formulated an inverse statistical mechanical method to infer the distribution of particle numbers from the data of protein abundance as chemical potentials for a grand canonical multicomponent mixture. Using grand canonical Monte Carlo simulations, we captured a distribution of high-order clusters in a protein complex of succinate dehydrogenase with four known subunits. The complexity of hierarchical clusters varies with the relative protein abundance of each subunit in distinctive cell types such as lung, heart, and brain. When the crowding content increases, we observed that crowding stabilizes emergent clusters that do not exist in dilute conditions. We, therefore, proposed a testable hypothesis that the hierarchical complexity of protein clusters on a molecular scale is a plausible biomarker of predicting the phenotypes of a cell.

59 BASIC BIOLOGICAL SCIENCES↗

Understanding the Adsorption of Rare-Earth Elements in Oligo-Grafted Mesoporous Carbon

Rare-earth elements (REEs) are 17 elements of the periodic table primarily consisting of lanthanides. In modern society, the usage of REEs is ubiquitous in almost all modern gadgets and therefore efficient recovery and separation of REEs are of high importance. Selective adsorption and chelation of REEs in solid sorbents is a unique and sustainable process for their recovery. In this work, single-stranded oligos with 100 units of thymine were grafted onto carboxylated mesoporous carbon to synthesize a sorbent with phosphorus and oxygen functionalities. Additionally, the sorbent was characterized by X-ray photoelectron spectroscopy, Fourier-transform infrared spectroscopy, and scanning electron microscopy-energy-dispersive X-ray spectroscopy. Three different REEs with varying atomic radii and densities, Lu, Dy, and La, were adsorbed onto the carbon from aqueous solutions. It was observed that the adsorbed amounts increased with the increase in the atomic radius or decrease in the atomic density. Calculation of the distribution coefficients for all the equilibrium adsorption amounts suggested that adsorption is more effective in the lower concentration region. The L 3 -edge X-ray absorption near-edge structure confirmed a 3+ oxidation state of REEs in the adsorbed phase. Extended X-ray absorption fine structure (EXAFS) confirmed the binding of REEs with oxygen functionalities in the adsorbed phase. The radial distribution functions calculated from the EXAFS data suggest a longer RE–O distance for La compared to those for Lu and Dy. The coordination numbers and Debye–Waller factors have typical values of about 8–9 atoms and 0.01–0.02 Å 2 , respectively.

36 MATERIALS SCIENCE↗

Bayesian Adaptive Polynomial Chaos Expansions

Polynomial chaos expansions (PCEs) are widely used for uncertainty quantification (UQ) tasks, particularly in the applied mathematics community. However, PCE has received comparatively less attention in the statistics literature, and fully Bayesian formulations remain rare—especially with implementations in R. Motivated by the success of adaptive Bayesian machine learning models such as BART, BASS and BPPR, we develop a new fully Bayesian adaptive PCE method with an efficient and accessible R implementation: khaos. Our approach includes a novel proposal distribution that enables data-driven interaction selection and supports a modified g-prior tailored to PCE structure. Through simulation studies and real-world UQ applications, we demonstrate that the Bayesian adaptive PCE provides competitive performance for surrogate modeling, global sensitivity analysis and ordinal regression tasks.

97 MATHEMATICS AND COMPUTING↗

Fission fragment decay simulations with the CGMF code

The CGMF code implements the Hauser-Feshbach statistical nuclear reaction model to follow the de-excitation of fission fragments by successive emissions of prompt neutrons and γ rays. The Monte Carlo technique is used to facilitate the analysis of complex distributions and correlations among the prompt fission observables. Starting from initial configurations for the fission fragments in mass, charge, kinetic energy, excitation energy, spin, and parity, Y(A,Z,KE,U,J,π), CGMF samples neutron and γ-ray probability distributions at each stage of the decay process, conserving energy, spin and parity. Nuclear structure and reaction input data from the RIPL3 library are used to describe fission fragment properties and decay probabilities. Characteristics of prompt fission neutrons, prompt fission γ rays, and independent fission yields can be studied consistently. Correlations in energy, angle and multiplicity among the emitted neutrons and γ rays can be easily analyzed as a function of the emitting fragments. Finally, in its current version, CGMF is limited to certain spontaneous fission reactions on Pu and Cf isotopes, as well as neutron-induced fission reactions up to 20 MeV on U, Np, and Pu isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Neural chaos: A spectral stochastic neural operator

Building surrogate models for operators with uncertainty quantification capabilities is essential for many engineering applications where randomness–such as variability in material properties, boundary conditions, and initial conditions–is unavoidable. Polynomial Chaos Expansion (PCE) is widely recognized as a go-to method for constructing stochastic surrogates in both intrusive and non-intrusive ways, and it has recently been used in the context of operator learning. However, its application becomes challenging for complex or high-dimensional processes, as achieving accuracy requires higher-order polynomials, which can increase computational demand and/or the risk of overfitting. Furthermore, PCE requires specialized treatments to manage random variables that are not independent, and these treatments may be problem-dependent or may fail with increasing complexity. Here, in this work, we adopt the same formalism as the spectral expansion used in PCE; however, we replace the classical polynomial basis functions with neural network (NN) basis functions to leverage their expressivity. To achieve this, we propose an algorithm that identifies NN-parameterized basis functions in a purely data-driven manner, without any prior assumptions about the joint distribution of the random variables involved, whether independent or dependent, or about their marginal distributions. The proposed algorithm identifies each NN-parameterized basis function sequentially, ensuring they are orthogonal with respect to the data distribution. The basis functions are constructed directly on the joint stochastic variables without requiring a tensor product structure or assuming independence of the random variables. This approach may offer greater flexibility for complex stochastic models, while simplifying implementation compared to the tensor product structures typically used in PCE to handle random vectors. This is particularly advantageous given the current state of open-source packages, where building and training neural networks can be done with just a few lines of code and extensive community support. We demonstrate the effectiveness of the proposed scheme through several numerical examples of varying complexity and provide comparisons with classical PCE.

Polynomial chaos expansion↗

Analysis of the MHD stability and energetic particles effects on EIC events in LHD plasma using a Landau-closure model

The aim of this study is to perform a theoretical analysis of the magnetohydrodynamic (MHD) stability and energetic particle effects on a LHD equilibria, calculated during a discharge where energetic-ion-driven resistive interchange mode (EIC) events were triggered. Here, we use the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particles species, including the effect of the acoustic modes, multiple energetic particles (EP) species, helical couplings and helically trapped EP. We add the Landau damping and resonant destabilization effects using a closure relation. The simulations suggest that the helically trapped EP driven by the perpendicular neutral beam injector (NBI) further destabilizes the 1/1 MHD-like mode located at the plasma periphery ( r / a = 0.88). If the β of the EP driven by the perpendicular NBI is larger than 0.0025 a 1/1 EIC with a frequency around 3 kHz is destabilized. If the effect of the passing EP driven by the tangential NBI is included on the model, any enhancement of the injection intensity of the tangential NBI below β leads to a decrease of the instability growth rate. The simulations indicate that the perpendicular NBI EP is the main driver of the EIC events, as it was observed in the experiment. If the effect of the helical couplings are added in the model, an 11/13 EIC is destabilized with a frequency around 9 kHz, inward shifted ( r / a = 0.81) compared to the 1/1 EIC. Thus, one possible explanation for the EIC frequency chirping down from 9 to 3 kHz is a transition between the 11/13 to the 1/1 EIC due to a weakening of the destabilizing effect of the high n modes, caused by a decrease of the EP drive due to a loss of helically trapped EP or a change in the EP distribution function after the EIC burst. The experimental data during the EIC bursting phase shows a complex mode structure and an inward shift of the instability, although no direct evidence of the proposed transition has been observed yet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Code Analysis and Refactoring with Clang Tools, Version 0.1

Code Analysis and Refactoring with Clang Tools is a small set of example code that demonstrates techniques for applying tools distributed with the open source Clang compiler. Examples include analyzing where variables are used and replacing old data structures with standard structures.

Kelley, TimothyM.↗