Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “instability criteria”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Destabilization of binary mixing layer in supercritical conditions

Compressible mixing layer instabilities are of importance to a wide range of environmental and industrial applications. Past studies have focused on either ideal-gas or real-fluid thermodynamic regimes of single-species mixing layers. However, mixing layers of binary mixtures at supercritical conditions, commonly encountered in fuel injection systems, introduce additional complexities due to the added compositional degree-of-freedom. Moreover, the effect of strong variations in thermodynamic response functions across the Widom line on the binary mixing layer stability remains poorly understood. Thus, the objective of this study is to examine the coupling between the hydrodynamic instability and the real-fluid thermodynamics across the Widom line and its effects on the overall binary mixing layer dynamics. To this end, we develop a linear stability analysis of the full binary-species compressible transport equations coupled with the PC-SAFT equation of state. Analysis shows the existence of a novel instability mechanism that arises from juxtapositioning of the Widom-line transition and the hydrodynamic inflection point. This novel thermodynamically induced instability mechanism has the net effect of destabilizing the binary mixing layer at lowering supercritical conditions towards the critical pressure point. This is in contrast to previous stability analyses of supercritical single-species mixing layers, where increasing pressure destabilizes the flow due to its effect on reducing the density stratification. The discovered thermodynamically induced instability mechanism of binary mixing flows highlights the need for an extension of classical instability criteria to incorporate the effect of strong variations in the thermodynamic response functions across the Widom line on mixing layer instability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Theory of the magnetothermal instability in coronal plasma flows

The theory of the magnetothermal instability (MTI) is revisited through the lens of the stability of uniform systems. The linear stability analysis includes flow advection and Nernst transport. The instability criteria derived distinguish between the convective and the absolute nature of the perturbation growth. It is proven that, in the region where the Nernst and plasma blowoff velocities cancel, the MTI can be absolute and wave-packet perturbations grow in situ. This instability is mediated by the internal feedback between the Biermann battery and Righi–Leduc terms. Here, the analysis is extended to derive the dispersion relation for short-wavelength perturbations developing in nonuniform profiles with the application to coronal plasmas. It is found that the condition for MTI requires the net B-field convection velocity to be small at the isothermal sonic section, and the plasma conditions in this section govern the dynamics of the instability. Analysis of hydro-equivalent implosions suggests that unstable perturbations undergo more e-foldings of growth in larger-size targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λ c , that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λ c -based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Nonlinear harmonics coupled by parallel wave propagations in a time-dependent plasma flow

In a time-dependent flow, nonlinear harmonics can be excited by coupling between linear waves and flow-induced harmonic waves. Examining the dispersion relations and selection rules for the coupling, we investigate nonlinearly coupled harmonics for waves propagating along the magnetic field line in a magnetized plasma, as well as waves in an unmagnetized plasma. The coupled harmonics in a plasma flow are described by analytic dispersion relations and selection rules. This nonlinear coupling is corroborated by the particle-in-cell simulation. Here, the coupled-harmonics model describes a mechanism for the excitation of nonlinear harmonics from linear waves in a time-dependent flow. The spectral analysis of the dispersion relation provides a useful way to evaluate the spatiotemporal behavior of a plasma flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Dispersion relation and instability for an anisotropic nonuniform flowing plasma

A generalized formula for wave instability is developed for an anisotropic nonuniform plasma with finite flows and temperatures. Six-moment fluid equations are solved to give the analytic expression for wave instability in arbitrarily nonuniform plasmas. The analytic formula explicitly states the dependence of wave instability on the nonuniformities of number density, flow velocity, and anisotropic or isotropic pressure. The accuracy of the formalism is verified by a numerical calculation of implicit dispersion relations in complex Fourier space. The analysis shows that nonuniformity plays a critical role in plasma instability, while the flow velocity and anisotropic pressures determine the growth rate of the instability. Lastly, the instability diagram and associated instability criterion for anisotropy-driven instability are introduced as applications of the formalism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Weibel-like instability in magnetohydrodynamics

In magnetohydrodynamics (MHD), a density perturbation perpendicular to an electron temperature gradient generates a magnetic field around itself that acts to increase the perturbation, which can lead to instability. An MHD dispersion relation is obtained for perturbations perpendicular to a fixed electron temperature gradient with an initial in-plane magnetic field, including resistivity, viscosity, and the electrothermal coefficient. Instability occurs for sufficiently small electron temperature-gradient scale lengths determined by the ion collisionless skin depth. Both viscosity and resistivity are required to prevent growth at arbitrarily small spatial scales and to give a physical result for the fastest growing mode. The perpendicular electrothermal coefficient is only significant for a narrow range of low electron Hall parameters, causing a modest reduction in magnetic field growth and modifying the criteria for instability in the presence of viscosity. If the definition of the Weibel instabilityis extended to include all instabilities due to anisotropy in the electron velocity distribution, then this is a Weibel-like instability because an electron temperature gradient implies an anisotropic electron velocity distribution. Here, the implications for the formation of filaments in laser-produced plasmas and for the verification of MHD codes are considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Magnetic reconnection in 3D fusion devices: non-linear reduced equations and linear current-driven instabilities

Abstract Magnetic reconnection in 3D fusion devices is investigated. With the use of Boozer co-ordinates, we reduce the non-linear resistive magnetohydrodynamic equations in the limit of large aspect ratio and finite pressure fluctuations, to obtain a set of non-linear equations suitable for magnetic reconnection studies in stellarators. Magnetic flux unfreezing due to a finite electron mass is also considered. Equations that govern the linear regime and some of their general properties are given. We emphasise the role of magnetic geometry and identify how some aspects of stellarator optimisation could have an impact on reconnecting instabilities, in particular by exacerbating those enabled by electron inertia. The effect of 3D coupling on the linear reconnection rates and the mode structure is quantitatively addressed in the case in which the equilibrium rotational transform has one specific resonant location for which one mode can reconnect while coupled to an arbitrary number of non-resonant harmonics. The full problem is rigorously reduced to an equivalent cylindrical one, by introducing some geometrically modified plasma inertial and dissipative scales. The 3D scalings for the growth rates of reconnection instabilities and their destabilisation criteria are given.

Physics↗
DOE OSTI2021

Cooling and instabilities in colliding flows

ABSTRACT Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power-law cooling function. The characteristic length and time-scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently and rapidly, the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one- and three-dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Gyrokinetic simulations compared with magnetic fluctuations diagnosed with a Faraday-effect radial interferometer-polarimeter in the DIII-D pedestal

Abstract Experimental data on electromagnetic fluctuations in DIII-D, made available by the Faraday-effect radial interferometer-polarimeter (RIP) diagnostic Chen (2016 Rev. Sci. Instrum. 87 11E108), is examined in comparison with detailed gyrokinetic simulations using gyrokinetic electromagnetic numerical experiment (GENE). The diagnostic has the unique capability of making internal measurements of fluctuating magnetic fields ∫ n e δ B r d R ∫ n e d R . Local linear simulations identify microtearing modes (MTMs) over a substantial range of toroidal mode numbers (peaking at n = 15) with frequencies in good agreement with the experimental data. Local nonlinear simulations reinforce this result by producing a magnetic frequency spectrum in good agreement with that diagnosed by RIP. Simulated heat fluxes are in the range of experimental expectations. However, magnetic fluctuation amplitudes are substantially lower than the experimental expectations. Possible sources of this discrepancy are discussed, notably the fact that the diagnostics are localized at the mid-plane—the poloidal location where the simulations predict the fluctuation amplitudes to be smallest. Despite some discrepancies, several connections between simulations and experiments, combined with general criteria discriminating between potential pedestal instabilities, strongly point to MTMs as the source of the observed magnetic fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Behavior of solitary waves of coupled nonlinear Schrödinger equations subjected to complex external periodic potentials with odd-$\mathcal{PT}$ symmetry

In this work, we discuss the response of both moving and trapped solitary wave solutions of a two-component nonlinear Schrödinger system in 1 + 1 dimensions to an odd-$\mathcal{PT}$ external periodic complex potential. The dynamical behavior of perturbed solitary waves is explored by conducting numerical simulations of the nonlinear system and using a collective coordinate variational approximation. We present case examples corresponding to choices of parameter values and initial conditions involved therein. The results of the collective coordinate approximation are compared against numerical simulations where we observe qualitatively good agreement between the two. Unlike the case for a single-component solitary wave in a complex periodic $\mathcal{PT}$-symmetric potential, the collective coordinate equations do not have a small oscillation regime, and initially the height of the two components changes in opposite directions often causing instability. We find that the dynamic stability criteria we have used in the one-component case are a good indicator for the onset of dynamic instabilities in the present setup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

First-Principles Examination of Multiple Criteria of Organic Solvent Oxidative Stability in Batteries

Oxidative instability of the liquid electrolyte at or near battery cathode oxide surfaces has significant detrimental effects on batteries. Organic solvent molecules are often the fuel and precursors of such degradation processes, releasing electrons and protons that react with cathode oxides and electrolyte anions. These reactions contribute to cathode–electrolyte interphase (CEI) film formation, transition-metal ion dissolution, and phase transformation of the surface regions of the cathode. Here we apply density functional theory calculations to examine four criteria of oxidative stability (oxidation potential, hydrogen removal energies, and initial reactivity on two types of oxide facets) using four different solvent/additive molecules (ethylene carbonate, fluoroethylene carbonate, 1,3-dioxolane, and dimethyl ether). The ranking of molecular stability differs with each criterion. Surprisingly, the all-oxygen-terminated basal planes of layered oxides exhibit lower reaction barriers than spinel surface facets with exposed transition-metal cations, especially for ether solvents; the calculations also suggest basal planes contribute to the dissolution of transition-metal ions. Further, the structure–degradation relation complexity underscores the challenge of understanding the function of the CEI but also offers a guide to future degradation-mitigation strategies including facet engineering. Our predictions and models help establish a framework for future studies relevant to high-voltage conditions.

25 ENERGY STORAGE↗
DOE OSTI2025

Structurally Frustrated Remeika Phases: Tolerance Factor and Non-Kramers’ Ion Effects

Quantum materials that are found on the verge of structural, magnetic, and electronic instabilities are deep reservoirs for exotic phenomena. Our criteria to study the reported Ln 5 Ru 6 Sn 18 (Ln = Gd, Tb) are guided by the tolerance factor for Remeika phases─a class of intermetallic compounds with the general formula A 3 M 4 X 13 (where A = rare-earth element, M = transition metal, and X = tetrel), which can be viewed as structural analogues of pseudoperovskites (4 × ABX 3 ), with one of the tetrel atoms occupying the A-site. We have grown single crystals, up to 0.5 cm on the largest facet, of Gd 5 Ru 6 Sn 18 (a = 13.8052(12) Å) and Tb 5 Ru 6 Sn 18 (a = 13.7825(9) Å), which we expect to be on the verge of a formation instability. Bulk magnetic, heat capacity, and electrical transport measurements reveal unusual behavior that originates from the trivalent Gd and Tb ions, with deviations from typical metallic behavior and magnetic ordering that is short-range or disrupted in some other way. Furthermore, to better understand these phenomena, single-crystal neutron diffraction is investigated, which reveals short-range correlations among the partially occupied Tb and Sn sublattices, with 3D-ΔPDF analysis supporting direction-dependent local order. Taken together, these measurements show the utility of using tolerance factors in identifying materials that are likely to exhibit complex phenomena and uncover environments where a formation instability may strongly impact emergent bulk magnetic and electronic phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A passive model for the evolution of subgrid-scale instabilities in turbulent flow regimes

In this work, a zero-dimensional model for the development of turbulent mixing layers driven by Rayleigh–Taylor (RT) or Richtmyer–Meshkov (RM) instabilities is presented, for application in Reynolds-Averaged Navier–Stokes (RANS) simulations. The model approximates the behavior of the turbulent mixing layer at a passively advected tracer particle that sits on the material interface of interest in the simulation, and is used to provide initial conditions for the RANS model once the mixing layer width has developed to a size that may be resolved on the computational mesh. The model relies on approximations to the Besnard–Harlow–Rauenzahn (BHR) RANS closure equations and is thus consistent the turbulent development of the mixing layer, unlike previous models of this type which have focused on the early-time laminar development of RT or RM instabilities. These instabilities will often grow through an initially laminar state before developing into turbulence, and criteria for switching from an existing laminar model to the proposed turbulent model are considered. The proposed model reproduces the mixing layer development of RT and RM instabilities seen in mixed width resolving BHR RANS simulations to a reasonable degree of accuracy, and relaxes the need for the RANS simulation to resolve the mixed width as soon as the mixing layer becomes turbulent, which may be impractical in cases where the mixing layer grows through a wide range of scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Discoveries in Blast-Driven Turbulence of Astrophysical Relevance

The fluid mixing caused by variable-density instabilities is important in a wide variety of scenarios from ocean mixing and astrophysical phenomena to nuclear fusion techniques and atomic weapons. This thesis explores the mixing resulting from a specific instability known as the Blast Driven Instability (BDI). This work investigates the variable density mixing in an explosively driven environment due to the fluid instabilities at the material interfaces. Specifically, diverging Richtmyer-Meshkov (impulsive-acceleration environment) and Rayleigh-Taylor (variable-acceleration environment) instabilities (present in supernova and inertial confinement fusion) are studied using advanced high-speed diagnostics in carefully designed laboratory experiments. The BDI morphology is presented through a time development of Mie scattering images, and steps through the parameter space (varying density ratio and driver speed), highlighting the development of the structures that form during mixing. A scaling criterion is used to relate the two systems of vastly different spatiotemporal scales. Velocity fields in the BDI have been captured for the first time using the high temporal resolution PIV technique. Subsequent analysis of the dynamics of the instability from the velocity fields illustrates the distribution of kinetic energy, the transition to turbulence, and the characteristic growth of the instability are discussed. This study furthers understanding of how blast-driven instability pertains to supernova and inertial confinement fusion science. The morphology of the BDI has been characterized for the first time. This work steps through the parameter space covered in Mie scattering experiments, and how the different parameters contribute to development of structures and mixing. It also examines a scaling of the Atwood number for expanding predictive capabilities to other experimental conditions and simulations. The first collection of velocity fields acquired for the BDI are recorded, and subsequent analysis evaluating the distribution of kinetic energy throughout space and time for two density ratios from the overall parameter space, as well as the transition to turbulence, estimated from a Reynolds number calculated based on momentum mixing are all presented. This information is useful in advancing the development of models to predict physics of high energy density applications where experiments are not always readily available. This research has successfully demonstrated understanding for the time criteria defining regimes where the shock driven (Richtmyer-Meshkov instability) and the buoyancy driven (Rayleigh-Taylor instability) dominates through a parametric study of density variation (Atwood number) and driver speed (Mach number). All this furthers understanding of how the BDI pertains to SN and ICF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Description of longitudinal space charge effects in beams and plasma through dielectric permittivity

We develop a universal framework which allows quickly solve a wide class of problems for longitudinal space charge effects in beams and plasmas in cylindrical geometry. We introduce the longitudinal dielectric permittivity for the beam of charged particles, which describes its collective space charge response. The analysis yields an effective plasma frequency, which depends on the transverse geometry of the system. This dielectric permittivity mirrors the dielectric permittivity of plasma and matches the one dimensional (1D) expression once the transverse size of the beam is large. Several particle species can be included as additive terms describing susceptibility of each specie. The developed approach allows to study stability criteria for collective beam-beam and beam-plasma instabilities for arbitrary transverse distributions in particle densities.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Uniform-density Bose-Einstein condensates of the Gross-Pitaevskii equation found by solving the inverse problem for the confining potential

Here, in this work, we consider a “reverse-engineering” approach to construct confining potentials that support exact, constant density kovaton solutions to the classical Gross-Pitaevskii equation (GPE) also known as the nonlinear Schr¨odinger equation (NLSE). In the one-dimensional case, the exact solution is the sum of stationary kink and anti-kink solutions, i.e. a kovaton, and in the overlapping region, the density is constant. In higher dimensions, the exact solutions are generalizations of this wave function. In the absence of self-interactions, the confining potential is similar to a smoothed out finite square well with minima also at the edges. When self-interactions are added, a term proportional to ±gψ*ψ gets added to the confining potential and ±gM, where M is the norm, gets added to the total energy. In the realm of stability analysis, we find (linearly) stable solutions in the case with repulsive self-interactions which also are stable to self-similar deformations. For attractive interactions, however, the minima at the edges of the potential get deeper and a barrier in the center forms as we increase the norm. This leads to instabilities at a critical value of M (related to the number of particles in the BEC). Comparing the stability criteria from Derrick’s theorem and Bogoliubov-de Gennes analysis stability results, we find that both predict stability for repulsive self-interactions and instability at a critical mass M for attractive interactions. However, the numerical analysis gives a much lower critical mass. The numerical analysis shows further that the initial instabilities violate the symmetry x → -x assumed by Derrick’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Phonon thermodynamics and elastic behavior of GaAs at high temperatures and pressures

The phonons of wurtzite and zinc blende GaAs were calculated at simultaneously elevated temperature and pressure, and elastic constants were calculated as functions of pressure. Pressure caused instabilities of shorter-wavelength transverse acoustic modes in both wurtzite and zinc blende GaAs, causing them to fall to zero at 18 and 20GPa, respectively. Here, the Born stability criteria, which depend on elastic constants and only long wavelength phonons, therefore overestimated the pressure needed to induce instability at 0 K. At elevated temperatures, explicit anharmonicity pushes the onset of instability to higher pressures in both wurtzite and zincblende GaAs. Phonon linewidth and densities of states data showed that the quasiharmonic approximation failed to account for temperature-induced phonon frequency shifts, and the quasiharmonic approximation became less reliable at elevated pressure. In general, the number of three-phonon processes increased with pressure, thereby increasing the temperature-driven broadening of phonon spectral lineshapes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Method for Selecting Rotor Suspension Design Criteria

A method for selecting rotor suspension design criteria that allow rotors with internal friction to negotiate critical speeds and avoid subsynchronous internal friction instabilities is described. A rotor model of a test rotor with internal friction is created and used to calculate the test rotor’s critical speeds, mode shapes, Campbell plots, and the rotational speed corresponding to the stability boundary for a nominal set of suspension parameters. The test rotor model is then used to demonstrate the suspension design criteria selection method.

Damiano, Brian↗