Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed computation”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Computational framework for behind-the-meter DER techno-economic modeling and optimization: REopt Lite

The energy system is undergoing a major transformation with the global emphasis on decarbonization. Distributed generation is projected to play a significant role in the new energy system, and energy models are informing how distributed generation can be integrated reliably and economically. In this work, we present an end-to-end computational framework for distributed energy resource (DER) modeling, REopt Lite™, which captures the interface of technology, economics, and policy in the energy modeling process. We describe the problem space, the building blocks of the model, the scaling capabilities of the design, the optimization formulation, and the extensibility of the model. We present a framework for accelerating the techno-economic analysis of behind-the-meter distributed energy resources to enable rapid planning and decision-making, thereby enabling greater renewable energy deployment. This computation framework is open-sourced to facilitate transparency, flexibility, and wider collaboration opportunities within the worldwide energy modeling community.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Tumor Microenvironment Targeting Nano–Bio Emulsion for Synergistic Combinational X-Ray PDT with Oncolytic Bacteria Therapy

Various cancer therapies have been developed, but tumor recurrence with incomplete tumor killing and remaining tumor cells/tissues is frequent in monotherapies. Herein, a nano–bio therapeutic emulsion formulated with multifunctional nanoscintillators and anaerobic Clostridium novyi-NT spores for synergistic image-guided combinational cancer therapy is reported. MRI visible nanoscintillators (NSs) are synthesized with a NaGdF 4 :Tb,Ce@NaGdF 4 core/shell structure for an image-guided X-ray photodynamic therapy (PDT) of the normoxic peripheral tumor. An anaerobic oncolytic bacterium (C. novyi-NT) therapy is combined to treat the hypoxic central tumor tissues. Photosensitizer-coated NSs (PS-NSs) and C. novyi-NT spores are emulsified with clinically available ethiodized oil (Lipiodol) to be the nano–bio therapeutic emulsion and injected into the tumor with computed tomography image guidance. The distribution of nano–bio therapeutic emulsion, including PS-NSs and anaerobic C. novyi-NT spores in the tumor site, is confirmed by both X-ray and T1-weighted magnetic resonance imaging. Following the image-guided X-ray PDT and anaerobic C. novyi- NT combination treatment, apoptotic cell death in cancer tissues, including both peripheral and central tumor regions, is significantly higher than in the control groups. Finally, this combination therapy approach using a nano–bio therapeutic emulsion is expected to overcome the limitations of conventional cancer therapy, resulting in increased cancer-therapeutic efficacy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Miscanthus × giganteus increases soil maximum water holding capacity compared to maize

Soil ecosystem services, like the ability to store water, have been depleted after a century of conventional, annual cropping, and perennial crops offer a solution to this and other agricultural environmental issues. We assessed the impact of Miscanthus × giganteus (miscanthus), a perennial biomass crop, on soil water holding capacity and structure compared to continuous maize (Zea mays L.) at two sites in Iowa. After three growing seasons, we measured the following: (1) maximum water holding capacity (MWHC) with and without soil structure, and (2) total porosity and pore size distribution (PSD) via micro-computed tomography (microCT). Miscanthus increased MWHC by 14.7% across both sites relative to maize (p = 0.002), and we attributed this to structural changes due to the lack of a crop effect when measured on structureless soils. No significant changes were detected in soil organic matter, texture, total porosity, or PSD that could explain the increase in MWHC under miscanthus. Our findings suggest that the increases in MWHC are primarily due to structural changes rather than increases in soil organic matter or porosity (at least porosity detectable by microCT). This study highlights miscanthus' potential to enhance soil water storage and underscores the need for further investigation to clarify the mechanisms through which this biomass crop influences soil structural properties.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Ability of x‐ray computed tomography to resolve critical flaw size in laser‐based, paste stereolithography ceramic printing of alumina

Abstract Complex alumina parts were printed using vat photopolymerization (VPP), which is a stereolithography‐based additive manufacturing (AM) technique used to shape ceramic preforms, or green parts. The critical flaw size was determined using classical fracture mechanics techniques. The strength and fracture toughness were measured and compared to flaws detected in x‐ray computed tomography (XCT or CT) distributions as well as the fracture surfaces. The strength was lower compared traditionally made alumina, and that is due to layering effects, slurry defects, and printing defects. The critical flaw size from fracture mechanics was 206 µm. XCT has high enough resolution to detect the critical flaw size and much smaller features, where the average flaw size observed in CT scans was around 80–100 µm. The fracture surfaces indicate that flaws causing failure are larger than that of the critical flaw size (∼300 µm), but fracture surfaces do not show definitive features compared to traditionally made ceramics. Since XCT can observe flaws smaller than the critical flaw size, this method can be used as a screening technique.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Jet substructure from dark sector showers

We examine the robustness of collider phenomenology predictions for a dark sector scenario with QCD-like properties. Pair production of dark quarks at the LHC can result in a wide variety of signatures, depending on the details of the new physics model. A particularly challenging signal results when prompt production induces a parton shower that yields a high multiplicity of collimated dark hadrons with subsequent decays to Standard Model hadrons. The final states contain jets whose substructure encodes their non-QCD origin. This is a relatively subtle signature of strongly coupled beyond the Standard Model dynamics, and thus it is crucial that analyses incorporate systematic errors to account for the approximations that are being made when modeling the signal. We estimate theoretical uncertainties for a canonical substructure observable designed to be sensitive to the gauge structure of the underlying object, the two-point energy correlator $ {e}_2^{\left(\beta \right)} $ , by computing envelopes between resummed analytic distributions and numerical results from Pythia. We explore the separability against the QCD background as the confinement scale, number of colors, number of flavors, and dark quark masses are varied. Additionally, we investigate the uncertainties inherent to modeling dark sector hadronization. Simple estimates are provided that quantify one’s ability to distinguish these dark sector jets from the overwhelming QCD background. Such a search would benefit from theory advances to improve the predictions, and the increase in statistics using the data to be collected at the high luminosity LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Quantifying Equiaxed vs Epitaxial Solidification in Laser Melting of CMSX-4 Single Crystal Superalloy

The competition between epitaxial vs. equiaxed solidification has been investigated in CMSX-4 single crystal superalloy during laser melting as practiced in additive manufacturing. Single-track laser scans were performed on a powder-free surface of directionally solidified CMSX-4 alloy with several combinations of laser power and scanning velocity. Electron backscattered diffraction (EBSD) mapping facilitated identification of new orientations, i.e., “stray grains” that nucleated within the fusion zone along with their area fraction and spatial distribution. Using high-fidelity computational fluid dynamics simulations, both the temperature and fluid velocity fields within the melt pool were estimated. This information was combined with a nucleation model to determine locations where nucleation has the highest probability to occur in melt pools. In conformance with general experience in metals additive manufacturing, the as-solidified microstructure of the laser-melted tracks is dominated by epitaxial grain growth; nevertheless, stray grains were evident in elongated melt pools. It was found that, though a higher laser scanning velocity and lower power are generally helpful in the reduction of stray grains, the combination of a stable keyhole and minimal fluid velocity further mitigates stray grains in laser single tracks.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Parallel sorting algorithm classification: is manual instrumentation necessary?

Understanding parallel algorithms is crucial for accelerating scientific simulations on complex, distributed memory, high-performance computers. Modern algorithm classification approaches learn semantics directly from source code to differentiate between algorithms, however, accessing source code is not always possible. We can learn about parallel algorithms from observing their performance, as programs running the same algorithms and using the same hardware should exhibit similar performance characteristics. We present an approach to learn algorithm classes from parallel performance data directly in order to classify algorithms without access to the source code. We extend previous work to enable classifying parallel sorting algorithms using automatic instrumentation instead of requiring manual region annotations in the source code. In this work, we design and demonstrate a study for classification of parallel sorting algorithms using parallel performance data collected from automatic instrumentation, and evaluate the performance of our new methodology on classification. We leverage Caliper to collect the performance data, Thicket for our exploratory data analysis (EDA), and PyTorch and Scikit-learn to evaluate the effectiveness of random forests, support vector machines (SVMs), decision trees, neural networks, and logistic regressions on parallel performance data. Additionally, we study noise in parallel performance data, whether the removal of noise and pre-processing of the data is necessary to accurately classify parallel sorting algorithms, and determine the effectiveness of features created from performance data. In conclusion, we demonstrate classification accuracy for these five different models of up to 97.7% across four different parallel algorithm classes.

Algorithm Classification↗
DOE OSTI2026

Scalable Generation of High-fidelity Synthetic Population Ensembles

Used within social simulations, synthetic population ensembles enable uncertainty quantification (UQ) methods for obtaining more robust model inference and prediction. A synthetic population ensemble is a series of plausible virtual reconstructions of an area’s population at the granularity of people and residences, generated stochastically to preserve privacy of the source population survey’s respondents. In this paper, we demonstrate the production of large synthetic population ensembles for the U.S. via Oak Ridge National Laboratory’s UrbanPop framework to support modeling of high spatial resolution energy affordability metrics from nationwide social surveys in collaboration with the fusionACS project. The study involves two scenarios: creating ensembles for (1) 17 U.S. metropolitan areas in 2019 and (2) full U.S. Census Divisions in 2023, with each scenario consisting of 41 population instances (a base realization and 40 replicates). To accomplish this task at scale, we configured an integrated system within a research cloud, comprised of virtual containerizations, GPU-enhanced functionality, and orchestrated deployments of UrbanPop’s maturing Likeness Python ecosystem. Results demonstrate we maintained high-fidelity approximations of residential totals by areas of interest and the demographic characteristics of neighborhoods while reducing manual workflow burdens. Finally, we discuss plans to fine-tune and further develop our automated workflows for truly distributed job orchestration to increase computational efficiency, as well as provide an outlook for broadening applications of the ensembles.

Cluster computing↗
DOE OSTI2023

Characterization of the In Vivo Deuteration of Native Phospholipids by Mass Spectrometry Yields Guidelines for Their Regiospecific Customization

Customization of deuterated biomolecules is vital for many advanced biological experiments including neutron scattering. However, because it is challenging to control the proportion and regiospecificity of deuterium incorporation in live systems, often only two or three synthetic lipids are mixed together to form simplistic model membranes. This limits the applicability and biological accuracy of the results generated with these synthetic membranes. Despite some limited prior examination of deuterating Escherichia coli lipids in vivo, this approach has not been widely implemented. In this report an extensive mass spectrometry-based profiling of E. coli phospholipid deuteration states with several different growth media was performed, and a computational method to describe deuterium distributions with a one-number summary is introduced. The deuteration states of 36 lipid species were quantitatively profiled in 15 different growth conditions, and tandem mass spectrometry was used to reveal deuterium localization. Regressions were employed to enable the prediction of lipid deuteration for untested conditions. Small-angle neutron scattering was performed on select deuterated lipid samples, which validated the deuteration states calculated from the mass spectral data. Based on these experiments, guidelines for the design of specifically deuterated phospholipids are described. This unlocks even greater capabilities from neutron-based techniques, enabling experiments that were formerly impossible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Self-consistent modeling of tokamak edge plasma transport with lithium sources

Magnetic confinement fusion devices require effective heat and particle exhaust solutions on the divertor plates to operate sustainably, especially under reactor-relevant conditions. Liquid lithium divertors have been proposed to address two major challenges: control of excessive heat flux to plasma-facing components through vapor shielding and minimization of core plasma contamination from impurities. The National Spherical Torus Experiment-Upgrade (NSTX-U) will explore lithium as a divertor material due to its potential to meet both objectives. We present a self-consistent coupling framework between the plasma boundary transport code UEDGE and the lithium wall transport code Wall–Li to evaluate the feasibility and operational limits of lithium-based divertors. The model aims to optimize lithium sourcing levels to prevent core plasma contamination via fuel dilution while ensuring divertor protection through vapor shielding. This integrated framework, applicable to any tokamak with lithium sources, dynamically adjusts lithium sourcing based on local plasma conditions and surface temperature. The coupled model is tested using NSTX-like geometry and plasma conditions to assess its performance and reliability. Wall–Li calculates lithium fluxes from plasma-facing components, incorporating physical sputtering, thermally enhanced sputtering, and evaporation driven by surface temperature and ion flux. These fluxes are reintroduced into UEDGE as neutral lithium atoms, enabling simulation of their transport and distribution within the plasma. UEDGE computes plasma and neutral transport, surface heat flux, and iteratively feeds this information back to Wall–Li. A small time step is employed to ensure numerical stability and convergence, enabling accurate simulations over typical tokamak discharge durations. This integrated modeling approach provides a robust tool for identifying operational regimes that balance effective lithium sourcing with minimal core plasma contamination, offering critical insights for optimizing lithium-based divertor systems in current and future fusion devices.

Magnetic confinement fusion↗
DOE OSTI2022

Atomic structure and Mott nature of the insulating charge density wave phase of 1T-TaS 2

Here, using x-ray pair distribution function (PDF) analysis and computer modeling, we explore structure models for the complex charge density wave (CDW) phases of layered 1T-TaS 2 that both well capture their atomic-level features and are amenable to electronic structure calculations. The models give the most probable position of constituent atoms in terms of 3D repetitive unit cells comprising a minimum number of Ta–S layers. Structure modeling results confirm the emergence of star-of-David (SD) like clusters of Ta atoms in the high-temperature incommensurate (IC) CDW phase and show that, contrary to the suggestions of recent studies, the low-temperature commensurate (C) CDW phase expands upon cooling thus reducing lattice strain. The C-CDW phase is also found to preserve the stacking sequence of Ta–S layers found in the room temperature, nearly commensurate (NC) CDW phase to a large extent. DFT based on the PDF refined model shows that bulk C-CDW 1T-TaS 2 also preserves the insulating state of individual layers of SD clusters, favoring the Mott physics description of the metal-to-insulator (NC-CDW to C-CDW) phase transition in 1T-TaS 2 . Our work highlights the importance of using precise crystal structure models in determining the nature of electronic phases in complex materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Neoclassical toroidal viscosity torque prediction via deep learning

GPECnet is a densely connected neural network that has been trained on GPEC data, to predict the plasma stability, neoclassical toroidal viscosity (NTV) torque, and optimized 3D coil current distributions for desired NTV torque profiles. Using NTV torque, driven by non-axisymmetric field perturbations in a tokamak, can be vital in optimizing pedestal performance by controlling the rotation profile in both the core, to ensure tearing stability, and the edge, to avoid edge localized modes (ELMs). The generalized perturbed equilibrium code (GPEC) software package can be used to calculate the plasma stability to 3D perturbations and the NTV torque profile generated by applied 3D magnetic fields. These calculations, however, involve complex integrations over space and energy distributions, which takes time to compute. Initially, GPECnet has been trained solely on data representative of the quiescent H-mode (QH) scenario, in which neutral beams are often balanced and toroidal rotation is low across the plasma profile. Lastly, this work provides the foundation for active control of the rotation shear using a combination of beams and 3D fields for robust and high performance QH mode operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

New method for detecting fast neutrino flavor conversions in core-collapse supernova models with two-moment neutrino transport

Fast-pairwise neutrino oscillations potentially affect many aspects of core-collapse supernovae (CCSNe): the explosion mechanism, neutrino signals, and nucleosynthesis in the ejecta. This particular mode of collective neutrino oscillations has a deep connection to the angular structure of neutrinos in momentum space; for instance, the appearance of electron-neutrino lepton number (ELN) angular crossings in momentum space is a good indicator of the occurrence of flavor conversions. However, many multidimensional CCSN simulations are carried out with approximate neutrino transport (such as two-moment methods), which limits access to the angular distributions of neutrinos, i.e., inhibits ELNcrossing searches. In this paper, we develop a new method of searching for ELN crossing in these CCSN simulations. The required data is the zeroth and first angular moments of neutrinos and the matter profile, all of which are available in CCSN models with the two-moment method. One of the novelties of our new method is the use of a ray-tracing neutrino transport to determine ELNs in the direction of the stellar center. It is designed to compensate for shortcomings of the crossing searches with only two angular moments. We assess the capability of the method by carrying out a detailed comparison with results of full Boltzmann neutrino transport in 1D and 2D CCSN models. In this work, we find that the ray-tracing neutrino transport improves the accuracy of crossing searches; indeed, the appearance/disappearance of the crossings is accurately detected even in the region of forward-peaked angular distributions. The new method is computationally cheap and has the benefit of efficient parallelization; hence, it will be useful for ELN-crossing searches in any CCSN models that employ two-moment neutrino transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Semicomputational calculation of Bragg shift in stratified materials

The fiber Bragg grating (FBG) may be viewed as a one dimensional photonic band-gap crystal by virtue of the periodic spatial perturbation imposed on the fiber core dielectric material. Similar to media supporting Bloch waves, the engraved weak index modulation, presenting a periodic “potential” to an incoming guided mode photon of the fiber, yields useful spectral properties that have been the basis for sensing applications and emerging quantum squeezing and solitons. The response of an FBG sensor to arbitrary external stimuli represents a multiphysics problem without a known analytical solution despite the growing use of FBGs in classical and quantum sensing and metrology. Here, we study this problem by first presenting a solid mechanics model for the thermal and elastic states of a stratified material. Additionally, the model considers an embedded optical material domain that represents the Bragg grating, here in the form of an FBG. Using the output of this model, we then compute the optical modes and their temperature- and stress-induced behavior. The developed model is applicable to media of arbitrary shape and composition, including soft matter and materials with nonlinear elasticity and geometric nonlinearity. Finally, we employ the computed surface stress and temperature distributions along the grating to analytically calculate the Bragg shift, which is found to be in reasonable agreement with our experimental measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Cosmic structure strikes back: The elimination of vector-mediated nonstandard interaction models as a mechanism for sterile neutrino dark matter production

We revisit sterile neutrino production enabled by nonstandard interactions (NSIs) among active neutrinos mediated by new bosons. We focus on vector mediators, including neutrinophilic, gauged 𝐿 𝜇 −𝐿 𝜏 , and 𝐵−𝐿 realizations, that modify in-medium dispersion and scattering, thereby altering the active-sterile conversion history. Building on a novel production framework with NSI thermal potentials and collision integrals, we compute nonthermal phase-space distributions across sterile neutrino mixing and NSI parameters and map each point to an equivalent thermal warm dark matter particle mass 𝑚 th via linear theory transfer function fitting with the cosmological structure formation Boltzmann solver. This enables a direct reinterpretation of state-of-the-art structure formation limits from Milky Way satellites, strong lensing, and the Lyman-𝛼 forest. These limits, in conjunction with x-ray decay searches, as well as results from a wide variety of particle physics experiments allow for a more complete examination of these models. We find that these vector-mediated models are ruled out when the full combination of current constraints, listed above, are taken into account. NSI scalar-mediated models and models with low reheating temperatures remain viable.

cosmology↗
DOE OSTI2020

Bayesian inference of Stochastic reaction networks using Multifidelity Sequential Tempered Markov Chain Monte Carlo

Stochastic reaction network models are often used to explain and predict the dynamics of gene regulation in single cells.These models usually involve several parameters, such as the kinetic rates of chemical reactions, that are not directly measurable and must be inferred from experimental data. Bayesian inference provides a rigorous probabilistic frame-work for identifying these parameters by finding a posterior parameter distribution that captures their uncertainty.Traditional computational methods for solving inference problems such as Markov Chain Monte Carlo methods based on classical Metropolis-Hastings algorithm involve numerous serial evaluations of the likelihood function, which in turn requires expensive forward solutions of the chemical master equation (CME). We propose an alternate approach based on a multifidelity extension of the Sequential Tempered Markov Chain Monte Carlo (ST-MCMC) sampler. This algorithm is built upon Sequential Monte Carlo and solves the Bayesian inference problem by decomposing it into a sequence of efficiently solved subproblems that gradually increase both model fidelity and the influence of the observed data. We reformulate the finite state projection (FSP) algorithm, a well-known method for solving the CME, to produce a hierarchy of surrogate master equations to be used in this multifidelity scheme. To determine the appropriate fidelity, we introduce a novel information-theoretic criteria that seeks to extract the most information about the ultimate Bayesian posterior from each model in the hierarchy without inducing significant bias. This novel sampling scheme is tested with high performance computing resources using biologically relevant problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Dispatch Manager for NEML2 Constitutive Model Calculations Embedded in MOOSE

This report describes the extended capabilities of the NEML2 constitutive modeling library, including a flexible and efficient work dispatching system designed to leverage both CPU and GPU resources. This enhancement addresses one of the primary computational challenges in large-scale simulations: the ability to distribute and execute batches of material model evaluations across heterogeneous computing devices. The new dispatch system introduces a modular set of dispatcher and scheduler classes that coordinate the flow of data and execution between devices. The dispatcher is responsible for efficiently packaging work, managing device-specific memory operations, and synchronizing results. This modularity allows for extensibility, making it straightforward to integrate additional computing backends in the future. From an implementation standpoint, the dispatcher system interfaces seamlessly with NEML2's existing models. They handle device-aware tensor operations, optimize memory transfers, and support asynchronous execution when applicable. This design ensures that batches of material points can be evaluated concurrently, substantially improving throughput compared to previous single-device or serial implementations. These improvements not only enhance the raw performance of NEML2 but also improve its usability in multiscale and high-fidelity simulations, where the simultaneous evaluation of large material point batches is critical. Benchmarks included in the report demonstrate the system’s scalability, highlighting its effectiveness when leveraging modern GPU architectures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

ExaTN: Scalable GPU-Accelerated High-Performance Processing of General Tensor Networks at Exascale

We present ExaTN (Exascale Tensor Networks), a scalable GPU-accelerated C++ library which can express and process tensor networks on shared- as well as distributed-memory high-performance computing platforms, including those equipped with GPU accelerators. Specifically, ExaTN provides the ability to build, transform, and numerically evaluate tensor networks with arbitrary graph structures and complexity. It also provides algorithmic primitives for the optimization of tensor factors inside a given tensor network in order to find an extremum of a chosen tensor network functional, which is one of the key numerical procedures in quantum many-body theory and quantum-inspired machine learning. Numerical primitives exposed by ExaTN provide the foundation for composing rather complex tensor network algorithms. We enumerate multiple application domains which can benefit from the capabilities of our library, including condensed matter physics, quantum chemistry, quantum circuit simulations, as well as quantum and classical machine learning, for some of which we provide preliminary demonstrations and performance benchmarks just to emphasize a broad utility of our library.

97 MATHEMATICS AND COMPUTING↗