Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA LAYER”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Analytical model of plasma response to external magnetic perturbation in absence of no-slip condition

Recent simulation and experimental results suggest that the magnetic island and flow on a resonant surface often do not satisfy the “no-slip” condition in the steady state. A new theoretical model on a nonlinear plasma response to external magnetic perturbation in the absence of a no-slip condition is proposed. The model is composed of the equations for the evolution of both the width and phase of the magnetic island due to forced reconnection driven by the external magnetic perturbation, and the force-balance equation for the plasma flow. When the island width is much less than the resistive layer width, the island growth is governed by the linear Hahm–Kulsrud–Taylor solution in the presence of time-dependent plasma flow. In the other regime when the island width is much larger than the resistive layer width, the evolution of both the island width and phase can be described using the Rutherford theory. The island solution is used to construct the quasilinear electromagnetic force, which, together with the viscous one, contributes to the nonlinear variation in plasma flow. The no-slip condition assumed in the conventional error field theory is not imposed here, where the island oscillation frequency depends on but does not necessarily equal to the plasma flow frequency at the rational surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

In-Situ Laser-Induced-Desorption-Spectroscopy Diagnostic Deployment on the Linear Plasma Device CSDX

A calibrated laser induced desorption diagnostic (LIDS) concept is developed under this award for application to fusion relevant plasma materials interactions studies. Small areal amounts of hydrogenic release from a deuterium filled PMI target are measured and quantified when irradiated by a LID flash heat event on the target surface. Good success was demonstrated by this project with the PISCES-E rf plasma device at UCSD by fully quantified measurement and removal of D from W-D co-deposit layers in experiments. However, adapting the LIDS technique to D plasma exposed W targets in the linear plasma CSDX device, now labelled PISCES-Rf, was prone to several technical challenges ranging from target line of sight access, laser focusing, and loss of optical spectroscopic signal strength in the larger PISCES-Rf vacuum system. The results are summarized in this report.

08 HYDROGEN↗
DOE OSTI2026

Upgrades to the Thomson Scattering Diagnostic on the Lithium Tokamak eXperiment-Beta (LTX-$β$)

We describe the upgrades to the Thomson scattering (TS) diagnostic on the Lithium Tokamak eXperiment-beta (LTX- β ) to provide complete radial coverage of the plasma. LTX- β is a spherical tokamak with plasma limited at a major radius of R = 14 cm and R = 66 cm. Flat temperature profiles with peaked density were observed in LTX- β using an existing multipoint TS diagnostic with radial coverage of the plasma core. There is a need for more comprehensive core and edge TS measurements to thoroughly examine beam fueling, beam heating, and transport under low recycling conditions with liquid lithium walls and a lower collisionality scrape-off layer (SOL). The core TS set up provides plasma profiles from R = 40–60 cm. To extend the TS diagnostic coverage: 1) we added a TS viewing optics to cover the SOL region at five radial points from R = 57–64 cm, and 2) another TS viewing optics in the HFS region from R = 15–45 cm to obtain TS data from ten radial positions while making use of the original laser beam. The main focus of this article is on the new HFS TS setup that is equipped with better-matched collection optics located in the horizontal midplane where the laser polarization lies and views from a scattering angle of 130°–150° with respect to the laser path. A set of unequal-length Y-branched fiber bundles transports the TS light from two radial points to each polychromator. The new HFS view is better matched in spot size, scattering angle, and throughput to the laser and polychromators and has new V-grooved style viewing dumps. The new HFS TS view, along with the core and SOL TS views, will provide a full radial measurement coverage from R = 15–64 cm to explore the flat temperature profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Machine learning-based real-time kinetic profile reconstruction in DIII-D

Abstract Kinetic equilibrium reconstruction plays a vital role in the physical analysis of plasma stability and control in fusion tokamaks. However, the traditional approach is subjective and prone to human biases. To address this, the consistent automatic kinetic equilibrium reconstruction (CAKE) method was introduced, providing objective results. Nonetheless, its offline nature limits its application in real-time plasma control systems (PCSs). To address this limitation, we present RTCAKENN, a machine learning model that approximates 7 CAKE-level output profiles, namely pressure, inverse q , toroidal current density, electron temperature and density, carbon ion impurity temperature and rotation profiles, using real-time available inputs. The deep neural network consists of an encoder layer, where the scalars and interdependent inputs such as plasma boundary coordinates and motional Stark effect data are encoded using multi-layer perceptrons (MLPs), while profile inputs are encoded by 1D convolutional layers. The encoded data is passed through a MLP for latent feature extraction, before being decoded in the decoding layers, which consist of upsampling and convolutional layers. RTCAKENN has been implemented in the DIII-D PCS and our model achieves accuracy comparable to CAKE and surpasses existing real-time alternatives. Through clever dropout training, RTCAKENN exhibits robustness and can operate even in the absence of Thomson scattering data or charge exchange recombination data. It executes in under 8 ms in the real-time environment, enabling future application in real-time control and analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Epitaxially Grown Single-Crystalline SrTiO 3 Membranes Using a Solution-Processed, Amorphous SrCa 2 Al 2 O 6 Sacrificial Layer

Water-soluble sacrificial layers based on epitaxially-grown, single crystalline (Ca, Sr, Ba) 3 Al 2 O 6 layer are widely used for creating free-standing perovskite oxide membranes. However, obtaining these sacrificial layers with intricate stoichiometry remains a challenge, especially for molecular beam epitaxy (MBE). In this study, we demonstrate the hybrid MBE growth of epitaxial, single crystalline SrTiO 3 films using a solution processed, amorphous SrCa 2 Al 2 O 6 sacrificial layer onto SrTiO 3 (001) substrates. Prior to the growth, the oxygen plasma exposure was used to first create the crystalline SrCa 2 Al 2 O 6 layer with well-defined surface crystallinity. Utilizing reflection high energy electron diffraction, x-ray diffraction, and atomic force microscopy, we observe an atomic layer-by-layer growth of epitaxial, single crystalline SrTiO 3 film on the SrCa 2 Al 2 O 6 layer with atomically smooth surfaces. The SrCa 2 Al 2 O6 layer was subsequently dissolved in de-ionized water to create free-standing SrTiO3 membranes that were transferred onto a metal-coated Si wafer. Membranes created with Sr-deficiency revealed ferroelectric-like behavior measured using piezo force microscopy whereas stoichiometric films remained paraelectric-like. Furthermore, these findings underscore the viability of using ex-situ deposited amorphous SrCa2Al2O6 for epitaxial, single crystalline growth, as well as the importance of point defects in determining the ferroic properties in membranes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Suppression of Midinfrared Plasma Resonance Due to Quantum Confinement in δ-Doped Silicon

The classical Drude model provides an accurate description of the plasma resonance of three-dimensional materials, but only partially explains two-dimensional systems where quantum mechanical effects dominate such as P:δ layers—atomically thin sheets of phosphorus dopants in silicon that induce electronic properties beyond traditional doping. Previously it was shown that P:δ layers produce a distinct Drude tail feature in ellipsometry measurements. However, the ellipsometric spectra could not be properly fit by modeling the δ layer as a discrete layer of classical Drude metal. In particular, even for large broadening corresponding to extremely short relaxation times, a plasma resonance feature was anticipated but not evident in the experimental data. In this work, we develop a physically accurate description of this system, which reveals a general approach to designing thin films with intentionally suppressed plasma resonances. Our model takes into account the strong charge-density confinement and resulting quantum mechanical description of a P:δ layer. Here, we show that the absence of a plasma resonance feature results from a combination of two factors: (i) the sharply varying charge-density profile due to strong confinement in the direction of growth; and (ii) the effective mass and relaxation time anisotropy due to valley degeneracy. The plasma resonance reappears when the atoms composing the δ layer are allowed to diffuse out from the plane of the layer, destroying its well-confined two-dimensional character that is critical to its distinctive electronic properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Deuterium retention in pre-lithiated samples and Li–D co-deposits in the DIII-D tokamak

Divertor designs involving liquid lithium have been proposed as an alternative to solid designs and wall conditioning techniques. However, Li affinity with tritium poses a risk for the fuel cycle. This study investigates deuterium retention in pre-lithiated samples and Li–D co-deposits in the DIII-D tokamak, making for the first time a direct comparison between Li–D co-deposits and pre-deposited Li films. Samples were exposed to H-mode plasmas in the far scrape-off layer (SOL), and Li powder was injected in-situ with the impurity powder dropper to study the uniformity of Li coatings, and the dependence of fuel retention on Li thickness. The results show that at temperatures below the melting point of lithium, deuterium retention is independent of the thickness of pre-deposited Li layers, with Li–D co-deposits being the primary factor for fuel retention. Both pre-deposited and in-situ deposited Li showed lower erosion than predicted by sputtering yield calculations. These results suggest that fuel retention in fusion reactors using lithium in the divertor will likely be dominated by co-deposits rather than in the divertor itself. If one desires to use Li to achieve flatter temperature profiles, operando Li injection is advantageous over pre-deposited Li films, at least at temperatures below the melting point of lithium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Reimagining full wave rf quasilinear theory in a tokamak

The velocity dependent resonant interaction of particles with applied radiofrequency (rf) waves during heating and current drive in the presence of pitch angle scattering collisions gives rise to narrow collisional velocity space boundary layers that dramatically enhance the role of collisions as recently shown by Catto ( J. Plasma Phys. , vol. 86, 2020, 815860302). The behaviour is a generalization of the narrow collisional boundary layer that forms during Landau damping as found by Johnston ( Phys. Fluids , vol. 14, 1971, pp. 2719–2726) and Auerbach ( Phys. Fluids , vol. 20, 1977, pp. 1836–1844). For a wave of parallel wave number ${k_{||}}$ interacting with weakly collisional plasma species of collision frequency $\nu$ and thermal speed ${v_{\textrm{th}}}$ , the effective collision frequency becomes of order $\nu {({k_{||}}{v_{th}}/\nu )^{2/3}} \gg \nu $ . The narrow boundary layers that arise because of the diffusive nature of the collisions allow a physically meaningful wave–particle interaction time to be defined that is the inverse of this effective collision frequency. The collisionality implied by the narrow boundary layer results in changes in the standard quasilinear treatment of applied rf fields in tokamaks while remaining consistent with causality. These changes occur because successive poloidal interactions with the rf are correlated in tokamak geometry and because the resonant velocity space dependent interactions are controlled by the spatial and temporal behaviour of the perturbed full wave fields rather than just the spatially local Landau and Doppler shifted cyclotron wave–particle resonance condition associated with unperturbed motion of the particles. The correlation of successive poloidal circuits of the tokamak leads to the appearance in the quasilinear operator of transit averaged resonance conditions localized in velocity space boundary layers that maintain negative definite entropy production.

Physics↗
DOE OSTI2021

A pushered capsule implosion as an alternate approach to the ignition regime for inertial confinement fusion

We report in inertial confinement fusion, the threshold for ignition is a highly dynamic quantity as the sources and sinks of power in the hot spot can vary rapidly. In this article, we consider the ignition condition as a race between heating and disassembly rates and make use of a prior solution to the fusion hot-spot thermodynamics to develop a Lawson-like ignition criteria for pressure × confinement time (p-τ) vs temperature. Low-Z capsule designs reach the temperature for this threshold using as much of the shell as feasible as ablator but then are limited in τ by low stagnated mass. An alternate approach, the pushered single shell (PSS) design [D. D.-M. Ho, S. MacLaren, and Y. Wang, “High-yield implosions via radiation trapping and high rho-R,” paper presented at the 60th Annual Meeting of the APS Division of Plasma Physics, 2018], introduces a dense inner layer of Mo-Be alloy that is smoothly graded outward to pure Be, increasing the confinement time at stagnation and lowering the temperature requirement at the ignition threshold. Here, we describe a PSS ignition design for the National Ignition Facility and use the theory as well as simulations to compare it with the low-Z capsule approach. Additionally, we show how an adjustment to the design is used to anticipate the effects of mixing at the fuel–ablator interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Role of the edge stochastic layer in density pump-out by resonant magnetic perturbations

The role of the edge stochastic layer on particle transport is addressed in DIII-D plasmas with applied resonant magnetic perturbations (RMPs) causing density pump-out (a term used to denote the density reduction at the top of the pedestal due to the applied RMPs). Using an analytical model that adds the stochastic parallel transport of electrons (Harvey et al 1981 Phys. Rev. Lett. 47 102) to the fluid equations, the ambipolar radial electric field and particle flux are calculated simultaneously. In this model the nonambipolar electron particle flux, driven by the stochastic magnetic field, is predominantly balanced by the nonambipolar perpendicular ion flux, driven by toroidal viscosity, across a narrow stochastic layer of order two percent of the plasma radius (0.98 < ψn < 1). The model reproduces the level of density pump-out and its dependence on RMP amplitude near the separatrix for three plasma discharges for which the pedestal foot is at medium to high collisionality $(2\leqslant {\nu }_{\text{sep}}^{\ast }\leqslant 30)$. In this work, the experimentally observed and numerically tested inverse density dependence of density pump-out (Hu et al 2020 Nucl. Fusion 60 076001) is accurately captured by the model: an increase in collisionality with density in the pedestal foot results in a decrease in stochastic diffusivity, and hence a decrease in the level of pump-out.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The effects of multispecies Hohlraum walls on stimulated Brillouin scattering, Hohlraum dynamics, and beam propagation

We report experiments and simulations have been conducted to investigate the efficacy of Ta 2 O 5 -lined Hohlraum walls at reducing stimulated Brillouin backscattering (SBS) as well as any subsequent effects on the Hohlraum dynamics and capsule implosions in indirect drive experiments at the National Ignition Facility. Using a 1.1 MJ 400 TW, 351 nm, shaped laser pulse, we measure a 5× reduction in SBS power in the peak of the pulse from the wall on the outer 50° cone beams. The SBS spectrum indicates a reduction in the high-Z spectral signature when using multispecies wall materials. Detailed hydrodynamic simulations were performed using different heat conduction models with flux limiters. Additional simulations were performed on the plasma maps using the 3D parallel paraxial code pF3D to compare backscatter powers between the pure Au and Ta 2 O 5 -lined Hohlraums. Further analysis, using hydrodynamically equivalent plasmas, shows that the SBS reduction is clearly a result of the added ion Landau damping caused by the oxygen ions and not from differences in plasma conditions. The experimental and simulation results also show an increase in the wall plasma expansion when using the Ta 2 O 5 liner leading to a 70% more oblate implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Inferring the scrape-off layer heat flux width in a divertor with a low degree of axisymmetry

Plasma facing components (PFCs) in the next generation of tokamak devices will operate in challenging environments, with heat loads predicted to exceed 10 MWm -2 . The magnitude of these heat loads is set by the width of the channel, the ‘scrape-off layer’ (SOL), into which heat is exhausted, and can be characterised by an e-folding length scale for the decay of heat flux across the channel. It is expected this channel will narrow as tokamaks move towards reactor relevant conditions. Understanding the processes involved in setting the SOL heat flux width is imperative to be able to predict the heat loads PFCs must handle in future devices. Measurements of the SOL width are performed on the high-field spherical tokamak, ST40, using a newly commissioned infrared thermography system. With its high on-axis toroidal magnetic field (≥1.5 T) ST40 is uniquely positioned to investigate the influence of toroidal field on the heat flux width in spherical tokamaks, whilst also extending measurements of the SOL width in spherical tokamaks to increased poloidal field (≥0.3 T). Due to the divertor on ST40 having a low degree of axisymmetry, it is necessary for a set of radial measurements of the heat flux to be taken across the divertor, made possible using an automated toolchain that fully incorporates its 3D geometry. These radial profiles are combined with the magnetic geometry of the plasma to infer the width of the SOL, with both Eich and double exponential profiles of heat flux observed. A reduction in the heat flux is observed toroidally across part of the divertor, along with increased heat loads observed locally around the edges of the tiles. Future work in characterising the impact of tile misalignment and uncertainties in the reconstructed divertor magnetic geometry is required in order to further understand the observed heat flux patterns, as are additional investigations into the role potentially being played by an inhomogeneous sheath electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Nitrogen plasma passivated niobium resonators for superconducting quantum circuits

Microwave loss in niobium metallic structures used for superconducting quantum circuits is limited by a native surface oxide layer formed over a timescale of minutes when exposed to an ambient environment. Here, in this work, we show that nitrogen plasma treatment forms a niobium nitride layer at the metal–air interface, which prevents such oxidation. X-ray photoelectron spectroscopy confirms the doping of nitrogen more than 5 nm into the surface and a suppressed oxygen presence. This passivation remains stable after aging for 15 days in an ambient environment. Cryogenic microwave characterization shows an average filling-factor-adjusted two-level-system loss tangent Fδ TLS of (2.9±0.5)·10 –7 for resonators with a 3 μm center strip and (1.0±0.3)·10 –7 for a 20 μm center strip, exceeding the performance of unpassivated samples by a factor of four.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Correlation between the relative blob fraction and plasma parameters in NSTX

Here, we present a simple analysis of the fraction of blobs observed in the scrape-off layer of NSTX and calculates the correlation between this blob fraction and various plasma parameters. The measurements were made using the gas puff imaging (GPI) diagnostic on a database of 103 shots from 2010 with neutral beam power P NBI = 0–6 MW. The blob fraction ranged from f b = 0.1% to 4.8%, where f b is defined as the fraction of time the GPI signal was larger than three times its average value within 0–7.5 cm outside the separatrix. The blob fraction generally decreased with increasing neutral beam injection power and was lower on average in H-mode plasmas than Ohmic and L-mode plasmas. The blob fraction had its highest correlation with the measured poloidal turbulence velocity and its radial gradient just inside the separatrix. This is in part consistent with a model in which the blob generation rate increases with the poloidal velocity shear.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Magnetic mesh generation and field line reconstruction for scrape-off layer and divertor modeling in stellarators

The design of divertor targets and baffles for optimal heat and particle exhaust from magnetically confined fusion plasmas requires a combination of fast, low-fidelity models (such as EMC3-Lite (Feng 2022 Plasma Phys. Control. Fusion 64 125012)) for scoping studies and high-fidelity ones (such as EMC3-EIRENE (Feng et al 2014 Contrib. Plasma Phys. 54 426–31)) for verification. Both of those approaches benefit from a magnetic flux tube mesh for fast interpolation and mapping of field line segments (Feng et al 2005 Phys. Plasmas 12 052505). A new automated mesh generator for unstructured quadrilateral flux tubes with adaptive refinement is presented and integrated into FLARE (Frerichs 2024 Nucl. Fusion 64 106034). For HSX with an extended first wall, it is found that several layers of flux tubes can span the entire half field period before splitting is required. This is an advantage over the traditional setup of the EMC3-EIRENE mesh where careful construction of several sub- domains is required already for the much tighter present first wall. In particular, there is no longer the need to manually construct a suitable outer boundary for the mesh. The divide and conquer paradigm with unstructured mesh layout offers a powerful alternative for fast head load approximation that is suitable for integration into optimization workflows. Further examples for W7-X and CTH demonstrate the versatile application range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

AlGaN/GaN metal–insulator–semiconductor high electron mobility transistors (MISHEMTs) using plasma deposited BN as gate dielectric

AlGaN/GaN metal–insulator–semiconductor high electron mobility transistors (MISHEMTs) were fabricated on Si substrates with a 10 nm boron nitride (BN) layer as a gate dielectric deposited by electron cyclotron resonance microwave plasma chemical vapor deposition. The material characterization of the BN/GaN interface was investigated by X-ray photoelectric spectroscopy (XPS) and UV photoelectron spectroscopy. The BN bandgap from the B1s XPS energy loss is ~5 eV consistent with sp2 bonding. The MISHEMTs exhibit a low off-state current of 1 × 10 –8 mA/mm, a high on/off current ratio of 10 9 , a threshold voltage of –2.76 V, a maximum transconductance of 32 mS/mm at a gate voltage of –2.1 V and a drain voltage of 1 V, a subthreshold swing of 69.1 mV/dec, and an on-resistance of 12.75 Ω·mm. The interface state density (D it ) is estimated to be less than 8.49 × 10 11 cm –2 eV –1 . Gate leakage current mechanisms were investigated by temperature-dependent current–voltage measurements from 300 K to 500 K. The maximum breakdown electric field is no less than 8.4 MV/cm. Poole–Frenkel emission and Fowler–Nordheim tunneling are indicated as the dominant mechanisms of the gate leakage through the BN gate dielectric at low and high electric fields, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Proton Selective Nanoporous Atomically Thin Graphene Membranes for Vanadium Redox Flow Batteries

Angstrom-scale proton-selective pores in atomically thin 2D materials present fundamentally new opportunities for advancing proton exchange membranes (PEMs). Vanadium Redox Flow Batteries (VRFBs) for grid-scale energy storage require PEMs with high areal proton conductance (>1 S cm −2 ) and minimal vanadium ion (VO 2+ ) crossover. However, state-of-the-art Nafion 212 membranes (N212 ≈50 µm thick), suffer from persistent VO 2+ crossover reducing performance and efficiency. Here, a layered PEM is demonstrated, comprising monolayer CVD graphene with Angstrom-scale proton-selective pores introduced via Ar plasma, integrated with an ultra-thin ≈300 nm polybenzimidazole (PBI) layer and sandwiched between two Nafion 211 (25 µm thick) layers. The layered architecture facilitates scalable membrane fabrication by mitigating defects while processing and facile stacking of graphene layers allows stochastic non-selective defect isolation enabling exceptionally low VO 2+ crossover (selectivity (H + areal conductance / VO 2+ permeability) ≈6709 × 10 6 S min cm −4 ), with proton conductance >8 S cm −2 . Systematic transport experiments supported by resistance-based transport modelling elucidate the role of defect size, defect isolation, and sealing, as well as layering/stacking, to enable orders of magnitude (>671× over N212) improvements in selectivity, along with areal proton conductance >8 S cm −2 . This work highlights the potential of atomic-scale proton-selective defect engineering in 2D materials, in conjunction with facile stacking and layering of materials as strategies for scalable, high-performance advances in PEMs for energy, electrochemical, and separation applications beyond VRFBs.

Chaturvedi, Pavan [Vanderbilt Univ., Nashville, TN↗
DOE OSTI2021

Computational general relativistic force-free electrodynamics: II. Characterization of numerical diffusivity

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe’s most extreme flows of energy. The discretization of Maxwell’s equations – needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling – introduces numerical diffusion. It acts as a source of dissipation independent of the system’s physical constituents. Understanding the numerical diffusion of scientific codes is the key to classifying their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field (j ∥ ), in combination with the breakdown of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g., the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of j ∥ .

79 ASTRONOMY AND ASTROPHYSICS↗