Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D MHD stability”

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 OSTI2019

Non-axisymmetric equilibrium and stability using the ELITE stability code

Linear perturbation theory is used to model the ideal magnetohydrodynamic stability of tokamak equilibria under the application of external 3D magnetic perturbations. We use the ELITE code to produce both a linear plasma response, as well as the linear axisymmetric toroidal eigenmodes which are used as basis functions for the 3D perturbative stability analysis. In this work, the symmetry breaking produces coupling of modes with different toroidal mode number, n, leading to modification of the linear growth rates, and poloidal localisation of the nonaxisymmetric peeling–ballooning mode at sufficiently large applied field.

3D MHD stability↗
DOE OSTI2025

Stability evaluation and mitigation strategies in advanced tokamaks using 3D MHD spectroscopy

Multi-modal, active 3D MHD spectroscopy is applied in high-performance advanced tokamak scenarios to study their stability time evolution, revealing an intriguing dependence on both and . A tailored applied 3D field provides a 3D plasma response to extract the growth rate of the least stable mode. The estimated growth rate finds a decrease in stability when the minimum in the safety factor (q) passes through 2.0 and reveals inherent risks of crossing an additional rational surface at integer q min , even above the usual q = 1 sawtooth condition. Based on this result, the potential scenario in which q min ~ 2 can be safely crossed during a more stable lower β N phase was investigated, and the improved stability of this scenario is confirmed by the estimated growth rate. This shows that 3D MHD spectroscopy can offer insights into strategies for improving stability by identifying the vulnerable aspects of such scenarios. In addition, the method highlights its potential for instability avoidance by enabling early detection of multiple modes, even before magnetic coils can measure them. The measured growth rate by the 3D MHD spectroscopy shows its reliability by exhibiting a correlation with the programmed rises in plasma beta across various high β N and high q min discharges. In addition, this method is successfully applied during rapidly evolving I p ramp-up phases, a key part of the scenario development. By achieving reasonable growth rate measurements at high-performance scenario developments, this technique contributes to the development of advanced diagnostic tools for tokamak scenario stability, which will help identify an effective pathway to stable, high-performance scenarios.

Yang, S. M. [Princeton Plasma Physics Laboratory (↗
NASA NTRS2004

3D MHD Models of Active Region Loops

Present imaging and spectroscopic observations of active region loops allow to determine many physical parameters of the coronal loops, such as the density, temperature, velocity of flows in loops, and the magnetic field. However, due to projection effects many of these parameters remain ambiguous. Three dimensional imaging in EUV by the STEREO spacecraft will help to resolve the projection ambiguities, and the observations could be used to setup 3D MHD models of active region loops to study the dynamics and stability of active regions. Here the results of 3D MHD models of active region loops are presented, and the progress towards more realistic 3D MHD models of active regions. In particular the effects of impulsive events on the excitation of active region loop oscillations, and the generation, propagations and reflection of EIT waves are shown. It is shown how 3D MHD models together with 3D EUV observations can be used as a diagnostic tool for active region loop physical parameters, and to advance the science of the sources of solar coronal activity.

Ofman, Leon↗
DOE OSTI2024

MHD stability of spherical tokamak equilibria with non-monotonic q-profiles

We use the 3D magnetohydrodynamic (MHD) code M3D-C1 [Jardin et al., Comput. Sci. Discovery 5, 014002 (2012)] to examine the MHD stability and subsequent evolution of NSTX shot 129169. This discharge had a period with a non-monotonic safety factor profile, q (reversed shear), which was terminated by a MHD event that abruptly lowered the central safety factor, q0, and greatly reduced the peakedness of the pressure profile. We show that the equilibrium just before the MHD event occurred was linearly unstable to many pressure-driven infernal modes. Modes with toroidal mode number n≥3 all had rational surfaces very close to the minimum value of q. However, a non-resonant pressure-driven (1, 1) mode was also present, and this dominated the nonlinear evolution. The final state in the simulation, after the MHD activity subsided, had a reduced and flattened pressure profile and a nearly monotonic q-profile, in qualitative agreement with experimental results. The initial state was also unstable to the resistive interchange criteria in the reversed-shear region, but the final state was stable everywhere. The “double tearing mode” (DTM) does not appear to play a role in the MHD activity of this discharge. In Appendix A, we show that in a torus, the DTM is strongly stabilized by pressure, but it is destabilized in cylindrical geometry (which has been the most extensively analyzed in the literature).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2006

Relativistic MHD Simulations of Precessed Jets

Relativistic jets are considered to be generated by magnetic fields in a rotating black hole with accretion disk. Consequently, resulting outflows contain magnetic fields in them and control the propagation of jets. We have performed 3D relativistic MHD simulations to investigate the stability and structure of precessed MHD jets with large Lorentz factor by using a newly developed 3D GRMHD code. We have performed simulations of supermagnetosonic jets surrounded by a fast wind. The simulation results reveal complex pressure structure inside the RMHD jet. The structure is produced by a combination of the helical surface and body modes excited by the precession as predicted theoretically. The wavelength of the body mode which occurs in an internal helical twist is much shorter than that of the helical twist surface mode. We will present some comparisons between the RMHD simulations and theoretical predictions, and potential observables and discuss the effect of wind.

Mizuno, Y.↗
DOE OSTI2019

Resistive stability of cylindrical MHD equilibria with radially localized pressure gradients

As a step toward understanding 3D magnetohydrodynamic (MHD) equilibria, for which smooth solutions may not exist, we develop a simple cylindrical model to investigate the resistive stability of MHD equilibria with alternating regions of constant and nonuniform pressure, producing states with continuous total pressure (i.e., no singular current sheets) but discontinuities in the parallel current density. We examine how the resistive stability characteristics of the model change as we increase the localization of pressure gradients at fixed radii, which approaches a discontinuous pressure profile in the zero-width limit. Equilibria with continuous pressure are found to be unstable to moderate/high-m modes and apparently tend toward ideal instability in some cases. We propose that additional geometric degrees of freedom or symmetry breaking via island formation may increase the parameter space on which equilibria of our model are physically realizable, while preserving the radial localization of pressure gradients. This is consistent with the possibility of realizing, in practice, 3D MHD equilibria which support both continuously nested flux surfaces (where ∇p ≠ 0) and chaotic field regions (where ∇p = 0)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Identification of multiple eigenmode growth rates towards real time detection in DIII-D and KSTAR tokamak plasmas

Abstract The successful application of three-dimensional (3D) magnetohydrodynamic (MHD) spectroscopy enables us to identify the multi-mode eigenvalues in DIII-D and KSTAR tokamak experiments with stable plasmas. The temporal evolution of the multi-modes’ stabilities have been detected. The new method is numerically efficient allowing the real time detection of MHD modes’ stabilities during the discharge. The method performs active detection of the plasma stability by utilizing the upper and lower rows of internal non-axisymmetric coils to apply a wide variety of 3D fields. Multi-mode eigenvalues are extracted using subspace system identification of the plasma response measured by 3D-field magnetic sensors distributed at different poloidal locations. The equivalence of this new method with the one introduced by Wang (2019 Nucl. Fusion 59 024001) has been numerically corroborated. The more robust and efficient calculation developed here will enable real time monitoring of the plasma stability based on the extracted eigenvalues of stable modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ideal MHD induced temperature flattening in spherical tokamaks

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Study of the Alfven eigenmodes stability in CFQS plasma using a Landau closure model

The aim of this study is to analyze the stability of the Alfven eigenmodes (AE) in the Chinese First Quasi-axisymmetric Stellarator (CFQS). The AE stability is calculated using the code FAR3d that solves the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moment for the energetic particles (EP) species including the effect of the helical couplings and acoustic modes. The Landau damping and resonant destabilization effects are added in the model by a given closure relation. The simulation results indicate the destabilization of n = 1 to 4 AEs by EP during the slowing down process, particularly n = 1 and n = 2 toroidal AEs (TAE), n = 3 elliptical AE (EAE) and n = 4 non circular AE (NAE). If the resonance is caused by EPs with an energy above 17 keV (weakly thermalized EP), n = 2 EAEs and n = 3 NAEs are unstable. On the other hand, EPs with an energy below 17 keV (late thermalization stage) lead to the destabilization of n = 3 and n = 4 TAEs. The simulations for an off-axis NBI injection indicate the further destabilization of n = 2 to 4 AEs although the growth rate of the n = 1 AEs slightly decreases, so no clear optimization trend with respect to the NBI deposition region is identified. In addition, n = 2, 4 helical AE (HAE) are unstable above an EP β threshold. Also, if the thermal β of the simulation increases (higher thermal plasma density) the AE stability of the plasma improves. Overall, the simulations including the effect of the finite Larmor radius and electron-ion Landau damping show the stabilization of the n = 1 to 4 EAE/NAEs as well as a decrease of the growth rate and frequency of the n = 1 to 4 BAE/TAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Inputs and Outputs for paper "The DESC stellarator code suite. Part 1. Quick and accurate equilibria computations"

Contains the DESC and VMEC input and output files used in the paper, as well as the plotting scripts used to create the figures in the paper (Current with the arXiv Mar 31 2022 version, https://arxiv.org/abs/2203.17173v1): 3D equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two 3D equilibrium codes, VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the MHD force error in real space directly. Accuracy as measured by final plasma energy and satisfaction of MHD force balance, as well as other metrics, will be presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near-axis. DESC's global Fourier-Zernike basis also yields the solution everywhere in the plasma volume, not just on discrete flux surfaces. Further, DESC can compute the same accuracy solution as VMEC in an order of magnitude less time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

FES 2023 Theory & Simulation Performance Target

The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Modeling of the ECCD injection effect on the Heliotron J and LHD plasma stability

The aim of the study is to analyze the stability of the energetic particle modes (EPM) and Alfven Eigenmodes (AE) in Helitron J and LHD plasma if the electron cyclotron current drive (ECCD) is applied. Additionally, the analysis is performed using the code FAR3d that solves the reduced MHD equations describing the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particle (EP) species, including the effect of the acoustic modes. The Landau damping and resonant destabilization effects are added via the closure relation. The simulation results show that the n = 1 EPM and n = 2 global AE (GAE) in Heliotron J plasma can be stabilized if the magnetic shear is enhanced at the plasma periphery by an increase (co-ECCD injection) or decrease (ctr-ECCD injection) of the rotational transform at the magnetic axis ($\rlap{-} \iota_{0}$). In the ctr-ECCD simulations, the EPM/AE growth rate decreases only below a given $\rlap{-} \iota_{0}$, similar to the ECCD intensity threshold observed in the experiments. In addition, ctr-ECCD simulations show an enhancement of the continuum damping. The simulations of the LHD discharges with ctr-ECCD injection indicate the stabilization of the n = 1 EPM, n = 2 toroidal AE (TAE) and n = 3 TAE, caused by an enhancement of the continuum damping in the inner plasma leading to a higher EP β threshold with respect to the co- and no-ECCD simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Analysis of Alfven eigenmode destabilization in ITER using a Landau closure model

Alfvén eigenmodes (AE) can be destabilized during ITER discharges driven by neutral beam injection (NBI) energetic particles (EP) and alpha particles. The aim of the present study is to analyze the AE stability of different ITER operation scenarios considering multiple energetic particle species. We use the reduced magneto-hydrodynamic (MHD) equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the EP species including the effect of the acoustic modes. Furthermore, the AEs driven by the NBI EP and alpha particles are stable in the configurations analyzed, only MHD-like modes with large toroidal couplings are unstable, although both can be destabilized if the EP increases above a threshold. The threshold is two times the model value for the NBI EP and alpha particles in the reverse shear (RS) case, leading to the destabilization of Beta induced AE (BAE) near the magnetic axis with a frequency of kHz and toroidal or elliptical AE (TAE/EAE) in the RS region with a frequency of kHz, respectively. On the other hand, the hybrid and steady state configurations show a threshold 3 times larger with respect to the model for the alpha particle and 40 times for the NBI EP, also destabilizing BAE and TAE between the inner and middle plasma region. In addition, a extended analysis of the RS scenario where the of both alpha particles and NBI EP are above the AE threshold, multiple EP damping effects are also identified as well as optimization trends regarding the resonance properties of the alpha particle and NBI EP with the bulk plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Pressure-driven tearing and thermal transport in finite-beta reversed field pinch computations

In this work, nonlinear resistive-magnetohydrodynamics (MHD) computation with heating and anisotropic transport is applied to examine the interaction between thermal energy and magnetic fluctuations in inductively driven reversed-field pinches (RFPs). The magnetic fluctuations underlie magnetic field reversal through dynamo-like correlations, and they enhance thermal energy transport through fluctuations of parallel heat flux density. With the unfavorable magnetic curvature that exists across the RFP profile, thermal energy also affects the magnetic fluctuations. Computations with the NIMROD code [Sovinec et al., J. Comput. Phys. 195, 355–386 (2004)] integrate nonlinear MHD dynamics with energy transport and reproduce an RFP state with experimentally relevant values of plasma-β. Equilibria constructed from results of the 3D computations are analyzed to assess the sources of free energy in the saturated nonlinear state. Linear computations for these profiles show unstable modes of tearing parity. Their eigenfunctions are used to evaluate and compare stabilizing and destabilizing contributions to the kinetic energy integral. An assessment of the drives in the integral reveals that the pressure gradient drive is of comparable magnitude to the parallel current drive, and only the sum of the two surpasses the stabilizing contributions. Correlation of magnetic and parallel heat flux density fluctuations in the nonlinear computations shows that fluctuation-induced thermal conduction is the dominant mode of energy loss, as expected from experimental evidence. Decomposition of the fluctuating heat flux density shows that second-order correlations, alone, do not explain the total energy transport. Higher-order correlations are also important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Analysis of fusion fuel conditions and instability development in shocklessly compressed MagLIF implosion simulations

Magnetized liner inertial fusion (MagLIF) implosions on the Z accelerator have almost exclusively been driven by ∼100-ns rise time current pulses. The rise time is selected to be as short as achievable on Z partially to minimize the time during which deleterious implosion instabilities can develop. Modifying the shape of the current pulse could provide benefits for MagLIF, including more efficient compression of the fusion fuel and the magnetic flux inside the liner cavity. Quasi-isentropic compression of the liner prevents formation of shocks in the liner material and reduces the amount of entropy generation within the liner. This allows for more final compression of the liner and fuel assembly. We present results from one-dimensional (1D) radiation-magnetohydrodynamic (rad-MHD) simulations comparing thermonuclear fuel conditions in MagLIF implosions driven with two different current pulses: a ∼100-ns rise time, ∼21.5 MA peak current “short pulse” and a ∼200-ns rise time, ∼21.5 MA peak current “shockless” pulse. We also quantify and compare the instability development in three-dimensional (3D) MHD implosion simulations driven by these two different pulse shapes. Our 1D simulations indicate that the shocklessly compressed MagLIF implosion performs better than the short pulse driven implosion with a >50% higher thermonuclear neutron yield, and 3D simulations indicate comparable implosion instability development, suggesting that pulse shaping could enable improvements to MagLIF performance on Z without compromising implosion stability.

Shipley, G. A. (ORCID:0000000205364001)↗
NASA NTRS2007

3D Relativistic MHD Simulations of Magnetized Spine-Sheath Relativistic Jets

We have performed numerical simulations of weakly and strongly magnetized relativistic jets embedded in a weakly and strongly magnetized stationary or mildly relativistic (0.5c) sheath using the RAISHIN code. In the numerical simulations a jet with Lorentz factor gamma=2.5 is precessed to break the initial equilibrium configuration. Results of the numerical simulations are compared to theoretical predictions from a normal mode-analysis of the linearized RMHD equations describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized moving sheath. The prediction of increased stability of a weakly-magnetized system with mildly relativistic sheath flow to Kelvin-Helmholtz instabilities and the stabilization of a strongly-magnetized system with mildly relativistic sheath flow is confirmed by the numerical simulations.

Mizuno, Yosuke↗