Search NASA⌕ Search

SEARCH · Search NASA

Results for “MHD (MAGNETOHYDRODYNAMICS)”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Hints of primordial magnetic fields at recombination and implications for the Hubble tension

Primordial Magnetic Fields (PMFs), long studied as relics of the early Universe, accelerate recombination and have been proposed as a way to relieve the Hubble tension. However, previous studies relied on simplified toy models. Here, in this study, we use recent evaluations of recombination with PMFs, incorporating full magnetohydrodynamic (MHD) simulations and detailed Lyman-alpha radiative transfer, to test PMF-enhanced recombination (bΛCDM) against observational data from the cosmic microwave background (CMB), baryon acoustic oscillations (BAO), and Type Ia supernovae (SN). Focusing on non-helical PMFs with a Batchelor spectrum, we find a preference for presentday total field strengths of approximately 5-10 pico-Gauss. Depending on the dataset combination, this preference ranges from mild (∼ 1.8σ with Planck + DESI) to moderate (∼ 3σ with Planck + DESI + SH0ES-calibrated SN) significance. The bΛCDM has Planck + DESI χ 2 values equal to or better than ΛCDM while predicting a higher Hubble constant. Future high-resolution CMB temperature and polarization measurements will be crucial for confirming or further constraining PMFs at recombination. Field strengths of 5-10 pico-Gauss align closely with those required for cluster magnetic fields to originate entirely from primordial sources, without the need for additional dynamo amplification.

Jedamzik, Karsten [Centre National de la Recherche↗
DOE OSTI2026

Nature of transonic sub-Alfvénic turbulence and density fluctuations in the near-Sun solar wind

Context. Recent Parker Solar Probe (PSP) measurements have revealed that solar wind (SW) turbulence transits from a subsonic to a transonic regime near the Sun, while remaining sub-Alfvénic. These observations call for a revision of the existing SW models, where turbulence is considered to be both subsonic and sub-Alfvénic. Aims. In this work, we introduce a new magnetohydrodynamic (MHD) model of transonic sub-Alfvénic turbulence (TsAT). Methods. We used 3D MHD simulations initialized with parameters measured by PSP to investigate the properties of the new near-Sun SW transonic turbulent regime. We then derived a reduced set of MHD equations in the transonic sub-Alfvénic limit to interpret our numerical results. Results. Our TsAT model shows that turbulence is effectively nearly incompressible (NI) and has a 2D + slab (quasi-2D) geometry not only in the subsonic limit, but also in the transonic regime, as long as it remains sub-Alfvénic, a condition essentially enforced everywhere in the heliosphere by the strong local magnetic field. These predictions are consistent with 3D MHD simulations, showing that transonic turbulence is dominated by low-frequency quasi-2D incompressible structures, while compressible fluctuations are a minor component corresponding to low-frequency slow modes and high-frequency fast modes. Conclusions. Our new TsAT model extends existing NI theories of turbulence, and is potentially relevant for the theoretical and numerical modeling of space and astrophysical plasmas, including the near-Sun SW, the solar corona, and the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Low latency optical-based mode tracking with machine learning deployed on FPGAs on a tokamak

Active feedback control in magnetic confinement fusion devices is desirable to mitigate plasma instabilities and enable robust operation. Optical high-speed cameras provide a powerful, non-invasive diagnostic and can be suitable for these applications. Here, in this study, we process high-speed camera data, at rates exceeding 100 kfps, on in situ field-programmable gate array (FPGA) hardware to track magnetohydrodynamic (MHD) mode evolution and generate control signals in real time. Our system utilizes a convolutional neural network (CNN) model, which predicts the n = 1 MHD mode amplitude and phase using camera images with better accuracy than other tested non-deep-learning-based methods. By implementing this model directly within the standard FPGA readout hardware of the high-speed camera diagnostic, our mode tracking system achieves a total trigger-to-output latency of 17.6 μs and a throughput of up to 120 kfps. This study at the High Beta Tokamak-Extended Pulse (HBT-EP) experiment demonstrates an FPGA-based high-speed camera data acquisition and processing system, enabling application in real-time machine-learning-based tokamak diagnostic and control as well as potential applications in other scientific domains.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Investigating nonlinear magnetohydrodynamics in an optimized, reactor-scale quasi-axisymmetric stellarator

We use the new simulation capabilities of the extended-magnetohydrodynamic (MHD) code, M3D-C1, to investigate the nonlinear MHD properties of a reactor-scale quasisymmetric stellarator equilibrium. Our model captures the self-consistent evolution of the magnetic field, temperature, density, and flow profiles without imposing restrictions on the structure of the first. We include the effects of resistivity using a realistic temperature-dependent Spitzer model, along with a model for heat transport that captures the key physical characteristic, namely, strongly anisotropic diffusion in directions perpendicular and parallel to the magnetic field. We consider a quasi-axisymmetric, finite-pressure equilibrium that was optimized for self-consistent bootstrap current, quasi-symmetry, and energetic particle confinement. Our assessment finds that the equilibrium is highly unstable to interchange-like pressure-driven instabilities near the plasma edge. The initially unstable modes rapidly destabilize other modes in the direction of the N-fold rotational symmetry (toroidal, in this case). For this equilibrium, N = 2, meaning destabilization of a large number of even-numbered toroidal Fourier modes. Thus, field-periodicity is likely to be an important factor in the nonlinear MHD stability characteristics of optimized stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Three-dimensional plasmoid-mediated reconnection and turbulence in Hall magnetohydrodynamics

Plasmoid instability accelerates reconnection in collisional plasmas by transforming a laminar reconnection layer into numerous plasmoids connected by secondary current sheets in two dimensions (2D) and by fostering self-generated turbulent reconnection in three dimensions (3D). In large-scale astrophysical and space systems, plasmoid instability likely initiates in the collisional regime but may transition into the collisionless regime as the fragmentation of the current sheet progresses toward kinetic scales. Hall magnetohydrodynamics (MHD) models are widely regarded as a simplified yet effective representation of the transition from collisional to collisionless reconnection. However, plasmoid instability in 2D Hall MHD simulations often leads to a single-X-line reconnection configuration, which significantly differs from fully kinetic particle-in-cell simulation results. This study shows that single-X-line reconnection is less likely to occur in 3D compared to 2D. Moreover, depending on the Lundquist number and the ratio between the system size and the kinetic scale, Hall MHD can also realize 3D self-generated turbulent reconnection. We analyze the features of the self-generated turbulent state, including the energy power spectra and the scale dependence of turbulent eddy anisotropy.

Physics↗
DOE OSTI2024

Orbit-following simulations of fast-ion transport and losses due to the Alfvén eigenmode burst in the Large Helical Device

Orbit-following simulations of fast-ion transport and losses with time-dependent electromagnetic perturbations are performed to clarify the roles of Alfvén eigenmodes (AEs) and the low-frequency magnetohydrodynamic (MHD) mode observed in the kinetic-MHD hybrid simulation of AE bursts in the Large Helical Device. Fast-ion pressure profile flattening in the kinetic-MHD hybrid simulation can be reproduced by an orbit-following simulation with only the primary single AE of the time-dependent amplitude following the kinetic-MHD hybrid simulation result, while orbit-following simulations with constant AE amplitude of average level during AE burst cannot reproduce the fast-ion pressure profile flattening observed. The effects of other modes are negligible on the fast-ion pressure profile flattening. The fast-ion losses in kinetic-MHD hybrid simulation can be reproduced by an orbit-following simulation with time-dependent amplitude when the low-frequency MHD mode is considered in addition to multiple AEs. This indicates the synergetic effect of multiple AEs and the low-frequency MHD mode on fast-ion losses.

Seki, Ryohsuke↗
DOE OSTI2024

Forming a database to study reversed magnetic shear from the National Spherical Torus eXperiment using machine learning

Achieving a long-lived reversed magnetic shear (RMS) target plasma in the National Spherical Torus eXperiment Upgrade will require developing various sustainment scenarios. To help with the ongoing plasma control efforts, the development of a new analysis for the motional Stark effect (MSE) diagnostic using a machine learning algorithm, namely, MSE-ML, is described. MSE-ML will be used to identify patterns during RMS discharges, some of which suffer magnetohydrodynamic (MHD) events resulting in current redistribution and monotonic q-profiles. A database consisting of q and magnetic shear profiles is being constructed primarily based on the existing National Spherical Torus eXperiment data with equilibrium reconstructions constrained by the magnetic field pitch angle profile measured using the multi-channel MSE diagnostic. An unsupervised k-means clustering of the data is developed to study the RMS formation as a function of time. The initial clustering from the q-profiles shows significant differences in both amplitude and the duration of the RMS period. As a goal, the clustering results that detect and distinguish shots with substantial and sustained RMS are to be used as a preprocessing step in a supervised algorithm to identify the underlying conditions that lead to long-lasting improved confinement with RMS. Another aim of the MSE-ML study is to identify precursors of RMS-destroying MHD events in either derived data such as the q-profile or directly measured data such as the magnetic field pitch angle profile.

Uzun-Kaymak, I. U. (ORCID:0000000276251493)↗
DOE OSTI2024

Effects of zero and reversed magnetic shear on resistive wall modes in a straight limiter tokamak with zero beta

Advanced tokamak scenarios often feature equilibriums with zero and reversed magnetic shear. To isolate and investigate their impacts on the resistive wall mode (RWM) instability analytically, we construct a series of cylindrical limiter equilibriums with different magnetic shear in the core and toward plasma edge, as a prototype of the configurations in advanced tokamak scenarios. Equilibriums with zero β are considered so that the analysis focuses on the current-driven RWMs. Analytical solutions for the n=1 resistive wall mode, derived from the reduced ideal magnetohydrodynamics (MHD) equations, consistently demonstrate key insights for both the core and the edge: the growth rate of RWMs is inversely proportional to the averaged magnetic shear when the safety factor at the edge (qa) is held constant. Specifically, an increase in positive averaged magnetic shear results in a decreased growth rate, whereas an increase in negative averaged magnetic shear leads to an increased growth rate. Furthermore, the influence of nqa on the mode's stability significantly surpasses the effects of magnetic shear changes previously discussed. These analytical calculation results provide a verification of NIMROD simulations.

Wan, Sui (ORCID:0009000464121677)↗
DOE OSTI2025

Burn propagation in magnetized high-yield inertial fusion

Recent experiments at the National Ignition Facility (NIF) have demonstrated ignition for the first time in an inertial confinement fusion (ICF) experiment, a major milestone allowing the possibility of high energy gain through burn propagation. Use of external magnetic fields, applied primarily to reduce thermal losses, could increase hotspot temperature and ease requirements for ignition, opening up the capsule design space for high energy gain. However, this same restriction of thermal transport has the potential to inhibit burn propagation, which is vital in the attainment of high gain. In this work, radiation-magnetohydrodynamics (MHD) simulations carried out using the code Chimera are used to investigate the effect of a pre-imposed magnetic field on ignition and burn propagation. This paper studies the propagation of burn using both an idealized planar model and in fully integrated 2D MHD simulations of an igniting NIF capsule. A study of magnetized burn propagation in the idealized planar model identifies three regimes of magnetized burn propagation: (1) thermal conduction driven; (2) alpha transport driven; and (3) fully suppressed burn. Simulations of NIF shot N210808 with an applied 40 T axial field show clear indication of burn suppression perpendicular to field lines, with rapid burn observed along field lines. Implosion shape is altered by the field, and anisotropic conduction causes significant modification to the rate of ablation during stagnation. These results highlight the fundamental changes to implosion dynamics in high-yield magnetized ICF and motivate further study to better optimize future magnetized target designs for high gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Nonlinear magnetohydrodynamic modeling of ideal ballooning modes in high- β Wendelstein 7-X plasmas

We present nonlinear magnetohydrodynamic (MHD) simulations of high- β Wendelstein 7-X plasmas using the stellarator extension of the M3D- C 1 code, building on the recent work that shows benign saturation of ideal ballooning modes above the designed β limit in the standard configuration [Zhou et al., Phys. Rev. Lett. 133, 135102 (2024)]. First, we examine the results' sensitivity to the parallel thermal conductivity. It is found that while an increased parallel conductivity reduces the linear growth rate, the saturated pressure profile is barely affected. Second, we consider the dependence on the profile shape. It is shown that an equilibrium with a peaked pressure profile and lower β is subject to more significant change than a broad profile with higher β and a larger growth rate, suggesting that benign saturation, or nonlinear stability, is not guaranteed and not dictated by linear growth. Third, we study the influence of the magnetic configuration, with the equilibrium rotational transform varied by adjusting the planar coil current. With similar growth rates, similar magnitudes of profile change are found regardless of the presence of a low-order resonance, which implies that the saturation mechanism is not specific to a resonant or non-resonant mode. These results indicate that MHD stability should still be treated seriously in stellarator operation and design, for which nonlinear modeling using tools like M3D- C 1 can play an instrumental role.

Zhou, Yao [Shanghai Jiao Tong University (China)] ↗
DOE OSTI2026

Expansion-driven Weibel instability in magnetohydrodynamics: Linear theory for static systems

The Weibel dispersion relation is obtained in the magnetohydrodynamic (MHD) approximation by including a tensor expression for temperature. MHD gives an upper cutoff wavenumber for Weibel growth identical to kinetic theory if the electron drift velocity is included in the evolution of the electron temperature. MHD overestimates growth rates compared to kinetic theory and gives maximum growth at larger wavenumbers, but it only leads to divergent results as temperature anisotropy tends to infinity. Thermal conduction in MHD lowers the growth rates and shifts maximum growth to smaller wavenumbers, but flux-limited or nonlocal thermal conduction is found to significantly limit these effects. The results are used as a first step in evaluating the potential of MHD to simulate expansion-driven Weibel instability, which occurs due to the temperature tensor only cooling in the direction of expansion. Finally, the small temperature anisotropy found in expansion-driven Weibel instability means that MHD could be an adequate model because growth rates are much less than the electron plasma frequency, and wavelengths are much greater than the Debye length.

Davies, J. R. [Univ. of Rochester, NY (United Stat↗
DOE OSTI2026

Multiscale and multidimensional modeling of particle acceleration and transport in solar flares

Multi-messenger, multi-viewpoint, and time-resolved observations of solar flares are now providing unprecedented constraints on particle acceleration sites, energy conversion, and energy transport. The interpretation of current observations, including microwave imaging spectroscopy from EOVSA, hard x-ray (HXR) imaging from Solar Orbiter/STIX, gamma-ray diagnostics from Fermi, and in situ measurements from Parker Solar Probe and Solar Orbiter, collectively demands modeling frameworks that go beyond traditional spatially unresolved, one-zone models or single-mechanism descriptions. This review surveys multiscale and multidimensional modeling approaches, including kinetic, magnetohydrodynamic (MHD), and macroscopic particle models, that are being developed to meet the need. Kinetic simulations reveal that three-dimensional (3D) effects, including field-line chaos and self-generated turbulence, are essential for sustained power-law particle acceleration. MHD simulations now capture flux-rope eruptions, plasmoid-unstable current sheets, and turbulent flare regions in realistic magnetic topologies. Macroscopic models coupling MHD with energetic-particle models produce spatially resolved electron distributions and synthetic HXR and microwave emissions for direct comparison with observations. Despite these advances, outstanding challenges remain in bridging kinetic and global scales, improving MHD simulations and macroscopic particle models, and achieving quantitative model-observation closure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design and Analysis of the Open-Surface Slow Li Flow Divertor and Comparison to the Fast Li Flow Divertor

In the ongoing U.S. project, “Liquid Metal Plasma Facing Components,” sponsored by the U.S. Department of Energy, efforts have been taken to develop two open-surface divertor designs for the Fusion Nuclear Science Facility using liquid lithium (Li) as a heat and particle flux removal media. The main focus of this study is the design and analysis of a slow (~1 mm/s) and thin (<1 mm) open-surface Li flow divertor with a Li-cooled substrate, which is then compared with an earlier design of a fast (up to 10 m/s) and thick (~0.5 cm) Li flow divertor with the substrate cooled with helium. The slow Li flow divertor design is based on the original LiWall concept developed at the Princeton Plasma Physics Laboratory. Such a thin and slow Li layer can remove the particle flux by reducing the recycling flux, while the heat flux is removed mainly through the heat sink located beneath. In the present study, the heat sink is provided through a Li cooling flow inside the substrate of reduced activation ferritic/martensitic steel. By performing a multiphysics analysis with COMSOL that included liquid-metal magnetohydrodynamics (MHD), heat transfer, and structural mechanics, the impact of various factors on the divertor heat removal capability, such as Li flow velocity, MHD effects, and inlet velocity boundary condition, were examined. Based on comparisons of the two divertor designs, it was shown that the fast-flow divertor significantly outperformed the slow-flow design, whose heat removal capability was limited to ~1 to 2 MW/m2.

Jiang, Yuchen [ORNL] (ORCID:0000000283280284)↗
DOE OSTI2025

Hybrid Simulations of FRC Merging and Compression

An improved understanding of field-reversed configuration (FRC) merging and stability in high acceleration and compression magnetic fields is needed to speed up the development of the pulsed fusion concept developed at Helion Energy. All previous theoretical and simulation work on FRC merging and compression was performed using two-dimensional (2D) magnetohydrodynamic (MHD) models. The results of novel 2D hybrid simulations (fluid electrons and full-orbit kinetic ions) of FRC merging and compression are presented. Results of kinetic and MHD simulations, computed using the HYM code, are compared and analyzed. In cases without axial magnetic compression, both the MHD and hybrid simulations show a high sensitivity to the initial parameters (i.e. FRC separation, velocity, normalized separatrix radius, and plasma viscosity), showing that FRCs with large elongation and separatrix radius either do not merge or merge partially, forming a doublet FRC. In conclusion, application of a mirror coil field at the FRC ends with increasing strength is shown to lead to fast and complete merging of the FRCs in MHD and kinetic simulations.

FRC↗
DOE OSTI2024

Accuracy of kinetic equilibrium reconstruction of NSTX and NSTX-U plasmas and its impact on the transport and stability analysis

An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

ELM-free H-mode phase and decoupling of peeling–ballooning stability boundary in the MAST Upgrade tokamak

A linear magnetohydrodynamic (MHD) peeling–ballooning stability analysis of the edge-localized mode (ELM)-free phase of a MAST Upgrade (MAST-U) H-mode plasma is presented. In contrast to other similar discharges, #47018 is found to have a significantly higher and wider pedestal during its ELM-free H-mode phase that lasts for approximately 80 ms; this is made possible by the reduced core MHD mode activity on the q = 2 surface. During this period, there is sustained decoupling of peeling and ballooning branches of the stability boundary on J–α space, opening an access channel to the second stability regime with higher peaks in pedestal current density $J_{\mathrm{N,ped}}$ and pressure gradient (α). Such decoupling of the stability boundary has not previously been observed in MAST-U H-modes, and if such a condition can be readily reproduced, it opens a wide range of opportunities for MAST-U to explore low-collisionality peeling-limited pedestal regimes as well as advanced scenarios such as quiescent H-modes that are relevant to future reactors such as STEP and ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Numerical investigation of two-dimensional buoyancy-driven eddies in liquid metal magnetohydrodynamic flows in breeding blankets *

Abstract Lead–lithium flows are key features in the design of tokamak breeding blanket concepts such as the dual-coolant lead–lithium (DCLL). Since they flow under magnetic fields, they are affected by magnetohydrodynamic (MHD) effects. The neutron flux originating in the tokamak plasma heats the breeding blanket channels in a non-uniform manner, inducing buoyancy forces in the liquid metal. Buoyancy may become a source of quasi-two-dimensional (Q2D) turbulence, and the appearance of eddies may affect the transport of heat and tritium across the blanket. Blankets characterized by high-speed liquid metal flows (such as DCLL) will need ceramic insulating walls to reduce the MHD-related pressure drop in the channels. In our simulations, we have used the Q2D model proposed by Sommeria and Moreau (SM82) which is especially suitable for modeling electrically insulating channel flows. On top of that, we have modelled buoyancy forces in the momentum equation using the Oberbeck–Boussinesq approximation. In this work, we include a validation of the implemented Q2D model in buoyancy-driven cases and identify a thermohydraulic configuration that promotes the generation of eddies and the accumulation of tritium. We also show the results provided by our post-processing tool based on the bi-dimensional fast Fourier transform for eddy detection and characterization. We complete our investigation by performing an initial assessment of how eddies can accumulate tritium in the breeding blanket. We conclude by discussing the relevance of preventing accumulation of tritium within the blanket and suggest a possible solution.

Suarez, Daniel (ORCID:0000000251740100)↗
DOE OSTI2025

Low- n stability and plasma response to RMP in various STEP scenarios

The low-n (n is the toroidal mode number) magnetohydrodynamic (MHD) stability and plasma response are numerically investigated for various scenarios designed for STEP, that are relevant for the H-mode pedestal analysis. Control of the edge-localized modes (ELMs) with externally applied resonant magnetic perturbations (RMPs) is considered. Optimization of the ELM control coil current configuration, based on the computed plasma MHD response and well-established figures of merit validated on present-day experiments, finds reasonable robustness of a fixed coil phasing (for a given n-number) to control ELMs in all five STEP plasmas considered. Based on certain semi-empirical criteria, the required coil current to achieve ELM suppression is estimated to be about 10–20 kAt with the n = 1 or 2 RMP configuration and about 100–200 kAt for the n = 4 RMP. Systematic linear stability calculations are used to map out stability windows for the low-n kink-peeling modes, in terms of the ideal-wall location and variation of the edge safety factor q 95 with respect to the target design. The kink-peeling stability boundary is found to be generally sensitive to the q 95 variation, which has implications for achieving the quiescent H-mode regime in STEP. Full toroidal quasilinear initial-value simulations for these STEP plasmas find that generation of the edge-harmonic oscillations (EHOs) depends sensitively on the plasma scenario, the initial linear stability of the kink-peeling modes, the initial plasma toroidal flow and q 95 . In general, it is easier (more robust) to access the EHO-regime for two of the cases considered with smaller plasma volume and higher on-axis safety factor. Finally, quasilinear simulations find robust density pumpout due to applied RMPs in these STEP plasmas, but the effect on the plasma toroidal flow varies among different cases.

EHO↗