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,477 records · Page 82

Some Effects of Leading-Edge Sweep on Boundary-Layer Transition at Supersonic Speeds

The effects of crossflow and shock strength on transition of the laminar boundary layer behind a swept leading edge have been investigated analytically and with the aid of available experimental data. An approximate method of determining the crossflow Reynolds number on a leading edge of circular cross section at supersonic speeds is presented. The applicability of the critical crossflow criterion described by Owen and Randall for transition on swept wings in subsonic flow was examined for the case of supersonic flow over swept circular cylinders. A wide range of applicability of the subsonic critical values is indicated. The corresponding magnitude of crossflow velocity necessary to cause instability on the surface of a swept wing at supersonic speeds was also calculated and found to be small. The effects of shock strength on transition caused by Tollmien-Schlichting type of instability are discussed briefly. Changes in local Reynolds number, due to shock strength, were found analytically to have considerably more effect on transition caused by Tollmien-Schlichting instability than on transition caused by crossflow instability. Changes in the mechanism controlling transition from Tollmien-Schlichting instability to crossflow instability were found to be possible as a wing is swept back and to result in large reductions in the length of laminar flow.

Chapman, Gray T.↗

Secondary instability mechanisms in compressible, axisymmetric boundary layers

Secondary instability mechanisms in compressible, axisymmetric boundary layers are analyzed using spectrally-accurate mean-flow and stability codes. Results show that subharmonic disturbances are the most dangerous secondary disturbances in an environment with a low to moderate intensity of the primary disturbance. The relation between spatial and temporal analyses of the secondary disturbance is explored at Mach 1.6 along a flat plate and Mach 6.8 along a cone. Spatial direct numerical simulations are utilized to confirm the quantitative predictions from spatial secondary instability theory.

Ng, Lian L.↗

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, the viscosity of the fluid is taken to be a linear function of temperature and real-gas effects are ignored. It is shown that the flow is susceptible to Goertler vortices and that they are trapped in the logarithmically thin adjustment layer in which the temperature of the basic flow changes rapidly to its free stream value. The vortices decay exponentially in both directions away from this layer and are most unstable when their wavelength is comparable with the depth of the adjustment layer. 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. It is shown that in the hypersonic limit, wall-cooling affects the leading-order term in the expansion of the Goertler number which is independent of the wavenumber, and which is due to the curvature of the basic state but it has a negligible effect on the other terms which are dependent on the wavenumber.

Hall, Philip↗

Nonlinear equilibration of localized instabilities on a baroclinic jet

Dynamical mechanisms underlying the equilibration of absolute instability are examined in a nonlinear, quasigeostrophic, two-layer model. The key to understanding the nonlinear equilibration is in recognizing that linear absolute instabilities can be stabilized both by a reduction of the vertical shear and by enhancement of the mean barotropic velocity. In a localized domain, the equilibration process proceeds with the creation of locally convectively unstable regions downstream, which encroach onto the locally absolutely unstable region until the local instability is suppressed. That local instabilities exist only if absolutely unstable regions span a minimum size is verified by eigenvalue calculations of three-dimensional flows. Numerical examples suggest that this critical size is at least 9000 km for a wide range of parameter values chosen to investigate the midlatitude storm tracks. Fluctuations arising from local absolute instability obtain maximum amplitude in the downstream convectively unstable regions rather than in the absolutely unstable regions themselves. Together, these results suggest that if an equilibrated absolute instability were to occur in midlatitudes, a zonal band of surface easterlies exceeding 900 km would be required and the associated enhanced variances would not be found coincident with the regions of absolute instability.

Delsole, Timothy M.↗

Hypersonic panel flutter in a rarefied atmosphere

Panel flutter is a form of dynamic aeroelastic instability resulting from the interaction between motion of an aircraft structural panel and the aerodynamic loads exerted on that panel by air flowing past one of the faces. It differs from lifting surface flutter in the sense that it is not usually catastrophic, the panel's motion being limited by nonlinear membrane stresses produced by the transverse displacement. Above some critical airflow condition, the linear instability grows to a limit cycle . The present investigation studies panel flutter in an aerodynamic regime known as 'free molecule flow', wherein intermolecular collisions can be neglected and loads are caused by interactions between individual molecules and the bounding surface. After collision with the panel, molecules may be reflected specularly or reemitted in diffuse fashion. Two parameters characterize this process: the 'momentum accommodation coefficient', which is the fraction of the specularly reflected molecules; and the ratio between the panel temperature and that of the free airstream. This model is relevant to the case of hypersonic flight vehicles traveling at very high altitudes and especially for panels oriented parallel to the airstream or in the vehicle's lee. Under these conditions the aerodynamic shear stress turns out to be considerably larger than the surface pressures, and shear effects must be included in the model. This is accomplished by means of distributed longitudinal and bending loads. The former can cause the panel to buckle. In the example of a simply-supported panel, it turns out that the second mode of free vibration tends to dominate the flutter solution, which is carried out by a Galerkin analysis. Several parametric studies are presented. They include the effects of (1) temperature ratio; (2) momentum accommodation coefficient; (3) spring parameters, which are associated with how the panel is connected to adjacent structures; (4) a parameter which relates compressive end load to its value which would cause classical column buckling; (5) a parameter proportional to the pressure differential between the front and back faces; and (6) initial curvature. The research is completed by an investigation into the possibility of accounting for molecular collisions, which proves to be infeasible given the speeds of current mainframe supercomputers.

Resende, Hugo B.↗

Aeroelastic Modeling and Full-Scale Loads Measurements for Investigation of Single-Axis PV Tracker Wind-Driven Dynamic Instabilities

While wind tunnel testing and proprietary industry modeling tools have been used for years to design and develop PV tracker systems, wind induced dynamic failures are becoming more prevalent and high visibility. Designers and manufactures have reacted by developing new systems or add-on products to improve system dynamics adding to overall system costs. To better understand the physics and aeroelastic behavior that leads to dynamic instabilities NREL researches have taken up a first of a kind study to both develop open source aeroelastic modeling tools and measure tracker loads on a single axis full-scale tracker system at a high wind site. These modeling tools can simulate the fluid-structure interaction driving torsional instabilities under a wide range of turbulent inflow conditions and stow angles. The results of these simulations reveal the flow features associated with panel rotation and the induced loads in the structure. The experimental load measurements can be used to validate models and provide a quantifiable understanding of tracker dynamics, critical loads paths, identify instability markers, inform resilient design, and suggest the most favorable stow approach.

14 SOLAR ENERGY↗

Hydrodynamic Dryout in Two-Phase Flows: Observations of Low Bond Number Systems

Dryout occurs readily in certain slug and annular two-phase flows for systems that exhibit partial wetting. The mechanism for the ultimate rupture of the film is attributed to van der Waals forces, but the pace towards rupture is quickened by the surface tension instability (Rayleigh-type) of the annular film left by the advancing slug and by the many perturbations of the free surface present in the Re(sub g) approximately 0(10(exp 3)), Re(sub l) approximately 0(10(exp 4)), and Ca approximately 0(10(exp -1) flows. Results from low-gravity experiments using three different test fluids are presented and discussed. For the range of tests conducted, the effect of increasing viscosity is shown to eliminate the film rupture while the decrease of surface tension via a surfactant additive is shown to dramatically enhance it. Laboratory measurements using capillary tubes are presented which reveal the sensitivity of the dryout phenomena to particulate and surfactant contamination. Rom such observations, dryout due to the hydrodynamic-van der Waals instability can be expected in a certain range of flow parameters in the absence of heat transfer. The addition of heat transfer may only exacerbate the problem by producing thermal transport lines replete with "hot spots." A caution to this effect is issued to future space systems designers concerning the use of partially wetting working fluids.

Weislogel, Mark M.↗

Low-Density Aerodynamics of the Stardust Sample Return Capsule

The aerodynamics of the Stardust Sample Return Capsule are analyzed in the low- density, transitional flow regime using free-molecular, Direct Simulation Monte Carlo, Navier-Stokes, and Newtonian methods to provide inputs for constructing a transitional flow bridging relation. The accuracy of this bridging relation in reconstructing the aero- dynamic coefficients given by the more exact methods is presented for a range of flight conditions and vehicle attitudes. There is good agreement between the various prediction methods, and a simple sine-squared bridging relation is shown to provide a reasonably good description of the axial force, normal force, and pitching moment over a range of Knudsen numbers from 0.001 to 10. The predictions show a static instability of the Star- dust capsule in the free-molecular regime that persists well into the transitional flow. The addition of a thin disk to the base of the capsule is shown to remove this static instability. However, the extremely high entry velocity of 12.6 km/s for the proposed trajectory introduces difficult design issues for incorporating this disk caused by the high aerothermal loads that occur even under relatively rarefied conditions.

Wilmoth, Richard G.↗

Radiative cooling effects on reverse shocks formed by magnetized supersonic plasma flows

We study the structure of reverse shocks formed by the collision of supersonic, magnetized plasma flows driven by an inverse (or exploding) wire array with a planar conducting obstacle. We observe that the structure of these reverse shocks varies dramatically with wire material, despite the similar upstream flow velocities and mass densities. For aluminum wire arrays, the shock is sharp and well-defined, consistent with magneto-hydrodynamic theory. In contrast, we do not observe a well-defined shock using tungsten wires, and instead we see a broad region dominated by density fluctuations on a wide range of spatial scales. We diagnose these two very different interactions using interferometry, Thomson scattering, shadowgraphy, and a newly developed imaging refractometer that is sensitive to small deflections of the probing laser corresponding to small-scale density perturbations. We conclude that the differences in shock structure are most likely due to radiative cooling instabilities, which create small-scale density perturbations elongated along magnetic field lines in the tungsten plasma. These instabilities grow more slowly and are smoothed by thermal conduction in the aluminum plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the Gortler vortex instability mechanism at hypersonic speeds

The linear instability of the hypersonic boundary layer on a curved wall is considered. As a starting point the viscosity of the fluid is taken to be a linear function of temperature and real-gas effects are ignored. It is shown that the flow is susceptible to Gortler vortices and that they are trapped in the logarithmically thin adjustment layer in which the temperature of the basic flow changes rapidly to its free stream value. The vortices decay exponentially in both directions away from this layer and are most unstable when their wavelength is comparable with the depth of the adjustment layer. The non-uniqueness of the neutral stability curve associated with incompressible Gortler 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. It is shown that in the hypersonic limit wall-cooling has a negligible effect on the stability of a fluid with a given value of the Chapman constant.

Hall, Philip↗

Instabilities and subharmonic resonances of subsonic heated round jets, volume 2

When a jet is perturbed by a periodic excitation of suitable frequency, a large-scale coherent structure develops and grows in amplitude as it propagates downstream. The structure eventually rolls up into vortices at some downstream location. The wavy flow associated with the roll-up of a coherent structure is approximated by a parallel mean flow and a small, spatially periodic, axisymmetric wave whose phase velocity and mode shape are given by classical (primary) stability theory. The periodic wave acts as a parametric excitation in the differential equations governing the secondary instability of a subharmonic disturbance. The (resonant) conditions for which the periodic flow can strongly destabilize a subharmonic disturbance are derived. When the resonant conditions are met, the periodic wave plays a catalytic role to enhance the growth rate of the subharmonic. The stability characteristics of the subharmonic disturbance, as a function of jet Mach number, jet heating, mode number and the amplitude of the periodic wave, are studied via a secondary instability analysis using two independent but complementary methods: (1) method of multiple scales, and (2) normal mode analysis. It is found that the growth rates of the subharmonic waves with azimuthal numbers beta = 0 and beta = 1 are enhanced strongly, but comparably, when the amplitude of the periodic wave is increased. Furthermore, compressibility at subsonic Mach numbers has a moderate stabilizing influence on the subharmonic instability modes. Heating suppresses moderately the subharmonic growth rate of an axisymmetric mode, and it reduces more significantly the corresponding growth rate for the first spinning mode. Calculations also indicate that while the presence of a finite-amplitude periodic wave enhances the growth rates of subharmonic instability modes, it minimally distorts the mode shapes of the subharmonic waves.

Ng, Lian Lai↗

Magnetic Measurements of Arc Instabilities in Segmented Arc Heater Column

The Aerodynamic Heating Facility (AHF) at NASA Ames Research Center (ARC) generates high enthalpy, supersonic flow environments for ground testing of NASAheat shield materials. This arc jet heater operates by passing a DC arc (<3200A,<2500V) through a gas column (~10cm diameter, 2m length) and ejecting the heatedgas through a supersonic nozzle into a test chamber. The arc plasma is a partially ionized air mixture at low temperature (< 2eV) but needs to be stable for long test times (up to 30min) at high power (< 20MW). Currently, the operational envelope of the facility is limited by arc instabilities including the kink instability. Consequently, it is of interest to characterize these instabilities to better understand facility limits and inform methods to extend stable operational ranges. Due to the high enthalpy and long timescales (i.e., minutes) diagnostics must be non-invasive. This work will cover the development of two sets of magnetic sensors: differential B-dot sensors and tri-axial magnetic hall sensors and will describe initial measurements of arc motion usingthese sensors.

Kink instability↗

On the role of mode resonances in regulating zonal-flow-moderated plasma microturbulence

Abstract The onset of turbulent heat transport at a higher temperature gradient than the critical gradient of linear instability, known as the Dimits shift, is a recurring feature of nonlinear simulations for magnetically confined fusion plasmas. Resonance in the nonlinear coupling between the modes that dominate energy transfer can lead to suppression of turbulence and transport above the linear critical gradient. As an expression of this resonance, gyrokinetic simulations show a quasi-coherent interaction between streamers and sidebands coupled through the zonal flow within the Dimits regime. This mechanism is further confirmed by use of artificial complex frequencies which break the resonance. By incorporating corresponding saturation physics, the standard quasilinear model for rapid head flux prediction is improved, which can now predict reduced heat flux in the Dimits regime. In particular, the triplet correlation time, the lifetime of the nonlinear interaction, is shown to be well-approximated by combinations of linear eigenvalues, and yields good representations of the heat flux variation both in and above the Dimits regime. Thus, a reduced but predictive model for transport near the critical gradient of zonal-flow saturated turbulence now exists.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Report on Year-4 of Water NSTF Matrix Testing: Facility Maintenance and Accident Testing

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report serves as a summary of maintenance and experimental activities during the program’s fourth year of water-based operation. A planned six-month maintenance period began in August 2021, during which major inspections, repairs, cleaning, and installation of new instrumentation and data acquisition hardware were conducted. Most significantly, two heaters that faulted during Year-3 were repaired, allowing the facility to resume use of the full heated section area and full range of available electric power. The remainder of the year consisted of eight months of test operations, during which the facility logged 211 hours of active heating across one bake-out (following the maintenance period) and seven matrix tests; five classifieds as Accepted per NQA-1, one as Trending, and one as Failed. Testing began by performing two repeat cases to confirm expected facility response and behavior during both single- and two-phase flow conditions, ensuring no changes were introduced during the maintenance period that might have altered the thermal-hydraulic characteristics of the facility. In continuation of the power parametric series initiated in previous years, a high-power test case was then performed examining heat removal performance at a decay heat load equivalent to 2.4 MWt, full-scale, a level exceeding maximum design targets. Additional testing then introduced various blockage scenarios along the network piping, examining the effects of partial and complete blockages of the flow paths on the system behavior and heat removal performance. A study of static boiling tests directed at understanding the geysering two-phase instability was also conducted. The loop was filled only to the bottom of the tank outlet, creating an open loop configuration that prevents any natural circulation flow from occurring, and the heaters were powered on until the facility reached saturation conditions. Following, a series of quasi-steady-state conditions were introduced by adjusting the inventory level in the adiabatic chimney piping at decreasingly lower elevations above the heated region. A strong correlation of geysering characteristics and loop level was observed, with flow and temperature excursions decreasing in intensity, but increasing in frequency, as the fill were reduced to lower elevations along the chimney piping. Once the level fell very low in the chimney, at points near the top of the heated section, the system reached a stable state of continuous boiling without any occurrence of geysering eruptions. A final significant testing accomplishment this year was successful completion of an “accident scenario” test, whose operating conditions were based on a prototypic decay heat curve provided by Framatome and scaled for the NSTF. This test began by establishing steady-state, single-phase “normal operation” conditions, before simulating an accident trip where the availability of active cooling systems was lost. Loop temperatures gradually increased until reaching saturation and subsequent two-phase boiling flow. Over the course of an extended operational period along the defined decay heat curve, steam boil-off caused gradual but continued depletion of liquid inventory until reaching a critically low level causing flow stagnation and cessation of natural circulation heat removal. At this point, after nearly 72 hours of continuous operation, a cold refill was performed to replenish the system inventory and allow the facility to re-establish closed loop natural circulation flow and return to a safe operational state.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

On the stability of self-consistent large amplitude waves in a cold plasma. I - Transverse circularly polarized waves in the absence of a large scale magnetic field

It is demonstrated that a self-consistent circularly polarized wave in an otherwise field-free homogeneous cold plasma is unstable to small amplitude perturbations. For either an electron-positron plasma or an electron-proton plasma the instability rate is at least about the order of the effective plasma frequency when the bulk flow speed is zero. For finite bulk flow speeds of the plasma, it is shown that the electron-positron plasma is unstable, again with a growth rate of the order of the effective plasma frequency; it is also shown that the electron-proton plasma is unstable (at least at small wave numbers, k) with a growth rate proportional to k. The calculated instability rates are conservative, for other modes not investigated here may be more unstable. The results of these calculations bear directly on the understanding of plasma systems thought to be driven by large amplitude waves.

Lee, M. A.↗

Eigenmode analysis of the sheared-flow Z-pinch

Experiments have demonstrated that a Z-pinch can persist for thousands of times longer than the growth time of global magnetohydrodynamic (MHD) instabilities such as the m=0 sausage and m=1 kink modes. These modes have growth times on the order of ta=a/vi, where vi is the ion thermal speed and a is the pinch radius. Axial flows with duz/dr ≲ vi/a have been measured during the stable period, and the commonly accepted theory is that this amount of shear is sufficient to stabilize these modes as predicted by numerical studies using the ideal MHD equations. However, these studies only consider specific equilibrium profiles that typically have a modest magnitude for the logarithmic pressure gradient, qP≡d ln P/d ln r, and may not represent experimental conditions. Linear stability of the sheared-flow Z-pinch is studied here via a direct eigen-decomposition of the matrix operator obtained from the linear ideal MHD equations. Several equilibrium profiles with a large variation of qP are examined. Considering a practical range of k, 1/3 ≲ ka ≲ 10, it is shown that the shear required to stabilize m=0 modes can be expressed as duz/dr≥Cγ0/(ka)α. Here, γ0=γ0(ka) is the profile-specific growth rate in the absence of shear, which scales approximately with |qP|. Both C and α are profile-specific constants, but C is order unity and α≈1. It is further demonstrated that even a large value of shear, duz/dr=3vi/a, is not sufficient to provide linear stabilization of the m=1 kink mode for all profiles considered. This result is in contrast to the currently accepted theory predicting stabilization at much lower shear, duz/dr=0.1vi/a, and suggests that the experimentally observed stability cannot be explained within the linear ideal-MHD model.

Angus, J. R. (ORCID:0000000314740002)↗

General Theory of Aerodynamic Instability and the Mechanism of Flutter

The aerodynamic forces on an oscillating airfoil or airfoil-aileron combination of three independent degrees of freedom have been determined. The problem resolves itself into the solution of certain definite integrals, which have been identified as Bessel functions of the first and second kind and of zero and first order. The theory, being based on potential flow and the Kutta condition, is fundamentally equivalent to the conventional wing-section theory relating to the steady case. The air forces being known, the mechanism of aerodynamic instability has been analyzed in detail. An exact solution, involving potential flow and the adoption of the Kutta condition, has been analyzed in detail. An exact solution, involving potential flow and the adoption of the Kutta condition, has been arrived at. The solution is of a simple form and is expressed by means of an auxiliary parameter K.

Theodorsen, Theodore↗

Hysteresis and Wavenumber Vacillation in Unstable Baroclinic Flows

Hysteresis and wavenumber vacillation are studied numerically in a weakly stratified quasigeostrophic model. In general, the amplitude of the most unstable wave increases, as the flow becomes more unstable. When the wave becomes saturated, the next longer wave will grow at the expanse of the most unstable wave and becomes the dominant wave. However, once the longwave state is established, it may remain in that regime even as the instability is decreased beyond the threshold where it first developed, thus constituting a hysteresis loop. In a highly unstable case, the flow may not show a preference for any single wave. Instead, the dominant wave aperiodically varies among several long waves. This phenomenon is known as wavenumber vacillation. Hysteresis is further examined in terms of eddy heat flux. It is shown that total eddy heat flux increases as the flow becomes more unstable, but displays a sharp drop when transition to a longer wave occurs. However, in a longwave state, the heat flux always decreases with decreasing instability even pass the threshold when wave transition first occurs.

Chou, Shih-Hung↗