Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence Stabilization”

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 937 records · Page 52

Experiments on the flow and acoustic properties of a moderate-Reynolds-number supersonic jet

Flow and acoustic properties of a jet at Reynolds number of 70,000 were studied at Mach 2.1. Measurements in a free jet test facility were made with pitot tubes and hot-wire anemometry. Center-line Mach number distributions for natural and excited jets were obtained. A slow initial growth rate was in the potential core region of the jet, indicating a transition from laminar to turbulent flow in moderate Reynolds number jets. The transition occurred within the first 2-3 diameters. Spectral components were calculated for the fluctuating flowfield, and sound pressure levels were measured for the overall near-field noise. The centroid of noise was located about 8 nozzle diameters downstream. The growth rates of instabilities were determined to be in agreement with linear stability theory predictions over a broad frequency range.

Troutt, T. R.↗

A quiet flow Ludwieg tube for study of transition in compressible boundary layers: Design and feasibility

Laminar-turbulent transition in high speed boundary layers is a complicated problem which is still poorly understood, partly because of experimental ambiguities caused by operating in noisy wind tunnels. The NASA Langley experience with quiet tunnel design has been used to design a quiet flow tunnel which can be constructed less expensively. Fabrication techniques have been investigated, and inviscid, boundary layer, and stability computer codes have been adapted for use in the nozzle design. Construction of such a facility seems feasible, at a reasonable cost. Two facilities have been proposed: a large one, with a quiet flow region large enough to study the end of transition, and a smaller and less expensive one, capable of studying low Reynolds number issues such as receptivity. Funding for either facility remains to be obtained, although key facility elements have been obtained and are being integrated into the existing Purdue supersonic facilities.

Schneider, Steven P.↗

High speed transition prediction

The main objective of this work period was to develop, acquire and apply state-of-the-art tools for the prediction of transition at high speeds at NASA Ames. Although various stability codes as well as basic state codes were acquired, the development of a new Parabolized Stability Equation (PSE) code was minimal. The time that was initially allocated for development was used on other tasks, in particular for the Leading Edge Suction problem, in acquiring proficiency in various graphics tools, and in applying these tools to evaluate various Navier-Stokes and Euler solutions. The second objective of this work period was to attend the Transition and Turbulence Workshop at NASA Langley in July and August, 1991. A report on the Workshop follows. From July 8, 1991 to August 2, 1991, the author participated in the Transition and Turbulence Workshop at NASA Langley. For purposes of interest here, analysis can be said to consist of solving simplified governing equations by various analytical methods, such as asymptotic methods, or by use of very meager computer resources. From the composition of the various groups at the Workshop, it can be seen that analytical methods are generally more popular in Great Britain than they are in the U.S., possibly due to historical factors and the lack of computer resources. Experimenters at the Workshop were mostly concerned with subsonic flows, and a number of demonstrations were provided, among which were a hot-wire experiment to probe the boundary layer on a rotating disc, a hot-wire rake to map a free shear layer behind a cylinder, and the use of heating strips on a flat plate to control instability waves and consequent transition. A highpoint of the demonstrations was the opportunity to observe the rather noisy 'quiet' supersonic pilot tunnel in operation.

Gasperas, Gediminis↗

A quiet-flow Ludwieg tube for experimental study of high speed boundary layer transition

Laminar-turbulent transition in high speed boundary layers is a complicated problem which is still poorly understood, partly because of experimental ambiguities caused by operating in noisy wind tunnels. The NASA Langley experience with quiet tunnel design has been used to design a quiet flow tunnel which can be constructed less expensively. Fabrication techniques have been investigated, and inviscid, boundary layer, and stability computer codes have been adapted for use in the nozzle design. Construction of such a facility seems feasible, at a reasonable cost. Two facilities have been proposed: a large one, with a quiet flow region large enough to study the end of transition, and a smaller and less expensive one, capable of studying low Reynolds number issues such as receptivity. Funding for either facility remains to be obtained, although key facility elements have been obtained and are being integrated into the existing Purdue supersonic facilities.

Schneider, Steven P.↗

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗

The Influence of Low Wall Temperature on Boundary-Layer Transition and Local Heat Transfer on 2-Inch-Diameter Hemispheres at a Mach Number of 4.95 and a Reynolds Number per Foot of 73.2 x 10(exp 6)

Measurements of the location of boundary-layer transition and the local heat transfer have been made on 2-inch-diameter hemispheres in the Langley gas dynamics laboratory at a Mach number of 4.95, a Reynolds number per foot of 73.2 x 10(exp 6), and a stagnation temperature of approximately 400 F. The transient-heating thin-skin calorimeter technique was used, and the initial values of the wall-to-stream stagnation- temperature ratios were 0.16 (cold-model tests) and 0.65 (hot-model test). During two of the four cold tests, the boundary-layer flow changed from turbulent to laminar over large regions of the hemisphere as the model heated. On the basis of a detailed consideration of the magnitude of roughness possibly present during these two cold tests, it appears that this destabilizing effect of low wall temperatures (cooling) was not caused by roughness as a dominant influence. This idea of a decrease in boundary-layer stability with cooling has been previously suggested. (See, for example, NASA Memorandum 10-8-58E.) For the laminar data obtained during the early part of the hot test, the correlation of the local-heating data with laminar theory was excellent.

Cooper, Morton↗

Asymptotic Behavior of Ensemble-Averaged Linear Disturbances in Homogeneous Shear Flow

In order to expand the predictive capability of single-point turbulence closure models too account for the early-stage transition regime, a methodology for the formulation and calibration of model equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is presented. The calibration is based on homogeneous shear flow where disturbances can be described by rapid distort,ion theory (RDT). The relationship between RDT and linear stability theory is exploit,c d in order to obtain a closed set, of modeled equations. The linear disturbance equations are solved directly so that, the numerical simulation yields a database from which the closure coefficient,s in the ensemble-averaged disturbance equations can he determined.

Thacker, W. D.↗

Numerical Simulation of Receptivity for a Transition Experiment

The cost of fuel to overcome turbulence induced viscous drag on a commercial airplane constitutes a significant fraction of the operating cost of an airline. Achieving laminar flow and maintaining it over a large portion of the wing can significantly reduce the viscous drag, and hence the cost. Design of such laminar-flow-control wings and their practical operation requires the ability to accurately and reliably predict the transition from laminar to turbulent flow. The transition process begins with the conversion of environmental and surface disturbances into the instability waves of the flow by a process called receptivity. The goal of the current research project has been to improve the prediction of transition through a better understanding of the physics of receptivity. The initial objective of this work was to investigate the specific stability and receptivity characteristics of a particular experimental investigation of boundary layer receptivity at NASA Langley. Some simulation results using direct solutions of the linearized Navier-Stokes equations which modeled this experiment where presented in the 1999 APS DFD meeting. However, based on these initial investigations, it became clear that to cover the vast receptivity parameter space required for a practical transition prediction tool, more efficient methods would be required. Thus, the focus of this research was shifted from modeling this particular experiment to formulating and developing new techniques that could efficiently yet accurately predict receptivity for a wide range of disturbance conditions.

Collis, S. Scott↗

Verification of Anisotropic Mesh Adaptation for Turbulent Simulations over ONERA M6 Wing

Unstructured anisotropic mesh adaptation is known to be an efficient way to control discretization errors in Computational Fluid Dynamics (CFD) simulations. Method verification is required to provide the confidence for routine use in production analysis. The current work aims at verification of anisotropic mesh adaptation for RANS simulations over the ONERA M6 wing. The present verification study is performed using four different flow solvers, three different implementations of the metric field, and three mesh mechanics packages. Two of the flow solvers use stabilized finite-element discretizations (FUN3D-SFE and GGNS), one uses finite-volume discretization (FUN3D-FV), and the last one uses mixed finite-volume and finite element discretizations (Wolf). The mesh adaptation is based on an error estimator that aims to control the quadratic error term in the linear interpolation of Mach number. Two sets of adaptations were performed; the first one controls the interpolation error in L2 norm and the second one controls the interpolation error in L4 norm. Convergence studies were performed on the forces and the pitching moment using all four solvers, and the results are compared with previously verified convergence studies on fixed (nonadapted) meshes. Both forces and pitching moment on adapted meshes are found to be converging to the fine mesh values faster than those on fixed meshes. In addition to forces and moments, convergence of surface pressure and skin friction coefficients at various measurement locations on the wing are also presented. Adapted-mesh surface pressure distributions agree with the fine fixed mesh pressure distributions. Adapted-mesh skin friction distributions contain high frequency noise with mean values approaching the fixed mesh pressure skin friction distributions.

Aravind Balan↗

A scheme for computing surface layer turbulent fluxes from mean flow surface observations

A physical model and computational scheme are developed for generating turbulent surface stress, sensible heat flux and humidity flux from mean velocity, temperature and humidity at some fixed height in the atmospheric surface layer, where conditions at this reference level are presumed known from observations or the evolving state of a numerical atmospheric circulation model. The method is based on coupling the Monin-Obukov surface layer similarity profiles which include buoyant stability effects on mean velocity, temperature and humidity to a force-restore formulation for the evolution of surface soil temperature to yield the local values of shear stress, heat flux and surface temperature. A self-contained formulation is presented including parameterizations for solar and infrared radiant fluxes at the surface. Additional parameters needed to implement the scheme are the thermal heat capacity of the soil per unit surface area, surface aerodynamic roughness, latitude, solar declination, surface albedo, surface emissivity and atmospheric transmissivity to solar radiation.

Hoffert, M. I.↗

Modern developments in shear flow control with swirl

Passive and active control of swirling turbulent jets is experimentally investigated. Initial swirl distribution is shown to dominate the free jet evolution in the passive mode. Vortex breakdown, a manifestation of high intensity swirl, was achieved at below critical swirl number (S = 0.48) by reducing the vortex core diameter. The response of a swirling turbulent jet to single frequency, plane wave acoustic excitation was shown to depend strongly on the swirl number, excitation Strouhal number, amplitude of the excitation wave, and core turbulence in a low speed cold jet. A 10 percent reduction of the mean centerline velocity at x/D = 9.0 (and a corresponding increase in the shear layer momentum thickness) was achieved by large amplitude internal plane wave acoustic excitation. Helical instability waves of negative azimuthal wave numbers exhibit larger amplification rates than the plane waves in swirling free jets, according to hydrodynamic stability theory. Consequently, an active swirling shear layer control is proposed to include the generation of helical instability waves of arbitrary helicity and the promotion of modal interaction, through multifrequency forcing.

Farokhi, Saeed↗

Modern developments in shear flow control with swirl

Passive and active control of swirling turbulent jets is experimentally investigated. Initial swirl distribution is shown to dominate the free jet evolution in the passive mode. Vortex breakdown, a manifestation of high intensity swirl, was achieved at below critical swirl number (S = 0.48) by reducing the vortex core diameter. The response of a swirling turbulent jet to single frequency, plane wave acoustic excitation was shown to depend strongly on the swirl number, excitation Strouhal number, amplitude of the excitation wave, and core turbulence in a low speed cold jet. A 10 percent reduction of the mean centerline velocity at x/D = 9.0 (and a corresponding increase in the shear layer momentum thickness) was achieved by large amplitude internal plane wave acoustic excitation. Helical instability waves of negative azimuthal wave numbers exhibit larger amplification rates than the plane waves in swirling free jets, according to hydrodynamic stability theory. Consequently, an active swirling shear layer control is proposed to include the generation of helical instability waves of arbitrary helicity and the promotion of modal interaction, through multifrequency forcing.

Farokhi, S.↗

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the nonlocal, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, such as those in the vicinity of roughness elements, steps, gaps, or corners. To compute the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 10^(7) to 10^(8) even in relatively simple flows. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear systems. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE. Additionally, inflow-resolvent analysis based on the HLNSE for discrete-roughness-induced disturbances in the nose tip of a blunt cone at Mach 6 demonstrates the importance of including the disturbance amplification along the near vicinity of the roughness element and separation region.

Boundary Layer Stability↗

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the non-local, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, as roughness elements, steps, gaps, or corners. To solve the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 107 to 108. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear system. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE.

Boundary Layer Stability↗

Impact of boundary layer simulation on predicting radioactive pollutant dispersion: a case study for HANARO Research Reactor using the WRF-MMIF-CALPUFF modeling system

Wind plays an important role in cases of unexpected radioactive pollutant dispersion, deciding distribution and concentration of the leaked substance. The accurate prediction of wind has been challenging in numerical weather prediction models, especially near the surface because of the complex interaction between turbulent flow and topographic effect. As such, in this study, we investigated the characteristics of atmospheric dispersion of radioactive material (i.e. 137 Cs) according to the simulated boundary layer around the HANARO research nuclear reactor in Korea using the Weather Research and Forecasting (WRF)-Mesoscale Model Interface (MMIF)-California Puff (CALPUFF) model system. We examined the impacts of orographic drag on wind field, stability calculation methods, and planetary boundary layer parameterizations on the dispersion of radioactive material under a radioactive leaking scenario. We found that inclusion of the orographic drag effect in the WRF model improved the wind prediction most significantly over the complex terrain area, leading the model system to estimate the radioactive concentration near the reactor more conservatively. We also emphasized the importance of the stability calculation method and employing the skillful boundary layer parameterization to ensure more accurate low atmospheric conditions, in order to simulate more feasible spatial distribution of the radioactive dispersion in leaking scenarios.

54 ENVIRONMENTAL SCIENCES↗

Simultaneous 10 kHz three-dimensional CH 2 O and tomographic PIV measurements in a lifted partially-premixed jet flame

High-speed, three-dimensional (3D) scalar-velocity field measurements are demonstrated in a lifted partially-premixed dimethyl-ether/air jet flame using simultaneous laser-induced fluorescence (LIF) of formaldehyde (CH 2 O) and tomographic particle image velocimetry (TPIV). The 3D LIF measurements are performed by raster scanning the laser beam from a 100 kHz pulse-burst laser across the probe volume using an acousto-optic deflector. The volumetric reconstruction of the LIF signal from ten parallel planes provides quasi-instantaneous 3D LIF measurements that are synchronized with 10 kHz TPIV measurements. Here, the temporally resolved CH 2 O-LIF and velocity field data are used to analyze Lagrangian particle trajectories and displacement speeds at the base of the lifted flame. The particle trajectories reveal flow structures that are difficult to observe in an Eulerian reference frame. Positive and negative displacement speeds are observed at the CH 2 O-LIF surfaces at the inner and outer regions of the jet flame with a maximum displacement speed of approximately eight times the laminar flame speed of a stoichiometric dimethyl-ether/air mixture. The observations in present study support formation of lean premixed flames near the flame base whose role in flame stabilization should be investigated in more detail.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

RF current condensation in the presence of turbulent enhanced transport

Sharp temperature gradients in a magnetically confined plasma can lead to turbulent motion of the plasma. This turbulence in turn enhances the transport of heat across magnetic field lines. The enhanced transport impacts the temperature differential that can be sustained in magnetic islands between the island center and its periphery. It is shown in this work that, by limiting this temperature differential, this enhanced transport can have a profound influence on the extent to which the RF current condensation effect stabilizes the island growth. Interestingly, because the heat transport is no longer simply linear in the temperature gradient, the RF current condensation effect also exhibits entirely new hysteresis phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of Surface Roughness on Dynamic Stall in Pitching Motion

Dynamic stall plays a critical role in determining the performance and stability of a wide range of fluid-dynamics systems in various engineering applications. This unsteady aerodynamic phenomenon is particularly significant for maneuvering aircraft wings, jet aircraft subjected to gust encounters, helicopter rotor blades, and wind turbine blades [1–3]. The prediction of the dynamic stall vortex (DSV) is challenging due to factors such as unsteady aerodynamics, three-dimensional (3-D) effects, turbulence and flow separation, incoming gust, and surface impact effects [4–6]. Hence, advanced computational techniques and modeling approaches in computational fluid dynamics (CFD) would be required to enhance the accuracy and reliability of DSV predictions in dynamic stall scenarios. In the past, Batther and Lee [7] employed delayed detached eddy simulations (DDES) to understand the flow physics associated with the onset of dynamic stall. Their approach demonstrated that DDES achieves results comparable to those obtained from large-eddy simulations at a reduced computational cost. In another study, Khalifa et al. [8] examined the 3-D aspects of dynamic stall on a NACA 0012 airfoil using DES solvers. In conclusion, the findings underlined the superiority of 3-D simulations over two-dimensional approaches, particularly in predicting the lift coefficient values and capturing dynamic stall stages more precisely.

97 MATHEMATICS AND COMPUTING↗