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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Asymptotic coefficients of the attached-eddy model derived from an adiabatic atmosphere

The attached-eddy model (AEM) predicts that the mean streamwise velocity and streamwise velocity variance profiles follow a logarithmic shape, while the vertical velocity variance remains invariant with height in the overlap region of high Reynolds number wall-bounded turbulent flows. Moreover, the AEM coefficients are presumed to attain asymptotically constant values at very high Reynolds numbers. Here, the AEM predictions are examined using sonic anemometer measurements in the near-neutral atmospheric surface layer, with a focus on the logarithmic behaviour of the streamwise velocity variance. Utilizing an extensive 210-day dataset collected from a 62 m meteorological tower located in the Eastern Snake River Plain, Idaho, USA, the inertial sublayer is first identified by analysing the measured momentum flux and mean velocity profiles. The logarithmic behaviour of the streamwise velocity variance and the associated ‘−1’ scaling of the streamwise velocity energy spectra are then investigated. The findings indicate that the Townsend–Perry coefficient (A 1 ) is influenced by mild non-stationarity that manifests itself as a Reynolds number dependence. After excluding non-stationary runs, and requiring the bulk Reynolds number defined using the atmospheric boundary layer height to be larger than 4 × 10 7 , the inferred A 1 converges to values ranging between 1 and 1.25, consistent with laboratory experiments. Furthermore, nine benchmark cases selected through a restrictive quality control reveal a close relation between the ‘−1’ scaling in the streamwise velocity energy spectrum and the logarithmic behaviour of streamwise velocity variance. However, additional data are required to determine whether the plateau value of the pre-multiplied streamwise velocity energy spectrum is identical to A 1 .

atmospheric flows↗
DOE OSTI2026

Modulation induced by very-large-scale motions on the inclination angle of wall-attached eddies: an atmospheric surface layer experiment

The forward leaning inclination angle, 𝛾, of coherent turbulent structures is a well-known feature of wall-bounded turbulent flows. Although invariant across friction Reynolds numbers within the range 𝑅𝑒 𝜏 =10 3 −10 6 , 𝛾 can vary significantly across turbulent scales within a high-Reynolds-number flow. Very-large-scale motions (VLSMs) are known to induce significant changes in the instantaneous shear profile, which is a conditioning event that could trigger variability in the inclination angle of smaller coherent turbulent structures. Although this aspect has been extensively studied via numerical and laboratory experiments, few studies have explored this feature for a very-high-Reynolds-number atmospheric flow. In this work, the inclination angle of turbulent structures within the atmospheric surface layer at a very high Reynolds number (𝑅𝑒 𝜏 =7.9 ×10 5 ) is investigated by deploying a scanning Doppler light detection and ranging and a super large particle image velocimetry (SLPIV) apparatus. The inclination angle of wall-attached eddies is inferred either from the two-point correlation of streamwise velocity (𝛾 =41.1∘) or with a scale-dependent approach through the spectral linear stochastic estimator (SLSE). The SLSE (and, thus, the scale-dependent inclination angle) is conditionally evaluated based on the high- and low-momentum events induced by VLSMs, both in the streamwise (𝑢′𝑉𝐿𝑆𝑀) and in the vertical (𝑤′ 𝑉𝐿𝑆𝑀) velocity components. As a result, lower inclination angles (𝛾 =30° −50°) are found for 𝑢$^{'}_{VLSM}$ >0 (𝑤$^{'}_{VLSM}$ <0), while higher values (50° −85°) are ascribed to 𝑢$^{'}_{VLSM}$ <0 (𝑤$^{'}_{VLSM}$ >0). This result emphasises the primary role that VLSMs play in shaping the wall-attached eddy geometry, which, in turn, is crucial to determine the Reynolds stress balance within the wall-attached eddy range.

42 ENGINEERING↗
DOE OSTI2021

Validation of the HyRAM+ physics models for use with propane

The Hydrogen Risk Assessment Model Plus (HyRAM+) toolkit combines quantitative risk assessment with simulations of unignited dispersion, ignited turbulent diffusion flames, and indoor accumulation with delayed ignition of fuels. HyRAM+ is differentiated from HyRAM in that it includes models and leak data for other alternate fuels. The models of the physical phenomena need to be validated for each of the fuels in the toolkit. This report shows the validation for propane which is being used as a surrogate for autogas, which is a mixture of propane and butane and used in internal combustion engines in vehicles. For flame length comparisons, five previously published experiments from peer reviewed journals were used to validate our models. The validation looked at flame lengths and flame widths with respect to different leak diameters, mass flow rates, and source pressures. Most of the sources included more than one set of experimental data, which were collected using different methods (CCD cameras, IR visualization etc.). In general, HyRAM+ overpredicts the flame lengths by around 65%. For heat and radiation models, we compared the heat flux and radiation data reported from two different sources to the values calculated by HyRAM+. For higher mass flow rates, the HyRAM+ calculated flame length results gave a better estimate of what is found in the experiments (65% error), but a higher error (85%) is observed between the HyRAM+ calculated lengths and the experimental flame lengthsfor lower mass flows. Some differences can be attributed to outdoor environmental effects (i.e. wind speed) and uncertainties in jet flame shapes. The propane flame trajectory is predicted for a high Reynolds number case with Re = 12,500 and a low Reynolds number case where Re = 2,000. The Re=12,500 case which is momentum dominated matches well with the experimental flame trajectory, but the agreement for the bouancy driven low Reynolds number case is not as good. Dispersion modeling for unignited propane was also analyzed. We compared the mole fraction, mixture fraction, mean velocity, concentration half width, and inverse mass concentration over an axial distance from different credible journals to the values calculated by HyRAM+. The results display good agreement but generally, HyRAM+ predicts a wider profile for mole fraction and mixture fraction experiments. Overall, HyRAM+’s results are reasonable for predicting the flame length, heat flux, flame trajectory, and dispersion for propane and can be used in risk analyses

08 HYDROGEN↗
DOE OSTI2025

Assessment and Improvement of the SST-Gamma Transition Model in Nalu-Wind

We conduct laminar–turbulent boundary-layer transition simulations using a local correlation-based transition model for two-dimensional incompressible flow and present enhancements to improve the accuracy of transition predictions. Menter’s Galilean-invariant 𝛾 transition model is implemented in the incompressible, unstructured-grid flow solver Nalu-Wind and is validated against experimental data and results from NASA’s flow solvers. The test cases of the AIAA Transition Prediction and Modeling Workshop are investigated, namely, the T3A/T3B flat plates and the NLF(1)-0416 and S809 airfoils. Based on the results, best practices for transition simulations, particularly for an unstructured-grid flow solver, are identified. Additional airfoil simulations are conducted for two wind turbine airfoils, S822 at Reynolds numbers of 𝑂⁡(10 5 ) and DU00-W-212 at Reynolds numbers of 𝑂⁡(10 7 ), to assess the model at low and high Reynolds numbers. Furthermore, through this work, we propose several approaches to enhance transition simulations, including 1) enforcing positivity of the implicit operator for the source terms of the transition model, 2) employing a constant turbulence intensity in stationary external flow simulations, and 3) recommending meshing for unstructured-grid flow solvers. Finally, we provide detailed documentation of the validation and data for the canonical cases to the transition modeling community.

17 WIND ENERGY↗
DOE OSTI2024

LES simulations of a vacuum membrane distillation channel with geometric alterations

3D LES simulations were carried out to study the performance of a vacuum membrane distillation module. A wiggly wall profile with/without embedded stiffeners was considered to alleviate axial and radial temperature polarization, the cause of performance decrease in membrane distillation. Results of the flow field show the wiggles create unsteady vortex shedding inducing intense mixing in the channel. Vortex shedding intensity increases as the Reynolds number increases or stiffeners are added. 3D results show a 56% improvement in flux when moving from a flat sheet membrane to a wiggly membrane with stiffeners at a constant Reynolds number, corresponding to an improvement from 11.7 to 41.7 in the Nusselt number, showing that the alterations improved the flux performance by enhancing the heat transfer along the membrane surface and therefore alleviating temperature polarization, a critical bottleneck in membrane distillation systems. Further, a merit criterion was defined based on the Nusselt number and friction factor, and a 40% increase in merit was shown switching from a flat to a wiggly channel, while a 97% merit increase was seen going from flat to wiggly with stiffeners. A 42% enhancement in the flux was also seen moving from a straight channel to a wiggly channel at a higher Reynolds number which highlights the importance of strategically choosing the mass flow rate, as well as inducing flow separation and vortex shedding in the channel to promote mixing and dissipate the thermal boundary layer. The variation in results between the wiggly module with stiffeners and the wiggly module at high Reynolds numbers suggests enhancing mixing structures can be more impactful on flux for this geometry than increasing the flow rate to a turbulent/transitional regime. However, both are preferable for peak system flux performance. Furthermore, a 2D approximation was used to perform simulations on more extended channels to examine the length degradation. The modules with wiggly channels performed at the same flux level with a doubling of the length. On the other hand, the flat sheet modules experienced length degradation by temperature polarization and dropped in flux yield by around 13% of the short-channel value. This work illustrates that modeling the system and understanding how the performance decreases as the membrane surface area increases are critical for a larger module (scaling up from a lab to a prototype module) to maintain high flux performance.

42 ENGINEERING↗
DOE OSTI2020

Heat Transfer Coefficients of Additively Manufactured Tubes with Internal Pin Fins for Supercritical Carbon Dioxide Cycle Recuperators

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO 2 flowing in additively manufactured tubes with internal pin fins at the US DOE’s National Energy Technology Laboratory in Morgantown, WV. The measurement procedures 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 Reynolds numbers (ReD) ranging from 5 × 10 4 to 2.5 × 10 5 , measured Nu was within 5% of the Dittus-Boelter correlation and measured friction factors were within 5% of the McAdams correlation 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 and were printed at a 30⁰ angle relative to the inner wall (to meet AM constraints). The pin arrangement was helical to promote enhanced heat transfer due to swirl flow. Pin length to diameter aspect ratio was 1.33, 2, and 8, while the pin diameter to tube diameter ratio was 0.188, 0.125, and 0.063. Tests were performed for ReD varying from 6.9 × 10 4 to 2.2 × 10 5 and at conditions equivalent to the low pressure (LP) outlet (8.69 MPa, 361 K) and the high pressure (HP) inlet (20.7 MPa, 350 K) of the low temperature recuperator (LTR) in an indirect sCO 2 cycle. The Wilson plot technique was utilized to measure the bulk heat transfer coefficients. For the best performing design (tube A, pin length to tube diameter ratio: 1.33, pin diameter to tube diameter ratio: 0.19), the local heat transfer coefficient increased by 136% relative to the Dittus-Boelter correlation at the LTR low pressure outlet and 194% at the LTR high pressure inlet. These correspond to a 282% and a 271% increase in the product of the heat transfer coefficient and surface area (adjusted for fin efficiency) product, respectively. Large pressure drops across the test articles were observed. For Tube Design A, the average friction factor, across the range of ReD considered, was significantly larger than the McAdams correlation at both the LTR LP outlet and the LTR HP inlet. A thermal performance factor was utilized to express the ratio of material required to build a finned heat exchanger relative to a finless heat exchanger with the same heat duty and pumping power. Tube Design A was estimated to decrease the required heat exchanger material by 13%.

Searle, Matthew↗
DOE OSTI2022

Large eddy simulation of flow through an axisymmetric sudden expansion

Here, this study aims at investigating the inlet flow conditions of flow through an axisymmetric sudden expansion with an expansion ratio of 2.0. A series of large eddy simulations with the WALE model were conducted for different inlet Reynolds numbers ( Re) and turbulence intensities ($u_{rms}/\bar{U}_m$). The reattachment length, defined as the length measured downstream of the expansion where the flow direction is reversed adjacent to the wall ( Lr), was measured for each case. For widely studied inlet turbulence intensity values (TI), the simulation results are in good agreement with the experimental and numerical results reported in the literature. Parametric studies revealed that turbulence intensity affects the critical Reynolds number, marking the transition between the laminar and transition regions and the reattachment length. The critical Reynolds number was found to decrease with increasing turbulence intensity. A correlation expression is proposed. Additional analysis with proper orthogonal decomposition was performed to enhance the understanding of complex flow structures downstream of the expansion. Finally, an overall correlation expression for the reattachment length was obtained for 500 ≤ Re ≤15 000 and 0.2 ≤ TI (%) ≤ 20. For a given turbulence intensity, the reattachment length can be expressed for laminar and turbulent regions as a function of the Reynolds number. The reattachment length in the transition region can be expressed as a fractional average of reattachment lengths for laminar and turbulent flows.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

On the use of a multigrid-reduction-in-time algorithm for multiscale convergence of turbulence simulations

Simulations of turbulent flow present challenges in terms of accuracy and affordability on modern highly-parallel computer architectures. A multigrid-reduction-in-time algorithm is used to provide a framework for separately evolving different scales of turbulence and for parallelizing the temporal domain, thereby increasing the concurrency. It is hypothesized that the space–time locality of the small scales of turbulence can be used to circumvent difficulties in applying temporal multigrid to flows dominated by inertial physics. For algorithms that fall well short of spectral accuracy (fourth-order is used in this work) attention must be paid to the accuracy of features on scales transferred between multigrid levels. Numerical experiments were performed using implicit large-eddy simulation. Results from applying the approach to an infinite-Reynolds number Taylor–Green flow and a double-shear flow at a Reynolds number of 11650 provide strong evidence that the approach has merit. The multigrid-reduction-in-time framework can be used to parallelize the temporal domain of a high-Reynolds-number turbulent flow and permit independent convergence of different scales. Establishing this foundation allows for future research in reducing the wall-clock time to solve turbulent flows while retaining the same accuracy as sequential solvers. In conclusion, current performance results from parallelizing the temporal domain are not competitive with those from sequential-in-time methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Grid resolution requirement for resolving rare and high intensity wall-shear stress events in direct numerical simulations

Turbulent signals are intermittent with large instantaneous fluctuations. Such large fluctuations lead to small Kolmogorov scales that are hard to resolve in numerical simulations [P. K. Yeung, K. R. Sreenivasan, and S. B. Pope, Effects of finite spatial and temporal resolution in direct numerical simulations of incompressible isotropic turbulence, Phys. Rev. Fluids 3, 064603 (2018)]. Here, the present paper follows the above basic logic, but instead of dissipation events in isotropic turbulence, we study wall-shear stress events in plane channel flow. Wall-shear stress fluctuations are increasingly more intermittent as the Reynolds number increases. Hence, one has to employ higher grid resolutions as the Reynolds number increases in order to resolve a given percentage of wall-shear stress events. The objective of this paper is to quantify effects of the grid resolutions on the rare and high intensity wall-shear stress events. We find that the standard grid resolution resolves about 99% of the wall-shear stress events at Re τ = 180. A slightly higher grid resolution has to be employed in order to resolve 99% of the wall-shear stress events at higher Reynolds numbers, and if the standard grid resolution is used for, e.g., a Re τ = 10000 channel flow, one resolves about 90%–95% wall-shear stress events.

42 ENGINEERING↗
DOE OSTI2020

Data-driven variational multiscale reduced order models

We propose a new data-driven reduced order model (ROM) framework that centers around the hierarchical structure of the variational multiscale (VMS) methodology and utilizes data to increase the ROM accuracy at a modest computational cost. The VMS methodology is a natural fit for the hierarchical structure of the ROM basis: In the first step, we use the ROM projection to separate the scales into three categories: (i) resolved large scales, (ii) resolved small scales, and (iii) unresolved scales. In the second step, we explicitly identify the VMS–ROM closure terms, i.e., the terms representing the interactions among the three types of scales. In the third step, we use available data to model the VMS–ROM closure terms. Thus, instead of phenomenological models used in VMS for standard numerical discretizations (e.g., eddy viscosity models), we utilize available data to construct new structural VMS–ROM closure models. Specifically, we build ROM operators (vectors, matrices, and tensors) that are closest to the true ROM closure terms evaluated with the available data. We test the new data-driven VMS–ROM in the numerical simulation of four test cases: (i) the 1D Burgers equation with viscosity coefficient v = 10 -3 ; (ii) a 2D flow past a circular cylinder at Reynolds numbers, Re = 100, Re = 500, and Re= 1000; (iii) the quasi-geostrophic equations at Reynolds number Re = 450 and Rossby number Ro = 0.0036; and (iv) a 2D flow over a backward facing step at Reynolds number Re = 1000. The numerical results show that the data-driven VMS–ROM is significantly more accurate than standard ROMs.

42 ENGINEERING↗
DOE OSTI2021

Nondimensional convection numbers modeling thermally stratified storage tanks: Richardson's number and hot-water tanks

Thermally stratified storage tank studies have spanned over 50 years to increase the thermal storage efficiency and accurate prediction of the outlet temperature particularly for solar applications. The studies have reviewed and modeled the jet and plume flow phenomena inside the tank due to the inlet mixing and stratification level. Kelvin–Helmholtz and Rayleigh–Taylor instabilities are the major drivers of the mixing in these tanks. Momentum jets deflecting off walls at the bottom of the tank also create significant mixing. Reviewing Richardson models shows that the categorization was based on the range of Reynolds numbers at the inlet. Unfortunately, the use of superficial velocity in calculating the Richardson number results in critical values in the literature ranging from below 0.25 to 100. The most used length scale associated with these flows is an inertial scale based on the tank height or diameter although the mixing can occur at a relatively smaller scale. The various inlet devices and a large span of flow rates experienced in thermally stratified storage tanks requisite the use of the Reynolds number in combination with a convection number for accurate one dimensional models that predict performance over the long-term. The evaluation of peak shifting of electric loads leveraging renewable sources for applications, including residential hot-water tanks, commercial water tanks, and large-scale chilled water storage tanks, require these models. Finally, this paper is focused on establishing the significance of the convection numbers in conjunction with the Reynolds number for modeling the thermal stratification in storage tanks.

14 SOLAR ENERGY↗
DOE OSTI2024

Direct Numerical Simulation Database of High-Speed Flow over Parameterized Curved Walls

This study presents a direct numerical simulation (DNS) database of high-speed turbulent boundary layers (TBLs) subject to pressure gradients due to parametrically varied backward-facing and forward-facing wall curvatures, with an inflow Mach number of 4.9 and a friction Reynolds number of [Formula: see text] immediately before the onset of wall curvature. The Mach and Reynolds numbers are significantly higher than those reported in the literature for the DNS of pressure-gradient TBLs. The flow conditions and baseline wall geometries are representative of experiments in the high-speed blowdown wind tunnel at the National Aerothermochemistry Laboratory at Texas A&M University. The wall steepness of the baseline geometry for both the backward-facing and forward-facing walls was systematically varied to cause attached, incipiently separated, and fully separated flows. Precomputed flow statistics, including turbulent kinetic energy budgets, are available on the website of the Turbulence Modeling Resource of the NASA Langley Research Center, allowing other investigators to query any property of interest.

Engineering↗
DOE OSTI2023

Taylor–Couette flow for astrophysical purposes

A concise review is given of astrophysically motivated experimental and theoretical research on Taylor–Couette flow. The flows of interest rotate differentially with the inner cylinder faster than the outer, but are linearly stable against Rayleigh’s inviscid centrifugal instability. At shear Reynolds numbers as large as 10 6 , hydrodynamic flows of this type (quasi-Keplerian) appear to be nonlinearly stable: no turbulence is seen that cannot be attributed to interaction with the axial boundaries, rather than the radial shear itself. Direct numerical simulations agree, although they cannot yet reach such high Reynolds numbers. This result indicates that accretion-disc turbulence is not purely hydrodynamic in origin, at least insofar as it is driven by radial shear. Theory, however, predicts linear magnetohydrodynamic (MHD) instabilities in astrophysical discs: in particular, the standard magnetorotational instability (SMRI). MHD Taylor–Couette experiments aimed at SMRI are challenged by the low magnetic Prandtl numbers of liquid metals. High fluid Reynolds numbers and careful control of the axial boundaries are required. The quest for laboratory SMRI has been rewarded with the discovery of some interesting inductionless cousins of SMRI, and with the recently reported success in demonstrating SMRI itself using conducting axial boundaries. Some outstanding questions and near-future prospects are discussed, especially in connection with astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Performance evaluation of finned tube heat exchanger using curved wavy delta winglet vortex generators with circular perforations

Vortex generation is recognized as an effective passive approach to improve the heat transfer rate in fin and tube heat exchangers (FTHEs). The current study proposed innovative designs of curved wavy delta winglet vortex generators (CWDWVGs), both without and with circular perforations, to improve the heat transfer efficiency of FTHEs. There is potential to increase heat transfer performance further through various CWDWVG designs. Here, this study explores seven unique CWDWVG configurations, from 1-wave to 7-wave. A 3-D computational numerical model is utilized to evaluate the Thermo-hydraulic performance of FTHEs fitted with these different CWDWVG configurations across Reynolds numbers from 400 to 2000. This comparative analysis of the Thermo-hydraulic performance of FTHEs featuring four parallel circular tube layouts assesses configurations both with and without vortex generators (VGs) and various hole configurations. The evaluation of Thermo-hydraulic performance involves different parameters, including the London area goodness factor (LAGF), Colburn factor (j), friction factor (f), pressure drop (?P), and Nusselt number (Nu). Results demonstrate that the various CWDWVG configurations and the number of holes in them substantially affect the efficiency, as evaluated by the dimensionless Performance Evaluation Criteria (PEC). Notably, the 7-wave CWDWVGs surpassed other configurations, and integrating circular punched perforations further improved the thermal-hydraulic performance of FTHE. Specifically, the 7-wave CWDWVGs without holes demonstrated superior performance over other configurations, showing a significant increase in Nusselt number by 75.18% and 85.16% at Reynolds numbers of 2000 and 400, respectively, alongside an increase in pressure drop by 216.38% to 224.96%. Meanwhile, the 7-wave CWDWVGs with eight holes, in comparison to those without holes, exhibited a Nusselt number increase of 0.85%, a pressure drop decrease of 7.31%, and a reduction in the friction factor by 5.82%.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Implications of inertial subrange scaling for stably stratified mixing

In this work, we investigate the effects of the turbulent dynamic range on active scalar mixing in stably stratified turbulence by adapting the theoretical passive scalar modelling and demonstrating their usefulness through consideration of the results of direct numerical simulations of statistically stationary homogeneous stratified and sheared turbulence. By analysis of inertial and inertial–convective subrange scalings, we show that the relationship between the active scalar and turbulence time scales is predicted by the ratio of the Kolmogorov and Obukhov–Corrsin constants, provided mean flow parameters permit the two subrange scalings to be appropriate approximations. We use the resulting relationship between time scales to parameterise an appropriate turbulent mixing coefficient $\varGamma \equiv \chi /\epsilon$, defined here as the ratio of available potential energy ($E_p$) and turbulent kinetic energy ($E_k$) dissipation rates. With the analysis presented here, we show that $\varGamma$ can be estimated by $E_p,E_k$ and a universal constant provided an appropriate Reynolds number is sufficiently high. This large Reynolds number regime appears here to occur at $ {{Re_b}} \equiv \epsilon / \nu N^{2} \gtrapprox 300$ where $\nu$ is the kinematic viscosity and $N$ is the characteristic buoyancy frequency. We propose a model framework for irreversible diapycnal mixing with robust theoretical parametrisation and asymptotic behaviour in this high-$ {{Re_b}}$ limit.

42 ENGINEERING↗
DOE OSTI2020

Experimental Study on the Effect of Localized Blockages on the Friction Factor of a 61-Pin Wire-Wrapped Bundle

Abstract The thermal-hydraulic behavior of the flow in rod bundles has motivated numerous experimental and computational investigations. Previous studies have identified potential for accumulation of debris within the small subchannels of typical wire-wrapped assemblies with subsequent total or partial blockage of subchannel coolant flow. A test campaign was conducted to study the effects of localized blockages on the bundle averaged friction factor of a tightly packed wire-wrapped rod bundle. Blockages were installed within the bundle, and fluid pressure drop was measured across one wire pitch for a Reynolds number range of 500–17,200. The Darcy–Weisbach friction factor of the perturbed rod bundle geometry was compared with that of the unblocked bundle, as well as with the predictions of a well-established friction factor correlation. Differing effects based on blockage size and location for various flow regimes were studied. A number of conclusions can be made about the effects of the blockages on the friction factor, such as an increasing effect of the blockage on friction factor with an increase in Reynolds number, a change in flow behavior in the turbulent transition flow regime near Reynolds number 3000, differences in effect on friction factor for different types of subchannel blockage, and a nonlinear trend in friction factor variation with flow area impeded for edge subchannels. To this end, all data and quantified uncertainty produced in this study are made available for comparison and validation of advanced computational tools.

Engineering↗
DOE OSTI2021

On the shear-current effect: toward understanding why theories and simulations have mutually and separately conflicted

ABSTRACT The shear-current effect (SCE) of mean-field dynamo theory refers to the combination of a shear flow and a turbulent coefficient β21 with a favourable negative sign for exponential mean-field growth, rather than positive for diffusion. There have been long-standing disagreements among theoretical calculations and comparisons of theory with numerical experiments as to the sign of kinetic ($\beta ^u_{21}$) and magnetic ($\beta ^b_{21}$) contributions. To resolve these discrepancies, we combine an analytical approach with simulations, and show that unlike $\beta ^b_{21}$, the kinetic SCE $\beta ^u_{21}$ has a strong dependence on the kinetic energy spectral index and can transit from positive to negative values at $\mathcal {O}(10)$ Reynolds numbers if the spectrum is not too steep. Conversely, $\beta ^b_{21}$ is always negative regardless of the spectral index and Reynolds numbers. For very steep energy spectra, the positive $\beta ^u_{21}$ can dominate even at energy equipartition urms ≃ brms, resulting in a positive total β21 even though $\beta ^b_{21}\lt 0$. Our findings bridge the gap between the seemingly contradictory results from the second-order-correlation approximation versus the spectral-τ closure, for which opposite signs for $\beta ^u_{21}$ have been reported, with the same sign for $\beta ^b_{21}\lt 0$. The results also offer an explanation for the simulations that find $\beta ^u_{21}\gt 0$ and an inconclusive overall sign of β21 for $\mathcal {O}(10)$ Reynolds numbers. The transient behaviour of $\beta ^u_{21}$ is demonstrated using the kinematic test-field method. We compute dynamo growth rates for cases with or without rotation, and discuss opportunities for further work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

An LES study of secondary motion and wall shear stresses in a pipe bend

In this study, large-eddy simulations (LES) of a single-phase, turbulent flow in a 90 ° pipe bend are performed at three Reynolds numbers (5300, 27 000, and 45 000) to investigate the correlation between secondary flow motion and wall shear stresses, which is suspected to be a potential mechanism responsible for material erosion. The isothermal flows are validated against available experimental and numerical data first. The snapshot proper orthogonal decomposition (POD) is applied for the medium and high Reynolds number flows to identify the secondary flow motions and the oscillation of the Dean vortices that are found to cause swirl-switching. Distinguished frequencies of the POD time coefficients at Strouhal numbers of 0.25 and 0.28 are identified for Reynolds numbers at 27 000 and 45 000, respectively. Moreover, shear stress on the pipe wall and the associated power spectral density are obtained and shown to have the same oscillating frequency as the swirl-switching.

42 ENGINEERING↗