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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High Fidelity Computational Model for Fluidized Bed Experiments (Final Report)

MFiX was integrated with preconditioners and linear solver packages in Trilinos via MFIX, Fortran, C, and C++ wrappers. The MFIX wrapper interpreted the structure of the matrix and vector from MFiX while the Fortran wrapper transfers this information to a low-level language, C. C-wrapper transforms the memory semantic between Fortran and C++ language. C++ wrapper passes the matrix and vector to a Trilinos package and gets the solution. This solution was transferred to MFiX via C++, C, Fortran, MFIX wrapper. The framework is used to integrate first, second and third generation linear solvers in Trilinos with MFiX.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Computational Modeling and Simulation for Nonproliferation: The History (and Present) of Monte Carlo and the MCNP(R) Code at Los Alamos [Slides]

The emergence of the Monte Carlo method as a research tool springs from work done at Los Alamos in the 1940s.Monte Carlo and the MCNP Code have been and continue to be developed at Los Alamos for many decades. From basic science in support of understanding nuclear interaction physics to global security applications, application uses of the code are extensive. Recent R&D projects and code modernization efforts make the MCNP code a great tool for nuclear nonproliferation applications. In collaboration with nuclear safeguards experts, new training has just recently been developed to help new practitioners learn how to use the code for nuclear safeguards applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Augmenting a Simulation Campaign for Hybrid Computer Model and Field Data Experiments

The Kennedy and O’Hagan (KOH) calibration framework uses coupled Gaussian processes (GPs) to meta-model an expensive simulator (first GP), tune its “knobs” (calibration inputs) to best match observations from a real physical/field experiment and correct for any modeling bias (second GP) when predicting under new field conditions (design inputs). There are well-established methods for placement of design inputs for data-efficient planning of a simulation campaign in isolation, that is, without field data: space-filling, or via criterion like minimum integrated mean-squared prediction error (IMSPE). Analogues within the coupled GP KOH framework are mostly absent from the literature. Here, in this study, we derive a closed form IMSPE criterion for sequentially acquiring new simulator data for KOH. We illustrate how acquisitions space-fill in design space, but concentrate in calibration space. Closed form IMSPE precipitates a closed-form gradient for efficient numerical optimization. We demonstrate that our KOH-IMSPE strategy leads to a more efficient simulation campaign on benchmark problems, and conclude with a showcase on an application to equilibrium concentrations of rare earth elements for a liquid–liquid extraction reaction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Computational Modeling of Water Infiltration From Ditches Into Roadway Embankments

Ditches are built along roadways to convey the stormwater runoff. They collect the runoff from the pavement and adjacent drainage areas. Water can move freely in a ditch away from the road surface, thus making it safer for the vehicular traffic. South Carolina DOT expressed interest in studying infiltration of water from ditches into roadway embankments. As part of the current design criteria, it is recommended that the road subgrades be 1.0 foot above the design event water surface elevation in the ditch. The distance between the water surface and bottom of the subbase is referred to as “freeboard”. The purpose of this requirement is to lower the risk of pavement damage by moving the flow away from the pavement, containing the captured flow, and conveying flow received from adjacent properties to outfalls. With changes in roadway design standards, larger shoulders and flatter cross-slopes have moved ditches further away from the pavement and the subbase. These changes could allow for a lower freeboard standard without reducing safety margins due to the increased lateral distance from the subbase to the ditch. The objective of this study is to evaluate whether the current SCDOT criterion requiring the roadway subgrade to be 1 ft above the design-event ditch water surface can be relaxed under selected geometric and soil conditions, particularly where larger shoulder widths and flatter side slopes increase horizontal separation between ditch water and pavement structure.

42 ENGINEERING↗
DOE OSTI2026

CRCNS21 Computational Models of Multisensory Integration by Upper Limb in Humanoids and Amputees

This international collaborative research project between Johns Hopkins University (JHU) and the Technical University of Munich (TUM) investigated how the human brain processes and integrates multiple types of sensory information, such as touch and force, with the goal of improving prosthetic limbs for amputees and advancing sensory capabilities in humanoid robots. The research advanced our understanding of how the brain responds to sensory feedback in upper-limb amputees. Through experiments in which amputees received electrical stimulation while performing phantom hand movements, we demonstrated that sensory feedback activates the cortical sensorimotor and multisensory regions, and that these regions communicate dynamically during stimulation. Experiments with intact-limb participants explored the integration of visual, haptic, and force feedback, as well as in virtual reality motor training, further showing how the brain processes multimodal sensory information. In addition, this research inspired work on examining the reliability of where amputees perceive sensations over time, which contributed to a successful doctoral fellowship for continued investigation. Our collaborators at TUM improved multimodal sensor technology combining tactile and thermal feedback for humanoid robots, demonstrating the feasibility of integrating multiple sensor types into a unified system for detecting and responding to environmental stimuli. The experimental methods and analysis techniques developed across both teams, including functional network analysis and multimodal sensor integration, provide a foundation for future research in prosthetics and robotics. This research benefits the public by generating knowledge about how amputees process restored sensory information. Advances in humanoid sensing contribute to safer human-robot interaction. The project also fostered international collaboration and cross-disciplinary training: one TUM doctoral student spent a summer at JHU working on multimodal sensor integration, while two JHU students traveled to TUM to host workshops on neuromorphic sensory encoding and sensory integration.

42 ENGINEERING↗
DOE OSTI2024

Computational modeling of microalgal biofilm growth in heterogeneous rotating algal biofilm reactors (RABRs) for wastewater treatment

Rotating algal biofilm reactors (RABRs) are innovative systems designed to cultivate microalgae biofilms efficiently. In this paper, we have developed a novel mathematical model to accurately capture the growth dynamics of algae biofilms within RABR. By considering the spatial heterogeneity of the RABR, we introduce a PDE-based model that addresses the spatial variations across the substratum, enabling a more accurate simulation of biofilm growth in RABRs. The photosynthesis process is modeled through reactive kinetics, driving the growth of the algae biofilm. To analyze the system's behavior, we employ finite difference numerical methods to solve the complex PDE model. We then conduct extensive numerical simulations to understand algae biofilm growth in the RABR environment under various operational factors and environmental conditions. One primary focus in these simulations is to investigate the impact of various harvesting strategies, harvesting frequencies, light intensity, and light exposure on the overall biomass productivity of the algae biofilm. The numerical results provide valuable insights into optimizing algae biofilm growth and designing harvesting techniques in RABR systems. Our proposed novel mathematical model provides an effective platform for the theoretical investigation and design of RABRs for wastewater treatment.

09 BIOMASS FUELS↗
DOE OSTI2026

Computational modeling of graphite degradation in molten salt reactors: Role of infiltration

Molten salt reactors (MSRs) often employ graphite as a moderator and reflector. An important challenge for deploying graphite in these reactors is that, due to limited experimental data, our understanding of graphite’s structural integrity in molten salt environments remains incomplete. Here, this study addresses heat generation from fuel-bearing salt that has infiltrated open pores in the graphite, driven primarily by pressure differentials. This is one of multiple identified physical and chemical mechanisms through which molten salt could potentially degrade graphite. Thermally driven stresses are quantified using the Molten-Salt Reactor Experiment (MSRE) graphite moderator elements as a case study. Finite element simulations predict stress distributions at varying infiltration levels, indicating that thermal stresses increase with higher infiltration. Rare-event simulations using the parallel subset simulation framework identify the combinations and corresponding ranges of input parameters that lead to stresses above a specified threshold. In particular, combinations involving high infiltration amounts, high power density, and low thermal conductivity tend to induce the highest stresses. Under the inputs and assumptions considered in this work, the magnitudes of the thermally driven stresses are quite low, with a very low likelihood of causing failure due to exceeding the graphite’s tensile strength. Additionally, rare-event simulations were performed for two more scenarios: a scaled-up moderator geometry and a localized hotspot in the original geometry. Both cases resulted in increased susceptibility to failure, though not to a detrimental extent. Furthermore, the combined effects of irradiation and infiltration-induced thermal stresses were evaluated. The results showed that thermal stresses from infiltration were negligible compared to those caused by irradiation. The findings of such a study are inherently component-specific, but the methodology presented here could be used for similar assessments of salt-infiltration effects in other graphite components.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Micro-structural features and material properties impact on adhesive metal joints via computational modeling and machine learning

The quality of structural bonding in practical applications depends on various factors arising from materials, pre-processing conditions, and manufacturing. Understanding how these factors influence bonding performance and determining their relative importance are of significant interest. Thus, this study evaluates the effects of microstructural features and material properties on the structural strength of adhesively-bonded metal joints at the submillimeter scale, utilizing a combination of Finite Element Modeling (FEM) and Machine Learning (ML) with Gradient Boosting Regression (GBR). The microstructural features include adhesive thickness, internal voids within the adhesive, adherend-adhesive interfacial voids, void size and volume fraction, and surface roughness. The material properties include the constitutive behavior of the adhesive, as well as the adherend-adhesive interfacial strength and fracture energy. The changes in structural strength and morphologies of the bonded metal structures with respect to different microstructural features and material properties were clarified by FEM. By further leveraging ML-GBR, the sequence of importance of these factors affecting bonding performance across various scenarios was summarized. This work provides valuable insights into the development of improved structural bonding for adhesive joints in industries such as automotive , aerospace, and beyond.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamics of episodic supershear in the 2023 M7.8 Kahramanmaraş/Pazarcik earthquake, revealed by near-field records and computational modeling

Abstract The 2023 M7.8 Kahramanmaraş/Pazarcik earthquake was larger and more destructive than what had been expected. Here we analyzed nearfield seismic records and developed a dynamic rupture model that reconciles different currently conflicting inversion results and reveals spatially non-uniform propagation speeds in this earthquake, with predominantly supershear speeds observed along the Narli fault and at the southwest (SW) end of the East Anatolian Fault (EAF). The model highlights the critical role of geometric complexity and heterogeneous frictional conditions in facilitating continued propagation and influencing rupture speed. We also constrained the conditions that allowed for the rupture to jump from the Narli fault to EAF and to generate the delayed backpropagating rupture towards the SW. Our findings have important implications for understanding earthquake hazards and guiding future response efforts and demonstrate the value of physics based dynamic modeling fused with near-field data in enhancing our understanding of earthquake mechanisms and improving risk assessment.

Environmental Sciences & Ecology↗
DOE OSTI2020

Computational modeling predicts the stability of both Pd + and Pd2 + ion-exchanged into H-CHA

Passive NO x adsorbers (PNA) using Pd/zeolites have emerged as a promising solution for the reduction of cold-start emissions from vehicle exhaust. However, the nature of the active sites and the mechanisms underlying NO x adsorption in Pd/zeolites remain a subject of ongoing investigation. In this study, we employ quantum chemical simulations to investigate the structure of Pd species in cation-exchange sites at isolated Al and Al pairs in the 6-ring and 8-ring of the CHA framework before the introduction of NO x . Our calculations show that the speciation of Pd in these exchange sites strongly depends on the precise Al arrangement within the framework, as well as the operating conditions. Ionically dispersed Pd is found to be the most favorable species over a wide range of oxidizing and reducing conditions. Small oligomers of PdO and metallic Pd do not appear to be competitive at either isolated Al or Al pairs. Notably, our calculations show that ion exchange sites other than next–next–nearest neighbor Al pairs in the 6-ring will be preferentially occupied by Pd + instead of Pd 2+ . The stability of Pd + in the zeolite environment is an interesting contrast with its rareness in molecular Pd compounds. Nonetheless, a detailed analysis of the electronic structure shows that predicted Pd oxidation states are consistent with chemical intuition for all complexes investigated in this study. We also discuss the potential ambiguity in Pd characterization provided by typical experimental techniques such as XANES, EXAFS and UV-Vis, and highlight the need for additional EPR spectroscopy studies to further elucidate the initial Pd speciation in zeolites for PNA applications.

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

Computational Modeling of Multi-Pass Rolling Parameters Effect on Resulting Fuel Foil Shape

A focus of the U.S. Department of Energy is to improve production yield and reduce the cost of Low Enriched Uranium (LEU)-molybdenum alloy (U-10Mo) monolithic fuel plates that will be replacing High Enriched Uranium (HEU) oxide dispersion fuels used currently in the United States High Performance Research Reactors (USHPRR). One area of improvement is lowering the high transverse waviness and longitudinal waviness currently present within rolled foils prior to cladding to produce fuel plates. Traditional rolling manufacturing techniques for other metal foils use winders to pull and straighten the foil as it is rolled back and forth to the final thickness. This approach cannot be used to roll thin U-10Mo foils (0.008-0.025” thick) because only small castings can be rolled due to nuclear criticality safety concerns. As a result, the fuel foils are too short (1 m in length) to use traditional winders. Therefore, it is crucial to identify other rolling parameters (i.e., roller friction, axial tension load, roller diameter, and roll pass reduction percent) that might reduce transverse waviness and longitudinal waviness in the rolled fuel foil and develop a high-yield, low-cost multi-pass rolling manufacturing process. This report documents a systematic finite element modeling study to investigate the effects of numerous rolling parameters to reduce resulting transverse waviness and longitudinal waviness in the fuel foil during multi-pass rolling of U-10Mo foils. The rolling of a U-10Mo plate with initial dimensions of 1”x1”x 0.048” is modeled using Abaqus CAE. This rolling is modeled to undergo eight 20% reduction roll passes to a final fuel foil thickness of 0.01”. The elastic-plastic constitutive model of the U-10Mo alloy was input to the fuel foil rolling model. The rollers were modeled as rigid bodies. A comparison of rolling friction coefficients of 0.3 and 0.7 over a wide range of applied axial tension loads were investigated in order to evaluate the effect of using a lubricant during rolling. The effect of roller diameter on the resulting transverse waviness and longitudinal waviness of the fuel foil over a wide range of axial tension loads were also investigated by modeling rollers 7/8” and 3.75” in diameter. The results of this systematic finite element method study will aid manufacturers in producing low transverse waviness and reduced longitudinal waviness in U-10Mo fuel foils.

U-10Mo, FEA, Rolling, Residual Stress, Fuel Foil↗