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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The Impact of Radiative Transfer at Reduced Spectral Resolution in Large‐Eddy Simulations of Convective Clouds

Abstract Many radiative transfer schemes approximate the spectral integration over ∼10 5 to ∼10 6 wavelengths with correlated k ‐distributions methods that typically require only 10 1 –10 2 spectral integration points ( g ‐points). The exact number of g ‐points is then chosen as an optimal balance between computational costs and accuracy, normally assessed in terms of a number of radiative quantities. How this radiative accuracy propagates to simulation accuracy, however, is not straightforward. In this study, we therefore explore the sensitivity of cloud properties in large‐eddy simulations (LES) to the accuracy of radiative fluxes and heating rates. We first generate smaller sets of g ‐points from existing k ‐distributions by repeatedly combining adjacent g ‐points while maintaining the highest possible accuracy on a chosen set of radiative metrics. Next, we perform three sets of LES with varying cloud—radiation coupling pathways, and therefore different requirements for the accuracy of the radiative transfer computations, to investigate how these smaller and thus less accurate k ‐distributions affect simulation characteristics. The decrease in radiative accuracy with 3–4 times smaller k ‐distributions results in biases in cloud properties that are relative small compared to their temporal fluctuations. These results show potential for speeding up radiative transfer computations in cloud‐resolving models by reducing the resolved spectral detail. However, more statistically converged simulations and a wider set of case studies is required to fully assess the robustness of our results.

Meteorology & Atmospheric Sciences↗
DOE OSTI2020

Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented here in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N -body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Modeling strength and failure variability due to porosity in additively manufactured metals

To model and quantify the variability in plasticity and failure of additively manufactured metals due to imperfections in their microstructure, we have developed uncertainty quantification methodology based on pseudo marginal likelihood and embedded variability techniques. We account for both the porosity resolvable in computed tomography scans of the initial material and the sub-threshold distribution of voids through a physically motivated model. Calibration of the model indicates that the sub-threshold population of defects dominates the yield and failure response. Finally, the technique also allows us to quantify the distribution of material parameters connected to microstructural variability created by the manufacturing process, and, thereby, make assessments of material quality and process control.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bayesian operator inference for data-driven reduced-order modeling

This work proposes a Bayesian inference method for the reduced-order modeling of time-dependent systems. Informed by the structure of the governing equations, the task of learning a reduced-order model from data is posed as a Bayesian inverse problem with Gaussian prior and likelihood. The resulting posterior distribution characterizes the operators defining the reduced-order model, hence the predictions subsequently issued by the reduced-order model are endowed with uncertainty. The statistical moments of these predictions are estimated via a Monte Carlo sampling of the posterior distribution. Since the reduced models are fast to solve, this sampling is computationally efficient. Furthermore, the proposed Bayesian framework provides a statistical interpretation of the regularization term that is present in the deterministic operator inference problem, and the empirical Bayes approach of maximum marginal likelihood suggests a selection algorithm for the regularization hyperparameters. The proposed method is demonstrated on two examples: the compressible Euler equations with noise-corrupted observations, and a single-injector combustion process.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Parametric reduced order models for graded lattice structures

Graded lattice structures, characterized by smoothly varying mechanical properties, hold significant promise for optimizing material distribution in advanced engineering applications. However, accurately modeling these structures poses substantial computational challenges due to the continuous geometric variations within their unit cells. Here, to address these challenges, this paper introduces a novel Efficient Reduced Order Model (EROM) that integrates the Matrix Discrete Empirical Interpolation Method (MDEIM) and Discrete Empirical Interpolation Method (DEIM) with polynomial regression to efficiently manage geometric parametrization in lattice structures. Unlike traditional reduced order models (ROMs) that require extensive precomputed libraries for each geometric configuration, our approach enables continuous geometric variations through a flexible algebraic formulation, significantly reducing computational costs while preserving high accuracy. The method constructs projection matrices for individual unit cells that can be efficiently assembled into global systems, leveraging the repetitive nature of lattice structures. Numerical studies demonstrate that our EROM achieves displacement errors below 1% and von Mises stress prediction errors below 4%, coupled with computational speedups exceeding two orders of magnitude compared to full-order simulations. The proposed method's modularity and scalability make it particularly suitable for design optimization and real-time simulation of functionally graded lattice structures, with applications spanning aerospace to biomedical engineering.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Regen: An object layout regenerator on large-scale production HPC systems

This article proposes an object layout regenerator called Regen which regenerates and removes the object layout dynamically to improve the read performance of applications. Regen first detects frequent access patterns from the I/O requests of the applications. Second, Regen reorganizes the objects and regenerates or preallocates new object layouts according to the identified access patterns. Finally, Regen removes or reuses the obsolete or regenerated object layouts as necessary. As a result, Regen accelerates access to objects by providing a flexible object layout. We implement Regen as a framework on top of Proactive Data Container (PDC) and evaluate it on Cori supercomputer, a production-scale HPC system, by using realistic HPC I/O benchmarks. The experimental results show that Regen improves the I/O performance by up to 16.92 × compared with an existing system.

Distributed file system↗
DOE OSTI2023

Fast approximation to supershell partition functions: Explicit forms of the coefficients

In a previous work, we derived a formula for supershell partition functions, which are the cornerstone of the Super-Transition-Array approach to radiative-opacity calculations. The new expression takes the form of a functional of the distribution of energies within the supershell and allows for fast and accurate computations, truncating the number of terms in the expansion. The latter involves coefficients (denoted Γ $k$ ) for which we obtained a recursion relation. In the present short paper, as a complement of the previous one, we give an explicit formula for the coefficients Γ $k$ . Another recursion relation is also provided, as well as an alternative expression involving Bell polynomials. In conclusion, the connections with cycle indexes of permutation groups and with elementary symmetric polynomials are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Validation of two-layer model for underexpanded hydrogen jets

Previous studies have shown that the two-layer model more accurately predicts hydrogen dispersion than the conventional notional nozzle models without significantly increasing the computational expense. However, the model was only validated for predicting the concentration distribution and has not been adequately validated for predicting the velocity distributions. In the present study, particle imaging velocimetry (PIV) was used to measure the velocity field of an underexpanded hydrogen jet released at 10 bar from a 1.5 mm diameter orifice. The two-layer model was the used to calculate the inlet conditions for a two-dimensional axisymmetric CFD model to simulate the hydrogen jet downstream of the Mach disk. The predicted velocity spreading and centerline decay rates agreed well with the PIV measurements. The predicted concentration distribution was consistent with data from previous planar Rayleigh scattering measurements used to verify the concentration distribution predictions in an earlier study. The jet spreading was also simulated using several widely used notional nozzle models combined with the integral plume model for comparison. These results show that the velocity and concentration distributions are both better predicted by the two-layer model than the notional nozzle models to complement previous studies verifying only the predicted concentration profiles. Thus, this study shows that the two-layer model can accurately predict the jet velocity distributions as well as the concentration distributions as verified earlier. Though more validation studies are needed to improve confidence in the model and increase the range of validity, the present work indicates that the two-layer model is a promising tool for fast, accurate predictions of the flow fields of underexpanded hydrogen jets.

08 HYDROGEN↗
DOE OSTI2025

Cross-Scale Catalyst Modeling Applied to H 2 Storage and Release via Formic Acid

Here, we propose the Systems-to-Atoms (S2A) modeling framework that integrates the kinetics of reaction chemistry and structural configurations across various length scales with the aim of establishing a versatile template for multiscale modeling of reactive flow problems and to predict the operando activity of catalyst materials. The approach encompasses a microkinetic model to analyze surface reactions on individual facets of catalyst nanoparticles coupled with the computation of average surface reaction rates for catalyst nanoparticles of specific size distributions. Macro-homogeneous surface reaction kinetics are derived as a function of catalyst loading and used as input parameters for the continuum-scale reactor model. The cross-scale framework enables the optimization of catalyst utilization through reactor design and operating strategy. To demonstrate the framework, we studied the storage and release of hydrogen from formic acid, a promising liquid organic hydrogen carrier (LOHC), over Pd, Pt, and Cu catalysts. The framework predicts observed trends in formic acid dehydrogenation activity for catalysts with comparable weight loadings and metal particle diameters, demonstrating satisfactory quantitative alignment. Finally, the seamless transmission of parameter uncertainties between scales is also discussed.

08 HYDROGEN↗
DOE OSTI2024

Oxygen Vacancies Alter Methanol Oxidation Pathways on NiOOH

A thorough comprehension of the mechanism underlying the methanol oxidation reaction (MOR) on Ni-based catalysts is critical for future electrocatalytic design and development. However, the mechanism of MOR on these materials remains a matter of controversy. Herein, we combine in situ surface-enhanced infrared absorption spectroscopy (SEIRAS) and density functional theory (DFT) calculations to identify the active sites and determine the mechanism of MOR on monometallic Ni-based catalysts in alkaline media. The SEIRAS results show that formate and (bi)carbonate are formed after the commencement of the MOR with potential-dependent relative distributions. These spectroscopic results are in good agreement with the DFT-computed reaction profiles over an oxygen vacancy, suggesting that the MOR mainly proceeds through the formate-involving pathway, in which the early consumption of methanol yields formate as the major product, while increasing potential drives further oxidation of formate to (bi)carbonate. We also find a parallel pathway for the generation of (bi)carbonate at high potentials that bypasses the formation of formate. The two main pathways are thermodynamically more feasible than the one predominantly reported in the literature for MOR on NiOOH that involves CHO and/or CO as key intermediates. These DFT results are supported by spectroscopic evidence showing that no band associated with CHO or CO can be detected by SEIRAS, which is attributed to the nature of the oxygen vacancies as the active sites, suppressing deep dehydrogenation of CH 2 O to CHO. Furthermore, this work thus shows the promising role of defect engineering in promoting the electrocatalytic MOR activity and selectivity.

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

A high-throughput and low-waste viability assay for microbes

Abstract Counting viable cells is a universal practice in microbiology. The colony-forming unit (CFU) assay has remained the gold standard to measure viability across disciplines, but it is time-intensive and resource-consuming. Here we describe the geometric viability assay (GVA) that replicates CFU measurements over 6 orders of magnitude while reducing over 10-fold the time and consumables required. GVA computes a sample’s viable cell count on the basis of the distribution of embedded colonies growing inside a pipette tip. GVA is compatible with Gram-positive and Gram-negative planktonic bacteria ( Escherichia coli, Pseudomonas aeruginosa and Bacillus subtilis ), biofilms and fungi ( Saccharomyces cerevisiae ). Laborious CFU experiments such as checkerboard assays, treatment time-courses and drug screens against slow-growing cells are simplified by GVA. The ease and low cost of GVA evinces that it can replace existing viability assays and enable viability measurements at previously impractical scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Physics informed deep neural network embedded in a chemical transport model for the Amazon rainforest

Secondary organic aerosols (SOA) are fine particles in the atmosphere, which interact with clouds, radiation and affect the Earth’s energy budget. SOA formation involves chemistry in gas phase, aqueous aerosols, and clouds. Simulating these chemical processes involve solving a stiff set of differential equations, which are computationally expensive steps for three-dimensional chemical transport models. Deep neural networks (DNNs) are universal function approximators that could be used to represent the complex nonlinear changes in aerosol physical and chemical processes; however, key challenges such as generalizability to extended time periods, preservation of mass balance, simulating sparse model outputs, and maintaining physical constraints have limited their use in atmospheric chemistry. Here, we develop an approach of using a physics-informed DNN that overcomes previous such challenges and demonstrates its applicability for the chemical formation processes of isoprene epoxydiol SOA (IEPOX-SOA) over the Amazon rainforest. The DNN is trained with data generated by simulating IEPOX-SOA over the entire atmospheric column, using the Weather Research and Forecasting Model coupled with Chemistry (WRF-Chem). The trained DNN is then embedded within WRF-Chem to replace the computationally expensive default solver of IEPOX-SOA formation. The trained DNN predictions generalizes well with the default model simulation of the IEPOX-SOA mass concentrations and its size distribution (20 size bins) over several days of simulations in both dry and wet seasons. The embedded DNN reduces the computational expense of WRF-Chem by a factor of 2. Our approach shows promise in terms of application to other computationally expensive chemistry solvers in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Interplay of molecular dynamics and radiative decay of a TADF emitter in a glass-forming liquid

Herein we investigate the role of molecular dynamics in the luminescent properties of a prototypical thermally activated delayed fluorescence (TADF) emitter, NAI-DMAC, in solution using a combination of temperature dependent time-resolved photoluminescence and absorption spectroscopies. We use a glass forming liquid, 2-methylfuran, to introduce an abrupt change in the temperature dependent diffusion dynamics of the solvent and examine the influence this has on the emission intensity of NAI-DMAC molecules. Comparison of experiment with first principles molecular dynamics simulations reveals that the emission intensity of NAI-DMAC molecules follows the temperature-dependent self-diffusion dynamics of the solvent. A marked reduction of emission intensity is observed as the temperature decreases toward the glass transition because the rate at which NAI-DMAC molecules can access emissive molecular conformations is greatly reduced. Below the glass transition, the diffusion dynamics of the solvent changes more slowly with temperature, which causes the emission intensity to decrease more slowly as well. The combination of experiment and computation suggests a pathway by which TADF emitters may transiently access a distribution of conformational states and avoid the need for an average conformation that strikes a balance between lower singlet–triplet energy splittings versus higher emission probabilities.

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

Efficiently predicting pressure-composition-temperature diagrams to discover low-stability metal hydrides

Quantitatively accurate computational predictions of metal hydride thermodynamics are challenging but critical for alloy performance optimization across a multitude of technological domains, including hydrogen storage, compression, purification, and getters. Recent machine learning approaches have demonstrated great success in this area, but can potentially suffer from several shortcomings since they rely on imbalanced experimental training data and can have poor out-of-distribution (ood) test performance. Here, in this study, we circumvent such pitfalls by developing a computationally efficient, first principles-based workflow for direct prediction of metal hydride phase equilibrium, i.e., the pressure-composition-temperature (PCT) diagram. We then demonstrate its utility on predicting low stability hydrides derived from compositionally complex C14 Laves phase AB2 alloys. Specifically, we computationally predict and then experimentally validate an AB 2 alloy series (z < 0.6 for Ti 2−z Zr z CrMnFeNi) with ideal hydriding thermodynamics for a two-stage metal hydride-based compressor for pressurizing boil off from liquefied hydrogen. Importantly, this study lays the groundwork for accurate and efficient discovery/optimization of ood, low-stability hydrides for which purely data-driven approaches lack sufficient accuracy.

08 HYDROGEN↗
DOE OSTI2021

Explaining an unusual electromigration behavior—A comprehensive experimental and theoretical analysis using finite element method

In metallic interconnects, it is generally assumed that electromigration (EM) failure location is independent of the applied electrical current and always occurs at the highest-current-density area. Our experiments show otherwise. We designed an Al interconnect that alters its failure location by only varying the applied current density. The failure occurs near the high for a current above 2 × 10 7 A/cm 2 , but at a location with 59% of the maximum for lower current densities. Thermoreflectance thermal imaging is employed to gather time-dependent high-resolution spatial temperature distributions of the Al interconnect during EM. More importantly, we propose a computationally inexpensive 2D finite element method that tracks EM evolution in time and matches well with the observations from different experimental conditions. A detailed analysis covering the major driving forces of EM is carried out to understand the complex physics behind EM. The atomic depletion rate contributed by each force is quantitatively studied. By examining the results from every tested experimental condition, the model reveals that the temperature gradient is the key reason causing atomic depletion near the failure location. Graphical illustrations and qualitative analysis are provided to intuitively show the key findings of our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Formation and optical properties of indium nanoparticle arrays for deep-UV plasmonics

We utilize a combined computational-experimental approach to examine the influence of indium nanoparticle (NP) array distributions on deep-ultraviolet (UV) plasmon resonances. For photon energies < 5.7 eV, analysis of ellipsometric spectra reveals an increase in silicon reflectance induced by indium NP arrays on silicon. For various energies in the range 5.7–7.0 eV, a decrease in reflectance is induced by the NP arrays. Similar trends in reflectance are predicted from finite-difference time-domain (FDTD) simulations using NP size distributions extracted from atomic-force micrographs as input. In addition, in the energy range of 7.4–9.2 eV, the FDTD simulations reveal reflectance minima, characteristic of localized surface plasmon resonances. Here, electron energy-loss spectroscopy collected from individual indium NPs reveals the presence of LSPR at ≈ 8 eV, further supporting the promise of indium NP arrays on silicon for deep-UV plasmonics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Inclusive breakup cross sections in reactions induced by the nuclides 6 He and 6,7 Li in the two-body cluster model

Calculations of inclusive α particle production are performed to interpret reaction experimental data induced by the two-neutron halo 6He and the isotopes 6,7 Li on different target nuclei. In this work, we have implemented the zero-range post-form DWBA to compute the elastic and nonelastic breakup cross sections. Integrated cross sections, angular distributions and spectra are presented. Projectiles are approximated as two-body clusters, namely 6 He as α+dineutron and 6 Li ( 7 Li) are interpreted as α+deuteron(triton). The São Paulo optical potential is employed in the calculation of the distorted wave functions for the incident channels, while standard phenomenological interactions are taken for the fragment-target interactions. Calculations for 6 Li and 7 Li furnish a good description of the data, while the reaction involving 6 He is found to be more complicated, where interpretation of the data must still be considered incomplete. The present analysis identifies an overall large contribution from inclusive breakup emission for all the cases studied.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗