Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow instability”

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

Pretest Computational Assessment of Boundary Layer Transition in the NASA Juncture Flow Model with an NACA 0015-Based Wing

The first two phases of the NASA Juncture Flow experiment were carried out on a DLR-F6 swept-wing model and were designed to provide “CFD validation-quality” data toward the assessment and improvement of existing CFD turbulence models in predicting onset and extent of three-dimensional separated flow near the wing-juncture trailing-edge region. The next phase of experiments will involve an NACA 0015-based swept wing, as prior risk reduction experiments had indicated that this wing shape resulted in reduced separation near the juncture region than the DLR-F6 wing, thus providing a better option to evaluate the ability of CFD models to predict incipient turbulent separation. The NACA 0015 measurements will also include IR thermography to infer the variation of transition front with respect to an increasing angle of attack. The primary objective of this work is to computationally make a preliminary assessment of the transition front on both surfaces of the NACA 0015 wing at a crank-chord-based Reynolds number of 2.4 x 106 for four different angles of attack, (0°, 2.5°, 5°, and 7.5°) and to determine the dominant mechanisms responsible for transition. This assessment includes both RANS-based transition models from NASA’s OVERFLOW 2.3b flow solver and linear parabolized stability equations (PSE) stability analysis based on the Langley Stability and Transition Analysis code, LASTRAC. Linear PSE results indicate that the upper surface of the wing is dominated by Tollmien- Schlichting (TS) instabilities, and that the laminar flow region shrinks from about 50% chord to a very small region just downstream of the attachment line as the angle of attack is increased from 0° to 7.5°. Consequently, the transition fronts predicted by the Spalart- Allmaras-based amplification factor transport (AFT-2017b) equation model (which accounts for the TS instabilities alone) and the Menter’s shear-stress transport equation (SST2003)- based Langtry-Menter transition model with ability to account for both TS and crossflow effects (LM2015) compare well with those predicted using linear PSE. On the lower surface of the wing, stationary crossflow (CF) instabilities begin to appear on the inboard portion of the wing in addition to the TS-instabilities for the larger angles of attack (5° and 7.5°), further reducing the laminar flow extent within the inboard region. The LM2015 model that accounts for CF effects is able to replicate this trend but appears to predict a slightly earlier transition. The outcome of this effort will inform the experiment and, when the actual experimental data become available, provide further opportunity to assess and improve the various transition models.

CFD modeling↗

Pretest Computational Assessment of Boundary Layer Transition in the NASA Juncture Flow Model with an NACA 0015-Based Wing

The first two phases of the NASA Juncture Flow experiment were carried out on a DLR-F6 swept-wing model and were designed to provide “CFD validation-quality” data toward the assessment and improvement of existing CFD turbulence models in predicting onset and extent of three-dimensional separated flow near the wing-juncture trailing-edge region. The next phase of experiments will involve an NACA 0015-based swept wing, as prior risk reduction experiments had indicated that this wing shape resulted in reduced separation near the juncture region than the DLR-F6 wing, thus providing a better option to evaluate the ability of CFD models to predict incipient turbulent separation. The NACA 0015 measurements will also include IR thermography to infer the variation of transition front with respect to an increasing angle of attack. The primary objective of this work is to computationally make a preliminary assessment of the transition front on both surfaces of the NACA 0015 wing at a crank-chord-based Reynolds number of 2.4 x 106 for four different angles of attack, (0°, 2.5°, 5°, and 7.5°) and to determine the dominant mechanisms responsible for transition. This assessment includes both RANS-based transition models from NASA’s OVERFLOW 2.3b flow solver and linear parabolized stability equations (PSE) stability analysis based on the Langley Stability and Transition Analysis code, LASTRAC. Linear PSE results indicate that the upper surface of the wing is dominated by Tollmien- Schlichting (TS) instabilities, and that the laminar flow region shrinks from about 50% chord to a very small region just downstream of the attachment line as the angle of attack is increased from 0° to 7.5°. Consequently, the transition fronts predicted by the Spalart- Allmaras-based amplification factor transport (AFT-2017b) equation model (which accounts for the TS instabilities alone) and the Menter’s shear-stress transport equation (SST2003)- based Langtry-Menter transition model with ability to account for both TS and crossflow effects (LM2015) compare well with those predicted using linear PSE. On the lower surface of the wing, stationary crossflow (CF) instabilities begin to appear on the inboard portion of the wing in addition to the TS-instabilities for the larger angles of attack (5° and 7.5°), further reducing the laminar flow extent within the inboard region. The LM2015 model that accounts for CF effects is able to replicate this trend but appears to predict a slightly earlier transition. The outcome of this effort will inform the experiment and, when the actual experimental data become available, provide further opportunity to assess and improve the various transition models.

CFD modeling↗

Streamwise vortices in heated boundary layers

The nonlinear instability of the boundary layer on a heated flat plate placed in an oncoming flow is investigated. Such flows are unstable to stationary vortex instabilities and inviscid traveling wave disturbances governed by the Taylor-Goldstein equation. For small temperature differences the Taylor-Goldstein equation reduces to Rayleigh's equation. When the temperature difference between the wall and free stream is small the preferred mode of instability is a streamwise vortex. It is shown in this case that the vortex, assumed to be of small wavelength, restructures the underlying mean flow to produce a profile which can be massively unstable to inviscid traveling waves. The mean state is shown to be destabilized or stabilized to inviscid waves depending on whether the Prandtl number is less or greater than unity.

Hall, Philip↗

Streamwise vortices in heated boundary layers

The nonlinear instability of the boundary layer on a heated flat plate placed in an oncoming flow is investigated. Such flows are unstable to stationary vortex instabilities and inviscid travelling wave disturbances governed by the Taylor-Goldstein equation. For small temperature differences the Taylor-Goldstein equation reduces to Rayleigh's equation. When the temperature difference between the wall and free stream is small the preferred mode of instability is a streamwise vortex. It is shown in this case that the vortex, assumed to be of small wavelength, restructures the underlying mean flow to produce a profile which can be massively unstable to inviscid travelling waves. The mean state is shown to be destabilized if the Prandtl number is less than unity.

Hall, Philip↗

Studies on Equatorial Shock Formation During Plasmaspheric Refilling

During the grant period from August 1, 1994 to October 31, 1995 we have continued to investigate the effects of plasma wave instabilities on the early stage plasmaspheric refilling. Since ion beams are the primary feature of the interhemispheric plasma flows during the early stage refilling, ion-beam driven instabilities and associated waves are of primary interest. The major findings of this research are briefly summarized here. After a systematic examination of the relevant plasma instabilities, we realized that when the interhemispheric plasma flows begin to interpenetrate at the equator, the most relevant plasma instability is the electrostatic ion cyclotron wave instability. Only at later stages the ion-acoustic instability may be affecting the plasma flow. An interesting property of the electrostatic ion cyclotron wave is that it heats ions perpendicular to the magnetic field. When the ions in the field-aligned flows are transversely heated, they are trapped in the magnetic flux tube, thus affecting the refilling process. The eic wave instability is a microprocess with scale length of the order of ion Larmor radius and the corresponding time scale is the ion cyclotron period. We have attempted to tackle the problem for the plasmaspheric refilling by incorporating the effects of eic wave instability on the mesoscale plasma flow when the properties of the latter exceeds the critical conditions for the former. We have compared the results on refilling from the model with and without the eic instability effects.

Singh, Nagendra↗

An Experimental Investigation of Incompressible Richtmyer-Meshkov Instability

Richtmyer-Meshkov (RM) instability occurs when two different density fluids are impulsively accelerated in the direction normal to their nearly planar interface. The instability causes small perturbations on the interface to grow and eventually become a turbulent flow. It is closely related to Rayleigh-Taylor instability, which is the instability of a planar interface undergoing constant acceleration, such as caused by the suspension of a heavy fluid over a lighter one in the earth's gravitational field. Like the well-known Kelvin-Helmholtz instability, RM instability is a fundamental hydrodynamic instability which exhibits many of the nonlinear complexities that transform simple initial conditions into a complex turbulent flow. Furthermore, the simplicity of RM instability (in that it requires very few defining parameters), and the fact that it can be generated in a closed container, makes it an excellent test bed to study nonlinear stability theory as well as turbulent transport in a heterogeneous system. However, the fact that RM instability involves fluids of unequal densities which experience negligible gravitational force, except during the impulsive acceleration, requires RM instability experiments to be carried out under conditions of microgravity. This experimental study investigates the instability of an interface between incompressible, miscible liquids with an initial sinusoidal perturbation. The impulsive acceleration is generated by bouncing a rectangular tank containing two different density liquids off a retractable vertical spring. The initial perturbation is produced prior to release by oscillating the tank in the horizontal direction to produce a standing wave. The instability evolves in microgravity as the tank travels up and then down the vertical rails of a drop tower until hitting a shock absorber at the bottom. Planar Laser Induced Fluorescence (PLIF) is employed to visualize the flow. PLIF images are captured by a video camera that travels with the tank. Figure 1 is as sequence of images showing the development of the instability from the initial sinusoidal disturbance far into the nonlinear regime which is characterized by the appearance of mushroom structures resulting from the coalescence of baroclinic vorticity produced by the impulsive acceleration. At later times in this sequence the vortex cores are observed to become unstable showing the beginnings of the transition to turbulence in this flow. The amplitude of the growing disturbance after the impulsive acceleration is measured and found to agree well with theoretical predictions. The effects of Reynolds number (based on circulation) on the development of the vortices and the transition to turbulence are also determined.

Jacobs, J. W.↗

The evolution of instabilities in the axisymmetric jet. I - The linear growth of disturbances near the nozzle. II - The flow resulting from the interaction between two waves

The modal distribution of coherent structures evolving near the nozzle of a circular jet was studied experimentally and theoretically, with particular attention given to the effects produced on the instability modes by transverse curvature, flow divergence, inhomogeneous inflow conditions, and the detailed shape of the mean velocity profile. Experiments were performed using a specially constructed air-jet facility; hot-wire anemometers were used in conjunction with Disa Model 55P11 sensors for flow measurements. The linear model used as a transfer function is capable of predicting the spectral distribution of the velocity perturbations in a jet. Consideration was also given to studies of leading nonlinear interactions generated by waves externally superimposed on an axisymmetric jet; theoretical predictions were verified experimentally.

Cohen, J.↗

Visualization studies of turbulent transition flows in a porous medium

Results are reported for flow-visualization studies of the flow regimes of water passing through a porous medium consisting of cylindrical glass and plexiglas rods arranged in a complex and fixed three-dimensional geometry. The Reynolds number (Re) varied from 50 to 700; the flow was visualized by injecting a 5% potassium permanganate dye solution into the pores and photographing the resulting dye streaklines with both a still camera and a movie camera. The results indicate that four distinct flow regimes exist in the porous medium: (1) Darcy or creeping flow up to Re = 3; (2) steady inertia-dominated laminar flow for Re = 3-150; (3) unsteady transitional laminar flow for Re = 150-250; and (4) fully turbulent flow for Re greater than 250. It is concluded that a laminar wake instability mechanism typical of the external flow about bluff bodies may be responsible for the overall transition from laminar to turbulent flow in porous media.

Bilardo, V. J.↗

On a class of unsteady three-dimensional Navier Stokes solutions relevant to rotating disc flows: Threshold amplitudes and finite time singularities

A class of exact steady and unsteady solutions of the Navier Stokes equations in cylindrical polar coordinates is given. The flows correspond to the motion induced by an infinite disc rotating with constant angular velocity about the z-axis in a fluid occupying a semi-infinite region which, at large distances from the disc, has velocity field proportional to (x,-y,O) with respect to a Cartesian coordinate system. It is shown that when the rate of rotation is large, Karman's exact solution for a disc rotating in an otherwise motionless fluid is recovered. In the limit of zero rotation rate a particular form of Howarth's exact solution for three-dimensional stagnation point flow is obtained. The unsteady form of the partial differential system describing this class of flow may be generalized to time-periodic equilibrium flows. In addition the unsteady equations are shown to describe a strongly nonlinear instability of Karman's rotating disc flow. It is shown that sufficiently large perturbations lead to a finite time breakdown of that flow whilst smaller disturbances decay to zero. If the stagnation point flow at infinity is sufficiently strong, the steady basic states become linearly unstable. In fact there is then a continuous spectrum of unstable eigenvalues of the stability equations but, if the initial value problem is considered, it is found that, at large values of time, the continuous spectrum leads to a velocity field growing exponentially in time with an amplitude decaying algebraically in time.

Hall, Philip↗

Heat Entrapment Effects Within Liquid Acquisition Devices

We introduce a model problem to address heat entrapment effects or the local accumulation of thermal energy within liquid acquisition devices. We show that the parametric space consists of six parameters, namely the Rayleigh and Prandtl numbers, the aspect ratio, and heat flux ratios for the bottom, side, and top boundaries of the enclosure. For the range of Ra considered 1 to 10(sup 9), beyond Ra on the order of 10(sup 5), convective instability is the dominant mode of convection in comparison to natural convection. The flow field transitions to asymmetric modes at Ra on the order of 10(sup 7). Direct numerical simulation of a large geometric length scale prototype for Ra on the order of 10(sup 9) shows that the flow field evolves from small wavelength instability which gives rise to nonlinear growth of thermals, propagation of the instability occurs via growth of secondary and tertiary modes, and a travelling wave mode occurs prior to asymmetry. The effect of a large aspect ratio is to increase the number of modes in the vertical direction. Due to the slow diffusion of heat in the prototype, asymptotic states are not readily attained, we show that dynamical similarity can be used for a model which allows the attainment of asymptotic states and that transition to a chaotic state occurs for Ra on the order of 10(sup 9) via a broadband power spectrum. These dynamical events show that for the baseline condition in which heat is absorbed from background laboratory environment, higher heat flux is absorbed at the top and bottom boundaries of the enclosure than a nominal value of 34.9 ergs per square centimeter -second.

Duval, W. M. B.↗

Nonlinear Development and Secondary Instability of Traveling Crossflow Vortices

Transition research under NASA's Aeronautical Sciences Project seeks to develop a validated set of variable fidelity prediction tools with known strengths and limitations, so as to enable "sufficiently" accurate transition prediction and practical transition control for future vehicle concepts. This paper builds upon prior effort targeting the laminar breakdown mechanisms associated with stationary crossflow instability over a swept-wing configuration relevant to subsonic aircraft with laminar flow technology. Specifically, transition via secondary instability of traveling crossflow modes is investigated as an alternate scenario for transition. Results show that, for the parameter range investigated herein, secondary instability of traveling crossflow modes becomes insignificant in relation to the secondary instability of the stationary modes when the relative initial amplitudes of the traveling crossflow instability are lower than those of the stationary modes by approximately two orders of magnitudes or more. Linear growth predictions based on the secondary instability theory are found to agree well with those based on PSE and DNS, with the most significant discrepancies being limited to spatial regions of relatively weak secondary growth, i.e., regions where the primary disturbance amplitudes are smaller in comparison to its peak amplitude. Nonlinear effects on secondary instability evolution is also investigated and found to be initially stabilizing, prior to breakdown.

Li, Fei↗

Spatial Direct Numerical Simulation of Boundary-Layer Transition Mechanisms: Validation of PSE Theory

A study of instabilities in incompressible boundary-layer flow on a flat plate is conducted by spatial direct numerical simulation (DNS) of the Navier-Stokes equations. Here, the DNS results are used to critically evaluate the results obtained using parabolized stability equations (PSE) theory and to study mechanisms associated with breakdown from laminar to turbulent flow. Three test cases are considered: two-dimensional Tollmien-Schlichting wave propagation, subharmonic instability breakdown, and oblique-wave break-down. The instability modes predicted by PSE theory are in good quantitative agreement with the DNS results, except a small discrepancy is evident in the mean-flow distortion component of the 2-D test problem. This discrepancy is attributed to far-field boundary- condition differences. Both DNS and PSE theory results show several modal discrepancies when compared with the experiments of subharmonic breakdown. Computations that allow for a small adverse pressure gradient in the basic flow and a variation of the disturbance frequency result in better agreement with the experiments.

Joslin, R. D.↗

Research highlights: June 1990 - May 1991

Linear instability calculations at MSFC have suggested that the Geophysical Fluid Flow Cell (GFFC) should exhibit classic baroclinic instability at accessible parameter settings. Interest was in the mechanisms of transition to temporal chaos and the evolution of spatio-temporal chaos. In order to understand more about such transitions, high resolution numerical experiments for the physically simplest model of two layer baroclinic instability were conducted. This model has the advantage that the numerical code is exponentially convergent and can be efficiently run for very long times, enabling the study of chaotic attractors without the often devastating effects of low-order trunction found in many previous studies. Numerical algorithms for implementing an empirical orthogonal function (EOF) analysis of the high resolution numerical results were completed. Under conditions of rapid rotation and relatively low differential heating, convection in a spherical shell takes place as columnar banana cells wrapped around the annular gap, but with axes oriented along the axis of rotation; these were clearly evident in the GFFC experiments. The results of recent numerical simulations of columnar convection and future research plans are presented.

Source record↗

The strong nonlinear interaction of Tollmien-Schlichting waves and Taylor-Goertler vortices in curved channel flow

Viscous fluid flows with curved streamlines can support both centrifugal and viscous traveling wave instabilities. Here the interaction of these instabilities in the context of the fully developed flow in a curved channel is discussed. The viscous (Tollmein-Schlichting) instability is described asymptotically at high Reynolds numbers and it is found that it can induce a Taylor-Goertler flow even at extremely small amplitudes. In this interaction, the Tollmein-Schlichting wave can drive a vortex state with wavelength either comparable with the channel width or the wavelength of lower branch viscous modes. The nonlinear equations which describe these interactions are solved for nonlinear equilibrium states.

Bennett, J.↗

Catalyst Bed Instability Within the USFE H2O2/JP-8 Rocket Engine

Orbital Sciences Corporation has been awarded a contract by NASA's Marshall Space Flight Center, in cooperation with the U.S. Air Force Research Laboratory's Military Space Plane Technology Program Office, for the Upper Stage Flight Experiment (USFE) program. Orbital is designing, developing, and will flight test a new low-cost, 10,000 lbf hydrogen peroxide/ JP-8 pressure fed liquid rocket. During combustion chamber tests at NASA Stennis Space Center (SSC) of the USFE engine, the catalyst bed showed a low frequency instability occurring as the H202 flow reached about 1/3 its design rate. This paper reviews the USFE catalyst bed and combustion chamber and its operation, then discusses the dynamics of the instability. Next the paper describes the dynamic computer model used to recreate the instability. The model was correlated to the SSC test data, and used to investigate possible solutions to the problem. The combustion chamber configuration which solved the instability is shown, and the subsequent stable operation presented.

Johnson, Curtis W.↗

Flow Visualization Study of a 1/48-Scale AFTI/F111 Model to Investigate Horizontal Tail Flow Disturbances

During flight testing of the AFTI/F111 aircraft, horizontal tail buffet was observed. Flutter analysis ruled out any aeroelastic instability, so a water-tunnel flow visualization study was conducted to investigate possible flow disturbances on the horizontal tail which might cause buffet. For this study, a 1/48-scale model was used. Four different wing cambers and one horizontal tail setting were tested between 0 and 20 deg angle of attack. These wing cambers corresponded to the following leading training edge deflections: 0/2, 10/10, 10/2, and 0/10. Flow visualization results in the form of still photographs are presented for each of the four wing cambers between 8 and 12 deg angle of attack. In general, the horizontal tail experiences flow disturbances which become more pronounced with angle of attack or wing trailing-edge deflection.

Bjarke, Lisa J.↗

Direct numerical study of crossflow instability

Disturbances in the swept Hiemenz flow are calculated by solving the Navier-Stokes equations. The spatially-evolving base flow is treated exactly, allowing a check of the 'local' stability theories. Different types of disturbances such as random noise, waves, and wave packets, are input near the attachment line, develop in space, and exit through an outflow boundary. They all generate streamwise vortices. The effect of the Reynolds number, of the time-dependence of the noise, and of nonlinearity, are investigated.

Spalart, P. R.↗

Crossflow Stability and Transition Experiments in a Swept-Wing Flow

An experimental examination of crossflow instability and transition on a 45 degree swept wing is conducted in the Arizona State University Unsteady Wind Tunnel. The stationary-vortex pattern and transition location are visualized using both sublimating-chemical and liquid-crystal coatings. Extensive hot-wire measurements are conducted at several measurement stations across a single vortex track. The mean and travelling-wave disturbances are measured simultaneously. Stationary-crossflow disturbance profiles are determined by subtracting either a reference or a span-averaged velocity profile from the mean-velocity data. Mean, stationary-crossflow, and travelling-wave velocity data are presented as local boundary-layer profiles and as contour plots across a single stationary-crossflow vortex track. Disturbance-mode profiles and growth rates are determined. The experimental data are compared to predictions from linear stability theory.

Dagenhart, John Ray↗