Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 91 records · Page 5

Performance and Wake Characterization of a Model Hydrokinetic Turbine: The Reference Model 1 (RM1) Dual Rotor Tidal Energy Converter

The mechanical power and wake flow field of a 1:40 scale model of the US Department of Energy’s Reference Model 1 (RM1) dual rotor tidal energy converter are characterized in an open-channel flume to evaluate power performance and wake flow recovery. The NACA-63(4)-24 hydrofoil profile in the original RM1 design is replaced with a NACA-4415 profile to minimize the Reynolds dependency of lift and drag characteristics at the test chord Reynolds number. Precise blade angular position and torque measurements were synchronized with three acoustic Doppler velocimeters (ADV) aligned with each rotor centerline and the midpoint between the rotor axes. Flow conditions for each case were controlled to maintain a hub height velocity, uhub= 1.04 ms−1, a flow Reynolds number, ReD= 4.4 × 105, and a blade chord length Reynolds number, Rec= 3.1 × 105. Performance was measured for a range of tip-speed ratios by varying rotor angular velocity. Peak power coefficients, CP= 0.48 (right rotor) and CP= 0.43 (left rotor), were observed at a tip speed ratio, λ= 5.1. Vertical velocity profiles collected in the wake of each rotor between 1 and 10 rotor diameters are used to estimate the turbulent flow recovery in the wake, as well as the interaction of the counter-rotating rotor wakes. The observed performance characteristics of the dual rotor configuration in the present study are found to be similar to those for single rotor investigations in other studies. Similarities between dual and single rotor far-wake characteristics are also observed.

30 DIRECT ENERGY CONVERSION↗

On the enhancement of boundary layer skin friction by turbulence: an angular momentum approach

Turbulence enhances the wall shear stress in boundary layers, significantly increasing the drag on streamlined bodies. Other flow features such as free stream pressure gradients and streamwise boundary layer growth also strongly influence the local skin friction. In this paper, an angular momentum integral (AMI) equation is introduced to quantify these effects by representing them as torques that alter the shape of the mean velocity profile. This approach uniquely isolates the skin friction of a Blasius boundary layer in a single term that depends only on the Reynolds number most relevant to the flow's engineering context, so that other torques are interpreted as augmentations relative to the laminar case having the same Reynolds number. The AMI equation for external flows shares this key property with the so-called FIK relation for internal flows (Fukagata et al. , Phys. Fluids , vol. 14, 2002, pp. L73–L76). Without a geometrically imposed boundary layer thickness, the length scale in the Reynolds number for the AMI equation may be chosen freely. After a brief demonstration using Falkner–Skan boundary layers, the AMI equation is applied as a diagnostic tool on four transitional and turbulent boundary layer direct numerical simulation datasets. Regions of negative wall-normal velocity are shown to play a key role in limiting the peak skin friction during the late stages of transition, and the relative strengths of terms in the AMI equation become independent of the transition mechanism a very short distance into the fully turbulent regime. The AMI equation establishes an intuitive, extensible framework for interpreting the impact of turbulence and flow control strategies on boundary layer skin friction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Characterization of the fractal dimension of helium plumes in laminar to turbulent transition

Abstract Schlieren imaging in conjunction with a high-speed camera was used to observe the behavior of metered helium plumes as they transition from laminar to turbulent flow in an air environment. The plumes were visualized at twelve jet Reynolds numbers ranging from 200 to 2980. The fractal dimension of the flows was obtained by applying a box counting algorithm to the recorded schlieren images. The results were analyzed to determine the correlation between the Reynolds number of the flow and the fractal dimension of the observed turbulence. A trend of increasing fractal dimension with increasing Reynolds number was observed for several different types of schlieren cutoffs including horizontal cutoff, vertical cutoff, circular cutoff, focused shadowgraphy and de-focused shadowgraphy. The vertical cutoff and focused shadowgraphy imaging methods showed the most consistent results for the fractal dimension characterization during the laminar to turbulent transition. For transitional plumes, it was observed that fractal dimension increased with distance from the jet outlet.

42 ENGINEERING↗

Bifurcation of equilibrium positions for ellipsoidal particles in inertial shear flows between two walls

We conducted a systematic numerical investigation of spherical, prolate and oblate particles in an inertial shear flow between two parallel walls, using smoothed particle hydrodynamics (SPH). It was previously shown that above a critical Reynolds number, spherical particles experience a supercritical pitchfork bifurcation of the equilibrium position in shear flow between two parallel walls, namely that the central equilibrium position becomes unstable, leading to the emergence of two new off-centre stable positions (Foxet al.,J. Fluid Mech., vol. 915, 2021). This phenomenon was unexpected given the symmetry of the system. In addition to confirming this finding, we found, surprisingly, that ellipsoidal particles can also return to the centre position from the off-centre positions when the particle Reynolds number is further increased, while spherical particles become unstable under this increased Reynolds number. By utilizing both SPH and the finite element method for flow visualization, we explained the underlining mechanism of this reverse of bifurcation by altered streamwise vorticity and symmetry breaking of pressure. Furthermore, we expanded our investigation to include asymmetric particles, a novel aspect that had not been previously modelled, and we observed similar trends in particle dynamics for both symmetric and asymmetric ellipsoidal particles. While further validation through laboratory experiments is necessary, our research paves the road for development of new focusing and separation methods for shaped particles.

Mechanics↗

Experimental Study on Flow Behavior and Heat Transfer Enhancement With Distinct Dimpled Gas Turbine Blade Internal Cooling Channel

Abstract The study presents the investigation on heat transfer distribution along a gas turbine blade internal cooling channel. Six different cases were considered in this study, using the smooth surface channel as a baseline. Three different dimples depth-to-diameter ratios with 0.1, 0.25, and 0.50 were considered. Different combinations of partial spherical and leaf dimples were also studied with the Reynolds numbers of 6,000, 20,000, 30,000, 40,000, and 50,000. In addition to the experimental investigation, the numerical study was conducted using a large eddy simulation (LES) to validate the data. It was found that the highest depth-to-diameter ratio showed the highest heat transfer rate. However, the deepest dimpled channel shows the highest pressure drop, which affects the overall thermal performance of the cooling channel. The results showed that the leaf dimpled surface is the best cooling channel based on the highest Reynolds number's heat transfer enhancement and friction factor. However, at the lowest Reynolds number, partial spherical dimples with a 0.25 depth-to-diameter ratio showed the highest thermal performance.

Energy & Fuels↗

The Madison plasma dynamo experiment: A facility for studying laboratory plasma astrophysics

The Madison plasma dynamo experiment (MPDX) is a novel, versatile, basic plasma research device designed to investigate flow driven magnetohydrodynamic instabilities and other high-β phenomena with astrophysically relevant parameters. A 3 m diameter vacuum vessel is lined with 36 rings of alternately oriented 4000 G samarium cobalt magnets, which create an axisymmetric multicusp that contains ∼14 m3of nearly magnetic field free plasma that is well confined and highly ionized (>50%). At present, 8 lanthanum hexaboride (LaB6) cathodes and 10 molybdenum anodes are inserted into the vessel and biased up to 500 V, drawing 40 A each cathode, ionizing a low pressure Ar or He fill gas and heating it. Up to 100 kW of electron cyclotron heating power is planned for additional electron heating. The LaB6cathodes are positioned in the magnetized edge to drive toroidal rotation through J × B torques that propagate into the unmagnetized core plasma. Dynamo studies on MPDX require a high magnetic Reynolds number Rm > 1000, and an adjustable fluid Reynolds number 10 1). Initial results from MPDX are presented along with a 0-dimensional power and particle balance model to predict the viscosity and resistivity to achieve dynamo action.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PERFORMANCE CHARACTERIZATION OF THE DIRECT CONTACT MEMBRANE DISTILLATION MODULE OPERATED WITH HIGHLY CONCENTRATED SOLUTION

The experimental study was conducted to characterize the performance of the flat sheet Direct Contact Membrane Distillation. The length degradation was studied in a short and long module with and without spacers. The experimental study is validated using the predicted flux for a wide range of operating conditions using deionized water. Three different membranes are used in the parametric study for feed Reynolds number varied between 500 and 2000, and inlet feed temperature varied between 45°C and 80°C. PTFE membrane was shown to be more effective than the PVDF due to its higher permeability coefficient and higher hydrophobicity. Fouling study was conducted on the selected PTFE membrane with various feed Reynolds numbers between 500 and 2000 and feed inlet temperature between 50°C and 80°C. Membranes used in the fouling study were characterized using scanning electron microscopy (SEM) and X-ray photoelectron spectroscopy (XPS) for the morphology and chemical composition of the fouling layer. 36 hours of operation with 10 wt% NaCl (100,000 ppm) at feed inlet temperatures 80°C and 60°C and fixed Reynolds 2000 demonstrated a drop in flux by 31% and 27%, respectively. Low flow rate and high temperature are favorable conditions for fouling formation and reduced permeate flux for highly concentrated solutions. The spacers improved the flux by 30.8% and 27.6% at the operating Reynolds numbers 2000 the short and long channel modules. This study helps to design and optimize DCMD modules for industrial-size desalination and water treatment applications.

Kim, Deliya↗

Experimental and computational investigation of interactive flow induced oscillations of two tandem rough cylinders at 3×10 4 ≤Re≤1.2×10 5

The interactive flow induced oscillations (FIO) of two adjacent, elastically mounted, rigid, tandem, locally-rough cylinders in transverse flow direction are analyzed utilizing two-dimensional Unsteady Reynolds-Averaged Navier-Stokes (2-D URANS) equations and verified experimentally in the proximity-wake region. Three sets of experiments and simulations (K = 600N/m, ζ = 0.14; K = 600N/m, ζ = 0.26; K = 1,200N/m, ζ = 0.26; K is spring stiffness, ζ is damping ratio) of two tandem cylinders with turbulence stimulation are tested and simulated for Reynolds number in the range of 30,000 = Re ≤ 120,000. The reduced velocity range is U* = 2.48–14.22, the mass ratio is m* = 1.343, and the center-to-center in-flow spacing to diameter ratio is d/D = 2.57. The characteristics of amplitude response, frequency response, lift force, and interactive wake patterns are presented and discussed. The trends of the amplitude and frequency responses from numerical simulations are in good agreement with experimental results. The main conclusions of Reynolds number effect on interactive flow induced oscillations are: (1) Five significant flow patterns between two tandem locally-rough cylinders for different Reynolds numbers are observed through analyzing the complex but stable interactions between vortices and cylinders. (2) In the initial and upper VIV branches, the downstream cylinder's FIO is seriously interfered by the wake of the upstream cylinder. (3) The downstream cylinder is strongly impinged by the vortices shed from the upstream cylinder resulting in nearly 180O out-of-phase oscillations in transition from VIV to galloping and in-phase oscillations in galloping.

42 ENGINEERING↗

Shear-induced lift force on spheres in a viscous linear shear flow at finite volume fractions

Several studies have shown a significant increase in drag on a distribution of solid spherical particles within a fluid with increasing particle volume fraction. As a result, many empirical drag laws accounting for the dependence on the Reynolds number and volume fraction can be found in the literature. This study investigates the possibility of a similar effect of the particle volume fraction on the mean hydrodynamic lift force on randomly distributed spherical particles in a linear shear flow. Particle-resolved direct numerical simulations are performed to evaluate the mean lift force, and the results are compared with the case of an isolated particle in a linear shear flow for the same Reynolds number and shear rate. The mean lift force acting on the particles appears to remain nearly the same as that on an isolated particle. However, due to the influence of neighboring particles, there is a substantial force variation in transverse directions on each individual particle, whose magnitude is comparable to the mean drag force. The distribution of drag force in a linear shear flow is shown to be nearly the same as in a uniform flow at the same volume fraction and Reynolds number. A simple stochastic model based on a Gaussian distribution is presented for the lift force variation, and its performance is compared to the prediction of the deterministic pairwise interaction extended point-particle model.

42 ENGINEERING↗

Critical Shear Stress for Erosion Under Laminar Jet Flow

Particle erosion is an important process in both natural and manmade environments. For example, in natural environments particle or sediment erosion rates determine the course of rivers, depth of streams, and fate of deltas. The interplay between erosion rates and flow patterns can have dramatic impact on shorelines. Where hydrodynamic jet flows or rip currents change direction, sediments may accumulate or disperse altering the shape of coastlines and the ecosystems and economies that rely on them. In man-made systems, jet flows may be used to mobilize particles from surfaces or within industrial scale mixing tanks. For example, in the nuclear industry, impinging jet mixers may be used to mobilize and suspend sludges at the bottom of waste tanks. In each of these examples, the Shields diagram is useful to determine the conditions under which erosion occurs. This diagram plots the scaled critical shear stress for erosion as a function of the particle size as embedded in particle Reynolds numbers, Archimedes numbers, or their equivalents. Although the data, that carries a large spread, is largely in agreement with mathematical models, for the smallest particles, model predictions diverge from each other and from experimental data. Indeed, there is substantial uncertainty about the Shields diagrams at relatively small particle sizes and flow rates. Here we evaluate critical shear stresses for erosion under laminar jet flow. To date critical shear stresses for erosion, as embodied in Shields diagrams, have largely been developed for turbulent flow both experimental and mathematically. However, curves on the Shields diagrams diverge among various models at lower particle Reynolds numbers. Indeed, how critical shear stresses for erosion develop under fully laminar jet flow conditions (not simply laminar boundary layers under turbulent flow) remains unclear. Here we address this gap, which is particularly important for mobilization of the smallest particles. We find the Shields parameter to be inversely proportional to the particle Reynolds numbers under laminar jet flow conditions similar to turbulent conditions.

critical shear stresses, erosion rates, effective ↗

Large-Eddy Simulation of Flow Over Boeing Gaussian Bump Using Multiagent Reinforcement Learning Wall Model: Preprint

We develop a wall model for large-eddy simulation (LES) that takes into account various pressure-gradient effects using multi-agent reinforcement learning. The model is trained using low-Reynolds-number flow over periodic hills with agents distributed on the wall at various computational grid points. It utilizes a wall eddy-viscosity formulation as the boundary condition to apply the modeled wall shear stress. Each agent receives states based on local instantaneous flow quantities at an off-wall location, computes a reward based on the estimated wall-shear stress, and provides an action to update the wall eddy viscosity at each time step. The trained wall model is validated in wall-modeled LES of flow over periodic hills at higher Reynolds numbers, and the results show the effectiveness of the model on flow with pressure gradients. The analysis of the trained model indicates that the model is capable of distinguishing between the various pressure gradient regimes present in the flow. To further assess the robustness of the developed wall model, simulations of flow over the Boeing Gaussian bump are conducted at a Reynolds number of 2 x 10^6, based on the free-stream velocity and the bump width. The results of mean skin friction and pressure on the bump surface, as well as the velocity statistics of the flow field, are compared to those obtained from equilibrium wall model (EQWM) simulations and published experimental data sets. The developed wall model is found to successfully capture the acceleration and deceleration of the turbulent boundary layer on the bump surface, providing better predictions of skin friction near the bump peak and exhibiting comparable performance to the EQWM with respect to the wall pressure and velocity field. We also conclude that the subgrid-scale model is crucial to the accurate prediction of the flow field, in particular the prediction of separation.

boundary layer↗

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation) [SWR-26-095]

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation): Multifidelity aerodynamic polar data generation for hydrofoil/tidal-turbine airfoil sections. Foilpolars ties together three pieces: *AeroSandbox supplies the baseline airfoil coordinates (UIUC database). *G2Aero parameterizes those shapes on a Grassmannian manifold (Karcher mean + PGA basis) and samples new perturbed shapes around that basis. *XFoil (panel method) and NeuralFoil (neural-network surrogate, shipped with AeroSandbox) each solve the resulting shapes for lift, drag, moment, and pressure at the swept angles of attack, Reynolds numbers, and n_crit values. Design optimization of foil shapes in a computationally efficient way requires polars data across many candidate shapes, not just a handful of baseline foils. However, high-fidelity CFD at that scale is too costly, and naive shape perturbation strays from realistic geometries. FOILPOLARS addresses this by loading baseline airfoils (via AeroSandbox) and mapping them onto a Grassmannian manifold (via G2Aero), computing a Karcher mean and principal geodesic analysis (PGA) basis. New shapes are sampled by perturbing PGA coefficients, keeping them close to the manifold of realistic foils. Each sampled shape is evaluated across a configurable sweep of angle of attack, Reynolds number, and critical amplification factor using two solvers: XFoil (panel method) and NeuralFoil (neural-network surrogate), producing a paired dataset of lift, drag, moment, pressure, convergence, and confidence, indexed alongside each shape's PGA coefficients and shared Grassmannian basis in a single xarray dataset. From this, FOILPOLARS produces convergence summaries and comparison plots per shape, Reynolds number, and n_crit. A command-line interface exposes each pipeline stage independently, supporting data-driven design, optimization, and machine-learning workflows for foils.

Sandhu, Rimple [National Laboratory of the Rockies↗

Reynolds stress scaling in the near-wall region of wall-bounded flows

A new scaling is derived that yields a Reynolds-number-independent profile for all components of the Reynolds stress in the near-wall region of wall-bounded flows, including channel, pipe and boundary layer flows. The scaling demonstrates the important role played by the wall shear stress fluctuations and how the large eddies determine the Reynolds number dependence of the near-wall turbulence behaviour.

42 ENGINEERING↗

Near-wall model for compressible turbulent boundary layers based on an inverse velocity transformation

In this work, a near-wall model, which couples the inverse of a recently developed compressible velocity transformation (Griffin et al., Proc. Natl Acad. Sci., vol. 118, 2021, p. 34) and an algebraic temperature–velocity relation, is developed for high-speed turbulent boundary layers. As input, the model requires the mean flow state at one wall-normal height in the inner layer of the boundary layer and at the boundary-layer edge. As output, the model can predict mean temperature and velocity profiles across the entire inner layer, as well as the wall shear stress and heat flux. The model is tested in an a priori sense using a wide database of direct numerical simulation high-Mach-number turbulent channel flows, pipe flows and boundary layers (48 cases, with edge Mach numbers in the range 0.77–11, and semi-local friction Reynolds numbers in the range 170–5700). The present model is significantly more accurate than the classical ordinary differential equation (ODE) model for all cases tested. The model is deployed as a wall model for large-eddy simulations in channel flows with bulk Mach numbers in the range 0.7–4 and friction Reynolds numbers in the range 320–1800. When compared to the classical framework, in the a posteriori sense, the present method greatly improves the predicted heat flux, wall stress, and temperature and velocity profiles, especially in cases with strong heat transfer. In addition, the present model solves one ODE instead of two, and has a computational cost and implementation complexity similar to that of the commonly used ODE model.

42 ENGINEERING↗

Optimization of a High Pressure Turbine Blade and Sector-Based Annular Rig Design for Supercritical CO2 Power Cycle Representative Testing

Abstract As part of the ongoing research into the design of hardware for zero emission cycles, a first-stage high-pressure turbine (HPT) blade is optimized for a 300 MWe supercritical CO2 (sCO2) power cycle using the surrogate-assisted genetic algorithm optimizer in Numeca FINE/Design three-dimensional with objectives of increasing efficiency and decreasing heat load to the blade. Supercritical CO2 property tables are constructed from NIST REFPROP data for the condensable gas simulation in FINE/Turbo. A detailed mesh sensitivity study is performed for a baseline design to identify the proper-grid refinement and efficiently allocate resources for the optimization. Seventy design variables are selected for the initial population generation. Self-organizing maps are then used to focus the design variables on the most important ones affecting the objective functions. The optimization results in approximately 3000 three-dimensional Reynolds Averaged Navier Stokes simulations of different blade shapes with increases in efficiency of up to 0.85% and decreases in heat load of 14%. Families of blade shapes are identified for experimental testing in an annular rig at the Purdue Experimental Turbine Aerothermal Laboratory. A design to adapt the annular cascade for testing optimized geometries is introduced, which features eccentric radius sectors allowing for scaled-up geometries of sCO2 optimized blade profiles to be tested at design cycle representative conditions at high Reynolds numbers in dry air. Analysis into the effects of Reynolds number, working fluid, and geometric relations are presented to prove the efficacy of the test method.

Engineering↗

Physics-constrained deep learning of nonlinear normal modes of spatiotemporal fluid flow dynamics

In this study, we present a physics-constrained deep learning method to discover and visualize from data the invariant nonlinear normal modes (NNMs) which contain the spatiotemporal dynamics of the fluid flow potentially containing strong nonlinearity. Specifically, we develop a NNM-physics-constrained convolutional autoencoder (NNM-CNN-AE) integrated with a multi-temporal-step dynamics prediction block to learn the nonlinear modal transformation, the NNMs containing the spatiotemporal dynamics of the flow, and reduced-order reconstruction and long-time future-state prediction of the flow fields, simultaneously. In test cases, we apply the developed method to analyze different flow regimes past a cylinder, including laminar flows with low Reynolds number in transient and steady states (RD = 100) and high Reynolds number flow (RD = 1000), respectively. The results indicate that the identified NNMs are able to reveal the nonlinear spatiotemporal dynamics of these flows, and the NNMs-based reduced-order modeling consistently achieves better accuracy with orders of magnitudes smaller errors in construction and prediction of the nonlinear velocity and vorticity fields, compared to the linear proper orthogonal decomposition (POD) method and the Koopman-constrained-CNN-AE using the same number or dimension of modes. We perform an analysis of the modal energy distribution of NNMs and find that compared to POD modes, the few fundamental NNMs capture a very high level of total energy of the flow, which is advantageous for reduced-order modeling and representation of the complex flows. Finally, we discuss the potentials and limitations of the presented method.

Mechanics↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Direct Numerical Simulation of Flow Over a Wall-Mounted Cube with the Nek5000 Spectral Element Code: DNS at Re h = 3900

The direct numerical simulation (DNS) of the canonical wall-mounted unit cube subjected to two distinct incident velocity profiles is performed at a Reynolds number where the bulk flow characteristics are known to become relatively Reynolds number insensitive. The aim of this work is to highlight the sensitivity of such bluff-body flows to mean shear and to provide a representative set of data for such flow scenarios where common turbulence modeling techniques often fail. Simple boundary conditions and a small domain are selected to allow for cost-effective and easy comparisons for model development purposes. In addition to mean velocity profiles, select turbulence statistics are presented in detail. Further, the basic efficacy of eddy viscosity models is evaluated and found to be adequate only for the shear stress components for the bluff-body flows examined here. Other failure mechanisms for RANS are proposed.

42 ENGINEERING↗