Search NASA⌕ Search

SEARCH · Search NASA

Results for “Three dimensional computational fluid dynamics model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effect of Surface Roughness on Dynamic Stall in Pitching Motion

Dynamic stall plays a critical role in determining the performance and stability of a wide range of fluid-dynamics systems in various engineering applications. This unsteady aerodynamic phenomenon is particularly significant for maneuvering aircraft wings, jet aircraft subjected to gust encounters, helicopter rotor blades, and wind turbine blades [1–3]. The prediction of the dynamic stall vortex (DSV) is challenging due to factors such as unsteady aerodynamics, three-dimensional (3-D) effects, turbulence and flow separation, incoming gust, and surface impact effects [4–6]. Hence, advanced computational techniques and modeling approaches in computational fluid dynamics (CFD) would be required to enhance the accuracy and reliability of DSV predictions in dynamic stall scenarios. In the past, Batther and Lee [7] employed delayed detached eddy simulations (DDES) to understand the flow physics associated with the onset of dynamic stall. Their approach demonstrated that DDES achieves results comparable to those obtained from large-eddy simulations at a reduced computational cost. In another study, Khalifa et al. [8] examined the 3-D aspects of dynamic stall on a NACA 0012 airfoil using DES solvers. In conclusion, the findings underlined the superiority of 3-D simulations over two-dimensional approaches, particularly in predicting the lift coefficient values and capturing dynamic stall stages more precisely.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

ADVANCING TURBULENT SPRAY AND COMBUSTION MODELS FOR COMPRESSION IGNITION ENGINE SIMULATIONS

This thesis seeks to investigate the turbulent mixing influence on spray atomization and combustion processes encountered in compression ignition diesel engines. Despite greater thermal efficiency of diesel engine than spark ignition engine, the nature of stratified air-fuel mixture and non-premixed flame gives rise to unacceptable levels of nitrogen oxides (NOx) and particulate matter (PM), thus the use of diesel engines has often been limited to heavy-duty vehicle and industrial power sources. However, recent advancement in diesel engine combustion strategies, e.g. low temperature combustion (LTC), has demonstrated promising pathways towards improvement in the engine-out pollutants. Therefore, particularly in the effort of computer-aided engine design tasks, such a new engine design concept requires more accurate modeling techniques applicable over a broader range of engine operating conditions than those of conventional engine strategies. In the notion of challenges in new engine operating conditions, this thesis aims to present successful implementation of improvement in numerical modeling techniques in high-pressure spray atomization and resulting turbulent spray flame of interest. Three-dimensional Computational Fluid Dynamics (CFD) in in-cylinder turbulent combustion is considered an integral part of engine design progress, but rather a cost-prohibitive to apply over a broad range of engine relevant conditions. In spite of successful use of existing spray atomization modeling, prior researchers have pointed out some degree of failure in LTC targeted injection strategies. Furthermore, finite rate and strong nonlinearity of chemistry influenced by local turbulent mixing still remain in challenges to account for in cost-efficient CFD analysis. In this context, a new attempt of hybrid spray primary breakup modeling is presented and demonstrated in successful application aimed at LTC technique. In addition, the Representative Interactive Flamelets (RIF) model with aid of multi-flamelets approach is extensively assessed in terms of predictive capability against classical combustion model. The combustion model employed in this study are fully examined in the general diesel combustion metric, e.g., ignition delay and flame lift-off length as well as newly suggested test metric, combustion recession. The combustion recession has been recently idenfied, but still remain largely unknown. Since the governing physics of this phenomenon is characterized by turbulent mixing coupled with finite rate chemistry, this can be considered as a relevant test metric for turbulent combustion models. In addition, very recent experimental studies have introduced a new non-sooting diesel combustion technique by manipulating direct injection method. The ducted fuel injection (DFI) has thus been demonstrated with its potential of low soot emissions. Knowing that the duct equipped ahead of injector nozzle was identifed to enhance turbulent mixing, investigations of DFI combustion may prove the effectiveness of turbulece-chemistry interaction modeling. This thesis presents comprehensive understandings of aformentioned diesel combustion techniques in terms of several important physics keywords, e.g., turbulent mixing and detailed chemistry.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2019

Advancing Turbulent Spray and Combustion Models for Compression Ignition Engine Simulations

This thesis seeks to investigate the turbulent mixing influence on spray atomization and combustion processes encountered in compression ignition diesel engines. Despite greater thermal efficiency of diesel engine than spark ignition engine, the nature of stratified air-fuel mixture and non-premixed flame gives rise to unacceptable levels of nitrogen oxides (NOx) and particulate matter (PM), thus the use of diesel engines has often been limited to heavy-duty vehicle and industrial power sources. However, recent advancement in diesel engine combustion strategies, e.g. low temperature com-bustion (LTC), has demonstrated promising pathways towards improvement in the engine-out pollutants. Therefore, particularly in the effort of computer-aided en-gine design tasks, such a new engine design concept requires more accurate modeling techniques applicable over a broader range of engine operating conditions than those of conventional engine strategies. In the notion of challenges in new engine operat-ing conditions, this thesis aims to present successful implementation of improvement in numerical modeling techniques in high-pressure spray atomization and resulting turbulent spray flame of interest. Three-dimensional Computational Fluid Dynamics (CFD) in in-cylinder turbu-lent combustion is considered an integral part of engine design progress, but rather a cost-prohibitive to apply over a broad range of engine relevant conditions. In spite of successful use of existing spray atomization modeling, prior researchers have pointed out some degree of failure in LTC targeted injection strategies. Furthermore, finite rate and strong nonlinearity of chemistry influenced by local turbulent mixing still re-main in challenges to account for in cost-efficient CFD analysis. In this context, a new attempt of hybrid spray primary breakup modeling is presented and demonstrated in successful application aimed at LTC technique. In addition, the Representative Interactive Flamelets (RIF) model with aid of multi-flamelets approach is extensively assessed in terms of predictive capability against classical combustion model. The combustion model employed in this study are fully examined in the general diesel combustion metric, e.g., ignition delay and flame lift-off length as well as newly sug-gested test metric, combustion recession. The combustion recession has been recently idenfied, but still remain largely unknown. Since the governing physics of this phe-nomenon is characterized by turbulent mixing coupled with finite rate chemistry, this can be considered as a relevant test metric for turbulent combustion models. In ad-dition, very recent experimental studies have introduced a new non-sooting diesel combustion technique by manipulating direct injection method. The ducted fuel in-jection (DFI) has thus been demonstrated with its potential of low soot emissions. Knowing that the duct equipped ahead of injector nozzle was identifed to enhance turbulent mixing, investigations of DFI combustion may prove the effectiveness of turbulece-chemistry interaction modeling. This thesis presents comprehensive under-standings of aformentioned diesel combustion techniques in terms of several important physics keywords, e.g., turbulent mixing and detailed chemistry

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

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↗
NASA NTRS1996

The AFDD International Dynamic Stall Workshop on Correlation of Dynamic Stall Models with 3-D Dynamic Stall Data

A variety of empirical and computational fluid dynamics two-dimensional (2-D) dynamic stall models were compared to recently obtained three-dimensional (3-D) dynamic stall data in a workshop on modeling of 3-D dynamic stall of an unswept, rectangular wing, of aspect ratio 10. Dynamic stall test data both below and above the static stall angle-of-attack were supplied to the participants, along with a 'blind' case where only the test conditions were supplied in advance, with results being compared to experimental data at the workshop itself. Detailed graphical comparisons are presented in the report, which also includes discussion of the methods and the results. The primary conclusion of the workshop was that the 3-D effects of dynamic stall on the oscillating wing studied in the workshop can be reasonably reproduced by existing semi-empirical models once 2-D dynamic stall data have been obtained. The participants also emphasized the need for improved quantification of 2-D dynamic stall.

Tan, C. M.↗
NASA NTRS2003

Flutter Calculations for an Experimental Fan

During testing, an experimental forward-swept fan encountered flutter at part-speed conditions. A three-dimensional propulsion aeroelasticity code, based on a computational fluid dynamics (CFD) approach, was used to model the aeroelastic behavior of this fan. This paper describes the flutter calculations and compares the results to the experimental measurements. Results of sensitivity studies are also presented that show the relative importance of different aspects of aeroelastic modeling.

Bakhle, Milind A.↗
NASA NTRS2002

Progress in the Development of a Nozzle Design Methodology for Pulsed Detonation Engines

The Johns Hopkins University Applied Physics Laboratory (JHU/APL), in support of the NASA Glenn Research Center (NASA GRC), is investigating performance methodologies and system integration issues related to Pulsed Detonation Engine (PDE) nozzles. The primary goal of this ongoing effort is to develop design and performance assessment methodologies applicable to PDE exit nozzle(s). APL is currently focusing its efforts on a common plenum chamber design that collects the exhaust products from multiple PDE tubes prior to expansion in a single converging-diverging exit nozzle. To accomplish this goal, a time-dependent, quasi-one-dimensional analysis for determining the flow properties in and through a single plenum and exhaust nozzle is underway. In support of these design activities, parallel modeling efforts using commercial Computational Fluid Dynamics (CFD) software are on-going. These efforts include both two and three-dimensional as well as steady and time-dependent computations to assess the flow in and through these devices. This paper discusses the progress in developing this nozzle design methodology.

Leary, B. A.↗
DOE OSTI2021

Three-Dimensional CFD Investigation of Pre-Spark Heat Release in a Boosted SI Engine

Low-temperature heat release (LTHR) in spark-ignited internal combustion engines is a critical step toward the occurrence of auto-ignition, which can lead to an undesirable phenomenon known as engine knock. As such, correct predictions of LTHR are of utmost importance to improve the understanding of knock and enable techniques aimed at controlling it. While LTHR is typically obscured by the deflagration following the spark ignition, extremely late ignition timings can lead to LTHR occurrence prior to the spark, i.e., pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection SI engine was numerically investigated using three-dimensional computational fluid dynamics (CFD). A hybrid approach was used, based on the G-equation model for representing the turbulent flame front and the multi-zone well-stirred reactor model for tracking the chemical reactions within the unburnt region. A recently developed best practice was also employed which keeps the well-stirred reactor model active throughout the entire simulation. This allowed for correct predictions of the previous cycle trapped residuals which have a considerable effect on the onset of LTHR. Multi-cycle simulations were conducted using Co-Optima alkylate and E30 fuels. The predicted in-cylinder pressure and heat release rate agreed well with the experimental data and served as validation for the CFD model. Following the initial validation, the dynamics of PSHR was discussed and a series of parameters of interest were assessed. First, the effect of exhaust valve temperature was investigated, qualitatively highlighting the importance of boundary conditions uncertainty. Further analyses were carried out on the effects of fuel properties, including laminar flame speed (LFS) and heat of vaporization (HOV). The results indicate that PSHR phasing is slightly advanced with lower LFS as more trapped unburnt fuel is made available for the next cycle. A similar trend was observed with lower HOV as less intense spray cooling led to higher mixture temperatures, i.e., higher mixture reactivity. Finally, a comparison of Co-Optima alkylate and E30 fuels was made using the pressure-temperature trajectory framework. It was shown that the differences between Co-Optima alkylate and E30 in PSHR tendency are correlated with both HOV and chemical effects.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS2024

Turbulence Model Validation through Joint Experimental / Computational Studies of Separated Flow Over a Three-Dimensional Tapered Bump: Part II - Computational Investigation

This report serves as documentation of activities related to Computational Fluid Dynamics (CFD) analysis applied to the Smooth Body Separation Experiment designed to produce validation-quality experimental data from wind-tunnel tests performed at the University of Notre Dame. All CFD analysis utilized the Boeing CFD (BCFD) flow solver. BCFD was used to analyze flow over the Boeing Speed Bump with second-order accuracy in space and time. Inviscid fluxes were computed using the HLLE scheme for steady-state simulations, and the Low-Dissipation Roe scheme for time-accurate simulations. The working fluid was modeled as calorically perfect air using Sutherland's law for viscosity. Thermal diffusivity was determined by assuming a constant molecular Prandtl number of 0.72, and a turbulent Prandtl number of 0.9. Reynolds Averaged Navier-Stokes (RANS) simulations were conducted using the one-equation Spalart-Allmaras turbulence model with rotation correction3 and the quadratic constitutive relation (SA-RC-QCR), as well as Menter's two-equation, Shear-Stress Transport (SST) model. Time-accurate computations were conducted using Delayed Detached eddy Simulation (DDES) with the one-equation RANS model, and Large Eddy Stress balance (LESb) with the two-equation RANS model. Implicit dual-time stepping was used to achieve second-order temporal accuracy. The experimental investigation conducted during the course of this study is described in detail in an accompanying report.

Flow separation↗
NASA NTRSDate not supplied

Modeling of Particle Transport and Impact in Transient, Multiphase Oxidizer-Rich Environments Using Computational Fluid Dynamics and Lagrangian Particle Tracking

Particle impacts from foreign object debris (FOD) in an oxidizer-rich environment have been shown to be a credible and significant risk for ignition and failure of components and vehicles. The risk of particle ignition is dependent on many factors, including particle impact velocity/energy and fluid composition (i.e. fluid phase, species presence, and state). Prediction of particle transport and particle impacts in relevant environments, including transient, multiphase flows, is important for assessing the risk and avoiding repeated expensive system level tests. The capability to model both the fluid dynamics and particle transport in complex three-dimensional geometries and multiphase environments has been developed at Marshall Space Flight Center (MSFC) in the propulsion fluid dynamic branch (ER42) by utilizing computational fluid dynamics (CFD) simulations and Lagrangian particle tracking. Using CFD tools, including density-based Loci/CHEM with a cavitation model and pressure-based Loci/STREAM with a volume of fluids model, MSFC/ER42 has provided particle impact predictions to multiple projects for oxidizer compatibility assessments. Simulations of transient, multiphase environments such as liquid slugs flowing through tubing and cavitating flow in valves have been completed alongside one-way coupled Lagrangian particle tracking to predict particle impacts.

M.V. Fischels↗
NASA NTRS2023

Turbulence Model Validation Through Joint Experimental /Computational Studies of Separated Flow Over A Three-Dimensional Tapered Bump: Part I - Experimental Investigation

The results of a three-year experimental campaign aimed at comprehensively documenting the separated flow over a three-dimensional bump are presented with the purpose of generating a benchmark experimental database useful in validating computational fluid dynamics (CFD) flow simulations. The bump model geometry was designed to provide well-defined and repeatable smooth-body flow separation conditions that were suitable for both experiments and simulations. The bump had a Gaussian streamwise profile with a constant height equal to 8.5% of its width over the central 60% of the test section width. The remaining 40% were outboard spanwise portions that gradually taper to zero using an error function profile to minimize side-wall boundary layer effects. The model was immersed in a turbulent boundary layer that was developed on a suspended flat plate in the Notre Dame Mach 0.6 Wind Tunnel. In order to document the effect of the incoming boundary layer thickness on the flow separation, the bump model could be located at two streamwise positions. The mean velocity and turbulence intensity of the wind tunnel freestream flow field and approaching turbulent boundary layer were fully documented. The measurements of the flow separation region included surface visualization, wall shear stress using oil-film interferometry, mean and dynamic surface pressure, and planar and stereoscopic particle image velocimetry. The experiments were conducted over a range of Mach numbers from 0.05 to 0.2 corresponding to a range of Reynolds numbers based on the test section spanwise dimension (0.914m) of 1.0 × 10 6 ≤ Re L = U ∞ L/ν ≤ 4.0 × 10 6 . The bulk of the results are presented for the higher Mach number conditions of 0.1 and 0.2 with Re L = 2.0 × 10 6 and 4.0 × 10 6 , respectively. Extensive uncertainty analysis of the data was performed. The data is archived in the NASA Langley Turbulence Modeling Resource website at https://turbmodels.larc.nasa.gov/Other_exp_Data/speedbump_sep_exp.html. In addition to the experiments, a computational effort was made in parallel by the CFD group at Boeing Research & Technology highlighting the usefulness of the data set, which is outlined in the accompanying CFD report.

Flow separation↗
DOE OSTI2017

3D-FE Modeling of 316 SS under Strain-Controlled Fatigue Loading and CFD Simulation of PWR Surge Line

In financial year 2017, we are focusing on developing a mechanistic fatigue model of surge line pipes for pressurized water reactors (PWRs). To that end, we plan to perform the following tasks: (1) conduct stress- and strain-controlled fatigue testing of surge-line base metal such as 316 stainless steel (SS) under constant, variable, and random fatigue loading, (2) develop cyclic plasticity material models of 316 SS, (3) develop one-dimensional (1D) analytical or closed-form model to validate the material models and to understand the mechanics associated with 316 SS cyclic hardening and/or softening, (4) develop three-dimensional (3D) finite element (FE) models with implementation of evolutionary cyclic plasticity, and (5) develop computational fluid dynamics (CFD) model for thermal stratification, thermal-mechanical stress, and fatigue of example reactor components, such as a PWR surge line under plant heat-up, cool-down, and normal operation with/without grid-load-following. This semi-annual progress report presents the work completed on the above tasks for a 316 SS laboratory-scale specimen subjected to strain-controlled cyclic loading with constant, variable, and random amplitude. This is the first time that the accurate 3D-FE modeling of the specimen for its entire fatigue life, including the hardening and softening behavior, has been achieved. We anticipate that this work will pave the way for the development of a fully mechanistic-computer model that can be used for fatigue evaluation of safety-critical metallic components, which are traditionally evaluated by heavy reliance on time-consuming and costly test-based approaches. This basic research will not only help the nuclear reactor industry for fatigue evaluation of reactor components in a cost effective and less time-consuming way, but will also help other safety-related industries, such as aerospace, which is heavily dependent on test-based approaches, where a single full-scale fatigue test can cost millions of dollars and require years of effort to conduct. Toward our goal of demonstration of fully mechanistic fatigue evaluation of reactor components, we also started work on developing a component-level computer model of reactor components, such as 316 SS surge line pipe. This requires developing a thermal-mechanical stress analysis model of the reactor surge line, which, in turn, requires time-dependent temperature and stratification information along the boundary of the pipe. Toward that goal, CFD models of surge lines are being developed. In this report, we also present some preliminary results showing the temperature conditions along the surge line wall under reactor heat-up, cool-down, and steady-state power operation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Numerical Modeling and Parametric Study of a Dual Purpose Underground Thermal Battery

Buildings’ thermal energy systems for providing space heating/cooling consume a considerable amount of electricity nationwide. Integrating thermal energy storage with buildings’ thermal systems has the potential to shave the peak electric demand and to overcome the mismatch between the intermittent renewable power and the fluctuating daily demand for electricity. A novel dual-purpose underground thermal battery (DPUTB) has been developed to provide diurnal thermal energy storage in conjunction with a dual-source heat pump (DSHP)system. A two-dimensional numerical model is developed to simulate the transient performance of the DPUTB. This model has similar accuracy as a detailed three-dimensional model developed with commercial heat transfer and fluid dynamics simulation program but with 1000 times faster computational speed. A parametric study is performed to identify a design of the DPUTB that can provide 14 kWh thermal storage capacity for each charging and discharging cycle. A system simulation has been developed to assess the performance of the designed DPUTB when it is integrated with a DSHP for conditioning a typical residential building in summer and winter. Simulation results indicate that the DPUTB can meet all the design requirements in both summer and winter with the simple rule-based control strategy and the overall integrated DPUTB and DSHP system can shift the electric load for meeting the thermal demand from the peak hours to off-peak hours during a day, which would result in considerable energy cost savings if a Time of Use electricity tariff is applicable.

25 ENERGY STORAGE↗
NASA NTRS1996

Endwall Heat Transfer Measurements in a Transonic Turbine Cascade

Turbine blade endwall heat transfer measurements are given for a range of Reynolds and Mach numbers. Data were obtained for Reynolds numbers based on inlet conditions of 0.5 and 1.0 x 106, for isentropic exit Mach numbers of 1.0 and 1.3, and for freestream turbulence intensities of 0.25% and 7.0%. Tests were conducted in a linear cascade at the NASA Lewis Transonic Turbine Blade Cascade Facility. The test article was a turbine rotor with 136' of turning and an axial chord of 12.7 cm. The large scale allowed for very detailed measurements of both flow field and surface phenomena. The intent of the work is to provide benchmark quality data for computational fluid dynamics (CFD) code and model verification. The flow field in the cascade is highly three-dimensional as a result of thick boundary layers at the test section inlet. Endwall heat transfer data were obtained using a steady-state liquid crystal technique.

Giel, P. W.↗
DOE OSTI2023

An Efficient Three-Dimensional CFD-Based Numerical Wave Tank for a Wave Energy Converter in Extreme Irregular Waves

A numerical wave tank approach for computational fluid dynamics (CFD) modelling of an extreme irregular seastate is presented. The technique couples a potential flow solution with a CFD solver for more efficient numerical predictions. This method has recently become attractive both for the research community and the industry working with offshore structures. The model is used to determine the response of a submerged pressure differential wave energy converter (WEC) in a fully nonlinear irregular wave condition using the high-fidelity CFD code, STAR-CCM+. Potential flow based numerical models are commonly used to predict motions and performance of wave energy converters. Wave kinematics can deviate from potential flow predictions for extreme wave conditions; the excitation loads on an absorber can also be increasingly influenced by viscous effects, not predicted by potential flow engineering level models. In these extreme conditions, a Reynolds-averaged Navier-Stokes CFD model can better predict motions and loads for a WEC. Long time series with varying random seed numbers can be used to identify singular extreme wave events from a stochastic irregular sea state. This approach simulates a more realistic wave series for a given sea state than a regular wave or a focused wave. However, it is computationally infeasible to run these long time series for three-dimensional (3D) CFD simulations. In this work, two-dimensional (2D) CFD simulations with a long domain allow the full development of an extreme nonlinear wave condition. The results are used to identify extreme events from a 50-year storm condition for the PacWave site off the coast of Oregon. A relatively short time window including this extreme event is then mapped to a 3D simulation using a user defined wave methodology. Convergence studies for domain length, wave forcing lengths, and time before the extreme event were conducted.

CFD↗
DOE OSTI2022

Coupled chemical–CFD modeling of unsealed dry storage of advanced test reactor spent fuel

We report given the current lack of long-term disposal for spent nuclear fuel storage in the United States, the current inventory of advanced test reactor (ATR) fuel at Idaho National Laboratory (INL) is expected to be stored in a dry storage facility for an extended storage period. As part of an effort to understand the conditions during a 50-year period, a three-dimensional computational fluid dynamics (CFD) of the unsealed canisters packed with spent fuel in the storage facility was constructed, this CFD model is coupled with radiolytic chemical reactions associated with water vapor and other species. The coupled model is simulated over a range of sensitivity parameters for the fuel decay heat, relative humidity, and oxyhydroxide film thickness for undried and fully-dried fuel. The most important parameter was the initial fuel decay heat. Over the range of cases, the hydrogen and nitric acid concentrations are low enough no significant impact is expected.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Modeling Approach for the Aluminum-clad Dry Storage Pilot using HFIR Fuel

To confirm that the dry storage of aluminum-clad research reactor spent nuclear fuel (ASNF) will remain within the safety envelope after applied drying schemes and that the resulting evolution of the gas space composition, temperature, and pressure conditions are understood, a dry storage pilot project is being established. The pilot will incorporate an instrumented lid for discrete interval or for on-demand gas composition and temperature monitoring of two DOE Standard Canisters (DSCs) loaded with three High Flux Isotope Reactor (HFIR) inner cores per DSC. Each DSC would be subjected to a separate alternative candidate drying scheme. Canisters will undergo 1 to 5 years of monitoring, including internal temperature and gas sampling to track pressure and composition changes. This report outlines the approach for modeling the ASNF-in-canister behavior in terms of evolving gas space conditions for the ASNF dry storage pilot using HFIR fuel. The ASNF has an adherent surface oxyhydroxide layer comprised of boehmite/bayerite that generates hydrogen when subjected to irradiation. Three-dimensional multi-physics computational fluid dynamics simulations will be executed to compute the thermal field within the DSC and provide inputs to a chemical model employed to compute pressure buildup as hydrogen is generated in the system. Implemented in Cantera, the chemical model solves gas phase and aluminum oxyhydroxide surface-mediated radiolysis reactions. Gas phase reactions are sourced from Wittman and Hanson (2015), whereas surface-mediated reactions are incorporated by fitting experimental data using an optimization algorithm (Abboud, 2023). Water radiolysis reactions from Wren and Ball (2001) are adopted with modifications as described in Abboud (2023c). Understanding the effect of the hydrogen buildup over time is important for long-term storage safety considerations. Modeling results will include the canister pressure, temperature, and composition evolution from the initial helium backfill with the addition of radiolytically-evolved chemical species (e.g., hydrogen and oxygen). The specific HFIR cores for the pilot program have not yet been selected, and the overall design is still in development. The CFD-chemical model used for this work will be based on prior models with necessary updates to allow for improved accuracy and efficiency. The experimental data obtained from the HFIR demonstration will be used to improve and validate the computational models to predict the ASNF-in-canister behavior.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Modeling of streamflow in a 30 km long reach spanning 5 years using OpenFOAM 5.x

Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.

58 GEOSCIENCES↗