Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle-in-cell method”

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 OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗
NASA NTRS2010

Particle-in-Cell Simulation of Collisionless Driven Reconnection with Open Boundaries

First results are discussed from an ongoing study of driven collisionless reconnection using a 2 1/2-dimensional electromagnetic particle-in-cell simulation model with open inflow and outflow boundaries. An extended electron diffusion region (EEDR) is defined as that region surrounding a reconnecting neutral line in which the out-of-plane nonideal electric field is positive. It is shown that the boundaries of this region in the directions of the outflow jets are at the positions where the electrons make the transition from unfrozen meandering motion in the current sheet to outward drifting with the magnetic field in the outflow jets; a turning length scale is defined to mark these positions, The initial width of the EEDR in the inflow directions is comparable to the electron bounce width. Later. as shoulders develop to form a two-scale structure. thc EEDR width expands to the ion bounce width scale. The inner portion of the EEDR or the electron diffusion region proper remains at the electron bounce width. Two methods are introduced for predicting the reconnection electric field using the dimensions of the EEDR. These results are interpreted as further evidence that the EEDR is the region that is relevant to understanding the electron role in the neutral line vicinity.

Kimas, Alex↗
DOE OSTI2025

Kinetic Plasma Simulation Capabilities in the MOOSE Framework: Verification of Particle-Particle Collisions

High-fidelity simulations of complex plasma systems allow researchers to gain key insights into and understanding of these systems. To facilitate massively parallel high-fidelity plasma simulations, finite-element-based particle-in-cell capabilities are being developed within the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) based framework called Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER). While SALAMANDER’s primary objective is modeling edge plasmas and plasma-facing components in fusion devices, the particle-in-cell capabilities being developed are general and will support modeling low-temperature plasmas as well. Previously, collisionless magnetostatic simulation capabilities have been verified with the two-stream and Dorey-Guest-Harris instabilities, and single particle motion. Collisions were implemented using the direct simulation Monte Carlo method, and verification of this capability will be presented here several verification problems: relaxation of a randomly initialized gas to a Maxwellian distribution, Fourier heat flow, and comparison of reaction rates to both analytic calculations and those calculated using a multi-term Boltzmann solver.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Critical magnetic Reynolds number of the turbulent dynamo in collisionless plasmas

The intracluster medium of galaxy clusters is an extremely hot and diffuse, nearly collisionless plasma, which hosts dynamically important magnetic fields of ~μG strength. Seed magnetic fields of much weaker strength of astrophysical or primordial origin can be present in the intracluster medium. In collisional plasmas, which can be approximated in the magnetohydrodynamical (MHD) limit, the turbulent dynamo mechanism can amplify weak seed fields to strong dynamical levels efficiently by converting turbulent kinetic energy into magnetic energy. However, the viability of this mechanism in weakly collisional or completely collisionless plasma is much less understood. In this study, we explore the properties of the collisionless turbulent dynamo using three-dimensional hybrid-kinetic particle-in-cell simulations. We explore the properties of the collisionless turbulent dynamo in the kinematic regime for different values of the magnetic Reynolds number, Rm, initial magnetic-to-kinetic energy ratio, (E mag /E kin ) i , and initial Larmor ratio, (r Larmor /L box ) i , i.e. the ratio of the Larmor radius to the size of the turbulent system. We find that in the ‘un-magnetized’ regime, (r Larmor /L box ) i > 1, the critical magnetic Reynolds number for the dynamo action Rm crit ≈ 107 ± 3. In the ‘magnetized’ regime, (r Larmor /L box ) i ≲ 1, we find a marginally higher Rm crit = 124 ± 8. We find that the growth rate of the magnetic energy does not depend on the strength of the seed magnetic field when the initial magnetization is fixed. We also study the distribution and evolution of the pressure anisotropy in the collisionless plasma and compare our results with the MHD turbulent dynamo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Explicit dispersion relations for warm fluid waves in a uniform plasma (invited)

Classical dispersion relations for waves in a fluid plasma are expressed as implicit functions of wave frequency and wave- number. The implicit dispersion relation must in general be solved numerically. This work introduces an astute method for obtaining an explicit dispersion relation for waves in a fluid plasma. The explicit expression has an advantage of providing eigenmodes without numerical computations and giving a dispersion relation more detailed than the implicit expression. As in the case of cold waves, the wavenumber for an arbitrary frequency can be directly obtained from the explicit formula. The explicit formula also enables an accurate evaluation of finite-temperature effects on a dispersion relation even near the cyclotron resonance, where a numerical analysis of the implicit relation is impractical. As a result, the analytic formula can be used to investigate temperature effects on electromagnetic waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

EZ: An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations

Here, we present EZ, a novel current deposition algorithm for particle-in-cell (PIC) simulations. EZ calculates the current density on the electromagnetic grid due to macro-particle motion within a time step by solving the continuity equation of electrodynamics. Being a charge conserving hybridization of Esirkepov's method and ZigZag, we refer to it as “EZ” as shorthand for “Esirkepov meets ZigZag”. Simulations of a warm, relativistic plasma with PIConGPU show that EZ achieves the same level of charge conservation as the commonly used method by Esirkepov, yet reaches higher performance for macro-particle assignment-functions up to third-order. In addition to a detailed description of the functioning of EZ, reasons for the expected and observed performance increase are given, and guidelines for its implementation aiming at highest performance on GPUs are provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Time-dependent saturation and physics-based nonlinear model of cross-beam energy transfer

The nonlinear physics of cross-beam energy transfer (CBET) for multi-speckled laser beams is examined using large-scale particle-in-cell simulations for a range of laser and plasma conditions relevant to indirect-drive inertial confinement fusion (ICF) experiments. The time-dependent growth and saturation of CBET involve complex, nonlinear ion and electron dynamics, including ion trapping-induced enhancement and detuning, ion acoustic wave (IAW) nonlinearity, oblique forward stimulated Raman scattering (FSRS), and backward stimulated Brillouin scattering (BSBS) in a CBET-amplified seed beam. Ion-trapping-induced detuning of CBET is captured in the kinetic linear response by a new δf-Gaussian-mixture algorithm, enabling an accurate characterization of trapping-induced non-Maxwellian distributions. Ion trapping induces nonlinear processes, such as changes to the IAW dispersion and nonlinearities (e.g., bowing and self-focusing), which, together with pump depletion, FSRS, and BSBS, determine the time-dependent nature and level of CBET gain as the system approaches a steady state. Using VPIC simulations at intensities at and above the onset threshold for ion trapping and the insight from the time-dependent saturation analyses, we construct a nonlinear CBET model from local laser and plasma conditions that predicts the CBET gain and the energy deposition into the plasma. This model is intended to provide a more accurate, physics-based description of CBET saturation over a wide range of conditions encountered in ICF hohlraums compared with linear CBET gain models with ad hoc saturation clamps often used in laser ray-based methods in multi-physics codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An implicit particle code with exact energy and charge conservation for electromagnetic studies of dense plasmas

A collisional particle code based on implicit energy- and charge-conserving methods is presented. A modified version of the particle-suppressed Jacobian-Free Newton-Krylov method that can enhance the solver efficiency is introduced. Mathematically, it is shown that this new approach can be viewed as a fixed-point iteration method for the particle positions. The model can exactly conserve global energy and local charge and can efficiently use time steps larger than the plasma period. In conclusion, the algorithm's ability to simulate dense plasmas accurately and efficiently is quantified by simulating the dynamic compression of a plasma slab via a magnetic piston in 1D planar geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Isolation and phase-space energization analysis of the instabilities in collisionless shocks

We analyse the generation of kinetic instabilities and their effect on the energization of ions in non-relativistic, oblique collisionless shocks using a 3D-3V (three spatial with three velocity components) simulation by dHybridR, a hybrid particle-in-cell code. At sufficiently high Mach number, quasi-perpendicular and oblique shocks can experience rippling of the shock surface caused by kinetic instabilities arising from free energy in the ion velocity distribution due to the combination of the incoming ion beam and the population of ions reflected at the shock front. To understand the role of the ripple on particle energization, we devise a new instability isolation method to identify the unstable modes underlying the ripple and interpret the results in terms of the governing kinetic instability. We generate velocity-space signatures using the field–particle correlation technique to look at energy transfer in phase space from the isolated instability driving the shock ripple, providing a viewpoint on the different dynamics of distinct populations of ions in phase space. Together, the field–particle correlation technique and our new instability isolation method provide a unique viewpoint on the different dynamics of distinct populations of ions in phase space and allow us to completely characterize the energetics of the collisionless shock under investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effect of magnetic field inclination on black hole jet power and particle acceleration

Rotating black holes are known to launch relativistic jets and accelerate particles if they accrete a magnetized plasma. It remains unclear, however, how the global magnetic field orientation affects the jet powering efficiency. We propose the first kinetic study of a collisionless plasma around a Kerr black hole embedded in a magnetic field that is inclined with respect to the black hole spin axis. Using three-dimensional general relativistic particle-in-cell simulations, we show that while oblique magnetic field configurations significantly reduce the jet power, particle acceleration still remains highly efficient. This suggests that black holes producing a weak jet might still be bright sources of nonthermal radiation and cosmic rays.

acceleration of particles↗
DOE OSTI2026

Demonstration of a Novel Phase-Space Painting Method in a Coupled Lattice to Mitigate Space Charge in High-Intensity Hadron Beams

Multiturn charge-exchange injection is the primary method of creating high-intensity hadron beams in circular accelerators, and phase space painting during injection enables tailoring of the accumulated phase space distribution. A technique we call eigenpainting allows injection of particles into a single mode of a coupled ring, providing full four-dimensional control of the phase space distribution. Under ideal conditions, uniform eigenpainting generates a linear-force equilibrium distribution in the transverse plane, with zero volume in four-dimensional transverse phase space, even including space charge. Here, we have implemented eigenpainting for the first time in the spallation neutron source accumulator ring. Injecting 8.8 μ⁢C of an 800 MeV beam, we obtain a final ratio of intrinsic transverse emittances of ≈2.4. We analyze the effect of space charge on the final distribution through comparison of the reconstructed phase space to particle-in-cell simulations.

Evans, Nicholas J. [Oak Ridge National Laboratory ↗
DOE OSTI2022

Absorption of charged particles in perfectly matched layers by optimal damping of the deposited current

Perfectly matched layers (PMLs) are widely used in particle-in-cell simulations, in order to absorb electromagnetic waves that propagate out of the simulation domain. However, when charged particles cross the interface between the simulation domain and the PMLs, a number of numerical artifacts can arise. Here, in order to mitigate these artifacts, we introduce a PML algorithm whereby the current deposited by the macroparticles in the PML is damped by an analytically derived optimal coefficient. The benefits of this algorithm are illustrated in practical simulations. In particular, it is shown that this algorithm is well suited for particles exiting the box in near-normal incidence, in the sense that the fields behave as if the exiting particle is propagating in an infinite vacuum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Full 3D Kinetic Modeling and Quantification of Positive Streamer Evolution in an Azimuthally Swept Pin-to-Plane Wedge Geometry

Cathode-directed streamer evolution in near atmospheric air is modeled in 3D pin-to-plane geometries using a 3D kinetic Particle-In-Cell (PIC) code that simulates particle-particle collisions via the Direct Simulation Monte Carlo (DSMC) method. Due to the computational challenges associated with a complete 360° volumetric domain, a practical alternative was achieved using a wedge domain and a range of azimuthal angles was explored (5°, 15°, 30°, and 45°) to study possible effects on the streamer growth and propagation due to the finite wedge angle. A DC voltage of 6 kV is administered to a hemispherical anode of radius 100 μm, with a planar cathode held at ground potential, generating an over-volted state with an electric field of 4 MV/m across a 1500 μm gap. The domain is seeded with an initial ion and electron density of 10 18 m -3 at 1 eV temperature confined to a spherical region of radius 100 μm centered at the tip of the anode. The air chemistry model [1] includes standard Townsend breakdown mechanisms (electron-neutral elastic, excitation, ionization, attachment, and detachment collision chemistry and secondary electron emission) as well as streamer mechanisms (photoionization and ion-neutral collisions) via tracking excited state neutrals which can then either quench via collisions or spontaneously emit a photon based on specific Einstein-A coefficients [2, 3]. In this work, positive streamer dynamics are formally quantified for each wedge angle in terms of electron velocity and density as temporal functions of coordinates r, Φ, and z. Applying a random plasma seed for each simulation, particles of interest are tracked with near femtosecond temporal resolution out to 1.4 ns and spatially binned. This process is repeated six times and results are averaged. Prior 2D studies have shown that the reduced electric field, E/n, can significantly impact streamer evolution [4]. We extend the analysis to 3D wedge geometries, to limit computational costs, and examine the wedge angle’s effect on streamer branching, propagation, and velocity. Results indicate that the smallest wedge angle that produced an acceptably converged solution is 30°. The potential effects that a mesh, when under-resolved with respect to the Debye length, can impart on streamer dynamics and numerical heating were not investigated, and we explicitly state here that the smallest cell size was approximately 10 times the minimum λ D in the streamer channel at late times. This constraint on cell size was the result of computational limitations on total mesh count.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Fluid-Kinetic Coupling: Advanced Discretizations for Simulations on Emerging Heterogeneous Architectures (LDRD FY20-0643)

Plasma physics simulations are vital for a host of Sandia mission concerns, for fundamental science, and for clean energy in the form of fusion power. Sandia's most mature plasma physics simulation capabilities come in the form of particle-in-cell (PIC) models and magnetohydrodynamics (MHD) models. MHD models for a plasma work well in denser plasma regimes when there is enough material that the plasma approximates a fluid. PIC models, on the other hand, work well in lower-density regimes, in which there is not too much to simulate; error in PIC scales as the square root of the number of particles, making high-accuracy simulations expensive. Real-world applications, however, almost always involve a transition region between the high-density regimes where MHD is appropriate, and the low-density regimes for PIC. In such a transition region, a direct discretization of Vlasov is appropriate. Such discretizations come with their own computational costs, however; the phase-space mesh for Vlasov can involve up to six dimensions (seven if time is included), and to apply appropriate homogeneous boundary conditions in velocity space requires meshing a substantial padding region to ensure that the distribution remains sufficiently close to zero at the velocity boundaries. Moreover, for collisional plasmas, the right-hand side of the Vlasov equation is a collision operator, which is non-local in velocity space, and which may dominate the cost of the Vlasov solver. The present LDRD project endeavors to develop modern, foundational tools for the development of continuum-kinetic Vlasov solvers, using the discontinuous Petrov-Galerkin (DPG) methodology, for discretization of Vlasov, and machine-learning (ML) models to enable efficient evaluation of collision operators. DPG affords several key advantages. First, it has a built-in, robust error indicator, allowing us to adapt the mesh in a very natural way, enabling a coarse velocity-space mesh near the homogeneous boundaries, and a fine mesh where the solution has fine features. Second, it is an inherently high-order, high-intensity method, requiring extra local computations to determine so-called optimal test functions, which makes it particularly suited to modern hardware in which floating-point throughput is increasing at a faster rate than memory bandwidth. Finally, DPG is a residual-minimizing method, which enables high-accuracy computation: in typical cases, the method delivers something very close to the $L^2$ projection of the exact solution. Meanwhile, the ML-based collision model we adopt affords a cost structure that scales as the square root of a standard direct evaluation. Moreover, we design our model to conserve mass, momentum, and energy by construction, and our approach to training is highly flexible, in that it can incorporate not only synthetic data from direct-simulation Monte Carlo (DSMC) codes, but also experimental data. We have developed two DPG formulations for Vlasov-Poisson: a time-marching, backward-Euler discretization and a space-time discretization. We have conducted a number of numerical experiments to verify the approach in a 1D1V setting. In this report, we detail these formulations and experiments. We also summarize some new theoretical results developed as part of this project (published as papers previously): some new analysis of DPG for the convection-reaction problem (of which the Vlasov equation is an instance), a new exponential integrator for DPG, and some numerical exploration of various DPG-based time-marching approaches to the heat equation. As part of this work, we have contributed extensively to the Camellia open-source library; we also describe the new capabilities and their usage. We have also developed a well-documented methodology for single-species collision operators, which we applied to argon and demonstrated with numerical experiments. We summarize those results here, as well as describing at a high level a design extending the methodology to multi-species operators. We have released a new open-source library, MLC, under a BSD license; we include a summary of its capabilities as well.

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 OSTI2026

Onset of Fluidization in MP-PIC Simulations using MFIX-Exa

Fluidized bed reactors are used across a variety of industries, including for energy processes like pyrolysis that result in low-cost energy products. Design and scale-up of fluidized beds is de-risked by modeling and simulation, utilizing tools like NETL’s MFIX-Exa High-Performance Computing (HPC) code for reacting multiphase flow. This report summarizes an investigation into the breadth of problems to which MFIX-Exa may be applied, specifically with regard to low fluid velocities and the onset of fluidization. A simple fluidization study is conducted both experimentally and numerically for particles of interest, then reactor simulations are compared to cold flow experiments for uniform distributor plates. Approaches for modeling bubble caps are also presented.

discrete particle method↗
DOE OSTI2025

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗
DOE OSTI2023

Exact Calculation of Nonideal Fields Demonstrates Their Dominance of Injection in Relativistic Reconnection

Magnetic reconnection is an important source of energetic particles in systems ranging from astrophysics to the laboratory. The large separation of spatiotemporal scales involved makes it critical to determine the minimum physical model containing the necessary physics for modeling particle acceleration. Here, by resolving the energy gain from ideal and nonideal magnetohydrodynamic electric fields self-consistently in kinetic particle-in-cell simulations of reconnection, we conclusively show the dominant role of the nonideal field for the early stage of energization known as injection. The importance of the nonideal field increases with magnetization, guide field, and in three dimensions, indicating its general importance for reconnection in natural astrophysical systems. We obtain the statistical properties of the injection process from the simulations, paving the way for the development of extended MHD models capable of accurately modeling particle acceleration in large-scale systems. The novel analysis method developed in this study can be applied broadly to give new insight into a wide range of processes in plasma physics.

79 ASTRONOMY AND ASTROPHYSICS↗