SEARCH · Engineering Papers
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.
Turbulent, nonisothermal, viscous, compressible flow in inlet sections of axisymmetric flat channels
A time-accurate finite-difference Navier-Stokes code has been used to calculate the viscous flow over a severely indented blunt body in a supersonic stream. An algebraic turbulence model is used and the results are compared with experimental data from wind-tunnel tests. Qualitative agreement is obtained for the surface pressure distribution and flow-field structure, including the separated bubble in the indented region. However, uncertainties still exist in the heating calculations, which are attributed to the turbulence model. For both laminar and turbulent calculations, the flow exhibits a fundamental unsteady character at a frequency of about 50 kHz.
A computationally-viable model describing the interaction between fluid-mechanical turbulence and finite-rate combustion reactions, principally in high-speed flows was developed. Chemical kinetic mechanisms, complete and global, were developed describing the finite rate reaction of fuels of interest to NASA. These fuels included principally hydrogen and silane, although a limited amount of work involved hydrocarbon fuels as well.
A new gauge for skin-friction measurements in turbulent flows involving a wire close to the surface and a second, flush-mounted film which serves as a constant temperature guard heater is described and evaluated. It is shown that a laminar calibration of the gauge applies to measurements in turbulent flow and that the frequency response of the gauge is appropriate for time-resolved data in such a flow. Initial results from a turbulent channel flow are presented.
Deep neural networks (DNNs) are developed from a data set obtained from the dynamic Smagorinsky model to emulate the subgrid-scale (SGS) viscosity (ν sgs ) and diffusivity (κ sgs ) for turbulent stratified shear flows encountered in the oceans and the atmosphere. These DNNs predict ν sgs and κ sgs from velocities, strain rates, and density gradients such that the evolution of the kinetic energy budget and density variance budget terms is similar to the corresponding values obtained from the original dynamic Smagorinsky model. These DNNs also compute ν sgs and κ sgs ~2–4 times quicker than the dynamic Smagorinsky model resulting in a ~2–2.5 times acceleration of the entire simulation. This study demonstrates the feasibility of deep learning in emulating the subgrid-scale (SGS) phenomenon in geophysical flows accurately in a cost-effective manner. In a broader perspective, deep learning-based surrogate models can present a promising alternative to the traditional parameterizations of the subgrid-scale processes in climate models.
Hot-wire anemometer and its application to measurements in turbulent flows
There exists currently considerable interest in the development of numerical procedures for solving the equations governing three-dimensional fluid flows. In particular, improved calculations of viscous flows on rotating bodies such a turbomachinery rotors are of great importance. The present investigation is mainly concerned with the analysis of the effect of rotation on the dissipation rate equation and the Reynolds stress equations. Attention is given to governing equations and modelling, algebraic Reynolds stress equations, a procedure for implementing the turbulence closure scheme, results for a model representative of centrifugal turbomachinery, and results for axial turbomachinery. The proposed turbulence model is valid for three-dimensional, rotating, and curved flows. The model has all the essential physics of turbulence, mean flow, rotation, and the curvature.
A transport equation for the root-mean-square pressure fluctuations of turbulent flow is derived from the time-dependent momentum equation for incompressible flow. Approximate modeling of this transport equation is included to relate terms with higher order correlations to the mean quantities of turbulent flow. Three empirical constants are introduced in the model. Two of the empirical constants are estimated from homogeneous turbulence data and wall pressure fluctuations measurements. The third constant is determined by comparing the results of large eddy simulations for a plane channel flow and an annulus flow.
Nonintrusive measurements of density and temperature and their turbulent fluctuation levels have been obtained in the boundary layer of an unseeded, Mach 2 wind tunnel flow. The spectroscopic technique that was used to make the measurements is based on the combination of laser-induced oxygen fluorescence and Raman scattering by oxygen and nitrogen from the same laser pulse. Results from this demonstration experiment compare favorably with previous measurements obtained in the same facility from conventional probes and an earlier spectroscopic technique.
It is known that the simulations of solid rocket motor internal flow field with AL-based propellants require complex multi-phase turbulent flow model. The objective of this study is to develop an advanced particulate multi-phase flow model which includes the effects of particle dynamics, chemical reaction and hot gas flow turbulence. The inclusion of particle agglomeration, particle/gas reaction and mass transfer, particle collision, coalescence and breakup mechanisms in modeling the particle dynamics will allow the proposed model to realistically simulate the flowfield inside a solid rocket motor. The Finite Difference Navier-Stokes numerical code FDNS is used to simulate the steady-state multi-phase particulate flow field for a 3-zone 2-D axisymmetric ASRM model and a 6-zone 3-D ASRM model at launch conditions. The 2-D model includes aft-end cavity and submerged nozzle. The 3-D model represents the whole ASRM geometry, including additional grain port area in the gas cavity and two inhibitors. FDNS is a pressure based finite difference Navier-Stokes flow solver with time-accurate adaptive second-order upwind schemes, standard and extended k-epsilon models with compressibility corrections, multi zone body-fitted formulations, and turbulence particle interaction model. Eulerian/Lagrangian multi-phase solution method is applied for multi-zone mesh. To simulate the chemical reaction, penalty function corrected efficient finite-rate chemistry integration method is used in FDNS. For the AL particle combustion rate, the Hermsen correlation is employed. To simulate the turbulent dispersion of particles, the Gaussian probability distribution with standard deviation equal to (2k/3)(exp 1/2) is used for the random turbulent velocity components. The computational results reveal that the flow field near the juncture of aft-end cavity and the submerged nozzle is very complex. The effects of the turbulent particles affect the flow field significantly and provide better prediction of the ASRM performance. The multi-phase flow analysis using the FDNS code in the present research can be used as a design tool for solid rocket motor applications.
A tip section of a delta wing having an NACA 0004-65 airfoil section and a 600 leading-edge sweepback was equipped with a porous leading-edge section through which hot gas was 'bled for anti-icing. Heating rates for anti-icing were determined for a wide range of icing conditions. The effects of gas flow through the porous leading-edge section on airfoil pressure distribution and drag in dry air were investigated. The drag increase caused by an ice formation on the unheated airfoil was measured for several icing conditions. Experimental porous surface- to free-stream convective heat-transfer coefficients were obtained in dry air and compared with theory. Adequate icing protection was obtained at all icing conditions investigated. Savings in total gas-flow rate up to 42 percent may be obtained with no loss in anti-icing effectiveness by sealing half the upper-surface porous area. Gas flow through the leading-edge section had no appreciable effect on airfoil pressure distribution. The airfoil section drag increased slightly (5-percent average) with gas flow through the porous surface. A heavy glaze-ice formation produced after 10 minutes of icing caused an increase in section drag coefficient of 240 percent. Experimental convective heat-transfer coefficients obtained with hot-gas flow through the porous area in dry air and turbulent flow were 20 to 30 percent lower than the theoretical values for a solid surface under similar conditions. The transition region from laminar to turbulent flow moved forward as the ratio of gas velocity through the porous surface to air-stream velocity was increased.
One of the main themes in fluid dynamics at present and in the future is going to be computational fluid dynamics with the primary focus on the determination of drag, flow separation, vortex flows, and unsteady flows. A computation of the flow of a viscous fluid requires an understanding and consideration of the physical aspects of the flow. This is done by identifying the flow regimes and the scales of fluid motion, and the sources of vorticity. Discussions of flow regimes deal with conditions of incompressibility, transitional and turbulent flows, Navier-Stokes and non-Navier-Stokes regimes, shock waves, and strain fields. Discussions of the scales of fluid motion consider transitional and turbulent flows, thin- and slender-shear layers, triple- and four-deck regions, viscous-inviscid interactions, shock waves, strain rates, and temporal scales. In addition, the significance and generation of vorticity are discussed. These physical aspects mainly guide computations of the flow of a viscous fluid.
A three dimensional finite element computer code was developed to analyze ejector and axisymmetric fluted mixer systems whose flow fields are not significantly influenced by streamwise diffusion effects. A two equation turbulence model was used to make comparisons between theory and data for various flow fields which are components of the ejector system, i.e., (1) turbulent boundary layer in a duct; (2) rectangular nozzle (free jet); (3) axisymmetric nozzle (free jet); (4) hypermixing nozzle (free jet); and (5) plane wall jet. Likewise, comparisons of the code with analytical results and/or other numerical solutions were made for components of the axisymmetric fluted mixer system. These included: (1) developing pipe flow; (2) developing flow in an annular pipe; (3) developing flow in an axisymmetric pipe with conical center body and no fluting and (4) developing fluted pipe flow. Finally, two demonstration cases are presented which show the code's ability to analyze both the ejector and axisymmetric fluted mixers.
In nuclear sodium fast reactors, bundles of rods are tightly packed into a triangular lattice, enclosed in a hexagonal duct, and each pin is spirally wrapped with a thin wire. Flow blockages can potentially impact the local flow characteristics and heat transfer mechanisms in the bundle due to its small subchannel size. The effects of the blockage on the flow structures and heat transfer mechanisms are important aspects that require an accurate investigation. In this study, the flow-field characteristics in the vicinity of a blockage located in the exterior subchannel of rod bundles with helically wrapped wires were experimentally investigated. The velocity fields in the exterior subchannel were acquired by applying matched-index-of-refraction and time-resolved particle image velocimetry (TR-PIV) techniques for Reynolds numbers of Re 1 = 4000 and Re 2 = 17 000, i.e., equivalent to Rew 1 = 19 600 and Rew 2 = 83 200, respectively, based on the blockage width. The results from the TR-PIV measurements revealed an arch-shaped vortex with a large flow recirculation and a pair of counter-rotating vortices in the wake region downstream of the blockage, which is commonly observed in the wake flow of bluff bodies. The relative lateral distance and angle between the two vortices decreased when the Reynolds numbers increased. Profiles of maximum turbulence intensity along the shear layers illustrated the transition process including the growth, peak, and decay along the flow direction. From the spectral analysis of the turbulent velocities extracted at points along the shear layer, the Strouhal numbers ( St ) representing the vortex shedding frequency were found to be St = 0.25 and St = 0.56 for the left and right shear layers, respectively. Characteristics of shear layers generated by the blockage in the exterior subchannel were investigated via the two-point cross correlation of fluctuating velocities. The patiotemporal cross correlations of turbulent velocities, computed at points in the region where the left shear layer exhibited rolling effects and vortex breakdowns, were considerably wider and longer. The convection velocity U c was estimated to be ~0.82 U m to 0.93 U m . Proper orthogonal decomposition (POD) analysis was applied to the instantaneous velocity fields to extract the statistically dominant flow structures. It was found that POD modes 2–3 and 4–5 formed the pair modes when the corresponding POD temporal coefficients depicted sinusoidal shapes and exhibited nearly circular orbits in the phase space. Spectral analysis of the POD temporal coefficients confirmed the vortex shedding frequencies detected in the analysis of turbulent velocities.
The results of laser anemometer and flow visualization based fluid mechanics studies of porous media are used to obtain heat transfer coefficients for porous materials. Average pore flow Re ranging from 0.16-700 were examined. Darcy, inertial steady laminar, unsteady laminar and turbulent flow regimes were detected. A passage length model was devised to derive the heat transfer coefficient. Sample data from flows through porous metals composed of powders and fibers validated the passage length for Darcy and inertial flow regimes. Unsteady laminar and turbulent flow coefficients require the identification of new parameters.
Rapid-distortion theory is adapted to introduce a truly unsteady closure into a simple phenomenological turbulence model in order to describe the unsteady response of a turbulent wall layer exposed to a temporarily oscillating pressure gradient. The closure model is built by taking the ratio of turbulent shear stress to turbulent kinetic energy to be a function of the effective strain. The latter accounts for the history of the flow. The computed unsteady velocity fluctuations and modulated turbulent stresses compare favorably in the 'non-quasi-steady' frequency range, where quasi-steady assumptions would fail. This suggests that the concept of rapid distortion is especially appropriate for unsteady flows. This paper forms the basis for acoustical studies of the problem to be reported elsewhere.
A cylinder and a thin plate were placed close together in a flow adjacent to a wall to study the effects on the turbulent boundary layer. Different spacings of the cylinder and plate within the shear flow were investigated to assess the possibility of lowering the production of fluctuating vorticity in the boundary layer by generating a fluctuating vorticity of opposite sign. Streakline photographs visualized changes in the flow induced by alterations in the cylinder/plate separation distance, the flow velocity and the angle of attack of the thin plate. Drag data were also acquired with varying thicknesses of the thin plate and diameters of the cylinder. Downstream skin friction reductions were obtained with the production of unsteady control vortices with a low turbulence boundary layer. Up to 4 percent drag reduction was also observed when the cylinder was sufficiently far from the wall.
Measurements of streamwise velocity fluctuation and surface pressure fluctuation spectra and wavespeeds are reported for a well-documented separating turbulent boundary layer. Because a portion of the acoustic pressure fluctuations is the same across the nominally two-dimensional turbulent flow, it is possible to decompose two microphone signals and obtain directly the turbulent flow contributions to the surface pressure spectra. The rms surface pressure fluctuation p' and spectra phi(omega) increase through the adverse pressure gradient attached flow region and the detached flow zone and scale on the maximum turbulent shearing stress tau(M); p'/tau(M) increases to the detachment location and decreases downstream due to the rapid movement of the pressure-fluctuation-producing motions away from the wall after the beginning of intermittent backflow. At lower frequencies for the attached flow phi(omega) is approximately omega to the -0.7 while phi(omega) is approximately omega to the -3 at higher frequencies. After the beginning of intermittent backflow, phi(omega) varies with omega at low frequencies and omega to the -3 at high frequencies; farther downstream the lower frequency range varies with omega to the 2.4. The surface pressure fluctuation celerity for the attached flow increases with frequency and agrees with the semi-logarithmic overlap equation of Panton and Linebarger. After the beginning of the separation process, the wavespeed decreases because of the oscillation of the instantaneous wavespeed direction and the streamwise coherence decreases drastically.