Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation and Modeling Tools”

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,081 records · Page 60

Assessment and Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY25 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., the 200 DV-1 Operable Unit), amendment particles need to be delivered a reasonable distance away from an injection well to provide cost effective in situ treatment. Field-scale particle transport models can estimate spatial deposition of amendment particles in the subsurface, which is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capability and a lack of experimental data to validate and parameterize particle transport models. During this fiscal year, the following progress has been made toward a field-scale particle transport modeling evaluation: (1) in addition to the two particle transport models implemented last FY, four additional particle transport models have been implemented within PFLOTRAN; (2) a Python-based pre-screening tool was finalized, enabling users to quickly estimate the particle radius of influence (ROI) for any given particle-amendment system; and (3) an initial compatibility assessment was completed using both particle transport simulations deployed through the pre-screening tool, in conjunction with general guidelines to (a) identify the most suitable amendment particle sizes for various Hanford sediments and (b) evaluate amendment-delivery fluid compatibility. The preliminary compatibility assessment revealed that for amendment delivery success to the various Hanford target formations, amendment particle sizes will likely need to be smaller than the amendment sizes tested in the DV-1 treatability study. It is recommended that amendment particles be decreased in size, or alternative smaller size amendments be obtained from the manufacture, prior to any further experimental testing. Also, preliminary testing suggests that xanthan gum may be the most broadly compatible delivery fluid. Planned laboratory experiments will be instrumental in validating and refining the PFLOTRAN particulate transport model formulations, ultimately enabling predictive capabilities to facilitate the design of field-scale amendment delivery systems. This work consists of acquiring new theoretical or experimental knowledge. The information associated with this report should not be used as design input or operating parameters without additional qualification.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

3D Theoretical and simulation tools for microbunched cooling

The concept of microbunched electron cooling (MBEC, introduced in Ref. [1]) is currently under investigation as a candidate for the strong hadron cooling (SHC) component of EIC. Most of the theoretical treatments of this scheme are based on a simplified model for the longitudinal space charge [2–4]. In this model, the particle beams are treated as collections of thin charged disks, with the interaction function for two such disks being calculated by averaging the forces between constituent point charges. Although the disk model is flexible and captures a large part of the underlying physics, it lacks the fidelity and robustness of an analysis that explicitly treats the beams as collections of point charges. In what follows, we explore such a model, both from a theory-based and a simulation perspective.

43 PARTICLE ACCELERATORS↗

Evaluation of distributed process-based hydrologic model performance using only a priori information to define model inputs

Fully distributed, integrated surface–subsurface hydrological models (ISSHMs) have seen renewed interest due to availability of better software, high performance computing facilities, and high-resolution, spatially extensive data products. ISSHMs are valuable as tools for advancing system understanding as they can resolve multiple processes defined on the plot scale including three-dimensional interaction of surface water and groundwater. Here, we evaluated the performance of an ISSHM, the Advanced Terrestrial Simulator (ATS), on seven diverse catchments across the continental US using widely available data products to define model inputs without calibration. We compare the ATS-simulated streamflow and evapotranspiration with gauge observations and MODIS-derived evapotranspiration, respectively. Using the Kling-Gupta Efficiency (KGE) as metric, ATS with default data products performed reasonably well at 6 of 7 catchments for streamflow. However, in one of those 6 catchments ATS had poor performance on baseflow and ATS’s overall performance was thus judged to be inadequate despite the acceptable KGE. ATS performance for evapotranspiration was good in all 7 catchments using default data products. In the two catchments where ATS streamflow performance using default data products was not acceptable, the performance was significantly improved by using local information on subsurface properties below the soil. We also compare the model-simulated streamflow and evapotranspiration with the Sacramento soil moisture accounting (SAC-SMA) model, a semi-distributed model that was calibrated on a catchment-by-catchment basis. Uncalibrated ATS performance is comparable to the calibrated SAC-SMA model in terms of streamflow while ATS performance is similar to or better (much better in certain catchments) in reproducing MODIS-derived evapotranspiration. Reasonably good performance of ATS without catchment-specific calibration provides new confidence in the ISSHM class of models and community data products as tools for advancing understanding of watershed function in a changing environment.

54 ENVIRONMENTAL SCIENCES↗

srlife: A Fast Tool for High Temperature Receiver Design and Analysis

This report describes a tool for estimating the structural service life of tubular, panel solar receivers operating at high temperatures. A complete version of the tool is available as open source software at https://github.com/Argonne-National-Laboratory/srlife and can be installed through the PyPi (https://pypi.org) package manager. Given the basic receiver geometry and the thermal loads on the receiver, the tool provides 1D, 2D, or 3D thermal and structural (single tube and simplified system) analysis and creep-fatigue service life prediction for six metallic alloys – 316H, 800H, Alloy 617, Alloy 740H, and Alloy 230. With the exception of Alloy 230 where the material data is preliminary, the software included a detailed set of models for the materials, well-supported by high temperature experimental test data. The tool is designed for easy integration with a software stack, including solar field and thermohydraulic simulations, for optimizing receiver designs to meet service life and economic targets. The report describes several heuristics that can be applied in srlife to reduce the analysis time by several orders of magnitude but with fairly accurate life estimation when compared with full analysis. The report provides several examples demonstrating the utility of srlife in receiver design. Finally, the report discusses high temperature tests on Alloy 282, collected as part of this project, used to develop and support the material model for that alloy.

47 OTHER INSTRUMENTATION↗

Simulating energy performance of buildings: a study using eQUEST and Energy Star ® portfolio manager

Commercial buildings consume significant energy in the United States and exhibit high potential for energy use reduction through retrofits. Benchmarking and energy simulation are well established tools in the industry to identify potential improvements and measure performance. Analysis to identify most sensitive retrofit parameters to energy performance can optimize investment and available energy savings. Presented study demonstrates methodology using a static model to determine sensitivity of building design and retrofit parameters with respect to energy performance. Calibrated simulation energy models (eQUEST) of two distribution centers (A, B) are presented. A fractional factorial analysis is conducted on retrofit parameters of efficiency measures targeting the highest energy consumers, and the results are benchmarked using Energy Star® Portfolio Manager. A custom Microsoft Excel® based simulation model is created to simulate occupancy levels, lighting, plug loads, and other equipment used in various spaces throughout the day. For Building A, efficient lighting was the most influential parameter for energy savings, carbon savings and benchmarking score; whereas, for Building B, HVAC efficiency was most influential for energy and demand controlled ventilation and economizers was most influential for benchmarking score. While retrofit projects can save energy and carbon emissions, variation in source-site ratios and state grid emissions, benchmarking scores may not always reflect equivalent improvement. State grid emissions factors, natural gas composition are difficult to model and hence not considered in this study. In conclusion, the synergistic analysis presented, emphasizes the importance of benchmarking and efficiency retrofits in promoting sustainable building practices to reduce energy consumption.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

An Open-Source Parallel EMT Simulation Framework

As the integration level of inverter-based resources (IBRs) increases, ensuring the reliable operation of the bulk power systems requires the use of electromagnetic transient (EMT) simulation tools to identify and mitigate system-wide stability risks. Conducting EMT studies for large-scale, IBR-rich grids, however, is challenging due to the inherent computational bottleneck caused by the underlying high-fidelity models and required small time steps. This paper introduces ParaEMT: an open-source, generic EMT simulation framework designed to accelerate simulations by leveraging advanced parallel computational technologies, such as high-performance computers. This paper presents a comprehensive exposition of ParaEMT, covering its modeling library, simulation strategy, framework structure, operational procedures, and auxiliary features, alongside its extensible parallel computational architecture. Notably, ParaEMT is a publicly accessible and modularized framework written in Python, thereby facilitating future development and the integration of new models and algorithms. The accuracy and efficiency of ParaEMT are demonstrated by rigorous validations via multiple case studies.

electromagnetic transient simulation↗

An Open-Source Parallel EMT Simulation Framework: Preprint

As the integration level of inverter-based resources (IBR) increases, ensuring the reliable operation of the bulk power systems requires the use of electromagnetic transient (EMT) simulation tools to identify and mitigate system-wide stability risks. Conducting EMT studies for large-scale, IBR-rich grids, however, is challenging due to the inherent computational bottleneck caused by the underlying high-fidelity models and required small time steps. This paper introduces ParaEMT: an open-source, generic EMT simulation framework designed to accelerate simulations by leveraging advanced parallel computational technologies, such as high-performance computers. This paper presents a comprehensive exposition of ParaEMT, covering its modeling library, simulation strategy, framework structure, operational procedures, and auxiliary features, alongside its extensible parallel computational architecture. Notably, ParaEMT is a publicly accessible and modularized framework written in Python, thereby facilitating future development and the integration of new models and algorithms. The accuracy and efficiency of ParaEMT are demonstrated by rigorous validations via multiple case studies.

electromagnetic transient simulation↗

Tool-workpiece stick-slip conditions and their effects on torque and heat generation rate in the friction stir welding

Friction stir welding (FSW) has found increased applications in automotive and aerospace industries due to its advantages of solid-state bonding, no fusion and melting, and versatility in various working conditions and material combinations. However, the relationship among processing parameters, material properties, and bonding extent and fidelity remains largely empirical, primarily because of the lack of the mechanistic understanding of the tool-workpiece frictional behavior that affects our subsequent understanding of microstructural evolution and interface bonding formation. While the tool-workpiece stick-slip condition is believed to dictate the resulting torque and heat generation rate during the welding process, it remains rare and elusive to conduct a quantitative experimental measurement of such interfacial field. On the other hand, numerical simulations based on Computational Fluid Dynamics (CFD) rely on ad hoc assumptions of interfacial pressure and shear-stress conditions, but predictions can only be validated via the medium- and far-range temperature field which is known to be insensitive to the interfacial frictional behavior. This work first presents a comparison among two CFD-based simulation methodologies and the Coupled Eulerian Lagrangian (CEL) model in finite element method, the last of which uses the Coulomb friction so that the stick-slip is naturally developed. Based on the Hill-Bower similarity relationship in the contact analysis, an analytical model is developed here to prove why a constant stick-slip fraction will be developed in the steady state, to correlate the stick-slip fraction to processing parameters such as the tool spin rate, and further to derive dimensionless functions for torque and heat-generation-rate predictions. Pros and cons of various numerical approaches in predicting stick-slip are discussed, and our analytical model has been found to agree well with our numerical simulation and literature experimental results. These analyses provide the critical strain-rate and temperature fields that are needed for the bonding analysis in our future work.

36 MATERIALS SCIENCE↗

Exploring the Use of Novel Spatial Accelerators in Scientific Applications

Driven by the need to find alternative accelerators which can viably replace GPUs in next-generation Supercomputing systems, this paper proposes a methodology to enable agile application/hardware co-design. The application-first methodology provides the ability to come up with design of accelerators while working with real-world workloads, available accelerators, and system software. The iterative design process targets a set of kernels in a workload for performance estimates that can prune the design space for later phases of detailed architectural evaluations. To this effect, in this paper, a novel data-parallel device model is introduced that simulates the latency of performance-sensitive operations in an accelerator including data transfers and kernel computation using multi-core CPUs. The use of off-the-shelf simulators, such as pre-RTL simulator Aladdin or multiple tools available for exploring the design of deep neural network accelerators (e.g., Timeloop) is demonstrated for evaluation of various accelerator designs using applications with realistic inputs. Examples of multiple device configurations that are instantiable in a system are explored to evaluate the performance benefit of deploying novel accelerators. The proposed device is integrated with a programming model and system software to potentially explore the impacts of high-level programming languages/compilers and low-level effects such as task scheduling on multiple accelerators. We analyze our methodology for a set of applications that represent high-performance computing (HPC) and graph analytics. The applications include a computational chemistry kernel realized using tensor contractions, triangle counting, GraphSAGE and Breadth-first Search. These applications include kernels such as dense matrix-dense matrix multiplication, sparse matrix-spare matrix multiplication, and sparse matrix-dense vector multiplication. Our results indicate potential performance benefits and insights for system design by including accelerators that realize these kernels along-side general purpose accelerators.

AI, codesign, Accelerated Computing, Modeling and ↗

Field Testing of a Hierarchical Model-Free Transactive Control Strategy in a Residential House

Demand response plays an important role in addressing the uncertain, intermittent, and variable nature of renewable energy sources. However, to be effective, it requires a significant amount of demand load flexibility, in which buildings are important flexible assets for electricity systems. In this regard, many transactive control tools have been developed to harness this demand flexibility. Yet, such tools are mostly simulation-based or deal with a single load, which may not be sufficient to demonstrate their effectiveness. Towards this end, this paper presents a field implementation and testing of a hierarchical model-free transactive control strategy using residential building thermostatically controlled loads. The transactive control implementation incorporates elements of virtual battery, game theory, and model-free control mechanisms. Experimental field-testing results are provided to demonstrate the effectiveness of the proposed transactive control strategy.

Olama, Mohammed M.↗

Investigation of irradiation damage and heat deposition: a comparative analysis for HEU-to-LEU conversion in HFIR

The planned conversion of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel requires detailed evaluation of experiment-relevant parameters to ensure continued performance for materials testing and isotope production. Here, this study presents the first comprehensive assessment of displacements per atom (dpa) and heat deposition rates in target materials within the HFIR flux trap with both HEU and candidate LEU core configurations. Seven analyses were conducted to evaluate key performance metrics, including fast neutron flux distribution, cross section response functions, cross section data, and local dpa and heat deposition rates using mesh- and cell-based tallies. Simulations employed Shift, Monte Carlo N-Particle (MCNP), and the HIFR Controller (HFIRCON) tool suite for high-fidelity transport and depletion modeling. The LEU designs—using U 3 Si 2 -Al dispersion fuel and operating at 95 MW—were compared to the current 85 MW HEU configuration. Results show that while the candidate LEU cores exhibit higher dpa rates due to a harder spectrum and extended cycle lengths, they also demonstrate reduced heat deposition rates in irradiation experiments, primarily due to increased gamma self-shielding from higher 238 U content in the core. These findings confirm that LEU conversion can maintain HFIR’s materials irradiation capabilities but may require redesigning existing experimental hardware.

HEU↗

Extreme variation in fatigue: Fatigue life prediction and dependence on build volume location in laser powder bed fusion of 17-4 stainless steel

Laser powder bed fusion (LPBF), a metal additive manufacturing technology, is well-suited for design optimization but fatigue life is limited by manufacturing defects. In this work, 17-4 stainless steel components were manufactured in densely populated build volumes, simulating at-scale LPBF production. Tests revealed extreme variability in fatigue life data, analyzed via rigorous statistical tools. The El-Haddad model, modified for finite-life, enabled defect-based life prediction. Specimen location within the build volume correlated to life, which was heteroscedastic. Investigating defect concentration over the build volume explained typical life and scatter. Finally, these findings argue for qualification approaches which acknowledge high material lot variability.

36 MATERIALS SCIENCE↗

Unraveling small-scale defects in irradiated ThO 2 using kinetic Monte Carlo simulations

Point defects and their clusters generated through irradiation can have significant impact on the physical and mechanical properties of materials. However, direct experimental visualization of these small-scale defects using high-resolution scanning transmission electron microscopy remains a challenging task. In this study, using thorium dioxide (ThO 2 ) with the fluorite structure as a model system, we demonstrate the use of ab initio basin-hopping simulations in synergy with object kinetic Monte Carlo simulations as a powerful tool for identifying small defect complexes in irradiated materials. In addition to providing quantitative insights into defect evolution in ThO 2 under irradiation, our study reveals an unexpected role of bound anti-Schottky defect clusters in mediating defect transport. Remarkably, despite their poor thermal stability against dissociation at high temperatures, the transient formation of bound anti-Schottky defects under irradiation and their subsequent migration provide the dominant mechanism for the growth of large interstitial loops that have been experimentally observed in ThO 2 .

36 MATERIALS SCIENCE↗

Role of the Terrestrial Biosphere in Atmospheric Chemistry and Climate

The terrestrial biosphere–atmosphere interface provides a key chemical, biological, and physical lower boundary for the atmosphere. The presence of vegetation itself modifies the physical boundary, or the biogeophysical aspects of the system, by controlling important climate drivers such as soil moisture, light environment, and temperature. The leaf surface area of the terrestrial biosphere provides additional surface area for emissions, and it can be up to 55% of the total Earth’s surface area during the boreal summer. Vegetation also influences the biogeochemical aspects of the system by emitting a broad suite of reactive trace gases such as biogenic volatile organic compound (BVOC) emissions and climate-relevant primary biological aerosol particles (PBAP). Many of these emissions are a function of meteorological and climatological conditions at the surface, including temperature, light environment, soil moisture, and winds. Once emitted, they can be processed in the troposphere through a suite of chemical reactions. BVOC can contribute to the formation of ozone and secondary organic aerosols (SOA), and PBAP can rupture to form smaller particles with climatic relevance. These emissions and subsequent aerosol products can influence atmospheric processes that affect the surface climate, such as the attenuation of radiation, the formation of greenhouse gases such as ozone that can feedback to surface air temperature, and the alteration of clouds and subsequent precipitation. These atmospheric changes can then feedback to the land surface and emissions themselves, creating positive or negative feedback loops that can dampen or amplify the emission response. For the dominant BVOC isoprene, the feedback response to temperature can be positive or negative depending on ambient temperatures that drive isoprene emissions. The feedback response to soil moisture and precipitation can be positive, negative, or uncoupled depending on the moisture content of the soil and the total atmospheric aerosol loading. For light, the isoprene response can be positive or negative depending on the role of diffuse light. Overall, these feedbacks highlight the dynamical response of the biosphere to changing atmospheric conditions across a range of time scales, from minutes for trace gases and aerosols, to months for phenological changes, to years for land cover and land use change. Here, the dynamic aspect of this system requires us to understand, simulate, and predict the complex feedbacks between the biosphere and atmosphere and understand their role in the simulation and understanding of climate and global change. From the observational perspective, these feedbacks are challenging to identify in observations, and predictive modeling tools provide a crucial link for understanding how these feedbacks will change under warming climate scenarios.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Penetration through slots in cylindrical cavities with cavity modes overlapping with the first slot resonance

In this work, we analyze the coupling into a slotted cylindrical cavity operating at fundamental cavity modal frequencies overlapping with the slot’s first resonance frequency through an unmatched formulation that accounts for the slot’s absorption and radiation processes. The model is validated through full-wave simulations and experimental data. We then couple the unmatched formulation to a perturbation theory model to investigate an absorber within the cavity to reduce the interior field strength, also validated with full-wave simulations and experiments. These models are pivotal to understanding the physical processes involved in the electromagnetic penetration through slots, and may constitute design tools to mitigate electromagnetic interference effects within cavities.

42 ENGINEERING↗

Towards Redefining the Reproducibility in Quantum Computing: A Data Analysis Approach on NISQ Devices

Although the building of quantum computers has kept making rapid progress in recent years, noise is still the main challenge for any application to leverage the power of quantum computing. Existing works addressing noise in quantum devices proposed noise reduction when deploying a quantum algorithm to a specified quantum computer. The reproducibility issue of quantum algorithms has been raised since the noise levels vary on different quantum computers. Importantly, existing works largely ignore the fact that the noise of quantum devices varies as time goes by. Therefore, reproducing the results on the same hardware will even become a problem. We analyze the reproducibility of quantum machine learning (QML) algorithms based on daily model training and execution data collection. Our analysis shows a correlation between our QML models’ test accuracy and quantum computer hardware’s calibration features. We also demonstrate that noisy simulators for quantum computers are not a reliable tool for quantum machine learning applications.

Senapati, Priyabrata↗

DOE-ICoM/RIFT

Rapid Infrastructure Flood Tool (RIFT) is a two-dimensional hydrodynamic model based on the complete shallow water equations. RIFT has specifically been designed with rapid simulation in mind by utilizing commodity high performance computing technology and best-available nation-wide data. RIFT is used to predict the movement of water over land and resolve the spatial and temporal variability of flood depths, extent, and velocity. RIFT can be applied to many flood situations and has primarily been used to quantify flood extents from dam/levee failure or inland rainfall flooding.

Perkins, Bill [Pacific Northwest National Laborato↗

SANDWake3D

SAND2026-23958O SANDWake3D models atmospheric surface layers and turbine wakes using a 3D parabolic RANS solver. This Python-based tool uses an efficient method to simulate wakes, capturing complex wake behaviors such as wind shear, veer, and atmospheric stratification, naturally. It has been calibrated and validated against higher-fidelity large eddy simulations (LES), providing reliable results for wind energy and atmospheric research applications. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗