Search NASA⌕ Search

SEARCH · Search NASA

Results for “magnetohydrodynamics (MHD)”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Magnetic Reconnection in Solar Flares [Slides]

Magnetic reconnection is the underlying mechanism in solar flares. Magnetic reconnection is an astrophysical process in an electrically conducting plasma where neighboring magnetic field lines running anti-parallel break and reconnection therefore exchanging the direction of their paths. The result of this process is the magnetic field energy converting to kinetic and thermal energy which is explosively released. The scale of a solar is extremely large compared to Earth so kinetic simulations are unable to portray an accurate model. In this project, we are using an astrophysical magnetohydrodynamic (MHD) simulation code known as Athena++ with the reconnection problem generator file created by Li et al. 2018. To further study this simulation, we began testing a heating/cooling term. The results of the original file will be compared with the modified file with an implemented heating/cooling term to better understand how solar flares function.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Physics of the thermal quench in MST tokamak plasmas (Final Technical Report)

This project at the University of Wisconsin-Madison targeted measurements in tokamak plasmas in the Madison Symmetric Torus (MST) to help elucidate the physics of the tokamak disruption. We intended to focus on the thermal quench, during which energy is lost from the plasma in a very short time. Understanding and ultimately mitigating the thermal quench are among the highest priority research thrusts for ITER, a large international fusion experiment presently under construction. The measurements in MST exploited an advanced set of diagnostics for internal measurements with high time resolution, targeting linear and nonlinear magnetohydrodynamic (MHD) activity as well as the electron temperature and density profiles. We also exploited other features of the MST device, including very flexible power supplies for the toroidal magnetic field and toroidal plasma current that are capable of producing almost arbitrary waveforms. While the work funded by this project was purely experimental, we continued and expanded an existing collaboration with UW-Madison experts in nonlinear MHD computational modeling of tokamak plasmas. That modeling incorporated realistic, experimentally measured parameters. The results of the modeling were anticipated to deepen considerably our understanding of the experimental measurements, and in turn, the experimental measurements were anticipated to help validate the modeling results. The science of transients such as disruptions is one of the highest priority research areas in the U.S. and international fusion programs. The work carried out in this project was targeted to help advance the science of tokamak disruptions which will in turn contribute to the success of ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulating of magnetic reconnection in solar flares

Magnetic reconnection is an astrophysical process where neighboring magnetic field lines running anti-parallel are reconfigured. Observations of solar flares are what pushed for serious research regarding magnetic reconnection as it seemed to be the underlying mechanism. The reconfiguration of the field lines results in an explosive release of energy as magnetic field energy is converted into plasma kinetic and thermal energies. Our work begins with Athena++, an astrophysical magnetohydrodynamic (MHD) simulation code, and an existing magnetic reconnection setup.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Optimal Chloride Salt Mixture for a Fusion Blanket

Deuterium-tritium fusion reactors cannot operate for a significant period without a closed tritium fuel cycle, according to a recent National Academies of Sciences (NAS) report on bringing fusion reactors to the US electrical grid. This fact places breeder blankets as one of the foundational systems for self-sustained fusion reactor operation. The typical functional requirements for a breeder blanket system include producing tritium, absorbing kinetic energy, transporting thermal energy, and being environmentally attractive. State-of-the-art research on liquid blankets has converged to primarily focus on (LiF) 2 and BeF 2 (FLiBe) molten salts and a metallic eutectic of lead and lithium (PbLi), but “virtually all of the technologies related to the tritium fuel cycle are at a low technological readiness level”. This work sought to explore optimum blanket configurations as it aligned with the Oak Ridge National Laboratory (ORNL) FY 2023 Laboratory Directed Research and Development Program’s research priority of developing and expanding the current understanding of fusion blanket science and technology. This purpose of this work was to address ORNL research priorities and NAS recommendations by investigating novel liquid blanket materials that could provide self-sustaining operation and draw on experience from research on molten salts used for advanced fission reactors, concentrated solar, and thermal energy storage. The hypothesis when proposing this research was that there could be chloride-based blanket designs that can exceed the tritium breeding ratios of (FLiBe) molten salt blankets while reducing the use of Be (FLiBe), avoiding the generation of HF (FLiBe), and minimizing magnetohydrodynamic (MHD)-perturbed flow fields (PbLi). The fastest and most cost-effective path to deploying liquid fusion breeder blankets could be from maximizing the synergistic technological overlap between fusion, fission, concentrated solar, and thermal energy storage industries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulating ‘Two Ribbon’ Type Solar Flares [Slides]

The underlying mechanism in solar flares in a process known as magnetic reconnection. This astrophysical process is when magnetic field lines, running anti-parallel, break and reconnect. The change in configuration of field lines results in an explosive release of energy as the leftover magnetic energy is converted to kinetic and thermal energies. In this project, we are using an astrophysical magnetohydrodynamic (MHD) simulation code known as Athena++ with a reconnection problem generator file created by Li et al. 2018. Our current research is to successfully implement a radiative cooling term into the MHD equations as it could have important physical effects on the plasma. We use Klimchuk et al. 2008 and SPEX_DM as our chosen cooling functions. The results from the SPEX_DM function show there is a condensation forming that we had predicted but had not seen before with our simulation code. We will continue to analyze the SPEX_DM function and potentially implement thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Exploration of Rayleigh-Taylor seeding mechanisms in laser- and pulsed-power-driven inertial fusion

The ubiquitous Rayleigh-Taylor instability (RTI) has been observed in fluids and plasmas in a wide range of parameter regimes from astrophysical to laboratory plasmas. In high-energy-density (HED) laboratory plasmas, such as laser-driven inertial confinement fusion (ICF) plasmas, the RTI can have a detrimental impact on achieving fusion ignition through the generation of hydrodynamic mix. Hence, an understanding of the seeding mechanisms that produce the RTI and identifying ways to mitigate hydrodynamic mix is of critical importance. Appropriately aligned magnetic fields have been known to stabilize short-wavelength RTI. A number of numerical and experimental studies have demonstrated the benefits of using imposed magnetic fields in laser-driven ICF to achieve higher ion temeratures, higher neutron yields, and a relaxation of the ignition criteria. This work addresses seeding mechanisms in laser-driven implosions to understand critical early-stage physics that ultimately leads to substantial growth of the RTI along with mechanisms for mitigation of this growth. Surface perturbations due to machining tolerances and single-feature seeds (for example due to fill tubes) can produce substantial RTI growth in the ignition-relevant high-convergence ratio targets for laser- driven ICF implosions. RTI growth from single-feature seeding in laser-based implosions has the potential to be mitigated through appropriately aligned externally applied magnetic fields and the goal of this work has been to quantify that numerically and experimentally. Resistive-magnetohydrodynamic (MHD) simu- lations are used to study the seeding and evolution of the RTI leveraging previous support from the DOE HEDLP program. Furthermore, experimental data has been obtained and applied towards code validation of unmagnetized and magnetized evolution of single-feature seeded RTI growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulating Magnetic Reconnection in ‘Two Ribbon’ Type Solar Flares

Magnetic reconnection is an astrophysical process where neighboring magnetic field lines, facing anti-parallel, are reconfigured. This reconfiguration results in built up magnetic energy being explosively released as it is being converted to plasma kinetic and thermal energies. There are various kinds of simulations used to simulation reconnection; our work begins with Athena++, a magnetohydrodynamic (MHD) simulation code typically used for astrophysical problems, and a reconnection specific code file. The resistive MHD equations are solved with Riemann solvers. There was an initial test run without modifying the code to understand the dynamics of the simulation. We expand on the original reconnection problem file by implementing a radiative cooling term specific to the corona. The radiative cooling is theorized to have an effect on solar coronal plasma and magnetic reconnection dynamics. The cooling term will be tested with various parameters and compared to the case without cooling to study these dynamics. The condensation found in only the with cooling case emphasizes the importance of implementing this feature and will be later tested with a thermal conduction term. We want to determine the parameter regime where non-equilibrium cooling will be important for the reconnection dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Identification of the weak-to-strong transition in Alfvénic turbulence from space plasma

Plasma turbulence is a ubiquitous dynamical process that transfers energy across many spatial and temporal scales in astrophysical and space plasma systems. Although the theory of anisotropic magnetohydrodynamic (MHD) turbulence has successfully described natural phenomena, its core prediction of an Alfvénic transition from weak to strong MHD turbulence when energy cascades from large to small scales has not been observationally confirmed. Here we report evidence for the Alfvénic weak-to-strong transition in small-amplitude, turbulent MHD fluctuations in Earth’s magnetosheath using data from four Cluster spacecraft. Our observations demonstrate the universal existence of strong turbulence accompanied by weak turbulent fluctuations on large scales. Moreover, we find that the nonlinear interactions of MHD turbulence are crucial to the energy cascade, as they broaden the cascade directions and fluctuating frequencies. The observed connection between weak and strong MHD turbulence systems may be present in many astrophysical environments, such as star formation, energetic particle transport, turbulent dynamos, and solar corona or solar wind heating.

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 OSTI2024

A hybrid-kinetic simulation tool for non-thermal warm x-ray z-pinch sources, with gas-puff and wire array exemplars

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Euler equations and the Sod shock tube problem

The Euler equations are a subset of the magnetohydrodynamic (MHD) equations in the infinitely collisional, unmagnetized limit. MHD modeling is central to many areas of plasma physics, ranging from low-temperature glow discharges to inertial confinement fusion. An important aspect of the Euler equations is their ability to describe states with discontinuities, such as shock waves. A standard benchmark test for numerical implementation of the Euler equations is the Sod shock tube. In this test, the system is initialized at rest with a pressure and density discontinuity, which results in a shock wave traveling into the low-pressure region and a rarefaction wave traveling into the high-pressure region. Starting with the presentation of the Euler equations, a numerical algorithm is presented here to solve these equations in one dimension. This is followed by an overview of the Sod shock tube problem that includes the precise initial setup and the analytic solution. Finally, the analytic solution is compared with results from numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

UNR Grant Proposal Contribution

The work of Carrier et. al. (citation needed for Matt’s effect of surface roughness on phase transitions/ETI paper) demonstrates that 1D and 2D resistive magnetohydrodynamic (MHD) simulations can reliably model exploding aluminum rods driven by megaampere currents. Close agreement between simulation results and photonic Doppler velocimetry data from the Mykonos electrothermal instability (METI-II) campaign builds confidence in predictive modelling capabilities for pulsed-power HED experiments. Furthermore, 2D MHD simulations show how machined features and micro-scale surface roughness can seed the electrothermal instability (ETI). Surface roughness was observed to reduce the time of melt by 19%, with ETI growth driving the enhanced heating of the rod surface.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Structure Formation Through Magnetohydrodynamical Instabilities in Primordial Disks

The shear flow instabilities under the presence of magnetic fields in the primordial disk can greatly facilitate the formation of density structures that serve as seeds prior to the onset of the gravitational Jeans instability. We evaluate the effects of the Parker, magnetorotational and kinematic dynamo instabilities by comparing the properties of these instabilities. We calculate the mass spectra of coagulated density structures by the above mechanism in the radial direction for an axisymmetric magnetohydrodynamic (MHD) torus equilibrium and power density profile models. Our local three-dimensional MHD simulation indicates that the coupling of the Parker and magnetorotational instability creates spiral arms and gas blobs in an accretion disk, reinforcing the theory and model. Such a mechanism for the early structure formation may be tested in a laboratory. The recent progress in experiments involving shear flows in rotating tokamak, field reversed configuration (FRC) and laser plasmas may become a key element to advance in nonlinear studies.

Noguchi, Koichi↗
DOE OSTI2022

Slow Shock Formation Upstream of Reconnecting Current Sheets

The formation, development, and impact of slow shocks in the upstream regions of reconnecting current layers are explored. Slow shocks have been documented in the upstream regions of magnetohydrodynamic (MHD) simulations of magnetic reconnection as well as in similar simulations with the kglobal kinetic macroscale simulation model. They are therefore a candidate mechanism for preheating the plasma that is injected into the current layers that facilitate magnetic energy release in solar flares. Of particular interest is their potential role in producing the hot thermal component of electrons in flares. During multi-island reconnection, the formation and merging of flux ropes in the reconnecting current layer drives plasma flows and pressure disturbances in the upstream region. These pressure disturbances steepen into slow shocks that propagate along the reconnecting component of the magnetic field and satisfy the expected Rankine–Hugoniot jump conditions. Plasma heating arises from both compression across the shock and the parallel electric field that develops to maintain charge neutrality in a kinetic system. Shocks are weaker at lower plasma β, where shock steepening is slow. While these upstream slow shocks are intrinsic to the dynamics of multi-island reconnection, their contribution to electron heating remains relatively minor compared with that from Fermi reflection and the parallel electric fields that bound the reconnection outflow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗