Search NASA⌕ Search

SEARCH · Search NASA

Results for “kink interaction”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Nonlinear MHD modeling of n = 1 RMP -induced pedestal transport and mode coupling effects on ELM suppression in KSTAR

We report fully suppressing edge-localized modes (ELMs), e.g., with resonant magnetic perturbations (RMPs), is essential to reach and sustain high-performance steady-state H-mode plasmas because large ELMs can significantly reduce the lifetime of divertor components in future tokamak reactors. RMP-driven ELM suppression in KSTAR has been modeled by coupling the neoclassical transport code PENTRC to the nonlinear 3D MHD code JOREK. We have found that the radial transport from the combined effects of the kink-peeling, tearing response, and neoclassical toroidal viscosity can explain the pedestal degradation observed in experiments. In addition, it has been found that the RMP response can increase the inter-ELM heat flux on the lower outer divertor by redistributing the heat transport between the divertor plates. In addition to the degraded pedestal, ELM suppression is also attributable to the RMP-induced mode interactions. While the linear stability of peeling-ballooning mode (PBMs) improves owing to the degraded pedestal, the PBM and RMP interaction increases the spectral transfer between edge harmonics, preventing catastrophic growth and the crash of unstable modes. Here, it turns out that the magnetic islands near the pedestal top can play a vital role in mediating the mode interactions.

mode coupling↗
DOE OSTI2023

Characterization and dynamics of the peculiar stream Jhelum: A tentative role for the Sagittarius dwarf galaxy

Stellar streams are a promising tool to study the Milky Way’s dark matter subhalo population, as interactions with subhalos are expected to leave visible imprints in the streams in the form of substructure. However, there may be other causes of substructure. Here we studied the kinematics and the unusual morphology of the stellar stream Jhelum. Using a combination of ground-based photometry and Gaia EDR3 astrometry, we characterized the morphology of Jhelum. We combined this new data with radial velocities from the literature to perform orbit integrations of the stream in static Galactic potentials. We also carried out N-body simulations in the presence of the Sagittarius dwarf galaxy. The new data reveal a previously unreported tertiary component in the stream, as well as several gaps and a kink-like feature in its narrow component. We find that for a range of realistic Galactic potentials, no single orbit is able to reproduce Jhelum’s radial velocity data entirely. A generic property of the orbital solutions is that they share a similar orbital plane to Sagittarius and this leads to repeated encounters with the stream. Using N-body simulations that include a massive Sagittarius, we explored its effect on Jhelum, and we show that these encounters are able to qualitatively reproduce the narrow and broad components in Jhelum, as well as create a tertiary component in some cases. We also find evidence that such encounters can result in an apparent increase in the velocity dispersion of the stream by a factor up to four due to overlapping narrow and broad components. Our findings suggest that the Jhelum stream is even more complex than once thought; however, its morphology and kinematics can tentatively be explained via the interactions with Sagittarius. In this scenario, the formation of Jhelum’s narrow and broad components occurs naturally, yet some of the smaller gap-like features remain to be explained.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2002

Surface Processes of Faceted Growth

Atomic force microscopy and high precision optical interferometry allow to analyze the processes that are in the core of our present understanding of faceted crystal growth. Some of these processes will be reviewed. Recent experiments suggest that the Gibbs-Thomson Law (GTL) may not be applicable to the weakly fluctuating strongly polygonized steps, with low kink density. Propagation rate of a straight short step segment changes with the segment length much steeper than predicted by GTL. If the step fluctuations are still well developed, the typical transition length may be determined just by the radius of rounded comers of the polygonized step. If fluctuations do not have enough time to develop, the kinetic effects may be essential. Indeed, 'communications' between comers of a short straight step segment may occur and, thus, the GTL may be implemented only via kink exchange. E.g., dissolutions of this short segment at macroscopic equilibrium occurs as follows: The kinks are split from one comer as a result of fluctuations, and annihilate with kinks of the opposite sign generated by another comer. Similarly, an exchange via chain of kinks on a rough step is the mechanism to implement the GT Law. If there is a supersaturation with respect to a large crystal, the kinks not only diffuse along the step, but move back the generating comers. This 'wind' prevents communication between the comers and places upper limit supersaturation only below which GTL is applicable. Steps replace kinks in the 3D case, though, of course, another dimensionality brings about another physics. Another specific problem to think about is if fluctuations are always fast enough to provide enough kinks for linear dependence of step rate on supersaturation. Increasing azimuthal polygonization of growth hillocks is a sign of insufficient fluctuation rate. Morphological stability of vicinal faces is usually considered on the basis of mutually parallel steps. Azimuthal anisotropy and interlacing of step bunches inducing the well known bunch splitting is another challenge. Among other challenging problems are: existence of surface diffusion on the crystal/so1ution interface, step interaction in solutions, achievement of step flow modes from liquids, impurity effects. New phenomena essential in biomacromolecular crystallization is another interesting area.

Chernov, A.↗
DOE OSTI2022

Kinematic Evidence for an Embedded Planet in the IM Lupi Disk

Abstract We test the hypothesis that an embedded giant planet in the IM Lupi protostellar disk can produce velocity kinks seen in CO line observations as well as the spiral arms seen in scattered light and continuum emission. We inject planets into 3D hydrodynamics simulations of IM Lupi, generating synthetic observations using Monte Carlo radiative transfer. We find that an embedded planet of 2–3 M Jup can reproduce non-Keplerian velocity perturbations, or “kinks”, in the 12 CO J = 2–1 channel maps. Such a planet can also explain the spiral arms seen in 1.25 mm dust continuum emission and 1.6 μ m scattered-light images. We show that the wake of the planet can be traced in the observed peak velocity map, which appears to closely follow the morphology expected from our simulations and from analytic models of planet–disk interaction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Discovery of multimechanisms of screw dislocation interaction in bcc iron from open-ended saddle point searches

Dislocation motion and interactions determine mechanical properties in body-centered cubic (bcc) metallic materials. However, studying mechanisms for the screw dislocation interaction is fundamentally challenging since many underlying processes involve mesotimescales and atomistic resolution, currently inaccessible by either experimental techniques or continuum theoretical methods. Here, in this paper, we develop a computational capability based on self-evolving atomistic kinetic Monte Carlo (SEAKMC) to sample the critical events and saddle point energies related to screw dislocations and their junctions. The method is first validated by calculating the stress dependence of Peierls barriers and formation energies of kink pairs and cross-slip kink pairs on a single screw dislocation in bcc iron. Then the method is applied to a binary junction of a pair of intersecting screw dislocations, the structure of which is crucial for low-temperature plastic deformation. We identify three important mechanisms: coplanar cross-slipping, jog-pinning, and a previously unknown unzipping mechanism during the evolution of the binary junction. The mechanisms are then further validated using classical molecular dynamics simulations. The computational capability developed in this paper provides an effective tool to evaluate screw dislocation related thermally activated events in complex stress conditions. The mechanisms discovered in this paper provide critical insights into temperature dependence of the anomalous slip, a breakdown of the Schmidt law, during the plastic deformation in bcc iron and can be generalized to other bcc metals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Non-linear MHD modelling of transients in tokamaks: a review of recent advances with the JOREK code

Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Global, quasi-steady dynamics of the distant solar wind 2: Deformation of the heliospheric current sheet

A three-dimensional (3-D) magnetohydrodynamic (MHD) numerical model is used to trace global deformations of the heliospheric current sheet (HCS) caused by large-scale dynamical interactions associated with corotating solar wind flows. Configurations incorporating the tilted-dipole geometry are investigated out past 30 AU for a variety of dipole tilt angles, alpha. Inclusion of the full, 3-D interplanetary dynamics allows north-south displacements and the east-west warping of the HCS by advective corotational effects to be accurately assessed for the first time. It is found that large-scale spatial correlations between velocity and density imposed at the coronal source (i.e., the geometric arrangement whereby regions of high-speed low-density material lying adjacent to areas of slow, dense flow interact obliquely under the influence of solar rotation) result in a distinctive pattern of deformation of the HCS. For an alpha = 30 deg tilted-dipole example, it is shown that typical zonal variations in radial velocity lead to significant folding of the HCS within about 5 AU of the Sun. By 10 AU, additional sharp bends appear near the latitudinal extremes of the HCS surface, where it is overtaken by shock fronts driven by 3-D corotating interaction regions (CIRs). Moreover, the model suggests that inside about 20 AU, major plasma structures are systematically organized, about the HCS, such that the greatest concentrations of material and magnetic field (the centroids of the 3-D CIR structures) are coincident with the folded crests of the HCS (near heliographic latitudes lambda = +/- alpha). Thus, in these circumstances many of the more interesting dynamical features inconveniently lie well away from the heliographic equator. At larger heliocentric distances, where neighboring CIRs begin to interact strongly, the warping of the HCS abates dramatically and the association between folds in the HCS and major field and density concentrations is weakened and ultimately breaks down. On the basis of the model calculations, quasi-steady deformations and related phenomena should figure prominently in the interplanetary medium, except during those periods when the magnetic dipole and solar spin axes are nearly aligned (alpha is approximately less than 10 deg). Even then, they may be of significance if there are any substantial local warps or kinks in the streamer belt.

Pizzo, V. J.↗
DOE OSTI2025

Geometric entropies and their Hamiltonian flows

In holographic theories, the Hubeny-Rangamani-Takayanagi (HRT) area operator plays a key role in our understanding of the emergence of semiclassical Einstein-Hilbert gravity. When higher derivative corrections are included, the role of the area is instead played by a more general functional known as the geometric entropy. It is thus of interest to understand the flow generated by the geometric entropy on the classical phase space. In particular, the fact that the associated flow in Einstein-Hilbert or Jackiw-Teitelboim (JT) gravity induces a relative boost between the left and right entanglement wedges is deeply related to the fact that gravitational dressing promotes the von Neumann algebra of local fields in each wedge to type II. This relative boost is known as a boundary-condition-preserving (BCP) kink-transformation. In a general theory of gravity (with arbitrary higher-derivative terms), it is straightforward to show that the flow continues to take the above geometric form when acting on a spacetime where the HRT surface is the bifurcation surface of a Killing horizon. However, the form of the flow on other spacetimes is less clear. In this paper, we use the manifestly-covariant Peierls bracket to explore such flows in two-dimensional theories of JT gravity coupled to matter fields with higher derivative interactions. The results no longer take a purely geometric form and, instead, demonstrate new features that should be expected of such flows in general higher derivative theories. We also show how to obtain the above flows using Poisson brackets.

AdS-CFT correspondence↗
NASA NTRS1987

Simulation studies for surfaces and materials strength

Three different investigations based on computer simulations were carried out. Simulation calculations were performed using model potentials with two- and three-body interactions which were represented by the Mie and the Axilrod-Teller potentials, respectively. In the first part, to analyze stability criteria for simple crystals, a parametrical investigation was performed and stability regions for five different crystalline structures were determined with respect to potential energy parameters. To simulate some of the basic surface processes such as diffusion and nucleation, in the second and third parts, silicon surfaces were considered as a specific example. In these studies diffusion of adatoms on stepped surfaces were simulated and energetics associated with the formation of kink sites were calculated.

Halicioglu, Timur↗
DOE OSTI2020

Anisotropic Motions of Fibrils Dictated by Their Orientations in the Lamella: A Coarse-Grained Model of a Plant Cell Wall

Plant cell walls are complex systems that exhibit the characteristics of both rigid and soft material depending on their external perturbations. The three main polymeric components in a plant primary cell wall are cellulose fibrils, hemicellulose, and pectins. These components interact in a hierarchical fashion giving rise to mesoscale structural features such as cellulose bundles, lamella stacking, and so on. Although several studies have focused on understanding these unique structural features, a clear picture linking them to cell wall mechanics is still lacking. As a first step toward this goal, a phenomenological model of plant cell wall has been developed in this work by using available experimental data to investigate the underlying connections between mesoscale structural features and the motions of fibrils during deformation. In this model cellulose fibrils exhibit motions such as angular reorientations and kinking upon forced stretching. Furthermore, these motions are dependent on the orientation of fibrils with respect to the stretch direction, i.e., fibrils that are at an angle to the stretch direction exhibit predominant angular reorientations, while fibrils transverse to the stretch direction undergo kinking as a result of transverse compression. Varying the chain length of pectin had negligible effects on these motions. One of the main contributions from this work is the development of a simple model that can be easily fine-tuned to test other hypotheses and extended to include additional experimental knowledge about the structural aspects of cell walls in the future.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Diffusion on Cu surfaces

Understanding surface diffusion is essential in understanding surface phenomena, such as crystal growth, thin film growth, corrosion, physisorption, and chemisorption. Because of its importance, various experimental and theoretical efforts have been directed to understand this phenomena. The Field Ion Microscope (FIM) has been the major experimental tool for studying surface diffusion. FIM have been employed by various research groups to study surface diffusion of adatoms. Because of limitations of the FIM, such studies are only limited to a few surfaces: nickel, platinum, aluminum, iridium, tungsten, and rhodium. From the theoretical standpoint, various atomistic simulations are performed to study surface diffusion. In most of these calculations the Embedded Atom Method (EAM) along with the molecular static (MS) simulation are utilized. The EAM is a semi-empirical approach for modeling the interatomic interactions. The MS simulation is a technique for minimizing the total energy of a system of particles with respect to the positions of its particles. One of the objectives of this work is to develop the EAM functions for Cu and use them in conjunction with the molecular static (MS) simulation to study diffusion of a Cu atom on a perfect as well as stepped Cu(100) surfaces. This will provide a test of the validity of the EAM functions on Cu(100) surface and near the stepped environments. In particular, we construct a terrace-ledge-kink (TLK) model and calculate the migration energies of an atom on a terrace, near a ledge site, near a kink site, and going over a descending step. We have also calculated formation energies of an atom on the bare surface, a vacancy in the surface, a stepped surface, and a stepped-kink surface. Our results are compared with the available experimental and theoretical results.

Karimi, Majid↗
DOE OSTI2020

Area-selective atomic layer deposition enabled by competitive adsorption

In this work, we investigate the atomic layer deposition (ALD) of ZrO2 thin films on Cu and SiO2 substrates, using Zr[N(C2H5CH3)]4 as the thin-film precursor, and H2O or O2 as the coreactants. Here, we introduce 3-hexyne as a coadsorbate molecule during the thin-film precursor half-cycle and examine its effect on the selectivity of growth. We find that 3-hexyne strongly inhibits growth on Cu, while having essentially no effect on the growth on SiO2. Calculations using dispersion-inclusive density functional theory verify that 3-hexyne undergoes sp → sp2 rehybridization on Cu, which results in strong chemisorption on the metal surface, while only binding weakly to SiO2 via nonbonded van der Waals/dispersion interactions. After 10 cycles of ALD using 3-hexyne as the coadsorbate, we observed the deposition of ∼1.5 nm of ZrO2 on SiO2. On a Cu substrate, we only detected <0.15 nm of ZrO2 after the same number of cycles of ALD. At this point in the process, we find evidence of the formation of cuprous oxide (Cu2O) from in situ x-ray photoelectron spectroscopy and a significant increase in the roughness of the Cu substrate. We conclude that both factors likely contribute to the loss of selectivity due to the formation of sites (e.g., Cu2O) that bind 3-hexyne less strongly and/or an increase in the density of highly reactive sites (e.g., steps, kinks) that promote dissociative chemisorption of the thin-film precursor.

Suh, Taewon (ORCID:0000000163823811)↗
DOE OSTI2021

Interfacial electron-phonon coupling and quantum confinement in ultrathin Yb films on graphite

Interfacial electron-phonon coupling in ultrathin films has attracted much interest recently. Here, by combining angle-resolved photoemission spectroscopy and scanning tunneling microscopy, we report quantized electronic states and strong interfacial electron-phonon coupling in ultrathin Yb films on graphite. We observed clear kinks in the energy-momentum dispersion of quantum well states, and the kink positions agree well with the energies of optical phonons of graphite. The extracted coupling strength λ is largest for the thinnest film with a preferred (“magic”) thickness of four monolayers and exhibits a strong band dependence, which can be qualitatively accounted for by a simple model. The interfacial electron-phonon coupling also gives rise to characteristic steplike structures in the $dI/dV$ spectra, implying dominant coupling with the phonons with zero in-plane momentum. A Lifshitz transition occurs at higher coverage, where quantum well states derived mainly from 5d electrons dominate near the Fermi level and possess large effective mass (up to ~ 19 m e ). Here our results highlight the potentially important role of interfacial electron-phonon interaction for ultrathin films and provide spectroscopic insight to understand this cross-interface fermion-boson interaction.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Three interaction energy scales in the single-layer high-$T_c$ cuprate HgBa 2 CuO 4+δ

The lamellar cuprate superconductors exhibit the highest ambient-pressure superconducting transition temperatures ($T_c$), and after more than three decades of extraordinary research activity, continue to pose formidable scientific challenges. A major experimental obstacle has been to distinguish universal phenomena from materials- or technique-dependent ones. Angle-resolved photoemission spectroscopy (ARPES) measures momentum-dependent single-particle electronic excitations and has been invaluable in the endeavor to determine the anisotropic momentum-space properties of the cuprates. HgBa 2 CuO 4+δ (Hg1201) is a single-CuO 2 -layer cuprate with a particularly high optimal Tc and a simple crystal structure, yet there exists little information from ARPES about the electronic properties of this model system. In this paper we present an ARPES study of doping-, temperature-, and momentum-dependent systematics of near-nodal dispersion anomalies in Hg1201. The data reveal a hierarchy of three distinct energy scales: a subgap low-energy kink, an intermediate-energy kink near 55 meV, and a peak-dip-hump structure. The first two features are attributed to the coupling of electrons to Ba-derived optical phonons and in-plane bond-stretching phonons, respectively. The nodal peak-dip-hump structure appears to have a common doping dependence in several single-layer cuprates and is interpreted as a manifestation of pseudogap physics at the node. These results establish several universal phenomena, both in terms of connecting multiple experimental techniques for a single material and in terms of connecting comparable spectral features in multiple structurally similar cuprates.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Enhanced Pedestal H-mode at low edge ion collisionality on NSTX

The Enhanced Pedestal (EP) H-mode regime is an attractive wide-pedestal ELM-free high-betap scenario for NSTX-U and next-step devices as it achieves enhanced energy confinement (H98y,2 > 1.5), large normalized pressure (betaN > 5) and significant bootstrap fraction (f_BS > 0.6) at I_p/B_T = 2 MA/T. This regime is realized when the edge ion collisionality becomes sufficiently small that a positive feedback interaction occurs between a reduction in the ion neoclassical energy transport and an increase in the particle transport from pressure-driven edge instabilities. EP H-mode was most often observed as a transition following a large ELM in conditions with low edge neutral recycling. It is hypothesized that the onset of pressure-driven instabilities prior to the full recovery of the neutral density leads to a temporary period with elevated ion temperature gradient that triggers the transition to EP H-mode. Linear CGYRO and M3D-C1 calculations are compared to beam emission spectroscopy (BES) and magnetic spectroscopy in order to describe the evolution of the edge particle transport mechanisms during the ELM recovery and the saturated EP H-mode state. The observations are consistent with the hypothesis that the onset of pressure-driven edge instabilities, such as the KBM and kink-peeling, can be responsible for the increased particle transport in EP H-mode.

edge and boundary physics↗
DOE OSTI2020

Enhanced pedestal H-mode at low edge ion collisionality on NSTX

The Enhanced Pedestal (EP) H-mode regime is an attractive wide-pedestal high-βp scenario for the National Spherical Torus Experiment Upgrade (NSTX-U) and next-step devices as it achieves enhanced energy confinement (H98y,2 > 1.5), large normalized pressure (βN > 5), and significant bootstrap fraction (fBS > 0.6) at Ip/BT = 2 MA/T. This regime is realized when the edge ion collisionality becomes sufficiently small that a positive feedback interaction occurs between a reduction in the ion neoclassical energy transport and an increase in the particle transport from pressure-driven edge instabilities. The EP H-mode was most often observed as a transition following a large edge-localized mode in conditions with low edge neutral recycling. It is hypothesized that the onset of pressure-driven instabilities prior to the full recovery of the neutral density leads to a temporary period with an elevated ion temperature gradient that triggers the transition to EP H-mode. Linear CGYRO and M3D-C1 calculations are compared to beam emission spectroscopy and magnetic spectroscopy in order to describe the evolution of the edge particle transport mechanisms during the ELM recovery and the saturated EP H-mode state. The observations are consistent with the hypothesis that the onset of pressure-driven edge instabilities, such as the kinetic ballooning mode and kink-peeling mode, can be responsible for the increased particle transport in EP H-mode.

Battaglia, D. J. (ORCID:0000000188979740)↗
DOE OSTI2023

Anomalous Hall effect in conical helimagnetic crystals

Spin-spiral texture can substantially change charge transport properties in helimagnets. Here we find the anomalous Hall effect (AHE) exhibiting the dramatic behavior with respect to chemical potential μ in conical magnetic structures. The direct conductivity demonstrates kinks, and the anomalous Hall current exhibits minima and maxima changing the sign. We analytically derive the expression for energy bands and eigenstates in the most general case. Because of the conical potential, the energy bands are split into two nonparabolic bands where the lower band can have one- or two-minima shapes in the k z direction ($\hat{z}$ is a direction of the spiral axis). We prove that the origin of the anomalous Hall effect is not topological and is due to the interplay between the asymmetry of energy bands in the x and z directions and spin restrictions in the phase space due to the conical potential. We also investigate the dependence of transport properties on cone half- angle θ, and find that the effects are most pronounced at θ=π/2 (a helical state). Electric current is calculated using the Boltzmann equation where the relaxation is caused by electron-acoustic phonon interaction. The transition probability is found to be a 2×2 matrix with nonvanishing off-diagonal elements indicating the strong interband transitions. The origin of interband transitions is because of the nature of the conical potential where conduction electron spins interact with localized magnetic moments. To verify the proposed theory, we calculate the temperature dependence of resistivity for MnSi crystals and find the discontinuity at the phase transition between conical and paramagnetic phases. The calculations are in the excellent agreement with the experimental data. In addition, we predict the discontinuity behavior for the anomalous Hall resistivity at the phase transition where the resistivity exhibits the abrupt change at T=T C , (a) to zero if the relativistic effects for the conduction electrons are small or (b) to a nonzero value if Rashba/Dresselhaus effects are taken into account.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗