Search NASA⌕ Search

SEARCH · Search NASA

Results for “destabilization”

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 991 records · Page 55

A study of resistive peeling–ballooning modes across spherical tokamaks

We investigate how non-ideal-magnetohydrodynamics (MHD) effects, in particular plasma resistivity, impact the peeling–ballooning stability thresholds in spherical tokamaks. This analysis follows the discovery of resistive kink-peeling modes in ELMing National Spherical Torus Experiment (NSTX) discharges. In the present study we extend this modeling to ELMing pulses in the Mega Ampere Spherical Tokamak (MAST) and MAST—Upgrade (MAST-U), where we find a clear resistive scaling for peeling–ballooning modes. While in NSTX ideal-MHD predicts stability for ELMing discharges, in MAST-U we find that the plasma is slightly unstable to peeling–ballooning modes, but is fully stabilized once diamagnetic effects are considered in terms of a growth rate normalization. A resistive power law scaling is calculated for these modes on MAST-U, which lies in between that of tearing modes and resistive interchange modes. A comparison between M3D-C1 and NIMROD shows reasonable agreement for this scaling. Resistivity destabilizes the modes and the peeling–ballooning unstable domain is considerably expanded in both, MAST and MAST-U. In addition to the MAST/-U pulses we also analyze resistive PB stability in a NSTX-similarity discharge on DIII-D. While having a different aspect ratio from NSTX, this discharge uses NSTX-like shaping parameters, toroidal field and plasma current. By considering these discharges alongside NSTX cases, we identify conditions influencing the onset of resistive peeling–ballooning modes. Furthermore, our findings indicate that magnetic shear in the pedestal region is closely linked to the emergence of resistive edge modes.

MAST-U↗

Effects beyond ideal MHD on stability of wide and enhanced pedestal regimes in NSTX

Stability of edge-localized modes (ELMs) in spherical tokamaks is explored using the extended MHD model. Linear NIMROD simulations have been performed for three NSTX discharges 132543, 132588, and 141133, to investigate the role of resistivity, diffusivity, and shear flows on the onset of ELMs. The first discharge represents the wide pedestal regime and the later two discharges represent ELM-free enhanced pedestal H-mode. We first present the effect of toroidal rotation shear and find a flow shear destabilizing effect in these NSTX discharges. Simulations are also extended to include the two-fluid and ion gyroviscosity effects. Simulations show that the flow shear can shift the mode spectrum and alter the critical condition of ELM onset. We also uncover that ELM onset prediction in spherical tokamaks requires effects beyond MHD, in particular gyroviscosity and diamagnetic terms could stabilize Peeling-Ballooning modes consistent with the experimental observation of ELM-free regimes in NSTX. The findings give new insight into the nature of the interplay between resistivity, flows, and diamagnetic stabilization in ELM suppression and have potential applications to ELM control schemes in NSTX-U and next-generation spherical tokamaks. This study identifies the essential physical effects that must be included in future predictive and validation simulations.

NSTX↗

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗

Magnetic phase diagram of the quantum spin chain compound SrCo 2 V 2 O 8 : a single-crystal neutron diffraction study

We explore magnetic order in the quantum spin chain compound SrCo 2 V 2 O 8 up to 14.9 T and down to 50 mK, using single-crystal neutron diffraction. Upon cooling in zero-field, commensurate antiferromagnetic (C-AFM) order with modulation vector ${{\boldsymbol{k}}}_{{\rm{C}}}$ = (0, 0, 1) develops below T N ≃ 5.0 K. Applying an external magnetic field (H∥c axis) destabilizes this C-AFM order, leading to an order-disorder transition between T N and ~1.5 K. Below 1.5 K, a commensurate to incommensurate (IC-AFM) transition occurs at 3.9 T, above which the magnetic reflections can be indexed by ${{\boldsymbol{k}}}_{{\rm{IC}}}$ = (0, 0, 1 ± δl). The incommensurability δl scales monotonically with H until the IC-AFM order disappears around 7.0 T. Magnetic reflections modulated by ${{\boldsymbol{k}}}_{{\rm{C}}}$ emerge again at higher fields. While the characters of the C-AFM, IC-AFM and the emergent AFM order in SrCo 2 V 2 O 8 appear to fit the descriptions of the Néel, longitudinal spin density wave and transverse AFM order observed in the related compound BaCo 2 V 2 O 8 , our results also reveal several unique signatures that are not present in the latter, highlighting the inadequacy of mean-field theory in addressing the complex magnetic order in systems of this class.

36 MATERIALS SCIENCE↗

Shape-driven, emergent behavior in active particle mixtures

Abstract Active particle systems can vary greatly from one-component systems of spheres to mixtures of particle shapes at different composition ratios. We investigate computationally the combined effect of anisotropy and stoichiometry on the collective behavior of two-dimensional active colloidal mixtures of polygons. We uncover three emergent phenomena not yet reported in active Brownian particle systems. First, we find that mixtures containing hexagons exhibit micro-phase separation with large grains of hexagonal symmetry. We quantify a measurable, implicit ‘steric attraction’ between the active particles as a result of shape anisotropy and activity. This calculation provides further evidence that implicit interactions in active systems, even without explicit attraction, can lead to an effective preferential attraction between particles. Next, we report stable fluid clusters in mixtures containing one triangle or square component. We attribute the fluidization of the dense cluster to the interplay of cluster destabilizing particles, which introduce grain boundaries and slip planes into the system, causing solid-like clusters to break up into fluid clusters. Third, we show that fluid clusters can coexist with solid clusters within a sparse gas of particles in a steady state of three coexisting phases. Our results highlight the potential for a wide variety of behavior to be accessible to active matter systems and establish a route to control active colloidal systems through simple parameter designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Low-scale supersymmetry from inflation

We investigate an inflation model with the inflaton being identified with a Higgs boson responsible for the breaking of U(1){sub B−L} symmetry. We show that supersymmetry must remain a good symmetry at scales one order of magnitude below the inflation scale, in order for the inflation model to solve the horizon and flatness problems, as well as to account for the observed density perturbation. The upper bound on the soft supersymmetry breaking mass lies between 1 TeV and 10{sup 3} TeV. Interestingly, our finding opens up a possibility that universes with the low-scale supersymmetry are realized by the inflationary selection. Our inflation model has rich implications; non-thermal leptogenesis naturally works, and the gravitino and moduli problems as well as the moduli destabilization problem can be solved or ameliorated; the standard-model higgs boson receives a sizable radiative correction if the supersymmertry breaking takes a value on the high side ∼ 10{sup 3} TeV.

79 ASTRONOMY AND ASTROPHYSICS↗

Bubble collision with gravitation

In this paper, we study vacuum bubble collisions with various potentials including gravitation, assuming spherical, planar, and hyperbolic symmetry. We use numerical calculations from double-null formalism. Spherical symmetry can mimic the formation of a black hole via multiple bubble collisions. Planar and especially hyperbolic symmetry describes two bubble collisions. We study both cases, when two true vacuum regions have the same field value or different field values, by varying tensions. For the latter case, we also test symmetric and asymmetric bubble collisions, and see details of causal structures. If the colliding energy is sufficient, then the vacuum can be destabilized, and it is also demonstrated. This double-null formalism can be a complementary approach in the context of bubble collisions.

79 ASTRONOMY AND ASTROPHYSICS↗

Stability analysis of inflation with an SU(2) gauge field

We study anisotropic cosmologies of a scalar field interacting with an SU(2) gauge field via a gauge-kinetic coupling. We analyze Bianchi class A models, which include Bianchi type I, II, VI{sub 0}, VII{sub 0}, VIII and IX. The linear stability of isotropic inflationary solution with background magnetic field is shown, which generalizes the known results for U(1) gauge fields. We also study anisotropic inflationary solutions, all of which turn out to be unstable. Then nonlinear stability for the isotropic inflationary solution is examined by numerically investigating the dependence of the late-time behaviour on the initial conditions. We present a number of novel features that may well affect physical predictions and viability of the models. First, in the absence of spatial curvature, strong initial anisotropy leads to a rapid oscillation of gauge field, thwarting convergence to the inflationary attractor. Secondly, the inclusion of spatial curvature destabilizes the oscillatory attractor and the global stability of the isotropic inflation with gauge field is restored. Finally, based on the numerical evidence combined with the knowledge of the eigenvalues for various inflationary solutions, we give a generic lower-bound for the duration of transient anisotropic inflation, which is inversely proportional to the slow-roll parameter.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining dark energy through the stability of cosmic structures

For a general dark-energy equation of state, we estimate the maximum possible radius of massive structures that are not destabilized by the acceleration of the cosmological expansion. A comparison with known stable structures constrains the equation of state. The robustness of the constraint can be enhanced through the accumulation of additional astrophysical data and a better understanding of the dynamics of bound cosmic structures.

79 ASTRONOMY AND ASTROPHYSICS↗

Stability analysis and future singularity of the m{sup 2} R □{sup -2} R model of non-local gravity

We analyse the classical stability of the model proposed by Maggiore and Mancarella, where gravity is modified by a term ∼ m{sup 2} R □{sup -2} R to produce the late-time acceleration of the expansion of the universe. Our study takes into account all excitations of the metric that can potentially drive an instability. There are some subtleties in identifying these modes, as a non-local field theory contains dynamical fields which yet do not correspond to degrees of freedom. Since some of them are ghost-like, we clarify the impact of such modes on the stability of the solutions of interest that are the flat space-time and cosmological solutions. We then find that flat space-time is unstable under scalar perturbations, but the instability manifests itself only at cosmological scales, i.e. out of the region of validity of this solution. It is therefore the stability of the FLRW solution which is relevant there, in which case the scalar perturbations are known to be well-behaved by numerical studies. By finding the analytic solution for the late-time behaviour of the scale factor, which leads to a big rip singularity, we argue that the linear perturbations are bounded in the future because of the domination of Hubble friction. In particular, this effect damps the scalar ghost perturbations which were responsible for destabilizing Minkowski space-time. Thus, the model remains phenomenologically viable.

79 ASTRONOMY AND ASTROPHYSICS↗

Fate of electroweak vacuum during preheating

Our electroweak vacuum may be metastable in light of the current experimental data of the Higgs/top quark mass. If this is really the case, high-scale inflation models require a stabilization mechanism of our vacuum during inflation. A possible candidate is the Higgs-inflaton/-curvature coupling because it induces an additional mass term to the Higgs during the slow roll regime. However, after inflation, the additional mass term oscillates, and it can destabilize our electroweak vacuum via production of large Higgs fluctuations during the inflaton oscillation era. In this paper, we study whether or not the Higgs-inflaton/-curvature coupling can save our vacuum by properly taking account of Higgs production during the preheating stage. We put upper bounds on the Higgs-inflaton and -curvature couplings, and discuss possible dynamics that might relax them.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Postinflationary vacuum instability and Higgs-inflaton couplings

The Higgs-inflaton coupling plays an important role in the Higgs field dynamics in the early Universe. Even a tiny coupling generated at loop level can have a dramatic effect on the fate of the electroweak vacuum. Such Higgs-inflaton interaction is present both at the trilinear and quartic levels in realistic reheating models. In this work, we examine the Higgs dynamics during the preheating epoch, focusing on the effects of the parametric and tachyonic resonances. We use lattice simulations and other numerical tools in our studies. We find that the resonances can induce large fluctuations of the Higgs field which destabilize the electroweak vacuum. Our considerations thus provide an upper bound on quartic and trilinear interactions between the Higgs and the inflaton. We conclude that there exists a favorable range of the couplings within which the Higgs field is stabilized during both inflation and preheating epochs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Decoherence can relax cosmic acceleration

In this work we investigate the semi-classical backreaction for a quantised conformal scalar field and classical vacuum energy. In contrast to the usual approximation of a closed system, our analysis includes an environmental sector such that a quantum-to-classical transition can take place. We show that when the system decoheres into a mixed state with particle number as the classical observable de Sitter space is destabilized, which is observable as a gradually decreasing Hubble rate. In particular we show that at late times this mechanism can drive the curvature of the Universe to zero and has an interpretation as the decay of the vacuum energy demonstrating that quantum effects can be relevant for the fate of the Universe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Astrophysical constraints on singlet scalars at LHC

We consider the viability of new heavy gauge singlet scalar particles at colliders such as the LHC . Our original motivation for this study came from the possibility of a new heavy particle of mass ∼ TeV decaying significantly into two photons at colliders, such as LHC, but our analysis applies more broadly. We show that there are significant constraints from astrophysics and cosmology on the simplest UV complete models that incorporate such new particles and its associated collider signal. The simplest and most obvious UV complete model that incorporates such signals is that it arises from a new singlet scalar (or pseudo-scalar) coupled to a new electrically charged and colored heavy fermion. Here we show that these new fermions (and anti-fermions) would be produced in the early universe, then form new color singlet heavy mesons with light quarks, obtain a non-negligible freeze-out abundance, and remain in kinetic equilibrium until decoupling. These heavy mesons possess interesting phenomenology, dependent on their charge, including forming new bound states with electrons and protons. We show that a significant number of these heavy states would survive for the age of the universe and an appreciable number would eventually be contained within the earth and solar system. We show that this leads to detectable consequences, including the production of highly energetic events from annihilations on earth, new spectral lines, and, spectacularly, the destabilization of stars. The lack of detection of these consequences rules out such simple UV completions, putting pressure on the viability of such new particles at LHC . To incorporate such a scalar would require either much more complicated UV completions or even further new physics that provides a decay channel for the associated fermion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling of lithium granule injection in NSTX using M3D-C1

In this paper, we present simulations of pedestal control by lithium granule injection (LGI) in NSTX. A model for small granule ablation has been implemented in the M3D-C1 code (Jardin et al 2012 Comput. Sci. Discovery 5 014002), allowing the simulation of realistic lithium granule injections. 2D and 3D simulations of Li injections in NSTX H-mode plasmas are performed and the effect of granule size, injection angle and velocity on the pedestal gradient increase is studied. The amplitude of the local pressure perturbation caused by the granules is found to be highly dependent on the solid granule size. Adjusting the granule injection velocity allows one to inject more particles at the pedestal top. 3D simulations show the destabilization of high order MHD modes whose amplitude is directly linked to the localized pressure perturbation, which is found to depend on the toroidal localization of the granule density source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of the n =1 resistive wall modes in the ITER high-mode confinement

The n = 1 resistive wall mode (RWM) stability of ITER high-mode confinement is investigated with bootstrap current included for equilibrium, together with the rotation and diamagnetic drift effects for stability. Here, n is the toroidal mode number. We use the CORSICA code for computing the free boundary equilibrium and AEGIS code for stability. We find that the inclusion of bootstrap current for equilibrium is critical. It can reduce the local magnetic shear in the pedestal, so that the infernal mode branches can develop. Consequently, the n = 1 modes become unstable without a stabilizing wall at a considerably lower beta limit, driven by the steep pressure gradient in the pedestal. Typical values of the wall position stabilize the ideal mode, but give rise to the 'pedestal' resistive wall modes. We find that the rotation can contribute a stabilizing effect on RWMs and the diamagnetic drift effects can further improve the stability in the co-current rotation case. But, generally speaking, the rotation stabilization effects are not as effective as the case without including the bootstrap current effects on equilibrium. We also find that the diamagnetic drift effects are actually destabilizing when there is a counter-current rotation.

Zheng, L. J.↗

Dependence of alpha-particle-driven Alfvén eigenmode linear stability on device magnetic field strength and consequences for next-generation tokamaks

Recently-proposed tokamak concepts use magnetic fields up to 12 T, far higher than in conventional devices, to reduce size and cost. Theoretical and computational study of trends in plasma behavior with increasing field strength is critical to such proposed devices. This paper considers trends in Alfvén eigenmode (AE) stability. Energetic particles, including alphas from deuterium–tritium fusion, can destabilize AEs, possibly causing loss of alpha heat and damage to the device. AEs are sensitive to device magnetic field via the field dependence of resonances, alpha particle beta, and alpha orbit width. We describe the origin and effect of these dependences analytically and by using recently-developed numerical techniques (Rodrigues et al 2015 Nucl. Fusion 55 083003). The work suggests high-field machines where fusion-born alphas are sub-Alfvénic or nearly sub-Alfvénic may partially cut off AE resonances, reducing growth rates of AEs and the energy of alphas interacting with them. High-field burning plasma regimes have non-negligible alpha particle beta and AE drive, but faster slowing down time, provided by high electron density, and higher field strength reduces this drive relative to low-field machines with similar power densities. The toroidal mode number of the most unstable modes will tend to be higher in high magnetic field devices. Finally, the work suggests that high magnetic field devices have unique, and potentially advantageous, AE instability properties at both low and high densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Prediction of Alfvén eigenmode energetic particle transport in ITER scenarios with a critical gradient model

A reduced 1D, local, critical-gradient model of energetic particle (EP) transport by Alfvén eigenmodes (AEs)—the TGLF-EP+Alpha model—is applied to a much-studied ITER base case and variations with lower plasma current and lower current penetration. The TGLF-EP+Alpha model is a greatly reduced and computationally inexpensive model of EP transport. Such a reduced critical gradient model, while inapplicable to transport driven by strongly nonlinear or non-local abrupt events, is a valuable tool for scoping studies needed in scenario optimization for ITER and beyond. Furthermore, it relies on the assumption of critical-gradient AE transport with the critical EP density gradient determined by linear AE stability calculations in the TGLF gyro-Landau fluid code automated with the parallel-processed TGLF-EP wrapper. EP transport is treated with simultaneous drive of AEs by fusion-born alpha particles and fast ions born from a 1 MeV neutral beam injection (NBI) heating. The effect of simultaneous drive creates about 50% increased particle transport in both EP channels. High magnetic safety factor q and low shear $\hat{s}$ are generally destabilizing to AEs, but low shear tends to be more important. A tailored q-profile, steady-staterelevant scenario can reduce AE-induced EP redistribution by better than 25% over the ITER base case despite having half the total current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗