Search NASA⌕ Search

SEARCH · Search NASA

Results for “effective kinetics”

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 145 records · Page 8

Structures of strong shocks in low-density helium and neon gases

Strong shocks are essential components in many high-energy-density environments such as inertial confinement fusion implosions. However, the experimental measurements of the spatial structures of such shocks are sparse. In this paper, the soft x-ray emission of a shock front in a helium gas mixture (90% helium, 10% neon) and a pure neon gas was spatially resolved using an imaging spectrometer. We observe that the shock width in the helium mixture gas is about twice as large as in the pure neon gas. Moreover, they exhibit different precursor layers, where electron temperature greatly exceeds ion temperature, extending for more than ~350 µm with the helium gas mixture but less than 30 µm in the pure neon. At the shock front, calculations show that the electrons are strongly collisional with mean-free path two orders of magnitude shorter than the characteristic length of the shock. However, the helium ions can reach a kinetic regime as a consequence of their mean-free path being comparable to the scale of the shock. A radiation-hydrodynamic simulation demonstrates the impact of thermal conduction on the formation of the precursors with charge state, Z, playing a major role in heat flow and the precursor formation in both the helium mixture and the pure neon gases. Particle-in-cell simulations are also performed to study the ion kinetic effects on the formation of the observed precursors. A group of fast-streaming ions is observed leading the shock only in the helium gas mixture. Additionally, both effects explain the longer precursor layer in the helium shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Species Separation and Hydrogen Streaming upon Shock Release from Polystyrene under Inertial Confinement Fusion Conditions

Shock release from inertial confinement fusion (ICF) shells poses a great challenge to single-fluid hydrodynamic equations, especially for describing materials composed of different ion species. This has been evidenced by a recent experiment [Haberberger et al., Phys. Rev. Lett. 123, 235001 (2019)], in which low-density plasmas ( 10 19 to 10 20 cm – 3 ) are measured to move far ahead of what radiation-hydrodynamic simulations predict. To understand such experimental observations, in this work we have performed large-scale nonequilibrium molecular-dynamics simulations of shock release in polystyrene (CH) at experimental conditions. These simulations revealed that upon shock releasing from the back surface of a CH foil, hydrogen can stream out of the bulk of the foil due to its mass being lighter than carbon. This released hydrogen, exhibiting a much broader velocity distribution than carbon, forms low-density plasmas moving in nearly constant velocities ahead of the in-flight shell, which is in quantitative agreement with the experimental measurements. Such kinetic effect of species separation is currently missing in single-fluid radiation-hydrodynamics codes for ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Equilibrium solute segregation to matrix- θ' precipitate interfaces in Al-Cu alloys from first principles

Particular combinations of solute atoms segregated to the interface of the Al and θ ' - Al 2 Cu in Al-Cu alloys can help stabilize θ ' precipitates at high temperatures. Stabilization of such precipitates is determined by a combination of thermodynamics (including driving forces for coarsening and transformation and solute segregation tendencies) and kinetic effects (including solute diffusion and interfacial mobility in the presence of interfacial solute segregation). For some alloys such as recent Al-Cu-based alloys, multiple solutes segregate in significant quantities to interfaces, and solute-solute interactions at the interface are important, with multiple types of solutes competing for similar interfacial sites. To treat this situation, we develop and apply a statistical mechanics approach to calculate the temperature-dependent equilibrium solute atoms distribution near the coherent and semicoherent interfaces between the Al matrix and the θ ' - Al 2 Cu precipitates. The developed approach is applied to the investigation of Si, Mn, and Zr segregation at the interface, as particular combinations of these elements affect the thermal stability of the θ ' precipitates. We demonstrate that because Si and Mn atoms segregate on the same semicoherent interface, the presence of Si reduces the concentration of Mn solutes at the interface. Si atoms preferably occupy the first layer of the interface and compete with Mn atoms for one type of particular sites in the layer. Mn atoms preferably occupy the second layer of the semicoherent interface, and the Mn-Mn interaction plays an important role in their distribution. Zr atoms mostly segregate on one of the two nonequivalent sites of the second layer of the coherent interface. Due to symmetry properties of the coherent interface, the calculations show that the segregation Zr of atoms to this interface will likely lead to the formation of L 1 2 ordered Al 3 Zr layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Phase transition kinetics revealed by in situ x-ray diffraction in laser-heated dynamic diamond anvil cells

We report successful coupling of dynamic loading in a diamond anvil cell and stable laser heating, which enables compression rates up to 500 GPa/s along high-temperature isotherms. Dynamic loading in a diamond-anvil cell allows exploration of a wider range of pathways in the pressure-temperature space compared to conventional dynamic compression techniques. By x-ray diffraction, we are able to characterize and monitor the structural transitions with the appropriate time resolution i.e., millisecond timescales. Using this method, we investigate the γ − ε phase transition of iron under dynamic compression, reaching compression rates of hundreds of GPa/s and temperatures of 2000 K. Our results demonstrate a distinct response of the γ − ε and α − ε transitions to the high compression rates achieved, possibly due to the different transition mechanisms. These findings open up new avenues to study tailored dynamic compression pathways in the pressure-temperature space and highlight the potential of this platform to capture kinetic effects (over ms time scales) in a diamond anvil cell. Published by the American Physical Society 2024

Ricks, Matthew (ORCID:000000015315612X)↗

Particle-in-cell simulations of burning inertial confinement fusion capsule implosions

Anomalies observed in the neutron spectral shift of high-yield shots at the National Ignition Facility (NIF) suggest the presence of suprathermal ions [E. P. Hartouni et al., Nat. Phys. 19, 72 (2023)], implying that kinetic effects play a significant role in burning inertial confinement fusion (ICF) plasmas. Furthermore, recent measurements of reaction-in-flight (RIF) neutrons offer a direct probe of the stopping power in the burning fuel region of high-energy alpha particles and up-scattered fuel ions. We have developed the particle-in-cell code PICNIC, an exactly energy-conserving particle-in-cell Monte-Carlo collision (PIC-MCC) code to simulate the burn stage in ICF. We present results from 1D spherical simulations of NIF shot N210808, which was the first to exceed the Lawson criterion for laser fusion. We find that the suprathermal ions generated by large-angle Rutherford and nuclear elastic scattering (NES) with fusion alphas produce an alpha knock-on neutron (AKN) signal consistent with the extent and relative yield of the AKN spectrum identified in ignition experiments at the NIF. Furthermore, we also find that the inclusion of large-angle scattering physics does not explain the anomalously large spectral shift observed in the experiment.

High-energy-density plasmas↗

Chemical vapor infiltration of additively manufactured preforms: Pore-resolved simulations and experimental validation

The densification of additively manufactured porous preforms by chemical vapor infiltration (CVI) is studied using pore-resolved simulations and experiments. Experimentally, 3D printed silicon carbide (SiC) preforms are subject to CVI synthesis using methyltrichlorosilane (MTS) precursor to obtain high purity SiC/SiC composites. Optical images of the cross sections of the processed preforms are analyzed to obtain the spatial porosity distribution. The numerical method is based on a level set formulation to capture the spatial distribution and time evolution of the pore scale microstructural characteristics. The coupled transport and kinetic effects are represented using a dimensionless Thiele modulus. Simulations are initialized using representative synthetic preform geometries comprising of packed particles based on the size distribution of the powder used for 3D printing. The simulation results are validated against the experimental observations in terms of total density and the distribution of residual porosity. The densification characteristics, porosity classification, concentration profiles, and structure functions are analyzed as functions of processing temperature and Thiele modulus.

36 MATERIALS SCIENCE↗

Supramolecular bending and twisting in the hierarchical self-assembly of monodisperse mesogenic oligomers

Understanding how different forms of supramolecular curvature arise during assembly is crucial to designing and tuning the microstructure of hierarchically self-assembled materials. Here, we show that in crystalline phases of mesogenic oligomers, the oligomer length is a critical parameter that determines the type of curvature (Gaussian or cylindrical) exhibited by the self-assembled structures. We use iterative exponential growth to synthesize monodisperse mesogenic oligomers ranging from dimer to octamer. By analyzing their phase behavior and microstructure, we elucidate how length-dependent thermodynamic and kinetic effects tune their hierarchical degree of ordering. The oligomers’ length-dependent crystalline order drives the formation of scrolled sheets in shorter oligomers and twisted ribbons in longer oligomers. These studies highlight how oligomer length interplays with mesogen geometry and crystalline packing to drive self-assembly, introducing oligomer length as a powerful design parameter toward tailored applications of mesogenic systems.

36 MATERIALS SCIENCE↗

Biochemical properties of glycerol kinase from the hypersaline-adapted archaeon Haloferax volcanii

ABSTRACT Extremophilic microorganisms are promising candidates for industrial and analytical biocatalysis.Haloferax volcanii, a halophilic archaeon that prefers glycerol over glucose, channels this substrate into central metabolism through glycerol kinase (GK). Here, we report the biochemical properties ofH. volcaniiGK and its potential for biotechnological applications. An N-terminal His-tagged GK was functionalin vivoand yielded 3 mg/L culture—4.5 times more enzyme than a C-terminal StrepII-tagged version. Size exclusion chromatography revealed a glycerol-induced oligomeric shift from homodimer to a dimer-dominant state with detectable tetramer. The purified enzyme showed robust activity across broad pH and salinity ranges, with optimal activity at 100 mM NaCl and 50°C–60°C. It retained catalytic activity in 5%–10% dimethyl sulfoxide (DMSO) and crude glycerol containing methanol. His-GK was freeze-thaw stable and thermotolerant in 2 M NaCl buffers. In the absence of ligands, the enzyme’s melting temperature (T m ) was 80°C. Glycerol increased the T m to 85°C, and combinations with MgCl₂ (84°C) or ATP (88°C) provided further stabilization. The highest T m (89°C) occurred with all three ligands, suggesting a cumulative stabilizing effect. Kinetic analyses revealed positive cooperativity for glycerol, ATP, and Mg² + ; Mn² + and Co² + also supported the activity.H. volcaniiGK is the first known GK to exhibit positive cooperativity with glycerol and ATP. Its high stability and substrate flexibility support its use in biodiesel waste valorization,in vitrobiocatalysis, and biosensor development—applications demanding robust, specific, and stable enzymes. IMPORTANCE This study reveals thatH. volcaniiGK exhibits positive cooperativity for glycerol, ATP, and Mg² + , a kinetic feature not previously reported for glycerol kinases. This behavior enables steep, switch-like responses to small substrate changes, offering unique advantages for biosensor design. Importantly,H. volcaniiGK also maintains high activity under extreme salinity, temperature, broad pH, and solvent conditions that typically limit enzyme use in industrial and environmental applications. These traits make this GK an ideal candidate for enzyme-based biosensors, which often suffer from poor tolerance to pH, solvent, and thermal stress. Its robustness supports its use in cross-linked enzyme crystals, an immobilization method that enhances enzyme stability and reusability under harsh conditions. Moreover, GKs are already employed in Mg² + detection kits; however,H. volcaniiGK’s ability to tolerate and respond to diverse divalent cations (e.g., Co² + , Mn² + ) broadens their potential for pollutant detection and environmental monitoring. These features collectively positionH. volcaniiGK as a valuable biocatalyst for biosensing,in vitrodiagnostics, and biotechnological applications requiring both precision and durability.

Biotechnology & Applied Microbiology↗

A Thermodynamic Reassessment of Lithium-Ion Battery Cathode Calorimetry

This work demonstrates how staged heat release from layered metal oxide cathodes in the presence of organic electrolytes can be predicted from basic thermodynamic properties. These prediction methods for heat release are an advancement compared to typical modeling approaches for thermal runaway in lithium-ion batteries, which tend to rely exclusively on calorimetry measurements of battery components. These calculations generate useful new insights when compared to calorimetry measurements for lithium cobalt oxide (LCO) as well as the most common varieties of nickel manganese cobalt oxide (NMC) and nickel cobalt aluminum oxide (NCA). Accurate trends in heat release with varying state of charge are predicted for all of these cathode materials. These results suggest that thermodynamic calculations utilizing a recently published database of properties are broadly applicable for predicting decomposition behavior of layered metal oxide cathodes. Aspects of literature calorimetry measurements relevant to thermal runaway modeling are identified and classified as thermodynamic or kinetic effects. The calorimetry measurements reviewed in this work will be useful for development of a new generation of thermal runaway models targeting applications where accurate maximum cell temperatures are required to predict cascading cell-to-cell propagation rates.

25 ENERGY STORAGE↗

Research three-dimensional hierarchical graphite architectures for anodes for fast charging

With current lithium ion batteries optimized for performance under relatively low charge rate conditions, implementation of XFC has been hindered by drawbacks including Li plating, kinetic polarization, and heat dissipation. This project will utilize model-informed design of 3-D hierarchical electrodes to tune key XFCrelated variables like 1) bulk porosity/tortuosity 2) vertical pore diameter, spacing, and lattice 3) crystallographic orientation of graphite particles relative to exposed surfaces 4) interfacial chemistry of the graphite surfaces through "artificial sEr' formation using ALD 5) current collector surface roughness (aspect ratio, roughness factor, etc.). A key aspect of implementing novel electrodes is characterizing them in relevant settings. For this project, ultimately led out of University of Michigan by Neil Dasgupta, that includes both coin cell and 2+ Ah pouch cell testing, as well as comparison testing against baselines. Sandia National Labs will be conducting detailed cell characterization on iterative versions/improvements of the model-based hierarchical electrodes, as well as COTS cells for baseline comparisons. Key metrics include performance under fast charge conditions, as well as the absence or degree of lithium plating. Sandia will use their unique high precision cycling and rapid EIS capabilities to accurately characterize performance and any lithium plating during 6C charging and beyond, coupling electrochemical observations with cell teardown. Sandia will also design custom fixturing to cool cells during rapid charge, to decouple any kinetic effects brought about by cell heating and allow comparisons between different cells and charge rates. Using these techniques, Sandia will assess HOH electrodes from the University of Michigan, as well as aiding in iterative model and electrode design.

25 ENERGY STORAGE↗

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Second Year Report/ Final Report)

This research will focus on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and ow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal ow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aim to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Final Report)

This research has focused on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and flow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal flow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aimed to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CRADA Final Report: CRADA Number NFE-22-09330 with General Fusion

General Fusion is developing a magnetized target fusion (MTF) approach that involves compressing an initial magnetically confined plasma inside a cavity formed in liquid metal. This approach builds from concepts initially developed under the Linus program at the U.S. Naval Research Laboratory and combines it with advances from compact toroid experiment (CTX) and sustained spheromak plasma experiment (SSPX) in compact toroid plasmas and coaxial Marshall gun systems. Modeling the tokamak during compression is central to designing a successful MTF device. The plasma is formed by coaxial helicity injection in the General Fusion device. Immediately after formation, the plasma has a diverted tokamak configuration with a single null. As the wall moves inwards, the plasma is repelled from the conducting surface and driven inwards by currents induced by its magnetic field in the liquid metal wall. As the liquid metal closes (or bridges) the opening of the coaxial plasma injector, the magnetic field topology alters to remove the null. Due to this, the plasma moves from a diverted to a wall-limited configuration. The liquid metal liner continues to close in and change shape, reducing in radius by a factor of ten at the peak of plasma compression. A model of the MTF plasma must be able to handle this continually varying geometry, and to be predictive, it must faithfully include the real imperfections arising in the process. In this project, we pursued a Monte Carlo approach to closures for MHD by computing kinetic electron trajectories in an MHD plasma background from simulations of GF devices. This requires enhancing the capabilities of the KORC-T code for running large ensembles of kinetic trajectories by porting it to GPU architectures and enabling workflows for large ensembles on OLCF machines. With these capabilities, it is possible to produce a large library of kinetic calculations of electron orbits evolving in plasma configurations spanning the magnetic configurations and plasma density profiles, including non-axisymmetry, arising in the General Fusion’s existing PI3 spherical tokamak device. Using ensembles will capture particles passing a single point in space in a given magnetic configuration, and the entire dataset will cover a range of global magnetic field geometries. By sampling around many starting points, this dataset will capture the spatial dependence of the plasma parameters. From this large dataset, it is possible to produce a reduced model for the kinetic effects not captured in MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Addressing key physics problems in high-energy-density plasmas with a novel kinetic simulation capability

Many important physical processes in inertial confinement fusion (ICF) and dense Z-pinch (DZP) experiments require a kinetic (velocity-space-dependent) description. Conventional particle-in-cell (PIC) methods are poorly suited for high-energy-density (HED) plasmas, due to restrictive time-step constraints and the inability to conserve energy. In a previous LDRD (21-FS-048), we demonstrated that a fully implicit PIC formulation overcomes these limitations: it conserves energy even when coupled with Coulomb collision models and can be solved efficiently with large grid cells and large time steps. Thus, it is feasible to use this method to study kinetic effects in ICF and DZP plasmas on hydro-like time and spatial scales. In this follow-on LDRD, we advanced this methodology into a high-fidelity tool for production-scale simulations and used it to answer key questions relevant to ICF and DZP experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FullWave — A Full Wave Parallel Code for Modeling RF Fields in Hot Tokamak Plasma

FullWave is a computer code that simulates how radio-frequency (RF) waves travel and deposit energy in the hot plasma inside a fusion reactor. RF waves are used to heat the plasma and drive electrical current, which is essential for sustaining fusion reactions. The code uses a new algorithm that can handle much finer spatial detail than previous codes — more than 100 times finer — while running efficiently on national supercomputers. It incorporates a detailed physics model that captures subtle kinetic effects important for accurate prediction of wave behavior. Under this project, FullWave was extended to cover multiple RF frequency ranges relevant to present and future tokamaks, and validated against experimental parameters from the DIII-D tokamak at General Atomics. Results were published in peer-reviewed journal articles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding the Impact of Hydrogen Activation by SrCe 0.8 Zr 0.2 O 3-δ Perovskite Membrane Material on Direct Non-Oxidative Methane Conversion

Direct non-oxidative methane conversion (DNMC) converts methane (CH 4 ) in one step to olefin and aromatic hydrocarbons and hydrogen (H 2 ) co-product. Membrane reactors comprising methane activation catalysts and H 2 -permeable membranes can enhance methane conversion by in situ H 2 removal via Le Chatelier's principle. Rigorous description of H 2 kinetic effects on both membrane and catalyst materials in the membrane reactor, however, has been rarely studied. In this work, we report the impact of hydrogen activation by hydrogen-permeable SrCe 0.8 Zr 0.2 O 3-δ (SCZO) perovskite oxide material on DNMC over an iron/silica catalyst. The SCZO oxide has mixed ionic and electronic conductivity and is capable of H 2 activation into protons and electrons for H 2 permeation. In the fixed-bed reactor packed with a mixture of SCZO oxide and iron/silica catalyst, stable and high methane conversion and low coke selectivity in DNMC was achieved by co-feeding of H 2 in methane stream. The characterizations show that SCZO activates H 2 to favor “soft coke” formation on the catalyst. The SCZO could absorb H 2 in situ to lower its local concentration to mitigate the reverse reaction of DNMC in the tested conditions. The co-existence of H 2 co-feed, SCZO oxide, and DNMC catalyst in the present study mimics the conditions of DNMC in the H 2 -permeable SCZO membrane reactor. The findings in this work offer the mechanistic understanding of and guidance for the design of H 2 -permeable membrane reactors for DNMC and other alkane dehydrogenation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laboratory Study of Magnetic Reconnection in Lunar-relevant Mini-magnetospheres

Abstract Mini-magnetospheres are small ion-scale structures that are well suited to studying kinetic-scale physics of collisionless space plasmas. Such ion-scale magnetospheres can be found on local regions of the Moon, associated with the lunar crustal magnetic field. In this paper, we report on the laboratory experimental study of magnetic reconnection in laser-driven, lunar-like ion-scale magnetospheres on the Large Plasma Device at the University of California, Los Angeles. In the experiment, a high-repetition rate (1 Hz), nanosecond laser is used to drive a fast-moving, collisionless plasma that expands into the field generated by a pulsed magnetic dipole embedded into a background plasma and magnetic field. The high-repetition rate enables the acquisition of time-resolved volumetric data of the magnetic and electric fields to characterize magnetic reconnection and calculate the reconnection rate. We notably observe the formation of Hall fields associated with reconnection. Particle-in-cell simulations reproducing the experimental results were performed to study the microphysics of the interaction. By analyzing the generalized Ohm’s law terms, we find that the electron-only reconnection is driven by kinetic effects through the electron pressure anisotropy. These results are compared to recent satellite measurements that found evidence of magnetic reconnection near the lunar surface.

Astronomy & Astrophysics↗

First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole

In this paper we provide a first physical interpretation for the Event Horizon Telescope's (EHT) 2017 observations of Sgr A*. Our main approach is to compare resolved EHT data at 230 GHz and unresolved non-EHT observations from radio to X-ray wavelengths to predictions from a library of models based on time-dependent general relativistic magnetohydrodynamics simulations, including aligned, tilted, and stellar-wind-fed simulations; radiative transfer is performed assuming both thermal and nonthermal electron distribution functions. We test the models against 11 constraints drawn from EHT 230 GHz data and observations at 86 GHz, 2.2μm, and in the X-ray. All models fail at least one constraint. Light-curve variability provides a particularly severe constraint, failing nearly all strongly magnetized (magnetically arrested disk (MAD)) models and a large fraction of weakly magnetized models. A number of models fail only the variability constraints. We identify a promising cluster of these models, which are MAD and have inclination i ≤ 30°. They have accretion rate (5.2–9.5) × 10$^{−9}$ M $_{⊙}$ yr$^{−1}$, bolometric luminosity (6.8–9.2) × 10$^{35}$ erg s$^{−1}$, and outflow power (1.3–4.8) × 10$^{38}$ erg s$^{−1}$. We also find that all models with i ≥ 70° fail at least two constraints, as do all models with equal ion and electron temperature; exploratory, nonthermal model sets tend to have higher 2.2μm flux density; and the population of cold electrons is limited by X-ray constraints due to the risk of bremsstrahlung overproduction. Finally, we discuss physical and numerical limitations of the models, highlighting the possible importance of kinetic effects and duration of the simulations.

79 ASTRONOMY AND ASTROPHYSICS↗