Search NASA⌕ Search

SEARCH · Search NASA

Results for “global 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Conservative Finite-Volume Based Interface-Tracking Algorithm Using the Signed Distance Function

Methods for tracking an interface between two fluid phases are developed to ensure desired fluid properties, conservation, and stability are preserved in a finitevolume (FV) discretization. Here, the interface is tracked using a level set method where the signed distance function implicitly defines the interface. Marching methods are used to evaluate the value of the signed distance function, including a novel initialization method to initialize any implicit function to the signed distance function around sharp corners in the level set. Global conservation and consistency with a set of governing equations is enforced by a compression coefficient that measures the volumetric compression or expansion due to inaccuracies in the level set evaluation. A redistribution method is integrated into the volume correction to eliminate the small-cell instability while maintaining global conservation. This suite of methods is implemented and tested using static uniform velocity, and potential flow cases with multiple interface geometries. Results show these methods achieve up to second order accuracy, and are conservative. The application for these methods is intended to track the interface of a 3D printing filament in a finite-volume discretization of the all-speed Navier-Stokes equations.

42 ENGINEERING↗
DOE OSTI2020

Eigenmode analysis of the sheared-flow Z-pinch

Experiments have demonstrated that a Z-pinch can persist for thousands of times longer than the growth time of global magnetohydrodynamic (MHD) instabilities such as the m=0 sausage and m=1 kink modes. These modes have growth times on the order of ta=a/vi, where vi is the ion thermal speed and a is the pinch radius. Axial flows with duz/dr ≲ vi/a have been measured during the stable period, and the commonly accepted theory is that this amount of shear is sufficient to stabilize these modes as predicted by numerical studies using the ideal MHD equations. However, these studies only consider specific equilibrium profiles that typically have a modest magnitude for the logarithmic pressure gradient, qP≡d ln P/d ln r, and may not represent experimental conditions. Linear stability of the sheared-flow Z-pinch is studied here via a direct eigen-decomposition of the matrix operator obtained from the linear ideal MHD equations. Several equilibrium profiles with a large variation of qP are examined. Considering a practical range of k, 1/3 ≲ ka ≲ 10, it is shown that the shear required to stabilize m=0 modes can be expressed as duz/dr≥Cγ0/(ka)α. Here, γ0=γ0(ka) is the profile-specific growth rate in the absence of shear, which scales approximately with |qP|. Both C and α are profile-specific constants, but C is order unity and α≈1. It is further demonstrated that even a large value of shear, duz/dr=3vi/a, is not sufficient to provide linear stabilization of the m=1 kink mode for all profiles considered. This result is in contrast to the currently accepted theory predicting stabilization at much lower shear, duz/dr=0.1vi/a, and suggests that the experimentally observed stability cannot be explained within the linear ideal-MHD model.

Angus, J. R. (ORCID:0000000314740002)↗
DOE OSTI2024

Global linear and nonlinear gyrokinetic simulations of tearing modes

Abstract To better understand multi-scale interactions between global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code Gene is modified to describe global tearing mode instability via a shifted Maxwellian (SM) distribution consistent with experimental equilibria. The implementation of the SM is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes ( n = 5, n = 6) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the m = 2 tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.

Physics↗
DOE OSTI2023

Disruption Mitigation Solutions for Long Pulse Tokamaks (Final Technical Report)

The project "Disruption Mitigation Solutions for Long Pulse Tokamaks" was a collaborative effort between ORNL, GA, and UCSD. The main goal of the project was to study the injection of shattered cryogenic pellets (SPI) as a method of reducing wall damage from global tokamak plasma instabilities called disruptions. This work was primarily carried out in the JET tokamak in the UK. The team already had extensive experience studying PSI in the medium sized DIII-D tokamak in the US, but adding data from the larger JET tokamak enabled more confident extrapolation to the even larger planned ITER tokamak in France. Each of the collaborating institutions focused on different aspects of SPI disruption mitigation. The UCSD collaboration focused on the volume recombination of post-disruption runaway electron (RE) beams by massive hydrogen injection. Disruptions can sometimes result in the formation of toroidal (ring shaped) beams of highly energetic (relativistic) electron beams, which can cause serious local wall damage, even in present-day tokamaks. A promising method to reduce wall damage from the RE beams is to inject massive amounts of H 2 or D 2 into the beams using SPI. although the precise method by which the injected H 2 reducses RE wall damage is still being investigated, it seems clear that a necessary and required step for the reduction to occur is for the H 2 injection to cause volume recombination (free electron density reduction) of the background cold plasma surrounding the fast electron beam. The UCSD collaboration focused on answering the question of if this recombination could be achievable in the future reactor-size tokamak ITER using realistic levels of H 2 injection. A model capturing the most essential physics was developed and was constrained using both DIII-D and JET data. This model was then used to extrapolate to ITER. It was predicted that RE beam recombination should be achievable in ITER, although in the worst-case scenari (a 10 MA Re beam in ITER), recombination is barely achieved within the maximum allowed levels of H 2 injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hydrostatic and Non‐Hydrostatic Baroclinic Instability in the Dynamical Core of the DOE Global Climate Model

Abstract The dynamical core of the Department of Energy global climate model is used to understand the role of non‐hydrostatic dynamics in the simulation of dry and moist baroclinic waves. To reduce computational cost, the Diabatic Acceleration and REscaling approach is adopted. Scale analysis and numerical simulations suggest that the model solution is not distorted by the change of spherical metric terms due to the change of Earth radius, neither are the inter‐scale interactions strongly altered as indicated by the analysis of spectral flux. Compared with hydrostatic simulations, the onset of baroclinic instability is delayed under non‐hydrostatic dynamics as the associated weaker vertical motions tend to increase the critical wavelength and narrow the range of waves that can be baroclinically unstable. During the development of baroclinic waves, non‐hydrostatic dynamics tends to induce vertical motions in the upper troposphere, accelerating the eastward propagation of upper‐level ridges through their impact on local vorticity tendency. These processes reduce the westward tilt of the vertical ridge axes and suppress the conversion of mean flow available potential energy to eddy kinetic energy, leading to weaker baroclinic eddies than in hydrostatic simulations. The contrasts between hydrostatic and non‐hydrostatic settings hold in supplemental experiments that vary the background flow Rossby number, the amount of water vapor content and vertical resolution. We also find that mesoscale and smaller‐scale activities can be considerably under‐represented when the vertical resolution is limited to that typically used in global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Disruption event characterization and forecasting in tokamaks

Disruption prediction and avoidance is a critical need for next-step tokamaks, such as ITER. Disruption Event Characterization and Forecasting (DECAF) research fully automates analysis of tokamak data to determine chains of events that lead to disruptions and to forecast their evolution allowing sufficient time for mitigation or complete avoidance of the disruption. Disruption event chains related to local rotating or global magnetohydrodynamic (MHD) modes and vertical instability are examined with warnings issued for many off-normal physics events, including density limits, plasma dynamics, confinement transitions, and profile variations. Along with Greenwald density limit evaluation, a local radiative island power balance theory is evaluated and compared to the observation of island growth. Automated decomposition and analysis of rotating tearing modes produce physical event chains leading to disruptions. A total MHD state warning model comprised of 15 separate criteria produces a disruption forecast about 180 ms before a standard locked mode detector warning. Single DECAF event analyses have begun on KSTAR, MAST, and NSTX/-U databases with thousands of shot seconds of device operation using from 0.5 to 1 × 10 6 tested sample times per device. An initial multi-device database comparison illustrates a highly important result that plasma disruptivity does not need to increase as βN increases. Global MHD instabilities, such as resistive wall modes (RWMs), can give the briefest time period of warning before disruption compared to other physics events. In an NSTX database with unstable RWMs, the mode onset, loss of boundary and current control, and disruption event warnings are found in all cases and vertical displacement events are found in 91% of cases. An initial time-dependent reduced physics model of kinetic RWM stabilization created to forecast the disruption chain predicts instability 84% of the time for experimentally unstable cases with a relatively low false positive rate. Instances of the disruption event chain analysis illustrate dynamics including H–L back transitions for rotating MHD and global RWM triggering events. Furthermore, disruption warnings are issued with sufficient time before the disruption (on transport timescales) to potentially allow active profile control for disruption avoidance, active mode control, or mitigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A reduced model for edge localized mode control by supersonic molecular beam injection and pellet injection

We develop a diffusive, bistable, tri-unstable cellular automata (CA) model to study the dynamics of H-mode pedestal with edge localized modes (ELMs) and their control by supersonic molecular beam injection (SMBI) and pellet injection (PI). It is shown that the new CA model can reproduce the key features of H-mode pedestals with various types of ELM, including Type-I ELM. SMBI and PI are modeled as additional grain injections into pedestal with varying degrees of injected materials and profiles. It is found that H-mode pedestal responds to SMBI differently depending on the baseline fueling. If the baseline fueling is large enough to allow Type-I ELM, SMBI enhances large transport avalanches caused by ballooning instabilities. These avalanches prevent the total pedestal current from reaching the boundary for peeling instability. On the other hand, if the baseline fueling is low to avoid Type-I ELM, SMBI enhances small scale avalanches, which prevent the pedestal from growing to profiles globally vulnerable to ballooning instabilities. These imply that SMBI can mitigate different types of ELM by converting them to more benign types. From CA modeling of pellet injection, it is shown that Type-I ELM can be triggered by pellet injection with sufficient strength and depth. Scanning the frequency of pellet injection, it is found that a maximum efficiency of pellet pacing is achieved when the injection frequency is approximately ten times the natural frequency of Type-I ELM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A generalized effective potential for differentially rotating plasmas

Global stability of differentially rotating plasma is investigated using a generalized effective potential. We first, for a current-free system, obtain a general form of an effective potential in terms of the free energies of global curvature and gradients of rotation for non-axisymmetric disturbances. We then examine the stability of differentially rotating disks for several rotation profiles and present the associated effective potential for the onset of these instabilities in the MHD regime. In particular, results for global axisymmetric magnetorotational instability as well as local and global non-axisymmetric modes are presented. The latter constitute two distinct non-axisymmetric modes, a high frequency local MRI and a global low-frequency non-axisymmetric mode (the magneto-curvature mode, introduced in Ebrahimi and Pharr [Astrophys. J. 936, 145 (2022)]), confined either between two Alfvénic resonances or an Alfvénic resonance and a boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nonlinear control of safety factor gradient in tokamaks using spatially variable electron cyclotron current drives

Active control of plasma properties may be necessary to achieve stable operation of next-generation tokamaks over large time scales. Such control algorithms can regulate the plasma properties to avoid the onset of magnetohydrodynamic (MHD) instabilities. For instance, the global and local properties of the safety factor profile are linked to the onset of neoclassical tearing modes (NTMs). This work proposes a model-based control approach for deterring/delaying NTMs through active modulation of one of the safety factor properties - the gradient of the safety factor profile at a particular rational safety factor surface. In particular, a novel control-oriented model for the local safety factor gradient is developed. The nonlinear control model is governed by a nonautonomous ordinary differential equation that accounts for a given rational safety factor surface’s spatial variation over time. Further, to improve the controllability of the spatially evolving parameter, the control model treats ECH&CD positions, along with noninductive powers, as controllable variables. A nonlinear control algorithm based on feedback linearization with optimization is synthesized to achieve the objective of regulating the safety factor gradient around a given target. The proposed algorithm allocates optimal ECH&CD positions, in addition to auxiliary powers, at each time instant as the rational safety factor surface drifts to locations with a low control authority. Stability guarantees of the proposed control law are also discussed in this work. The proposed algorithm is tested for a DIII-D tokamak scenario in nonlinear simulations carried out using the Control Oriented Transport SIMulator (COTSIM). Both fixed and moving ECH&CD cases are studied, and their outcomes are compared. Simulation results demonstrate that enthusiastic regulation of the safety factor gradient can be achieved during the ramp-up and flat-top phases of tokamak operation in both fixed and moving ECH&CD cases. However, real-time updates of ECH&CD positions can prevent the saturation of auxiliary powers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Measurement of pressure gradients near the interface in the viscous fingering instability

The viscous fingering instability, which forms when a less-viscous fluid invades a more-viscous one within a confined geometry, is an iconic system for studying pattern formation. For both miscible and immiscible fluid pairs the growth dynamics change after the initial instability onset and the global structures, typical of late-time growth, are governed by the viscosity ratio. Here we introduce an experimental technique to measure flow throughout the inner and outer fluids. This probes the existence of a new length scale associated with the local pressure gradients around the interface and allows us to compare our results to the predictions of a previously proposed model for late-time finger growth. Published by the American Physical Society 2024

Physics↗
DOE OSTI2022

Troubles with global monopoles in quantum gravity

A theory with a global [Formula: see text] symmetry broken at a scale [Formula: see text] admits topological configurations: global monopoles realized by Goldstone fields winding around the core of false vacuum. One may expect them to behave as heavy, big composite objects, difficult to make and therefore mostly harmless. However, after gravity is turned on, as long as Equivalence Principle holds, one finds that global monopoles have a negative mass, [Formula: see text], where [Formula: see text] is the field theory coupling. This could catalyze an instability in the space of a global monopole, by the production of additional global monopole–antimonopole pairs along with normal particles, leading to energy production ab nihilo, if the pair energy is dominated by their negative rest masses. In a theory with unbroken local Poincaré symmetry, this could lead to a divergent ‘decay rate’ of the global monopole configurations.

Astronomy & Astrophysics↗
DOE OSTI2021

Normal and anomalous self-healing mechanism of crystalline calcium silicate hydrates

The origin of different stability of crystalline calcium silicate hydrates was investigated. The tobermorite crystal has been used as an analog of cement hydrate that is being mostly manufactured material on earth. Normal tobermorite is thermally unstable and transforms to amorphous at low pressure. Meanwhile, anomalous tobermorite with high Al content does not significantly transform under high pressure or high temperature. Conducted X-ray absorption spectroscopy explains the weak stability of normal tobermorite which was originally hypothesized by the role of zeolitic Ca ions in the cavities of silicate chains. Atomic simulations reproduced the experimentally observed trend of pressure behavior once the ideal structures were modified to account for the Al content as well as the chain defects. The simulations also suggested that the stability of tobermorite under stress could be rationalized as a self-healing mechanism in which the structural instabilities were accommodated by a global sliding of the CaO layer.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Artificial Intelligence/Deep Learning FRNN Software for Prediction & Real-Time Control of DIII-D Plasma Control System (PCS)

This collaborative project integrated an improved version of the Artificial Intelligence/Deep Learning FRNN prediction and control software into the real-time DIII-D PCS (plasma control system). A key AI/DL software challenge is to build a modern high-performance computing (HPC) enabled “synthetic plasma simulator” capable of carrying out HPC-driven real-time plasma control applications. This involves development of a deep learning framework to train the surrogate model for a first-principles-based instability analysis simulator (“SGTC”) derived from the global gyrokinetic code GTC. The role of SGTC is to provide accurate and detailed plasma instability information from a real-time AI-based simulator capability to complement the deep learning prediction and control from experimentally-measured signals, such as ECE Imaging, supplemented by synthetic SGTC-ECEI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Magnetogenesis from early structure formation due to Yukawa forces

Yukawa interactions can mediate relatively long-range attractive forces between fermions in the early universe. Such a globally attractive interaction creates an instability that can result in the growth of structure in the affected species even during the radiation dominated era. The formation and collapse of fermionic microhalos can create hot fireballs at the sites of the collapsing halos which inject energy into the cosmic plasma. In this paper we study a new phenomena which can take place in such models. We show that the injected energy can be partially converted into primordial magnetic fields and we estimate the correlation scale and the power spectrum of these fields. We show that they may be the seeds of the observed astrophysical magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Advection versus diffusion in Richtmyer-Meshkov mixing

The Richtmyer-Meshkov (RM) instability is one of the most severe degradation mechanisms for inertial confinement fusion (ICF), and mitigating it has been a priority for the global ICF effort. In this Letter, the instability's ability to atomically mix is linked to its background decay of residual turbulent energy. Here, we show how recently derived inequalities from the mathematical theory of PDEs constrain the evolution. A model RM process at leading order may diffusively mix or retain imprints of its initial structures indefinitely, depending on initial conditions, and there exists a theoretical range of zero-mixing for certain values of parameters. The results may apply to other systems resembling scalar transport in decaying turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Nonideal stability analysis of differentially rotating plasmas with global curvature effects

The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.

Haywood, Alexander [Princeton Univ., NJ (United St↗