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 577 records · Page 32

Computation of viscous flows over airfoils, including separation, with a coupling approach

Viscous incompressible flows over single or multiple airfoils, with or without separation, were computed using an inviscid flow calculation, with modified boundary conditions, and by a method providing calculation and coupling for boundary layers and wakes, within conditions of strong viscous interaction. The inviscid flow is calculated with a method of singularities, the numerics of which were improved by using both source and vortex distributions over profiles, associated with regularity conditions for the fictitious flows inside of the airfoils. The viscous calculation estimates the difference between viscous flow and inviscid interacting flow, with a direct or inverse integral method, laminar or turbulent, with or without reverse flow. The numerical method for coupling determines iteratively the boundary conditions for the inviscid flow. For attached viscous layers regions, an underrelaxation is locally calculated to insure stability. For separated or separating regions, a special semi-inverse algorithm is used. Comparisons with experiments are presented.

Leballeur, J. C.↗

Meteorological and aeronomical requirements for MST radar networks (keynote paper), part 1

Mesosphere - stratosphere - troposphere (MST) radar are phase coherent radars that measure the amplitude and Doppler shift of radio waves that are scattered back to the receiving antennas. For a monostatic system, the line-of-sight projection of the wind vector is obtained if one assumes that the atmospheric scatterers are being swept along with the wind velocity. The three-dimensional wind is then derived either by using multiple beams or by beam swinging. The turbulence intensity is derived either by measuring the backscattered power or by deriving the width of the autocorrelation function for the wind. Furthermore, some information on sharp changes in the atmospheric static stability (e.g., at the tropopause) can be obtained by looking for specular reflections. The discussion addresses the question of how these MST measurement capabilities can contribute to various meteorological and aeronomical research areas.

Geller, M. A.↗

Laminar-Turbulent Transition in High Mach-Number Flows

Accurate prediction of laminar-turbulent transition is critical for the design and characterization of high-speed flight systems. Presently, laminar-turbulent transition prediction in computational fluid dynamics (CFD) is accomplished using PDE-based models coupled to the governing equations of the flow system. Current industrial models suffer significant performance issues when analyzing high-speed systems because they were developed for lower speed flows. Current research focuses on multiple solutions to this problem. The coupling of physics-based stability analysis tools which produce quality transition predictions for high-speed flows to NASA’s OVERFLOW CFD solver is being pursued, as well as the implementation and improvement of emerging models developed specifically for high-speed flow.

High-Speed Flow↗

A 7.5-Year Dataset of SSM/I-Derived Surface Turbulent Fluxes Over Global Oceans

The surface turbulent fluxes of momentum, latent heat, and sensible heat over global oceans are essential to weather, climate and ocean problems. Wind stress is the major forcing for driving the oceanic circulation, while Evaporation is a key component of hydrological cycle and surface heat budget. We have produced a 7.5-year (July 1987-December 1994) dataset of daily, individual monthly-mean and climatological (1988-94) monthly-mean surface turbulent fluxes over the global oceans from measurements of the Special Sensor Microwave/Imager (SSM/I) on board the US Defense Meteorological Satellite Program F8, F10, and F11 satellites. It has a spatial resolution of 2.0x2.5 latitude-longitude. Daily turbulent fluxes are derived from daily data of SSM/I surface winds and specific humidity, National Centers for Environmental Prediction (NCEP) sea surface temperatures, and European Centre for Medium-Range Weather Forecasts (ECMWF) air-sea temperature differences, using a stability-dependent bulk scheme. The retrieved instantaneous surface air humidity (with a 25-km resolution) IS found to be generally accurate as compared to the collocated radiosonde observations over global oceans. The surface wind speed and specific humidity (latent heat flux) derived from the F10 SSM/I are found to be -encrally smaller (larger) than those retrieved from the F11 SSM/I. The F11 SSM/I appears to have slightly better retrieval accuracy for surface wind speed and humidity as compared to the F10 SSM/I. This difference may be due to the orbital drift of the F10 satellite. The daily wind stresses and latent heat fluxes retrieved from F10 and F11 SSM/Is show useful accuracy as verified against the research quality in si -neasurerrients (IMET buoy, RV Moana Wave, and RV Wecoma) in the western Pacific warm pool during the TOGA COARE Intensive observing period (November 1992-February 1993). The 1988-94 seasonal-mean turbulent fluxes and input variables derived from FS and F11 SSM/Is show reasonable patterns related to seasonal variations of atmospheric general circulation. This dataset of SSM/I-derived turbulent fluxes is useful for climate studies, forcing of ocean models, and validation of coupled ocean-atmosphere global models and can be accessed through the NASA/GSFC Distributed Active Archive Center.

Chou, Shu-Hsien↗

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.

forward-facing step↗

Monte Carlo turbulence simulation using rational approximations to von Karman spectra

Turbulence simulation is computationally much simpler using rational spectra, but turbulence falls off as f exp -5/3 in frequency ranges of interest to aircraft response and as predicted by von Karman's model. Rational approximations to von Karman spectra should satisfy three requirements: (1) the rational spectra should provide a good approximation to the von Karman spectra in the frequency range of interest; (2) for stability, the resulting rational transfer function should have all its poles in the left half-plane; and (3) at high frequencies, the rational spectra must fall off as an integer power of frequency, and since the -2 power is closest to the -5/3 power, the rational approximation should roll off as the -2 power at high frequencies. Rational approximations to von Karman spectra that satisfy these three criteria are presented, along with spectra from simulated turbulence. Agreement between the spectra of the simulated turbulence and von Karman spectra is excellent.

Campbell, C. W.↗

Stability of compressible Taylor-Couette flow

The objectives of this paper are to: (1) develop both analytical and numerical tools that can be used to predict the onset of instability and subsequently to simulate the transition process by which the originally laminar flow evolves into a turbulent flow; and (2) conduct the preliminary investigations with the purpose of understanding the mechanisms of the vortical structures of the compressible flow between tow concentric cylinders.

Kao, K.↗

Further results of natural laminar flow flight test experiments

Flight test experiments were conducted to measure the extent and nature of natural laminar flow on a smoothed test region of a swept-wing business jet wing. Surface hot film anemometry and sublimating chemicals were used for transition detection. Surface pressure distributions were measured using pressure belts. Engine noise was monitored by a microphone attached to the wing surface to study possible acoustic effects on stability of the laminar boundary layer. Side-slip conditions were flown to simulate changes in effective wing sweep. Flight instrumentation and ground data analysis techniques and a method for measuring intermittency of turbulence are described. Correlation was obtained between the hot film gage signals and chemicals for transition detection. Cross-flow vortices were observed for some flight conditions. Results of spectral and statistical analysis of the hot film signals for various flight test conditions are presented.

Wentz, W. H., Jr.↗

Satellite-Based Surface Heat Budgets and Sea Surface Temperature Tendency in the Tropical Eastern Indian and Western Pacific Oceans for the 1997/98 El Nino and 1998/99 La Nina

The 1997/98 is a strong El Nino warm event, while the 1998/99 is a moderate La Nina cold event. We have investigated surface heat budgets and sea surface temperature (SST) tendency for these two events in the tropical western Pacific and eastern Indian Oceans using satellite-retrieved surface radiative and turbulent fluxes. The radiative fluxes are taken from the Goddard Satellite-retrieved Surface Radiation Budget (GSSRB), derived from radiance measurements of the Japanese Geostationary Meteorological Satellite 5. The GSSRB covers the domain 40 deg S - 4 deg N, 90 deg E-17 deg W and a period from October 1997 to December 2000. The spatial resolution is 0.5 deg x 0.5 deg lat-long and the temporal resolution is 1 day. The turbulent fluxes are taken from Version 2 of the Goddard Satellite-based Surface Turbulent Fluxes (GSSTF-2). The GSSTF-2 has a spatial resolution of 1 deg x 1 deg lat-long over global Oceans and a temporal resolution of 1 day covering the period July 1987-December 2000. Daily turbulent fluxes are derived from the S S M (Special Sensor Microwave/Imager) surface wind and surface air humidity, and the SST and 2-m air temperature of the NCEP/NCAR reanalysis, using a stability-dependent bulk flux algorithm. The changes of surface heat budgets, SST and tendency, cloudiness, wind speed, and zonal wind stress of the 1997/98 El Nino relative to the1998/99 La Nina for the northern winter and spring seasons are analyzed. The relative changes of surface heat budgets and SST tendency of the two events are quite different between the tropical eastern Indian and western Pacific Oceans. For the tropical western Pacific, reduced solar heating (more clouds) is generally associated with decreased evaporative cooling (weaker winds), and vise versa. The changes in evaporative cooling over-compensate that of solar heating and dominate the spatial variability of the changes in net surface heating. Both solar heating and evaporative cooling offset each other to reduce interannual variability of net surface heating. In addition, the area of increased SST tendency is larger than that of increased net surface heating, due to less solar radiation penetration through the bottom of deeper ocean mixed layer (stronger winds). For the tropical eastern Indian Ocean, enhanced solar heating (less clouds) is generally associated with reduced evaporative cooling (weaker winds). Both solar heating and evaporative cooling reinforce each other to increase interannual variability of net surface heating. In addition, the area of increased SST tendency is smaller than that of increased net surface heating in the southern domain. The relative changes in wind and zonal wind stress indicate more solar radiation penetration through the ocean mixed layer and more northward heat transport by Ocean current from the south to the north Indian Ocean for the El Nino than for the La Nina.

Chou, Shu-Hsien↗

Studying Turbulence Using Numerical Simulation Databases - X Proceedings of the 2004 Summer Program

This Proceedings volume contains 32 papers that span a wide range of topics that reflect the ubiquity of turbulence. The papers have been divided into six groups: 1) Solar Simulations; 2) Magnetohydrodynamics (MHD); 3) Large Eddy Simulation (LES) and Numerical Simulations; 4) Reynolds Averaged Navier Stokes (RANS) Modeling and Simulations; 5) Stability and Acoustics; 6) Combustion and Multi-Phase Flow.

Moin, Parviz↗

A computational examination of directional stability for smooth and chined forebodies at high-alpha

Computational Fluid Dynamics (CFD) has been used to study aircraft forebody flowfields at low-speed, angle-of-attack conditions with sideslip. The purpose is to define forebody geometries which provide good directional stability characteristics under these conditions. The flows over the experimentally investigated F-5A forebody and chine type configuration, previously computed by the authors, were recomputed with better grid topology and resolution. The results were obtained using a modified version of CFL3D (developed at NASA Langley) to solve either the Euler equations or the Reynolds equations employing the Baldwin-Lomax turbulence model with the Degani-Schiff modification to account for massive crossflow separation. Based on the results, it is concluded that current CFD methods can be used to investigate the aerodynamic characteristics of forebodies to achieve desirable high angle-of-attack characteristics. An analytically defined generic forebody model is described, and a parametric study of various forebody shapes was then conducted to determine which shapes promote a positive contribution to directional stability at high angle-of-attack. An unconventional approach for presenting the results is used to illustrate how the positive contribution arises. Based on the results of this initial parametric study, some guidelines for aerodynamic design to promote positive directional stability are presented.

Ravi, Ramakrishnan↗

Transition in Turbine Flows

We have further developed our capabilities to analyze transition in turbine boundary layers from first principles by integrating the nonlinear parabolized stability equations (PSE) with improved initial and boundary conditions. With modified iteration schemes, we are able to proceed deeper into the transition region where skin friction coefficient and heat transfer coefficient significantly increase. Initial and boundary conditions at elevated turbulence levels can be derived by receptivity analysis. Test runs for ERCOFTAC test case T3A at 2.4\% turbulence level provide results in good agreement with the experimental data. The sharper minimum of the skin coefficient also shown by DNS results is likely due to the missing intermittency. The method has been applied to various experimentally studied turbine blades (UTRC, VKI, Zierke, Langston, Hippensteele, and others). The PSE results, though physically reasonable, do not agree as well as expected with the experimental findings. We have, therefore, performed an extensive search for the reasons of the seemingly systematic deviations. A first source of uncertainty has been found in the often insufficient documentation of the experiments (e.g. on blockage by end-wall boundary layers). However, variation of the relevant parameters does not lead to more satisfactory agreement. A second reason has been found in the "standard procedure" which considers a 2D flow at midspan and uses a panel code and subsequent boundary-layer code to obtain the laminar basic flow for the transition analysis. Comparison with the pressure distribution obtained with a 3D design code (RVC3D) shows significant three-dimensionality of the flow (e.g. in the UTRC experiments). The spanwise variation has been neglected in our original PSE code. To overcome this problem, we have developed the PSE/3D for fully 3D boundary layers to account for streamwise and spanwise variations. Since the design code does not provide the boundary-layer flow with sufficient resolution, we have generated the Euler solution and employed a 3D boundary-layer code to obtain the viscous basic flow. Although only the linear stability level of PSE/3D has been implemented so far, the discrepancies with the experiments change but do not disappear. We still find deviations between the computed and experimental variations of C(sub f), and St along the blade for laminar flow. The main reason can be seen by comparing the solution of the boundary-layer code with the viscous results of the design code. The conventional boundary-layer solution exhibits an asymptotic behavior appropriate in external aerodynamics but does not match the steep gradients of the inviscid flow through the passage and consequently provides biased results for C(sub f), and St. An attempt is currently being made to correct this deficiency. Before attempting to perform the transition analysis for the viscous flow provided by the design code, we have analyzed the implementation and "best possible" results. Code and results exhibit flaws that may negatively affect the design and are intolerable for transition analysis. Therefore, we have decided to develop a new code to obtain a reliable basis for stability and transition studies. We expect to report improved results by the time of the meeting.

Herbert, Thorwald↗

Hypersonic Boundary-Layer Instabilities over Ogive-Cylinder Models

Computational investigations of an ogive-cylinder geometry with varying nosetips at zero degrees angle of attack are presented. The model geometry and conditions are selected to match experiments conducted in the Air Force Research Laboratory (AFRL) Mach 6 Ludwieg Tube. Five nosetips of interest were selected for the computational studies herein: two sharp ogives, two blunt ogives, and one hemispherical nosetip. Computations are performed at a freestream Reynolds number of 7.01 × 10 6 m -1 . Each sharp and blunt tip ogives had a 14 and 28 degree version for the tip angles. A cylindrical section follows the nosetip, resulting in a meter long model such as in the experiments. The laminar flow solutions are analyzed. The boundary-layer-edge properties and the velocity and temperature profiles are compared across streamwise locations aft of the ogive-cylinder junction. The blunt nosetips induce an entropy layer that envelopes the boundary-layer profiles. Modal stability analysis identifies most amplified frequencies corresponding to Mack’s second modes that agree with experimental results for the sharp tips. Similar to the experimental measurements based on wall-mounted pressure sensors, no unstable modes are found for the blunt models. Nonmodal analysis revealed a broadband set of disturbances present for the blunter tips, in agreement with experimental observations. Flow perturbation contours of most amplified planar and oblique disturbances are shown to qualitatively match wind tunnel Schlieren images, with a switch from rope-like to elongated structures, i.e., from high frequency Mack’s second modes to low frequency Mack’s first modes, as the nosetip angle is increased for the sharp tip, and from boundary-layer to entropy-layer disturbances as the bluntness is increased.

Boundary layer transition↗

Optimization and simulation of flight control laws under parameter uncertainty and external disturbances

Several tasks pertinent to flight control in parameter uncertainty and wind-gust loading were successfully completed. Identification algorithms for extracting stability and control derivatives from flight data taking gust loading into account were developed. They were verified by simulation and evaluated throughly on actual flight data taken on a Lockheed Jet Star flying in turbulence. In particular the need for automatically generated dither-like inputs was studied. Criteria for performance evaluation using stochastic models were developed for gust alleviation as well as handling quantities. Algorithms for assessing degradation in performance due to parameter uncertainty were developed and evaluated using flight test data.

Source record↗

Laboratory Study of Magnetorotational Instability and Hydrodynamic Stability at Large Reynolds Numbers

Two plausible mechanisms have been proposed to explain rapid angular momentum transport during accretion processes in astrophysical disks: nonlinear hydrodynamic instabilities and magnetorotational instability (MRI). A laboratory experiment in a short Taylor-Couette flow geometry has been constructed in Princeton to study both mechanisms, with novel features for better controls of the boundary-driven secondary flows (Ekman circulation). Initial results on hydrodynamic stability have shown negligible angular momentum transport in Keplerian-like flows with Reynolds numbers approaching one million, casting strong doubt on the viability of nonlinear hydrodynamic instability as a source for accretion disk turbulence.

Ji, H.↗

Effects of confinement on partially premixed flames

Partially premixed combustion is an intermediate regime between the limiting cases of premixed and nonpremixed combustion. Although combustion problems are generally approached from one of these two limiting cases, there are many practical situations where flames cannot be considered as purely premixed or nonpremixed, and thus the partially premixed approach must be used. In partially premixed combustion, mechanisms from the premixed and nonpremixed regimes can coexist, and as a result some interesting new phenomena can arise. One such phenomenon is the flame stabilization in laminar mixing layers by triple flames. One of the first observations of triple flames was made by Phillips (1965), who investigated a triple flame propagating in a methane mixing layer. Kioni et al. (1993) also examined triple flames both experimentally and numerically. There have also been numerous analytical studies on the shape and propagation of triple flames under various assumptions by Dold (1989), Dold et al. (1991), and Hartley and Dold (1991). In terms of modeling, Muller et al. (1994) have combined the flamelet formulations for premixed and nonpremixed combustion in order to treat lifted diffusion flames. One common feature in the analytical and numerical studies mentioned above is the assumption of zero heat release, which is necessary to make the problem tractable. The effect of heat release on triple flames was investigated by Ruetsch et al. (1995), where for the unconfined case, flame speeds larger than their premixed counterparts were found. One of the most important practical situations in which these conditions arise is in lifted turbulent jet diffusion flames. At a critical velocity the burning zone of a fuel jet lifts off from the nozzle, moves to increasing distances as the jet velocity increases, and finally blows off. The mechanisms that control these phenomena, i.e. that determine the stability of these flames, are still not understood. In addition to regions where diffusion flame stabilization takes place, partially premixed conditions also exist during the ignition process in nonpremixed systems. Numerical simulations by Reveillon et al. (1994) of the ignition process in a weakly stirred mixture of fuel and oxidizer show that triple flames propagate along lines of stoichiometric mixture fraction throughout the fluid. In addition, Peters (1994) notes that NO(x) emissions are likely to be large in such transient cases, and therefore an understanding of triple flames can provide information concerning pollutant formation. This study extends the work previously done and examines the effects of lateral confinement on partially premixed flames. Once again, we study both the flame structure and propagation.

Ruetsch, G. R.↗

Simultaneous Water Vapor and Dry Air Optical Path Length Measurements and Compensation with the Large Binocular Telescope Interferometer

The Large Binocular Telescope Interferometer uses a near-infrared camera to measure the optical path length variations between the two AO-corrected apertures and provide high-angular resolution observations for all its science channels (1.5-13 microns). There is however a wavelength dependent component to the atmospheric turbulence, which can introduce optical path length errors when observing at a wavelength different from that of the fringe sensing camera. Water vapor in particular is highly dispersive and its effect must be taken into account for high-precision infrared interferometric observations as described previously for VLTI/MIDI or the Keck Interferometer Nuller. In this paper, we describe the new sensing approach that has been developed at the LBT to measure and monitor the optical path length fluctuations due to dry air and water vapor separately. After reviewing the current performance of the system for dry air seeing compensation, we present simultaneous H-, K-, and N-band observations that illustrate the feasibility of our feed forward approach to stabilize the path length fluctuations seen by the LBTI nuller uses a near-infrared camera to measure the optical path length variations between the two AO-corrected apertures and provide high-angular resolution observations for all its science channels (1.5-13 microns). There is however a wavelength dependent component to the atmospheric turbulence, which can introduce optical path length errors when observing at a wavelength different from that of the fringe sensing camera. Water vapor in particular is highly dispersive and its effect must be taken into account for high-precision infrared interferometric observations as described previously for VLTI MIDI or the Keck Interferometer Nuller. In this paper, we describe the new sensing approach that has been developed at the LBT to measure and monitor the optical path length fluctuations due to dry air and water vapor separately. After reviewing the current performance of the system for dry air seeing compensation, we present simultaneous H-, K-, and N-band observations that illustrate the feasibility of our feed forward approach to stabilize the path length fluctuations seen by the LBTI nuller.

which can introduce optical path length errors whe↗

Streak Instability Analysis for BOLT Configuration

The BOLT configuration designed for an upcoming flight experiment exemplifies the complexities of hypersonic boundary layer transition in the presence of multiple and potentially interacting instability mechanisms. A numerical study of the boundary layer instabilities over this configuration has been performed with the goal of applying state-of-the-art tools in multidimensional stability analysis that complement the other ongoing research related to the BOLT flight campaign. Specifically, the paper outlines modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane of a 1/3 rd scale BOLT configuration that was tested in the low-disturbance hypersonic wind tunnel at Purdue University and both quiet and conventional tunnels at Texas A&M University (TAMU). The basic state computations agree well with the heat flux measurements in the Purdue quiet tunnel as well as with the off-body velocity field described in a previous DNS study at the University of Minnesota. The encouraging agreement between the predicted frequency spectra of the streak instability in the region of increased boundary layer thickness and the surface pressure spectra measured during the TAMU experiments suggests that these streak instabilities may play an important role during the transition process near the symmetry plane. The modest value of the predicted peak N-factor (N max ≈ 5) for the streak instabilities provides a potential explanation for why transition is not observed under quiet flow conditions, but may occur in a noisy wind tunnel. Similar instability analyses are performed for the flight configuration at two selected flow conditions from the ascent portion of the planned flight trajectory. The predicted values of the peak N-factor suggest that transition onset in the vicinity of the center plane is possible at the ascent time of t = 26.98 seconds (N max ≈ 12.5 at X/L = 0.77), but the streak instabilities weaken sufficiently at the later time of t = 28.88 seconds (N max ≈ 11.5 at X/L ≈ 1.00) and are less likely to cause transition at that time.

Laminar turbulent transition↗