Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid 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 613 records · Page 34

Taylor–Couette flow for astrophysical purposes

A concise review is given of astrophysically motivated experimental and theoretical research on Taylor–Couette flow. The flows of interest rotate differentially with the inner cylinder faster than the outer, but are linearly stable against Rayleigh’s inviscid centrifugal instability. At shear Reynolds numbers as large as 10 6 , hydrodynamic flows of this type (quasi-Keplerian) appear to be nonlinearly stable: no turbulence is seen that cannot be attributed to interaction with the axial boundaries, rather than the radial shear itself. Direct numerical simulations agree, although they cannot yet reach such high Reynolds numbers. This result indicates that accretion-disc turbulence is not purely hydrodynamic in origin, at least insofar as it is driven by radial shear. Theory, however, predicts linear magnetohydrodynamic (MHD) instabilities in astrophysical discs: in particular, the standard magnetorotational instability (SMRI). MHD Taylor–Couette experiments aimed at SMRI are challenged by the low magnetic Prandtl numbers of liquid metals. High fluid Reynolds numbers and careful control of the axial boundaries are required. The quest for laboratory SMRI has been rewarded with the discovery of some interesting inductionless cousins of SMRI, and with the recently reported success in demonstrating SMRI itself using conducting axial boundaries. Some outstanding questions and near-future prospects are discussed, especially in connection with astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗

Summary of Multiphysics Modeling for Sublimation Mass Transfer

This report summarizes the effort to develop a multiphysics modeling framework for sublimation mass transfer processes. The aim of the project was to develop a multiphysics code capable of predicting phase change between solid and gas in a closed container, and the movement of material within such a container when exposed to various exterior environmental conditions. Two modeling frameworks were developed towards this aim: one with a high-fidelity computational fluid dynamics (CFD) structure, and the other as a fast reduced-order model. Ultimately, both attempts were unsuccessful due to instabilities in the code and physical processes for which there is no adequate numerical representation. Both modeling attempts are briefly detailed before providing a brief survey of recent updates in the literature, and a recommendation for future work on this subject.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Liquid Engine Jet Atomization Workshop report

The workshop was oriented to disclose information and unsettled problems to understand the fundamental physical mechanism of the droplet formation process. Based on presentation and discussion of results, recommendations were made which should lead to associated future activities. To accomplish this task, existing observations and experiments, contributing to the basic knowledge, providing data for analytical concept verification, and forming a basis for empirical correlations were solicited. Advanced analytical modeling methods or results from specific studies were requested as well as the experience and advice from injector designers. All effort is directed to advance current analytical techniques, simulating the flow behavior downstream of the injection elements in a liquid rocket combustion chamber. Such a tool can be used to optimize injector designs with respect to short length weight savings, wall material protection, or large heat energy transport to a regenerative cooling fluid, while simultaneously achieving the maximum specific impulse in performance. The liquid atomization process also forms a sound basis for combustion instability analysis.

Gross, Klaus W.↗

Generation of circumferential velocity contours associated with pulsed point suction on a rotating disk

Numerous experimental studies were conducted on the steady, three-dimensional boundary layer over a disk rotating at constant angular speed in an otherwise undisturbed fluid. The subject flow geometry is of interest because it provides a relatively simple way to study the cross-flow instability phenomenon which occurs in three-dimensional boundary layers, as on swept wings. This flow instability results in the formation of a stationary spiral vortex flow field over the disk, as shown by Wilkinson and Malik. Using a hot-wire probe, the spatial wave pattern of stationary vortices, which filled the entire circumference of the disk was mapped. The subject flow instability caused transition-to-turbulent flow as the periphery of the disk was approached. The effect on receptivity and transition of discrete disturbance modes, such as three-dimensional toughness elements and acoustic excitation was investigated. The present study (an extension of the work of Wilkinson and Malik) is focused on the effect of pulsed point suction on flow instability and transition, and consequently, on the classical stationary vortical flow pattern.

Selby, Gregory V.↗

Current Status of the Gasdynamic Mirror Fusion Propulsion Experiment

Nuclear fusion appears to be the most promising concept for producing extremely high specific impulse rocket engines. One particular fusion concept which seems to be particularly well suited for fusion propulsion applications is the gasdynamic mirror (GDM). An experimental GDM device has been constructed at the NASA Marshall Space Flight Center to provide an initial assessment of the feasibility of this type of propulsion system. An initial shakedown of the device is currently underway with initial experiments slated to occur in late 2001. This device would operate at much higher plasma densities and with much larger L/D ratios than previous mirror machines. The high L/D ratio minimizes to a large extent certain magnetic curvature effects which lead to plasma instabilities causing a loss of plasma confinement. The high plasma density results in the plasma behaving much more like a conventional fluid with a mean free path shorter than the length of the device. This characteristic helps reduce problems associated with 'loss cone' microinstabilities. The device has been constructed to allow a considerable degree of flexibility in its configuration thus permitting the experiment to grow over time without necessitating a great deal of additional fabrication.

Emrich, William J., Jr.↗

First Results of the Gasdynamic Mirror Experiment

An experimental Gasdynamic Mirror or GDM device has been constructed at the NASA Marshall Space Flight Center to provide an initial assessment of the applicability of this technology for propulsion systems. This paper presents the first experimental results obtained from the machine. It is intended that this device operate at higher plasma densities and with much larger L/D ratios than previous mirror machines. The high L/D ratio minimizes to a large extent certain magnetic curvature effects which lead to plasma instabilities causing a loss of plasma confinement. The high plasma density results in the plasma behaving much more like a conventional fluid with a mean free path shorter than the length of the device. This characteristic helps reduce problems associated with "loss cone" microinstabilities. The device has been constructed to allow a considerable degree of flexibility in its configuration thus permitting the experiment to grow over time without necessitating a great deal of additional fabrication.

Emrich, William J., Jr.↗

First Results of the Gasdynamic Mirror Fusion Propulsion Experiment

An experimental Gasdynamic Mirror or GDM device has been constructed at the NASA Marshall Space Flight Center to provide an initial assessment of the applicability of this technology for propulsion systems. This paper presents the first experimental results obtained from the machine and an analysis of the types of plasma instabilities likely to be encountered. It is intended that this device operate at higher plasma densities and with much larger L/D ratios than previous mirror machines. The high L/D ratio minimizes to a large extent certain magnetic curvature effects which lead to plasma instabilities causing a loss of plasma confinement. The high plasma density results in the plasma behaving much more like a conventional fluid with a mean free path shorter than the length of the device. This characteristic helps reduce problems associated with "loss cone" microinstabilities. The device has been constructed to allow a considerable degree of flexibility in its configuration thus permitting the experiment to grow over time without necessitating a great deal of additional fabrication.

Emrich, William J., Jr.↗

Stopping-power enhancement from discrete particle-wake correlations in high-energy-density plasmas

Three-dimensional (3D) simulations of electron beams propagating in high-energy-density plasmas using the quasistatic Particle-in-Cell (PIC) code QuickPIC demonstrate a significant increase in stopping power when beam electrons mutually interact via their wakes. Each beam electron excites a plasma wave wake of wavelength ~ 2 π c / ω p e , where c is the speed of light and ω p e is the background plasma frequency. Here we show that a discrete collection of electrons undergoes a beam-plasma-like instability caused by mutual particle-wake interactions that causes electrons to bunch in the beam, even for beam densities n b for which fluid theory breaks down. This bunching enhances the beam's stopping power, which we call “correlated stopping,” and the effect increases with the “correlation number” N b ≡ n b ( c / ω p e ) 3 . For example, a beam of monoenergetic 9.7 MeV electrons with N b = 1 / 8 , in a cold background plasma with n e = 10 26 cm – 3 (450 g cm – 3 DT), has a stopping power of 2.28 ± 0.04 times the single-electron value, which increases to 1220 ± 5 for N b = 64 . The beam also experiences transverse filamentation, which eventually limits the stopping enhancement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electromagnetic ion beam instabilities - Comparison of oneand two-dimensional simulations

Numerical simulations of electromagnetic instabilities, driven by a cool tenuous ion beam propagating along an ambient magnetic field, have been conducted in one and two spatial dimensions. The calculations employ particle ions, fluid electrons, and a predictor corrector scheme for solving the electromagnetic field in two dimensions that is described in some detail. While the principal features of the one-dimensional calculations (which reproduce previous work) are retained, the two-dimensional simulations show some reduction of the overall level of the magnetic field fluctuations. Enhancement of the heating of the beam ions at the expense of the core ions also occurs in the case where the beam density is sufficiently large that the right-hand nonresonant instability dominates. Implications of the results for modeling of the ion foreshock and quasi-parallel shocks are discussed.

Winske, D.↗

Convective Instability in Ice I: Application to Callisto and Ganymede

Laboratory experiments measuring ice rheology suggest that it deforms under the influence of several nonNewtonian creep mechanisms, where the viscosity depends on both strain rate and temperature [Goldsby & Kohlstedt, 2001]. Whether or not a fluid with a purely temperature-dependent viscosity convects can be determined by comparing the Rayleigh number of the system to the critical Rayleigh number (Racr), which depends on rheological, thermal, and physical parameters of the fluid layer. However, in a nonNewtonian fluid where viscosity depends on the strain rate (i.e. velocity), convection can only occur if a temperature or velocity perturbation is issued to the system to lower the viscosity and permit fluid motions. Therefore, whether convection occurs in an ice I layer depends on initial conditions, in addition to rheological, thermal, and physical properties of the layer. We show new results for a scaling between the critical Rayleigh number and perturbation amplitude for grain boundary sliding rheology. This scaling can be used to determine the conditions required to initiate convection in the ice I shell of a generic icy satellite. We use this scaling to judge the convective instability of Ganymede and Callisto's ice shells in the absence of tidal dissipation.

Barr, A. C.↗

Numerical analysis of Microgravity Bubble Separation in Open Channels with ISS Bemchmarks

Spacecraft fluid systems often require specific considerations to operate in a manner similar to those on Earth. For instance, an open wedge shaped channel can serve as a pathway for passively separating bubbles from a liquid in a two phase flow, analogous to buoyancy effects on Earth. As the bubbles merge within the channel, they become increasingly confined in the wedge, driving them upwards and outwards toward the free surface due to capillary forces. These same forces ensure coalescence and allow the bubbles to escape through the free surface. However, in low gravity environments, a challenge arises when bubbles are not sufficiently large. Under such conditions, the bubbles tend to follow trajectories near their inscribed elevations, never fully leaving the liquid. Such unseparated bubbles can lead to downstream pump failure, cavitation, flow instabilities, dryout, and other adverse effects aboard spacecraft. Understanding and exploiting the underlying mechanisms behind this behavior is crucial for enhancing the reliability of passive capillary fluids management in space applications, including cryogenic and storable fuels transport, thermal fluids circulation, water recycling for life support, plant watering systems, and more. The joint German Aerospace Center (DLR) and NASA Capillary Channel Flow (CCF) experiment conducted on the International Space Station has provided an extensive publicly available database documenting such phenomena (https://science.nasa.gov/physical sciences informatics psi). Leveraging this resource, we aim to establish benchmarks for our numerical investigation of zero gravity bubbly two phase flow in open wedge channels. Specifically, we seek to quantitatively understand the mechanisms governing ‘inscribed’ bubble motion, considering the apparent interplay of two competing forces: Saffman lift (causing upward/outward motion) and the Magnus effect (causing downward/inward motion). Our cross cutting findings directly inform the development of robust passive phase separation methods within spacecraft plumbing systems.

bubble↗

A study of baroclinic instability in a cylindrical annulus with the temperature gradient imposed on the lower surface

Results are presented from experimental and numerical-modeling studies carried out on a rotating thermally driven fluid system in a cylindrical annulus with the horizontal gradient imposed upon the lower horizontal surface. The rotation rate and the temperature difference were varied to construct a regime diagram in thermal Rossby-number/Taylor-number space. It is shown that the curve separating the axisymmetric flow from the wave flow is 'knee-shaped', similar to the conventional side-heated and side-cooled baroclinic annulus. The experimentally observed transition curve agrees well with numerical calculations and the agreement between the predicted and measured wavenumbers is good both near the transition and within the wave regime. Away from the transition curve and well within the wave regime, there is an evolution from larger to smaller wavenumbers as the flow equilibrates, which is well captured by the numerical model for these cases.

Miller, Timothy L.↗

Fluid-Dynamic Mechanisms Underlying Wind Turbine Wake Control with Strouhal-Timed Actuation

A reduction in wake effects in large wind farms through wake-aware control has considerable potential to improve farm efficiency. This work examines the success of several emerging, empirically derived control methods that modify wind turbine wakes (i.e., the pulse method, helix method, and related methods) based on Strouhal numbers on the O(0.3). Drawing on previous work in the literature for jet and bluff-body flows, the analyses leverage the normal-mode representation of wake instabilities to characterize the large-scale wake meandering observed in actuated wakes. Idealized large-eddy simulations (LES) using an actuator-line representation of the turbine blades indicate that the n=0 and ±1 modes, which correspond to the pulse and helix forcing strategies, respectively, have faster initial growth rates than higher-order modes, suggesting these lower-order modes are more appropriate for wake control. Exciting these lower-order modes with periodic pitching of the blades produces increased modal growth, higher entrainment into the wake, and faster wake recovery. Modal energy gain and the entrainment rate both increase with streamwise distance from the rotor until the intermediate wake. This suggests that the wake meandering dynamics, which share close ties with the relatively well-characterized meandering dynamics in jet and bluff-body flows, are an essential component of the success of wind turbine wake control methods. A spatial linear stability analysis is also performed on the wake flows and yields insights on the modal evolution. In the context of the normal-mode representation of wake instabilities, these findings represent the first literature examining the characteristics of the wake meandering stemming from intentional Strouhal-timed wake actuation, and they help guide the ongoing work to understand the fluid-dynamic origins of the success of the pulse, helix, and related methods.

17 WIND ENERGY↗

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↗

Enhanced shear stabilization of turbulence in NSTX

In studying a particular non-stationary NSTX L-mode plasma, we observed unexpectedly high levels of flux—first with quasilinear (TGLF) modeling and subsequently with nonlinear gyrokinetic simulations. Upon more detailed analysis, a novel confinement regime was discovered in which a modest increase in E x B shear (beyond baseline experimental estimates) rapidly reduced turbulent transport to levels consistent with power balance. This modest increase is plausible given the errors inherent to the estimation of shearing rates, and the added complexity of the non-stationary (time-dependent) power balance. Remarkably, an additional small increase in shear yields the familiar ion-neoclassical transport level with what appears to be the onset of high-k electron transport only. Although analyses using the TGLF-SAT2 model successfully capture numerous parametric dependencies of this plasma, TGLF does not reproduce the rapid E x B stabilization seen in CGYRO. We believe the results presented should help to better characterize the nonlinear physics of spherical tokamak confinement regimes, provide useful ST datasets for reduced model development, and motivate more accurate experimental diagnosis of E x B shearing rates.

Atomic and molecular collisions↗

Rotating spokes, potential hump and modulated ionization in radio frequency magnetron discharges

Here, in this work, the transition from the gradient drift instability (GDI) into an m = 1 rotating spoke in the radio frequency magnetron discharge was studied by means of the two-dimensional axial-azimuthal (z–y) particle-in-cell/Monte Carlo collision method. The kinetic model combined with the linear analysis of the perturbation revealed that the cathode sheath (axial) electric field E z triggers the GDI, deforming the local potential until the instability condition is not fulfilled and the fluctuation growth stops in which moment the instability becomes saturated. The potential deformation consequently leads to the formation of the potential hump, surrounding which the azimuthal electric field E y is present. The saturation level of E y was found to be synchronized with and proportional to the time-changing voltage applied at the cathode, resulting in the RF-modulation of the electron heating in the E y due to $\nabla B$ drift. In the saturated stage of the instability, it was shown that the rotation velocity and direction of the spoke present in the simulations agree well with the experimental observation (Panjan 2019 J. Appl. Phys.125 203303). In the instability linear stage, the instability mode wavelength and the growth rate were also found to be in good agreement with the prediction of the GDI linear fluid theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗