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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dynamics and chemical mode analysis of plasma thermal-chemical instability

The stability of the weakly ionized plasma and the transition from a stable homogeneous discharge to unstable filaments play an important role in gas laser physics, plasma-assisted combustion, chemical reforming, and material synthesis. Here, theoretical stability analysis and thermal-chemical mode analysis were performed to understand the mechanism of plasma thermal-chemical instability by using a zero-dimensional plasma system with both simplified and detailed chemical kinetics of H 2 /O 2 /N 2 mixtures. The plasma dynamic and kinetic models accounted for multiple physical mechanisms in the chemically-reactive weakly ionized plasma, including ionization, attachment/detachment, recombination, vibrational and electronic energy relaxation, convective and diffusive species/heat removal, Joule heating, and detailed chemical kinetics. An analytical criterion and the explosive mode species/temperature pointers were formulated while the representative active species were identified for different thermal-chemical modes. The results showed that in addition to the classical thermal-ionization mechanism, various chemical modes from chemical heat imbalance and elementary kinetics significantly modified the time dynamics and the stability of the weakly ionized plasma. The present analysis provides insights and guidance to control plasma instability using chemical kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Nonlocal Magneto-curvature Instability in a Differentially Rotating Disk

A global mode is shown to be unstable to nonaxisymmetric perturbations in a differentially rotating Keplerian disk containing either vertical or azimuthal magnetic fields. In an unstratified cylindrical disk model, using both global eigenvalue stability analysis and linear global initial-value simulations, it is demonstrated that this instability dominates at strong magnetic fields where local standard magnetorotational instability (MRI) becomes stable. Unlike the standard MRI mode, which is concentrated in the high flow shear region, these distinct global modes (with low azimuthal mode numbers) are extended in the global domain and are Alfvén-continuum-driven unstable modes. As its mode structure and relative dominance over MRI are inherently determined by the global spatial curvature as well as the flow shear in the presence of a magnetic field, we call it the magneto-curvature (magneto-spatial-curvature) instability. Consistent with the linear analysis, as the field strength is increased in the nonlinear simulations, a transition from MRI-driven turbulence to a state dominated by global nonaxisymmetric modes is obtained. This global instability could therefore be a source of nonlinear transport in accretion disks at a higher magnetic field than predicted by local models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

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↗
DOE OSTI2020

The Role of Multiscale Interaction in Tropical Cyclogenesis and Its Predictability in Near-Global Aquaplanet Cloud-Resolving Simulations

Tropical cyclogenesis (TCG) is a multiscale process that involves interactions between large-scale circulation and small-scale convection. A near-global aquaplanet cloud-resolving model (NGAqua) with 4-km horizontal grid spacing that produces tropical cyclones (TCs) is used to investigate TCG and its predictability. This study analyzes an ensemble of three 20-day NGAqua simulations, with initial white-noise perturbations of low-level humidity. TCs develop spontaneously from the northern edge of the intertropical convergence zone (ITCZ), where large-scale flows and tropical convection provide necessary conditions for barotropic instability. Zonal bands of positive low-level absolute vorticity organize into cyclonic vortices, some of which develop into TCs. A new algorithm is developed to track the cyclonic vortices. A vortex-following framework analysis of the low-level vorticity budget shows that vertical stretching of absolute vorticity due to convective heating contributes positively to the vorticity spinup of the TCs. A case study and composite analyses suggest that sufficient humidity is key for convective development. TCG in these three NGAqua simulations undergoes the same series of interactions. The locations of cyclonic vortices are broadly predetermined by planetary-scale circulation and humidity patterns associated with ITCZ breakdown, which are predictable up to 10 days. Whether and when the cyclonic vortices become TCs depend on the somewhat more random feedback between convection and vorticity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

First implementation of gyrokinetic exact linearized Landau collision operator and comparison with models

Gyrokinetic simulations are fundamental to understanding and predicting turbulent transport in magnetically confined fusion plasmas. Previous simulations have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius (FLR) effects. We have implemented the linearized Fokker–Planck collision operator with exact field-particle terms and full FLR effects in a gyrokinetic code (GENE). The new operator, referred to as “exact” in this paper, allows the accuracy of model collision operators to be assessed. The conservative Landau form is implemented because its symmetry underlies the conservation laws and the H-theorem, and enables numerical methods to preserve this conservation, independent of resolution. The implementation utilizes the finite-volume method recently employed to discretize the Sugama collision model in GENE, allowing direct comparison between the two operators. Results show that the Sugama model appears accurate for the growth rates of trapped electron modes (TEMs) driven only by density gradients, but appreciably underestimates the growth rates as the collisionality and electron temperature gradient increase. The TEM turbulent fluxes near the nonlinear threshold using the exact operator are similar to the Sugama model for the n e = d ln T e /d ln n e = 0 case, but substantially larger than the Sugama model for the n e = 1 case. The FLR effects reduce the growth rates increasingly with wavenumber deepening a “valley” at the intermediate binormal wavenumber as the unstable mode extends from the TEM regime to the electron temperature gradient instability regime. Application to the Hinton–Rosenbluth problem shows that zonal flows decay faster as the radial wavenumber increases and the exact operator yields weaker decay rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Observation of Nonaxisymmetric Standard Magnetorotational Instability Induced by a Free-Shear Layer

The standard magnetorotational instability (SMRI) with a magnetic field component parallel to the rotation axis is widely believed to be responsible for the fast accretion in astronomical disks. In conventional base flows with a Keplerian profile or an ideal Couette profile, most studies focus on axisymmetric SMRI, since excitation of nonaxisymmetric SMRI in such flows requires a magnetic Reynolds number (Rm) more than an order of magnitude larger. Here, we report that, in a magnetized Taylor-Couette flow, nonaxisymmetric SMRI with an azimuthal mode number m = 1 can be triggered by a free-shear layer in the base flow at Rm ≳1, the same threshold as for axisymmetric SMRI. Global linear analysis reveals that the free-shear layer reduces the required Rm, possibly by introducing an extremum in the vorticity of the base flow. Nonlinear simulations validate the results from linear analysis and confirm that a novel instability recently discovered experimentally [Wang et al., Nat. Commun. 13, 4679 (2022)] is the nonaxisymmetric m =1 SMRI. Further, our finding has astronomical implications as free-shear layers are ubiquitous in celestial systems, such as the disk-star boundary layer, the solar tachocline, and the edge of planet-opened gaps in protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Flow Channels and the Generation of Alfvénic Turbulence Along Storm‐Time Inner Magnetospheric Field‐Lines

Abstract A feature of Earth's storm‐time magnetosphere outside the plasmapause is the occurrence of broad‐spectrum Alfvénic fluctuations. In this letter observations from the Van Allen Probes are compared with 3‐D fluid‐kinetic simulations of an evolving convective flow channel to investigate the mechanisms generating the observed spectrum. It is shown how narrow channels of fast convection are unstable to the Kelvin‐Helmholtz instability which on closed field‐lines initiates a cascade to small scales. Sustained driving of the flow combined with reflection from the topside ionosphere leads to the generation of an intensified spectrum of electromagnetic structures having similar spectral and morphological characteristics to those observed. This process couples enhanced magnetospheric convection to kinetic scale electromagnetic fluctuations that drive particle transport, scattering and energization through the outer radiation belt and ring current during geomagnetic storms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Turbulence theories and statistical closure approaches

When discussing research in physics and in science more generally, it is common to ascribe equal importance to the three components of the scientific trinity: theoretical, experimental, and computational studies. This review will explore the future of modern turbulence theory by tracing its history, which began in earnest with Kolmogorov’s 1941 analysis of turbulence cascade and inertial range [A.N. Kolmogorov, Dokl. Akad. Nauk SSSR, 30, 299, (1941); 32, 19, (1941)]. The 80th Anniversary of Kolmogorov’s landmark study is a welcome opportunity to survey the achievements and evaluate the future of the theoretical approach of turbulence research. Over the years, turbulence theories have been critically important in laying the foundation of our understanding of the nature of turbulent flows. In particular, the Direct Interaction Approximation (DIA) [R.H. Kraichnan, J. Fluid Mech., 5, 497 (1959)] and its subsequent development, known as the statistical closure approach, can be identified as perhaps the most profound single advancement. The remarkable success of the statistical closure has furnished a platform to study such essential concepts as the energy transfer process and interacting scales, and the roles of the straining and sweeping motions. More recently, the quasi-Lagrangian formulation of V. L’vov & I. Procaccia and Kraichnan’s solvable passive scalar model provided powerful ways to explore another fundamental aspect of turbulent flows, the phenomena of intermittency, and the associated anomalous scaling exponents. In the meantime, the theory of fluid equilibria has been developed to describe the large-scale structures that can emerge from turbulent cascades of two-dimensional and geophysical flows at a later time. And yet, despite all these successes, analytical treatments suffer from mathematical complexities. As a result, the utility of theoretical approaches has been limited to relatively idealized flows. On the other hand, in recent decades, computational abilities and experimental facilities have reached an unprecedented scale. Looking beyond the horizon, the imminent deployment of exascale supercomputers will generate complete datasets of the entire flow field of key benchmark flows, allowing researchers to extract additional measurements concerning fully developed, complex turbulent flow fields far beyond those available from the statistical closure theories. Some other developments that could potentially influence the future course of turbulence theories include the advancement of machine learning, artificial intelligence, and data science; likely disruptions arising from the advent of quantum computation; and the increasingly prominent role of turbulence research in providing more accurate climate scientific data. Finally, turbulence theorists can leverage these developments by asking the right questions and developing advanced, sophisticated frameworks that will be able to predict and correlate vast amounts of data from the other two components of the trinity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Modeling Hydrodynamic Instabilities, Shocks, and Radiation Waves in High Energy Density Experiments [Dissertation]

This thesis presents the computational design, modeling, and analysis of three experiments in high-energy density physics (HEDP), all of them concerning fundamental radiation flows. The first experiment is a laboratory astrophysics experiment to investigate the role of the Kelvin-Helmholtz instability (KHI) in the process of galactic filaments supplying gas to galactic halos. The achieved goal was to provide a first study in which the role of the instability is maximal and predict behavior in future iterations of an experiment accessing a more radiative regime where the role of the instability is stifled. This experiment would help answer how certain galaxies are able to grow so rapidly and produce many stars as the KHI limits this process. The second experiment, COAX, is a radiation flow experiment with a novel spectroscopy diagnostic configuration, designed to spatially measure the temperature of a radiation wave as it travels down a doped foam. A key result of this work was the development of a synthetic spectroscopy application and application of modern spectroscopy comparison techniques to provide our first temperature reconstructions from the experimental data. This experimental platform serves as the launching ground for a number of new experiments that vary the basic premise and thus is foundational to our ongoing research. The final experiment is a full integration of modeling, design, and theoretical development for the Radishock experiment. This experimental platform studies the head-on collision of a radiation wave with a counter-propagating shock, and like COAX, uses spectroscopy to diagnose and detect the interaction. My research analyzes the successful shots, indicating aspects of successful detections and suggests improvements to future iterations of the design.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

High resolution numerical simulations of methane pool fires using adaptive mesh refinement

The ability to accurately predict the structure and dynamics of pool fires using computational simulations is of great interest in a wide variety of applications, including accidental and wildland fires. However, the presence of physical processes spanning a broad range of spatial and temporal scales poses a significant challenge for simulations of such fires, particularly at conditions near the transition between laminar and turbulent flow. Here, in this study, we examine the transition to turbulence in methane pool fires using high-resolution simulations with multi-step finite rate chemistry, where adaptive mesh refinement (AMR) is used to directly resolve small-scale flow phenomena. We perform three simulations of methane pool fires, each with increasing diameter, corresponding to increasing inlet Reynolds and Richardson numbers. As the diameter increases, the flow transitions from organized vortex roll-up via the puffing instability to much more chaotic mixing associated with finger formation along the shear layer and core collapse near the inlet. These effects combine to create additional mixing close to the inlet, thereby enhancing fuel consumption and causing more rapid acceleration of the fluid above the pool. We also make comparisons between the transition to turbulence and core collapse in the present pool fires and in inert helium plumes, which are often used as surrogates for the study of buoyant reacting flows.

42 ENGINEERING↗
DOE OSTI2022

Absolute vs Convective Instabilities and Front Propagation in Lipid Membrane Tubes

Here we analyze the stability of biological membrane tubes, with and without a base flow of lipids. Membrane dynamics are completely specified by two dimensionless numbers: the well-known Föppl-von Kármán number Γ and the recently introduced Scriven-Love number SL, respectively quantifying the base tension and base flow speed. For unstable tubes, the growth rate of a local perturbation depends only on Γ, whereas SL governs the absolute versus convective nature of the instability. Furthermore, nonlinear simulations of unstable tubes reveal an initially localized disturbance result in propagating fronts, which leave a thin atrophied tube in their wake. Depending on the value of Γ, the thin tube is connected to the unperturbed regions via oscillatory or monotonic shape transitions-reminiscent of recent experimental observations on the retraction and atrophy of axons. We elucidate our findings through a weakly nonlinear analysis, which shows membrane dynamics may be approximated by a model of the class of extended Fisher-Kolmogorov equations. Our study sheds light on the pattern selection mechanism in axonal shapes by recognizing the existence of two Lifshitz points, at which the front dynamics undergo steady-to-oscillatory bifurcations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Roles of Mn and Co in the Air Synthesizability of Layered Oxide Cathodes for Lithium-Based Batteries

High-nickel layered oxides (LiNi 1-x-y Mn x Co y O 2 ) are the prevailing cathode materials for high-energy-density lithium-based batteries, but they are plagued with deleterious surface air instabilities stemming from residual lithium formation. These issues severely hinder mass production as cathode calcination is limited to a flowing oxygen atmosphere, which entails high manufacturing costs as opposed to simpler and more economical air calcination. Here, while higher Ni contents are known to worsen air instabilities, the influence of Mn and Co contents on impacting these phenomena are less elucidated. We herein present the synthesis in ambient air and flowing oxygen atmospheres of three cathode variants with the same Ni contents, but varying Mn and Co contents: LiNi 0.7 Mn 0.3 O 2 , LiNi 0.7 Mn 0.15 Co 0.15 O 2 , and LiNi 0.7 Co 0.3 O 2 . It is found that the critical parameter influencing the air stability of the cathodes is the average Ni oxidation state, which is greatly dependent on the Mn and Co contents. Substitution of Mn for Ni drives down the Ni oxidation state as Mn exists as Mn4 + and reduces surface residual lithium formation, which vastly improves the overall air stability and, therefore, the synthesizability in air, but with a penalty of lowered capacity. In contrast, substitution of Co for Ni maintains Ni 3+ as Co exists as Co 3+ , offering increased initial capacity, but worsens the air stability and cyclability as the driving force for residual lithium formation and surface reactivity is increased.

25 ENERGY STORAGE↗
DOE OSTI2025

An experimental and computational study of thin-layer Rayleigh–Taylor instability development during deceleration with and without an externally applied magnetic field

The importance of mitigating the Rayleigh–Taylor instability (RTI) in inertial confinement fusion (ICF) is critical to successfully achieve high gain fusion yield. Consequently, understanding the seed mechanisms of RTI and the potential evolution of RTI in ICF relevant conditions is crucial. Single feature perturbations consistently demonstrate non-linear RTI evolution, for which an experimental platform on OMEGA-EP is developed. Manufacturing defects introduced into the target design require exploration of unanticipated changes to RTI development and an identification of targets that will still render quantifiable physics results. Consequently, it is presented that the inherent 3D nature of experimental targets necessitates 3D modeling for accurate design work and predictive modeling of experimental targets, especially when high resolution imaging diagnostics, like Fresnel Zone Plates, are utilized. A study of the morphology of the RTI evolution due to changing initial conditions and the presence of an externally applied magnetic field are also explored. Experimental data show thin-layer RTI morphology comparable to resistive magneto hydrodynamic 3D results. A discussion on the impacts of an externally applied magnetic field makes the case for continued efforts to observe a magnetic field's impact on RTI morphology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Experimental, numerical and analytical evaluation of j × B -thrust for fast-liquid-metal-flow divertor systems of nuclear fusion devices

Abstract Divertor systems of fusion devices are exposed to intense heat loads from plasmas, which degrade solid plasma-facing components. Fast liquid metal (LM) flow divertors may be more advantageous for this purpose but have risk of piling due to intense magnetohydrodynamic (MHD) drag. However, severe deceleration of the flow could be countered with the injection of currents that are transverse to external magnetic fields, allowing to thrust the flow with j × B (Lorentz) forces. Given that the injection of currents as an approach to propel LM-divertor flows has remained experimentally understudied, this article focuses on the evaluation of j × B -thrust and finding its drawbacks. j × B -thrust was experimentally tested with free-surface-LM flows, a vertical magnetic field and an externally applied current. Experiments were reviewed with a theoretical model, showing agreement in the trends of theory and experiments. Full 3D-MHD-free-surface-flow simulations were also performed with FreeMHD and confirmed the sensitivity to unstable flow behavior in LM systems when applying external currents. Furthermore, excessive power requirements are expected for the implementation of j × B -thrust at the reactor scale, making these systems inefficient for commercial devices. This paper evidences that the simple operation of a LM-flow divertor with j × B -thrust, without any of the instabilities caused from reactor plasmas or parasitic currents, already presents intrinsic challenges.

Magnetohydrodynamics↗
DOE OSTI2021

White paper on NIF Discovery Science and frontier regimes of HEDP

Over the past two decades, experimental capabilities at high energy density (HED) facilities such as the National Ignition Facility (NIF), Omega, and Omega EP lasers, and the Z pulsed power facility are able to make precision measurements of fundamental quantities such as equations of state (EOS), material phase, opacities, and ionization levels Z. Flow based quantities are also being studied, such as hydrodynamic instabilities, evolution into turbulence, turbulent dynamo magnetic field amplification, magnetic reconnection; plasma instabilities, plasma wakefield particle acceleration; and relativistic plasma phenomena, such as (e + ,e - ) pair creation. Experimental regimes can be created that reproduce conditions relevant to planetary and exoplanet interiors, stellar interiors, the interiors of brown dwarfs, and the envelope conditions of white dwarf stars. Diagnostic techniques have been developed hand in hand with these new experimental regimes that can be accessed. In this white paper, a brief summary of key experimental achievements from the NIF Discovery Science (basic science) program will be presented; followed by a discussion of current challenges, needs, and potential paths forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗