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 235 records · Page 13

An Experimental Investigation of Unsteady Surface Pressure on an Airfoil in Turbulence

Measurements of fluctuating surface pressure were made on a NACA 0015 airfoil immersed in grid generated turbulence. The airfoil model has a 2 ft chord and spans the 6 ft Virginia Tech Stability Wind Tunnel test section. Two grids were used to investigate the effects of turbulence length scale on the surface pressure response. A large grid which produced turbulence with an integral scale 13% of the chord and a smaller grid which produced turbulence with an integral scale 1.3% of the chord. Measurements were performed at angles of attack, alpha from 0 to 20 . An array of microphones mounted subsurface was used to measure the unsteady surface pressure. The goal of this measurement was to characterize the effects of angle of attack on the inviscid response. Lift spectra calculated from pressure measurements at each angle of attack revealed two distinct interaction regions; for omega(sub r) = omega b / U(sub infinity) is less than 10 a reduction in unsteady lift of up to 7 decibels (dB) occurs while an increase occurs for omega(sub r) is greater than 10 as the angle of attack is increased. The reduction in unsteady lift at low omega(sub r) with increasing angle of attack is a result that has never before been shown either experimentally or theoretically. The source of the reduction in lift spectral level appears to be closely related to the distortion of inflow turbulence based on analysis of surface pressure spanwise correlation length scales. Furthermore, while the distortion of the inflow appears to be critical in this experiment, this effect does not seem to be significant in larger integral scale (relative to the chord) flows based on the previous experimental work of McKeough suggesting the airfoils size relative to the inflow integral scale is critical in defining how the airfoil will respond under variation of angle of attack. A prediction scheme is developed that correctly accounts for the effects of distortion when the inflow integral scale is small relative to the airfoil chord. This scheme utilizes Rapid Distortion Theory to account for the distortion of the inflow with the distortion field modeled using a circular cylinder.

Mish, Patrick F.↗

Propagation of a Free Flame in a Turbulent Gas Stream

Effective flame speeds of free turbulent flames were measured by photographic, ionization-gap, and photomultiplier-tube methods, and were found to have a statistical distribution attributed to the nature of the turbulent field. The effective turbulent flame speeds for the free flame were less than those previously measured for flames stabilized on nozzle burners, Bunsen burners, and bluff bodies. The statistical spread of the effective turbulent flame speeds was markedly wider in the lean and rich fuel-air-ratio regions, which might be attributed to the greater sensitivity of laminar flame speed to flame temperature in those regions. Values calculated from the turbulent free-flame-speed analysis proposed by Tucker apparently form upper limits for the statistical spread of free-flame-speed data. Hot-wire anemometer measurements of the longitudinal velocity fluctuation intensity and longitudinal correlation coefficient were made and were employed in the comparison of data and in the theoretical calculation of turbulent flame speed.

Mickelsen, William R↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗

Concepts for a theoretical and experimental study of lifting rotor random loads and vibrations. Phase 5B: Analysis of gust alleviation methods and rotor dynamic stability

The effects of various gust alleviation methods on the random blade response in flapping were studied analytically, assuming a rigid rotor support. The analytical model assumes rigid flapping blades with elastic root restraints. Linearized equations which are approximately valid at low lift conditions were used. Because of the interblade coupling from the feedback devices, the method of multiblade generalized coordinates was most convenient and was extended to include coning, differential coning, and warping of the rotor. The numerical examples cover the dynamic stability characteristics as affected by feedback gains of three- to six-bladed rotors. The number of blades had large effects on stability limits and modal time functions at these limits. The random flapping response of the blades to atmospheric turbulence was determined at 1.6 rotor advance ratio using feedback gains below the stability limit. The most effective reduction of the flapping response per unit gain was achieved with a rotor coning angle feedback.

Hohenemser, K. H.↗

Wind-tunnel investigation of a full-scale general aviation airplane equipped with an advanced natural laminar flow wing

An investigation was conducted in the Langley 30- by 60-Foot Wind Tunnel to evaluate the performance, stability, and control characteristics of a full-scale general aviation airplane equipped with an advanced laminar flow wing. The study focused on the effects of natural laminar flow and advanced boundary layer transition on performance, stability, and control, and also on the effects of several wing leading edge modifications on the stall/departure resistance of the configuration. Data were measured over an angle-of-attack range from -6 to 40 deg and an angle-of-sideslip range from -6 to 20 deg. The Reynolds number was varied from 1.4 to 2.4 x 10 to the 6th power based on the mean aerodynamic chord. Additional measurements were made using hot-film and sublimating chemical techniques to determine the condition of the wing boundary layer, and wool tufts were used to study the wing stall characteristics. The investigation showed that large regions of natural laminar flow existed on the wing which would significantly enhance cruise performance. Also, because of the characteristics of the airfoil section, artificially tripping the wing boundary layer to a turbulent condition did not significantly effect the lift, stability, and control characteristics. The addition of a leading-edge droop arrangement was found to increase the stall angle of attack at the wingtips and, therefore, was considered to be effective in improving the stall/departure resistance of the configuration. Also the addition of the droop arrangement resulted in only minor increases in drag.

Murri, Daniel G.↗

Wave-turbulence dynamics in the stably stratified boundary layer

New data obtained at the Boulder Atmospheric Observatory are analyzed to obtain separation of wave, turbulence, and mean field necessary for a complete treatment of wave-turbulence interaction. The data were compared with a linear stability analysis of the background atmospheric state, showing good agreement between measured wave parameters (such as wavelength, period, and vector phase velocity) and the eigenvalues of the linear solution. The analysis of the budgets of wave heat flux and temperature variance revealed the essential role of wave-turbulence interaction in maintaining a large amplitude temperature wave and countergradient heat flux. A mechanism for the maintenance of turbulence by waves in strongly stratified boundary layers is described, which emphasizes that the time-mean Richardson number is an irrelevant parameter at such times.

Einaudi, F.↗

A turbulence module for the NPARC code

A turbulence module is developed for the 2D version of the NPARC code which is currently restricted to planar or axisymmetric flows without swirling. Four turbulence models have been built into the module: Baldwin-Lomax, Chien, Shih-Lumley and CMOTT models. The first is a mixing-length eddy-viscosity model which is mainly used for initialization of computational fields and the last three are the low Reynolds number two-equation models. Unlike chien's model, both the Shih-Lumley and CMOTT models do not involve the dimensionless wall distance y(sup +), an advantage for separated flow calculations. Contrary to the NPARC and most other compressible codes, the non-delta form of transport equations is used which leads to a simpler linearization and is more effective than using the delta form in ensuring the positiveness of the turbulent kinetic energy and its dissipation rate. To reduce numerical diffusion while maintaining necessary stability, a second-order accurate and bounded scheme is used for the convective terms of the turbulent transport equations. This scheme is implemented in a deferred correction manner so that the main coefficients of the resulting difference equations are always positive, thus making the numerical solutions process unconditionally stable. The system of equations is solved via a decoupled method and by the alternating direction TDMA of Thomas. The module can be easily linked to the NPARC code for turbulent flow calculations.

Zhu, J.↗

Turbulent Jet Flames Into a Vitiated Coflow

Examined is the vitiated coflow flame, an experimental condition that decouples the combustion processes of flows found in practical combustors from the associated recirculating fluid mechanics. The configuration consists of a 4.57 mm diameter fuel jet into a coaxial flow of hot combustion products from a lean premixed flame. The 210 mm diameter coflow isolates the jet flame from the cool ambient, providing a hot environment similar to the operating conditions of advanced combustors; this important high temperature element is lacking in the traditional laboratory experiments of jet flames into cool (room) air. A family of flows of increasing complexity is presented: 1) nonreacting flow, 2) all hydrogen flame (fuel jet and premixed coflow), and 3) set of methane flames. This sequence of experiments provides a convenient ordering of validation data for combustion models. Laser Raman-Rayleigh-LIF diagnostics at the Turbulent Diffusion Flame laboratory of Sandia National Laboratories produced instantaneous multiscalar point measurements. These results attest to the attractive features of the vitiated coflow burner and the well-defined boundary conditions provided by the coflow. The coflow is uniform and steady, isolating the jet flame from the laboratory air for a downstream distance ranging from z/d = 50-70. The statistical results show that differential diffusion effects in this highly turbulent flow are negligible. Complementing the comprehensive set of multiscalar measurements is a parametric study of lifted methane flames that was conducted to analyze flame sensitivity to jet and coflow velocity, as well as coflow temperature. The linear relationship found between the lift-off height and the jet velocity is consistent with previous experiments. New linear sensitivities were found correlating the lift-off height to coflow velocity and temperature. A blow-off study revealed that the methane flame blows off at a common coflow temperature (1260 K), regardless of coflow or jet velocity. An explanation for this phenomenon is that entrainment of ambient air at the high lift-off heights prevents autoignition. Analysis of the results suggests that flame stabilization occurs through a combination of flame propagation, autoignition, and localized extinction processes. Proposed is an expanded view of distributed reaction combustion based on analysis of the distributions of probe volume conditions at the stabilization region of the lifted hydrogen and methane flames. Turbulent eddies the size of the flame thickness mix fuel and hot coflow across the flame front, thereby enhancing the reaction zone with autoignition of reactants at elevated temperatures; this is the reverse effect of turbulent flames in ambient air, where intense turbulence in cool mixtures result in localized extinction. Each of the three processes (i.e., flame propagation, autoignition and localized extinction) contributes to flame stabilization in varying degrees, depending on flow conditions.

Holdeman, James D.↗

Stability of the Boundary Layer and the Spot

The similarity among turbulent spots observed in various transition experiments, and the rate in which they contaminate the surrounding laminar boundary layer is only cursory. The shape of the spot depends on the Reynolds number of the surrounding boundary layer and on the pressure gradient to which it and the surrounding laminar flow are exposed. The propagation speeds of the spot boundaries depend, in addition, on the location from which the spot originated and do not simply scale with the local free stream velocity. The understanding of the manner in which the turbulent manner in which the turbulent spot destabilizes the surrounding, vortical fluid is a key to the understanding of the transition process. We therefore turned to detailed observations near the spot boundaries in general and near the spanwise tip of the spot in particular.

Wygnanski, I.↗

A comprehensive model to determine the effects of temperature and species fluctuations on reaction rates in turbulent reaction flows

The use of silane (SiH4) as an effective ignitor and flame stabilizing pilot fuel is well documented. A reliable chemical kinetic mechanism for prediction of its behavior at the conditions encountered in the combustor of a SCRAMJET engine was calculated. The effects of hydrogen addition on hydrocarbon ignition and flame stabilization as a means for reduction of lengthy ignition delays and reaction times were studied. The ranges of applicability of chemical kinetic models of hydrogen-air combustors were also investigated. The CHARNAL computer code was applied to the turbulent reaction rate modeling.

Magnotti, F.↗

In Situ Detection of Strong Langmuir Turbulence Processes in Solar Type III Radio Bursts

The high time resolution observations obtained by the WAVES experiment of the STEREO spacecraft in solar type III radio bursts show that Langmuir waves often occur as intense localized wave packets. These wave packets are characterized by short durations of only a few ms and peak intensities, which well exceed the supersonic modulational instability (MI) thresholds. These timescales and peak intensities satisfy the criterion of the solitons collapsed to spatial scales of a few hundred Debye lengths. The spectra of these wave packets consist of primary spectral peaks corresponding to beam-resonant Langmuir waves, two or more sidebands corresponding to down-shifted and up-shifted daughter Langmuir waves, and low frequency enhancements below a few hundred Hz corresponding to daughter ion sound waves. The frequencies and wave numbers of these spectral components satisfy the resonance conditions of the modulational instability (MI). Moreover, the tricoherences, computed using trispectral analysis techniques show that these spectral components are coupled to each other with a high degree of coherency as expected of the MI type of four wave interactions. The high intensities, short scale lengths, sideband spectral structures and low frequency spectral enhancements and, high levels of tricoherences amongst the spectral components of these wave packets provide unambiguous evidence for the supersonic MI and related strong turbulence processes in type III radio bursts. The implication of these observations include: (1) the MI and related strong turbulence processes often occur in type III source regions, (2) the strong turbulence processes probably play very important roles in beam stabilization as well as conversion of Langmuir waves into escaping radiation at the fundamental and second harmonic of the electron plasma frequency, fpe, and (3) the Langmuir collapse probably follows the route of MI in type III radio bursts.

Langmuir↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗

Flight evaluation of highly augmented controls and electronic displays for precision approach and landing of powered-lift aircraft

Experiments were conducted on simulators and on the Quiet Short-Haul Research Aircraft to evaluate the effect of highly augmented control modes and electronic displays on the ability of pilots to execute precision approaches and landings on a short runway. It is found that the primary benefits of highly augmented flightpath and airspeed controls and electronic displays are realized when the pilot is required to execute precisely a complex transition and approach under instrument conditions and in the presence of a wide range of wind and turbulence conditions. A flightpath and airspeed command and stabilization system incorporating nonlinear, inverse system concepts produced fully satisfactory flightpath control throughout the aircraft's terminal operating envelope.

Franklin, J. A.↗

Three-Dimensional Navier-Stokes Simulation of Space Shuttle Main Propulsion 17-inch Disconnect Valves

A steady incompressible three-dimensional viscous flow analysis has been conducted for the Space Shuttle external tank/orbiter propellant feed line disconnect flapper valves with upstream elbows. The Navier-Stokes code, INS3D, is modified to handle interior obstacles and a simple turbulence model. The flow solver is tested for stability and convergence in the presence of interior flappers. An under-relaxation scheme has been incorporated to improve the solution stability. Important flow characteristics such as secondary flows, recirculation, vortex and wake regions, and separated flows are observed. Computed values for forces, moments, and pressure drop are in satisfactory agreement with water flow test data covering a maximum tube Reynolds number of 3.5 million. The predicted hydrodynamical stability of the flappers correlates well with the measurements.

Kandula, M.↗

Reliability enhancement of Navier-Stokes codes through convergence acceleration

Methods for enhancing the reliability of Navier-Stokes computer codes through improving convergence characteristics are presented. The improving of these characteristics decreases the likelihood of code unreliability and user interventions in a design environment. The problem referred to as a 'stiffness' in the governing equations for propulsion-related flowfields is investigated, particularly in regard to common sources of equation stiffness that lead to convergence degradation of CFD algorithms. Von Neumann stability theory is employed as a tool to study the convergence difficulties involved. Based on the stability results, improved algorithms are devised to ensure efficient convergence in different situations. A number of test cases are considered to confirm a correlation between stability theory and numerical convergence. The examples of turbulent and reacting flow are presented, and a generalized form of the preconditioning matrix is derived to handle these problems, i.e., the problems involving additional differential equations for describing the transport of turbulent kinetic energy, dissipation rate and chemical species. Algorithms for unsteady computations are considered. The extension of the preconditioning techniques and algorithms derived for Navier-Stokes computations to three-dimensional flow problems is discussed. New methods to accelerate the convergence of iterative schemes for the numerical integration of systems of partial differential equtions are developed, with a special emphasis on the acceleration of convergence on highly clustered grids.

Merkle, Charles L.↗

The Martian Surface Boundary Layer

The acquisition of meteorological data from the surface of Mars by the two Viking Landers and Mars Pathfinder make it possible to estimate atmospheric boundary layer parameters and surface properties at three different locations on the planet. Because the Martian atmosphere is so thin the majority of the solar radiance is converted to heat at the surface. The difference between surface and atmospheric temperature can also constraint surface albedo, thermal inertia, and infrared emissivity. The Mars Pathfinder Atmospheric Structure Instrument/Meteorological package (ASI/MET) was the most capable weather monitoring system ever sent to the surface of another planet to date. One of the prime objectives of the ASI/MET package is to characterize the surface boundary layer parameters, particularly the heat and momentum fluxes, scaling temperature and friction velocity, and estimate surface roughness. Other important boundary layer parameters, such as Richardson Number, Monin-Obukhov length, analysis of turbulence characteristics of wind and temperature, and atmospheric stability class can also be determined from these measurements.

Wilson, G. R.↗

Enhancements to Linear Stability-Based, CFD-integrated Transition Prediction for High-Speed Flows

Combining linear stability calculations with computational fluid dynamics (CFD) simulations has great potential for the automated modeling of high-speed flows, especially when adequate information about the configuration and the disturbance environment is available. However, a significant impediment to the applicability of this technique is the lack of an efficient method to calculate the crucial amplification ratio corresponding to the onset of transition in hypersonic flows. This ratio, also known as the "transition N-factor," is dependent upon the freestream disturbance environment as well as the surface properties of the test article. In response to the need for an engineering solution to predict the transition N-factor within conventional hypersonic wind tunnels, this paper presents a data-driven correlation that expands the existing correlations from straight circular cones with a narrow range of half angles to a broader array of axisymmetric configurations. Furthermore, when tested against a chosen dataset that was not used in its calibration, the suggested correlation shows good predictive accuracy with an RMS error of only 6.9%. Although similar accuracy may also be achieved via existing correlations based on similar datasets, predictions based on the proposed correlation have the advantage of not requiring an extensive amount of configuration-specific data. Practical applications often have access to the input parameters needed for this correlation, such as the freestream disturbance intensity, Mach number, and body-based slenderness Reynolds number. Additionally, this correlation outperforms the traditional assumption of a constant N-factor, particularly for configurations with blunted nose geometries. The development of this correlation is grounded in an extensive dataset encompassing conical models with body half-angles varying between 5 degrees and 16 degrees, Mach numbers ranging from 5 to 14, and nosetip-based Reynolds numbers approaching the transition reversal limit for blunt-nosed cones.

CFD↗

Enhancements to Linear Stability-Based, CFD-integrated Transition Prediction for High-Speed Flows

Combining linear stability calculations with computational fluid dynamics (CFD) simulations has great potential for the automated modeling of high-speed flows, especially when adequate information about the configuration and the disturbance environment is available. However, a significant impediment to the applicability of this technique is the lack of an efficient method to calculate the crucial amplification ratio corresponding to the onset of transition in hypersonic flows. This ratio, also known as the "transition N-factor," is dependent upon the freestream disturbance environment as well as the surface properties of the test article. In response to the need for an engineering solution to predict the transition N-factor within conventional hypersonic wind tunnels, this paper presents a data-driven correlation that expands the existing correlations from straight circular cones with a narrow range of half angles to a broader array of axisymmetric configurations. Furthermore, when tested against a chosen dataset that was not used in its calibration, the suggested correlation shows good predictive accuracy with an RMS error of only 6.9%. Although similar accuracy may also be achieved via existing correlations based on similar datasets, predictions based on the proposed correlation have the advantage of not requiring an extensive amount of configuration-specific data. Practical applications often have access to the input parameters needed for this correlation, such as the freestream disturbance intensity, Mach number, and body-based slenderness Reynolds number. Additionally, this correlation outperforms the traditional assumption of a constant N-factor, particularly for configurations with blunted nose geometries. The development of this correlation is grounded in an extensive dataset encompassing conical models with body half-angles varying between 5 degrees and 16 degrees, Mach numbers ranging from 5 to 14, and nosetip-based Reynolds numbers approaching the transition reversal limit for blunt-nosed cones.

CFD↗