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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Laminar-Turbulent Transition Behind Discrete Roughness Elements in a High-Speed Boundary Layer

Computations are performed to study the flow past an isolated roughness element in a Mach 3.5, laminar, flat plate boundary layer. To determine the effects of the roughness element on the location of laminar-turbulent transition inside the boundary layer, the instability characteristics of the stationary wake behind the roughness element are investigated over a range of roughness heights. The wake flow adjacent to the spanwise plane of symmetry is characterized by a narrow region of increased boundary layer thickness. Beyond the near wake region, the centerline streak is surrounded by a pair of high-speed streaks with reduced boundary layer thickness and a secondary, outer pair of lower-speed streaks. Similar to the spanwise periodic pattern of streaks behind an array of regularly spaced roughness elements, the above wake structure persists over large distances and can sustain strong enough convective instabilities to cause an earlier onset of transition when the roughness height is sufficiently large. Time accurate computations are performed to clarify additional issues such as the role of the nearfield of the roughness element during the generation of streak instabilities, as well as to reveal selected details of their nonlinear evolution. Effects of roughness element shape on the streak amplitudes and the interactions between multiple roughness elements aligned along the flow direction are also investigated.

Choudhari, Meelan M.↗
NASA NTRS1994

Turbulent mixing noise from supersonic jets

There is now a substantial body of theoretical and experimental evidence that the dominant part of the turbulent noise of supersonic jets is generated directly by the large turbulence structures/instability waves of the jet flow. Earlier, Tam and Burton provided a description of the physical mechanism by which supersonically traveling instability waves can generate sound efficiently. They used the method of matched asymptotic expansions to construct an instability wave solution which is valid in the far field. The present work is an extension of the theory of Tam and Burton. It is argued that the instability wave spectrum of the jet may be regarded as generated by stochastic white noise excitation at the nozzle lip region. The reason why the excitation has white noise characteristics is that near the nozzle lip region the flow in the jet mixing layer has no intrinsic length and time scales. The present stochastic wave model theory of supersonic jet noise contains a single unknown multiplicative constant. Comparisons between the calculated noise directivities at selected Strouhal numbers and experimental measurements of a Mach 2 jet at different jet temperatures have been carried out. Favorable agreements are found.

Tam, Christopher K. W.↗
DOE OSTI2020

A Comparison of Interface Growth Models Applied to Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

Sophisticated numerical models that contain fluid interfaces rely upon interface evolution models to approximate the transition to turbulence near interfaces, in the presence of Rayleigh–Taylor (RTI) or Richtmyer–Meshkov (RMI) instability. Semi-analytical models have been developed in recent decades to predict the interface growth from an initial state into the nonlinear regime. Two of these models are considered in this study. They are the Goncharov and the z-models. Both of these interface models have strengths and weaknesses, which are examined here. Both of them have been implemented in the xRAGE compressible flow solver, which models fluid interfaces. The flow solver provides the fluids acceleration as a body force to the interface model. The purpose of such interface model is to evolve the early times of interface position as a subgrid model within a compressible flow simulation in order to then initialize a turbulence model. In this work, the interface models are assessed and compared for their evolution of RTI and RMI. The z-model performed better than the Goncharov model for all cases that were explored.

42 ENGINEERING↗
DOE OSTI2022

3D Printed Carbon Aerogels for Polymer-Electrolyte Fuel Cells

Carbon aerogels (CA) are known for ultra-low density, very high surface area, fine porous structure, and excellent electrical conductivity. Recent advances in high resolution 3D printing at LLNL enable the manufacturing of complex structured CA electrodes which have been demonstrated in electrolysis, flow battery, and supercapacitor applications. Polymer-electrolyte fuel cells (PEFC) may be an excellent candidate for these materials/methods as designers seek to optimize membrane-electrode assemblies to achieve higher efficiencies. A major challenge to PEFC design is two-phase flow due to water production which causes pore flooding and leads to cell instability and degradation. To mitigate these effects, groups have investigated the creation of liquid specific pathways in the porous layers to alleviate liquid buildup and enhance gas diffusion. While some progress has been made, traditionally manufactured GDLs, which consist of randomly distributed carbon fibers, are limiting. 3D printed CAs are a promising solution as complex macroporous electrode structures can be manufactured out of micro and nano porous filaments with a high level of control. Additionally, the process may allow for the consolidation of multiple layers (MPL, GDL, & flow channels) into a single printed architecture. In this work various CA lattice structures were printed using a direct ink write (DIW) method and liquid breakthrough pressures were measured to demonstrate the ability to tune the water pathways to achieve specific pressure thresholds. Furthermore, a discussion of the prospects for printed CA electrodes in PEFC will focus on material properties, synthesis techniques, and current challenges.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Nonlinear interaction between a pair of oblique modes in a supersonic mixing layer: Long-wave limit

The nonlinear interaction between a pair of symmetric, oblique, and spatial instability modes is studied in the long-wave limit using asymptotic methods. The base flow is taken to be a supersonic mixing layer whose Mach number is such that the corresponding vortex sheet is marginally stable according to Miles' criterion. It is shown that the amplitude of the mode obeys a nonlinear integro-differential equation. Numerical solutions of this equation show that, when the obliqueness angle is less than pi/4, the effect of the nonlinearity is to enhance the growth rate of the instability. The solution terminates in a singularity at a finite streamwise location. This result is reminiscent of that obtained in the vicinity of the neutral point by other authors in several different types of flows. On the other hand, when the obliqueness angle is more than pi/4, the streamwise development of the amplitude is characterized by a series of modulations. This arises from the fact that the nonlinear term in the amplitude equation may be either stabilizing or destabilizing, depending on the value of the streamwise coordinate. However, even in this case the amplitude of the disturbance increases, though not as rapidly as in the case for which the angle is less than pi/4. Quite generally then, the nonlinear interaction between two oblique modes in a supersonic mixing layer enhances the growth of the disturbance.

Balsa, Thomas F.↗
NASA NTRS1990

On the Goertler vortex instability mechanism at hypersonic speeds

The linear instability of the hypersonic boundary layer on a curved wall is considered. As a starting point real-gas effects are ignored and the viscosity of the fluid is taken to be related to the temperature either by Chapman's Law or by Sutherland's Law. It is shown that the flow is susceptible to Goertler vortices. If Chapman's Law is used the vortices are trapped in the logarithmically thin adjustment layer in which the temperature of the basic flow changes rapidly to its free stream value and the nonuniqueness of the neutral stability curve associated with incompressible Goertler vortices is shown to disappear at high Mach numbers if the appropriate 'fast' streamwise dependence of the instability is built into the disturbance flow structure. If, on the other hand, Sutherland's Law is used, the vortices are found to spread into an O(1) region and the concept of a unique neutral stability curve is not tenable because of the nonparallel effects. For both laws the leading order terms in the expansions of the Goertler number are independent of the wave number and are due to the curvature of the basic state.

Hall, P.↗
NASA NTRS2000

Phase Shift Interferometer and Growth Set Up to Step Pattern Formation During Growth From Solutions. Influence of the Oscillatory solution Flow on Stability

We have assembled an experimental setup based on Michelson interferometry with the growing crystal surface as one of the reflective surfaces. The crystallization part of the device allows optical monitoring of a face of a crystal growing at temperature stable within 0.05 C in a flow of solution of controlled direction and speed. The reference arm of the interferometer contains a liquid crystal element that allows controlled shifts of the phase of the interferograms. We employ an image-processing algorithm, which combines five images with a pi/2 phase difference between each pair of images. The images are transferred to a computer by a camera capable of capturing 60 frames per second. The device allows data collection on surface morphology and kinetics during the face layers growth over a relatively large area (approximately 4 sq. mm) in situ and in real time during growth. The estimated depth resolution of the phase shifting interferometry is approximately 50 Angstroms. The data will be analyzed in order to reveal and monitor step bunching during the growth process. The crystal chosen as a model for study in this work is KH2PO4 (KDP). This optically non-linear material is widely used in frequency doubling applications. There have been a number of studies of the kinetics of KDP crystallization that can serve as a benchmark for our investigations. However, so far, systematic quantitative characteristics of step interaction and bunching are missing. We intend to present our first quantitative results on the onset, initial stages and development of instabilities in moving step trains on vicinal crystal surfaces at varying supersaturation, flow rate, and flow direction. Behavior of a vicinal face growing from solution flowing normal to the steps and periodically changing its direction in time was considered theoretically. It was found that this oscillating flow reduces both stabilization and destabilization effects resulted from the unidirectional solution flow directed up the step stream and down the step stream. This reduction of stabilization and destabilization comes from effective mixing which entangles the phase shifts between the spatially periodic interface perturbation and the concentration wave induced by this perturbation. Numerical results and simplified mixing criterion will be discussed.

Chernov, Alex A.↗
NASA NTRS1995

Contact and pressure balance structures in two-fluid cosmic-ray hydrodynamics

The role of cosmic-ray-modified contact discontinuities and pressure balance structures in two-fluid cosmic-ray hydrodynamics in one Cartesian space dimension are investigated by means of analytic and numerical solution examples, as well as by weakly nonlinear asymptotics. The fundamental wave modes of the two-fluid cosmic-ray hydrodynamic equations in the long-wavelength limit consist of the backward and forward propagating cosmic-ray-modified sound waves, with sound speed dependent on both the cosmic-ray and thermal gas pressures; the contact discontinuity; and a pressure balance mode in which the sum ofthe cosmic ray and thermal gas pressure perturbations is zero. The pressure balance mode, like the contact discontinuity is advected with the background flow. The interaction of the pressure balance mode with the contact discontinuity is investigated by means of the method of multiple scales. The thermal gas and cosmic-ray pressure perturbations satisfy a linear diffusion equation, and entropy perturbations arising from nonisentropic initial conditions for the thermal gas are frozen into the fluid. The contact discontinuity and pressure balance eigenmodes both admit nonzero perturbations in the thermal gas, whereas the cosmic-ray-modified sound waves are isentropic. The total entropy perturbation is shared between the contact discontinuity and pressure balance eigenmodes, and examples are given in which there is a transfer of entropy between the two modes. In particular, N-wave type density disturbances are obtained which arise as a result of the entropy transfer between the two modes. A weakly nonlinear geometric optics perturbation expansion is used to study the long timescale evolution of the short-wavelength entropy wave and the thermal gas sound waves in a slowly varying, large-scale background flow. The weakly nonlinear geometric optics expansion is also used to generalize previous studies of squeezing instability for short-wavelength sound waves in the two fluid model, by including a weakly nonlinear wave steepening term that leads to shock formation, as well as the effect of long time and space dependence of the background flow. Implications of cosmic-ray-modified pressure balance structures and contact discontinuities in models of the interaction of traveling interplanetary shocks and compression and rarefraction waves with the solar wind termination shock are briefly discussed.

Webb, G. M.↗
NASA NTRS1978

AE-C observations of low-energy particles and ionospheric temperatures in the turbulent polar cusp - Evidence for the Kelvin-Helmholtz instability

Particle observations at 283 km acquired with the AE-C spacecraft during the large geomagnetic storm of May 16, 1975 indicate that the polar cusp was displaced to 71 deg invariant latitude between 1020 and 1244 MLT. Three regions of low-energy particle fluxes were determined which may be indentified with regions of field-aligned current flow in the dayside auroral zone and cusp. It is suggested that the low-energy electrons are scattered and that their pitch angles are isotropized by magnetic fluctuations associated with the Kelvin-Helmholtz instability caused by shear in the proton flow into the cusp.

Potemra, T. A.↗
NASA NTRS1981

Turbulence in sound excited jets - Measurements and theory

The mechanisms by which the turbulent structure of an axisymmetric jet is modified by the presence of an acoustic excitation are examined. A model is described in which the excitation triggers instability waves at the jet exit. As these waves propagate downstream they extract energy from the mean flow and transfer it to the random turbulence. This results in an increase of the random turbulence levels and a more rapid mixing and spreading for the jet. Models are introduced for the Reynolds stress and the 'wave-induced stress'. It is shown that at high frequencies the presence of the instability wave may reduce the random turbulence levels. Numerical calculations are presented for both the radial and axial variation in the time-averaged properties of the flow as a function of excitation conditions. The calculations are compared with measurements of fluctuating velocity and pressure in a round jet with a Reynolds number of 375,000, based on jet diameter and exit velocity.

Morris, P. J.↗
NASA NTRS1987

On the onset of three-dimensionality and time-dependence in the Goertler vortex problem

The instability of large amplitude Goertler vortices in a growing boundary layer is discussed in the fully nonlinear regime. It is shown that a three-dimensional breakdown to a flow with wavy vortex boundaries similar to that which occurs in the Taylor vortex problem takes place. However, the instability is confined to the thin shear layers which have been shown to trap the region of vortex activity. The disturbance eigenfunctions decay exponentially away from the center of these layers so that the upper and lower shear layers can support independent modes of instability. The structure of the instability, in particular its location and speed of downstream propagation, is found to be entirely consistent with recent experimental results. Furthermore, it is shown that the upper and lower layers support wavy vortex instabilities with quite different frequencies. This result is again consistent with the available experimental observations.

Hall, Philip↗
NASA NTRS1994

Simulation of Unsteady Combustion in a Ramjet Engine Using a Highly Parallel Computer

Combustion instability in ramjets is a complex phenomenon that involve nonlinear interaction between acoustic waves, vortex motion and unsteady heat release in the combustor. To numerically simulate this 3-D, transient phenomenon, enormous computer resources (time, memory and disk storage) are required. Although current generation vector supercomputers are capable of providing adequate resources for simulations of this nature, their high cost and limited availability, makes such machines less than satisfactory for routine use. The primary focus of this study is to assess the feasibility of using highly parallel computer systems as a cost-effective alternative for conducting such unsteady flow simulations. Towards this end, a large-eddy simulation model for combustion instability was implemented on the Intel iPSC/860 and a careful study was conducted to determine the benefits and the problems associated with the use of such machines for transient simulations. Details of this study along with the results obtained from the unsteady combustion simulations carried out on the iPSC/860 are discussed in this paper.

Menon, Suresh↗
NASA NTRSDate not supplied

Supersonic Traveling Crossflow Wave Characteristics in Ground and Flight Tests

This paper continues analysis of data acquired in a series of supersonic crossflow transition experiments on a 65° swept- wing model tested in both wind tunnel and flight. Flow visualization data are summarized to provide context for the swept-wing instability mechanisms of interest. Unsteady pressure measurements obtained near the wing model leading edge (LE) are studied using 2nd- and 3rd-order statistical methods to extract traveling crossflow wave characteristics from these data. Comparisons between tunnel and flight data with boundary-layer (BL) Linear Stability Theory (LST) predictions provide insights into BL transition phenomena similarities and differences observed in tests using the same swept-wing model configurations. Comparisons of unsteady pressure data results for wing LE configurations with and without distributed-roughness-element (DRE) patterns installed, suggest traveling crossflow disturbance growth is reduced when this type of flow control is applied. A brief stud y of different background surface roughness levels on the wing LE suggests the potential for another flow control approach to promote laminar flow on swept wings using streamwise-biased surface finishes.

Lewis R Owens↗
NASA NTRS2010

Control of Vortex Breakdown in Critical Swirl Regime Using Azimuthal Forcing

We finally go back to the four swirl cases and see how the flow responds to either forcing m = -1 or m = -2. On the left we see the flow forced at m = -1 We see that the PVC locks onto the applied forcing also for lower swirl number causing this high TKE at the jet center. The amplification of this instability causes VB to occur at a lower swirl number. The opposite can be seen when forcing the flow at m=-2 which is basically growing in the outer shear layer causing VB to move downstream . There is no energy at the center of the vortex showing that the precessing has been damped. The mean flow is most altered at the swirl numbers were VB is unstable.

Oberleithner, Kilian↗
NASA NTRS1989

Some Observations on Transitory Stall in Conical Diffusers

Results from an experimental investigation on the flow through conical diffusers are presented. The mean and fluctuating velocity fields are compared for three diffusers with total diffusion angles of 16, 20 and 24 degrees, in the throat Mach number (M sub t) range of 0.05 to 0.95. Each of the diffusers were 14 cm long and had a 5.08 cm inlet diameter, and the flow exited into the ambient. The boundary layer at the throat was thin with the throat diameter (D sub t) to momentum thickness (O) ratio being as high as 800 at M(sub t) = 0.4. While the 16 deg diffuser flow exited with a top-hat mean velocity profile, increasing losses due to increasing separation resulted in fuller profiles for the 20 and 24 degree cases. A detailed flow field study was conducted for the 16 deg. diffuser. The u'-spectrum, measured at the exit plane, exhibited a peak apparently due to the ensuing jet column instability throughout the M(sub t) range covered. In addition, a much lower frequency spectral peak also occurred in the M(sub t) range of 0.3 to 0.7. Both of the spectral peaks were due to axisymmetric flow fluctuations. A self-sustaining flow oscillation occurred in the M(sub t) range of 0.6 to 0.85, emitting a loud tone, when the jet column instability frequency matched the resonance frequency of the diffuser. Limited data showed that artificial acoustic excitation was effective in reducing the flow fluctuations, with a resultant increase in the pressure recovery, at low M(sub t).

Zaman, K. B. M. Q.↗
NASA NTRS1990

Some observations on transitory stall in conical diffusers

Results from an experimental investigation on the flow through conical diffusers are presented. The mean and fluctuating velocity fields are compared for three diffusers with total diffusion angles of 16, 20 and 24 degrees, in the throat Mach number (M sub t) range of 0.05 to 0.95. Each of the diffusers were 14 cm long and had a 5.08 cm inlet diameter, and the flow exited into the ambient. The boundary layer at the throat was thin with the throat diameter (D sub t) to momentum thickness (O) ratio being as high as 800 at M(sub t) = 0.4. While the 16 deg diffuser flow exited with a top-hat mean velocity profile, increasing losses due to increasing separation resulted in fuller profiles for the 20 and 24 degrees cases. A detailed flow field study was conducted for the 16 deg diffuser. The u'-spectrum, measured at the exit plane, exhibited a peak apparently due to the ensuing jet column instability throughout the M(sub t) range covered. In addition, a much lower frequency spectral peak also occurred in the M(sub t) range of 0.3 to 0.7. Both of the spectral peaks were due to axisymmetric flow fluctuations. A self-sustaining flow oscillation occurred in the M(sub t) range of 0.6 to 0.85, emitting a loud tone, when the jet column instability frequency matched the resonance frequency of the diffuser. Limited data showed that artificial acoustic excitation was effective in reducing the flow fluctuations, with a resultant increase in the pressure recovery, at low M(sub t).

Zaman, K. B. M. Q.↗
NASA NTRS1973

Ion streaming instabilities with application to collisionless shock wave structure

The electromagnetic dispersion relation for two counterstreaming ion beams of arbitrary relative strength flowing parallel to a dc magnetic field is derived. The beams flow through a stationary electron background and the dispersion relation in the fluid approximation is unaffected by the electron thermal pressure. Magnetic effects on the ion beams are included, but the electrons are treated as a magnetized fluid. The dispersion relation is solved with a zero net current condition applied and the regions of instability in the k-U space (U is the relative velocity between the two ion beams) are presented. These results are extensions of Kovner's analysis for weak beams. The parameters are then chosen to be applicable for parallel shocks. It is found that unstable waves with zero group velocity in the shock frame can exist near the leading edge of the shock for upstream Alfven Mach numbers greater than 5.5.

Golden, K. I.↗
NASA NTRSDate not supplied

Aerodynamic Design Optimization for Natural Laminar Flow Airfoils

Natural laminar flow technology is a passive laminar flow control (LFC) strategy that seeks to delay the onset of boundary-layer transition (BLT) through shape optimization to reduce the drag of the aerodynamic vehicle. Adjoint-based design optimization for LFC is proposed in an integrated multidisciplinary framework, which includes the computational fluid dynamics (CFD), geometry and grid deformation, and linear stability analysis (LSA) for transition prediction. In particular, the BLT location is predicted using the dual N-factor method that is based on a linear stability theory (LST) eigenvalue problem. The dual N-factor criterion accounts for the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) boundary-layer instabilities to predict the transition location in three-dimensional boundary-layer flows. The adjoint-based shape optimization procedure is based on an iteratively coupled CFD and LSA methodology to converge the transition location and flow solutions, as well as to calculate the sensitivities of the aerodynamic metrics of interest with respect to the flow and shape design parameters. The RAE 2822 airfoil at 0 and 30 degrees yaw angles, an angle of attack of 0.72 degrees, and subsonic conditions (M∞ = 0.19, Rec = 5.6 × 106 ) are used as baseline configurations for design optimization. The angle of attack and the vertical displacement of free-form-deformation control points are used as design variables to reduce the drag coefficient while reaching a specified lift coefficient. The optimized unswept airfoil designs achieve a 30% drag reduction accompanied by a downstream shift of the transition locations over both suction and pressure sides of the airfoil. The initial design iterations for the swept case also show a favorable trend in the drag reduction with transition delay over both sides.

Transition↗