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,099 records · Page 61

The structure and dynamics of reacting plane mixing layers

Experimental studies of the combusting mixing layers formed over a rearward facing step and between two parallel streams are reported. The steady behavior of the mean properties of the reacting flows is similar to that of nonreacting flows except for the global effects of thermodynamic property changes. Vortices form by the two-dimensional K-H instability and grow by subharmonic combination until the layer begins to interact with the walls. Unsteady behavior occurs when there are resonator interactions with the K-H instability or the instability associated with the recirculation vortex in the rearward facing step flow.

Keller, J. O.↗

Pore-Scale Modeling of Mineral Growth and Nucleation in Reactive Flow

A fundamental understanding of mineral precipitation kinetics relies largely on microscopic observations of the dynamics of mineral surfaces exposed to supersaturated solutions. Deconvolution of tightly bound transport, surface reaction, and crystal nucleation phenomena still remains one of the main challenges. Particularly, the influence of these processes on texture and morphology of mineral precipitate remains unclear. This study presents a coupling of pore-scale reactive transport modeling with the Arbitrary Lagrangian-Eulerian approach for tracking evolution of explicit solid interface during mineral precipitation. It incorporates a heterogeneous nucleation mechanism according to Classical Nucleation Theory which can be turned “on” or “off.” This approach allows us to demonstrate the role of nucleation on precipitate texture with a focus at micrometer scale. In this work precipitate formation is modeled on a 10 micrometer radius particle in reactive flow. The evolution of explicit interface accounts for the surface curvature which is crucial at this scale in the regime of emerging instabilities. The results illustrate how the surface reaction and reactive fluid flow affect the shape of precipitate on a solid particle. It is shown that nucleation promotes the formation of irregularly shaped precipitate and diminishes the effect of the flow on the asymmetry of precipitation around the particle. The observed differences in precipitate structure are expected to be an important benchmark for reaction-driven precipitation in natural environments.

58 GEOSCIENCES↗

Fuel Injection Dynamics and Composition Effects on RDE Performance

Rotating detonation engines (RDEs) provide a promising route to substantially increasing cycle efficiency in stationary gas turbines. Much of this increase relies on the ability to achieve consistent pressure gain within the combustor. In particular, the design of injectors that feed fuel and air into the detonation channel plays a crucial role. Such injectors have to ensure proper mixing of fuel and oxidizer, while minimizing backflow of detonation products into the feed plenums, and reduce susceptibility to the complex wave structures that exist within the combustor. From a practical perspective, such RDEs also need to operate with variable fuel composition. When fuel mixtures with components that possess vastly different oxidation pathways and time-scales are used, there could be additional losses through deflagrative burning instead of detonation-driven heat release. Such sensitivity to the complete flow path is akin to the physics of thermoacoustic instabilities in conventional gas turbines. In this sense, RDEs pose a unique research challenge: the performance of the device relies on the small-scale heat release process, which is highly dependent on the flow interactions within the full-scale system. As a result, canonical flow configurations, instrumented with detailed diagnostics or modeled using high-fidelity tools, but only focus on the small-scale processes will not contain the key system-level interactions. At the same time, macroscopic measurements and models that only capture system-level performance will not provide insight into the key sources of pressure losses. These couplings and sensitivities provide a formidable challenge to both experimental and simulation studies of the effects Thus, a joint experimental/computational program designed specifically to address these challenges was undertaken in this program. The focus of this program was on two key topics: a) the interaction between injector flow and the overall wave dynamics within the combustor, and b) the deflagration/detonation structure in multi-component fuels that are of practical interest. Both topics involve interaction of small-scale heat release processes with the geometry-dependent wave structure. Studies focused on the study of full-scale RDE systems, based on a 6-inch conventional annular geometry. Experimentally RDEs were studied using a combination of diagnostics. A combination of optical diagnostics and aero-thermo-acoustic analysis based on a combination of spectral and mode decomposition analysis was used to identify the dynamics of the detonation wave and other secondary waves that exist in the system. These studies have helped the identification and investigation of injector and detonation dynamics arising from coupling, and how they affect RDE mixing, detonation structure, operability and performance. Performance of RDEs was investigated through thrust stand measurements, which was used to evaluate the effective pressure gain generated by the system through the concept of equivalent available pressure. Optical diagnostics were developed and implemented to investigate the distribution of heat release, across the detonation wave. Novel optical diagnostics of NIR imaging was also developed and applied to investigate the high temperature / high pressure distribution across the detonation wave. In order to complement the experiments, the computational tools were geared to simulate the full experimental setup. GPU-based acceleration of the models and computations were developed to enable rapid simulation of the full system. In addition, the use of adaptive mesh refinement, and unstructured grid formulation, enabled the investigation of realistic geometries studied in the laboratory. The simulations produced a wealth of detail on the structure of the detonation wave under different operating conditions. Emphasis was placed on quantifying mixture pre-burning and the impact on wave propagation and structure.

03 NATURAL GAS↗

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↗

Separation and Transition on a Cone-Cylinder-Flare: Computational Investigations

Base flow computation and stability analysis were conducted for hypersonic flow over a cone-cylinder-flare (CCF) geometry for conditions that correspond to the experimental runs carried out in three wind tunnels. Owing to the presence of an attached boundary layer, a separation bubble induced by a shock-boundary layer interaction, and a reattachment region, the chosen flow configuration is physically rich. The complexity of this flowfield encompasses a combination of convective instabilities developing on the cone, global instabilities in the separation bubble, and shear-layer modes and streaks in the reattachment region. Thus, the selected CCF configuration provides the opportunity for a comprehensive comparison of the currently available methodologies for analyzing boundary layer instabilities. Various tools are used for the analysis, including global stability codes as well as convective instability analyses based on a local theory, a weakly non-parallel analysis, and tools that are applicable to strongly non-parallel flows. The paper presents a comparison of the convective instability characteristics based on different methodologies, such as linear stability theory (LST), the harmonic form of linearized Navier-Stokes equations (HLNSE), and resolvent analysis. The CCF configuration provided an effective framework for conducting a detailed cross-validation of this type, which had not yet been addressed in existing literature. This document is accompanied by a companion paper that is focused on the experimental aspects of the CCF configuration. Both papers are being presented in a dedicated session that highlights the research activities of the NATO STO Research Task Group AVT-346.

Boundary Layer Transition↗

The onset and development of Kelvin-Helmholtz instability at the Venus ionopause

The paper investigates the Kelvin-Helmholtz instability at the Venus ionopause resulting from the flow of the (shocked) solar wind tangential to the ionopause for the case where the interplanetary field is oriented normal to the direction of flow. It is found that gravity stabilizes the long wavelength perturbations, and the finite thickness of the boundary layer stabilizes short wavelength modes. The magnetic 'gyroviscosity' due to finite Larmor radius effects either destabilizes the boundary or stabilizes it according to whether the solar wind electric field points away from or toward the ionosphere. For solar wind and ionosphere plasma parameters consistent with Pioneer Venus observations, it is found that the instabilities with the greatest growth rates (shortest growth times) have wavelengths of 50-150 km and growth times of 0.5 to several seconds. In addition, it is found how distortion of the ionopause by Kelvin-Helmholtz instability might lead to the formation of magnetic 'flux ropes' inside the ionosphere as well as ionospheric 'bubbles' embedded in the solar wind.

Wolff, R. S.↗

Development of Physics-Based Transition Models for Unstructured-Mesh CFD Codes Using Deep Learning Models

Predicting transition locations over a vehicle surface is of fundamental importance for many engineering applications. With the transition information, the Reynolds-averaged Navier-Stokes (RANS) computations can turn on the turbulence model at the right locations so that drag, lift and other aerodynamic quantities can be accurately predicted. In contrast to the popularity of RANS-based transition modeling in which transition onset is governed by the turbulence equations, physics-based transition models that account for instability waves within the boundary layer, thus more compliant to flow physics, only gained more attention in recent years. This paper describes the development of a new physics-based transition model based on either the linear stability theory (LST) or parabolized stability equations (PSE). The model is designed to communicate with a structured or unstructured-mesh RANS solver back and forth in order to more accurately compute transition fronts over a three-dimensional body. In the developed model, the Python suite of interface codes in conjunction with the LASTRAC software can be executed autonomously to produce transition onset locations for a given laminar or RANS-computed transitional state. In addition, as a proof of concept, the tool set consists of a deep learning neural network model that has been designed and trained to predict instability wave evolutions inside the boundary layer for various instability wave mechanisms across a selected speed range. A machine-learned intelligent profile interpolation model has also been devised to enable reliable instability-wave spectra predictions with just a few points in the mean flow profiles.

Transition Modeling↗

The effects of cold plasma on the Kelvin-Helmholtz instability

The results are presented of a detailed study that examines the effects of cold plasma on the Kelvin-Helmholtz instability. Both electrostatic and electromagnetic modes are considered. The interest in the Kelvin-Helmholtz instability stems from observations of low frequency waves in the polar cusp, magnetopause, in fast moving auroral forms, the plasma pause, and other regions of magnetospheric flow boundaries. It is pointed out that the Kelvin-Helmholtz instability is driven by a velocity shear and can excite a broad spectrum of wave frequencies below the ion gyrofrequency. A model is obtained of a flow boundary that includes a density gradient perpendicular to the velocity shear. The analysis shows a compressional component in Kelvin-Helmholtz waves.

Melander, B. G.↗

Dissipation in the magnetic turbulence of reversed field pinch plasmas

Reversed field pinch plasmas are subject to tearing instability that creates a broad spectrum of magnetic fluctuations. The dominant fluctuations have poloidal and toroidal mode numbers (m,n) = (1,6 - 10) and can grow to 2 - 3% of the mean magnetic field. Through non-linear coupling this growth culminates in a strong reconnection event and broadening of the magnetic spectrum extending to the ion gyroradius scale. Multiple developments occur during the reconnection stage: ions and electrons are energized, magnetic fluctuation amplitudes increase, plasma flow is halted, and the toroidal magnetic flux increases in a sawtooth-like fashion as the RFP dynamo becomes stronger. Magnetic fluctuations are measured in the plasma edge at multiple radial locations from r/a = 0.75 to 0.96 to assess and characterize the magnetic turbulence. The measured spectrum perpendicular to the mean field, S(k ⊥ ), can be fit to a model spectrum consisting of power law and exponential component with one free parameter that characterizes dissipation. The measured dissipation is much larger than estimated from classical viscous or resistive dissipation, but it is consistent with a flow damping measurement of anomalous viscosity. The measurements show an evolution of the spectrum during which fluctuation power builds up in the smallest wavenumbers and cascades to the larger wavenumber due to the non-linear coupling between the linear (m,n) = (1, > 6) and the non-linear (m,n) = (0, 1) tearing mo

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The collisionless macroscopic Kelvin-Helmholtz instabiity. I Transverse electrostatic mode

The electrostatic limit of the Kelvin-Helmholtz instability in a collisionless plasma is investigated for a magnetic field perpendicular to the velocity flow. A two-fluid eigenvalue analysis of the transverse electrostatic Kelvin-Helmholtz instability is performed using the Chew-Goldberger-Low hydromagnetic equations with first-order finite Larmer radius corrections, establishing that the transverse electrostatic mode reduces to the incompressible transverse MHD mode in the limit of zero ion gyroradius. The linear growth rates are determined and a strong stabilization of the instability for (k-perpendicular) (rho sub i) of 0.2 or greater is observed. The nonlinear stage of the instability displays large vortices whose size is many times the initial velocity shear length. At saturation the simulation is dominated by the longest-wavelength mode permitted in the system. The simulation results can be understood in terms of a simple model which assumes that saturation occurs when the angular frequency of the vortex motion reaches the linear growth rate of the instability.

Prichett, P. L.↗

Electromagnetic instabilities attributed to a cross-field ion drift

Instabilities due to a cross-field ion flow are reexamined by including the electromagnetic response of the ions, which has been ignored in existing discussions. It is found that this effect can lead to significant enhancement of the growth rate. Among the new results, a purely growing, electromagnetic unstable mode with a wave vector k parallel to the ambient magnetic field is found. The plasma configuration under consideration is similar to that used in the discussion of the well-known modified-two-stream instability. This instability has a growth rate faster than the ion cyclotron frequency, and is not susceptible to high-plasma-beta stabilization.

Chang, C. L.↗

Experiments on annular liquid jet instability and on the formation of liquid shells

An annular jet flow of liquid surrounding a flow of gas at its core is extremely unstable. Experiments are described in which such a flow is generated by an annular nozzle operated at fairly specific conditions. It is shown that periodic, axisymmetric oscillations arise spontaneously within the cylindrical sheet emerging from the nozzle and grow with such rapidity along the axial dimension that a sealing-off and encapsulation of the core gas occurs within a few jet diameters. This is closely followed by a pinchoff of the liquid between adjacent bubbles. The liquid shells set free thereby assume spherically symmetric form under capillary forces, and each contains a precisely uniform measure of gas and of liquid on account of the extremely high frequency-stability of the process. Description is given of the fluid dynamic processes by which the shells are formed, and mention is made of exploiting the instability for the production of rigid shells for technological applications.

Kendall, J. M.↗

Explicit structure-preserving geometric particle-in-cell algorithm in curvilinear orthogonal coordinate systems and its applications to whole-device 6D kinetic simulations of tokamak physics

Explicit structure-preserving geometric particle-in-cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km s -1 , agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2% level by the E × B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.

43 PARTICLE ACCELERATORS↗

The instability of a liquid jet in a compressible airstream

The instability of a liquid jet immersed in a coflowing compressible airstream is investigated, for both two-dimensional and circular jets, using a simplified flow model in which the effects associated with the cross-flow component of velocity are ignored. Results of the analysis agreed qualitatively with experimental results in that they suggest a more violent instability when the free-stream Mach number is close to unity. It was also found that a change from convective to absolute instability occurs in this transonic flow regime.

Li, H.-S.↗

Flow boiling in geometrically modified microchannels

Microchannels are a promising solution for high-heat-flux thermal management scenarios, including high-performance microelectronics cooling and power electronics cooling. However, thermohydraulic instabilities result from the rapid vapor bubble formation. The prior literature has examined several methods, including constricted inlet microchannels, expanding microchannels, and auxiliary jetting microchannels, to mitigate the effect of these instabilities. Computational fluid dynamics and heat transfer (CFD/HT) modeling of the flow boiling phenomena in these microchannel configurations has seen limited examination, and one-to-one numerical comparisons of the different mitigation strategies have not been performed. In the present investigation, CFD/HT analyses using a three-dimensional (3D) volume of fluid model coupled with a phase-change model for the interfacial heat and mass transfer were performed for multiple microchannel configurations (constricted inlet, expanding, and auxiliary jetting microchannels). A benchmark case of a rectangular microchannel was examined to quantify baseline thermohydraulic performance. Results demonstrated slight to significant thermal performance improvements for all cases, and significant pressure benefits for the expanding and jetting cases, consistent with experimental results in the literature. Bubble dynamics and visualization for the baseline and alternative configurations are provided to give insight into their underlying physics, and the differences in performance were investigated and compared with available literature.

Broughton, Justin (ORCID:000000026848457X)↗

Stability of compressible Taylor-Couette flow

The objectives of this paper are to: (1) develop both analytical and numerical tools that can be used to predict the onset of instability and subsequently to simulate the transition process by which the originally laminar flow evolves into a turbulent flow; and (2) conduct the preliminary investigations with the purpose of understanding the mechanisms of the vortical structures of the compressible flow between tow concentric cylinders.

Kao, K.↗

Shear banding, intermittency, jamming, and dynamic phases for skyrmions in inhomogeneous pinning arrays

In this work, we examine driven skyrmion dynamics in systems with inhomogeneous pinning where a strip of strong pinning coexists with a region containing no pinning. For driving parallel to the strip, we find that the initial skyrmion motion is confined to the unpinned region and the skyrmion Hall angle is zero. At larger drives, a transition occurs to a phase in which motion also appears within the pinned region, creating a shear band in the skyrmion velocity, while the skyrmion Hall angle is still zero. As the drive increases further, the flow becomes disordered and the skyrmion Hall angle increases with drive until saturating at the highest drives when the system transitions into a moving crystal phase. The different dynamic phases are associated with velocity and density gradients across the pinning boundaries. We map out the dynamic phases as a function of pinning strength, skyrmion density, and Magnus force strength, and correlate the phase boundaries with features in the velocity-force curves and changes in the local and global ordering of the skyrmion structure. For large Magnus forces, the shear banding instability is replaced by large-scale intermittent flow in the pinned region accompanied by simultaneous motion perpendicular to the direction of the drive, which appears as oscillations in the transport curves. We also examine the case of a drive applied perpendicular to the strip, where we find a jamming effect in which the skyrmion flow is blocked by skyrmion-skyrmion interactions until the drive is large enough to induce plastic flow.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effect of Gravity on the Near Field Flow Structure of Helium Jet in Air

Experiments have shown that a low-density jet injected into a high-density surrounding medium undergoes periodic oscillations in the near field. Although the flow oscillations in these jets at Richardson numbers about unity are attributed to the buoyancy, the direct physical evidence has not been acquired in the experiments. If the instability were indeed caused by buoyancy, the near-field flow structure would undergo drastic changes upon removal of gravity in the microgravity environment. The present study was conducted to investigate this effect by simulating microgravity environment in the 2.2-second drop tower at the NASA Glenn Research Center. The non-intrusive, rainbow schlieren deflectometry technique was used for quantitative measurements of helium concentrations in buoyant and non-buoyant jets. Results in a steady jet show that the radial growth of the jet shear layer in Earth gravity is hindered by the buoyant acceleration. The jet in microgravity was 30 to 70 percent wider than that in Earth gravity. The microgravity jet showed typical growth of a constant density jet shear layer. In case of a self-excited helium jet in Earth gravity, the flow oscillations continued as the jet flow adjusted to microgravity conditions in the drop tower. The flow oscillations were however not present at the end of the drop when steady microgravity conditions were reached.

Agrawal, Ajay K.↗