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 1,207 records · Page 67

MHz-Rate NO PLIF Imaging in a Mach 10 Hypersonic Wind Tunnel

NO PLIF imaging at repetition rates as high as 1 MHz is demonstrated in the NASA Langley 31 inch Mach 10 hypersonic wind tunnel. Approximately two hundred time correlated image sequences, of between ten and twenty individual frames, were obtained over eight days of wind tunnel testing spanning two entries in March and September of 2009. The majority of the image sequences were obtained from the boundary layer of a 20 flat plate model, in which transition was induced using a variety of cylindrical and triangular shaped protuberances. The high speed image sequences captured a variety of laminar and transitional flow phenomena, ranging from mostly laminar flow, typically at lower Reynolds number and/or in the near wall region of the model, to highly transitional flow in which the temporal evolution and progression of characteristic streak instabilities and/or corkscrew-shaped vortices could be clearly identified. A series of image sequences were also obtained from a 20 compression ramp at a 10 angle of attack in which the temporal dynamics of the characteristic separated flow was captured in a time correlated manner.

Jiang, N.↗

Dynamics and chemical mode analysis of plasma thermal-chemical instability

The stability of the weakly ionized plasma and the transition from a stable homogeneous discharge to unstable filaments play an important role in gas laser physics, plasma-assisted combustion, chemical reforming, and material synthesis. Here, theoretical stability analysis and thermal-chemical mode analysis were performed to understand the mechanism of plasma thermal-chemical instability by using a zero-dimensional plasma system with both simplified and detailed chemical kinetics of H 2 /O 2 /N 2 mixtures. The plasma dynamic and kinetic models accounted for multiple physical mechanisms in the chemically-reactive weakly ionized plasma, including ionization, attachment/detachment, recombination, vibrational and electronic energy relaxation, convective and diffusive species/heat removal, Joule heating, and detailed chemical kinetics. An analytical criterion and the explosive mode species/temperature pointers were formulated while the representative active species were identified for different thermal-chemical modes. The results showed that in addition to the classical thermal-ionization mechanism, various chemical modes from chemical heat imbalance and elementary kinetics significantly modified the time dynamics and the stability of the weakly ionized plasma. The present analysis provides insights and guidance to control plasma instability using chemical kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Nonlocal Magneto-curvature Instability in a Differentially Rotating Disk

A global mode is shown to be unstable to nonaxisymmetric perturbations in a differentially rotating Keplerian disk containing either vertical or azimuthal magnetic fields. In an unstratified cylindrical disk model, using both global eigenvalue stability analysis and linear global initial-value simulations, it is demonstrated that this instability dominates at strong magnetic fields where local standard magnetorotational instability (MRI) becomes stable. Unlike the standard MRI mode, which is concentrated in the high flow shear region, these distinct global modes (with low azimuthal mode numbers) are extended in the global domain and are Alfvén-continuum-driven unstable modes. As its mode structure and relative dominance over MRI are inherently determined by the global spatial curvature as well as the flow shear in the presence of a magnetic field, we call it the magneto-curvature (magneto-spatial-curvature) instability. Consistent with the linear analysis, as the field strength is increased in the nonlinear simulations, a transition from MRI-driven turbulence to a state dominated by global nonaxisymmetric modes is obtained. This global instability could therefore be a source of nonlinear transport in accretion disks at a higher magnetic field than predicted by local models.

79 ASTRONOMY AND ASTROPHYSICS↗

Transition in supersonic flow past axisymmetric bodies

The paper reviews some recent work on the stability of supersonic flow past axisymmetric bodies. Results indicate that tranverse curvature effect can both be stabilizing or destabilizing depending upon the particular instability modes involved. Small nose bluntness is found to have a stabilizing influence on the flow past a cone. The relevance of these results to supersonic boundary-layer transition is brought out and comparison with the experimental data is made where possible.

Malik, M. R.↗

Test results for rotordynamic coefficients of anti-swirl self-injection seals

Test results are presented for rotordynamic coefficients and leakage for three annular seals which use anti-swirl self-injection concept to yield significant improvement in whirl frequency ratios as compared to smooth and damper seals. A new anti-swirl self-inection mechanism is achieved by deliberately machining self-injection holes inside the seal stator mechanism which is used to achieve effective reduction of the tangential flow which is considered as a prime cause of rotor instability in high performance turbomachinery. Test results show that the self-injection mechanism significantly improves whirl frequency ratios; however, the leakage performance degrades due to the introduction of the self-injection mechanism. Through a series of the test program, an optimum anti-swirl self-injection seal which uses a labyrinth stator surface with anti-axial flow injections is selected to obtain a significant improvement in the whirl frequency ratio as compared to a damper seal, while showing moderate leakage performance. Best whirl frequency ratio is achieved by an anti-swirl self-injection seal of 12 holes anti-swirl and 6 degree anti-leakage injection with a labyrinth surface configuration. When compared to a damper seal, the optimum configuration outperforms the whirl frequency ratio by a factor of 2.

Kim, C. H.↗

Verification and Validation of MFiX-PIC

This work is aimed at verification and validation (V&V) of particle-in-cell (PIC) approach in the National Energy Technology Laboratory’s (NETL) Multiphase Flow with Interphase eXchanges (MFiX-PIC) which was officially released in April 2019. Under verification, the results from MFiX-PIC are compared with exact or analytical solutions, while experimental results are used for benchmarking under validation. The test suite consists of the following cases: Terminal velocity. Particle advection in time-varying flow field. Particle settling in a dense medium. Evaporation. Rayleigh-Taylor instability. Minimum fluidization. Mini-circulating fluidized bed. The cases reported in a format consistent with the existing Verification and Validation Manual (Musser et al., 2018) explore a wide range of flow velocities and concentration regimes. Only one case considered in this work verifies the implementation of heat and mass transfer routines in MFiX-PIC, while a few other studies are currently underway to further expand the validation space in reacting flows. This study highlights the need for a rigorous analysis to objectively quantify the effect of model parameters. In addition, there is a scarcity of high-quality experimental data in the literature suitable for validating PIC methodology. Future efforts would include performing experiments which provide results for a valuable database. While designing such experiments, particle count should be significant to have sufficient particles per parcel while maintaining an optimal parcel to cell volume ratio. Also, datasets should have objectively assessed experimental uncertainty whereby confidence interval in those measurements become more explicit.

42 ENGINEERING↗

Effect of Compliant Walls on Secondary Instabilities in Boundary-Layer Transition

For aerodynamic and hydrodynamic vehicles, it is highly desirable to reduce drag and noise levels. A reduction in drag leads to fuel savings. In particular for submersible vehicles, a decrease in noise levels inhibits detection. A suggested means to obtain these reduction goals is by delaying the transition from laminar to turbulent flow in external boundary layers. For hydrodynamic applications, a passive device which shows promise for transition delays is the compliant coating. In previous studies with a simple mechanical model representing the compliant wall, coatings were found that provided transition delays as predicted from the semi-empirical e(sup n) method. Those studies were concerned with the linear stage of transition where the instability of concern is referred to as the primary instability. For the flat-plate boundary layer, the Tollmien-Schlichting (TS) wave is the primary instability. In one of those studies, it was shown that three-dimensional (3-D) primary instabilities, or oblique waves, could dominate transition over the coatings considered. From the primary instability, the stretching and tilting of vorticity in the shear flow leads to a secondary instability mechanism. This has been theoretical described by Herbert based on Floquet theory. In the present study, Herbert's theory is used to predict the development of secondary instabilities over isotropic and non-isotropic compliant walls. Since oblique waves may be dominant over compliant walls, a secondary theory extention is made to allow for these 3-D primary instabilities. The effect of variations in primary amplitude, spanwise wavenumber, and Reynolds number on the secondary instabilities are examined. As in the rigid wall case, over compliant walls the subharmonic mode of secondary instability dominates for low-amplitude primary disturbances. Both isotropic and non-isotropic compliant walls lead to reduced secondary growth rates compared to the rigid wall results. For high frequencies, the non-isotropic wall suppresses the amplification of the secondary instabilities, while instabilities over the isotropic wall may grow with an explosive rate similar to the rigid wall results. For the more important lower frequencies, both isotropic and non-isotropic compliant walls suppress the amplification of secondary instabilities compared to the rigid wall results. The twofold major discovery and demonstration of the present investigation are: (1) the use of passive devices, such as compliant walls, can lead to significant reductions in the secondary instability growth rates and amplification; (2) suppressing the primary growth rates and subsequent amplification enable delays in the growth of the explosive secondary instability mechanism.

Joslin, Ronald D.↗

Boundary Layer Instabilities Over a Cone-Cylinder-Flare Model at Mach 6

Computations are performed to investigate the boundary-layer instabilities over a sharp cone-cylinder-flare model at zero degrees angle of attack. The model geometry and the flow conditions are selected to match the experiments conducted in the Boeing/AFOSR Mach 6 Quiet Tunnel (BAM6QT) at Purdue University. The geometry consists of a nominally sharp 5-degree half-angle cone, followed by a cylindrical segment and then a 10-degree flare. An axisymmetric separation bubble is generated as a result of the laminar shock/boundary-layer interaction in the cylinder-flare region. The comparison of the laminar flow solution and the schlieren images shows a remarkable agreement between the respective locations of both the boundarylayer edge and the reattachment shock. The predicted heat flux distribution is also in agreement with the measured values downstream of the reattachment location. The analysis of convective and global instabilities is performed for flare half angles equal to 8, 10, and 12 degrees and nosetip radii equal to 0.1, 1, and 5 mm. The linear amplification of first and second Mack mode instabilities that begin to amplify in the cone region are computed with a combination of the parabolized stability equations (PSE) and the harmonic linearized Navier-Stokes equations (HLNSE). The predicted frequency spectra of the surface pressure fluctuations associated with both planar and oblique instability waves are compared with the measured spectra at the various locations of the PCB and Kulite sensors. The comparison shows that the computational analysis captures the distinct lobes within the disturbance amplification spectra measured in the experiments, but some differences in amplification characteristics are noted at low frequencies. Overall, the oblique disturbances are found to be more amplified than the planar disturbances. To our knowledge, this represents the first successful comparison between convective instability analysis and measured surface pressure fluctuations for a hypersonic configuration with a separation bubble. Finally, the global instability analysis shows that the laminar flow becomes supercritical for flare half angles larger than 8 degrees. The unstable global mode for the experimental configuration of a 10 degrees flare and a sharp nosetip cone corresponds to a stationary three-dimensional disturbance that is concentrated in the recirculation region and achieves its maximum growth rate for an azimuthal wavenumber of 5.

Hydrodynamic instability↗

TF34 engine compression system computer study

The stability of the fan and the compressor components was examined individually using linearized and time dependent, one dimensional stability analysis techniques. The stability of the fan core integrated compression system was investigated using a two dimensional compression system model. The analytical equations on which this model was based satisfied the mass, axial momentum, radial momentum, and energy conservation equations for flow through a finite control volume. The results gave an accurate simulation of the flow through the compression system. The speed lines of the components were reproduced; the points of instability were accurately predicted; the locations where the instability was initiated in the fan and the core were indicated; and the variation of the bypass ratio during flow throttling was calculated. The validity of the analytical techniques was then established by comparing these results with test data and with results obtained from the steady state cycle deck.

Hosny, W. M.↗

Numerical simulation of nonlinear development of instability waves

The nonlinear interactions of high amplitude instability waves in turbulent jets are described. In plane shear layers Riley and Metcalf (1980) and Monkewitz (1987) have shown that these interactions are dependent, among other parameters, on the phase-difference between the two instability waves. Therefore, here researchers consider the nonlinear development of both the amplitudes and the phase of the instability waves. The development of these waves are also coupled with the development of the mean flow and the background turbulence. In formulating this model it is assumed that each of the flow components can be characterized by conservation equations supplemented by closure models. Results for the interactions between the two instability waves under high-amplitude forcing at fundamental and subharmonic frequencies are presented here. Qualitative agreements are found between the present predictions and available experimental data.

Mankbadi, Reda R.↗

Correlated Temporal and Spectral Variability in Neutron Star and Black Hole X-Ray Binaries

The variability of neutron star and black hole X-ray sources has several dimensions, because of the roles played by different important time-scales. The variations on time scales of hours, weeks, and months, ranging from 50% to orders of magnitude, arise out of changes in the flow in the disk. The most important driving forces for those changes are probably various possible instabilities in the disk, though there may be effects with other dominant causes. The changes in the rate of flow appear to be associated with changes in the flow's configuration, as the accreting material approaches the compact object, for there are generally correlated changes in both the X-ray spectra and the character of the faster temporal variability. There has been a lot of progress in tracking these correlations, both for Z and Atoll neutron star low-mass X-ray binaries, and for black hole binaries. I will discuss these correlations and what they tell us about the physical states of the systems.

Swank, J.↗

Correlated Temporal and Spectral Variability

The variability of neutron star and black hole X-ray sources has several dimensions, because of the roles played by different important time-scales. The variations on time scales of hours, weeks, and months, ranging from 50% to orders of magnitude, arise out of changes in the flow in the disk. The most important driving forces for those changes are probably various possible instabilities in the disk, though there may be effects with other dominant causes. The changes in the rate of flow appear to be associated with changes in the flow's configuration, as the accreting material approaches the compact object, for there are generally correlated changes in both the Xray spectra and the character of the faster temporal variability. There has been a lot of progress in tracking these correlations, both for Z and Atoll neutron star low-mass X-ray binaries, and for black hole binaries. I will discuss these correlations and review briefly what they tell us about the physical states of the systems.

Swank, Jean H.↗

Crystal growing from the melt

The mechanical and electrical properties of crystals produced by a unidirectional process depend strongly on the temperature and flow fields since these control the concentration of solute at the melt-crystal interface. The solute gradient there drives morphological instabilities that lead to cellular or dendritic interfaces. In the presentation several features of flow-solidification interactions will be discussed. These will include the effects of convection driven by density changes and buoyancy and the imposition of forced flow.

Davis, S. H.↗

Simple Sound Source in a 50- to 300- Foot-Per-Second Air Stream

In jet-engine combustor research, it may sometimes be desirable to excite some flow field within a combustor with a train of sound waves. This is because some of the flow fields are unstable, and the sound can trigger the instability. A sound source is thus used as an instrument to search out unstable flame, and as a mixer to increase heat release rates. When the combustor is small, it is possible to use electronic equipment to produce the desired sound. As the size of the combustor increases, it reaches a point where electronic drivers are no longer convenient. For such a cornbustor a sound source sustained by the approach flow is desired.

Donald R Boldman↗

Resonant instability of supersonic shear layers

A computer simulation of possible resonant instability of a supersonic shear layer is carried out. The resonance of two acoustic duct modes of the flow induced by periodic Mach waves generated by a wavy wall is sought. Results of the simulations are reported. Simulations are unable to document a resonant instability and the mixing characteristics remain unchanged. Possible weakness of the present simulations are discussed. A second set of simulations involving a mixing layer separating a supersonic and a subsonic stream were performed. A wavy wall placed adjacent to the supersonic stream to produce a set of periodic Mach waves terminating at the shear layer is modelled. The entire flow field is similar to that of an imperfectly expanded supersonic jet discharging into a subsonic coflowing stream for which enhanced mixing due to the onset of screech (feedback instability) is known to occur. The purpose of these simulations is to see if enhanced mixing and feedback instability would, indeed, take place. Some evidence of feedback oscillations is found in the simulated flow.

Tam, C. K. W.↗

Early time evolution of negative ion clouds and electron density depletions produced during electron attachment chemical release experiments

Two-dimensional electrostatic particle-in-cell simulations are used to study the early time evolution of electron depletions and negative ion clouds produced during electron attachment chemical releases in the ionosphere. The simulation model considers the evolution in the plane perpendicular to the magnetic field and a three-species plasma that contains electrons, positive ions, and also heavy negative ions that result as a by-product of the electron attachment reaction. The early time evolution (less than the negative ion cyclotron period) of the system shows that a negative charge surplus initially develops outside of the depletion boundary as the heavy negative ions move across the boundary. The electrons are initially restricted from moving into the depletion due to the magnetic field. An inhomogenous electric field develops across the boundary layer due to this charge separation. A highly sheared electron flow velocity develops in the depletion boundary due to E x B and Delta-N x B drifts that result from electron density gradients and this inhomogenous electric field. Structure eventually develops in the depletion boundary layer due to low-frequency electrostatic waves that have growth times shorter than the negative ion cyclotron period. It is proposed that these waves are most likely produced by the electron-ion hybrid instability that results from sufficiently large shears in the electron flow velocity.

Scales, W. A.↗

Aerodynamics of Stardust Sample Return Capsule

Successful return of interstellar dust and cometary material by the Stardust Sample Return Capsule requires an accurate description of the Earth entry vehicle's aerodynamics. This description must span the hypersonic-rarefied, hypersonic-continuum, supersonic, transonic, and subsonic flow regimes. Data from numerous sources are compiled to accomplish this objective. These include Direct Simulation Monte Carlo analyses, thermochemical nonequilibrium computational fluid dynamics, transonic computational fluid dynamics, existing wind tunnel data, and new wind tunnel data. Four observations are highlighted: 1) a static instability is revealed in the free-molecular and early transitional-flow regime due to aft location of the vehicle s center-of-gravity, 2) the aerodynamics across the hypersonic regime are compared with the Newtonian flow approximation and a correlation between the accuracy of the Newtonian flow assumption and the sonic line position is noted, 3) the primary effect of shape change due to ablation is shown to be a reduction in drag, and 4) a subsonic dynamic instability is revealed which will necessitate either a change in the vehicle s center-of-gravity location or the use of a stabilizing drogue parachute.

Mitcheltree, R. A.↗

Blunt-Body Paradox and Improved Application of Transient-Growth Framework

The “Reshotko-Tumin transition criterion" based on optimal transient growth successfully correlates laboratory measurements of roughness induced transition over blunt body configurations. Even though transient growth has not been conclusively linked to the measured onset of transition, the above correlation denotes the only available physics-based model for subcritical transition in blunt body flows, since the latter do not support any modal instabilities at typical experimental conditions. Unlike other established models based on empirical curve fits that are valid for a specific subclass of datasets, the optimal-growth-based transition criterion appears to provide a reasonable correlation with measurements in various wind tunnel and ballistic range facilities and for a broad range of surface temperature ratios. This paper is focused on optimal growth calculations that improve upon significant shortcomings of the computations underlying the Reshotko-Tumin correlation. The improved framework is applied to leeward transition over a spherical section forebody that was tested in the Mach 6 Adjustable Contour Expansion wind tunnel at Texas A&M University. The computed results highlight the significance of nonparallel basic state evolution, curvature terms, and an optimization procedure that varies both inflow and outflow locations of the transient growth interval. More important, the results indicate that the modified correlation is very close to its original form, and hence, that the accuracy of the transient-growth-based transition criterion is not compromised by using a more thorough theoretical framework. Yet the results also show that the optimal energy gain up to the predicted transition onset location can be rather small, highlighting the need to further investigate the optimal growth criterion for additional experimental configurations and to also uncover the in-depth physics underlying blunt body transition.

Freestream conditions↗