Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Assessment of the 3420 Building Filtered Exhaust Stack Sampling Probe Location: Stack Verification Following Fan and Air Blender Additions

The Pacific Northwest National Laboratory (PNNL) 3420 Building, which is part of the Physical Sciences Facility (PSF), houses radiological capabilities that results in the requirement that emissions monitoring must be conducted for potential radionuclides in the exhaust air discharge of this building. The air monitoring system is required to conform to Title 40 of the Code of Federal Regulations part 61 (40 CFR 61) Subpart H, which in turns requires a sampling probe in the exhaust stream to conform to the criteria of American National Standards Institute / Health Physics Society (ANSI/HPS) N13.1-2011, Sampling and Monitoring Releases of Airborne Radioactive Substances from the Stack and Ducts of Nuclear Facilities. To support the air emissions permit for the 3420 Building on the PNNL campus, stack testing that used computational fluid dynamics (CFD) modeling as a surrogate stack and verification tests of velocity uniformity and flow angle on the retrofitted facility stack was performed. The ANSI/HPS N13.1-2011 criteria for the air monitoring probe location are that velocity uniformity, gaseous tracer uniformity, and particulate tracer uniformity must be less than or equal to 20%COV. Furthermore, no point in the sampling location may have a gaseous tracer concentration that varies from the mean concentration by more than 30%. Additionally, the flow angle at the sampling location must not be more than 20°. The CFD modeling of the stack, as reported by Recknagle et al. (2018) demonstrated that the stack meets the criteria at the probe location. The velocity uniformity and flow angle results from the 3420 stack verification tests, performed in October 2020, demonstrated that the CFD model results may be used to support the qualification of the stack sampling location. The measured velocity uniformity verification test result was 1.4%COV. This value is well within the uniformity criterion, which is that the velocity uniformity be =20%COV. Additionally, this value is well within the criterion that the actual stack measurement must be within 5% of the surrogate stack result of 2.1%COV when all four fans were operating. Additionally, the measured average flow angle at the 3420 stack monitor location was 15.5 degrees. Although this is higher than expected based on the CFD model, the result is =20 degrees, so the criterion is met. Based on these stack verification test results, the reconfigured 3420 Building filtered exhaust stack meets the qualification criteria given in the ANSI/HPS N13.1-2011 standard. Further changes to the system configuration or operating conditions that are outside the bounds described in this and the CFD report (Recknagle et al., 2018) may require additional tests and additional analysis to determine compliance with the standard.

3420 Building↗
DOE OSTI2020

Accelerated Aging Humidity Chamber for Nuclear Grade HEPA Filter Media

A unique accelerated aging humidity chamber for simulating the natural aging process of nuclear grade high efficiency particulate air (HEPA) filter media is being designed, implemented, and tested. The nuclear industry currently lacks information regarding the shelf life and service life of HEPA filter media. The Department of Energy recommends disposing of HEPA filters 10 years after the manufacture date for filters operated under dry conditions. Studies have shown the tensile strength and water repellency of HEPA filter media diminish with age. The chamber's principal function is to sustain an elevated relative humidity condition for an extended duration specified by the Arrhenius aging model. To accomplish this, various design parameters were established: (1) consistently achieve uniform humidity conditions within the chamber to equally expose media sheets, (2) ensure the chamber is sufficiently sealed to mitigate the uncontrolled infiltration of air, (3) confidently secure the media sheets to avert the introduction of unnecessary stressing/creasing, (4) fabricate the chamber from clear polycarbonate material to allow inspection of the media sheets, (5) provide a safe and ergonomic design for personnel, (6) ensure the chamber design is replicable and relatively easy to fabricate/assemble. A humidity source, fan, PID controller, humidity and temperature transmitters, and data loggers are required for proper operation and control of the chamber. The humidity source must accommodate target values up to 95% relative humidity throughout the chamber. The three humidity and temperature transmitters are strategically spaced throughout the chamber to ensure uniform conditions. The 10 cubic feet per minute (CFM) fan selected intends to accomplish the desired air change rate inside the chamber of approximately 20 air changes per hour (ACH). The 2 plenums on the upper and lower portions of the chamber utilize a perforated design to uniformly distribute air, promote air mixing, and control the air velocity entering the chamber. Industrial blueprint hanging clamps are being retrofitted to be seated within the chamber drawers and firmly secure the nuclear grade filter media. The aged media will be evaluated using autopsy methodology, and the results obtained intend to help clarify the useful life of nuclear grade HEPA filters. Nuclear grade high efficiency particulate air (HEPA) filters are defined as disposable, extended-media, dry-type filters with a rigid casing enclosing the full depth of the pleats and have a minimum particle removal efficiency of 99.97%. [1] - HEPA filters are constantly exposed to stressors and subsequently have been shown to degrade as the filter ages. Properties such as tensile strength and water repellency, along with others, are negatively impacted. [2] - Due to the scarcity of naturally aged HEPA filters available, it is crucial to develop accelerated aging methods which effectively mimic the natural degradation effects. - After accelerated aging, the filters are subject to testing in the Axial Flow Large Scale Test Stand (ALSTS) and other analysis techniques such as the same qualification tests performed directly after being manufactured. - Quality data is necessary for future decisions regarding the lifetime of nuclear grade HEPA filters. - The objective is to develop a prototype accelerated aging humidity chamber for nuclear grade HEPA filter media. The chamber should expose the media to commonly occurring stressors to expedite the degradation process caused by aging. PID yields reliable and reproducible results: - CFD model correlates well with actual performance; - Further analysis is needed to complete characterization; - Future work for the chamber is still to be completed. The data shown in Figure 5 and Table 1 indicate the system is functioning as designed. As can be determined from the system response evaluation in Table 1, the PID controller is yielding reliable and reproducible results at each set point. As the set point increases, the time constant and 20% settling time increase accordingly. From the results in Table 1, it stands to reason that the CFD model correlates well with the actual performance of the system. Each system response parameter for the 90%RH* CFD model is slightly lower than the parameters collected from the physical model tested at an equivalent set point. This is to be expected considering the inlet in the CFD model is a constant supply of water vapor, whereas the physical model uses a PID controller to monitor the amount of water vapor needed and reacts accordingly. More analysis is still to be done on the prototype chamber. A pressure test will be conducted in order to determine the hourly leak rate of the chamber. This test will involve constant air flow into the chamber until the desired pressure is reached. A uniformity test will also be conducted in order to ensure each filter medium is being exposed to equivalent relative humidity levels. This test will involve an array of humidity meters strategically oriented inside the chamber to record any possible gradients. Future work for the chamber includes dehumidification abilities, tests while fully loaded with filters, and temperature control. Once completed, the autopsy team at Institute of Clean Energy Technology will conduct necessary accelerated aging studies as needed.

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

MULTI-DIMENSIONAL MODELING OF THE CFR ENGINE FOR THE INVESTIGATION OF SI NATURAL GAS COMBUSTION AND CONTROLLED END-GAS AUTOIGNITION

Engine knock and misfire are barriers to pathways leading to high-efficiency Spark-Ignited (SI) Natural Gas engines. The general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by low emission limits and the wide range of chemical reactivity in pipeline quality natural gas. Depending on the region and the source of the natural gas, its reactivity, described by its methane number (analogous to the octane number for liquid SI fuels) can span from 65 - 95. In order to realize diesel-like efficiencies, SI natural gas engines must be designed to operate at high BMEP near knock limits over a wide range of fuel reactivity. This requires a deep understanding regarding the combustion-engine interactions pertaining to flame propagation and end-gas autoignition (EGAI). However, EGAI, if controlled, provides an opportunity to increase SI natural gas engine efficiency by increasing combustion rate and the total burned fuel, mitigating the effects of the slow flame speeds of natural gas fuels which generally reduce BMEP and increase unburned hydrocarbon emissions. For this reason, in order to study EGAI phenomenon, the present work highlights multi-dimensional computational fluid dynamics (CFD) models of the Cooperative Fuel Research (CFR) engine. The CFR engine models are used to investigate fuel-engine interactions that lead to EGAI with natural gas, including effects of fuel reactivity, engine operating parameters, and exhaust gas recirculation (EGR). A Three-Pressure Analysis, performed with GT-Power, was used to estimate initial and boundary conditions for the three-dimensional CFD model. CONVERGE CFD v2.4 was used for the three-dimensional CFD modeling where the level set G-Equation model and SAGE detailed chemical kinetics solver were used. An assessment of the different modeling approaches is also provided to evaluate their limitations, advantages and disadvantages, and for which situations they are most applicable. Model validation was performed with experimental data taken with a CFR engine over varying compression ratio, CA50, EGR fraction, and IMEP and shows good agreement in Peak Cylinder Pressure (PCP), PCP crank angle, and the location of the 10%, 50%, and 90% mass fraction burned (CA10, CA50, and CA90, respectively). The models can predict the onset crank angle and pressure rise rate for light, medium, and heavy EGAI under a variety of fuel reactivities and engine operating conditions.

Bestel, Diego↗
DOE OSTI2023

Assessment of the 3430 Building Filtered Exhaust Stack Sampling Probe Location: Stack Verification Following Fan and Air Blender Additions

The velocity uniformity and flow angle results from the 3430 stack verification tests, performed in April 2023, demonstrated that the CFD model results may be used to support the qualification of the stack sampling location. The measured velocity uniformity verification test result was 2.1 %COV. This value is well within the uniformity criterion, which is that the velocity uniformity be =20 %COV. Additionally, this value is well within the criterion that the actual stack measurement must be within 5% of the surrogate stack (i.e., CFD modeled stack); in this case the CFD modeled average result of 2.85 %COV for the nominal operating range of 22,800 cfm to 62,400 cfm. Additionally, the measured average flow angle at the 3430 stack monitor location was 5.6 degrees. The result is =20 degrees, so the criterion is met.

3430 Building↗
DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗
DOE OSTI2025

A computational investigation of high-flux, plate-and-frame membrane modules for industrial carbon capture

In this work, we study the application of membrane-based separation systems for carbon capture, considering plate-and-frame membrane modules. The successful deployment of membrane CO 2 capture system relies on high-performing membranes as well as effective membrane modules that can fully exploit the developed membranes. A plate-and-frame membrane module is especially attractive for CO 2 capture from industrial flue gas due to its lower pressure drop compared to its counterparts such as spiral wound modules and hollow fiber modules. To design better plate-and-frame modules, we investigate their basic unit - a single membrane stack through a combination of computational modeling and experimental investigations. The modeling approach is based on Computational Fluid Dynamics (CFD) to represent a multiphysics problem, including the fluid flow and diffusion processes within a membrane module. We use experimental data collected under different operating conditions to validate the CFD model. Numerical results suggest a good agreement between experiments and model outputs for the CO 2 recovery, CO 2 mole fraction in the retentate and permeate, and stage-cut. The CFD model is able to predict accurately the flow behavior, providing valuable insights on the effects of fluid dynamics on mass transfer of CO 2 . We also carry out a sensitivity analysis to identify the effect of key parameters on the CO 2 recovery and the CO 2 purity of the outlet streams.

CFD simulation↗
DOE OSTI2022

THERMAL MODELING OF HANFORD LEAD CANISTER’S HEATER BENCH TESTS

A computational fluid dynamics (CFD) model was built to simulate planned testing of heater assemblies for the Hanford Lead Canister (HLC) project. The HLC is a canister storage system that will contain heaters to simulate the decay heat of nuclear material and provide the canister storage system with environmental conditions equivalent to the operating conditions on a dry storage pad. The HLC will be equipped with long-term data collection and monitoring systems to provide an early warning of corrosion, pitting, cracking, or other signs of canister degradation that might threaten the integrity of the containment boundary over the potentially long term of dry storage. An important part of the HLC development is to confirm the function and ability of the electric heater assemblies that were specially designed to provide heating similar to the decay heat of nuclear material contained within the canister storage system. Heater bench testing is planned for early 2022 in a test configuration that does not include the canister. The goal of the bench testing is to verify that the heaters can replicate the decay heat of a canister with nuclear material and to validate the thermal models, which are critical to understanding the HLC’s thermal environment, including the local air flow within the canister storage system. Testing of the heater assemblies inside the canister system are planned in the future to validate canister level thermal models, and rigorous pre-deployment testing of the complete HLC cask and canister system is intended to be completed before the HLC is deployed in the 2025-2026 timeframe. This study presents the pre-test temperature predictions of the bench testing. A description of the heater assemblies and planned bench testing is presented. The model was developed with the commercial CFD code STAR-CCM+. An uncertainty analysis was run with the CFD model to determine the uncertainty in the temperature predictions and provide a range over which the predicted temperatures are expected to vary. The uncertainty analysis was performed by coupling STAR-CCM+ with the software Dakota, which provides advanced parametric analyses, including quantification of margins and uncertainty with computational models. This work is expected to provide insight into SNF canister behavior.

Suffield, Sarah R.↗
DOE OSTI2026

Laser absorption measurements of temperature, pressure, CO, and CO 2 at near-MHz rates in post-detonation fireballs with comparison to synthetic measurements

A laser absorption spectroscopy (LAS) diagnostic was used to obtain measurements of temperature, pressure, CO, and CO 2 at 500 kHz or 1 MHz in post-detonation fireballs produced by hemispherical samples of pentaerythritol tetranitrate (PETN). A quantum-cascade laser was scanned over multiple CO absorption transitions near 2008.5 cm −1 at 1 MHz, while an interband-cascade laser was scanned over a CO 2 absorption transition near 2394.8 cm −1 at 500 kHz. Light from each laser was combined onto a single path and passed through a detonation chamber approximately 83 mm above the 12-mm diameter hemispherical PETN charge. The CO and CO 2 absorption signals were post-processed to obtain time histories of temperature, pressure, species column pressures (P CO L, P CO2 L), and species column mole fractions (X CO L, X CO2 L). Additionally, schlieren imaging was performed simultaneously at 500 kHz to aid interpretation of the LAS measurements. Experimental and synthetic (i.e., CFD based) LAS measurements were compared to evaluate the accuracy of the CFD model and its ability to model the turbulent afterburning of the detonation products in air. In general, the experimental measurements exhibit reasonable agreement with the synthetic measurements at early times; thereby supporting the accuracy of the CFD model. Periods of disagreement between experimental and synthetic measurements at later times are most likely due to a reflected shock and detonator cavity jetting, which are not accounted for in the CFD model.

Schwartz, Charles J. [Purdue Univ., West Lafayette↗
DOE OSTI2022

Assessment of Thermal Stratification in Versatile Test Reactor Transients

The Versatile Test Reactor (VTR) is a fast-spectrum test reactor currently being developed in the United States under the direction of the US Department of Energy. The conceptual design of the 300 MWth pool-type sodium-cooled fast reactor (SFR) has been led by the US National Laboratories in collaboration with General Electric-Hitachi and Bechtel National Inc. Safety performance analysis for the VTR conceptual design is being performed with the systems thermal-hydraulics (SYSTH) module of the SAS4A/SASSYS-1 liquid-metal reactor safety analysis code system. Since the current model of the VTR employs a simple perfect mixing model for large plena like the hot pool, it is not able to predict temperature variations that may develop during the transient. Prior work simulating the response of SFRs to postulated events like the Protected Station Blackout (PSBO) has shown that the phenomenon of thermal stratification, where stable thermal layers accumulate in the hot pool, may delay the transition to natural circulation and thus impact the predicted transient progression. Thus, an effort has begun to model this transient by integrating a Computational Fluid Dynamics (CFD) model of the hot pool into the SAS4A/SASSYS-1 model of the Primary Heat Transport System during the simulation of the PSBO event. A three-dimensional Volume-Of-Fluid CFD model of the VTR hot pool has been developed for the co-simulation of SAS4A/SASSYS-1 with CFD. In this work, the standalone SAS4A/SASSYS-1 calculation and the standalone CFD calculation based on the SAS4A/SASSYS-1 calculation result were produced on the Idaho National Laboratory High Performance Computing cluster, SAWTOOTH. At this time, only the standalone CFD and SAS4A/SASYS-1 simulation results are provided.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

A Numerical Study on the Energy Performance of a Novel Furnace With Acidic Gas Trap Absorbers

Natural gas furnaces are widely used in US residential and commercial building markets. An important issue for natural gas furnaces is serious corrosion and fouling problems caused by acidic gas, such as SOx. An advanced adsorption technology based on acidic gas trap (AGT) absorbers offers the possibility to remove SOx acidic gas from natural gas furnaces with high efficiency and low cost, thereby enabling the development of condensing furnaces without the use of expensive corrosion resistant materials in the heat exchanger. A three-dimensional (3D) computational fluid dynamics (CFD) model has been developed to evaluate the heat transfer performance of a furnace with AGT absorbers and to compare it with a baseline conventional furnace without the AGT. Moreover, an axisymmetric model has been built focusing on the absorbing process in the AGT. The baseline conventional furnace used for the study is a commercial condensing furnace (Rheem 92% AFUE 84,000 BTU Multi-Position Gas Furnace). This furnace was completely disassembled, and the dimensions of each part were carefully measured and used to build a detailed CFD model. A model representing the new furnace, incorporating the AGT absorbers, was developed by adding the AGT system to the conventional furnace model. For the CFD analysis, a mixture model was employed to characterize the heat and mass transfer during the condensing process in the furnace while considering three components—air, water vapor and liquid water. Condensation takes place in the condensing heat exchanger, where water vapor changes phase to liquid water, and the latent heat is thus used in the furnace for useful heating. The simulation results characterize the energy performance of both the conventional furnace and the novel furnace with AGT absorbers, as well as the reactive processing in the AGT. These results provide insightful guidance for the development of the AGT absorber-based furnace from the perspective of its energy performance and will be used to further optimize this novel furnace design.

Laclair, Tim↗
DOE OSTI2023

Assessment of Sodium Thermal Stratification Models Utilizing the TSTF Benchmark

As a result of certain transient scenarios, a thermally stratified layer of liquid sodium can develop in the bulk coolant volumes of a sodium-cooled fast reactor (SFR). In addition to the effects a stratification layer has on the temperature of the heat transport system, a stratification layer can also influence the transition to and establishment of natural circulation flow, which plays an important role in passive cooling and the inherent safety of a pool-type SFR. Therefore, the ability to accurately capture thermal stratification phenomena is important when demonstrating the safety basis of a pool-type SFR during transient sequences. The present work assesses various computational models with different fidelities in their ability to predict thermal stratification in the upper plenum of an SFR. Each computational model will be assessed using the data generated at the Thermal Stratification Test Facility (TSTF) located at the University of Wisconsin-Madison. Using measured flow rate and inlet temperature data, the measured temperature distributions of the tests are compared to the predictions of the lumped volume-based models in SAS4A/SASSYS-1, a 1D-based model in SAM, and a 3-D computational fluid dynamics (CFD) model using STAR-CCM+. The relative performance of the various computational methods is assessed with respect to key metrics such as bulk coolant temperature distribution and plenum exit temperature. A total of eight tests are analyzed, covering different combinations of flow rates (3 and 10 GPM) and upper internal structure (UIS) configurations (none, solid, porous, and open) The perfect mixing model of SAS4A/SASSYS-1 provides the highest accuracy when the flow rate is high and there is no UIS in the test vessel, as high flow rate injection promotes thermal mixing of the sodium in the test vessel. For most of the analyzed tests, the stratified volume model of SAS4A/SASSYS-1 is able to predict the delay in the outlet temperature drop and temperature distribution in the test vessel by a small number of layers to represent thermal stratification. However, the stratified volume model can only simulate a maximum of three temperature layers within a volume and when a layer approaches the elevation of the outlet, the predicted outlet temperature can demonstrate rapid, non-physical changes. The 1-D axial mixing model of SAM provides results that agree reasonably well with the measured data in the prediction of the temporal evolution of the outlet temperature with the exception of the case with a high flow rate and no UIS. The SAM 1-D model has a similar level of accuracy to CFD results when it comes to predicting the outlet temperature. CFD shows overall good agreement in predicting the temperature distribution in the test vessel and outlet temperature. As CFD can model the test vessel geometry in detail, it performs well in the cases of complex geometries such as tests that included a UIS and internal flow through the UIS resulting in active mixing of the coolant in the test vessel. Each of the models discussed in the present work has the potential to be useful during the various stages of reactor design, analysis, and licensing. The lumped-volume approach can be applied for fast turnaround safety calculations to obtain overall reactor behavior during transients. The 1-D models provide improved accuracy when stratification is expected for a relatively low increase in the computational cost. The CFD model can be utilized for confirmatory analysis of the 1-D model, when experimental measurements are not available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗
DOE OSTI2023

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↗
DOE OSTI2022

Sensitivity study of coupled chemical-CFD simulations for analyzing aluminum-clad spent nuclear fuel storage in sealed canisters

We report the United States Department of Energy (DOE) manages over 50 Metric Tons Heavy Metal (MTHM) of aluminum-clad spent nuclear fuel. One main source for DOE’s Aluminum-clad spent nuclear fuel (ASNF) inventory is the advanced test reactor (ATR) at the INL site, which makes this fuel of particular interest for storage scenarios. Road-ready and final disposition packaging configurations for the ATR fuel dictates storage within helium-backfilled, sealed DOE standard canisters. The conditions within these sealed canisters for extended (greater than50 year dry) storage is of interest. To further this goal, a three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model is developed of the sealed DOE standard canisters. This 3D CFD model is one-way coupled with bulk gas radiolysis reactions considering sealed canisters with inert gas and possible trace amounts of air and water vapor. This study looks at the evolution of the thermal history of the canisters over a 50 year time period with a coupling to the chemical reactions occurring from radiolytic breakdown of residual water. A sensitivity study is then carried out over the parameters of the model including the fuel decay heat, residual water content, sealed pressure, canister external temperature, and canister emissivity. In pure helium, hydrogen generation rates are low, under 10 ppm, but hydrogen generation rates are affected greatly by the presence of even 1% residual air, increasing by 50-plus-fold, and nitric acid generation with residual air also occurs ranging from 500 to 4000 ppm after 50 years. The fuel decay heat and the residual water content show the most importance in the generation of hydrogen gas in pure air, and for nitric acid with a residual air condition. External temperature, canister emissivity and sealed pressure all show minor sensitivity effects to the generation of potentially harmful species.

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

Study of effect of PWR cold leg temperature gradient on reactor core condition

Effects of temperature and flow gradients in Westinghouse designed three-loop Pressurized Water Reactor (PWR) cold legs, the piping between the main coolant pump and the reactor vessel, were evaluated using Computational Fluid Dynamics (CFD) code STAR CCM+ and coupled neutronic and thermal-hydraulic (T/H) code system VERA. In the parametric study, several symmetric and asymmetric temperature gradients that were significantly larger than those observed from plant measurements were applied to the cold leg inlets for comparison with the base case without any temperature gradient. A CFD model using the STAR-CCM+ code was developed for a portion of the RCS region between the Reactor Coolant Pump (RCP) and the core inlet based on previously validated modeling approach. The CFD simulation results were processed for the temperature and flow rate distributions at the core inlet as input to the VERA calculations. The VERA code system consists of COBRA-TF (CTF) for thermal-hydraulics, MPACT for reactor physics and neutron transport, and ORIGEN for isotopic depletion. The VERA model was for depletion calculations of a high-burnup loading pattern with the reactor core in pin-by-pin and subchannel resolution. The results of the study indicate that the postulated temperature gradients within the PWR cold legs do not result in any significant changes in the core inlet temperature distributions and the core power distributions during the reactor operation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

A Modeling-Based Flammable Risk Treatment of Refrigerant Leakage from a Commercial R-290 Refrigeration Machine

Because of serious concerns about global warming, manufacturers have started phasing out high global warming potential (GWP) refrigerants in commercial refrigeration equipment (e.g., R-134a). As a potential replacement, propane (R-290) is an environmentally friendly refrigerant for commercial refrigeration equipment because its GWP is only three. However, propane is flammable and is therefore classified as a Class A3 refrigerant per ASHRAE Standards, so safety is a very important consideration when propane-based equipment is designed and deployed in buildings. In the event of a refrigerant leak, flammability of the refrigerant depends on the refrigerant’s local concentration, which is highly affected by the indoor air environment, including temperature and air flow. In this study, a ventilation system attached to a commercial R-290 refrigeration device was designed to eliminate the flammability risk. Moreover, a computational fluid dynamics (CFD) model was developed to investigate the refrigerant leak, thereby evaluating effects of the ventilation system. The CFD model can visualize the flammable zones owing to the leak.

42 ENGINEERING↗
DOE OSTI2022

A whole-core steady-state thermal-hydraulic model for annular fuel type fluoride-salt-cooled reactors

A whole-core, steady-state thermal-hydraulic model is developed for the fluoride-salt-cooled small modular advanced high-temperature reactor (SmAHTR) that employs an annular fuel configuration. This pre-conceptual reactor design by Oak Ridge National Laboratory (ORNL) has the annular fuel and moderator pins arranged in a hexagonal layout. The FLiBe coolant flows from the bottom to the top of the core, parallel to the hexagonal bundle. The fuel and moderator pins in the core are discretized into finite volumes and the 3-D heat conduction equation is solved to obtain the temperature profile. Inter-fuel assembly conduction is also addressed. For this fuel assembly configuration, the coolant flows through two distinct regions – the hexagonal pin bundle and the annulus between the fuel pin and the tie rod. The fluid flow through the hexagonal bundles is modeled using the subchannel approach, in which the coolant region is discretized into corner, edge and interior subchannels and the resulting conservation equations are systematically solved. The 1-D mass, momentum and energy equations are solved for the annulus channels between the fuel pin and the tie rod. Pertinent closure models from the literature are employed to close the system of equations. We also performed a preliminary code-to-code comparison between the present model and a CFD model.. The resulting thermal-hydraulic model can provide temperature, flow rate and pressure drop profiles for the different solid and fluid regions throughout the entire core. Whole-core thermal-hydraulic results for a representative power profile are presented and discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Design of High-Deflection Foils MHK Applications - FEA models

The Ocean Renewable Power Company's (ORPC's) goal is to design, develop, and test hydrofoils with large deflections. The effects of the deflections on cross-flow turbine performance would be evaluated in order to inform design considerations for full-scale water turbines and other marine hydrokinetic devices. Finite element models - NASTRAN files Model scale turbines tested in UNH tow tank Model loads from CFD models

16 TIDAL AND WAVE POWER↗