Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Fluid Dynamic (CFD)”

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 433 records · Page 24

Analytical Study and Experimental Verification of Shear-Thinning Ink Flow in Direct Ink Writing Process

Abstract Direct ink writing (DIW) process is a facile additive manufacturing technology to fabricate three-dimensional (3D) objects with various materials. Its versatility has attracted considerable interest in academia and industry in recent years. As such, upsurging endeavors are invested in advancing the ink flow behaviors in order to optimize the process resolution and the printing quality. However, so far, the physical phenomena during the DIW process are not revealed in detail, leaving a research gap between the physical experiments and its underlying theories. Here, we present a comprehensive analytical study of non-Newtonian ink flow behavior during the DIW process. Different syringe-nozzle geometries are modeled for the comparative case studies. By using the computational fluid dynamics (CFD) simulation method, we reveal the shear-thinning property during the ink extrusion process. Besides, we study the viscosity, shear stress, and velocity fields, and analyze the advantages and drawbacks of each syringe-nozzle model. On the basis of these investigations and analyses, we propose an improved syringe-nozzle geometry for stable extrusion and high printing quality. A set of DIW printing experiments and rheological characterizations are carried out to verify the simulation studies. The results developed in this work offer an in-depth understanding of the ink flow behavior in the DIW process, providing valuable guidelines for optimizing the physical DIW configuration toward high-resolution printing and, consequently, improving the performance of DIW-printed objects.

Engineering↗

Large Eddy Simulation of Convective Heat Transfer in a Random Pebble Bed Using the Spectral Element Method

The development of fluoride-cooled high-temperature reactors has drastically increased the demand for an in-depth understanding of the heat transfer (HT) in packed beds cooled by liquid salts. The complex flow fields and space-dependent porosity found in a pebble bed require a detailed understanding to ensure the proper cooling of the reactor core during normal and accident conditions. As detailed experimental data are complicated to obtain for these configurations, high-fidelity simulation such as large eddy simulation and direct numerical simulation (DNS) can be employed to create a high-resolution heat transfer numerical database that can assist in addressing industrial-driven issues associated with the heat transfer behavior of fluoride-cooled high-temperature reactors. In this paper, we performed a series of large eddy simulation using computational fluid dynamics (CFD) code NekRS to investigate the heat transfer for a bed of 1741 pebbles. Further, the characteristics of the flow, such as average, rms, and time series of velocity and temperature, have been analyzed. Porous media averages have also been performed. The simulation results show a good agreement between non-conjugate heat transfer and conjugate heat transfer. The generated data will be used to benchmark heat transfer modeling methods and local maxima/minima of heat transfer parameters. It will also be used for supporting convective heat transfer quantification for Kairos Power and benchmarking lower fidelity models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Helical Fins for Concentrated Solar Receivers: Design Optimization and Entropy Analysis

Concentrated solar power (CSP) with thermal energy storage (TES) has the potential to achieve grid parity. This can be realized by operating CSP systems at temperatures above 700 °C with high-efficiency sCO 2 power cycles. However, operating CSP systems at such temperatures poses several challenges, among which the design of solar receivers to accommodate increased thermal loads is critical. To this end, this work explores and optimizes various swirl-inducing internal fin designs for solar receiver tubes. These fin designs not only improve the thermal performance of receiver tubes but also levelize temperature unevenness caused by non-uniform thermal loading. In this work, the geometric parameters of the fin designs are optimized to maximize the Nusselt number with a constraint on the friction factor. This optimization, however, is computationally intensive, requiring hundreds of simulation calls to computational fluid dynamics (CFD) models. Here, to circumvent this problem, surrogate models are used to approximate the simulation outputs needed during the optimization. In addition, this study also examines the fin designs from an entropy generation perspective. To this end, the entropy contributions from thermal and viscous effects are quantitatively compared while varying the operational Reynolds number.

14 SOLAR ENERGY↗

Identification of Dynamic Force Coefficients for an Additively Manufactured Hermetic Squeeze Film Bearing Support Damper Utilizing a Pass-Through Channel

Abstract The following paper presents breakthrough experimental results for a new hermetic squeeze film damper (HSFD) concept that is integrally designed within an externally pressurized tilting-pad radial gas bearing support. The flexibly damped gas bearing module was designed for a 7.2″ (183 mm) diameter shaft and fabricated using direct metal laser melting (DMLM); also known as additive manufacturing. The bearing and HSFD were sized based on ongoing studies for oil-free supercritical carbon dioxide (sCO2) power turbines in the 8.5 MW–10 MW power range. The development of the new damper concept was motivated by past dynamic testing on HSFD, which generated frequency-dependent stiffness and damping force coefficients. In efforts to eliminate the frequency dependency, a new HSFD architecture was conceived that adds accumulator volumes and a pass-through channel to previously conceived HSFD flow network designs. The other motivation for the work is the need to develop a cost-effective and reliable oil-free bearing technology that is scalable to large power turbomachinery applications. There were several objectives for the following work. The first objective was to successfully design and fabricate a single piece bearing-damper using additive manufacturing, while dimensionally controlling critical design features. The paper discusses the manufacturing steps and shows cut-ups that reveal adequate clearance control capability with internal damper clearances. The second objective was to perform experimental testing with the new HSFD design in efforts to extract stiffness and damping coefficients for excitation frequencies within 20–160 Hz and peak vibration amplitudes between 0.25 mils (6.35 microns) to 1 mil (25.4 microns). The test results for a single HSFD bearing module indicated that the design modifications to the HSFD architecture were successful in eliminating nearly all the frequency dependencies for the stiffness force coefficient. The dynamic tests yielded a stiffness coefficient that varied between 112 klb/in. (19.6 MN/m) and 96 klb/in. (16.8 MN/m). The damping force coefficient however, exhibited relatively more variation with frequency with values residing between 175 lb-s/in. (31 kN-s/m) to 214 lb-s/in. (37 kN-s/m). Finally, the paper advances a three-dimensional fluid-structure interaction (FSI) model using transient finite element analysis (FEA) coupled to a computational fluid dynamics (CFD) model. The FSI analysis performed between 20 Hz and 80 Hz was used to predict the stiffness and damping of the HSFD using a quarter-section model of the damper. The FSI analysis was able to support test results by showing only a 6–7.4% change in the magnitude of force coefficients. Stiffness predictions agree reasonably well with experiments whereas damping is underpredicted.

Engineering↗

Characterization of Flex-Fuel Prechamber Enabled Mixing-Controlled Combustion With Gasoline/Ethanol Blends at High Load

This numerical study focuses on the characterization of prechamber-enabled mixing-controlled combustion (PC-MCC) at ~18 bar brake mean effective pressure (BMEP) and 2200 rpm with 10% by volume ethanol-gasoline blend (E10) and pure ethanol (E100). Computational fluid dynamic (CFD) simulations of a stock and prechamber retrofitted single-cylinder Caterpillar C9.3B are carried out using CONVERGE. Prechamber equivalence ratio at spark timing, prechamber spark timing advance, and main chamber injection strategy are assessed with respect to their impact on ignition assistance performance and emissions characteristics relative to a diesel baseline at the same boundary conditions. Simulation results indicate that PC-MCC is flex-fuel capable and operates well for both E10 and E100 at the operating conditions considered. The results demonstrate that the use of a pilot-main injection strategy enables spark timing in the prechamber to be advanced and thus reduces spark plug firing pressure while maintaining robust ignition assistance. Results also indicate that the rich prechamber operation is favored for improved ignition assistance capabilities. Furthermore, the findings of this work suggest that a heavy-duty vehicle using a PC-MCC engine can utilize any blend of gasoline and ethanol, up to including pure ethanol, with no major sacrifices in performance relative to the diesel engine.

33 ADVANCED PROPULSION SYSTEMS↗

Flow and Heat Transfer Experimental Study for 3D-Printed Solar Receiving Tubes With Helical Fins at Internal Surface

3D-printing technology was applied to fabricate novel solar thermal collection tubes that have internal heat transfer enhancement fins and external surfaces with high solar absorptivity and low emissivity due to the ability to use different materials in one tube. Helical fins were selected to introduce circumferential flow and thus minimize the circumferential temperature difference of the tube that receives sunlight on one side. The structures of the helical fins were previously optimized from computational fluid dynamics (CFD) analysis with the objective of low entropy production rate by looking for high heat transfer coefficient and relatively lower pressure loss. High-temperature alloy, Inconel-718, was used to 3D print the tubes, which can resist corrosion for the potential application of molten chloride salts as heat transfer fluid. Experimental tests were carried out using water as the heat transfer fluid with the high heat flux provided by a tubular furnace heater. The tested Reynolds number ranges from 3.9 × 10 3 to 6.1 × 10 4 . Heat transfer coefficients of up to 2.8 times that of the smooth tube could be obtained with the expense of increased pressure loss compared to that of the smooth tube. The total system entropy generation can be significantly reduced due to the benefit of heat transfer enhancement that is greater than the expenses of the increased pressure loss. The experimental results of the 3D-printed heat transfer tubes confirmed the CFD-based results of fin optimization. Furthermore, the novel heat transfer tube is recommended for application in concentrating solar power systems.

14 SOLAR ENERGY↗

Thermomechanical Stress and Creep-Fatigue Analysis of a High-Temperature Prototype Receiver for Heating Particles

This work presents a three-dimensional (3D) thermomechanical model of a prototype-scale enclosed light trapping solar receiver for heating particles. Results of the thermoelastic model are used to estimate receiver lifetime under maximum flux conditions. A computational fluid dynamics (CFD) model is first developed to predict the temperature fields in a multi-panel assembly under steady operating conditions. Solar flux distributions on the receiver are obtained from the software package SolTrace and applied to the 3D thermal model. The subsequent particle heating is captured through a simplified 1D energy balance. Panel reradiation is considered through a surface-to-surface radiation model and natural convection loss to the surrounding air is captured in a representative fluid domain surrounding the receiver. The resulting temperature fields from the CFD analysis are used as inputs for a thermoelastic mechanical model with representative boundary conditions. With the resultant temperature and stress fields, a creep-fatigue damage and lifetime analysis is performed using the linear damage accumulation (LDA) theory. The Manson-Coffin formula and Larson Miller correlation are used to calculate the fatigue and creep, respectively. A maximum damage (corresponding to a 30-year service life) is defined for design assessment. The model was first developed and verified in detail by comparing with published results in the literature (temperature and stress profiles and distributions, and creep/fatigue damage fractions) for tubular solar receivers with supercritical carbon dioxide as the working fluid. It was then implemented to model a planar-cavity receiver with various design parameters. Specifically, three different design geometries are considered, and the results show that a maximum temperature of approximately 1200 K could be reached for each design with the given incident solar flux, with the main difference being the distribution of these temperatures. Preliminary resulting stresses for the small-scale prototype without design optimization vary from 20 MPa to 250 MPa for each design, with the maximum stresses occurring on the front face and concave geometry on the side of the panel. In future work, the developed methodology shown here will be applied to analyze a full-scale (50-150 MWth) receiver.

concentrated solar power↗

Robust Solution Verification Experiments on Nonuniform Meshes

The activities of verification, validation, and uncertainty quantification (VVUQ) provide a comprehensive means to assess the credibility of computational models. Within VVUQ, solution verification assesses numerical errors and evaluates whether the simulation is sufficiently accurate for its intended applications. As computational modeling gains traction in the development of complex, high-consequence systems, the need for robust solution verification intensifies, particularly because experimental data for these systems are often limited. This work examines improvements in the robustness of Richardson extrapolation (RE), a method commonly used in solution verification to study the discretization error of computational models using a power law. Nonuniform mesh refinement is discussed alongside other pollutants that affect the robustness of the power law model. Maximum likelihood estimation (MLE) is proposed as a robust strategy to address the uncertainty generated by nonuniform mesh refinement. An exploratory computational fluid dynamics (CFD) study of a 2D planar Poiseuille flow is conducted to determine if nonuniform mesh noise can be modeled with this MLE approach for more robust RE.

Weinmeister, Justin [ORNL] (ORCID:0000000160090237↗

Flashback of Hydrogen-Methane Mixtures in a Fixed-Vane, Low-Swirl Burner at Atmospheric Pressure

Equivalence ratio at flashback, 𝜙FB, of low swirl injectors for 60% – 95% (H2, by mole) hydrogen-methane mixtures was determined experimentally in a < 100 kW atmospheric pressure burner. Pre-mixer velocities equaled 5, 7.5, and 10 m/s, with the reactants entering the burner at 294 K. The swirler design was selected following a literature review, and the test articles were additively manufactured. The test articles had two turning angles, 26° and 33°, with center-body hole diameters equaling 1.08 and 1.16 mm. A test article was built with each parameter combination, which yielded four test articles with measured swirl numbers varying from 0.43 to 0.49. The results showed agreement with prior work. 𝜙FB increased with increasing pre-mixer velocity and decreasing percent hydrogen. 𝜙FB was also observed to increase with swirl angle and decrease with increasing center-body hole diameter. The data was correlated using a semi-empirical quenching approach, with 𝐶quench=0.01, where the chemical time scale was calculated using a perfectly stirred reactor (PSR) model. Two different chemical kinetics mechanisms independently yielded the same results. The correlation makes the present small-scale laboratory results relevant to large, high-pressure gas turbines. The flashback data set is available for computational fluid dynamics (CFD) validation.

Searle, Matthew↗

Design, Analysis, and Experimental Testing of Hydrogen Lean Direct Injection Nozzles at Elevated Pressure

Abstract There are many challenges of commissioning a hydrogen combustor into future gas turbine engines; especially regarding achieving emissions goals. Previously, Escudero et al. and Tran et al. conducted a study to adapt the liquid fuel Lean Direct Injection (LDI) concept from Jet-A to gaseous natural gas-hydrogen blends and pure hydrogen [1], [2]. Experimental data was collected at atmospheric conditions using a Box Behnken design of experiments. The design of experiments suggested that biasing the air split in favor of the inner air circuit and increasing the swirl strength of this inner air passage resulted in improved NOx emissions, while the inverse was true for stability, which was quantified by studying the lean blowoff point (LBO) [1], [2]. The trends revealed by the original experiment [1], [2] provided a design direction for further iterations of the experimental hardware. The study presented herein describes the further investigation of such LDI injectors through experimental methods and computational fluid dynamic (CFD) simulations at atmospheric conditions, which were used to identify potential flow behaviors driving enhanced emissions performance. Further evaluation of select injectors from both studies was then conducted at elevated pressures up to 6 atmospheres. The results from both experiments are presented in this study, which include flame observations, emissions measurements, and operational challenges. NOx emissions results are reported on a volume basis in ppmvd corrected to 15% O2 and corrected for fuel. A predictive model for relating NOx emissions to test conditions at atmospheric conditions show high significance to adiabatic flame temperature while little to no significance to fuel composition for the best performing configurations. The results illustrate the connection between atmospheric testing and testing elevated pressures. The design direction indicated by the initial tests and CFD results in promising configurations for implementation into a Multi-point LDI array.

08 HYDROGEN↗

A NUMERICAL STUDY OF REFRIGERANT LEAKAGE FROM A PROPANE-BASED REFRIGERATION SYSTEM

Low Global Warming Potential (GWP) refrigerants, e.g. Propane (R290) (GWP = 3) and R1234yf (GWP = 4), are becoming a popular choice among refrigeration and HVAC systems, such as vending machines and air conditioning. However, most of Low GWP refrigerants are flammable (ASHRAE class A2 or A3), so safety is a very important consideration when designing and deploying low GWP refrigerants-based equipment in buildings. In the event of a refrigerant leak, the flammability of low GWP refrigerants depends on the local concentration of the refrigerant within the vicinity of the leak. In addition, the low GWP refrigerants concentration is affected by the indoor air environment, such as air flow rate and temperature. In the present study, a computational fluid dynamics (CFD) model is developed to model the flow and temperature of air surrounding a propane-based refrigeration system, as well as the concentration of leaked low GWP refrigerants surrounding the refrigeration system. The model results reveal the concentration distribution of low GWP refrigerants within the building, and as a result, the flammable regions within the building can be identified. Moreover, different ventilation layouts will be tested using the model to improve the design of ventilation. The numerical model can assist in the design of ventilation systems to minimize flammable regions within buildings due to leakage of flammable refrigerants.

Zhang, Mingkan↗

NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE

In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.

detailed chemical kinetics↗

Numerical Investigation of a Central Fuel Property Hypothesis under Boosted Spark-ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel's performance irrespective of its chemical composition, was numerically investigated In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel's knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (= 98) and MON (= 90), were first generated using a non-linear regression model based on artificial neural network (ANN). Three unconventional bioblendstocks were included in the analysis: Di-isobutylene (DIB), Isobutanol and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed was tabulated for each blend as a function of pressure, temperature and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON and engine operating conditions, the TPRF-Anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (DI) expectation was shown to be caused by differences in laminar flame speed (LFS). However; it was found that relatively large fuel-specific differences in LFS (> 20%) would have to be present to cause any appreciable deviation from the 01 framework. Otherwise, RON and MON would still be robust enough to predict a fuel's knock-limited performance.

Pal, Pinaki↗

Performance and Emissions of an SI Engine Fueled With DME-Propane Blends

Dimethyl Ether (DME) is an alternative fuel that can be produced renewably and has the potential for lower CO and NOx emissions than conventional petroleum-based fuels. Blending DME with another gaseous fuel such as propane, which has a lower knock tendency than gasoline, can allow this fuel to be leveraged on SI engines. In this study, the use of DME-propane blends on a spark ignition (SI) engine was studied via computer simulations in order to understand the impact on engine performance and emissions and to identify the knock limitations of using such fuel blends. A 2L Hyundai SI engine was modeled in GT Power and the model was validated by comparing it with computational fluid dynamics (CFD) simulation results. Starting from pure propane, DME was added incrementally until knock was observed in the engine. Results indicate that it is feasible to run propane with DME percentages up to 35% before severe knock impacts were observed. The BTE was higher, but the BSFC also increased for DME-propane blends as compared to gasoline. Here, an increase in NOx emissions was detected along with a significant decrease in CO emissions. CO 2 emissions declined for propane as compared to gasoline but increased with the addition of DME.

alternative fuels↗

Multi-Scale Thermal Analysis for Design of SiC-Based Medium Voltage Motor Drive: Preprint

This paper describes a multi-scale thermal analysis approach for the design of an air-cooled 1.7-kV SiC MOSFET-based medium-voltage variable-speed motor drive. The National Renewable Energy Laboratory supported the project objectives by providing thermal management design support to the project partner The Ohio State University (OSU). The scope of the models and limited resources required efficient and flexible thermal models to be developed. Two modeling techniques are described that significantly reduced model run time and enabled more complex models to be run faster while retaining needed accuracy. The first technique uses the effectiveness-NTU method to extract convection boundary conditions from a computational fluid dynamics (CFD) model that can be applied to a fast-running finite element analysis model. The second is a porous media technique that enables system-level CFD simulations that incorporate effects from heat exchangers (e.g., pin fin heat sinks) that run in a fraction of the time required for fully resolved CFD simulations. The multi-scale approach to the thermal analysis enabled fast and accurate simulation for the converter design ranging from the die level up to the full system with 36 submodules. The modeling results were validated against experimental data from system tests performed by OSU.

DIRECT ENERGY CONVERSION↗

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial-flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that are neglected for wind turbines. These include buoyancy, added mass and inertial loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbines at a speed comparable to real time. While efficient for large sets of load cases and design studies, mid-fidelity codes make simplifying assumptions that may impact their accuracy. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system was modeled both with OpenFAST and with the commercial CFD code STAR-CCM+. The CFD model used a three-dimensional unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modeled with a custom coupling of the CFD solver with the open-source code MoorDyn. This improves upon the existing quasi-static catenary solver in STAR-CCM+, which lacks seabed contact or line-to-line connections. Simulation results for a floating marine turbine are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗

OVERVIEW OF PARTICLE DEPOSITION MODELS FOR SPENT NUCLEAR FUEL STORAGE SYSTEMS

Deposition models were built to evaluate contaminant deposition on spent nuclear fuel (SNF) canisters. The primary contaminant of concern is chloride, which is dispersed in the atmosphere and then deposits onto the canisters. During dry storage, the primary degradation process is likely to be Chloride Induced Stress Corrosion Cracking (CISCC) at the heat-affected zones of the canister welds. It is known that stainless steel canisters are susceptible to CISCC; however, the rate of chloride deposition onto the canisters is poorly known, based on sparse field data from a small number of sites. The models presented in this study could be useful for determining the rate of deposition on the canisters and the likelihood of CISCC to help with SNF canister ageing management. The deposition models were developed with the commercial computational fluid dynamics (CFD) code STAR-CCM+. Various deposition mechanisms were considered and incorporated into the models, and a sensitivity study was conducted to determine the most important mechanisms for deposition within a SNF storage system. The models included both a vertical and horizontal configuration storage system: NAC International’s Modular, Advanced Generation, Nuclear All-purpose STORage System (MAGNASTOR®) and a NUHOMS® horizontal storage module respectively. The resulting canister deposition on the horizontal canister is visually compared with inspection data taken onsite at the Calvert Cliffs Nuclear Power Plant. These models are preliminary, and development of the models will continue. Future validation exercises are currently being planned, including the Canister Deposition Field Demonstration (CDFD) effort funded by U.S. Department of Energy office of Nuclear Energy office of Spent Fuel Waste Science and Technology. The goal of the modeling presented is to demonstrate a potential modeling technique that could be used to plan and inform SNF canister ageing management programs with predictive models for the timing and occurrence of canister CISCC.

Suffield, Sarah R.↗

Computational Fluid Dynamic Modeling of Dry Cask Simulator with Crosswind

The purpose of this study is to create a STAR-CCM+ model of a Belowground Vertical Dry Cask Simulator (BVDCS) at Sandia National Laboratories (SNL) and validate the model with SNL’s experimental results. The BVDCS consists of a single boiling water reactor assembly fitted with electric heaters encompassed by a containment vessel and shell to represent a belowground spent nuclear fuel (SNF) dry storage system. Blowers are located near the inlet and outlet of the BVDCS to simulate crosswind conditions. In addition to the experimental results, the STAR-CCM+ model developed for this study is compared with a previous computational fluid dynamics (CFD) model in a different software program, which is used as a software-to-software benchmark. The experimental results provide a dataset to compare the STAR-CCM+ model results for a variety of different conditions. The main objective is to validate and improve STAR-CCM+ CFD models for spent nuclear fuel storage systems with explicitly modeled external environments and “wind driven” crossflows. These CFD models aide in the study of external particle deposition in spent nuclear fuel storage systems, which is important to predicting the significance of chloride induced stress corrosion cracking (CISCC). In addition to experimental comparison, a sensitivity analysis study is performed using the STAR-CCM+ model. The sensitivity analysis provides a quantitative assessment of the sensitivity of various parameters. This helps provide information on various parameters that are of particular importance to constructing a model representative of real life systems. The STAR-CCM+ model compared well to the experimental results showing similar responses to changes in cross wind flow, and a number of parameters are identified for model improvement.

Jensen, Ben J.↗