SEARCH · Engineering Papers

Results for “Low Reynolds Number”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗

NASA NTRS1988

Boundary-layer control for drag reduction

Although the number of possible applications of boundary-layer control is large, a discussion is given only of those that have received the most attention recently at NASA Langley Research Center to improve airfoil drag characteristics. This research concerns stabilizing the laminar boundary layer through geometric shaping (natural laminar flow, NLF) and active control involving the removal of a portion of the laminar boundary layer (laminar flow control, LFC) either through discrete slots or a perforated surface. At low Reynolds numbers, a combination of shaping and forced transition has been used to achieve the desired run of laminar flow and control of laminar separation. In the design of both natural laminar flow and laminar flow control airfoils and wings, boundary layer stability codes play an important role. A discussion of some recent stability calculations using both incompressible and compressible codes is given.

Harvey, William D.↗

NASA NTRS1987

Advanced experimental techniques for transonic testing

Problems common to wind tunnels, including low Reynolds numbers, wall interference, support interference, and flow unsteadiness are discussed, and it is noted that with increasing aircraft size, existing transonic tunnels are becoming more inadequate in terms of their Reynolds-number capability. Problems attributed to scale effect are covered, and techniques for increasing Reynolds number are presented. These include operating with a heavy gas or a mixture of gases, increasing the operating pressure in the tunnel, increasing the size of the model, and/or reducing the test temperature. Acceptable levels of flow unsteadiness are reviewed, along with the problems associated with the flow over wind-tunnel models influenced by flow unsteadiness. Sources of flow unsteadiness present in a continuous transonic tunnel driven by a compressor or a fan are listed, and rules for ensuring low levels of flow unsteadiness are given.

Kilgore, Robert A.↗

NASA NTRSDate not supplied

Wall-modeled LES of the Three-dimensional Speed Bump Experiment

We evaluate the performance of wall-modeled large eddy simulation (WMLES) in predicting smooth-body turbulent flow separation for a three-dimensional Gaussian-shaped speed bump geometry. The Reynolds number based on the bump length is 2 million, and an unstructured compressible finite-volume solver is used with an equilibrium wall model and dynamic subgrid scale model. Spanwise periodic simulations of the centerline two-dimensional bump were used to assess grid resolution requirements, and resolving the thin internal layer in the accelerating region was found to be necessary to correctly capture the downstream separated flow region. Based on these insights, an optimized grid that was smaller by a factor of two provided comparable accuracy to a finer grid that has been used by us and other researchers in past studies. Using insights gained from the spanwise periodic simulation, an unstructured polyhderal grid with about 250 million cells was used for the three-dimensional (3D) configuration with inviscid tunnel side and top walls. Detailed comparisons of wall skin-friction coefficient, wall pressure, velocity and turbulent stresses with available experimental data indicated excellent agreement. While the low-Reynolds Number Spalart-Allmaras RANS model with rotation/curvature correction gave qualitatively good agreement with experiments, WMLES showed significantly more accurate quantitative predictions of the flowfield. The 3D WMLES showed good agreement with experiments in the separated region, and the centerplane results indicated a different separation topology compared to the spanwise periodic simulation.

Computational Fluid Dynamics↗

NASA NTRS1999

Boundary Layer Transition in the NTF: HSR Experience and Plans

Efforts towards understanding boundary layer transition characteristics on a High Speed Civil Transport (HSCT)-class configuration in the National Transonic Facility (NTF) are ongoing. The majority of the High Speed Research (HSR) data base in the NTF has free transition on the wing, even at low Reynolds numbers (Rn) attainable in conventional facilities. Limited data has been obtained and is described herein showing the effects of a conventional, Braslow method based wing boundary-layer trip on drag. Comparisons are made using force data polars and surface flow visualization at selected angles-of-attack and Mach number. Minimum drag data obtained in this study suggest that boundary layer transition occurred very near the wing leading edge by a chord Rn of 30 million. Sublimating chemicals were used in the air mode of operation only at low Rn and low angles-of-attack with no flap deflections; sublimation results suggest that the forebody and outboard wing panel are the only regions with significant laminar flow. The process and issues related to the sublimating chemical technique as applied in the NTF are discussed. Beyond the existing experience, status of efforts to develop a production transition detection system applicable to both air and cryogenic nitrogen environments is presented.

Owens, Lewis R., Jr.↗

NASA NTRS1987

Studies of pressure-velocity coupling schemes for analysis of incompressible and compressible flows

Two pressure-velocity coupling schemes, both of which solve the fully implicit discretized equations governing the flow of fluids were examined, and the capability of performing large Reynolds number, low Mach number compressible flow calculations were assessed. The semi-implicit iterative SIMPLE algorithm is extended to handle transient compressible flow calculations. This extension takes into account a strong coupling between the pressure and temperature through a correction procedure, based on the equations of state. Results obtained from the extended SIMPLE algorithm are then compared to similar results obtained from the non-iterative PISO algorithm. Both time-dependent and steady state calculations were performed using an axisymmetric 2:1 pipe expansion geometry and laminar flow conditions corresponding to Reynolds number of 1000 and Mach number of 2.0. For calculations simulating a time-dependent compression/expansion process, both schemes exhibit transient features in excellent agreement with each other, and moreover, the PISO method shows a significant computational time reduction of 60 percent over the SIMPLE scheme, regardless of the time step size or grid size employed. The effects of numerical diffusion are shown to be significant in these calculations. For steady state compressible flows, however, the SIMPLE algorithm displays increasing computational efficiency over the PISO method as the time step sizes employed to reach steady state conditions are decreased.

Burgreen, Gregory Wayne↗

NASA NTRS1988

On the two-dimensional modelling of horizontal CVD reactors and its limitations

In order to clarify under what conditions two-dimensional theoretical models can realistically account for the actual growth rate distribution in horizontal CVD reactors, a systematic numerical study was performed, and the results were compared with published experimental data for the MOCVD of gallium arsenide. Three models of increasing complexity were studied. It is shown that two-dimensional models can produce realistic predictions only for reactors with large width-to-height (aspect) ratios, that are operated at subcritical Rayleigh number,i.e., under dominant forced flow conditions. But even then, depending on the carrier gas, thermal (Soret) diffusion must be included in the model. Furthermore, velocity corrections for finite aspect ratios must be made, and buoyancy effects can be significant in the entrance region. On the other hand, the results are relatively insensitive to axial diffusion, temperature dependence of transport properties (versus averaged properties), and an assumed Poiseuille velocity profile (versus a rigorous solution of the velocity field). Yet, boundary layer approximations are prone to give unrealistic results at the low Reynolds numbers typical for CVD operations.

Ouazzani, Jalil↗

NASA NTRS1987

Arcjet thruster research and technology, phase 1

The objectives of Phase 1 were to evaluate analytically and experimentally the operation, performance, and lifetime of arcjet thrusters operating between 0.5 and 3.0 kW with catalytically decomposed hydrazine (N2H4) and to begin development of the requisite power control unit (PCU) technology. Fundamental analyses were performed of the arcjet nozzle, the gas kinetic reaction effects, the thermal environment, and the arc stabilizing vortex. The VNAP2 flow code was used to analyze arcjet nozzle performance with non-uniform entrance profiles. Viscous losses become dominant beyond expansion ratios of 50:1 because of the low Reynolds numbers. A survey of vortex phenomena and analysis techniques identified viscous dissipation and vortex breakdown as two flow instabilities that could affect arcjet operation. The gas kinetics code CREK1D was used to study the gas kinetics of high temperature N2H4 decomposition products. The arc/gas energy transfer is a non-equilibrium process because of the reaction rate constants and the short gas residence times. A thermal analysis code was used to guide design work and to provide a means to back out power losses at the anode fall based on test thermocouple data. The low flow rate and large thermal masses made optimization of a regenerative heating scheme unnecessary.

Knowles, Steven C.↗

NASA NTRS2004

Effects of Elevated Free-Stream Turbulence on Active Control of a Separation Bubble

The Effects of elevated free-stream turbulence (FST) on the natural and periodically excited separation bubbles were studied experimentally, due to the relevance of this flow to low-pressure turbine blades at low Reynolds numbers. A bubble was formed at the leading edge of a flat plate and the FST level was altered by placing a grid across the flow at different locations upstream of the plate. The mixing across the separated shear-layer, forming the free boundary of the bubble, increased due to the elevated FST and due to nominally two-dimensional periodic excitation, both flattening and shortening the bubble. Periodic excitation at frequencies that were at least an order of magnitude lower than those associated with the initial shear-layer instability, were very effective at low FST, because the amplitudes of the excitation frequency and its harmonic were amplified over the bubble. High frequency excitation (F+ 3, based on the length of the baseline low FST bubble) had a major effect close to the separation location, while farther downstream the excited fluctuations rapidly decayed in the reattachment region. Low frequency excitation, that generated waves comparable to the length of the unperturbed bubble (F+ 1) were less effective and their magnitude decayed at a slower rate downstream of reattachment. An increase in the level of the FST reduced the net effect of the periodic excitation on the mixing enhancement and subsequent reattachment process, probably due to a destructive interference between the nominally 2D excitation and the random (in space and time) FST, reducing the spanwise coherence and therefore the effectiveness of the current control strategy. However, even at the reduced effectiveness of 2D periodic excitation at elevated FST, it accelerated the reattachment process and the recovery rate of the reattached boundary layer, enhancing the boundary layer resistance to repeat separation and reducing its momentum loss further downstream.

Ashpis, D.↗

NASA NTRS1992

Experimental studies of hypersonic shock-wave boundary-layer interactions

Two classes of shock-wave boundary-layer interactions were studied experimentally in a shock tunnel in which a low Reynolds number, turbulent flow at Mach 8 was developed on a cold, flat test surface. The two classes of interactions were: (1) a swept interaction generated by a wedge ('fin') mounted perpendicularly on the flat plate; and (2) a two-dimensional, unseparated interaction induced by a shock impinging near an expansion corner. The swept interaction, with wedge angles of 5-20 degrees, was separated and there was also indication that the strongest interactions prossessed secondary separation zones. The interaction spread out extensively from the inviscid shock location although no indication of quasi-conical symmetry was evident. The surface pressure from the upstream influence to the inviscid shock was relatively low compared to the inviscid downstream value but it rose rapidly past the inviscid shock location. However, the surface pressure did not reach the downstream inviscid value and reasons were proposed for this anomalous behavior compared to strongly separated, supersonic interactions. The second class of interactions involved weak shocks impinging near small expansion corners. As a prelude to studying this interaction, a hypersonic similarity parameter was identified for the pure, expansion corner flow. The expansion corner severely damped out surface pressure fluctuations. When a shock impinged upstream of the corner, no significant changes to the surface pressure were found as compared to the case when the shock impinged on a flat plate. But, when the shock impinged downstream of the corner, a close coupling existed between the two wave systems, unlike the supersonic case. This close coupling modified the upstream influence. Regardless of whether the shock impinged ahead or behind the corner, the downstream region was affected by the close coupling between the shock and the expansion. Not only was the mean pressure distribution modified but the unsteadiness in the surface pressure was reduced compared to the flat-plate case.

Lu, Frank K.↗

NASA NTRSDate not supplied

Design Exploration of a Transonic Cruise Slotted Airfoil

A knowledge-based aerodynamic design method, CDISC, has been leveraged for the computational design exploration of the cruise slotted airfoil as a drag-saving technology for single aisle, transonic transport aircraft. Aerodynamic predictions were generated using the NASA USM3D Reynolds-averaged Navier-Stokes flow solver, and laminar flow assessments were conducted using the NASA BLSTA3D boundary layer profile solver paired with the LASTRAC stability analysis and transition prediction software. The CDISC design method was used to parametrically vary several slotted-airfoil design variables to understand their impact on aerodynamic performance. In addition to establishing best practices for slotted airfoil design, it was observed that the historically noted skin-friction penalty for cruise slotted airfoils was attributed to the low-Reynolds-number boundary layer of the flap. At cruise, this skin-friction penalty negated any potential decrease in pressure drag enabled by the cruise slotted airfoil architecture. This observation led to the development of an Aft-Laminar Multi-element Airfoil concept that uses airfoil shaping to promote natural laminar flow on the flap. At cruise conditions, the concept is predicted to offer a 2.7% reduction in sectional drag relative to a fully turbulent supercritical airfoil. Off-design analyses showed that laminar flow could be maintained with near-cruise variations in angle of attack and Mach number, resulting in sustained aerodynamic efficiency improvements and a delay in the drag rise Mach number. A low-speed, high-lift analysis at takeoff conditions predicted a maximum lift coefficient of 2.4 could be achieved by use of a variable-camber leading edge on the main element and simple deflection of the flap. These computationally-predicted benefits have motivated future research toward an aft-laminar cruise slotted wing design as an incremental step toward the application of natural laminar flow technology on transonic commercial transports.

Brett R Hiller↗

DOE OSTI2020

Analysis of Laminar Convective Heat Transfer Over Structured Non-Wetting Surfaces

Structured non-wetting surfaces provide alternating no-slip and partial slip boundary conditions to the fluid flow which, in turn, affects the convective heat transfer performance over the surfaces. In this work, an analytical model is developed for the interfacial Nusselt number, the overall Nusselt number and a thermal hydraulic performance factor for fluid flow in a cylinder patterned with structured non-wetting surfaces, for the two cases of uniform wall heat flux and uniform wall temperature. In addition, by considering the stability of the Cassie state of wettability and its transition to the Wenzel state for flow over superhydrophobic surfaces, the present model overcomes certain limitations of the previously reported studies in the literature. Based on the analytical formulations and the stability constraints, the present paper provides optimum design maps for tailoring structured non-wetting surfaces for maximizing convective heat transfer and the combined thermal-hydraulic performance in applications. Use of the design maps on example cases is also discussed. It is shown that the use of structured non-wetting surfaces is most effective for low Reynolds numbers and/or small cylinder radius.

42 ENGINEERING↗

NASA NTRS2019

Airfoil Selection for Mars Rotor Applications

This short white paper is aimed at providing an overview of design considerations for airfoil choices for rotor applications in the Martian atmosphere, at very low chord-based Reynolds number flows, around Rec = O (103 -104). The low Reynolds number typical of rotorcraft operation in the Martian atmosphere reduces the rotor lifting force and efficiency, which is only partially compensated for by a lower gravity on Mars. Additionally, the low temperature and largely CO2 based atmosphere of Mars compound the overall aerodynamic problem by resulting in a lower speed of sound, further constraining rotor operation in the Martian atmosphere by limiting the maximum rotor tip speed possible so as not to exceed an acceptable tip Mach number. In light of the expected reduced rotor efficiency, evaluating airfoils for compressible, low-Reynolds number Mars rotor applications is key. Prior research on airfoil optimization and performance evaluation at low Reynolds numbers, especially in the compressible regime, is scarce and deserves further understanding. Specifically, the proposed goal stemming from this overview is to develop airfoils tailored to the unique demands of the second generation of Mars rotorcraft, i.e. the Mars Science Helicopter (MSH). This research focuses on the airfoil performance at low Reynolds numbers and hopes to add to the work performed by, amongst others, Kroo et al., Kunz and Kroo, Oyama and Fujii, Anyoji et al.

Koning, Witold J. F.↗

NASA NTRS1999

Prediction of Heat and Mass Transfer in a Rotating Ribbed Coolant Passage With a 180 Degree Turn

Numerical results are presented for flow in a rotating internal passage with a 180 degree turn and ribbed walls. Reynolds numbers ranging from 5200 to 7900, and Rotation numbers of 0.0 and 0.24 were considered. The straight sections of the channel have a square cross section, with square ribs spaced one hydraulic diameter (D) apart on two opposite sides. The ribs have a height of 0.1D and are not staggered from one side to the other. The full three dimensional Reynolds Averaged Navier-Stokes equations are solved combined with the Wilcox k-omega turbulence model. By solving an additional equation for mass transfer, it is possible to isolate the effect of buoyancy in the presence of rotation. That is, heat transfer induced buoyancy effects can be eliminated as in naphthalene sublimation experiments. Heat transfer, mass transfer and flow field results are presented with favorable agreement with available experimental data. It is shown that numerically predicting the reattachment between ribs is essential to achieving an accurate prediction of heat/mass transfer. For the low Reynolds numbers considered, the standard turbulence model did not produce reattachment between ribs. By modifying the wall boundary condition on omega, the turbulent specific dissipation rate, much better agreement with the flow structure and heat/ mass transfer was achieved. It is beyond the scope of the present work to make a general recommendation on the omega wall boundary condition. However, the present results suggest that the omega boundary condition should take into account the proximity to abrupt changes in geometry.

Rigby, David L.↗

DOE OSTI2026

Modeling Flow in a Microreactor Core: From High-Fidelity CFD to Subchannel Analysis

The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

NASA NTRS2000

SHIVA (Spaceflight Holography Investigation in a Virtual Apparatus)

This paper provides the description and results of a ground-based experiment designed to support the spaceflight program SHIVA, which started its flight definition in 1998. SHIVA will apply new diagnostic tools and experimental techniques to test the validity of a newly discovered analytical solution to the general equation of motion of a particle in a fluid. We have designed a spaceflight experiment with the help of the theoretical model that is optimized for testing the model, and for measuring g, g-jitter, and other microgravity phenomena. Our ongoing, ground-based particle/fluid experiment supports both the experimental and theoretical aspects of the project. The ultimate spaceflight experiment will be similar to the ground-based experiment. The "virtual spaceflight chamber" concept asserts that certain spaceflight experiments can be recorded in holograms in such a manner that having the holograms on earth is optically equivalent to being back in space with unlimited time to conduct the experiment. Properly exploited, this concept can save a significant amount of experiment time in space by effectively bringing the. experiment optically back to earth. SHIVA will accomplish the following: record a large number of holograms of particle fields in space under controlled conditions, extract the precise, three-dimensional position of all of the particles as a function of time, examine the effects of all parameters on the motion of the particles, and test these against predictions of the Coimbra-Range1 solution to the general equation of motion. Particle sizes will range from hundreds of microns up to about 2 mm in diameter and will cover a range of densities and fluid viscosities. Forcing functions will be introduced onto the particle field, including at least the following: a) No isolation in the Spacelab vibration environment. b) Isolation from Spacelab. c) Oscillatory motion from 1 to 100 Hz with amplitudes of a few millimeters. d) Convective fields to be introduced mechanically. The first objective of this research is to understand the physics of particle interactions with fluids and other particles in low Reynolds number flows in microgravity. Secondary objectives are to: (1) observe and quantify g-jitter effects and microconvection on particles in fluids, (2) validate an exact solution to the general equation of motion of a particle in a fluid, and (3) to characterize'the ability of isolating platforms to isolate experiments containing particles in liquids.

Trolinger, James D.↗

DOE OSTI2023

Hydrodynamic irreversibility of non-Brownian suspensions in highly confined duct flow

The irreversible behaviour of a highly confined non-Brownian suspension of spherical particles at low Reynolds number in a Newtonian fluid is studied experimentally and numerically. In the experiment, the suspension is confined in a thin rectangular channel that prevents complete particle overlap in the narrow dimension and is subjected to an oscillatory pressure-driven flow. In the small cross-sectional dimension, particles rapidly separate to the walls, whereas in the large dimension, features reminiscent of shear-induced migration in bulk suspensions are recovered. Furthermore, as a consequence of the channel geometry and the development and application of a single-camera particle tracking method, three-dimensional particle trajectories are obtained that allow us to directly associate relative particle proximity with the observed migration. Companion simulations of a steadily flowing suspension highly confined between parallel plates are conducted using the force coupling method, which also show rapid migration to the walls as well as other salient features observed in the experiment. While we consider relatively low volume fractions compared to most prior work in the area, we nevertheless observe significant and rapid migration, which we attribute to the high degree of confinement.

42 ENGINEERING↗

NASA NTRS2014

Simulation of VSPT Experimental Cascade Under High and Low Free-Stream Turbulence Conditions

Variable-Speed Power Turbines (VSPT) for rotorcraft applications operate in conditions of low Reynolds number and a wide range in incidence resulting from rotational speed variation. A comprehensive data set obtained in a linear cascade which includes the effects of Reynolds number, free-stream turbulence and incidence is now available and this paper concerns itself with the post-diction of boundary layer transitionseparation, blade pressure loading and total pressure loss pertaining to the conditions set for measurements in that data set. The distinction between the state of the measured data presented here and the earlier publications is the addition of high free-stream turbulence intensity. We will, for the purposes of the numerical post-diction, present some of the higher free stream turbulence data in this paper but defer a comprehensive presentation and treatment of the measured data will be done elsewhere.

power turbine↗