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.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Quantifying Uncertainty in Ice Particle Velocity–Dimension Relationships Using MC3E Observations

Abstract Ice particle terminal fall velocity ( V t ) is fundamental for determining microphysical processes, yet remains extremely challenging to measure. Current theoretical best estimates of V t are functions of Reynolds number. The Reynolds number is related to the Best number, which is a function of ice particle mass, area ratio ( A r ), and maximum dimension ( D max ). These estimates are not conducive for use in most models since model parameterizations often take the form , where ( α , β ) depend on habit and D max . A previously developed framework is used to determine surfaces of equally plausible ( α , β ) coefficients whereby ice particle size/shape distributions are combined with V t best estimates to determine mass- ( V M ) or reflectivity-weighted ( V Z ) velocities that closely match parameterized V M ,SD or V Z ,SD calculated using the ( α , β ) coefficients using two approaches. The first uses surfaces of equally plausible ( a , b ) coefficients describing mass ( M )–dimension relationships (i.e., ) to calculate mass- or reflectivity-weighted velocity from size/shape distributions that are then used to determine ( α , β ) coefficients. The second investigates how uncertainties in A r , D max , and size distribution N ( D ) affect V M or V Z . For seven of nine flight legs flown on 20 and 23 May 2011 during the Mesoscale Continental Convective Clouds Experiment (MC3E), uncertainty from natural parameter variability—namely, the variability in ice particle parameters in similar meteorological conditions—exceeds uncertainties arising from different A r assumptions or D max estimates. The combined uncertainty between A r , D max , and N ( D ) produced smaller variability in ( α , β ) compared to varying M ( D ), demonstrating M ( D ) must be accurately quantified for model fall velocities. Primary sources of uncertainty vary considerably depending on environmental conditions. Significance Statement Ice particle fall velocity is fundamental for numerous processes within clouds, and hence is a critical property that must be accurately represented in weather and climate models. Using aircraft observations of ice particle shapes and sizes obtained in clouds behind midlatitude thunderstorms, this work develops a new framework for estimating ice particle fall velocities and their uncertainty, including quantifying the importance of different uncertainty sources from cloud microphysics measurements. Natural parameter variability contributes the most uncertainty in ice particle fall velocity estimates, although other sources can also be important contributors to uncertainty in certain conditions. Additional work examining ice particle data is needed to further understand how dependent uncertainty in certain ice particle properties are to local environmental conditions.

Meteorology & Atmospheric Sciences↗

DOE OSTI2023

How do magnetized HED flows transition to collimated plasma jets? (Final Report)

We present results from a pulsed-power experiment designed to generate accretion-like flows in the laboratory that may be used to study the dynamics of astrophysical accretion-driven outflows. To emulate systems like Young Stellar Object (YSO) outflows, we create radially converging plasma flows in a short disk that transition into much longer bipolar outflows. Our experiment uses a ~1 MA pulsed-power driver to generate plasma from an array of thin aluminum wires mounted to a 3D-printed current path that controls the angular momentum and magnetic field topology. In this paper, we show that experimental data scale favorably to observed properties of YSO jets. This set of experiments tests the simplest version of our possible plasma dynamics, where outflows are generated without rotation or axial magnetic field present. The temperature, velocity and density of the disk and outflows are characterized using interferometry, gated optical and ultraviolet imaging, and Thomson scattering diagnostics. We discuss the Reynolds number, magnetic Reynolds number and Mach numbers calculated in three different stages of the experiment, and compare our results to 3D numerical simulations in the PERSEUS magnetohydrodynamics code. This report is a summary of the work done by the graduate students who were supported by the project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

A Study of the Transition to Turbulence in a Bed of 67 Spherical Pebbles

Packed beds are commonly found in many engineering systems and have been widely studied for decades. A relatively new packed bed system is the Pebble Bed Reactor, a type of generation-IV nuclear reactor. Unlike many of the packed beds encountered in chemical and process engineering applications, Pebble Bed Reactors are larger and operate at significantly higher Reynolds numbers. As a result of these differences, there is a very limited amount of information on the detailed flow physics that exist in these complex geometries. This work seeks to contribute to a growing database of flow data for Pebble Bed Reactor systems by performing Direct Numerical Simulations of the flow in an experimental bed of 67 pebbles for a range of conditions. Simulations are performed at a Prandtl number of 0.66 and Reynolds numbers from 300–600. These Reynolds numbers are chosen to gain additional knowledge on the spatial development of turbulence in these systems. Analysis of the Turbulent Kinetic Energy, turbulence anisotropy, and Turbulent Heat Flux is performed. Results demonstrate significant development of the TKE across the tested range of Reynolds numbers. Examination of both the TKE and THF reveal that development first occurs near the center of the bed and propagates radially as the flow moves further into the bed. Notable regions of negative production of turbulent kinetic energy are observed in regions where flow accelerates around pebble contact points. Furthermore, these regions are found to coincide with regions of 1-component turbulence.Kindly check and confirm, all authors email id is correctly identified.These are correct

Direct numberical simulation↗

DOE OSTI2022

Investigation of Gas Turbine Internal Cooling Using Supercritical CO 2 —Effect of Surface Roughness and Channel Aspect Ratio

In this paper, an experimental and numerical investigation of internal cooling channels with rib turbulators is presented with sCO 2 as the working fluid at process conditions (pressure-20.7 MPa and temperature up to 150 °C). The effect of channel aspect ratio up to 2:1 on thermal-hydraulic performance is explored in additively manufactured rectangular channels and square channels, both with and without 60 deg ribs on the top and bottom sides. The Wilson-plot method is employed to experimentally measure channel-averaged Nusselt number over a Reynolds number range up to 370,000. The friction factor is calculated from pressure drop and mass flow rate and additionally, the overall thermal performance factor (TPF) is reported. A companion computational fluid dynamics (CFD) simulation is performed for the rib turbulated cooling configurations reported in the experiments using the Reynolds average Navier–Stokes-based turbulence model. The objective of the numerical study is to gain insight into the local heat transfer augmentation in the ribbed channels as a result of varying the aspect ratio, channel configuration (square versus rectangular), operating conditions (Reynolds number) and the surface roughness, an inherent outcome of the additive manufacturing process. Surface roughness is simulated using sand grain roughness height (K S ) calculated from the experimental data, and a comparison is presented with the corresponding channel configuration with varying surface roughness heights starting from smooth surfaces (K S = 0). Experimental results indicate that the heat transfer augmentation is negligible in the rectangular channels with ribs on the long side compared to the square channel. However, it is enhanced by 60% in comparison to placing ribs on the shorter side. The TPF remains constant at around 1 for the entire range of Reynolds numbers consistent with prior work at the National Energy Technology Laboratory (NETL). Here, the simulation results highlight that increased surface roughness can have a favorable considerable influence on Nusselt number and overall thermal performance enhancement.

42 ENGINEERING↗

DOE OSTI2024

Fluid inertia controls mineral precipitation and clogging in pore to network-scale flows

Mineral precipitation caused by fluid mixing presents complex control and predictability challenges in a variety of natural and engineering processes, including carbon mineralization, geothermal energy, and microfluidics. Precipitation dynamics, particularly under the influence of fluid flow, remain poorly understood. Combining microfluidic experiments and three-dimensional reactive transport simulations, we demonstrate that fluid inertia controls mineral precipitation and clogging at flow intersections, even in laminar flows. We observe distinct precipitation regimes as a function of Reynolds number (Re). At low Reynolds numbers (Re < 10), precipitates form a thin, dense layer along the mixing interface, which shuts precipitation off, while at high Reynolds numbers (Re > 50), strong three-dimensional flows significantly enhance precipitation over the entire intersection, resulting in rapid clogging. When injection rates from two inlets are uneven, flow symmetry-breaking leads to unexpected flow bifurcation phenomena, which result in enhanced concurrent precipitation in both downstream channels. Finally, we extend our findings to rough channel networks and demonstrate that the identified inertial effects on precipitation at the intersection scale are also present and even more dramatic at the network scale. This study sheds light on the fundamental mechanisms underlying mixing-induced mineral precipitation and provides a framework for designing and optimizing processes involving mineral precipitation.

Science & Technology - Other Topics↗

DOE OSTI2022

Conditional analysis of temperature and strain rate effects on dissipation structure in turbulent non-premixed jet flames

Here, this work presents results from simultaneous high-resolution temperature and velocity measurements in a series of turbulent non-premixed jet flames. The filtered Rayleigh scattering (FRS)-based temperature measurements demonstrate sufficient signal-to-noise (SNR) and spatial resolution to estimate the smallest scalar length scales and accurately determine dissipation rate fields. A comprehensive set of conditional statistics are used to characterize the small-scale structure, including the dependence of dissipation layer widths on Reynolds number, temperature, and dissipation magnitude. In general, the dissipation layer thickness decrease with increasing Reynolds number and increase with increasing temperature. However, dissipation layer widths show two distinct behaviors with respect to dissipation magnitude. For small dissipation values, increases in magnitude results in broadening of the dissipation layer, while for larger magnitude values of dissipation, the layer widths are thinned, highlighting the complexity of small-scale turbulent mixing. Additionally, measured ratios of the dissipation layer width to the Batchelor length scale are consistent across all Reynolds numbers and agree with previous studies in non-reacting flows. The unique aspect about the current set of measurements is the ability to examine the interaction of dissipation structure with turbulent flow parameters for the first time in turbulent non-premixed flames. Particularly, the strain rate/dissipation relationship is examined and compared to previous studies in non-reacting flows. It is found that the dissipation layers tend to align normal to the principal compressive strain axis and this tendency increases with increasing Reynolds number. For the lowest Reynolds number case, no dependence of the dissipation layer width nor dissipation rate magnitude on strain rate is found. However, for higher Reynolds numbers, a strong dependence of the dissipation layer width and dissipation rate magnitude on the principal compressive strain rate is observed. These results indicate the direct role of the compressive strain rate field on small-scale mixing structure in reacting flows.

42 ENGINEERING↗

DOE OSTI2020

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Small-scale properties from exascale computations of turbulence on a $\mathbf{32\,768^3}$ periodic cube

To study the physics of small-scale properties of homogeneous isotropic turbulence at increasingly high Reynolds numbers, direct numerical simulation results have been obtained for forced isotropic turbulence at Taylor-scale Reynolds number R λ = 2500 on a 32 768 3 three-dimensional periodic domain using a GPU pseudo-spectral code on a 1.1 exaflop GPU supercomputer (Frontier). These simulations employ the multi-resolution independent simulation (MRIS) technique (Yeung & Ravikumar 2020, Phys. Rev. Fluids, vol. 5, 110517) where ensemble averaging is performed over multiple short segments initiated from velocity fields at modest resolution, and subsequently taken to higher resolution in both space and time. Reynolds numbers are increased by reducing the viscosity with the large-scale forcing parameters unchanged. Although MRIS segments at the highest resolution for each Reynolds number last for only a few Kolmogorov time scales, small-scale physics in the dissipation range is well captured – for instance, in the probability density functions and higher moments of the dissipation rate and enstrophy density, which appear to show monotonic trends persisting well beyond the Reynolds number range in prior works in the literature. Attainment of range of length and time scales consistent with classical scaling also reinforces the potential utility of the present high-resolution data for studies of short-time-scale turbulence physics at high Reynolds numbers where full-length simulations spanning many large-eddy time scales are still not accessible. A single snapshot of the 32 768 3 data is publicly available for further analyses via the Johns Hopkins Turbulence Database.

intermittency↗

DOE OSTI2024

On the stability of a pair of vortex rings

The growth of perturbations subject to the Crow instability along two vortex rings of equal and opposite circulation undergoing a head-on collision is examined. Unlike the planar case for semi-infinite line vortices, the zero-order geometry of the flow (i.e. the ring radius, core thickness and separation distance) and by extension the growth rates of perturbations vary in time. The governing equations are therefore temporally integrated to characterize the perturbation spectrum. The analysis, which considers the effects of ring curvature and the distribution of vorticity within the vortex cores, explains several key flow features observed in experiments. First, the zero-order motion of the rings is accurately reproduced. Next, the predicted emergent wavenumber, which sets the number of secondary vortex structures emerging after the cores come into contact, agrees with experiments, including the observed increase in the number of secondary structures with increasing Reynolds number. Finally, the analysis predicts an abrupt transition at a critical Reynolds number to a regime dominated by a higher-frequency, faster-growing instability mode that may be consistent with the experimentally observed rapid generation of a turbulent puff following the collision of rings at high Reynolds numbers.

42 ENGINEERING↗

DOE OSTI2024

Heat Transfer Experiments of Ribbed, Serpentine Cooling Passages with Supercritical CO2

The results of heat transfer experiments with supercritical CO2 (sCO2) to capture the performance of ribbed, serpentine cooling passages for Reynolds numbers in the range of 100,000 - 400,000 are presented. The oxy-combustion turbine operating in an Allam-Fetvedt cycle has a turbine inlet temperature exceeding 1100°C, necessitating internally cooled turbine blades due to material creep strength limits. Reynolds numbers for sCO2, with its combination of high fluid density and low viscosity, can be up to an order of magnitude higher than comparable internal cooling conditions for an air-breathing engine where the majority of experimental data exists. To improve the prediction of cooling performance through turbulated passages, a test rig was designed and constructed with interchangeable insert capability for operation in sCO2 with heat transfer from a heater outlet flow stream to a lower temperature recuperator outlet flow stream. The inserts included serpentine passage geometry with five passes, one of the inserts including chevron ribs and the other a plain wall. The multi-pass serpentine geometry allowed for the evaluation of the effects 180 deg. tip turns have on the return passage flow characteristics, as they can depart significantly from fully-developed flow characteristics. The test procedure included the use of a modified Wilson plot method, with hot flow stream conditions kept effectively constant so that monitored changes in the overall thermal resistance could be attributed to manipulations in cooling flow stream conditions. Conditions during testing consisted of test section pressures of approximately 180 bar, cooling flow stream temperature of 175°C and hot flow stream temperature of 410°C. Test points at discrete cool flow Reynolds numbers up to 400,000 demonstrated an expected trend of decreased fluid thermal resistance with increasing Reynolds number, while at comparable Reynolds number there was an exhibited lower thermal resistance and higher passage pressure drop for the ribbed insert compared to the plain-walled insert. Processing of test data included establishment of Nusselt number enhancement ratios, so that for different internal cooling passage dimensions in the actual turbine blade at comparable Reynolds number the ratio could be applied to predict cooling performance.

Marshall, Michael↗

DOE OSTI2024

Generating Buoyancy-Driven Convection in Membrane Distillation

Membrane distillation (MD) is a thermally-driven desalination process that can treat hypersaline brines. Considerable MD literature has focused on mitigating temperature and concentration polarization. This literature largely neglects that temperature and concentration polarization increase the feed density near the membrane. With gravity properly oriented, this increase in density could trigger buoyancy-driven convection and increase permeate production. Convection could also be strengthened by heating the feed channel wall opposite the membrane. To investigate that possibility, we perform a series of experiments using a plate-and-frame direct contact MD system with an active membrane area of 300 cm2 and a feed channel wall heated using a resistive heater. The experiments measure the average transmembrane permeate flux for two gravitational orientations, feed Reynolds numbers between 128 and 1128, and wall heat fluxes up to 12 kW/m2. The results confirm that with gravity properly oriented, wall-heating can trigger buoyancy-driven convection for a wide range of feed Reynolds numbers, and increase permeate production between roughly 20 and 130 %. We estimate, however, that at high Reynolds numbers (Re > 800), more than 70 % of the wall heat is carried out of the MD system by the feed flow, without contributing to permeate production. This suggests the need for longer membranes and heat recovery steps in any future practical implementation.

buoyancy-driven convection↗

DOE OSTI2021

The discrete Green's function paradigm for two-way coupled Euler–Lagrange simulation

We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within 6% at all wall-normal separations for particle Reynolds numbers up to 24. Furthermore, the discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.

42 ENGINEERING↗

DOE OSTI2020

The discrete Green's function paradigm for two-way coupled Euler-Lagrange simulation

We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.

42 ENGINEERING↗

DOE OSTI2021

A review of modelling techniques for floating offshore wind turbines

Modelling floating offshore wind turbines (FOWTs) is challenging due to the strong coupling between the aerodynamics of the turbine and the hydrodynamics of the floating platform. Physical testing at scale is faced with the additional challenge of the scaling mismatch between Froude number and Reynolds number due to working in the two fluid domains, air and water. In the drive for cost-reduction of floating wind energy, designers may be seeking to move towards high-fidelity numerical modelling as a substitute for physical testing. However, the numerical engineering tools typically used for FOWT modelling are considered as mid-fidelity to low-fidelity tools, and currently lack the level of accuracy required to do so. Furthermore, there is a lack of operational FOWT data available for further development and validation. High-fidelity tools, such as CFD, have greater accuracy but are cumbersome tools and still require validation. Physical scale model testing therefore continues to play an essential role in the development of FOWTs both as a source of validation data for numerical models and as an important development step along the path to commercialization of all platform concepts. The aim of this paper is to provide an overview of both numerical modelling and physical FOWT scale model testing approaches and to provide guidance on the selection of the most appropriate approach (or combination of approaches). The current state-of-the-art will be discussed along with current research trends and areas for further investigation.

17 WIND ENERGY↗

DOE OSTI2021

CFD Evaluation of Pressure Change Along Coolant Passages in Sodium-Cooled Fast Reactor with Nek5000

To support the design efforts of advanced sodium-cooled fast reactors (SFRs), a series of computational fluid dynamics (CFD) simulations are performed to investigate the pressure change along various flow passages in the proposed SFR system. The simulations are carried out with the state-of-the-art spectral element flow solver, Nek5000. Two specific case studies are presented in this paper: the flow exiting the axial neutron reflector channels and the flow entering the fuel pin bundle. Due to the high Reynolds numbers expected, a Reynolds-averaged Navier-Stokes (RANS) approach is necessary to model the turbulence. A newly developed regularized k-omega RANS model is adopted in the related CFD calculations. The first case study explores the effect of Reynolds number on the pressure change when flow exits the reflector channels. The pressure change in this case has two major contributors: the change due to wall friction and the Bernoulli effect. It is noted that the nondimensional pressure loss follows a log-linear trend up to Re = 10 5 , and then the trend is flattened. In the second case study, the advanced NekNek coupling capability is tested where an integral domain can be divided into multiple subdomains with coupling interfaces, which would greatly ease the meshing process of complex engineering geometries and potentially save computational resources. The preliminary results obtained so far confirm the consistency between the NekNek results and those produced by regular Nek5000 simulation. The presented work demonstrates the readiness and flexibility of the related CFD techniques, which is part of the broader effort to leverage cutting-edge CFD to inform the advanced nuclear reactor designs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Experimental Study of Eutectic Molten Salts NaCl/KCl/ZnCl2 Heat Transfer Inside a Smooth Tube for High-Temperature Application

Abstract Eutectic salts NaCl-KCl-ZnCl2 and NaCl-KCl-MgCl2 are two of the chloride salt systems that are promising for being used as high-temperature heat transfer fluid (HTF) and thermal energy storage (TES) materials in a wide range of temperatures from 230 °C to 750 °C in concentrated solar thermal power systems. To conduct measurement of the heat transfer coefficient of the molten salt at high temperatures, a circulation system and instrumentation of flow and heat transfer was constructed. Experimental measurement of the convective heat transfer coefficients of NaCl-KCl-ZnCl2 (molar fraction: 13.8–41.9–44.3%) inside tubes has been accomplished to find the most suitable heat transfer correlations. Experience of salt processing and operation of the high-temperature heat transfer test system was obtained. Two correlations, Dittus-Boelter equation and Gnielinski’s correlation for Nusselt number against Reynolds number and Prandtl number, are evaluated using the test results, and the latter correlation is recommended due to its better agreement of prediction against tested data.

Energy & Fuels↗

DOE OSTI2023

Large Eddy Simulation of Externally Induced Ingress about an Axial Seal by Stator Vanes

Turbine inlet temperatures in advanced gas turbines could be as high as 2000 °C. To prevent ingress of this hot gas into the wheelspace between the stator and rotor disks, whose metals can only handle temperatures up to 850 °C, rim seals and sealing flows are used. This study examines the abilities of large eddy simulation (LES) based on the WALE subgrid model and Reynolds-averaged Navier–Stokes (RANS) based on the SST model in predicting ingress in a rotor–stator configuration with vanes but no blades, a configuration with experimental data for validation. Results were obtained for an operating condition, where the ratio of the external Reynolds number to the rotational Reynolds number is 0.538. At this operating condition, both LES and RANS were found to correctly predict the coefficient of pressure, C p , located downstream of the vanes and upstream of the seal, but only LES was able to correctly predict the sealing effectiveness. This shows C p by itself is inadequate in quantifying externally induced ingress. RANS was unable to predict the sealing effectiveness because it significantly under predicted the pressure drop in the hot gas path along the axial direction, especially about the seal region. This affected the pressure difference across the seal in the radial direction, which ultimately drives ingress.

42 ENGINEERING↗

DOE OSTI2021

Simultaneous velocity and density measurements of fully developed Rayleigh-Taylor mixing

Here, the dynamics of molecular mixing and the energy transfer process in the Rayleigh-Taylor instability (RTI) are studied through the collection of simultaneous density-velocity measurements. These experiments provide simultaneous density-velocity field measurements, in contrast to previous point measurements. Statistically stationary experiments are performed in a “convective-type” gas tunnel facility, with density contrast achieved through the injection of helium into the bottom stream. Three experiments at Atwood number ≈ 0.1 are captured at three outer scale Reynolds numbers Re = 520, 2260, and 4050. Particle image velocimetry and laser induced fluorescence are employed simultaneously. Statistics of the density and velocity show self-similar collapse of RTI profiles at large Reynolds number Re >2000. Flat velocity profiles indicate homogeneous turbulence characteristics in the core of the mixing region. Significant anisotropy develops in the flow, with horizontal velocity fluctuations being only 60% of the vertical velocity fluctuations. The turbulent mass flux is found to be asymmetric about the centerline, with increased peak towards the spike. Measurements of the molecular mixing show that mixing is maximized at the core of the flow and increases with increased Reynolds number. The transport equation of density-specific-volume correlation b shows that it is mostly produced in the core of the mixing region, and that the spatial evolution of its profile is the result of transport by bulk motion of the bubble and spike. Energy transfer from gravitational potential energy to turbulent kinetic energy and viscous dissipation is observed to occur in the experiment with a ratio of dissipated energy to potential energy released of 38%. The analysis of the turbulent kinetic energy transport equation budget reveals that production is the primary mechanism towards the growth of turbulent kinetic energy in the core of the flow, and is asymmetrically slightly skewed towards the spike. However, it is through the transport that the strong advection at the edges of the mixing region is maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗