Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow instabilities”

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 541 records · Page 30

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.↗

Emergence of Unstable Focused Flow Induced by Variable‐Density Flows in Vertical Fractures

Abstract Fluids with different densities often coexist in subsurface fractures and lead to variable‐density flows that control subsurface processes such as seawater intrusion, contaminant transport, and geologic carbon sequestration. In nature, fractures have dip angles relative to gravity, and density effects are maximized in vertical fractures. However, most studies on flow and transport through fractures are often limited to horizontal fractures. Here, we study the mixing and transport of variable‐density fluids in vertical fractures by combining three‐dimensional (3D) pore‐scale numerical simulations and visual laboratory experiments. Two miscible fluids with different densities are injected through two inlets at the bottom of a fracture and exit from an outlet at the top of the fracture. Laboratory experiments show the emergence of an unstable focused flow path, which we term a “runlet.” We successfully reproduce the unstable runlet using 3D numerical simulations and elucidate the underlying mechanisms triggering the runlet. Dimensionless number analysis shows that the runlet instability arises due to the Rayleigh‐Taylor instability (RTI), and flow topology analysis is applied to identify 3D vortices that are caused by the RTI. Even under laminar flow regimes, fluid inertia is shown to control the runlet instability by affecting the size and movement of vortices. Finally, we confirm the emergence of a runlet in rough‐walled fractures. Since a runlet dramatically affects fluid distribution, residence time, and mixing, the findings in this study have direct implications for the management of groundwater resources and subsurface applications.

Environmental Sciences & Ecology↗

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.↗

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗

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↗

Control of internal transport barriers in magnetically confined tokamak burning plasmas

The initiation, termination, and control of internal transport barriers associated with E × B flow shear near local minima of magnetic shear are examined for burning plasmas to determine if the positive feedback loops between profiles, instability, transport, and flow shear operate in regimes with fusion self-heating. A five-field transport model for the evolution of profiles of density, ion and electron temperature, ion and electron fluctuations, and radial electric field is utilized to examine the efficacy of controls associated with external inputs of heat and particles, including neutral beam injection, RF, pellets, and gas puffing. The response of the plasma to these inputs is studied in the presence of self-heating. The latter is affected by the external inputs and their modification of profiles and is, therefore, not an external control. Provided sufficient external power is applied, internal transport barriers can be created and controlled, both in ion and electron channels. Barrier control is sensitive to the locations of power deposition and pellet ablation, as well as temporal sequencing of external inputs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Traffic Control via Connected and Automated Vehicles (CAVs): An Open-Road Field Experiment with 100 CAVs

The CIRCLES project aims to reduce instabilities in traffic flow, which are naturally occurring phenomena due to human driving behavior. Also called “phantom jams” or “stop-and-go waves,” these instabilities are a significant source of wasted energy. Toward this goal, the CIRCLES project designed a control system, referred to as the MegaController by the CIRCLES team, that could be deployed in real traffic. Our field experiment, the MegaVanderTest (MVT), leveraged a heterogeneous fleet of 100 longitudinally controlled vehicles as Lagrangian traffic actuators, each of which ran a controller with the architecture described in this article. The MegaController is a hierarchical control architecture that consists of two main layers. The upper layer is called the Speed Planner and is a centralized optimal control algorithm. It assigns speed targets to the vehicles, conveyed through the LTE cellular network. The lower layer is a control layer, running on each vehicle. It performs local actuation by overriding the stock adaptive cruise controller, using the stock onboard sensors. The Speed Planner ingests live data feeds provided by third parties as well as data from our own control vehicles and uses both to perform the speed assignment. The architecture of the Speed Planner allows for the modular use of standard control techniques, such as optimal control, model predictive control (MPC), kernel methods, and others. The architecture of the local controller allows for the flexible implementation of local controllers. Corresponding techniques include deep reinforcement learning (RL), MPC, and explicit controllers. Depending on the vehicle architecture, all onboard sensing data can be accessed by the local controllers or only some. Likewise, control inputs vary across different automakers, with inputs ranging from torque or acceleration requests for some cars to electronic selection of adaptive cruise control (ACC) setpoints in others. The proposed architecture technically allows for the combination of all possible settings proposed previously, that is {Speed Planner algorithms} × {local Vehicle Controller algorithms} × {full or partial sensing} × {torque or speed control}. As a result, most configurations were tested throughout the ramp up to the MegaVandertest (MVT).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Theory of instability and transition

The strongly nonlinear area of theory is discussed, as well as linear and weakly nonlinear (i.e., perturbation) theories, and it is noted that the weaknesses of the weakly nonlinear theory are essentially the inappropriate formulation in earlier works and the lack of guidance for the choice of the lowest-order basis. Attention is focused on the areas of theoretical/numerical development contributing to understanding the transition mechanism and new means for analyzing and predicting transition quantitatively. The nonlinear stability of nonparallel flows, linear secondary instability, and nonlinear wave interaction are analyzed. The incompressible flow over a flat plate with zero-pressure gradient is chosen as an example, while applications range to other shear flows including three-dimensional and compressible boundary layers.

Herbert, Thorwald↗

Microgravity Liquid-Gas Two-Phase Flow: Review of Pressure Drop and Heat Transfer Correlations and Guidelines for Equipment Operability

In this report we have catalogued the flow regimes observed in microgravity, summarized correlations for the pressure drop and rate of heat transfer that are commonly used, and discuss the validation of a few correlations from available experimental results. Two-phase flow through some specific components such as bends, tees, filters and pumps are discussed from a physical perspective to guide the designer on how reduced gravity might affect their performance. Phase separation in zero gravity is addressed through the behavior and basic design concepts for devices based on passive centrifugal action, capillary forces, gas extraction through a membrane installed in a channel wall and the use of a syringe with a perforated piston to remove bubbles from small liquid volumes. We address the common instabilities that develop in flow loops owing exclusively to the two-phase nature of the flow, e.g., Ledinegg instability and concentration waves. Finally we briefly review flow metering and gauging; two-phase flow through porous media, where pressure drop and flow regime map correlations in zero-g are a current research topic; and basic operation principles of heat pipes and capillary pumped loops.

liquid-gas↗