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 379 records · Page 21

Linear Disturbance Amplification Over Blunted Flat Plates in High-Speed Flows

Modal and nonmodal instability characteristics of cylindrically blunted flat plates with varying leading edge radii are described for Mach 4 and Mach 6 freestream conditions. The selection of leading edge radii and freestream parameters is informed by experimental conditions. The investigation of this 2D problem provides a slow entropy layer swallowing which allows for an isolated development of perturbations seeded upstream within different wall-normal regions of the flow. At both Mach numbers, a decrease in modal instability amplification was seen as the leading edge radius was increased. Nonmodal analysis reveals amplifying perturbations in the boundary layer as well as in the entropy layer. The medium bluntness regime exhibits the strongest amplification of nonmodal disturbances that is nonmonotonic in character. Small amplitude boundary forcing at the plate surface or volumetric forcing at various wall-normal heights was used to account for receptivity effects. While wall forcing effectively induced modal instabilities, only an actuation above the boundary layer captured disturbances that amplify within the entropy layer. The optimal nonmodal theory’s entropy-layer disturbance evolution exhibited outstanding agreement with controlled forcing, including receptivity effects. The evolution of entropy-layer disturbances from the optimal nonmodal theory showed excellent agreement with the results of receptivity to controlled forcing. Therefore, the nonmodal optimal growth analysis may provide a useful as well as efficient technique to identify the complete disturbance spectrum in blunt hypersonic configurations, where both modal and nonmodal disturbances can amplify in the boundary-layer and entropy-layer regions.

Boundary layer transition↗

Boiler loop instabilities.

Closed loop low pressure forced flow convective boiling system analyzed for instabilities due to feed system/boiler interaction

Grace, T. M.↗

Wave Number Selection for Incompressible Parallel Jet Flows Periodic in Space

The temporal instability of a spatially periodic parallel flow of an incompressible inviscid fluid for various jet velocity profiles is studied numerically using Floquet Analysis. The transition matrix at the end of a period is evaluated by direct numerical integration. For verification, a method based on approximating a continuous function by a series of step functions was used. Unstable solutions were found only over a limited range of wave numbers and have a band type structure. The results obtained are analogous to the behavior observed in systems exhibiting complexity at the edge of order and chaos.

Miles, Jeffrey Hilton↗

Deleterious Thermal Effects due to Randomized Flow Paths in Pebble Bed, and Particle Bed Style Reactors

Reactor fuel rod surface area that is perpendicular to coolant flow direction (+S) i.e. perpendicular to the P creates areas of coolant stagnation leading to increased coolant temperatures resulting in localized changes in fluid properties. Changes in coolant fluid properties caused by minor increases in temperature lead to localized reductions in coolant mass flow rates leading to localized thermal instabilities. Reductions in coolant mass flow rates result in further increases in local temperatures exacerbating changes to coolant fluid properties leading to localized thermal runaway. Unchecked localized thermal runaway leads to localized fuel melting. Reactor designs with randomized flow paths are vulnerable to localized thermal instabilities, localized thermal runaway, and localized fuel melting.

Moran, Robert P.↗

Linear Disturbance Amplification Over Blunted Flat Plates in High-Speed Flows

Modal and nonmodal instability characteristics of cylindrically blunted flat plates with varying leading edge radii are described for Mach 4 and Mach 6 freestream conditions. The selection of leading edge radii and freestream parameters is informed by experimental conditions. The investigation of this 2D problem provides a slow entropy layer swallowing which allows for an isolated development of perturbations seeded upstream within different wall-normal regions of the flow. At both Mach numbers, a decrease in modal instability amplification was seen as the leading edge radius was increased. Nonmodal analysis reveals amplifying perturbations in the boundary layer as well as in the entropy layer. The medium bluntness regime exhibits the strongest amplification of nonmodal disturbances that is nonmonotonic in character. Small amplitude boundary forcing at the plate surface or volumetric forcing at various wall-normal heights was used to account for receptivity effects. While wall forcing effectively induced modal instabilities, only an actuation above the boundary layer captured disturbances that amplify within the entropy layer. The optimal nonmodal theory’s entropy-layer disturbance evolution exhibited outstanding agreement with controlled forcing, including receptivity effects. The evolution of entropy-layer disturbances from the optimal nonmodal theory showed excellent agreement with the results of receptivity to controlled forcing. Therefore, the nonmodal optimal growth analysis may provide a useful as well as efficient technique to identify the complete disturbance spectrum in blunt hypersonic configurations, where both modal and nonmodal disturbances can amplify in the boundary-layer and entropy-layer regions.

Boundary layer transition↗

Nonlinear instability of accelerating shock waves with application to supernovae

We consider the stability of a planar accelerating shock front in an exponential atmosphere, a situation that can be described by a self-similar solution. A previous linear instability analysis showed that there are three regimes, depending on the wavelength of the perturbation along the direction of the shock front: at long wavelengths the shock is unstable, at intermediate wavelengths it is overstable, and at short wavelenghts it is stable, where the characteristic length is set by the initial density scale height. We have carried out numerical simulations of such an accelerating shock front and have confirmed the results of the semianalytic linear analysis in the linear regime. In the nonlinear regime, the evolution again depends on the three wavelength domains. At long wavelengths, the instability continuously grows but the flow remains smooth, and at short wavelengths, the flow is stable. The intermediate wavelength, overstable regime shows more complex evolution. The growing linear oscillations saturate soon after entering the nonlinear regime and the oscillations continue with the same approximate period. The shock front develops moving points of intersection, which generate weak shock fronts and density and pressure structure in the immediate postshock flow. The density fluctuations, with a contrast of a factor of 2-3, become frozen into the downstream flow. The growth of the instability is slow, so substantial initial perturbations are needed; thesemay be present in Type II supernovae with red supergiant progenitor stars which have outer convective regions. The clumping in the outer supernova atmosphere may affect spectral line formation and may play a role in the formation of fast knots.

Luo, Ding↗

Advanced stability analysis for laminar flow control

Five classes of problems are addressed: (1) the extension of the SALLY stability analysis code to the full eighth order compressible stability equations for three dimensional boundary layer; (2) a comparison of methods for prediction of transition using SALLY for incompressible flows; (3) a study of instability and transition in rotating disk flows in which the effects of Coriolis forces and streamline curvature are included; (4) a new linear three dimensional instability mechanism that predicts Reynolds numbers for transition to turbulence in planar shear flows in good agreement with experiment; and (5) a study of the stability of finite amplitude disturbances in axisymmetric pipe flow showing the stability of this flow to all nonlinear axisymmetric disturbances.

Orszag, S. A.↗

Instabilities of dynamic thermocapillary liquid layers. I Convective instabilities

A steady shear flow is set up in a planar liquid layer bounded below by a rigid plate, and above by an interface with a passive gas, by imposing a temperature gradient along the layer and driving the motion by thermocapillarity. The thermal-convective instabilities to which this dynamic state is susceptible include stationary longitudinal rolls involving the classical Marangoni instability, and unsteady hydrothermal waves, whose new instability mechanism derives its energy from the horizontal temperature gradients. The thermal stability characteristics for liquid layers with and without return flow profiles are presently given as functions of the Prandtl number for the liquid and the Biot number for the interface, and comparisons are made with available experimental observations.

Smith, M. K.↗

Stability of an oscillated fluid with a uniform density gradient

Instabilities in a fluid with a constant density gradient that is subject to arbitrarily oriented oscillatory accelerations are considered. With the Boussinesq approximation and for the case of an unbounded fluid, transformation to Lagrangian coordinates allows the reduction of the problem to an ordinary differential equation for each three-dimensional wavenumber. The problem has three parameters: the nondimensional amplitude R of the base-state oscillation, the nondimensional level of background steady acceleration, which for some cases can be represented in terms of a local (in time) Richardson number Ri, and the Prandtl number Pr. Some general bounds on stability are derived. For Pr = 1 closed-form solutions are found for impulse (delta function) accelerations and a general asymptotic solution is constructed for large R and general imposed accelerations. The asymptotic solution takes advantage of the fact that at large R wave growth is concentrated at 'zero points'. These are times when the effective vertical wavenumber passes through zero. Kelvin-Helmholtz instabilities are found to dominate at low R, while Rayleigh-Taylor instabilities dominate at high R. At high R, the uniform shear of the Kelvin-Helmholtz case tends to distort and weaken instability waves. With unsteady flows, Ri = 1/4 is no longer an instability limit. Significant instabilities have been found for sinusoidal forcing for Ri up to 0.6.

Jacqmin, David↗

Entry Dispersion Analysis for the Stardust Comet Sample Return Capsule

Stardust will be the first mission to return samples from beyond the Earth-Moon system. The sample return capsule, which is passively controlled during the fastest Earth entry ever, will land by parachute in Utah. The present study analyzes the entry, descent, and landing of the returning sample capsule. The effects of two aerodynamic instabilities are revealed (one in the high altitude free molecular regime and the other in the transonic/subsonic flow regime). These instabilities could lead to unacceptably large excursions in the angle-of-attack near peak heating and main parachute deployment, respectively. To reduce the excursions resulting from the high altitude instability, the entry spin rate of the capsule is increased. To stabilize the excursions from the transonic/subsonic instability, a drogue chute with deployment triggered by an accelerometer and timer is added prior to main parachute deployment. A Monte Carlo dispersion analysis of the modified entry (from which the impact of off-nominal conditions during the entry is ascertained) shows that the capsule attitude excursions near peak heating and drogue chute deployment are within Stardust program limits. Additionally, the size of the resulting 3-sigma landing ellipse is 83.5 km in downrange by 29.2 km in crossrange, which is within the Utah Test and Training Range boundaries.

Desai, Prasun N.↗

Entry Trajectory Issues for the Stardust Sample Return Capsule

The Stardust mission was successfully launched on February 7, 1999. It will be the first mission to return samples from a comet. The sample return capsule, which is passively controlled during the fastest Earth entry ever, will land by parachute in Utah. The present study describes the analysis of the entry, descent, and landing of the returning sample capsule utilizing the final, launch configuration capsule mass properties. The effects of two aerodynamic instabilities are revealed (one in the high altitude free molecular regime and the other in the transonic/subsonic flow regime). These instabilities could lead to unacceptably large excursions in the angle-of-attack near peak heating and main parachute deployment, respectively. To reduce the excursions resulting from the high altitude instability, the entry spin rate of the capsule is increased. To stabilize the excursions from the transonic/subsonic instability, a drogue chute with deployment triggered by a gravity-switch and timer is added prior to main parachute deployment. A Monte Carlo dispersion analysis of the modified entry (from which the impact of off-nominal conditions during the entry is ascertained) predicts that the capsule attitude excursions near peak heating and drogue chute deployment are within Stardust mission limits. Additionally, the size of the resulting 3-sigma landing ellipse is 60.8 km in downrange by 19.9 km in crossrange, which is within the Utah Test and Training Range boundaries.

Desai, Prasun N.↗

Entry Trajectory Issues for the Stardust Sample Return Capsule

The Stardust mission was successfully launched on February 7, 1999. It will be the first mission to return samples from a comet. The sample return capsule, which is passively controlled during the fastest Earth entry ever, will land by parachute in Utah. The present study describes the analysis of the entry, descent, and landing of the returning sample capsule utilizing the final, launch configuration capsule mass properties. The effects of two aerodynamic instabilities are revealed (one in the high altitude free molecular regime and the other in the transonic/subsonic flow regime). These instabilities could lead to unacceptably large excursions in the angle-of-attack near peak heating and main parachute deployment, respectively. To reduce the excursions resulting from the high altitude instability, the entry spin rate of the capsule is increased. To stabilize the excursions from the transonic/subsonic instability, a drogue chute with deployment triggered by a gravity-switch and timer is added prior to main parachute deployment. A Monte Carlo dispersion analysis of the modified entry (from which the impact of off-nominal conditions during the entry is ascertained) predicts that the capsule attitude excursions near peak heating and drogue chute deployment are within Stardust mission limits. Additionally, the size of the resulting 3-sigma landing ellipse is 60.8 km in downrange by 19.9 km in crossrange, which is within the Utah Test and Training Range boundaries.

Desai, Prasun N.↗

Morphological stability in the presence of fluid flow in the melt

Recent experiments have shown that the presence of a vertical buoyancy-driven flow adjacent to an initially cylindrical crystal-melt interface may produce a time-dependent helical deformation of the interface, with a rotation period ranging from several minutes to many hours, depending upon the width of the melt. The temperature distribution is such that the interface is expected to be morphologically stable in the absence of fluid flow. A linear stability analysis reveals that the instability is due to a coupling between a basic hydrodynamic instability in the buoyant flow and the deformable boundary separating the two phases. The crystal-melt interface lowers the critical Grashof number of an analogous rigid-walled system by an order of magnitude for succinonitrile with a Prandtl number P = 22.8; furthermore, the hydrodynamic mode that is actually destabilized by the interface is not the least stable mode in the rigid-walled system for P = 22.8. The results show that the instability may be regarded either as a rather large alteration of a basic hydrodynamic instability by the crystal-melt interface, or as a significant modification of the morphological stability of the interface by the presence of the buoyant flow.

Mcfadden, G. B.↗

Turbulence modeling in aerodynamic shear flows - Status and problems

This paper briefly summarizes the status and problems of turbulence modeling for aerodynamical applications. For complex flows the 'approach of choice' is (increasingly) full second-order (Reynolds stress equation) closure. These closures have not yet developed to anywhere near their full potential, significant further research is required especially regarding length-scale equations, representation of pressure-strain correlations, and wall region treatments. Recent developments in computer capability, algorithms, numerical simulations, theory and quantitative flow visualization should assist in and hasten this research. Several problem areas such as shock interaction and discrete dynamic instabilities of turbulent flows may require mega-to-large eddy simulation or theoretical adjuncts.

Bushnell, D. M.↗

On the thermal instability of galactic and cluster halos

The paper presents a detailed study of thermal instabilities in cooling flows associated with galaxies and clusters of galaxies. In the case of purely radiation-driven accretion onto a central object such as the cD galaxy M87, it is found that the gas is largely subject to overstability, rather than to monotonic instability. If thermal conductivity is taken into account, the flow is stabilized on scales of several kiloparsecs, even if the conductivity is appreciably reduced (e.g., about 1 percent) with respect to the Spitzer value. In no case are the globular perturbations (i.e., perturbations with comparable radial and azimuthal dimensions) found to be monotonically unstable. The paper presents numerical solutions of the local dispersion relation for the cooling flow in M87 and discusses the possible consequences of the results for a correct understanding of cooling flows.

Malagoli, A.↗

Direct numerical simulation of hot jets

The ultimate motivation of this work is to investigate the stability of two dimensional heated jets and its implications for aerodynamic sound generation from data obtained with direct numerical simulations (DNS). As pointed out in our last report, these flows undergo two types of instabilities, convective or absolute, depending on their temperature. We also described the limits of earlier experimental and theoretical studies and explained why a numerical investigation could give us new insight into the physics of these instabilities. The aeroacoustical interest of these flows was also underlined. In order to reach this goal, we first need to succeed in the DNS of heated jets. Our past efforts have been focused on this issue which encountered several difficulties. Our numerical difficulties are directly related to the physical problem we want to investigate since these absolutely or almost absolutely unstable flows are by definition very sensitive to the smallest disturbances and are very likely to reach nonlinear saturation through a numerical feedback mechanism. As a result, it is very difficult to compute a steady laminar solution using a spatial DNS. A steady state was reached only for strongly co-flowed jets, but these flows are almost equivalent to two independent mixing layers. Thus they are far from absolute instability and have much lower growth rates.

Jacob, Marc C.↗

Multi-Element Unstructured Analyses of Complex Valve Systems

The safe and reliable operation of high pressure test stands for rocket engine and component testing places an increased emphasis on the performance of control valves and flow metering devices. In this paper, we will present a series of high fidelity computational analyses of systems ranging from cryogenic control valves and pressure regulator systems to cavitating venturis that are used to support rocket engine and component testing at NASA Stennis Space Center. A generalized multi-element framework with sub-models for grid adaption, grid movement and multi-phase flow dynamics has been used to carry out the simulations. Such a framework provides the flexibility of resolving the structural and functional complexities that are typically associated with valve-based high pressure feed systems and have been difficult to deal with traditional CFD methods. Our simulations revealed a rich variety of flow phenomena such as secondary flow patterns, hydrodynamic instabilities, fluctuating vapor pockets etc. In the paper, we will discuss performance losses related to cryogenic control valves, and provide insight into the physics of the dominant multi-phase fluid transport phenomena that are responsible for the choking like behavior in cryogenic control elements. Additionally, we will provide detailed analyses of the modal instability that is observed in the operation of the dome pressure regulator valve. Such instabilities are usually not localized and manifest themselves as a system wide phenomena leading to an undesirable chatter at high flow conditions.

Sulyma, Peter↗