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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Flow visualizations of perpendicular blade vortex interactions

Helium bubble flow visualizations have been performed to study perpendicular interaction of a turbulent trailing vortex and a rectangular wing in the Virginia Tech Stability Tunnel. Many combinations of vortex strength, vortex-blade separation (Z(sub s)) and blade angle of attack were studied. Photographs of representative cases are presented. A range of phenomena were observed. For Z(sub s) greater than a few percent chord the vortex is deflected as it passes the blade under the influence of the local streamline curvature and its image in the blade. Initially the interaction appears to have no influence on the core. Downstream, however, the vortex core begins to diffuse and grow, presumably as a consequence of its interaction with the blade wake. The magnitude of these effects increases with reduction in Z(sub s). For Z(sub s) near zero the form of the interaction changes and becomes dependent on the vortex strength. For lower strengths the vortex appears to split into two filaments on the leading edge of the blade, one passing on the pressure and one passing on the suction side. At higher strengths the vortex bursts in the vicinity of the leading edge. In either case the core of its remnants then rapidly diffuse with distance downstream. Increase in Reynolds number did not qualitatively affect the flow apart from decreasing the amplitude of the small low-frequency wandering motions of the vortex. Changes in wing tip geometry and boundary layer trip had very little effect.

Rife, Michael C.↗
NASA NTRS1948

Theoretical Analysis of Stationary Potential Flows and Boundary Layers at High Speed

The present report consists of two parts. The first part deals with the two-dimensional stationary flow in the presence of local supersonic zones. A numerical method of integration of the equation of gas dynamics is developed. Proceeding from solutions at great distance from the body the flow pattern is calculated step by step. Accordingly the related body form is obtained at the end of the calculation. The second part treats the relationship between the displacement thickness of laminar and turbulent boundary layers and the pressure distribution at high speeds. The stability of the boundary layer is investigated, resulting in basic differences in the behavior of subsonic and supersonic flows. Lastly, the decisive importance of the boundary layer for the pressure distribution, particularly for thin profiles, is demonstrated.

Oswaititsch, K.↗
NASA NTRS2007

Analysis of Crossflow Transition Flight Experiment aboard the Pegasus Launch Vehicle

The Pegasus wing-glove flight experiment was designed to provide crossflow transition data at high Mach numbers, specifically to help validate stability based predictions for transition onset in a flight environment. This paper provides an analysis of the flight experiment, with emphasis on computational results for crossflow disturbances and the correlation of disturbance growth factors with in-flight transition locations via the e(sup N) method. Implications of the flight data for attachment line stability are also examined. Analysis of the thermocouple data reveals that transition (from turbulent to laminar flow) was first detected during the ascending flight of the rocket when the free stream Mach number exceeded about 4. Therefore, computations have been performed for flight Mach numbers of 4.13, 4.35, 4.56 and 4.99. Due to continually decreasing unit Reynolds number at higher altitudes, the entire wing-glove boundary layer became laminar at the highest flight Mach number computed above. In contrast, the boundary layer flow over the inboard tile region remained transitional up to and somewhat beyond the time of laminarization over the instrumented glove region. Linear stability predictions confirmed that the tile boundary layer is indeed more unstable to crossflow disturbances than the much colder stainless steel glove boundary layer. The transition locations based on thermocouple data from both the glove and the tile regions are found to correlate with stationary-crossflow N-factors within the range of 7 to 12.4 and with traveling mode N-factors between 7.6 and 14.1. Data from the thermocouples and hot film sensors indicates that transition from turbulent to laminar flow (i.e., laminarization) at a fixed point over the glove is generally completed within a flight time interval of 3 seconds. However, the times at which transition begins and ends as inferred from the hot film sensors are found to differ by about 2 seconds from the corresponding estimates based on the thermocouple data.

Malik, Mujeeb R.↗
NASA NTRS1984

Applications of spectral methods to turbulent magnetofluids in space and fusion research

Recent and potential applications of spectral method computation to incompressible, dissipative magnetohydrodynamics are surveyed. Linear stability problems for one dimensional, quasi-equilibria are approachable through a close analogue of the Orr-Sommerfeld equation. It is likely that for Reynolds-like numbers above certain as-yet-undetermined thresholds, all magnetofluids are turbulent. Four recent effects in MHD turbulence are remarked upon, as they have displayed themselves in spectral method computations: (1) inverse cascades; (2) small-scale intermittent dissipative structures; (3) selective decays of ideal global invariants relative to each other; and (4) anisotropy induced by a mean dc magnetic field. Two more conjectured applications are suggested. All the turbulent processes discussed are sometimes involved in current carrying confined fusion magnetoplasmas and in space plasmas.

Montgomery, D.↗
NASA NTRS2021

Coupling of the FUN3D Unstructured Flow Solver and the LASTRAC Stability Code to Model Transition

We develop an iterative automated method to predict transition locations in boundary-layer flows by using the FUN3D solver to perform flow simulations and the LASTRAC code for linear stability computations. The coupling of FUN3D and LASTRAC allows for a robust physics-based approach to model boundary-layer transition by analyzing the growth of different instability waves and then using that information to iteratively update the resulting transition location. There is no user involvement during the iterative computations. We apply this automated method to subsonic flow over a flat plate with a sharp leading edge. The final solution has regions of laminar and turbulent flow with a transition onset location that agrees with experiments and stability-based correlations. This iterative automated method is also applied to an NLF(1)-0416 airfoil for conditions with and without a separation bubble. Along with predicting transition locations, we compare the streamwise distributions of surface-pressure and skin-friction coefficients to a transport-equation-based model. We consider a 6:1 prolate spheroid at three angles of attack, namely, five, ten, and fifteen degrees, where mixed-mode transition occurs due to both Tollmien-Schlichting and crossflow instabilities. The skin-friction contours and transition fronts at every angle of attack from our iterative automated method show good agreement with past experimental and computational results in the literature for the 6:1 prolate spheroid.

Transition↗
NASA NTRS1994

Postinstability Models in Dynamics

This paper is devoted to the concept of instability in dynamical systems with the main emphasis on orbital, Hadamard, and Reynolds instabilities. It demonstrates that the requirement about differentiability in dynamics in some cases is not consistent with the physical nature of motions, and may lead to unrealistic solutions. Special attention is paid to the fact that instability is not an invariant of motion: it depends upon frames of reference, the metric of configuration space, and classes of functions selected for mathematical models of physical phenomena. This leads to the possibility of elimination of certain types of instabilities (in particular, those which lead to chaos and turbulence) by enlarging the class of functions using the Reynolds-type transformation in combination with the stabilization principle: the additional terms (the so-called Reynolds stresses) are found from the conditions that they suppress the original instability. Based upon these ideas, a new approach to chaos and turbulence as well as a new mathematical formalism for nonlinear dynamics are discussed.

Zak, Michail↗
NASA NTRS2017

Instability Wave-Streak Interactions in a Supersonic Boundary Layer

Optimal initial conditions for transient growth in a two-dimensional boundary layer flow correspond to stationary, counter-rotating vortices that subsequently develop into streamwise elongated streaks, which are characterized by an alternating pattern of low and high streamwise velocity. For incompressible flows, previous studies have shown that boundary layer modulation due to streaks below a threshold amplitude level can stabilize the Tollmien-Schlichting instability waves, resulting in a delay in the onset of laminar-turbulent transition. In the supersonic regime, the linearly, most-amplified waves become three-dimensional, corresponding to oblique, first-mode waves. This change in the character of dominant instabilities leads to an important change in the transition process, which is now dominated by oblique breakdown via nonlinear interactions between pairs of first-mode waves that propagate at equal but opposite angles with respect to the free stream. Because the oblique breakdown process is characterized by a rapid amplification of stationary streamwise streaks, artificial excitation of such streaks may be expected to promote transition in a supersonic boundary layer. Indeed, suppression of those streaks has been shown to delay the onset of transition in prior literature. Consistent with those findings, the present study shows that optimally growing stationary streaks indeed destabilize the first-mode waves, but only when the spanwise wavelength of the instability waves is equal to or smaller than twice the streak spacing. Transition in a benign disturbance environment typically involves first-mode waves with significantly longer spanwise wavelengths, and hence, these waves are stabilized by the optimal growth streaks. Thus, as long as the amplification factors for the destabilized, short wavelength instability waves remain below the threshold level for transition, a significant net stabilization is achieved, yielding a transition delay that is comparable to the length of the laminar region in the uncontrolled case.

Boundary layer transition↗
NASA NTRS1965

Flight Tests Related to Jet-Transport Upset and Turbulent-Air Penetration

A flight program, utilizing a Convair 880 and a Boeing 720 airplane, was conducted in conjunction with wind-tunnel and simulator programs to study problems related to jet-transport upsets and operation in a turbulent environment. During the handling-qualities portion of the program the basic static stability of the airplanes was considered to be satisfactory and the lateral-directional damping was considered to be marginal without damper augmentation. An evaluation of the longitudinal control system indicated that this system can become marginal in effectiveness in the high Mach number and high dynamic-pressure range of the flight envelope. From the upset and recovery phase of the program it was apparent that retrimming the stabilizer and spoiler deployment were valuable tools in effecting a positive recovery; however, if these devices are to be used safely, it appears that a suitable g-meter should be provided in the cockpit because the high control forces in recovery tend to reduce the pilot's sensitivity to the actual acceleration loads. During the turbulence penetrations the pilot noted that the measured vibrations of 4 to 6 cps in the cockpit considerably disrupted their normal scan pattern and suggested that an improvement should be made in the seat cushion and restraint system. Also it was observed that the indicator needles on the flight instruments were quite stable in the turbulent environment.

Andrews, William H.↗
NASA NTRS1983

Stability and Interaction of Coherent Structure in Supersonic Reactive Wakes

A theoretical formulation and analysis is presented for a study of the stability and interaction of coherent structure in reacting free shear layers. The physical problem under investigation is a premixed hydrogen-oxygen reacting shear layer in the wake of a thin flat plate. The coherent structure is modeled as a periodic disturbance and its stability is determined by the application of linearized hydrodynamic stability theory which results in a generalized eigenvalue problem for reactive flows. Detailed stability analysis of the reactive wake for neutral, symmetrical and antisymmetrical disturbance is presented. Reactive stability criteria is shown to be quite different from classical non-reactive stability. The interaction between the mean flow, coherent structure and fine-scale turbulence is theoretically formulated using the von-Kaman integral technique. Both time-averaging and conditional phase averaging are necessary to separate the three types of motion. The resulting integro-differential equations can then be solved subject to initial conditions with appropriate shape functions. In the laminar flow transition region of interest, the spatial interaction between the mean motion and coherent structure is calculated for both non-reactive and reactive conditions and compared with experimental data wherever available. The fine-scale turbulent motion determined by the application of integral analysis to the fluctuation equations. Since at present this turbulence model is still untested, turbulence is modeled in the interaction problem by a simple algebraic eddy viscosity model. The applicability of the integral turbulence model formulated here is studied parametrically by integrating these equations for the simple case of self-similar mean motion with assumed shape functions. The effect of the motion of the coherent structure is studied and very good agreement is obtained with previous experimental and theoretical works for non-reactive flow. For the reactive case, lack of experimental data made direct comparison difficult. It was determined that the growth rate of the disturbance amplitude is lower for reactive case. The results indicate that the reactive flow stability is in qualitative agreement with experimental observation.

Menon, Suresh↗
NASA NTRS1999

Ocean Turbulence: Two-Point Closure Model Momentum, Heat and Salt Vertical Diffusivities in the Presence of Shear - Paper 3

In papers 1 and 2 we have presented the results of the most updated 1-point closure model for the turbulent vertical diffusivities of momentum, heat and salt, K(sub m,h,s). In this paper, we derive the analytic expressions for K(sub m,h,s) using a new 2-point closure model that has recently been developed and successfully tested against some approx. 80 turbulence statistics for different flows. The new model has no free parameters. The expressions for K(sub m, h. s) are analytical functions of two stability parameters: the Turner number R(sub rho) (salinity gradient/temperature gradient) and the Richardson number R(sub i) (temperature gradient/shear). The turbulent kinetic energy K and its rate of dissipation may be taken local or non-local (K-epsilon model). Contrary to all previous models that to describe turbulent mixing below the mixed layer (ML) have adopted three adjustable "background diffusivities" for momentum. heat and salt, we propose a model that avoids such adjustable diffusivities. We assume that below the ML, K(sub m,h,s) have the same functional dependence on R(sub i) and R(sub rho) derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking measured by Gargett et al. (1981). The procedure frees the model from adjustable background diffusivities and indeed we use the same model throughout the entire vertical extent of the ocean. Using the new K(sub m,h, s), we run an O-GCM and present a variety of results that we compare with Levitus and the KPP model. Since the traditional 1-point (used in papers 1 and 2) and the new 2-point closure models used here represent different modeling philosophies and procedures, testing them in an O-GCM is indispensable. The basic motivation is to show that the new 2-point closure model gives results that are overall superior to the 1-point closure in spite of the fact that the latter rely on several adjustable parameters while the new 2-point closure has none. After the extensive comparisons presented in papers 1 and 2, we conclude that the new model presented here is overall superior for it not only is parameter free but also 2 because is part of a more general turbulence model that has been previously successfully tested on a wide variety of other types of turbulent flows.

Canuto, V. M.↗
NASA NTRS2007

The Many Worlds of Transition Research.

The transition from laminar to turbulent flow in a boundary layer is a complex phenomenon that may take different routes, each involving distinct stages governed by different, often not-yet unraveled dynamical principles. There are, surprisingly, questions concerning virtually every stage in the process, beginning with receptivity to external disturbances, the linear stability of spatially developing flows, different possible nonlinear end games, the formation and propagation of turbulent spots and the emergence of fully developed turbulent flow. There seems no doubt that the flow has to be seen as a forced, nonlinear spatio-temporal system, but the system is so complex that to extract simple insights is still very difficult.

Narasimha, R.↗
DOE OSTI2022

Large-eddy simulation of non-vaporizing sprays using the spectral-element method

Predictive simulations of high-pressure sprays require accurate representation of the turbulent gaseous flow field generated by liquid jet. Typically, the accuracy that can be obtained with low-order numerical methods (e.g. finite volume, finite element) is limited by stability issues in fine grids and the order of convergence of the method. In this work, we resolve the turbulent flow field in an Eulerian manner using the high-order spectral element method, coupled with a Lagrangian parcels approach to model the atomizing liquid jet. Large eddy simulations of single-hole sprays under non-evaporative conditions were conducted and compared against experimental data from Margot et al. (2008) and Spray A data from the Engine Combustion Network. The sensitivity of liquid penetration and droplet sizes to different breakup model parameters was studied. The effect of different numerical parameters, such as polynomial order of the solution (grid resolution), on liquid penetration was also analyzed. The method achieved grid-independent results using p-refinement, achieving finer resolution (by a factor of x 1.7 - x 3.5) in the gas-phase solution than in state-of-the-art simulations using the finite-volume method. Here the results showed good agreement with experimental data, demonstrating the ability of the current method to accurately capture liquid penetration and the shape of the spray.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Pseudospectral simulation of compressible turbulence using logarithmic variables

The direct numerical simulation of dissipative, highly compressible turbulent flow is performed using a pseudospectral Fourier technique. The governing equations are cast in a form where the important physical variables are the fluid velocity and the natural logarithms of the fluid density and temperature. Bulk viscosity is utilized to model polyatomic gases more accurately and to ensure numerical stability in the presence of strong shocks. Numerical examples include three-dimensional supersonic homogeneous turbulence and two-dimensional shock-turbulence interactions.

Shebalin, John V.↗
NASA NTRS1999

Ocean Turbulence: One-Point Closure Model Momentum, Heat and Salt Vertical Diffusivities in the Presence of Shear - Paper 2

We develop and test a 1-point closure turbulence model with the following features: 1) we include the salinity field and derive the expression for the vertical turbulent diffusivities of momentum K(sub m) , heat K(sub h) and salt K(sub s) as a function of two stability parameters: the Richardson number R(sub i) (stratification vs. shear) and the Turner number R(sub rho) (salinity gradient vs. temperature gradient). 2) to describe turbulent mixing below the mixed layer (ML), all previous models have adopted three adjustable "background diffusivities" for momentum, heat and salt. We propose a model that avoids such adjustable diffusivities. We assume that below the ML, the three diffusivities have the same functional dependence on R( sub i) and R(sub rho) as derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking.measured by Gargett et al. The procedure frees the model from adjustable background diffusivities and indeed we employ the same model throughout the entire vertical extent of the ocean. 3) in the local model, the turbulent diffusivities K(sub m,h,s) are given as analytical functions of R(sub i) and R(sub rho). 5) the model is used in an O-GCM and several results are presented to exhibit the effect of double diffusion processes. 6) the code is available upon request.

Canuto, V. M.↗
NASA NTRS1980

Experimental analysis of the boundary layer transition with zero and positive pressure gradient

The influence of a positive pressure gradient on the boundary layer transition is studied. The mean velocity and turbulence profiles of four cases are examined. As the intensity of the pressure gradient is increased, the Reynolds number of the transition onset and the length of the transition region are reduced. The Tollmein-Schlichting waves disturb the laminar regime; the amplification of these waves is in good agreement with the stability theory. The three dimensional deformation of the waves leads finally to the appearance of turbulence. In the case of zero pressure gradient, the properties of the turbulent spots are studied by conditional sampling of the hot-wire signal; in the case of positive pressure gradient, the turbulence appears in a progressive manner and the turbulent spots are much more difficult to characterize.

Arnal, D.↗
NASA NTRS1993

Effects of free-stream turbulence on boundary-layer transition

Traditional methods for the prediction of the transition location rest either on lump parameters like the momentum thickness or on N factors obtained from the approximate linear stability characteristics of the flow. These methods account for the effects of the disturbance environment by empirical correlation with observed transition locations. We have developed a highly efficient and accurate approach to stability analysis and transition simulation based on parabolized stability equations (PSE). Here we present a first application of this approach to realistic problems of practical interest. We study the effects of the turbulence level on transition and heat transfer in flows over flat plates, curved plates, and the stator blade of a gas-turbine model. We discuss the formulation of an input model for the PSE analysis in the light of the present knowledge of receptivity mechanisms. The input model is used to compute transition locations for various configurations and flow conditions. The results compare favorably with experimental data.

Herbert, TH.↗
DOE OSTI2020

Reynolds-Averaged Turbulence Modeling Using Deep Learning with Local Flow Features: An Empirical Approach

Reynolds-Averaged Navier-Stoke (RANS) models offer an alternative avenue in predicting flow characteristics when the corresponding experiments are difficult to achieve due to geometry complexity, limited budget, or knowledge. RANS models require the knowledge of subgrid scale physics to solve conservation equations for mass, energy, and momentum. Mechanistic turbulence models, such as k-epsilon, are generally evaluated and calibrated for specific flow conditions with various degrees of uncertainty. These models have limited capability to assimilate a substantial amount of data due to model form constraints. Meanwhile, deep learning (DL) has been proven to be universal approximators with the potential to assimilate available, relevant, and adequately evaluated data. Moreover, deep neural networks (DNNs) can create surrogate models without knowing function forms. Such a data-driven approach can be used in updating fluid models based on observations as opposed to hard-wiring models with precalibrated correlations. The paper presents progress in applying DNNs to model Reynolds stress using two machine learning (ML) frameworks. A novel flow feature coverage mapping is proposed to quantify the physics coverage of DL-based closures. It can be used to examine the sufficiency of training data and input flow features for data-driven turbulence models. The case of a backward-facing step is formulated to demonstrate that not only can DNNs discover underlying correlation behind fluid data but also they can be implemented in RANS to predict flow characteristics without numerical stability issues. Finally, the presented research is a crucial stepping-stone toward the data-driven turbulence modeling, which potentially benefits the design of data-driven experiments that can be used to validate fluid models with ML-based fluid closures.

42 ENGINEERING↗
NASA NTRS1973

An analytical study of the response of a constant-attitude aircraft to atmospheric turbulence

A light airplane equipped with an automatic control system which drives large wing flaps and the stabilator so as to produce a constant pitch attitude in all flight modes was analyzed for its response to a specific gust. The aircraft was also equipped with a bank-angle steering, zero sideslip automatic control system which was studied for its effectiveness in suppressing a specific lateral gust. The gusts were assumed to be comprised of 200 lateral and 400 vertical sinusoids. Each was used to excite the controlled aircraft and the time response to the sum of all sinusoids was plotted. The assumption was that the gust may be treated as stationary in space but variable in time rather than the reverse. Results indicate that such a control system can suppress vertical gusts up to the limit of control authority. Either the lateral accelerations or the yawing velocity response to lateral gusts can be suppressed with this system but not both simultaneously.

Smetana, F. O.↗