Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer model”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The effect of polymer stiffness on magnetization reversal of magnetorheological elastomers

Ultrasoft magnetorheological elastomers (MREs) offer convenient real-time magnetic field control of mechanical properties that provides a means to mimic mechanical cues and regulators of cells in vitro. Here, we systematically investigate the effect of polymer stiffness on magnetization reversal of MREs using a combination of magnetometry measurements and computational modeling. Poly-dimethylsiloxane-based MREs with Young’s moduli that range over two orders of magnitude were synthesized using commercial polymers Sylgard™ 527, Sylgard 184, and carbonyl iron powder. The magnetic hysteresis loops of the softer MREs exhibit a characteristic pinched loop shape with almost zero remanence and loop widening at intermediate fields that monotonically decreases with increasing polymer stiffness. A simple two-dipole model that incorporates magneto-mechanical coupling not only confirms that micrometer-scale particle motion along the applied magnetic field direction plays a defining role in the magnetic hysteresis of ultrasoft MREs but also reproduces the observed loop shapes and widening trends for MREs with varying polymer stiffnesses.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The SONATA data format for efficient description of large-scale network models

Increasing availability of comprehensive experimental datasets and of high-performance computing resources are driving rapid growth in scale, complexity, and biological realism of computational models in neuroscience. To support construction and simulation, as well as sharing of such large-scale models, a broadly applicable, flexible, and high-performance data format is necessary. To address this need, we have developed the Scalable Open Network Architecture TemplAte (SONATA) data format. It is designed for memory and computational efficiency and works across multiple platforms. The format represents neuronal circuits and simulation inputs and outputs via standardized files and provides much flexibility for adding new conventions or extensions. SONATA is used in multiple modeling and visualization tools, and we also provide reference Application Programming Interfaces and model examples to catalyze further adoption. SONATA format is free and open for the community to use and build upon with the goal of enabling efficient model building, sharing, and reproducibility.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Understanding Model Inadequacy in TRISO Nuclear Fuel Fission Products Release Models: Empirical and Mechanistic Approaches

The increasing use of tristructural isotropic (TRISO) particle fuel in both advanced and existing reactors necessitates a thorough evaluation of uncertainties and shortcomings in TRISO fission product release models. These inadequacies arise from the simplifications made in computational models compared to experimental data. Utilizing the BISON fuel performance code and experimental data from the Advanced Gas Reactor (AGR) program provides a unique chance to rigorously assess these inadequacies within a Bayesian uncertainty quantification (UQ) framework. This study contrasts the standard Bayesian framework with the Kennedy-O'Hagan (KOH) framework, which explicitly accounts for modeling inadequacies, in the context of UQ for TRISO silver release models. It examines both the traditional Arrhenius equation and a more advanced lower-length-scale (LLS)-informed model that incorporates microstructure information. The inverse UQ process applied to AGR-2 and AGR-3/4 datasets identified modeling inadequacy as the primary source of uncertainty, with experimental noise also being significant, while model parameter uncertainty was minimal. Both the Arrhenius and LLS-informed models showed similar levels of modeling inadequacy. For forward predictive UQ using the AGR-1 dataset, the KOH framework enhanced the accuracy and quality of quantified uncertainties by approximately 30% and 40%, respectively, compared to the standard Bayesian framework. This improvement was observed for both the Arrhenius and LLS-informed models. At the engineering scale, both models performed similarly, but the LLS-informed model outperformed the Arrhenius equation at the mesoscale. These findings underscore the importance of explicitly considering modeling inadequacy in the UQ process and highlight the need for ongoing refinement of physics-based models to address these shortcomings.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

High-Strength, High-Ductility, High Entropy Alloys with High-Efficiency Native Oxide Solar Absorbers for Concentrating Solar Power Systems

This EPSCoR Project has been investigating the synergy between the excellent high-temperature mechanical behavior of FeMnNiAlCr high entropy alloys (HEA) and the high solar absorptance of their native oxides for high efficiency concentrated solar thermal power (CSP) systems working at >700°C. While HEAs have attracted substantial interest in recent years, most investigations have focused on their applications as structural materials rather than functional materials. This EPSCoR project discovered that FeMnNiAlCr HEAs can potentially be applied synergistically as both a structural and functional material for high-efficiency concentrating solar thermal power (CSP) systems working at >700°C. The HEA itself would be used in high-temperature tubing to carry molten salts or supercritical CO 2 , while its surface oxide would act as a high-efficiency solar thermal absorber. With Fe and Mn being the major components in these HEAs (adding up to ~70 at.% of the alloy), these materials are much more cost-effective than the Ni-based superalloys currently being investigated for high-temperature CSP systems. Through this research, these Fe-Mn based HEAs have demonstrated yield strengths 2-3x greater than that of stainless steel at 700°C and a creep lifetime >800 h at 700ºC under a typical CSP tubing mechanical load of 35 MPa. Their Mn-rich surface oxides maintain a high optical-to-thermal conversion efficiency of ~87% under 1000x solar concentration ratio for 20 simulated day-night thermal cycles between 750ºC and room temperature. In preliminary corrosion studies, these HEAs have sustained immersion in unpurified bromide molten salts for 14 days at 750°C with <2% weight loss, in contrast to 70% weight loss from a 316 stainless steel reference. The simultaneous achievement of promising mechanical, optical, and thermochemical properties in this FeMnNiAlCr system opens the door to new applications of HEAs in solar energy harvesting. Partnerships with Ames Laboratory and Oak Ridge National Laboratory (ORNL) also advanced our understanding of the fundamental structure-property relationships through atomic scale material characterization and first-principles computational modeling. The key research results in this project can potentially be extended to other HEAs and their native oxides. In terms of applications, the proposed FeMnNiAlCr HEA/native oxide system could potentially exceed the mechanical and the optical performance of existing tubing and solar coating materials under EERE’s CSP program at lower cost, which also aligns well with the EPSCoR Science and Technology strategies of New Hampshire in boosting the deployment of renewable energy.

14 SOLAR ENERGY↗
DOE OSTI2020

Two-dimensional electronic Stark spectroscopy of a photosynthetic dimer

Stark spectroscopy, which measures changes in the linear absorption of a sample in the presence of an external DC electric field, is a powerful experimental tool for probing the existence of charge-transfer (CT) states in photosynthetic systems. CT states often have small transition dipole moments, making them insensitive to other spectroscopic methods, but are particularly sensitive to Stark spectroscopy due to their large permanent dipole moment. In a previous study, we demonstrated a new experimental method, two- dimensional electronic Stark spectroscopy (2DESS), which combines two-dimensional electronic spectroscopy (2DES) and Stark spec- troscopy. In order to understand how the presence of CT states manifest in 2DESS, here, we perform computational modeling and calculations of 2DESS as well as 2DES and Stark spectra, studying a photosynthetic dimer inspired by the photosystem II reaction center. Here, we identify specific cases where qualitatively different sets of system parameters produce similar Stark and 2DES spectra but significantly different 2DESS spectra, showing the potential for 2DESS to aid in identifying CT states and their coupling to excitonic states.

14 SOLAR ENERGY↗
DOE OSTI2022

Thermoelectric coolers for high-power-density 3D electronics heat management

Future advancements in three-dimensional (3D) electronics require robust thermal management methodology. Thermoelectric coolers (TECs) are reliable and solid-state heat pumping devices with high cooling capacity that can meet the requirements of emerging 3D microelectronic devices. In this work, we first provide the design of TECs for electronics cooling using a computational model and then experimentally validate the main predictions. Key device parameters such as device thickness, leg density, and contact resistance were studied to understand their influence on the performance of TECs. Our results show that it is possible to achieve high cooling power density through optimization of TE leg height and packing density. Scaling of TECs is shown to provide ultra-high cooling power density.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

RESOLVING THE ELECTROCHEMICAL EQUATIONS OF A SOLID OXIDE FUEL CELL FOR USE IN TRANSIENT SIMULATION AND INTEGRATION INTO CYBER-PHYSICAL SYSTEMS

A major challenge with complex cyber-physical systems stems from long model computational time that creates a mismatch between the model system and the physical system. The numerical modeling of solid oxide fuel cells (SOFCs) presents particular challenges due to the highly coupled nature of the underlying equations and the multiphysics needed to fully resolve their behavior during a transient event. To this end current approaches revolve around splitting the computational efforts into resolving temperature effects and resolving electrochemical effects. Current methods employed for the transient simulation of an SOFC for implementation in the Hybrid Performance (HyPer) facility cyberphysical plant at the National Energy Technology Laboratory reveal a distinct need for accelerated results with a high degree of stability. To this aim, an investigation into the computational time for the code reveals that the underlying electrochemical algorithm takes an order of magnitude more time than its thermal counterpart and has a tendency to vary in terms of iteration time and as such a rework of the underlying system is proposed. The primary method for accelerated electrochemical algorithm solutions is to employ higher order root finding recipes for the resolution of the highly coupled electrochemical equations. This is done with the intention to reduce the overall number of subiterations necessary for resolving voltage, current density, and species concentration, properties of the fuel cell that are all directly coupled and require nested iterative approaches. The overall objective of this approach is an order of magnitude reduction in calculation time without sacrificing stability and increasing accuracy. Specific approaches involve using both bounded and unbounded techniques, such as the False Position method and the Secant method (or if applicable Newton-Raphson) respectively, the drawbacks being slower convergence for False Position and instability for the Secant or Newton-Raphson methods. Current preliminary results on simplified versions of the parent functions involved for electrochemical calculations indicate a reduction in computational steps by a factor of two for the secant method and a factor of three for Newton-Raphson. When implemented into new modified electrochemical algorithms, the results indicate a possible order of magnitude reduction in calculation time.

Arias, Jesus↗
DOE OSTI2022

Variational multiscale reinforcement learning for discovering reduced order closure models of nonlinear spatiotemporal transport systems

Abstract A central challenge in the computational modeling and simulation of a multitude of science applications is to achieve robust and accurate closures for their coarse-grained representations due to underlying highly nonlinear multiscale interactions. These closure models are common in many nonlinear spatiotemporal systems to account for losses due to reduced order representations, including many transport phenomena in fluids. Previous data-driven closure modeling efforts have mostly focused on supervised learning approaches using high fidelity simulation data. On the other hand, reinforcement learning (RL) is a powerful yet relatively uncharted method in spatiotemporally extended systems. In this study, we put forth a modular dynamic closure modeling and discovery framework to stabilize the Galerkin projection based reduced order models that may arise in many nonlinear spatiotemporal dynamical systems with quadratic nonlinearity. However, a key element in creating a robust RL agent is to introduce a feasible reward function, which can be constituted of any difference metrics between the RL model and high fidelity simulation data. First, we introduce a multi-modal RL to discover mode-dependant closure policies that utilize the high fidelity data in rewarding our RL agent. We then formulate a variational multiscale RL (VMRL) approach to discover closure models without requiring access to the high fidelity data in designing the reward function. Specifically, our chief innovation is to leverage variational multiscale formalism to quantify the difference between modal interactions in Galerkin systems. Our results in simulating the viscous Burgers equation indicate that the proposed VMRL method leads to robust and accurate closure parameterizations, and it may potentially be used to discover scale-aware closure models for complex dynamical systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Data analysis and modeling pipelines for controlled networked social science experiments

There is large interest in networked social science experiments for understanding human behavior at-scale. Significant effort is required to perform data analytics on experimental outputs and for computational modeling of custom experiments. Moreover, experiments and modeling are often performed in a cycle, enabling iterative experimental refinement and data modeling to uncover interesting insights and to generate/refute hypotheses about social behaviors. The current practice for social analysts is to develop tailor-made computer programs and analytical scripts for experiments and modeling. This often leads to inefficiencies and duplication of effort. In this work, we propose a pipeline framework to take a significant step towards overcoming these challenges. Our contribution is to describe the design and implementation of a software system to automate many of the steps involved in analyzing social science experimental data, building models to capture the behavior of human subjects, and providing data to test hypotheses. The proposed pipeline framework consists of formal models, formal algorithms, and theoretical models as the basis for the design and implementation. We propose a formal data model, such that if an experiment can be described in terms of this model, then our pipeline software can be used to analyze data efficiently. The merits of the proposed pipeline framework is elaborated by several case studies of networked social science experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Population balance modeling of polyurethane foam formation with pressure‐dependent growth kernel

Abstract Polyurethane foams are widely used materials often chosen for their useful characteristics such as low thermal conductivity, ease of application, and high strength‐to‐weight ratios. Computational models are needed to predict the dynamics of the flow and expansion, and the resulting material properties, to improve manufacturing processes. In this paper, a model for PMDI, a water‐blown polyurethane foam, is presented. By extending a kinetics‐based approach by adding bubble‐scale information via a population balance equation (PBE) using the quadrature method of moments, we can track bubble size distributions during foaming. We present results from a three‐dimensional computational fluid dynamics model using arbitrary Lagrangian–Eulerian interface tracking implemented in finite element software. The model compares favorably with experimental data, including dynamics, bubble distributions measured by both camera and diffusion wave spectroscopy, and post‐test bubble size from scanning electron microscopy and density measurements from x‐ray computed tomography.

Ortiz, Weston↗
DOE OSTI2021

An efficient numerical model for predicting residual stress and strain in parts manufactured by laser powder bed fusion

Abstract Computational modeling of additively manufactured structures plays an increasingly important role in product design and optimization. For laser powder bed fusion processes, the accurate modeling of stress and distortion requires large amount of computational cost due to very localized heat input and evolving complex geometries. The current study takes advantage of a graphics processing unit accelerated explicit finite element analysis code and approximated heat conduction analysis to predict the macroscopic thermo-mechanical behavior in laser selective melting. Adjacent layers and tracks were lumped to reduce the number of time steps and elements in the finite element model. The effects of track and layer grouping on prediction accuracy and solution efficiency are investigated to provide a guidance for a cost-effective simulation. Thin-wall builds from Inconel alloy 625 (IN625) powders were simulated by applying the developed modeling approach to get the detailed residual stress and distortion at a computational speed 50 times higher than conventional approach. Under repeated heating and cooling cycles, a high tensile stress was produced near surfaces of a build due to a larger shrinkage on surface than that in central area. It is also shown that horizontal stresses concentrate near the root and top layers of the IN625 build. The predicted residual elastic strain distribution was validated by the experimental measurement using x-ray synchrotron diffraction.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Atomic layer deposited aluminum nitride for interface protection of lithium germanium thiophosphate solid electrolytes

Lithium germanium thiophosphate (LGPS) is an attractive solid-state electrolyte material due to its exceptionally high ionic conductivity (~1.2 × 10 –2 S cm –1 ), comparable to many organic liquid electrolytes commonly used in batteries. Despite the high conductivity of LGPS, the susceptibility of LGPS to deleterious degradation reactions has impeded its commercial adoption into solid-state batteries. In particular, the poor voltage stability of LGPS with high-voltage cathode or lithium metal potentials often results in dramatically increasing cell impedance during galvanic cycling. Here, we use aluminum nitride (AlN) as a protection layer for Li metal anode, applied directly to the LGPS at 250 °C using plasma-enhanced atomic layer deposition, to enhance cell performance by preventing LGPS-Li metal reactions. We compare the surface chemistry and electrochemical cycling performance of atomic layer-deposited AlN grown using both plasma N 2 and NH 3 precursors. Galvanostatic cycling and electrochemical impedance spectroscopy show that AlN-coated LGPS cells perform better than bare LGPS cells in contact with Li metal anodes, with the AlN able to improve cycling longevity by over a factor of 3 in certain cases. Finally, we utilize x-ray photoelectron spectroscopy (XPS) line scans to highlight the slow room-temperature reactivity between AlN and evaporated lithium metal, and a computational model is built to aid further XPS analysis.

25 ENERGY STORAGE↗
DOE OSTI2025

Establishing a process-structure-property-performance framework for SLS additive manufacturing through integrated multiscale modeling

This study presents a comprehensive suite of high-fidelity computational models that integrate multiscale and multiphysics simulations to capture the full Selective Laser Sintering (SLS) additive manufacturing process—from initial melting and solidification to mechanical response under external loads. Process simulations are linked with mechanical analysis through Representative Volume Elements (RVEs), establishing a process-structure–property-performance framework. The interaction between laser light and polyamide 12 (PA12) powder is modeled, accounting for laser characteristics and the optical, thermal, and geometrical properties of the powder. The heat source is incorporated into a heat transfer model, coupled with crystallization kinetics and densification models to predict material density and crystallinity. The porosity distribution from the densification model and crystallinity interpolated from experimental data are used to construct the RVEs. A multi-mechanism constitutive model is then calibrated using mechanical tests to predict the stress–strain response. Simulation results show good agreement with experimental data in terms of porosity, crystallinity, and mechanical performance when sufficient laser power (62 W or higher) is used. This research supports the inverse design of 3D-printed structures by introducing a high-fidelity framework that combines multiscale and multiphysics modeling with experimental calibration for predictive and performance-driven additive manufacturing.

SLS↗
DOE OSTI2026

Bryan Mound Geomechanical Model Update

This report summarizes the recent updates of Sandia's computational model for geomechanical analysis of salt caverns at the Bryan Mound site, one of the U.S. Strategic Petroleum Reserves.

02 PETROLEUM↗
DOE OSTI2023

High-Temperature Air Stability of Electrodeposited Copper Cobalt Oxide and Copper Manganese Oxide Absorber Coatings for Concentrating Solar Power

The solar absorptance of a thermal receiver surface significantly affects the photothermal conversion efficiency of concentrating solar power (CSP) plants. The stability of solar receivers and coatings applied to improve solar absorptance at high temperatures is critical for efficient and long-term operation. This study analyzes the high-temperature air stability of multiscale fractal textured electrodeposited copper manganese oxide and copper cobalt oxide solar receiver coatings over a range of temperatures. The study also investigates the effect of high temperature on phase stability, chemical and morphological evolution, and the optical properties of the two solar absorber materials. A computational model is developed to determine the absorptance of the textured absorbers and is shown to match the experimental measurements closely over a range of heat treatment temperatures from 450 °C to 850 °C and durations up to 100 h. Degradation curves and coating stability maps are developed from the study. Copper manganese oxide and copper cobalt oxide coatings exhibit good thermal stability up to ~815 °C and 630 °C, respectively, with minimal changes in their optical and structural properties with thermal efficiency greater than 90%. Furthermore, the study offers valuable insights into the high-temperature annealing effects on the optical and structural properties of multiscale copper alloyed coatings fabricated by electrodeposition and presents a means of designing the optimal temperature range for maximizing absorber efficiency in practical applications.

14 SOLAR ENERGY↗
DOE OSTI2022

Non-equilibrium plasma generation via nano-second multi-mode laser pulses

The formation and growth of plasma kernels generated via nano-second mode-beating laser pulses is investigated here via a non-equilibrium self-consistent computational model. Chemically reactive Navier–Stokes equations are used to describe the hydrodynamics, and non-equilibrium effects are taken into account with a two-temperature model. Inverse Bremsstrahlung and multiphoton ionization are included self-consistently in the model via a coupled solution of the plasma governing equations and the radiative transfer equation (that describes the laser beam propagation and attenuation). A self-consistent approach (despite carrying additional challenges) minimizes the empiricism and it allows for a more accurate description since it prevents both the utilization of artificial plasma seeds to trigger the breakdown and the implementation of tuning parameters to simulate the laser-energy deposition. The advantages of this approach are confirmed by the good agreement between the numerically predicted and the experimentally measured plasma boundary evolution and absorbed energy. This also holds true for the periodic plasma kernel structures that, as suggested by the experiments and confirmed by the simulations presented here, are connected to the modulating frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Hybrid (PDE+ML) models in the context of land ice modeling

Focal Area(s): Primary focal areas are: predictive modeling through the use of AI-derived model components; advanced methods including network design/optimization/deep learning. Science Challenge: The atmospheric, ocean and ice dynamics components of the Energy Exascale Earth System Model (E3SM) are governed by Partial Differential Equations (PDEs) and significant efforts have been made during the last decades to develop such computational models. Here we propose to fundamentally improve these PDE-based codes by enhancing them with Machine Learning (ML) sub-models for complex, poorly understood physical processes in the context of ice sheet modeling. We propose to train these models with a novel approach that allows the assimilation of the different sources of data available (direct/indirect observations and possibly simulation data), improving on existing simplified models. We also highlight computational challenges originating from the coexistence of PDE-based and ML-based models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Multi-scale dynamics of atmospheric-pressure discharges ignited over liquid electrodes

A multiphase computational model of atmospheric-pressure gas discharges ignited over liquid electrodes is developed. The model takes into account both the liquid electrode motion under applied electric fields and plasma generation in the gas phase. The influence of the applied voltage and liquid properties (density, viscosity, and surface tension) on the liquid surface deformation is analyzed. It is shown that the shape of dynamic cones formed on the liquid surface differs from the static Taylor's cones. The influence of the liquid surface protrusions on gas breakdown dynamics is demonstrated. It is shown that the breakdown develops in two stages: first, a fast ionization wave propagates from the cathode to the anode, and then the cathode sheath collapses due to secondary electron emission from the liquid surface. It is shown that dynamics of the multiphase system containing gas, plasma, and liquid states is characterized by three disparate time scales: the fast electron time scale, the intermediate ion time scale, and the slow liquid dynamics time scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗