Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational modeling and simulations”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Benchmarking of massively parallel phase-field codes for directional solidification

We present a detailed benchmark comparing two state-of-the-art phase-field implementations for simulating alloy solidification under experimentally relevant conditions. The study investigates the directional solidification of Al-3wt%Cu under high-velocity solidification conditions and SCN-0.46wt% camphor under microgravity conditions from National Aeronautics and Space Administration (NASA) DECLIC-DSI-R experiments. Both codes, one employing finite-difference discretization with uniform mesh and GPU-acceleration (GPU-PF) and the other one employing finite-element discretization with adaptive-mesh and CPU-parallelization (PRISMS-PF), solve the same quantitative phase-field formulation that incorporates an anti-trapping current for the solidification of dilute alloys. We evaluate the predictions of each code for dendritic morphology, primary spacing, and tip dynamics in both 2D and 3D, as well as their numerical convergence and computational performance. While existing benchmark problems have primarily focused on simplified or small-scale simulations, they do not reflect the computational and modeling challenges posed by employing experimentally relevant time and length scales. Our results provide a practical framework for assessing phase-field code performance as well as validating and facilitating their application in integrated computational materials engineering (ICME) workflows that require integration with realistic experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling and mitigating airborne pathogen risk factors in school buses

Computational fluid dynamics (CFD) models were developed to simulate the impact of different ventilation scenarios on airborne exposure risks in a 72-passenger school bus. Here, scenarios and factors that were investigated included a moving vs. stationary bus, impacts of a heating unit within the bus, and impacts of alternative ventilation scenarios with different combinations of openings (e.g., windows, door, emergency hatch). Results of the simulations showed that when the bus was stationary, use of the heater increased receptor concentrations unless there was another opening. When the bus was moving, simulations with at least two sets of openings separated from each other in the forward and aft directions produced a through-flow condition that reduced concentrations via dilution from outside air by a factor of ten or more. A single opening in a moving bus generally increased concentrations throughout the cabin due to increased mixing with minimal ventilation. The cumulative exposure risk (time-averaged concentrations) was found to be inversely correlated to the air exchange rate. Stationary and moving-bus scenarios that yielded above ~20 air changes per hour resulted in the lowest cumulative exposures. Recommendations from this study were implemented in new safety and operating procedures by the Albuquerque Public Schools Transportation Center.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Optimal CO 2 storage management considering safety constraints in multi-stakeholder multi-site GCS projects: A Markov game perspective

Geological carbon storage (GCS) projects could involve a diverse array of stakeholders or players from public, private, and regulatory sectors, each with different objectives and responsibilities. Given the complexity, scale, and long-term nature of GCS operations, determining whether individual stakeholders can independently optimize their interests — or whether collaborative coalition agreements are needed — remains a central question for effective GCS project planning and management. To access large, high-quality storage resources, future GCS deployment may increasingly occur in geologically connected sites, where shared geological features such as pressure space and reservoir pore capacity can lead to competitive behavior among stakeholders. In this work, we propose a paradigm based on Markov games to quantitatively investigate how different coalition structures affect the goals of stakeholders. We frame this multi-stakeholder multi-site problem as a multi-agent reinforcement learning problem with safety constraints. Our approach enables agents to learn optimal strategies while complying with safety regulations. We present an example where multiple operators are injecting CO 2 into their respective project areas in a geologically connected basin. To address the high computational cost of repeated simulations of high fidelity models, a previously developed surrogate model based on the Embed-to-Control (E2C) framework is employed. Our results demonstrate the effectiveness of the proposed framework in addressing optimal management of CO 2 storage when multiple stakeholders with different objectives and goals are involved.

58 GEOSCIENCES↗
DOE OSTI2023

Enabling rapid X-ray CT characterisation for additive manufacturing using CAD models and deep learning-based reconstruction

Metal additive manufacturing (AM) offers flexibility and cost-effectiveness for printing complex parts but is limited to few alloys. Qualifying new alloys requires process parameter optimisation to produce consistent, high-quality components. High-resolution X-ray computed tomography (XCT) has not been effective for this task due to artifacts, slow scan speed, and costs. We propose a deep learning-based approach for rapid XCT acquisition and reconstruction of metal AM parts, leveraging computer-aided design models and physics-based simulations of nonlinear interactions between X-ray radiation and metals. This significantly reduces beam hardening and common XCT artifacts. We demonstrate high-throughput characterisation of over a hundred AlCe alloy components, quantifying improvements in characterisation time and quality compared to high-resolution microscopy and pycnometry. Our approach facilitates investigating the impact of process parameters and their geometry dependence in metal AM.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modelling Airborne Transmission and Ventilation Impacts of a COVID-19 Outbreak in a Restaurant in Guangzhou, China

Computational fluid dynamics (CFD) modelling was performed to simulate spatial and temporal airborne pathogen concentrations during an observed COVID-19 outbreak in a restaurant in Guangzhou, China. The reported seating configuration, overlap durations, room ventilation, layout, and dimensions were modelled in the CFD simulations to determine relative exposures and probabilities of infection. Results showed that the trends in the simulated probabilities of infection were consistent with the observed rates of infection at each of the tables surrounding the index patient. Alternative configurations that investigated different boundary conditions and ventilation conditions were also simulated. Increasing the fresh-air percentage to 10%, 50%, and 100% of the supply air reduced the accumulated pathogen mass in the room by an average of ~30%, ~70%, and ~80%, respectively, over 73 min. Overall, the probability of infection was reduced by ~10%, 40%, and 50%, respectively.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Design Overview of a High-Pressure Helium Flow Visualization Apparatus for Blanket Cooling Studies

Cooling of the fusion blanket first wall remains a significant challenge given the adverse conditions of heat and particle flux encountered near the plasma. Helium emerges as an attractive cooling candidate because of its chemical and neutronic inertness and separability from hydrogenic species (e.g. tritium). Because of the low thermal mass of helium, optimization of these coolant channels is warranted to provide high heat transfer performance at low pumping costs. Increasingly, computational fluid dynamics (CFD) simulations are employed to model and optimize these flow channels, and accompanying experimental data are needed to validate the predictions of these models. To provide the aforementioned experimental data, a high-pressure helium flow visualization upgrade has been designed for the Helium Flow Loop Experiment facility. This apparatus was built to American Society of Mechanical Engineers boiler and pressure vessel standards to withstand operating pressure of 4 MPa and mated to high-pressure glass windows. Seedless flow visualization is performed via high-speed background oriented schlieren (BOS), with image correlation used for time-resolved two-dimensional velocimetry at frequencies in excess of 60 kHz. Rectangular flow channel test articles are additively manufactured via laser powder bed fusion and installed into this visualization apparatus, with one-sided heating supplied by resistive heaters. In conclusion, the chosen test geometries were informed by prior CFD simulations, and the helium flow structures observed via BOS (detachment, recirculation, etc.) will be used for the validation of these accompanying models, in support of the design and optimization of blanket cooling channel configurations.

Helium flow↗
DOE OSTI2023

Disentangling the effects of non-adiabatic interactions upon ion self-diffusion within warm dense hydrogen

Warm dense matter is a material state in the region of parameter space connecting condensed matter to classical plasma physics. In this intermediate regime, we investigate the significance of non-adiabatic electron-ion interactions upon ion dynamics. To disentangle non-adiabatic from adiabatic electron-ion interactions, we compare the ion self-diffusion coefficient from the non-adiabatic electron force field computational model with an adiabatic, classical molecular dynamics simulation. A classical pair potential developed through a force-matching algorithm ensures the only difference between the models is due to the electronic inertia. Here, we implement this new method to characterize non-adiabatic effects on the self-diffusion of warm dense hydrogen over a wide range of temperatures and densities. Ultimately we show that the impact of non-adiabatic effects is negligible for equilibrium ion dynamics in warm dense hydrogen.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Signaling and scrambling with strongly long-range interactions

Strongly long-range interacting quantum systems—those with interactions decaying as a power law 1/$r^α$ in the distance $\textit{r}$ on a $\textit{D}$-dimensional lattice for $\textit{α ≤ D}$—have received significant interest in recent years. They are present in leading experimental platforms for quantum computation and simulation, as well as in theoretical models of quantum-information scrambling and fast entanglement creation. Since no notion of locality is expected in such systems, a general understanding of their dynamics is lacking. In a step towards rectifying this problem, we prove two Lieb-Robinson-type bounds that constrain the time for signaling and scrambling in strongly long-range interacting systems, for which no tight bounds were previously known. In this work, our first bound applies to systems mappable to free-particle Hamiltonians with long-range hopping, and is saturable for $\textit{α ≤ D/2}$. Our second bound pertains to generic long-range interacting spin Hamiltonians and gives a tight lower bound for the signaling time to extensive subsets of the system for all $\textit{α < D}$. This many-site signaling time lower bounds the scrambling time in strongly long-range interacting systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Mass measurement of 27 P to constrain type-I x-ray burst models and validate the isobaric multiplet mass equation for the $A = 27, T = \frac{3}{2}$ isospin quartet

Background: Light curves are the primary observable of type-I x-ray bursts. Computational x-ray burst models must match simulations to observed light curves. Most of the error in simulated curves comes from uncertainties in rp process reaction rates, which can be reduced via precision mass measurements of neutron-deficient isotopes in the rp process path. Purpose: Perform a precise atomic mass measurement of 27 P. Use this new measurement to calculate rp process reaction rates and input these rates into an x-ray burst model to reduce simulated light curve uncertainty. Use the mass measurement of 27 P to validate the isobaric multiplet mass equation (IMME) for the A = 27 T = 3/2 isospin quartet which 27 P belongs to. Method: High-precision Penning trap mass spectrometry utilizing the time-of-flight ion cyclotron resonance technique was used to determine the atomic mass of 27 P. The MESA code (Modules for Experiments in Stellar Astrophysics) was then used to simulate x-ray bursts using a one-dimensional multizone model to produce updated light curves. Results: The mass excess of 27 P was measured to be –670.7(6) keV, a 14-fold precision increase over the mass reported in the 2020 Atomic Mass Evaluation (AME2020). The 26 Si(p, $γ$) 27 P– 27 P($γ$, p) 26 Si rate equilibrium has been determined to a higher precision based on the precision mass measurement of 27 P. x-ray burst light curves were produced with the MESA code using the new reaction rates. Changes in the mass of 27 P seem to have minimal effect on light curves, even in burster systems tailored to maximize impact. Conclusion: The mass of 27 P does not play a significant role in x-ray burst light curves. It is important to understand that more advanced models do not just provide more precise results, but often qualitatively different ones. This result brings us a step closer to being able to extract stellar parameters from individual x-ray burst observations. In addition, the IMME has been validated for the A = 27, T = 3/2 quartet. The normal quadratic form of the IMME using the latest data yields a reduced $χ$ 2 of 2.9. Finally, the cubic term required to generate an exact fit to the latest data matches theoretical attempts to predict this term.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Evaluation of the Performance and Exhaust Emissions of a 4 Cylinder CI Engine Operating With Dimethyl Ether (DME) and Propane Blends

In response to stringent emissions regulations and the need for higher efficiency engines, the utilization of DME and propane fuel blends in compression ignition (CI) engines has gained interest in the automotive industry. In this study, a range of DME-propane blends are explored in a CI combustion strategy at high injection pressures. A GT-Power model of a 2.2 L Hyundai CI engine was developed to facilitate evaluation of the impacts of variations of DME and propane blends at a light and medium engine speed-torque-load operating condition; specifically at 1500 rpm, 50 Nm and 2.84 bar brake mean effective pressure (BMEP); and 2000 rpm, 150 Nm, 8.53 bar BMEP speed-torque-load combination. Here, the GT-Power model was validated using Computational Fluid Dynamics (CFD) simulations. The results indicate that high diesel-like efficiencies can be achieved with a 100% DME mass substitution and up to 50% propane-DME blends could be implemented without a significant penalty on engine performance indicators. Significant brake specific nitrogen oxides (BSNOx) reductions were also observed along with reductions in carbon dioxide (CO2) and soot when leveraging these fuel blends.

computational fluid dynamics↗
DOE OSTI2020

Gas Generation Models for Thermal Runaway in Lithium Ion Batteries

A numerical model is implemented to study cell venting, internal pressure, and gas-phase dynamics behavior of Li-ion cells undergoing thermal runaway. A k-ɛ Reynolds-Averaged Navier-Stokes (RANS) model is adopted to describe the turbulent flow out of the cells, while the fluid dynamics inside the cells is described by Darcy-Forchheimer’s equation. Thermal abuse reactions and gas generation kinetics are described by a single-step lumped reaction model. A series of computational fluid dynamics (CFD) simulations are conducted on a single cell at various states-of-charge (100%, 50%, 25%) to demonstrate detailed flow and thermal behavior as a function of quantity of gas generated during cell venting. Venting events are categorized into two stages: i) breaching of the cell container and ii) thermal runaway reactions. It is shown that the cell response is dominated by the second stage since most of the gases are generated during thermal runaway. Also, the propensity for propagation is highly affected by state-of-charge (SOC). Cells at higher SOCs produce more heat and gas during the venting event, owing to higher mass and concentrations of reacting gases, and consequently reach higher internal cell pressures which increase the risk of side-wall breaching.

Kim, Jinyong↗
DOE OSTI2021

The Trinity High Explosive Implosion System: The Foundation for Precision Explosive Applications

This article is set during the 1944 and 1945 final push to complete Project Y—the Manhattan Project at Los Alamos—and focuses primarily on overcoming the challenge of creating and demonstrating a successful convergent explosive implosion to turn a subcritical quantity of plutonium into a critical mass. The critical mass would then efficiently yield kilotons of trinitrotoluene (TNT)-equivalent energy in about a microsecond, demonstrating the implosion atomic bomb concept. This work culminated in the Trinity atomic test near Alamogordo on July 16, 1945. This implosion effect demarcated the approach to explosive science and technology the Laboratory has followed ever since, including development of high-explosive synthesis and formulation, small and large test and diagnostic facilities, shock dynamics theory, high-explosive system design engineering, and three-dimensional implosion modeling and simulation using some of the fastest computers in the world. This work also ushered in a period of broader application of precision high explosives in conventional munitions, demolition, mining and oil exploration, and space travel.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Grain2Mesh: Mesh Generation for Grain-Scale Nonlinear Elasticity Modeling

The nonlinear hysteretic behavior of rocks under cyclic loading is a crucial area of study in geomechanics. The macroscopic response of a variety of materials has been found to be contingent upon the behavior of the micro-scale structure. This project aims to develop a functional and maintainable software package for generating a multi-phase numerical mesh and accompanying simulation files for finite element modeling used in computational mechanics solvers. Meshes generated from images often lack key preprocessing that reduces noise and prevents mesh element distortion that can increase computational cost. By incorporating user feedback throughout, grain2mesh ensures a high-fidelity mesh that can be used to model grain-scale interactions such as shearing, crack propagation, and interfacial material contrast. Scientific applications of this software include material fracturing, stress-strain analysis for natural and engineered materials, and nonlinear meso-scale analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Aerodynamic Characterization of 3D Scanned Wind Turbine Blades Using Experimental and Computational Methods

This study presents an aerodynamic characterization of 3D scanned wind turbine blades using both experimental and computational methods. The research was conducted by Gulf Wind Technology and Sandia National Laboratories. The primary objective was to investigate the aerodynamic impacts of leading-edge manufacturing defects on wind turbine blades. The study utilized the Stratasys NEO 800 3D Printer for high-precision manufacturing and the GWT Accelerator Wind Tunnel for experimental testing. Computational simulations were performed using COMSOL Multiphysics to model the wind tunnel and analyze flow characteristics and OpenFOAM to study the aerodynamic impacts of leading-edge defects. OpenFAST was used to estimate how these defects can lead to revenue losses for wind farm operators as high as 6%. The results demonstrated significant aerodynamic performance variations due to defects, with detailed analysis provided through wind tunnel and CFD data. The findings contribute to the understanding of defect impacts on wind turbine blade performance and offer insights for future design improvements.

17 WIND ENERGY↗
DOE OSTI2020

Increasing Compressed Gas Energy Storage Density Using CO2–N2 Gas Mixture

This paper demonstrates a new method by which the energy storage density of compressed air systems is increased by 56.8% by changing the composition of the compressed gas to include a condensable component. A higher storage density of 7.33 MJ/m3 is possible using a mixture of 88% CO2 and 12% N2 compared to 4.67 MJ/m3 using pure N2. This ratio of gases representing an optimum mixture was determined through computer simulations that considered a variety of different proportions from pure CO2 to pure N2. The computer simulations are based on a thermodynamic equilibrium model that predicts the mixture composition as a function of volume and pressure under progressive compression to ultimately identify the optimal mixture composition (88% CO2 + 12% N2). The model and simulations predict that the optimal gas mixture attains a higher energy storage density than using either of the pure gases.

25 ENERGY STORAGE↗
DOE OSTI2022

Lignin Utilization

Valorization of lignin has massive potential economic and sustainability benefits for the lignocellulosic biorefinery. However, challenges remain to realize lignin conversion to coproducts, especially related to selective, high yield depolymerization to monomers and quantitative analytics on lignin, the latter of which is critical for accurate process modeling. Towards these goals, the Lignin Utilization project focuses on catalytic lignin deconstruction chemistries for C-O and C-C bond cleavage, development of analytical chemistry techniques to quantitatively characterize lignin, and syntheses of requisite model compounds for understanding lignin transformations. This work is done with and supports multiple BETO projects, including the Biological Lignin Valorization (BLV) project, the Separations Consortium, and others. Outcomes of Lignin Utilization for analytics and synthesis include 1) development of new mass spectrometry (MS) methods to characterize lignin dimers/oligomers in process streams, 2) deployment of a computational-experimental tool (with Biochemical Process Modeling and Simulation) to identify lignin-derived compounds with high fidelity from MS, and 3) delivery of >40 unique compounds. From a catalysis perspective, we have developed new oxidative approaches to cleave C-O and C-C bonds and produce >50% bio-available aromatic monomers for biological funneling for the BLV project and have developed recoverable bases for base-catalyzed deconstruction of lignin.

BIOMASS FUELS↗
DOE OSTI2024

An open-source framework for balancing computational speed and fidelity in production cost models

Studies of bulk power system operations need to incorporate uncertainty and sensitivity analyses, especially around exposure to weather and climate variability and extremes, but this remains a computational modeling challenge. Commercial production cost models (PCMs) have shorter runtimes, but also important limitations (opacity, license restrictions) that do not fully support stochastic simulation. Open-source PCMs represent a potential solution. They allow for multiple, simultaneous runs in high-performance computing environments and offer flexibility in model parameterization. Yet, developers must balance computational speed (i.e. runtime) with model fidelity (i.e. accuracy). In this paper, we present Grid Operations (GO), a framework for instantiating open-source, scale-adaptive PCMs. GO allows users to search across parameter spaces to identify model versions that appropriately balance computational speed and fidelity based on experimental needs and resource limits. Results provide generalizable insights on how to navigate the fidelity and computational speed tradeoff through parameter selection. We show that models with coarser network topologies can accurately mimic market operations, sometimes better than higher-resolution models. It is thus possible to conduct large simulation experiments that characterize operational risks related to climate and weather extremes while maintaining sufficient model accuracy.

42 ENGINEERING↗
DOE OSTI2021

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗