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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗
NASA NTRS1980

Maximum likelihood method for estimating airplane stability and control parameters from flight data in frequency domain

A frequency domain maximum likelihood method is developed for the estimation of airplane stability and control parameters from measured data. The model of an airplane is represented by a discrete-type steady state Kalman filter with time variables replaced by their Fourier series expansions. The likelihood function of innovations is formulated, and by its maximization with respect to unknown parameters the estimation algorithm is obtained. This algorithm is then simplified to the output error estimation method with the data in the form of transformed time histories, frequency response curves, or spectral and cross-spectral densities. The development is followed by a discussion on the equivalence of the cost function in the time and frequency domains, and on advantages and disadvantages of the frequency domain approach. The algorithm developed is applied in four examples to the estimation of longitudinal parameters of a general aviation airplane using computer generated and measured data in turbulent and still air. The cost functions in the time and frequency domains are shown to be equivalent; therefore, both approaches are complementary and not contradictory. Despite some computational advantages of parameter estimation in the frequency domain, this approach is limited to linear equations of motion with constant coefficients.

Klein, V.↗
NASA NTRS1992

Code validation for the simulation of supersonic viscous flow about the F-16XL

The viewgraphs and discussion on code validation for the simulation of supersonic viscous flow about the F-16XL are provided. Because of the large potential gains related to laminar flow on the swept wings of supersonic aircraft, interest in the applications of laminar flow control (LFC) techniques in the supersonic regime has increased. A supersonic laminar flow control (SLFC) technology program is currently underway within NASA. The objective of this program is to develop the data base and design methods that are critical to the development of laminar flow control technology for application to supersonic transport aircraft design. Towards this end, the program integrates computational investigations underway at NASA Ames-Moffett and NASA Langley with flight-test investigations being conducted on the F-16XL at the NASA Ames-Dryden Research Facility in cooperation with Rockwell International. The computational goal at NASA Ames-Moffett is to integrate a thin-layer Reynolds averaged Navier-Stokes flow solver with a stability analysis code. The flow solver would provide boundary layer profiles to the stability analysis code which in turn would predict transition on the F-16XL wing. To utilize the stability analysis codes, reliable boundary layer data is necessary at off-design cases. Previously, much of the prediction of boundary layer transition has been accomplished through the coupling of boundary layer codes with stability theory. However, boundary layer codes may have difficulties at high Reynolds numbers, of the order of 100 million, and with the current complex geometry in question. Therefore, a reliable code which solves the thin-layer Reynolds averaged Navier-Stokes equations is needed. Two objectives are discussed, the first in greater depth. The first objective is method verification, via comparisons of computations with experiment, of the reliability and robustness of the code. To successfully implement LFC techniques to the F-16XL wing, the flow about the leading edge must be maintained as laminar flow. Therefore, the second objective is to focus on a series of numerical simulations with different values of angle of attack, alpha, and Reynolds numbers. The purpose of the simulations is to study their effects on the two main factors which precipitate transition to turbulence at leading edges of highly swept wings (e.g., 'spanwise contamination' and 'crossflow instability').

Flores, Jolen↗
NASA NTRS1994

Triple flame structure and diffusion flame stabilization

The stabilization of diffusion flames is studied using asymptotic techniques and numerical tools. The configuration studied corresponds to parallel streams of cold oxidizer and fuel initially separated by a splitter plate. It is shown that stabilization of a diffusion flame may only occur in this situation by two processes. First, the flame may be stabilized behind the flame holder in the wake of the splitter plate. For this case, numerical simulations confirm scalings previously predicted by asymptotic analysis. Second, the flame may be lifted. In this case a triple flame is found at longer distances downstream of the flame holder. The structure and propagation speed of this flame are studied by using an actively controlled numerical technique in which the triple flame is tracked in its own reference frame. It is then possible to investigate the triple flame structure and velocity. It is shown, as suggested from asymptotic analysis, that heat release may induce displacement speeds of the triple flame larger than the laminar flame speed corresponding to the stoichiometric conditions prevailing in the mixture approaching the triple flame. In addition to studying the characteristics of triple flames in a uniform flow, their resistance to turbulence is investigated by subjecting triple flames to different vortical configurations.

Veynante, D.↗
NASA NTRS1986

Analysis and interpretation of the shocklike electrostatic noise observed during the AMPTE solar wind lithium releases

During the AMPTE (Active Magnetospheric Particle Tracer Explorers) solar wind lithium release on September 11, and again on September 20, 1984, an intense burst of electrostatic noise was observed near the upstream edge of the ion cloud. Comparisons with measurements by the IMP-6 and ISEE-1 spacecraft show that the spectrum and overall features of this noise are very similar to electrostatic noise observed at the Earth's bow shock. A stability analysis using realistic parameters shows that the electrostatic noise can be accounted for by an ion beam-plasma instability caused by the solar wind proton beam streaming through the nearly stationary lithium cloud. The growth rate of this instability is largest when the ion density and solar wind proton density are similar, which explains why the noise only occurs near the outer edge of the ion cloud. The similarity to the noise in the Earth's bow shock suggest that a shock may exist in the solar wind plasma flow upstream of the ion cloud. If the noise is associated with a shock, then it must be an electrostatic shock, since the ion cyclotron radii are too small for the existence of a MHD shock. Since the electrostatic instability occurs at phase velocities near the lithium thermal velocity, the electrostatic turbulence may play a role in heating the lithium ions and transferring momentum from the solar wind to the ion cloud.

Gurnett, D. A.↗
NASA NTRS2021

Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition

We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.

Transition↗
NASA NTRS2011

Control of Stationary Cross-Flow Modes in a Mach 3.5 Boundary Layer Using Patterned Passive and Active Roughness

Spanwise-periodic roughness designed to excite selected wavelengths of stationary cross- ow modes was investigated in a 3-D boundary layer at Mach 3.5. The test model was a sharp-tipped 14deg right-circular cone. The model and integrated sensor traversing system were placed in the Mach 3.5 Supersonic Low Disturbance Tunnel (SLDT) equipped with a "quiet design" nozzle at the NASA Langley Research Center. The model was oriented at a 4:2deg angle of attack to produce a mean cross-fl ow velocity component in the boundary layer over the cone. Five removable cone tips have been investigated. One has a smooth surface that is used to document the baseline ("natural") conditions. Two had minute (20 - 40 micron) "dimples" that are equally spaced around the circumference, at a streamwise location that is just upstream of the linear stability neutral growth branch for cross- ow modes. The azimuthal mode numbers of the dimpled tips were selected to either enhance the most amplified wave numbers, or to suppress the growth of the most amplified wave numbers. Two of the cone tips had an array of plasma streamwise vortex generators that were designed to simulate the disturbances produced by the passive patterned roughness. The results indicate that the stationary cross-fl ow modes were highly receptive to the patterned roughness of both passive and active types. The patterned passive roughness that was designed to suppress the growth of the most amplified modes had an azimuthal wavelength that was 66% smaller that that of the most amplified stationary cross- ow mode. This had the effect to increase the transition Reynolds number from 25% to 50% depending on the measurement technique. The application of the research is on turbulent transition control on swept wings of supersonic aircraft. The plasma-based roughness has the advantage over the passive roughness of being able to be adaptable to different conditions that would occur during a flight mission.

Schuele, Chan Yong↗
NASA NTRS1990

Active turbulence control in wallbounded flow using direct numerical simulations

An exploratory study of concepts for active control of turbulent boundary layers using the direct numerical simulation technique was performed. Significant drag reduction was achieved when the surface boundary condition was modified such that it could suppress the large-scale structures present in the wall region. This was achieved by prescribing the normal component of velocity at the wall to be 180 deg out of phase with the normal velocity slightly above the wall at each instant. The drag reduction was accompanied with significant reduction in the intensity of the wall-layer structures and reductions in the magnitude of Reynolds stresses throughout the flow. Suitability of wall-pressure and shear-stress fluctuations for detection of flow structures above the wall was examined. A preliminary result obtained by applying the present control strategy to a transitional flow is also briefly described, from which one can infer a possible linkage between the control strategy and flow stability.

Kim, J.↗
NASA NTRS2011

Assessing Uncertainties in Boundary Layer Transition Predictions for HIFiRE-1 at Non-zero Angles of Attack

Boundary layer stability was analyzed for the HIFiRE-1 flight vehicle geometry for ground tests conducted at the CUBRC LENS I hypersonic shock test facility and the Langley Research Center (LaRC) 20- inch Mach 6 Tunnel. Boundary layer stability results were compared to transition onset location obtained from discrete heat transfer measurements from thin film gauges during the CUBRC test and spatially continuous heat transfer measurements from thermal phosphor paint data during the LaRC test. The focus of this analysis was on conditions at non-zero angles of attack as stability analysis has already been performed at zero degrees angle of attack. Also, the transition onset data obtained during flight testing was at nonzero angles of attack, so this analysis could be expanded in the future to include the results of the flight test data. Stability analysis was performed using the 2D parabolized stability software suite STABL (Stability and Transition Analysis for Hypersonic Boundary Layers) developed at the University of Minnesota and the mean flow solutions were computed using the DPLR finite volume Navier-Stokes computational fluid dynamics (CFD) solver. A center line slice of the 3D mean flow solution was used for the stability analysis to incorporate the angle of attack effects while still taking advantage of the 2D STABL software suite. The N-factors at transition onset and the value of Re(sub theta)/M(sub e), commonly used to predict boundary layer transition onset, were compared for all conditions analyzed. Ground test data was analyzed at Mach 7.2 and Mach 6.0 and angles of attack of 1deg, 3deg and 5deg. At these conditions, the flow was found to be second mode dominant for the HIFiRE-1 slender cone geometry. On the leeward side of the vehicle, a strong trend of transition onset location with angle of attack was observed as the boundary layer on the leeward side of the vehicle developed inflection points at streamwise positions on the vehicle that correlated to angle of attack. Inflection points are a strong instability mechanism that lead to rapid breakdown and transition to turbulence. The transition onset location on the windward side of the vehicle displayed no trend with angle of attack or freestream Reynolds number and transition was observed farther down the vehicle than observed on the leeward side of the vehicle. In analysis of both windward and leeward sides of the vehicle, use of the N factor methodology to develop trends to predict boundary layer transition onset showed improvements over the Re(sub theta)/M(sub e) empirical correlation methodology. Stronger correlations and less scatter in the data were observed when using the N factor method for these cases.

Marek, Lindsay C.↗
NASA NTRS1994

Wind tunnel test of a variable-diameter tiltrotor (VDTR) model

This report documents the results from a wind tunnel test of a 1/6th scale Variable Diameter Tiltrotor (VDTR). This test was a joint effort of NASA Ames and Sikorsky Aircraft. The objective was to evaluate the aeroelastic and performance characteristics of the VDTR in conversion, hover, and cruise. The rotor diameter and nacelle angle of the model were remotely changed to represent tiltrotor operating conditions. Data is presented showing the propulsive force required in conversion, blade loads, angle of attack stability and simulated gust response, and hover and cruise performance. This test represents the first wind tunnel test of a variable diameter rotor applied to a tiltrotor concept. The results confirm some of the potential advantages of the VDTR and establish the variable diameter rotor a viable candidate for an advanced tiltrotor. This wind tunnel test successfully demonstrated the feasibility of the Variable Diameter rotor for tilt rotor aircraft. A wide range of test points were taken in hover, conversion, and cruise modes. The concept was shown to have a number of advantages over conventional tiltrotors such as reduced hover downwash with lower disk loading and significantly reduced longitudinal gust response in cruise. In the conversion regime, a high propulsive force was demonstrated for sustained flight with acceptable blade loads. The VDTR demonstrated excellent gust response capabilities. The horizontal gust response correlated well with predictions revealing only half the response to turbulence of the conventional civil tiltrotor.

Matuska, David↗
NASA NTRS2002

Boundary Layer Transition on Slender Cones in Conventional and Low Disturbance Mach 6 Wind Tunnels

An experimental investigation was conducted on a 5-degree half-angle cone and a 5-degree half-angle flared cone in a conventional Mach 6 wind tunnel to examine the effects of facility noise on boundary layer transition. The influence of tunnel noise was inferred by comparing transition onset locations determined from the present test to that previously obtained in a Mach 6 low disturbance quiet tunnel. Together, the two sets of experiments are believed to represent the first direct comparison of transition onset between a conventional and a low disturbance wind tunnel using a common test model and transition detection technique. In the present conventional hypersonic tunnel experiment, separate measurements of heat transfer and adiabatic wall temperatures were obtained on the conical models at small angles of attack over a range of Reynolds numbers, which resulted in laminar, transitional, and turbulent flow. Smooth model turbulent heating distributions are compared to that obtained with transition forced via discrete surface roughness. The model nosetip radius was varied to examine the effects of bluntness on transition onset. Despite wall to total temperature differences between the transient heating measurements and the adiabatic wall temperature measurement, the two methods for determining sharp cone transition onset generally yielded equivalent locations. In the 'noisy' mode of the hypersonic low disturbance tunnel, transition onset occurred earlier than that measured in the conventional hypersonic tunnel, suggesting higher levels of freestream acoustic radiation relative to the conventional tunnel. At comparable freestream conditions, the transition onset Reynolds number under low disturbance conditions was a factor of 1.3 greater than that measured on flared cone in the LaRC conventional hypersonic tunnel and a factor of 1.6 greater that the flared cone run in the low disturbance tunnel run 'noisy'. Navier-Stokes mean flow computations and linear stability analysis were conducted to assess the experimental results and have indicated N factors associated with sharp flared cone transition onset to be approximately a factor of 2 lower than that inferred from the corresponding low disturbance tunnel measurements.

Horvath, Thomas J.↗
NASA NTRS1999

Numerical Solution of Incompressible Navier-Stokes Equations Using a Fractional-Step Approach

A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.

Kiris, Cetin↗
NASA NTRS2016

Unsteady Heat-Flux Measurements of Second-Mode Instability Waves in a Hypersonic Boundary Layer

In this paper we report on the application of the atomic layer thermopile (ALTP) heat- flux sensor to the measurement of laminar-to-turbulent transition in a hypersonic flat plate boundary layer. The centerline of the flat-plate model was instrumented with a streamwise array of ALTP sensors and the flat-plate model was exposed to a Mach 6 freestream over a range of unit Reynolds numbers. Here, we observed an unstable band of frequencies that are associated with second-mode instability waves in the laminar boundary layer that forms on the flat-plate surface. The measured frequencies, group velocities, phase speeds, and wavelengths of these instability waves are in agreement with data previously reported in the literature. Heat flux time series, and the Morlet-wavelet transforms of them, revealed the wave-packet nature of the second-mode instability waves. In addition, a laser-based radiative heating system was developed to measure the frequency response functions (FRF) of the ALTP sensors used in the wind tunnel test. These measurements were used to assess the stability of the sensor FRFs over time and to correct spectral estimates for any attenuation caused by the finite sensor bandwidth.

Kergerise, Michael A.↗
NASA NTRS1984

Design of infrasound-detection system via adaptive LMSTDE algorithm

A proposed solution to an aviation safety problem is based on passive detection of turbulent weather phenomena through their infrasonic emission. This thesis describes a system design that is adequate for detection and bearing evaluation of infrasounds. An array of four sensors, with the appropriate hardware, is used for the detection part. Bearing evaluation is based on estimates of time delays between sensor outputs. The generalized cross correlation (GCC), as the conventional time-delay estimation (TDE) method, is first reviewed. An adaptive TDE approach, using the least mean square (LMS) algorithm, is then discussed. A comparison between the two techniques is made and the advantages of the adaptive approach are listed. The behavior of the GCC, as a Roth processor, is examined for the anticipated signals. It is shown that the Roth processor has the desired effect of sharpening the peak of the correlation function. It is also shown that the LMSTDE technique is an equivalent implementation of the Roth processor in the time domain. A LMSTDE lead-lag model, with a variable stability coefficient and a convergence criterion, is designed.

Khalaf, C. S.↗
NASA NTRS1993

Convergence acceleration of implicit schemes in the presence of high aspect ratio grid cells

The performance of Navier-Stokes codes are influenced by several phenomena. For example, the robustness of the code may be compromised by the lack of grid resolution, by a need for more precise initial conditions or because all or part of the flowfield lies outside the flow regime in which the algorithm converges efficiently. A primary example of the latter effect is the presence of extended low Mach number and/or low Reynolds number regions which cause convergence deterioration of time marching algorithms. Recent research into this problem by several workers including the present authors has largely negated this difficulty through the introduction of time-derivative preconditioning. In the present paper, we employ the preconditioned algorithm to address convergence difficulties arising from sensitivity to grid stretching and high aspect ratio grid cells. Strong grid stretching is particularly characteristic of turbulent flow calculations where the grid must be refined very tightly in the dimension normal to the wall, without a similar refinement in the tangential direction. High aspect ratio grid cells also arise in problems that involve high aspect ratio domains such as combustor coolant channels. In both situations, the high aspect ratio cells can lead to extreme deterioration in convergence. It is the purpose of the present paper to address the reasons for this adverse response to grid stretching and to suggest methods for enhancing convergence under such circumstances. Numerical algorithms typically possess a maximum allowable or optimum value for the time step size, expressed in non-dimensional terms as a CFL number or vonNeumann number (VNN). In the presence of high aspect ratio cells, the smallest dimension of the grid cell controls the time step size causing it to be extremely small, which in turn results in the deterioration of convergence behavior. For explicit schemes, this time step limitation cannot be exceeded without violating stability restrictions of the scheme. On the other hand, for implicit schemes, which are typically unconditionally stable, there appears to be room for improvement through careful tailoring of the time step definition based on results of linear stability analyses. In the present paper, we focus on the central-differenced alternating direction implicit (ADI) scheme. The understanding garnered from this analyses can then be applied to other implicit schemes. In order to systematically study the effects of aspect ratio and the methods of mitigating the associated problems, we use a two pronged approach. We use stability analyses as a tool for predicting numerical convergence behavior and numerical experiments on simple model problems to verify predicted trends. Based on these analyses, we determine that efficient convergence may be obtained at all aspect ratios by getting a combination of things right. Primary among these are the proper definition of the time step size, proper selection of viscous preconditioner and the precise treatment of boundary conditions. These algorithmic improvements are then applied to a variety of test cases to demonstrate uniform convergence at all aspect ratios.

Buelow, B. E. O.↗
NASA NTRS2011

Effects of Cavities and Protuberances on Transition over Hypersonic Vehicles

Surface protuberances and cavities on a hypersonic vehicle are known to cause several aerodynamic or aerothermodynamic issues. Most important of all, premature transition due to these surface irregularities can lead to a significant rise in surface heating. To help understand laminar-turbulent transition induced by protuberances or cavities on a Crew Exploration Vehicle (CEV) surface, high-fidelity numerical simulations are carried out for both types of trips on a CEV wind tunnel model. Due to the large bluntness, these surface irregularities reside in an accelerating subsonic boundary layer. For the Mach 6 wind tunnel conditions with a roughness Reynolds number Re(sub kk) of 800, it was found that a protuberance with a height to boundary layer thickness ratio of 0.73 leads to strong wake instability and spontaneous vortex shedding, while a cavity with identical geometry only causes a rather weak flow unsteadiness. The same cavity with a larger Reynolds number also leads to similar spontaneous vortex shedding and wake instability. The wake development and the formation of hairpin vortices for both protuberance and cavity were found to be qualitatively similar to that observed for an isolated hemisphere submerged in a subsonic, low speed flat-plate boundary layer. However, the shed vortices and their accompanying instability waves were found to be slightly stabilized downstream by the accelerating boundary layer along the CEV surface. Despite this stabilizing influence, it was found that the wake instability spreads substantially in both wall-normal and azimuthal directions as the flow is evolving towards a transitional state. Similarities and differences between the wake instability behind a protuberance and a cavity are investigated. Computations for the Mach 6 boundary layer over a slender cylindrical roughness element with a height to the boundary layer thickness of about 1.1 also shows spontaneous vortex shedding and strong wake instability. Comparisons of detailed flow structures associated with protuberances at subsonic and supersonic edge Mach numbers indicate distinctively different instability mechanisms.

Chang, Chau-Lyan↗
NASA NTRS1957

A discussion of cone and flat-plate Reynolds numbers for equal ratios of the laminar shear to the shear caused by small velocity fluctuations in a laminar boundary layer

By use of the linear theory of boundary-layer stability and Schlichting's formula for the maximum amplification of a disturbance, an approximate relation is derived between the Reynolds number on a cone and the Reynolds number on a flat plate for equal closeness to transition. The indication is that the ratio of the cone Reynolds number for transition, based on the distance to the cone apex, to the plate Reynolds number for transition, based on the distance to the leading edge, is not in general equal to 3, as has been suggested by other investigators, but varies from 3 when transition occurs at the minimum critical Reynolds number to unity when transition occurs at a large multiple of the critical Reynolds number.

FLOW, TURBULENT↗
NASA NTRS1999

Hydrodynamic Stability and Magnetic Reconnection in Disks and Stars

The purpose of this grant is to study parametric instability. The simplest example of parametric instability is a harmonic oscillator with a periodic modulation of the spring constant. If the modulation frequency is close to twice the natural frequency of the oscillator, the amplitude of oscillation tends to grow exponentially. The growth rate is proportional to the strength of the modulation, but it also depends upon the closeness to resonance of the two frequencies, and upon natural damping rate or "Q" of the oscillator. Parametric instabilities are very common in physics. A familiar example is a jogger's ponytail--normally a very strongly damped pendulum, it can be destabilized by the variation in effective gravity during the jogger's stride. Observation confirms that the period of the pendulum is half that of the jogger's vertical motion. In astrophysics, parametric instability may occur by external tidal forcing, or by interaction among eigenmodes. In the latter case, an energetic eigenmode may destabilize modes of half its frequency, provided some weak nonlinearity exists to couple them. Under a previous Astrophysical Theory grant (NAGW-2419), the PI discovered a parametric instability of tidally forced disks such as the accretion disks in cataclysmic variables and X ray binaries [2]. The destabilized modes are tightly-wound, incompressible, three-dimensional waves analogous to g-modes and r-modes in stars. Later work has confirmed our analysis [4]. It was hoped that these modes might provide a source of turbulence and angular momentum transport in accretion disks. However, a follow-up investigation of this instability by local numerical simulations, although confirming the analytically estimated growth rates, found negligible angular momentum flux [3]. Other work, partly supported by the ATP, now strongly indicates that the transport mechanism in such disks is magnetohydrodynamic turbulence [6]. Nevertheless, the parametric mechanism may truncate the outer edges of disks in close binaries [2], and it may be important in disks of very low ionization such as protostellar disks, or even cataclysmic-variable disks in quiescence where the MHD mechanism may be ineffective [5]. All analyses up to 1996 were done in a local approximation where the orbital frequency, shear rate, and tidal field were treated as constants. The locally computed growth rate turns out to depend strongly on radius, and it was unclear how to average these local rates to obtain the correct global rate. This is a critical issue for accretion disks in close binaries, because the local growth rate is comparable to the orbital frequency towards the outer edge of the disk but decreases rapidly inwards. Paper #1 examined this issue in a simplified global model where the destabilizing terms vary with position. We found that the global growth rate is essentially equal to the maximum local rate, provided that the latter is smoothed over a radial range equal to the distance that the destabilized wave propagates at its group speed in one growth time. Thus, in an accretion disk, waves would grow rapidly in the outer parts but would propagate both inwards and outwards at a maximum group speed of order the disk thickness divided by the orbital period.

Goodman, Jeremy↗