Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Low Reynolds Number”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Improved natural convection heat transfer correlations for reactor cavity cooling systems of high-temperature gas-cooled reactors: From computational fluid dynamics to Pronghorn

The Reactor Cavity Cooling System (RCCS) is a common reactor safety system in High Temperature Gas Cooled Reactors (HTGR) that removes heat from the Reactor Pressure Vessel (RPV) by radiation ($\sim 80\%$) and natural convection ($\sim 20\%$). For simulation of accident scenarios of HTGRs, intermediate fidelity and system codes models must be employed for limiting the models' execution time. While accurate quantification of the radiative heat transfer is available in these models, quantification of natural convection must rely on correlations of questionable accuracy for the Nusselt number. Commonly used correlations are based in experiments performed at low Rayleigh numbers and/or using isothermal walls in simplified geometries. Here, this work improves on the accuracy of natural convection heat transfer correlations in support for HTGR designs. These correlations include both local and average Nusselt numbers as a function of the global Rayleigh number, the local Rayleigh number, and the temperature profile at the hot wall of the RCCS. In the absence of dedicated experiments and the difficulty of performing high-fidelity simulations at realistic Rayleigh numbers, the data to fit the correlations are generated with Computational Fluid Dynamics (CFD) using Reynolds Averaged Navier-Stokes (RANS) models. First, a careful selection of the RANS turbulence model is performed by comparing the results obtained with different RANS turbulence models against high fidelity simulations of natural convection at $Ra \ 1 \times 10^{11}$ in a rectangular cavity. Next, the selected model is used to perform simulations of an HTGR cavity at different high Rayleigh numbers $\in [6.1 \times 10^{11},2.9 \times 10^{13}]$ to encompass several HTGR designs, assuming an isothermal RPV wall. The results obtained are used to fit a correlation for the average and space-varying Nusselt number as a function of the global and local Rayleigh numbers via a sparsity-promoting least-squares method. The selected RANS model is then used to perform simulations of a PBMR 400 HTGR cavity with the temperature profiles at the RPV wall obtained during a PLOFC transient. We use the results obtained to fit a temperature-dependent correction to the space-varying Nusselt number with the sparsity-promoting least-squares method. The results obtained in this work, enable system-level codes, such as Pronghorn, to perform higher-fidelity simulations of the heat exchange process in the RCCS while still maintaining a low computational cost.

42 ENGINEERING↗
DOE OSTI2022

An improved pressure drop correlation for modeling localized effects in a pebble bed reactor

Advances in the development of pebble bed reactors (PBRs) has created a desire for accurate and cost-effective simulation tools for design scoping studies and safety analysis. The current state-of-the-art for these simulations is the use of porous media models, although these models rely on correlations to capture the effects of flow features that are not explicitly modeled. One of the areas where correlation accuracy is currently lacking is in the near-wall region of the bed. In this region, the presence of the wall causes the pebbles to pack more orderly, drastically changing the geometry and flow behavior in this region. This work presents a new generalized pressure drop correlation for PBRs based on the KTA equation. A high-to-low methodology is applied, where large eddy simulation (LES) is performed on two beds of 1568 and 1700 pebbles to generate a high-fidelity dataset. The flow fields are then averaged in time and separated into concentric rings of 0.05 D peb width. Average porosity, velocity, and pressure drop are extracted for each ring and the friction and form losses are calculated. The Reynolds number range for this study is 625–10,000, and thus the form losses are dominant over the friction losses. The form losses across the rings are investigated, and a correction term for the form loss calculation is determined and applied to the KTA equation to drastically improve the capability of modeling localized porosity effects in a porous media code. Finally, the improved correlation reduces near-wall velocity prediction error from over 50% with the KTA correlation to around 5%. Agreement in pressure drop prediction between LES and porous media simulations is also improved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Hanford Waste Treatment Plant Low Activity Waste Facility Stack Effluent Monitoring - Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant low activity waste (LAW) facility stack monitor locations were qualified using scale model stacks to mitigate the risk of identifying that sampling locations do not meet the qualification criteria on the full-scale stack. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling location were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the LAW stacks were performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing determines the range of conditions for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. Based on these DV values, the corresponding stack flow rates for each of the LAW stacks are 815–55,758 scfm for LV-S1, 980–112,078 scfm for LV-S2, 264–22,901 scfm for LV-S3, and 981–79,832 scfm for LV-C2. The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be =20°. Second, the velocity uniformity at the full-scale stack must be =20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack tests at the LAW facility. Flow angle results were primarily less than 10°, except for one LV-S2 Fan A result, which was 13.2°; all flow angle results were within the =20° criterion. The velocity uniformity results for each test condition ranged between 1.5 COV and 9.2% COV, which were all within the range of the target % COV values from the scale model tests. Based on these stack verification test results, the four LAW filtered exhaust stack sampling locations meet the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all fan operating configurations. This includes single-fan operating conditions for LV-S1 and LV-S2, dual-fan operations for LV-S3 at both the continuous air monitor and record sampler locations, and both the single-fan as well as the dual-fan operations for LV-C2. Further changes to the system configuration or operating conditions that are outside the bounds described in this report may require additional tests or analyses to determine compliance with the standard.

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

Hanford Waste Treatment Plant Low Activity Waste Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation (Rev.1)

The Hanford Tank Waste Treatment and Immobilization Plant Low Activity Waste (LAW) facility stack monitor locations were qualified using scale model stacks to mitigate the risk of discovering that sampling locations do not meet the qualification criteria on the full-scale stacks. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling location were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the LAW stacks were performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing, determines the range of stack flow rates for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. For this analysis, the range of qualified flow rates listed is conservatively based on the average DV through 6 DV for LV-S1, LV-S2, and LV-C2, and 1/3 DV to 3 DV for LV-S3. Table S1 lists the operating flow rates along with the conservative lower and upper qualified stack flow rates for each of the LAW facility stacks. For each stack, the operating flow is below the upper qualified stack.

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

Mycoparasites, Gut Dwellers, and Saprotrophs: Phylogenomic Reconstructions and Comparative Analyses of Kickxellomycotina Fungi

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

TEAMER: Experimental performance characterization of a shrouded axial-flow turbine

Sitkana has developed a shrouded hydrokinetic turbine with a modular, low-cost design that can be scaled to meet the needs of remote communities. With technical support from the University of Washington, Sitkana sought to experimentally characterize the mechanical power and structural loads of various 1:3.3 scale rotor geometries. In all, 11 different rotor geometries were characterized with variations in height-to-diameter ratio, blade number, and blade type (foiled versus flat). All tests were conducted in Reynolds-independent flow conditions in the Alice C. Tyler Flume at the University of Washington. Results allow Sitkana to (1) refine the optimal rotor geometry, (2) validate numerical models, and (3) predict power output for a full-scale system. This project is part of the TEAMER RFTS 8 (request for technical support) program.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Foaming prediction in pure liquids from dimensionless numbers inspired by the theory of fluid behavior for drops

Foaming prediction is critical for selecting materials and designing processes in industries such as bioprocessing and gas processing. Existing models lack the generality needed for a wide range of materials and overlook the foaming behavior in pure liquids. Here, this work presents a novel method for predicting foaming in pure liquids based on their density, surface tension, and viscosity, using Reynolds ( Re ) and Ohnesorge ( Oh ) numbers. A foaming prediction map, leveraging the theory of fluid drop behavior, was developed by plotting these numbers. This map delineates distinct non-foaming and foaming regions, functioning as a binary classifier for foaming predictions. The map was fitted and validated through shake test experiments on 46 liquids, demonstrating reliable predictions, except for a specific region characterized by small Oh and large Re numbers. This region corresponded to relatively low foam stability and high turbulence, making foaming predictions challenging for liquids in this category.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Numerical modeling of hydrogen mixing in a direct-injection engine fueled with gaseous hydrogen

Hydrogen is considered as one of the most promising options to achieve effective decarbonization of the energy and transportation sectors. As such, it has recently been receiving increasing attention because of its promising potential as an energy carrier for advanced energy and propulsion systems. With a focus on internal combustion engines, direct injection (DI) of gaseous hydrogen during the compression stroke offers great potential for high engine efficiency and specific power while reducing the risk of backfiring and pre-ignition issues. Therefore, many experimental and numerical efforts have recently been dedicated to understanding the physical and chemical behaviors of hydrogen in engine during mixing and combustion. This study focuses on computational fluid dynamics (CFD) modeling of the hydrogen DI process in a hydrogen optical research engine. Under the conditions studied, gaseous hydrogen is injected into the combustion chamber via a centrally located single-hole injector at a pressure of 100 bar. Two configurations, namely low-and high-tumble, are investigated to understand the impact of different in-cylinder flow patterns on the fuel-air mixture preparation. Simulations are carried out using the commercial CFD software CONVERGE. Here, the in-cylinder turbulence is modeled with an unsteady Reynolds-averaged Navier-Stokes (URANS) formulation closed by the renormalization group (RNG) k-ε model. Several numerical methods and model constants, including but not limited to turbulent Schmidt number, are evaluated. The numerical results are systematically compared against experimental measurements of velocity and hydrogen concentration fields on the vertical center plane to assess the performance of the CFD model, unveil the physics of hydrogen mixing, and establish best practices for modeling hydrogen DI under relatively high injection pressure conditions.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Hydrokinetic Power Conversion Using Vortex-Induced Oscillation with Cubic Restoring Force

A cubic-spring restoring function with high-deformation stiffening is introduced to passively improve the harnessed marine hydrokinetic power by using flow-induced oscillations/vibrations (FIO/V) of a cylinder. In these FIO/V experiments, a smooth, rigid, single-cylinder on elastic end-supports is tested at Reynolds numbers ranging from 24,000 < Re < 120,000. The parameters of the tested current energy converter (CEC) are cubic stiffness and linear damping. Using the second generation of digital virtual spring-damping (Vck) controller developed by the Marine Renewable Energy Laboratory (MRELab), the cubic modeling of the oscillator stiffness is tested. Experimental results show the influence of the parameter variation on the amplitude, frequency, energy conversion, energy efficiency, and power of the converter. All experiments are conducted in the low turbulence-free surface water (LTFSW) channel of the MRELab of the University of Michigan. The main conclusions are: (1) The nonlinearity in the cubic oscillator is an effective way to extend the vortex-induced vibration (VIV) upper branch, which results in higher harnessing power and efficiency compared to the linear stiffness cylinder converter. (2) Compared to the linear converter, the overall power increase is substantial. The nonlinear power optimum, occurring at the end of the VIV upper branch, is 63% higher than its linear counterpart. (3) The cubic stiffness converter with low harnessing damping achieves consistently good performance in all the VIV regions because of the hardening restoring force, especially at higher flow velocity.

13 HYDRO ENERGY↗
DOE OSTI2024

Preliminary Assessment of Additively Manufactured Cooling Channel Performance for Helium-Cooled Blanket Concepts

Sufficient cooling of the blanket first wall remains a critical challenge for the design and deployment of fusion power plants. Helium has been targeted as a potential blanket coolant due to its inertness and low neutron interactivity, among other advantages. However, the low thermal mass of helium creates a need for heat transfer enhancements in coolant channels to provide adequate cooling to the blanket’s first wall. Toward this end, a series of ribbed flow channels of various rib cross sections and configurations has been produced via additive manufacturing (AM) to study the efficacy of AM for first wall heat transfer enhancement and the optimization of heat transfer geometries. Helium cooling performance is studied in AM test articles at 4 MPa operating pressure, Reynolds numbers up to 197 000, and outer surface heat fluxes up to 42 kW/m 2 in the recently commissioned helium flow loop experiment (HFLE). Preliminary results of this study are presented herein. Heat transfer performance of nominally smooth (i.e., featureless) AM channels is quantified via measured Nusselt numbers and friction factors and compared to off-the-shelf smooth pipe experiments. Results are compared to existing correlations and used to discuss the effects of the AM processes on thermal-hydraulic performance. It is seen that the inherent roughness of the AM channels leads to an increase in both heat transfer coefficient and pressure drop when compared to the conventional pipe. In conclusion, recommendations are made for future studies based on these findings and additional considerations for the deployment of AM blanket cooling components.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Time-resolved particle image velocimetry measurements of the turbulent Richtmyer–Meshkov instability

Experiments are presented on the Richtmyer–Meshkov instability (RMI) with a three-dimensional, multi-mode initial perturbation. The experiments use a vertical shock tube, where a stably stratified interface is formed between air and sulphur hexafluoride (SF 6 ) via counterflow. A perturbation is imposed at the interface by vertical oscillation of the gas column, forming Faraday waves. The interface is accelerated by a Mach 1.17 (in air) shock wave, and the development of the mixing region between the gases is investigated using particle image velocimetry. Following shock acceleration, a reflected shock wave from the bottom of the shock tube interacts with the mixing layer a second time (reshock). The experiment is initialized with both high and low amplitude perturbations to examine the effect of the perturbation amplitude on measured quantities. The instability growth exponent (θ) is determined from the kinetic energy field using the width of the mixing layer and the decay of kinetic energy, which are found to be in agreement when the flow is most strongly excited. A growth exponent of θ ≈ 0.5 is found for all cases except the high-amplitude reshocked regime (where θ ≈ 0.33). High-amplitude experiments exhibit the transitional outer Reynolds number (Re ≡ hh˙/ν >10 4 ) required for mixing transition following the incident shock, and both experiments are elevated well above this threshold following reshock. However, neither set of experiments meet the more stringent requirements proposed by Zhou et al. (Phys. Rev.E, vol. 67, issue 5, 2003) which include the time dependent aspect of the RMI, an observation which is also made when examining the spectra.

42 ENGINEERING↗
DOE OSTI2024

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Measurement of Convective Heat Transfer Coefficients with Supercritical CO2 in Novel Additively Manufactured Helically Patterned Pin Fin Tubes Using the Wilson Plot Technique

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO2 flowing in pin-fin patterned pipes in the Heat Exchange and Experimental Testing (HEET) facility at the US DoE’s National Energy Technology Laboratory (NETL) in Morgantown, WV. The measurement procedures in the HEET rig were validated by conducting benchmark tests with smooth stainless-steel tube and comparing the results with published correlations for Nusselt number (Nu) and friction factor. Over typical Reynolds number range in sCO2, the measured Nu and friction factors were within 7% of classical correlations for smooth tube flow.The candidate pin fin patterned pipes were additively manufactured (AM) at the Oak Ridge National Laboratory. The pins were circular or elliptical in cross-section. Pin length to diameter aspect ratios were 1.33 and 2, while the pin diameter to tube diameter ratio was 0.188 and 0.125. Tests were performed for ReD varying from 6.9×104 to 2.2×105 and at conditions equivalent to the low pressure outlet (8.69 MPa, 361 K) of the low temperature recuperator (LTR) in an indirect sCO2 power cycle. The Wilson plot technique was utilized to measure the bulk heat transfer coefficients.For the better performing design (tube A, pin length to tube diameter ratio: 1.33, pin diameter to tube diameter ratio: 0.188), the local heat transfer coefficient increased by 112% relative to the Dittus-Boelter correlation at the LTR low pressure outlet. This corresponded to a 282% increase in the product of the heat transfer coefficient and the surface area. Large pressure drops across the test articles were observed.

Searle, Matthew↗
DOE OSTI2020

Effect of the numerical discretization scheme in Shock-Driven turbulent mixing simulations

In this work, we evaluate the effects of distinct numerical strategies in simulations of shock-driven turbulent mixing. An air-SF 6 -air gas curtain subjected to Mach 1.2 shock-waves is computed with Implicit Large-Eddy Simulation based on three numerical schemes: directional-split with Harten-Lax-van Leer (HLL) solver; directional-unsplit using HLL-Contact (HLLC) solver; and directional-unsplit with HLLC solver and a Low Mach number Correction (LMC). The results illustrate the importance of the numerical strategy to the accuracy of the predictions. Whereas both split and unsplit schemes result in a similar spatial development of the initial shock-driven instability, only the unsplit schemes can predict the turbulent mixing transition after reshock observed by the laboratory experiments. Such feature increases the mixing rate of the two fluids, this being particularly pronounced when the LMC is active due to i) the reduced flow Mach number, and ii) the larger effective Reynolds number. Since the selected mixing problem is driven by the deposition of vorticity at the fluids’ interface, the resultant flow physics is analyzed by investigating the contribution of distinct inviscid mechanisms to the production of vorticity: baroclinicity, stretching, and dilatation. As expected, it is observed that the production of vorticity is initially dominated by the baroclinicity mechanism. Yet, the relevance of the remainder mechanisms is enhanced after the reshock and may even surpass the baroclinicity term.

42 ENGINEERING↗
DOE OSTI2022

(Towards) DNS of a Laboratory Lean CH4/H2 Low-Swirl Flame Impinging on an Inclined Wall

Due to downsizing trends, flame-wall interaction (FWI) is increasingly prominent in gas turbines (GTs). FWI has direct consequences on flame stabilization and pollutant emissions, but it is not well understood in turbulent flows representative of GTs. We present results from a direct numerical simulation (DNS) of a turbulent CH4/H2 model GT low-swirl laboratory-scale flame interacting with an inclined wall. The results from the laboratory flame include simultaneous measurements of velocity using stereo particle imaging velocimetry and OHxCH2O planar laser induced fluorescence. The adaptive-mesh refinement solver PeleLMeX is used, with 24-species, 105-reaction reduced Aramco chemical kinetics mechanism. The premixed fuel-air mixture consists of hydrogen-enriched methane with 70% hydrogen volume fraction and 0.4 equivalence ratio. The inflow is prescribed to match experimental measurements at the burner exit. Karlovitz and turbulent Reynolds numbers are 300 and 400, respectively. The simulation and experimental results show excellent agreement. The flame features a bowl-shape stabilization, with a corrugated, continuous flame front at the leading edge, followed by fragmented reaction zones downstream. A large diffuse cloud of CH2O is formed downstream of the quenching point. The simulation results indicate that the cloud of CH2O is the result of incomplete methane combustion, with CH2O "leaking" from the locally quenched reaction zones.The DNS provides fine-grain resolution of turbulence-flame-wall interaction that cannot be captured with experimental measurements. With access to the entire solution vector at each cell of the computational domain, the local quenching.

flame-wall interactions↗
DOE OSTI2021

Aero-Thermal Characterization of Accelerating and Diffusing Passages Downstream of Rotating Detonation Combustors

Cycle benefits of rotating detonation engines show up to five percentage points of efficiency gain for low-pressure ratio engines. An optimal integration between the combustor and the turbine needs to be guaranteed to realize this potential gain. The rotating detonation combustor (RDC) exhausts transonic flow with shocks rotating at frequencies ranging from a few to tens of kilohertz depending on the number of present waves. Hence, the turbine design requires precise knowledge of the fluctuations and losses downstream of the combustor. This paper focuses on the quantification of fluctuations and losses for accelerating and diffusing passages. The analysis of the combustor is performed via reactive unsteady Reynolds Averaged Navier-Stokes (URANS) simulations. The unsteady RANS equations are solved via CFD++ from Metacomp with a one-step reaction mechanism for an H2-air mixture. The resolving of the boundary layer is achieved with a structured mesh of around 36 million cells. Inlet pressure of 10 bar and two different back pressures are applied to the combustor to model the interconnection with downstream turbines. Finally, we present and assess a methodology to reduce the computational time to model these passages ten times.

Braun, James↗
DOE OSTI2025

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗