SEARCH · Engineering Papers

Results for “Simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Simulation of Powder Layer Deposition in Additive Manufacturing Processes Using the Discrete Element Method

This document serves as a final report to a small effort where several improvements were added to a LLNL code GEODYN-L to develop Discrete Element Method (DEM) algorithms coupled to Lagrangian Finite Element (FE) solvers to investigate powder-bed formation problems for additive manufacturing. The results from these simulations will be assessed for inclusion as the initial conditions for Direct Metal Laser Sintering (DMLS) simulations performed with ALE3D. The algorithms were written and performed on parallel computing platforms at LLNL. The total funding level was 3-4 weeks of an FTE split amongst two staff scientists and one post-doc. The DEM simulations emulated, as much as was feasible, the physical process of depositing a new layer of powder over a bed of existing powder. The DEM simulations utilized truncated size distributions spanning realistic size ranges with a size distribution profile consistent with realistic sample set. A minimum simulation sample size on the order of 40-particles square by 10-particles deep was utilized in these scoping studies in order to evaluate the potential effects of size segregation variation with distance displaced in front of a screed blade. A reasonable method for evaluating the problem was developed and validated. Several simulations were performed to show the viability of the approach. Future investigations will focus on running various simulations investigating powder particle sizing and screen geometries.

36 MATERIALS SCIENCE↗

DOE OSTI2015

UCLA Final Technical Report for the "Community Petascale Project for Accelerator Science and Simulation”.

The UCLA Plasma Simulation Group is a major partner of the “Community Petascale Project for Accelerator Science and Simulation”. This is the final technical report. We include an overall summary, a list of publications, progress for the most recent year, and individual progress reports for each year. We have made tremendous progress during the three years. SciDAC funds have contributed to the development of a large number of skeleton codes that illustrate how to write PIC codes with a hierarchy of parallelism. These codes cover 2D and 3D as well as electrostatic solvers (which are used in beam dynamics codes and quasi-static codes) and electromagnetic solvers (which are used in plasma based accelerator codes). We also used these ideas to develop a GPU enabled version of OSIRIS. SciDAC funds were also contributed to the development of strategies to eliminate the Numerical Cerenkov Instability (NCI) which is an issue when carrying laser wakefield accelerator (LWFA) simulations in a boosted frame and when quantifying the emittance and energy spread of self-injected electron beams. This work included the development of a new code called UPIC-EMMA which is an FFT based electromagnetic PIC code and to new hybrid algorithms in OSIRIS. A new hybrid (PIC in r-z and gridless in φ) algorithm was implemented into OSIRIS. In this algorithm the fields and current are expanded into azimuthal harmonics and the complex amplitude for each harmonic is calculated separately. The contributions from each harmonic are summed and then used to push the particles. This algorithm permits modeling plasma based acceleration with some 3D effects but with the computational load of an 2D r-z PIC code. We developed a rigorously charge conserving current deposit for this algorithm. Very recently, we made progress in combining the speed up from the quasi-3D algorithm with that from the Lorentz boosted frame. SciDAC funds also contributed to the improvement and speed up of the quasi-static PIC code QuickPIC. We have also used our suite of PIC codes to make scientific discovery. Highlights include supporting FACET experiments which achieved the milestones of showing high beam loading and energy transfer efficiency from a drive electron beam to a witness electron beam and the discovery of a self-loading regime a for high gradient acceleration of a positron beam. Both of these experimental milestones were published in Nature together with supporting QuickPIC simulation results. Simulation results from QuickPIC were used on the cover of Nature in one case. We are also making progress on using highly resolved QuickPIC simulations to show that ion motion may not lead to catastrophic emittance growth for tightly focused electron bunches loaded into nonlinear wakefields. This could mean that fully self-consistent beam loading scenarios are possible. This work remains in progress. OSIRIS simulations were used to discover how 200 MeV electron rings are formed in LWFA experiments, on how to generate electrons that have a series of bunches on nanometer scale, and how to transport electron beams from (into) plasma sections into (from) conventional beam optic sections.

43 PARTICLE ACCELERATORS↗

DOE OSTI2017

Computational Approach for Improving Three-Dimensional Sub-Surface Earth Structure for Regional Earthquake Hazard Simulations in the San Francisco Bay Area

In our Exascale Computing Project (ECP) we seek to simulate earthquake ground motions at much higher frequency than is currently possible. Previous simulations in the SFBA were limited to 0.5-1 Hz or lower (Aagaard et al. 2008, 2010), while we have recently simulated the response to 5 Hz. In order to improve confidence in simulated ground motions, we must accurately represent the three-dimensional (3D) sub-surface material properties that govern seismic wave propagation over a broad region. We are currently focusing on the San Francisco Bay Area (SFBA) with a Cartesian domain of size 120 x 80 x 35 km, but this area will be expanded to cover a larger domain. Currently, the United States Geologic Survey (USGS) has a 3D model of the SFBA for seismic simulations. However, this model suffers from two serious shortcomings relative to our application: 1) it does not fit most of the available low frequency (< 1 Hz) seismic waveforms from moderate (magnitude M 3.5-5.0) earthquakes; and 2) it is represented with much lower resolution than necessary for the high frequency simulations (> 5 Hz) we seek to perform. The current model will serve as a starting model for full waveform tomography based on 3D sensitivity kernels. This report serves as the deliverable for our ECP FY2017 Quarter 4 milestone to FY 2018 “Computational approach to developing model updates”. We summarize the current state of 3D seismic simulations in the SFBA and demonstrate the performance of the USGS 3D model for a few selected paths. We show the available open-source waveform data sets for model updates, based on moderate earthquakes recorded in the region. We present a plan for improving the 3D model utilizing the available data and further development of our SW4 application. We project how the model could be improved and present options for further improvements focused on the shallow geotechnical layers using dense passive recordings of ambient and human-induced noise.

58 GEOSCIENCES↗

DOE OSTI2019

Supporting Simulation Families in Nuclear Reactor Workflows SBIR Phase I. Final Scientific/Technical Report

The US Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) Toolkit has accelerated both the solutions to existing problems and deployment of new designs for current and advanced reactors. Within the toolkit, there exist various families of simulation tools that target specific aspects of a multi-physics simulation for a nuclear reactor. Tools within a specific family vary in both complexity and fidelity as well as the types of resources they require. To complement the NEAMS effort, the Computational Model Builder (CMB) has been designed to provide a graphical user interface to design and execute nuclear reactor workflows. Previously, we had designed a workflow in CMB to design and prepare for execution a neutronics simulation using DIF3D, a diffusion and transport solver developed at Argonne National Laboratory. For this Phase I effort, we extended CMB’s reactor modeling tools to design and execute a workflow using Proteus-SN, a higher-fidelity transport solver developed at Argonne with more preparatory requirements from the workflow. Furthermore, this workflow was designed to maximize the common features of the two workflows, describing the analogous tasks between the workflows as belonging to a “simulation family”. This perspective facilitates the abstraction of the neutronics simulations to a single bifurcated workflow, allowing the engineer to design a reactor and then choose the appropriate simulation according to his or her needs for fidelity and computational cost.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2019

Particle/Ray Interaction Simulation Manager (PRISM) Assessment Report

The Particle/Ray Interaction Simulation Manager (PRISM) is a software tool developed by Triple Ring Technologies (TRT), http://www.tripleringtech.com/, under a Small Business Innovation Research Award from the Department of Homeland Security (DHS) to investigate the utility of simulations in X-ray-based explosive detection systems. The PRISM tool is intended to provide a standardized and intuitive user interface to control and configure realistic simulations of X-ray imaging systems. It includes a command-line interface that enables the user to configure simulated X-ray radiography and computed tomography (CT) systems. It allows the import of computer-aided design (CAD) models of objects and it supports Monte Carlo simulation engines such as GEANT4. PRISM is designed to simulate realistic X-ray projection images including scatter, noise, and various detector effects such as quantum efficiency, lag, and veiling glare. At the request of DHS, the authors assessed PRISM’s capabilities to simulate data generated by a Micro-CT test-bed system (MCT-TB) housed at the Lawrence Livermore National Laboratory (LLNL). This document describes the results of that assessment.

42 ENGINEERING↗

DOE OSTI2020

Cesium Removal Performance Comparisons of Crystalline Silicotitanate Media Batches with Savannah River Site Waste Simulant

The Tank Closure Cesium Removal (TCCR) system uses ion exchange columns filled with Crystalline Silicotitanate (CST) media to process radioactive waste solutions for the removal of Cs137. TCCR currently focuses on dissolving Savannah River Site (SRS) Tank 10H waste (primarily sodium salt cake solids) within the tank followed by at-tank ion exchange column treatment. Plans are underway to prepare and install a second TCCR unit at SRS. Capacity and particle size differences exist between archived (IE-911) and more recently-prepared CST media batches (9120- B and 9140-B). Side-by-side comparison testing was performed to evaluate the cesium removal performance of each batch to aid in selecting the preferred CST batch and media characteristics to load into the second TCCR unit. Batch contact equilibrium and flow-through column tests have been conducted with three CST batches using an SRS Average Simulant. Simulant batch contact equilibrium cesium loading results for the three CST batches (including two lots of one batch) are provided in Table ES-1. The 35 °C data indicates that the archived IE-911 batch has a higher cesium capacity than recently-prepared CST media and that the minor TCCR CST 9120-B media lot (2099000035) has similar cesium removal performance to the major TCCR lot (2099000034). Tests conducted at 25 ºC for the archived IE-911 CST batch indicated lower cesium removal performance with this simulant batch than was observed recently with a different SRS Average Simulant batch. A dilution factor (DF) is typically utilized when modeling engineered CST media cesium loading performance to account for mass contributions from the binder material. In cases where CST performance is lower than expected, this factor includes corrections for the binder and for low performance. The ZAM (Zheng, Anthony, and Miller) Isotherm Model DF values are provided in Table ES-1 for each simulant batch contact result. DF values near 0.5 were determined for the 9120-B and 9140-B CST batch contact tests while DF values near 0.6 (20% higher) were determined for the tests with IE-911. These dilution factors are lower than recently observed with a different SRS Average Simulant batch (9140-B DF = 0.68; IE-911 DF = 1.0). Based on these results, it appears that some component in the simulant solution used for equilibrium testing may have resulted in reduced cesium loading on the CST media.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Demonstration and performance testing of extreme-resolution simulations with static meshes on Summit (CPU & GPU) for a parked-turbine configuration and an actuator-line (mid-fidelity model) wind farm configuration (ECP-Q4 FY2020 Milestone Report)

17 WIND ENERGY↗

DOE OSTI2019

High-order methods for wind turbine simulations and moving to next-generation platforms

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. This milestone represents an effort to increase the fidelity of Nalu-Wind at a fixed resolution through the implementation of a tensor-product based, matrix-free high order scheme. High order finite element methods have increased local work per datum communicated and have the potential to provide significantly more accurate solutions at a fixed number of degrees of freedom. Previous to this milestone, Nalu-Wind had an arbitrary order Control Volume Finite Element Method discretization as a solver option, but it required too much memory and was too slow to be of practical use. The work in this milestone addresses these issues by first implementing an implicit, high order solver that only partially assembles the global system. This reduces the memory footprint of the high-order scheme by orders of magnitude for higher polynomial orders. Second, a faster, tensor-product based method for evaluating the action of the left-hand side was implemented. This reduces the amount of computational work required by the scheme and dramatically enhanced the time-to-solution on example problems. Finally, this milestone is an evaluation of the value of high order methods in the wind application space. With the enhancements to memory and computational cost, accuracy vs. time-to-solution was evaluated for several resolutions on an under-resolved Taylor Green vortex test case. Results show that the high order scheme is cost-competitive with the production low-order schemes in Nalu-Wind, being moderately more expensive than the production edge-based vertex centered finite volume scheme. The evaluation of accuracy on the test case shows a potential benefit to high order at the highest resolution while not deteriorating accuracy on the lowest tested resolution. More work is needed to show value in the wind application, but positive strides have been made.

17 WIND ENERGY↗

DOE OSTI2016

High Fidelity Dynamic Simulation of Cycling in a Natural Gas Combined Cycle Power Plant

The CRADA objective is to use a high-fidelity natural gas combined cycle (NGCC) dynamic simulator deployed at the National Energy Technology Laboratory (NETL) to produce practical engineering ideas that would either reduce the impact on equipment resulting from cycling, or improve certain plant operational modes. The goal is to demonstrate to National Rural Electric Cooperative Association (NRECA) that a high-fidelity dynamic simulator offers a relatively low-cost option for investigation and testing without direct plant experimentation. Initial operation of the simulator revealed the potential for saturation conditions in the final high pressure superheater as the attemperator tried to control temperature at the superheater outlet during gas turbine loading and unloading. Subsequent plant operational data confirmed simulation results. Thus, this particular area was chosen for further investigation. Multiple simulations were performed during loading and unloading of the gas turbine in order to determine operational strategies that prevented saturation and increased the approach to saturation temperature. The results of these simulations are presented in this report. The summary at the end of this report suggests trials that could be performed at the plant to confirm the benefits seen by simulator investigations, such as changes to the attemperator temperature control setpoints and strategic control of the steam turbine inlet pressure control valve.

03 NATURAL GAS↗

DOE OSTI2022

Numerical Simulations of Non-Proliferation Experiments

We have implemented a set of numerical simulations to numerically investigate the outcomes of non-proliferation experiments. The high explosive (HE) in the simulations is in the form of a cylinder. The two-dimensional simulations are carried out in the r-z cylindrical coordinate, and the three-dimensional simulation is in the Cartesian coordinate. From the simulation, the cavity formed after detonation stops its expansion after about 0.15 sec. At the time the cavity stops expanding, the cavity shape seems quite spherical although the initial HE is in the shape of cylinder, but the radius in the vertical direction, which is about 19.7m, is a little less than the one in the horizonal direction, which is about 21.5m. Although there is no layer of air between the HE and tuff in real experiments, in addition to the simulations without air, we add a layer of air in one simulation to assess any possible influence of the air on the shape of the cavity.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

DOE OSTI2022

Update on the Investigation of Commercial Drying Cycles Using the Advanced Drying Cycle Simulator

The purpose of this report is to document updates on the apparatus to simulate commercial vacuum drying procedures at the Nuclear Energy Work Complex at Sandia National Laboratories. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system during subsequent storage and disposal. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates well-designed investigations of drying process efficacy and water retention. Scaled tests that incorporate relevant physics and well-controlled boundary conditions are essential to provide insight and guidance to the simulation of prototypic systems undergoing drying processes. This report documents a new test apparatus, the Advanced Drying Cycle Simulator (ADCS). This apparatus was built to simulate commercial drying procedures and quantify the amount of residual water remaining in a pressurized water reactor (PWR) fuel assembly after drying. The ADCS was constructed with a prototypic 17×17 PWR fuel skeleton and waterproof heater rods to simulate decay heat. These waterproof heaters are the next generation design to heater rods developed and tested at Sandia National Laboratories in FY20. This report describes the ADCS vessel build that was completed late in FY22, including the receipt of the prototypic length waterproof heater rods and construction of the fuel basket and the pressure vessel components. In addition, installations of thermocouples, emissivity coupons, pressure and vacuum lines, pressure transducers, and electrical connections were completed. Preliminary power functionality testing was conducted to demonstrate the capabilities of the ADCS. In FY23, a test plan for the ADCS will be developed to implement a drying procedure based on measurements from the process used for the High Burnup Demonstration Project. While applying power to the simulated fuel rods, this procedure is expected to consist of filling the ADCS vessel with water, draining the water with applied pressure and multiple helium blowdowns, evacuating additional water with a vacuum drying sequence at successively lower pressures, and backfilling the vessel with helium. Additional investigations are expected to feature failed fuel rod simulators with engineered cladding defects and guide tubes with obstructed dashpots to challenge the drying system with multiple water retention sites.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗

DOE OSTI2022

Vitrification of Hanford Tank 241-AP-101 Waste and Simulant

Hanford tank 241-AP-101 (referred to herein as AP-101) is the second Hanford radioactive tank waste planned to be processed and vitrified. A simulant version of AP 101 waste was formulated from the best-basis inventory (BBI) for the Hanford Tank 241-AP-101 liquid with an assumed target dilution of the waste from the BBI sodium molarity of 8.61 M to the desired 5.5 M Na. After the addition of glass-forming chemicals (GFCs), the simulant melter feed was processed in a non-radioactive, continuous laboratory-scale melter (CLSM) system. The AP-101 simulant melter feed was charged into the CLSM for 6.11 h of processing, which produced 6.55 kg of glass, for an average glass production rate of 2275 kg m 2 d -1 . Since there were no processing issues with the AP-101 simulant melter feed, AP-101 melter feed made with actual waste was then processed in a CLSM system built into a contamination area in a radioactive environment. The melting behavior characteristics appeared similar for both the simulant and waste melter feeds. The AP-101 waste melter feed was charged into the CLSM for 12.14 h of processing, which produced 8.75 kg of glass, for an average glass production rate of 1530 kg m 2 d -1 . During the AP-101 waste melter feed charging, the pump used to move the feed reached a maximum and it is believed that if the pump had a greater capacity, a greater average glass production rate could have been achieved. A constituent of interest present in low quantities in the AP-101 waste is 99 Tc or its non-radioactive surrogate, Re, added to the AP-101 simulant. Analysis for the quantities of 99 Tc and Re in the AP-101 glass product resulted in an average single-pass retention from the melter feed during relative chemical steady state of 55 ± 2 % for 99 Tc and 45 ± 2 % for Re. Compared to the processing of other melter feeds, the retention of 99 Tc in the AP-101 glass was greater than in both AP-107 and AP-105 glass, while the retention of Re in the AP-101 was less than in the AP-107 glass, but greater than in the AP-105 glass. A spike of I was added into the AP-101 melter feed that could be detected above the analysis detection limits. However, the iodine was only detectable above the ~6 ppm limit in one glass pour: the pour immediately following the burn off of the cold cap, where the I level reached ~30 ppm. This event was significant because the glass was poured immediately after burn off and thus it is presumed that the iodine had yet to volatilize from the glass melt while idling. It is recommended to perform future tests with I spikes at greater levels so that it can be detected in additional glass pours to determine if the expected 50 % retention of I used in the Kim et al. glass models can be confirmed. Offgas liquid samples were analyzed for acetonitrile, which was present at greater concentrations in CLSM liquids than in other scaled melter systems. This result was expected based on unique conditions with the CLSM system including a small plenum space leading to low residence time for offgas and the rapidity of offgas cooling upon exiting the CLSM vessel due to the location and environment. About 90 % of the total acetonitrile captured during both the AP-101 simulant and waste CLSM runs was found in the offgas condensate and demister liquids, thus it is recommended that only those liquids be sent for analysis if future testing to study the presence of acetonitrile in offgas products is desired.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Mesoscale simulations of pressure-shear loading of granular tungsten carbide

Numerical simulations of pressure-shear loading of a granular material are performed using the shock physics code CTH. A simple mesoscale model for the granular material is used that consists of a randomly packed arrangement of solid circular or spherical grains of uniform size separated by vacuum. The grain material is described by a simple shock equation of state, elastic perfectly plastic strength model, and fracture model with baseline parameters for WC taken from previous mesoscale modeling work. Simulations using the baseline material parameters are performed at the same initial conditions of pressure-shear experiments on dry WC powders. Except for some localized flow regions appearing in simulations with an approximate treatment of sliding interfaces among grains, the samples respond elastically during shear, which is in contrast to experimental observations. By extending the simulations to higher shear wave amplitudes, macroscopic shear failure of the simulated samples is observed with the shear strength increasing with increasing stress confinement. The shear strength is also found to be strongly dependent on the grain interface treatment and on the fracture stress of the grains, though the variation in shear strength due to fracture stress decreases with increasing stress confinement. At partial compactions, the transverse velocity histories show strain-hardening behavior followed by formation of a shear interface that extends through the transverse dimensions of the sample. Near full compaction, no strain hardening is observed and, instead, the sample transitions sharply from an elastic response to formation of an internal shear interface. Agreement with experiment is shown to worsen with increasing confinement stress with simulations overpredicting the shear strengths measured in experiment. The source of the disagreement can be ultimately attributed to the Eulerian nature of the simulations, which do not treat contact and fracture realistically.

36 MATERIALS SCIENCE↗

DOE OSTI2023

RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models

This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.

42 ENGINEERING↗

DOE OSTI2024

Establishing capabilities for quantum computing and simulations for energy applications

Quantum information science (QIS) is creating potential transformative opportunities to exploit intricate quantum mechanical phenomena in new ways for obtaining and processing information to advance many areas of science and engineering. Since the National Quantum Initiative Act was signed into law in 2018, developing QIS capability and competency is one of the most urgent tasks of DOE to make sure the US win the quantum race. The QIS contains four pillars: quantum computing, quantum simulations, quantum sensing, and quantum networking. To apply QIS in energy related applications, the key is to develop the capability of quantum computing & simulation tools. In this project, we propose to develop the capability of quantum computing and simulations at NETL to target fossil energy related problems. We will install a simulator (e. g. IBM qiskit) on the NETL supercomputer to simulate the environments of quantum computer. Based on the current available quantum algorithms for quantum chemistry, we will develop quantum computing codes to perform simulation which will focus on fossil energy sector challenges, including CO2 capture & conversion, sensing, and fuel conversion. We then will seek opportunities to run our codes on real quantum computers (such as IBM-Q, google Sycamore, etc.). Through this project, the new capability of quantum computing and simulations will be established at NETL. In addition, the NETL workforce in this area will be trained ready to conduct more complicated tasks in line with NETL missions to enhance the nation’s energy foundation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

New Muon Campus Simulations for the Muon $g\textrm{-}2$ Experiment at Fermilab

The Muon $g\textrm{-}2$ Experiment (E989) at Fermilab conducted high-precision measurements of the muon anomalous magnetic moment $a_{\mu}$ using a storage ring from 2018 to 2023, achieving a remarkable precision of $200\:\mathrm{ppb}$ over Runs 1-3, with analyses for Runs 4-6 ongoing. A comprehensive understanding of the storage ring's beam dynamics and its accurate simulations are crucial for achieving the experiment's ambitious goals. One of the requirements for this effort is a very detailed knowledge of the phase space distribution of the beam. To address this requirement, we performed high-statistics simulations of the Muon $g\textrm{-}2$ Target Station (AP0) and the Muon Campus beamlines: M2 and M3, followed by the Delivery Ring, and then M4 and M5. The resulting muon distribution at the end of the M5 beamline from our previous $3\times10^{12}$ protons-on-target (PoT) simulation serves as an essential input for the storage ring simulations. In 2024, to facilitate the analyses of Runs 4-6, we have updated our Muon Campus models and re-optimised certain parameters to reflect the operational currents and wire chamber measurements of the beam. For these optimisations, we employed the heterogeneous island method, implemented in our evolutionary optimisation tool, glyfada. This key update addresses the need to use the best possible beam for the storage ring simulations and thus supports the experiment's overall precision. The Muon $g\textrm{-}2$ Experiment at Fermilab has successfully achieved its $70\:\mathrm{ppb}$ systematic uncertainty goal and collected 21 times more data than its predecessor at BNL. The updated and improved Muon Campus models and simulations not only facilitate the experiment's efforts to potentially resolve the current tension between experimental measurements and theoretical predictions of $a_{\mu}$, but also provide a basis for future simulations for the Mu2e Experiment (E-973), which utilises shared Muon Campus beamlines.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Simulating the CMS High Granularity Calorimeter with ML

Detector simulation is a key component of physics analysis and related activities in CMS. In the upcoming High Luminosity LHC era, simulation will be required to use a smaller fraction of computing in order to satisfy resource constraints. At the same time, CMS will be upgraded with the new High Granularity Calorimeter (HGCal), which requires significantly more resources to simulate than the existing CMS calorimeters. This computing challenge motivates the use of generative machine learning models as surrogates to replace full physics-based simulation. We study the application of state-of-the-art diffusion models to simulate particle showers in the CMS HGCal. We will discuss methods to overcome the challenges posed by the high-dimensional, irregular geometry of the HGCal. The quality of the showers produced by the diffusion model will be assessed by comparison to the full GEANT4-based simulation. The increase in simulation throughput will be quantified and methods to accelerate the diffusion model inference will also be discussed.

Amram, Oz↗