Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 775 records · Page 43

The impact of aerosol mixing state on immersion freezing: insights from classical nucleation theory and particle-resolved simulations

Immersion freezing, initiated by ice-nucleating particles (INPs) in supercooled aqueous droplets, plays an important role in the formation of ice crystals within clouds. The efficiency of immersion freezing depends strongly on INP composition and, crucially, on the mixing state – how chemical species are distributed across the particle population. Here, we quantify the impact of aerosol mixing state on immersion freezing using a combined theoretical and particle-resolved modeling approach. We derive analytical expressions for the frozen fraction of internally and externally mixed INP populations based on classical nucleation theory, showing that the frozen fraction is sensitive to whether ice-active species are present in all particles or only in a subset of the population. We introduce a multi-species immersion freezing scheme into the particle-resolved model PartMC, using the water activity-based immersion freezing model (ABIFM) to compute freezing probabilities for mixed-composition particles. To improve computational efficiency, we implement a Binned Tau-Leaping algorithm and demonstrate an order-of-magnitude speedup with minimal accuracy loss. Simulations reproduce the analytical trends in limiting cases and extend the analysis to more general aerosol populations, where mixing state continues to exert a substantial control on frozen fraction. Sensitivity analyses across particle size, species type, and cooling condition reveal that the mixing state effect is most pronounced when small amounts of highly efficient INPs are mixed with less efficient materials. These findings underscore the need to represent aerosol mixing state explicitly in models of heterogeneous ice nucleation to reduce uncertainty in cloud-phase partitioning.

54 ENVIRONMENTAL SCIENCES↗

Effective Lanthanide Diffusivity through U-Zr Metallic Fuel

Fuel-cladding chemical interaction (FCCI) is one of the main concerns for the performance of U-Zr metallic fuels. Although FCCI results in complex phase transformations and morphological changes in both the fuel and cladding sides of the interface, the formation of a brittle layer in the cladding side of the interface (sometimes referred to as wastage) has one of the most significant impacts on the cladding mechanical integrity. This brittle layer is associated with the formation of intermetallic compounds between the cladding constituents Fe, Cr and the lanthanide fission product species. The formation of this cladding wastage layer has been the focus of recent efforts to develop a mechanistic modeling framework, including production of lanthanides in the fuel, transport through the fuel to the fuel-cladding interface, and phase transformation to form intermetallic phases. To calculate the rate of transport of lanthanides through the fuel, a multi-scale computational approach has been used. Atomistic calculations have been performed to determine the diffusivities of Nd, the most prevalent lanthanide species, through the solid fuel matrix and along pore surfaces. These diffusivities have been used as input parameters for a mesoscale model to determine an effective diffusion coefficient that accounts for the formation of porosity, its interconnection, and infiltration with bond sodium from the fuel-cladding gap. The effective diffusivity has been used in engineering-scale simulations using the BISON fuel performance code, and BISON simulations have been compared to representative validation cases from the EBR-II reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An Insight-Centric Paradigm for Data Reduction and Inference Speed Improvement at the Scurry Area Canyon Reef Operator’s Committee (SACROC) Unit

The poster presents the work conducted under SMART focusing on using insight-centric approach to design a meaningful proxy for machine learning. Domain insights are critical not just in understanding the prediction results but also in designing the model. This study demonstrated that a single meaningful scaler (as an extreme case) can effectively replace full-size 3D geologic properties. The model's accuracies are on par with other models, and it is the fastest model to predict all test cases, 5,000 times faster than traditional simulations.

Shih, Chung Yan↗

An Agent-Based Modeling Approach for Spatiotemporal Optimization of Electric Vehicle Fast-Charging Station Demand

With increasing electric vehicle (EV) adoption, managing public fast-charging demand effectively is crucial to avoid grid strain. This study investigates the potential of using dynamic pricing schemes to address this challenge. Presented in this study is a scalable agent-based simulation framework, which is applied to a case study in Richmond, Virginia, that assumes a 50% EV adoption rate in 2040. Two pricing schemes are compared: (1) a dynamic-pricing scheme based on station utilization and (2) a dynamic-pricing scheme based on peak power at the station. These schemes are compared to two baseline scenarios: (1) unscheduled first-come, first-served and (2) scheduled with constant price. The study’s results suggest that dynamic pricing has the potential to influence EV charging behavior, inducing both spatial and temporal shifts, but does so at the cost of inducing inconvenience to EV drivers. The results suggest the peak-power dynamic pricing scheme has the potential to mitigate peak demand pressures on the grid with minimal inconvenience, offering a promising approach for sustainable EV charging infrastructure expansion.

33 - ADVANCED PROPULSION SYSTEMS↗

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗

Edge-localized mode mitigation enabled by active control of pedestal density gradient with new EAST tokamak divertor

Mitigation of large edge-localized modes (ELMs) has been achieved by actively reducing the pedestal density gradient with the EAST new right-angled lower divertor through changing the strike point position from the vertical target to the horizontal target. A series of dedicated experiments in the 2021–2024 EAST campaigns demonstrate that this ELM control solution is highly reproducible in a broad parameter space of edge safety factor q 95 = 4.7–7.1, heating power P total = 2.3–5 MW, and pedestal collisionality $ν_{e,ped}^{*}$ = 1–6, under both favorable and unfavorable magnetic configurations. Higher plasma density could facilitate the achievement of this ELM control solution. Statistical results indicate that the ELM mitigation effect can be observed at relatively larger Greenwald density fraction of f GW > 0.47. In addition, this ELM mitigation effect can be achieved with both lithium-coated and boronized metal walls. The pedestal density gradient is systematically lower in the horizontal target case than that of the vertical target case when the ELM mitigation effect can be observed. SOLPS-ITER simulation results indicate that the pedestal fueling from divertor recycling is significantly lower in the horizontal target case. This could contribute to the formation of a flattened pedestal density profile with small ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation budgeting for hybrid effective field theories

In this work, we forecast the number of, and requirements on, N-body simulations needed to train hybrid effective field theory (HEFT) emulators for a range of use cases, using a hybrid of HMcode and perturbation theory as a surrogate model. Our accuracy goals, determined with careful consideration of statistical and systematic uncertainties, are 1% accurate in the high-likelihood range of cosmological parameters, and 2% accurate over a broader parameter space volume for k < 1 h Mpc -1 and z < 3. Focusing in part on the 8-parameter w 0 w a CDM+ m ν cosmological model, we find that < 225 simulations are required to meet our error goals over our wide parameter space, including models with rapidly evolving dark energy, given our simulation and emulator recommendations. For a more restricted parameter space volume, as few as 80 simulations are sufficient. We additionally present simulation forecasts for example use cases, and make the code used in our analyses publicly available. These results offer practical guidance for efficient emulator design and simulation budgeting in future cosmological analyses.

cosmological parameters from LSS↗

Rare Earth and Transition Metal Containing Glasses

Transition metal (TM) and rare earth (RE) ions have been incorporated into many glass systems such as silicate, phosphate, and borosilicate-based oxide glasses, as well as in halide and chalcogenide glasses, that find applications ranging from optical, photonic, and magnetic devices, solid-state battery, to nuclear waste disposal. Understanding the structural role of RE and TM in these glasses can help to develop glass compositions for targeted applications with either high-optical emission efficiency, electrical conductivity, or chemical durability. In this chapter, we first provide a general introduction of the applications and structural features of RE and TM in glasses, then the critical aspects of molecular dynamics (MD) simulations of these glasses such as interatomic potentials, structural analysis tools to study RE and TM ions in glasses and their clustering behaviors, Quantitative Structure–Property Analysis (QSPR), diffusion and dynamic property calculations, and electronic structure calculations to understand electronic defects such as charge trapping and radiation effects are introduced. Three representative case studies are presented: the first one is on MD simulations of erbium- and europium-doped silica and silicate glasses, as well as cerium doped aluminophosphate glasses, that revealed the effect of glass composition on RE ion local structure and clustering behavior. Electronic structure calculations of cerium-doped glass show how the existence of multioxidation states help to mediate radiation-induced damages caused by excited electron–hole pairs was also discussed. The second one focuses on alkali vanadophosphate glasses where the existence of two vanadium oxidation states help to provide electronic conduction in the glasses while alkali ions provide ionic conduction. MD simulations were used to understand vanadium environments and other structural aspects in the phosphate glasses, as well as the ionic transport behaviors of alkali ions. The third case study is on zirconium-containing borosilicate and aluminosilicate glasses that find wide applications in nuclear waste disposal. MD simulations help to provide structural details of zirconium ions that are validated by diffraction and EXAFS spectra. The structural information was used to interpret changes of mechanical properties and chemical durability by using QSPR and other analyses-based MD-generated structure models.

Du, Jincheng↗

High-level fuel fabrication facility designs from discrete-event simulation

Like other industrial processes, the production of metallic nuclear fuels (MNF) requires that fabrication facilities be able to reliably meet production demands, operate efficiently, and adhere to federal and local safety regulations. In turn, the set of design variables employed by such a facility, such as operations policies, infrastructure and machinery purchased, and the type and number of staff hired, directly impact a facility’s ability to satisfy these goals. Therefore, facility designers must carefully determine which set of design variable values optimally satisfies these constraints. In this paper, we explore how values for these high-level design variables, namely hiring requirements, can be determined in the context of nuclear fuel manufacturing through the coupling of physics-based and discrete-event simulation technologies. Using the Versatile Test Reactor (VTR) program as a case study, we demonstrate how SCALE and MCNP nuclear physics model outputs can be integrated into ExtendSim discrete-event simulation (DES) models of the fuel fabrication process to determine the optimal number of staff hired to ensure fuel production goals are met, operations comply with effective dose limit regulations, and overall project costs are reduced.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Wavelet-based local mesh refinement for rainfall–runoff simulations

A wavelet-based local mesh refinement (wLMR) strategy is designed to generate multiresolution and unstructured triangular meshes from real digital elevation model (DEM) data for efficient hydrological simulations at the catchment scale. The wLMR strategy is studied considering slope- and curvature-based refinement criteria to analyze DEM inputs: the slope-based criterion uses bed elevation data as input to the wLMR strategy, whereas the curvature-based criterion feeds the bed slope data into it. The performance of the wLMR meshes generated by these two criteria is compared for hydrological simulations; first, using three analytical tests with the systematic variation in topography types and then by reproducing laboratory- and real-scale case studies. The bed elevation on the wLMR meshes and their simulation results are compared relative to those achieved on the finest uniform mesh. Repeated analytical tests show that the slope- and curvature-based criteria are equally effective with the wLMR strategy, and that it is easier to decide which criterion to take in relation to the (regular) shape of the topography. For the realistic case studies: (i) slope analysis provides a better metric to assess the correlation of a wLMR mesh to the fine uniform mesh and (ii) both criteria predict outlet hydrographs with a close predictive accuracy to that on the uniform mesh, but the curvature-based criterion is found to slightly better capture the channeling patterns of real DEM data.

54 ENVIRONMENTAL SCIENCES↗

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

tih_vqe [SWR-23-32]

This software supports the paper, "Exploring the scaling limitations of the variational quantum eigensolver with the bond dissociation of hydride diatomic molecules," published in the International Journal of Quantum Chemistry, whose abstract is as follows: Materials simulations involving strongly correlated electrons pose fundamental challenges to state-of-the-art electronic structure methods but are hypothesized to be the ideal use case for quantum computing. To date, no quantum computer has simulated a molecule of a size and complexity relevant to real-world applications, despite the fact that the variational quantum eigensolver (VQE) algorithm can predict chemically accurate total energies. Nevertheless, because of the many applications of moderately-sized, strongly correlated systems, such as molecular catalysts, the successful use of the VQE stands as an important waypoint in the advancement toward useful chemical modeling on near-term quantum processors. In this paper, we take a significant step in this direction. We lay out the steps, write, and run parallel code for an (emulated) quantum computer to compute the bond dissociation curves of the TiH, LiH, NaH, and KH diatomic hydride molecules using VQE. TiH was chosen as a relatively simple chemical system that incorporates d orbitals and strong electron correlation. Because current VQE implementations on existing quantum hardware are limited by qubit error rates, the number of qubits available, and the allowable gate depth, recent studies have focused on chemical systems involving s and p block elements. Through VQE + UCCSD calculations of TiH, we evaluate the near-term feasibility of modeling a molecule with d-orbitals on real quantum hardware. We demonstrate that the inclusion of d-orbitals and the use of the UCCSD ansatz, which are both necessary to capture the correct TiH physics, dramatically increase the cost of this problem. We estimate the approximate error rates necessary to model TiH on current quantum computing hardware using VQE+UCCSD and show them to likely be prohibitive until significant improvements in hardware and error correction algorithms are available.

Graf, Peter↗

Cost effectiveness of preemptive school closures to mitigate pandemic influenza outbreaks of differing severity in the United States

Background: Nonpharmaceutical interventions (NPIs) may be considered as part of national pandemic preparedness as a first line defense against influenza pandemics. Preemptive school closures (PSCs) are an NPI reserved for severe pandemics and are highly effective in slowing influenza spread but have unintended consequences. Methods: We used results of simulated PSC impacts for a 1957-like pandemic (i.e., an influenza pandemic with a high case fatality rate) to estimate population health impacts and quantify PSC costs at the national level using three geographical scales, four closure durations, and three dismissal decision criteria (i.e., the number of cases detected to trigger closures). At the Chicago regional level, we also used results from simulated 1957-like, 1968-like, and 2009-like pandemics. Our net estimated economic impacts resulted from educational productivity costs plus loss of income associated with providing childcare during closures after netting out productivity gains from averted influenza illness based on the number of cases and deaths for each mitigation strategy. Results: For the 1957-like, national-level model, estimated net PSC costs and averted cases ranged from $\$7.5$ billion (2016 USD) averting 14.5 million cases for two-week, community-level closures to $\$97$ billion averting 47 million cases for 12-week, county-level closures. We found that 2-week school-by-school PSCs had the lowest cost per discounted life-year gained compared to county-wide or school district–wide closures for both the national and Chicago regional-level analyses of all pandemics. The feasibility of spatiotemporally precise triggering is questionable for most locales. Theoretically, this would be an attractive early option to allow more time to assess transmissibility and severity of a novel influenza virus. However, we also found that county-wide PSCs of longer durations (8 to 12 weeks) could avert the most cases (31–47 million) and deaths (105,000–156,000); however, the net cost would be considerably greater ($\$88$-$\$103$ billion net of averted illness costs) for the national-level, 1957-like analysis. Conclusions: We found that the net costs per death averted ($\$180,000$-$\$4.2$ million) for the national-level, 1957-like scenarios were generally less than the range of values recommended for regulatory impact analyses ($\$4.6$ to 15.0 million). This suggests that the economic benefits of national-level PSC strategies could exceed the costs of these interventions during future pandemics with highly transmissible strains with high case fatality rates. In contrast, the PSC outcomes for regional models of the 1968-like and 2009-like pandemics were less likely to be cost effective; more targeted and shorter duration closures would be recommended for these pandemics.

60 APPLIED LIFE SCIENCES↗

Large-Scale Forcing Impact on the Development of Shallow Convective Clouds Revealed From LASSO Large-Eddy Simulations

Real-world large-eddy simulations (LES) are driven by time-varying large-scale forcings (LSF) - e.g., temperature advection, moisture advection, and subsidence - derived from large-scale weather models. This study investigates the impact of the uncertainty in LSF on real-world LES in terms of the development of shallow convection at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) atmospheric observatory for the 11 June 2016 case using LES provided by the U.S. Department of Energy's LES ARM Symbiotic Simulation and Observation (LASSO) activity. The LASSO dataset provides an ensemble of LES for the selected case, which consists of LES runs that were driven by different LSF. The two contrasting LES runs investigated here generate different types of convective clouds, i.e., nonprecipitating shallow clouds and precipitating cumulus congestus, mainly due to the difference of LSF in temperature advection in the free troposphere. The temperature advection modulates the strength of the capping inversion and therefore the buoyancy of the air parcels rising from the atmospheric boundary layer (ABL). The inversion, together with large-scale updrafts, controls the penetration of the ABL thermals into the free troposphere, leading to cumulus congestus in the case of a weaker inversion. In contrast, clouds remain shallow in the case of a strong inversion. Differences between the two simulations are amplified over time, as mixed-phase clouds are formed near the top of the congestus in the weaker inversion case. Furthermore, this high dependency of LES results to LSF stresses the importance of accurate LSF by large-scale models to real-world LES simulations.

54 ENVIRONMENTAL SCIENCES↗

An embedded 3D fracture modeling approach for simulating fracture-dominated fluid flow and heat transfer in geothermal reservoirs

In this paper, we describe an efficient modeling approach, named embedded discrete fracture method (EDFM), for incorporating arbitrary 3D, discrete fractures, such as hydraulic fractures or faults, into modeling fracture-dominated fluid flow and heat transfer in fractured geothermal reservoirs. This technique allows 3D discrete fractures to be discretized independently from surrounding rock volume and inserted explicitly into a primary fracture/matrix grid, generated without including 3D discrete fractures in prior. An effective computational algorithm is developed to discretize these 3D discrete fractures and construct local connections between 3D fractures and fracture/matrix grid blocks representing the surrounding rock volume. The constructed gridding information on 3D fractures is then added to the primary grid. This embedded fracture modeling approach can be directly implemented into a developed geothermal reservoir simulator via the integral finite difference (IFD) method or with TOUGH2 technology. This embedded fracture modeling approach is very promising and computationally efficient to handle realistic 3D discrete fractures with complicated geometries, connections, and spatial distributions. Compared with other fracture modeling approaches, it avoids cumbersome 3D unstructured, local refining procedures, and increases computational efficiency by simplifying Jacobian matrix size and sparsity, while maintaining enough accuracy. Several numeral simulations are presented to demonstrate the utility and robustness of the proposed technique. Our numerical experiments show that this approach captures all the key patterns about fluid flow and heat transfer dominated by fractures in these cases. Thus, this approach is readily available to the simulation of fractured geothermal reservoirs with both artificial and natural fractures.

15 GEOTHERMAL ENERGY↗

Phase Field Modeling of Crystallographic Corrosion Pits

The modeling of corrosion to understand and predict corrosion behavior is a topical issue. Here, a 3D phase field model is developed to simulate the pit morphology, primarily focusing on crystallographic pits. A crystallographic function is employed to incorporate different corrosion rates for different crystallographic planes. The model is benchmarked and validated against an analytical solution for a simple case. 3D crystallographic pits in a single crystal are simulated and the effect of substrate orientation on the pit morphology is studied. The crystallographic pit morphology changes significantly with the substrate orientation and these morphologies have a symmetry consistent with the substrate orientation. This first 3D phase field model of crystallographic pits will help in predicting the intricate shapes of pits, thereby, pushing the frontiers of pitting corrosion modeling.

36 MATERIALS SCIENCE↗

An Open-Source Virtual Testbed for a Real Net-Zero Energy Community

Net zero energy communities (NZECs) are critical to ensure sustainability and resilience of modernized power systems. System modeling helps overcome technical challenges in designing and operating NZECs. In this paper, we present the modeling work based on a real NZEC. Two sets of models are developed: higher-fidelity physics-based models considering the interaction between subsystems of the studied NZEC and capturing fast-dynamics; and lower-fidelity data-driven models requiring less resource to establish and/or run. All models are validated against measurements from this real NZEC. In addition, we create a simulation framework which streamlines the processes for simulation and thus allows using developed models to form a virtual testbed. To demonstrate the usage of the virtual testbed, a case study is conducted where a building-to-grid integration control is evaluated via simulation. The evaluation results suggest that the tested control significantly smooths the power draw of the studied community and doesn’t sacrifice the thermal comfort to a great extent.

Huang, Sen↗

Modeling Environmental Effects on Ventilated Spent Fuel Storage Systems

This report describes newly developed external environment wind effects models of spent nuclear fuel (SNF) dry storage systems. The primary purpose of these wind effects models is to better understand particle deposition on SNF canisters in the context of chloride-induced stress corrosion cracking. The goal of this effort is to further the understanding and improve the wind effects model of the Nuclear Horizontal Modular Storage (NUHOMS ® ) Advanced Horizontal Storage Module–High Seismic (AHSM-HS) storage system in support of the planned Canister Deposition Field Demonstration (CDFD) project (Durbin et al. 2021, Suffield et al. 2021). The steps in this wind effects investigation begin with a validation case of comparing experimental data with a STAR-CCM+ computational fluid dynamics (CFD) model of the Belowground Vertical Dry Cask Simulator (BVDCS) with the external environment explicitly modeled. Next, a test case is performed with the AHSM-HS models comparing solar-loading strategies for a standalone and wind effects model in STAR-CCM+ and a standalone model in ANSYS Mechanical Parametric Design Language (APDL). The final portion of this study compares results of the validation exercise with applications of wind effects models for two horizontal storage systems, a site specific NUHOMS ® horizontal storage module and a NUHOMS ® AHSM-HS.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗