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,783 records · Page 99

Spontaneous magnetization of collisionless plasma

Here we study within a fully kinetic framework the generation of “seed” magnetic fields through the Weibel instability, driven in an initially unmagnetized plasma by a large-scale shear force. We develop an analytical model that describes the development of thermal pressure anisotropy via phase mixing, the ensuing exponential growth of magnetic fields in the linear Weibel stage, and the saturation of the Weibel instability when the seed magnetic fields become strong enough to instigate gyromotion of particles and thereby inhibit their free-streaming. The predicted scaling dependencies of the saturated fields on key parameters (e.g., ratio of system scale to electron skin depth and forcing amplitude) are confirmed by two-dimensional and three-dimensional particle-in-cell simulations of an electron–positron plasma. This work demonstrates the spontaneous magnetization of a collisionless plasma through large-scale motions as simple as a shear flow and therefore has important implications for magnetogenesis in dilute astrophysical systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Boundary-Layer Transition Prediction Through Loose Coupling of OVERFLOW and LASTRAC

Transition prediction based on linear stability theory is expected to more accurately reflect the causality of transition onset than phenomenological transition models based on RANS-like transport equations. To help achieve the CFD vision 2030 aim of building a CFD tool chain with automated prediction of boundary layer transition, a technique to loosely tie the NASA OVERFLOW CFD solver with the LASTRAC stability analysis tool is described. The coupled solver is then used to compute transition over selected over a flat plate in a freestream with sufficiently low levels of turbulence, NLF(1)-0416 airfoil, the 6:1 prolate spheroid at an angle of attack, and a NASA juncture flow model with symmetric wing configuration. The findings show that the loosely coupled approach can reliably predict the transition location accurately in scenarios that are dominated by a single transition mechanism involving Tollmien-Schlichting instabilities, crossflow instabilities, or separation bubble-induced transition, or include a mixture of selected mechanisms. The toolset presents here appears to be robust to the prescription of the initial transition location, and it can lead to a converged solution in four or five rounds of the mean flow calculation and stability analysis, with minimal input from the user.

boundary layer transition↗

Z-pinch fusion

The Z-pinch configuration offers the promise of a compact fusion device owing to its simple geometry, unity beta, and absence of external magnetic field coils. Increasing the axial current compresses the plasma, resulting in a rapid rise of the fusion reaction rate. Historically, the Z pinch has been plagued by fast growing instabilities that limit plasma lifetimes. Recent progress has resulted from investigating approaches that provide stability. One approach exploits sheared-flow stabilization to produce an equilibrium Z pinch, which sustains the compressed plasma state for durations much longer than other plasma timescales. Recent experimental and simulation results encourage pursuit of the stabilized Z pinch to explore its fusion performance limits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Preface for the Sherwood Fusion Theory 2021–2022 special collection

The Sherwood Fusion Theory conference has a 60-year history and has been held nearly annually for at least five decades. With theoretical and computational research directed toward fusion energy as an organizing principle, the typical Sherwood program encompasses a wide range of plasma physics topics pertaining to a variety of confinement concepts along with advancements in algorithms and numerical methods for fusion plasma simulations. Additionally, the typical annual meeting comprises roughly a dozen invited talks plus a small number of plenary speakers and several poster sessions. The program committee presents several awards each year to graduate (or undergraduate) students with outstanding poster presentations. Beginning in 2021, both invited speakers and student poster award winners were asked to submit manuscripts for this special collection, which combines contributions from the 2021 and 2022 Sherwood meetings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Load Margin Constrained Moving Target Defense against False Data Injection Attacks

Cyber physical security of power systems with high penetration of renewable generation has attracted attention from researchers. One critical issue is that cyber-physical attacks, disguised as uncertain renewable generation, can target conventional power system state estimation (SE). Moving target defense (MTD) is a promising defense strategy to detect stealthy false data injection (FDI) attacks against SE. However, all existing studies myopically perturb the reactance of transmission lines equipped with distributed flexible AC transmission system (D-FACTS) devices without adequately considering the system voltage stability. Exacerbated by the renewable generation uncertainty, existing MTD may cause voltage instability when the power grid is under stress. To address this issue, we propose a novel MTD framework that explicitly considers system voltage stability by using continuation power flow. We utilize the sensitivity matrix of power injection to line impedance, on which an optimization problem for maximizing load margin is formulated. This framework is validated on the IEEE 14-bus system and the IEEE 118-bus system, in which net load redistribution attacks are launched by sophisticated attackers. Steady-state simulations and dynamic simulations on PSS/E show the effectiveness of the proposed framework in circumventing the voltage instability while maintaining the detection effectiveness of MTD. The impact of the proposed method on attack detection effectiveness is also revealed.

Zhang, Hang↗

Structure and Energetics of Medium-scale Atmospheric Waves in the Southern Hemisphere Summer

The NMC analyses were used to study medium-scale waves during three Southern Hemisphere (SH) summers. The SH summer medium-scale waves are observed to exhibit remarkably regular eastward phase progression, and wave maxima can often be traced continuously around the globe. Frequent downstream development of existing wave patterns is observed. The medium-scale waves at times appear to be longitudinally localized features; at other times they resemble truly global-scale models. The time-mean wave structure is found to be consistent with basic-state propagation characteristics and the conservation of wave activity. There is a tendency for background wind structure to restrict wave propagation. From energetics studies it is found that wave-zonal-mean exchange is a valid concept for describing the SH summer circulation, and that the flow vacillated between periods of highly perturbed and zonally symmetric states. The medium-scale waves result from nonlinear baroclinic instabilities, and exhibit a well-defined life cycle of baroclinic growth, maturity, and barotropic decay.

Randel, W. J.↗

Anti-swirling demonstration

Fluid in bearings and seals, involved in motion due to shaft rotation and generating fluid dynamic forces, represents the well recognized source of instability known as whirl and whip. It is characterized by shaft forward circular precession at subysnchronous frequency. Control of the flow pattern in the bearing/seal clearance may enhance the stability. The rotor system demonstrates the dynamic effect of perturbation of the flow pattern in a seal gap by blowing compressed air in the rotor tangential direction. This effect results in shifting the rotor threshold of stability to a higher value of rotative speed when the air jet direction opposes rotation, and to a lower value of rotative speed when the air jet direction coincides with rotation. This effect implies that an unstable rotor can be stabilized by blowing air in the direction opposite to rotation. This technique is known as anti-swirl.

Muszynska, A.↗

Straight cylindrical seal for high-performance turbomachines

A straight cylindrical seal configuration representing the seal for a high-performance turbopump (e.g., the space shuttle main engine fuel pump) was tested under static (nonrotating) conditions. The test data included critical mass flux and pressure profiles over a wide range of inlet temperatures and pressures for fluid nitrogen and fluid hydrogen with the seal in concentric and fully eccentric positions. The critical mass fluxes (or leakage rates) for the concentric and fully eccentric configurations were nearly the same when based on stagnation conditions upstream of the seal. The fully eccentric configuration pressure profiles of the gas and liquid were different. Further, the pressure differences between the maximum and the minimum clearance positions were highly dependent on the geometric conditions, the temperature, and the absolute pressure at both the inlet and the exit. The pressure differences were greatest in the inlet region. The results, although complex, tend to follow the corresponding-states principles for critical flows. Gaseous injection near the seal exit plane significantly altered the pressure profiles and could be used to control turbomachine instabilities.

Hendricks, Robert C.↗

Low-dimensional chaos in turbulence

Direct numerical simulations are being performed on two different fluid flows in an attempt to discover the mechanism underlying the transition to turbulence in each. The first system is Taylor-Couette flow; the second, two-dimensional flow over an airfoil. Both flows exhibit a gradual transition to high-dimensional turbulence through low-dimensional chaos. The hope is that the instabilities leading to chaos will be easier to relate to physical processes in this case, and that the understanding of these mechanisms can then be applied to a wider array of turbulent systems.

Vastano, John A.↗

Observations of the frequencies in a sphere wake and of drag increase by acoustic excitation

Vortex shedding and instability wave frequencies have been measured in the wakes of spheres in the Reynolds number range 500 less than Re less than 60,000. The effect of acoustic excitation was examined and an interaction between the two frequency modes was found at the lower Reynolds numbers; through this interaction, external forcing at the instability frequency could change the vortex shedding frequency. The development of the mean wake was manipulated by forcing near to the dominant shear layer instability frequency. With this forcing, the separated shear layer moved closer to the surface of the sphere and the reversed flow region of the wake was shortened. Concomitantly, the base pressure decreased and drag increased.

Kim, H. J.↗

Laminar Dust Flames: A Program of Microgravity and Ground Based Studies at McGill

Fundamental knowledge of heterogeneous combustion mechanisms is required to improve utilization of solid fuels (e.g. coal), safe handling of combustible dusts in industry, and solid propulsion systems. The objective of the McGill University research program on dust combustion is to obtain a reliable set of data on basic combustion parameters for dust suspensions (i.e. laminar burning velocity, flame structure, quenching distance, flammability limits, etc.) over a range of particle sizes, dust concentrations, and types of fuel. This set of data then permits theoretical models to be validated and, when necessary, new models to be developed to describe the detailed reaction mechanisms and transport processes. Microgravity is essential to the generation of a uniform dust suspension of arbitrary particle size and concentration. When particles with a characteristic size on the order of tens of microns are suspended, they rapidly settle in a gravitational field. To maintain a particulate in suspension for time duration adequate to carry out combustion experiments invariably requires continuous convective flow in excess of the gravitational settling velocity (which is comparable with and can even exceed the dust laminar burning velocity). This makes the experiments turbulent in nature and thus renders it impossible to study laminar dust flames. Even for small particle sizes on the order of microns, a stable laminar dust flow can be maintained only for relatively low dust concentrations at normal gravity conditions. High dust loading leads to gravitational instability of the dust cloud and to the formation of recirculation cells in the dust suspension in a confined volume, or to the rapid sedimentation of the dense dust cloud, as a whole, in an unconfined volume. Many important solid fuels such as carbon and boron also have low laminar flame speeds (of the order of several centimeters per second). Convection that occurs in combustion products due to buoyancy disrupts the low speed dust flames and makes observation of such flames at normal gravity difficult.

Goroshin, Sam↗

Fluid Physics in a Fluctuating Acceleration Environment

Our program of research aims at developing a stochastic description of the residual acceleration field onboard spacecraft (g-jitter) to describe in quantitative detail its effect on fluid motion. Our main premise is that such a statistical description is necessary in those cases in which the characteristic time scales of the process under investigation are long compared with the correlation time of g-jitter. Although a clear separation between time scales makes this approach feasible, there remain several difficulties of practical nature: (i), g-jitter time series are not statistically stationary but rather show definite dependences on factors such as active or rest crew periods; (ii), it is very difficult to extract reliably the low frequency range of the power spectrum of the acceleration field. This range controls the magnitude of diffusive processes; and (iii), models used to date are Gaussian, but there is evidence that large amplitude disturbances occur much more frequently than a Gaussian distribution would predict. The lack of stationarity does not constitute a severe limitation in practice, since the intensity of the stochastic components changes very slowly during space missions (perhaps over times of the order of hours). A separate analysis of large amplitude disturbances has not been undertaken yet, but it does not seem difficult a priori to devise models that may describe this range better than a Gaussian distribution. The effect of low frequency components, on the other hand, is more difficult to ascertain, partly due to the difficulty associated with measuring them, and partly because they may be indistinguishable from slowly changing averages. This latter effect is further complicated by the lack of statistical stationarity of the time series. Recent work has focused on the effect of stochastic modulation on the onset of oscillatory instabilities as an example of resonant interaction between the driving acceleration and normal modes of the system, and on cavity flow as an example of how an oscillatory response under periodic driving becomes diffusive if the forcing is random instead. This paper describes three different topics that illustrate behavior that is peculiar to a stochastic acceleration field. In the first case, we show that g-jitter can induce effective attractive or repulsive forces between a pair of spherical particles that are suspended in an incompressible fluid of different density provided that the momentum diffusion length is larger than the interparticle separation (as in the case in most colloidal suspensions). Second, a stochastic modulation of the control parameter in the vicinity of a pitchfork or supercritical bifurcation is known not to affect the location of the threshold. We show, however, that resonance between the modulation and linearly stable modes close to onset can lead to a shift in threshold. Finally, we discuss the classical problem of vorticity diffusion away from a plane boundary that is being vibrated along its own plane. Periodic motion with zero average vorticity production results in an exponential decay of the vorticity away from the boundary. Random vibration, on the other hand, results in power law decay away from the boundary even if vorticity production averages to zero.

Drolet, Francois↗

Highly Structured Slow Solar Wind Emerging From an Equatorial Coronal Hole

At solar minimum, the solar wind is observed at high solar latitudes as a predominantly fast (> 500 km/s), highly Alfvenic, rarefied stream of plasma originating deep within coronal holes, while near the ecliptic plane it is interspersed with a more variable slow (< 500 kms) wind. The precise origins of the slow wind streams are less certain, with theories and observations supporting sources from the tips of helmet streamers, interchange reconnection near coronal hole boundaries, and origins within coronal holes with highly diverging magnetic fields. The heating mechanism required to drive the solar wind is also an open question and candidate mechanisms include Alfven wave turbulence, heating by reconnection in nanoflares, ion cyclotron wave heating and acceleration by thermal gradients1. At 1 au, the wind is mixed and evolved and much of the diagnostic structure of these sources and processes has been lost. Here we present new measurements from Parker Solar Probe at 36 to 54 solar radii that show clear evidence of slow, Alfvenic solar wind emerging from a small equatorial coronal hole. The measured magnetic field exhibits patches of large, intermittent reversals associated with jets of plasma and enhanced Poynting flux and interspersed in a smoother and less turbulent flow with near-radial magnetic field. Furthermore, plasma wave measurements suggest electron and ion velocity-space micro-instabilities that have been identified with plasma heating and thermalization processes. Our measurements suggest an impulsive mechanism associated with solar wind energization and a heating role for micro-instabilities and provide strong evidence for low latitude coronal holes as a significant contribution to the source of the slow solar wind.

Solar physics↗

From Disturbances to Instabilities, to Breakdown to Turbulence: The Physics of Transition in Boundary Layers

In order to understand the end-stages of boundary layer transition in low as well as high disturbance environments it is desirable to establish a unified view of the sequences of physico-mathematical phenomena that lead from laminar flow to self-sustained "bursting" in wall turbulence. The dominant driving disturbances: oncoming free turbulence, unsteady pressure fields, inhomogeneous density fields, inhomogeneities in wall geometry, all force disturbed motions within the boundary layer via multiple competitive receptivity mechanisms. For small disturbances, a sequence of instabilities then leads to sporadic local bursting very near the wall which can sustain turbulence. The local seeds of turbulence then somehow propagate to engulf quite rapidly the surrounding disturbed but still laminar regions.

Morkovia, Mark V.↗

Unsteady flow in cavitating turbopumps

Unsteady flow in a cavitating axial inducer pump is analyzed with the help of a simple two-dimensional cascade model. This problem was motivated by a desire to study the effect of unsteady cavitation on the so-called POGO instability in the operation of liquid rocket engines. Here, an important feature is a closed loop coupling between several different modes of oscillation, one of which is due to the basic unsteady characteristics of the cavitation itself. The approaching and leaving flow velocities up- and downstream of the inducer oscillate, and the cavity-blade system participates dynamically with the basic pulsating flow. In the present work, attention is focused on finding a transfer matrix that relates the set of upstream variables to those downstream. This quantity, which is essentially equivalent to cavitation compliance in the quasi-static analyses, is found to be complex and frequency dependent.

Kim, J. H.↗

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work is limited to very simple models. In this project, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. The proposed work is based on previously funded NASA work that allowed us to develop a two-dimensional numerical method for boiling flows and to demonstrate the ability of the method to simulate film boiling. While numerical simulations of multi-fluid flows have been advanced in a major way during the last five years, or so, similar capability for flows with phase change are still in their infancy. Although the feasibility of the proposed approach has been demonstrated, it has yet to be extended and applied to fully three-dimensional simulations. Here, a fully three-dimensional, parallel, grid adaptive code will be developed. The numerical method will be used to study nucleate boiling in microgravity, with particular emphasis on two aspects of the problem: 1) Examination of the growth of bubbles at a wall nucleation site and the instabilities of rapidly growing bubbles. Particular emphasis will be put on accurately capturing the thin wall layer left behind as a bubble expands along a wall, on computing instabilities on bubble surfaces as bubbles grow, and on quantifying the effects of both these phenomena on heat transfer; and 2) Examination of the effect of shear flow on bubble growth and heat transfer.

Tryggvason, Gretar↗

Electrostatic ion-cyclotron instability caused by a nonuniform electric field perpendicular to the external magnetic field

A new mechanism that can destablize kinetic ion-cyclotron waves in the presence of a nonuniform electric field perpendicular to the uniform ambient magnetic field is given. In the absence of the electric field, the mode energy is positive, while in the presence of a uniform electric field the mode energy could be negative. However, when the electric field is nonuniform, it is possible for a finite region to be of negative wave energy surrounded by regions of positive wave energy. A nonlocal wave packet couples the two regions so that a flow of energy from the region of negative wave energy to the region of positive wave energy will cause the mode to grow. This gives rise to the instability.

Ganguli, G.↗

Sensitivity of mixing layers to three-dimensional forcing

It is well known that turbulent mixing layers are dominated by large scale, fairly coherent structures, and that these structures are related to the stability characteristics of the flow. These facts have led researchers to attempt controlling such flows by selectively forcing certain unstable modes, which can in addition have the effect of suppressing other modes. Much of the work on controlling the mixing layer has relied on forcing 2-D instabilities. The results of forcing 3-D instabilities are addressed. The objectives of the work are twofold: to understand how a mixing layer responds to 3-D perturbations, and to test the validity of an amplitude expansion in predicting the mixing layer development. The amplitude expansion could be very useful in understanding and predicting the 3-D response of the flow to a variety of initial conditions.

Riley, J. J.↗