Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnetohydrodynamical simulations”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Reynolds number dependence of Lagrangian dispersion in direct numerical simulations of anisotropic magnetohydrodynamic turbulence

Large-scale magnetic fields thread through the electrically conducting matter of the interplanetary and interstellar medium, stellar interiors and other astrophysical plasmas, producing anisotropic flows with regions of high-Reynolds-number turbulence. It is common to encounter turbulent flows structured by a magnetic field with a strength approximately equal to the root-mean-square magnetic fluctuations. In this work, direct numerical simulations of anisotropic magnetohydrodynamic (MHD) turbulence influenced by such a magnetic field are conducted for a series of cases that have identical resolution, and increasing grid sizes up to $2048^3$ . The result is a series of closely comparable simulations at Reynolds numbers ranging from 1400 up to 21 000. We investigate the influence of the Reynolds number from the Lagrangian viewpoint by tracking fluid particles and calculating single-particle and two-particle statistics. The influence of Alfvénic fluctuations and the fundamental anisotropy on the MHD turbulence in these statistics is discussed. Single-particle diffusion curves exhibit mildly superdiffusive behaviours that differ in the direction aligned with the magnetic field and the direction perpendicular to it. Competing alignment processes affect the dispersion of particle pairs, in particular at the beginning of the inertial subrange of time scales. Scalings for relative dispersion, which become clearer in the inertial subrange for a larger Reynolds number, can be observed that are steeper than indicated by the Richardson prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

As a Matter of State: The Role of Thermodynamics in Magnetohydrodynamic Turbulence

Turbulence simulations play a key role in advancing the general understanding of the physical properties of turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence simulations are often conducted in the isothermal limit. Given that the majority of astrophysical systems are not governed by isothermal dynamics, we aim to quantify the impact of thermodynamics on the physics of turbulence, through varying adiabatic index, γ, combined with a range of optically thin cooling functions. Here, we present a suite of ideal magnetohydrodynamics simulations of thermally balanced stationary turbulence in the subsonic, super-Alfvénic, high ${\beta }_{{\rm{p}}}$ (ratio of thermal to magnetic pressure) regime, where turbulent dissipation is balanced by two idealized cooling functions (approximating linear cooling and free–free emission) and examine the impact of the equation of state by considering cases that correspond to isothermal, monatomic, and diatomic gases. We find a strong anticorrelation between thermal and magnetic pressure independent of thermodynamics, whereas the strong anticorrelation between density and magnetic field found in the isothermal case weakens with increasing γ. Similarly, the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing γ. This suggests that there exists a degeneracy in these relations with respect to thermodynamics and Mach number in this regime, which is dominated by slow magnetosonic modes. These results have implications for attempts to infer (e.g.,) Mach numbers from (e.g.,) Faraday rotation measurements, without additional information regarding the thermodynamics of the plasma. However, our results suggest that this degeneracy can be broken by utilizing higher-order moments of observable distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Hybrid simulation of energetic particles interacting with magnetohydrodynamics using a slow manifold algorithm and GPU acceleration

The hybrid method combining particle-in-cell and magnetohydrodynamics can be used to study the interaction between energetic particles and global plasma modes. In this paper we introduce the M3D-C1-K code, which is developed based on the M3D-C1 finite element code solving the magnetohydrodynamics equations, with a newly developed kinetic module simulating energetic particles. The particle pushing is done using a new algorithm by applying the Boris pusher to the classical Pauli particles to simulate the slow-manifold of particle orbits, with long-term accuracy and fidelity. The particle pushing can be accelerated using GPUs with a significant speedup. The moments of the particles are calculated using the δƒ method, and are coupled into the magnetohydrodynamics simulation through pressure or current coupling schemes. Several linear simulations of magnetohydrodynamics modes driven by energetic particles have been conducted using M3D-C1-K with the δƒ method, including fishbone, toroidal Alfvén eigenmodes and reversed shear Alfvén eigenmodes. Good agreement with previous results from other eigenvalue, kinetic and hybrid codes have been achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Numerical Simulation of Electron Magnetohydrodynamics with Landau-quantized Electrons in Magnetar Crusts

Abstract In magnetar crusts, magnetic fields are sufficiently strong to confine electrons into a small to moderate number of quantized Landau levels. This can have a dramatic effect on the crust's thermodynamic properties, generating field-dependent de Haas–van Alphen oscillations. We previously argued that the large-amplitude oscillations of the magnetic susceptibility could enhance the ohmic dissipation of the magnetic field by continuously generating small-scale, rapidly dissipating field features. This could be important to magnetar field evolution and contribute to their observed higher temperatures. To study this, we performed quasi-3D numerical simulations of electron MHD in a representative volume of neutron star crust matter, for the first time including the magnetization and magnetic susceptibility resulting from Landau quantization. We find that the potential enhancement in the ohmic dissipation rate due to this effect can be a factor ∼3 for temperatures of the order of 10 8 K, and ∼4.5 for temperatures of the order of 5 × 10 7 K, depending on the magnetic field configuration. The nonlinear Hall term is crucial to this amplification: without it, the magnetic field decay is only enhanced by a factor ≲2 even at 5 × 10 7 K. These effects generate a high wavenumber plateau in the magnetic energy spectrum associated with the small-scale de Haas–van Alphen oscillations. Our results suggest that this mechanism could help explain the magnetar heating problem, though due to the effect's temperature-dependence, full magneto-thermal evolution simulations in a realistic stellar model are needed to judge whether it is viable explanation.

Rau, Peter B. (ORCID:0000000152209277)↗
DOE OSTI2022

Nonresonant particle acceleration in strong turbulence: Comparison to kinetic and MHD simulations

Collisionless, magnetized turbulence offers a promising framework for the generation of nonthermal high-energy particles in various astrophysical sites. Yet, the detailed mechanism that governs particle acceleration has remained subject to debate. By means of 2D and 3D particle-in-cell, as well as 3D (incompressible) magnetohydrodynamic (MHD) simulations, we test here a recent model of nonresonant particle acceleration in strongly magnetized turbulence, which ascribes the energization of particles to their continuous interaction with the random velocity flow of the turbulence, in the spirit of the original Fermi model. To do so, we compare, for a large number of particles that were tracked in the simulations, the predicted and the observed histories of particles momenta. The predicted history is that derived from the model, after extracting from the simulations, at each point along the particle trajectory, the three force terms that control acceleration: the acceleration of the field line velocity projected along the field line direction, its shear projected along the same direction, and its transverse compressive part. Overall, we find a clear correlation between the model predictions and the numerical experiments, indicating that this nonresonant model can successfully account for the bulk of particle energization through Fermi-type processes in strongly magnetized turbulence. Additionally we also observe that the parallel shear contribution tends to dominate the physics of energization in the particle-in-cell simulations, while in the magnetohydrodynamic incompressible simulation, both the parallel shear and the transverse compressive term provide about equal contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

First M87 Event Horizon Telescope Results. VIII. Magnetic Field Structure near The Event Horizon

Event Horizon Telescope (EHT) observations at 230 GHz have now imaged polarized emission around the supermassive black hole in M87 on event-horizon scales. This polarized synchrotron radiation probes the structure of magnetic fields and the plasma properties near the black hole. Here we compare the resolved polarization structure observed by the EHT, along with simultaneous unresolved observations with the Atacama Large Millimeter/submillimeter Array, to expectations from theoretical models. The low fractional linear polarization in the resolved image suggests that the polarization is scrambled on scales smaller than the EHT beam, which we attribute to Faraday rotation internal to the emission region. We estimate the average density n_e of order 10^4-7 cm-3, magnetic field strength B of order 1-30 G, and electron temperature Te of order (1-12) x 10^10 K of the radiating plasma in a simple one-zone emission model. We show that the net azimuthal linear polarization pattern may result from organized, poloidal magnetic fields in the emission region. In a quantitative comparison with a large library of simulated polarimetric images from general relativistic magnetohydrodynamic (GRMHD) simulations, we identify a subset of physical models that can explain critical features of the polarimetric EHT observations while producing a relativistic jet of sufficient power. The consistent GRMHD models are all of magnetically arrested accretion disks, where near-horizon magnetic fields are dynamically important. We use the models to infer a mass accretion rate onto the black hole in M87 of (3-20) x 10^-4 Msun yr-1.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Turbulent Regimes in Collisions of 3D Alfvén-wave Packets

Using three-dimensional gyrofluid simulations, we revisit the problem of Alfvén-wave (AW) collisions as building blocks of the Alfvénic turbulent cascade and their interplay with magnetic reconnection at magnetohydrodynamic (MHD) scales. Depending on the large-scale value of the nonlinearity parameter χ 0 (the ratio between the AW linear propagation time and nonlinear turnover time), different regimes are observed. For strong nonlinearities (χ 0 ~ 1), turbulence is consistent with a dynamically aligned, critically balanced cascade—fluctuations exhibit a scale-dependent alignment $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/4}$, resulting in a ${k}_{\perp }^{-3/2}$ spectrum and ${k}_{\parallel }\propto {k}_{\perp }^{1/2}$ spectral anisotropy. At weaker nonlinearities (small χ 0 ), a spectral break marking the transition between a large-scale weak regime and a small-scale ${k}_{\perp }^{-11/5}$ tearing-mediated range emerges, implying that dynamic alignment occurs also for weak nonlinearities. At χ 0 < 1 the alignment angle ${\theta }_{{k}_{\perp }}$ shows a stronger scale dependence than in the χ 0 ~ 1 regime, namely $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/2}$ at χ 0 ~ 0.5, and $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1}$ at χ 0 ~ 0.1. Dynamic alignment in the weak regime also modifies the large-scale spectrum, scaling approximately as ${k}_{\perp }^{-3/2}$ for χ 0 ~ 0.5 and as ${k}_{\perp }^{-1}$ for χ 0 ~ 0.1. A phenomenological theory of dynamically aligned turbulence at weak nonlinearities that can explain these spectra and the transition to the tearing-mediated regime is provided; at small χ 0 , the strong scale dependence of the alignment angle combines with the increased lifetime of turbulent eddies to allow tearing to onset and mediate the cascade at scales that can be larger than those predicted for a critically balanced cascade by several orders of magnitude. Such a transition to tearing-mediated turbulence may even supplant the usual weak-to-strong transition.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Bridging the Bondi and Event Horizon Scales: 3D GRMHD Simulations Reveal X-shaped Radio Galaxy Morphology

X-shaped radio galaxies (XRGs) produce misaligned X-shaped jet pairs and make up ≲10% of radio galaxies. XRGs are thought to emerge in galaxies featuring a binary supermassive black hole (SMBH), SMBH merger, or large-scale ambient medium asymmetry. We demonstrate that XRG morphology can naturally form without such special, preexisting conditions. Our 3D general-relativistic magnetohydrodynamic (GRMHD) simulation for the first time follows magnetized rotating gas from outside the SMBH sphere of influence of radius RB to the SMBH of gravitational radius Rg at the largest scale separation, RB/Rg = 103, to date. Initially, our axisymmetric system of constant-density hot gas contains a weak vertical magnetic field and rotates in the equatorial plane of a rapidly spinning SMBH. We seed the gas with small-scale 2% level pressure perturbations. Infalling gas forms an accretion disk, and the SMBH launches relativistically magnetized collimated jets reaching well outside RB. Under the pressure of the infalling gas, the jets intermittently turn on and off, erratically wobble, and inflate pairs of cavities in different directions, resembling an X-shaped jet morphology. Synthetic X-ray images reveal multiple pairs of jet-powered shocks and cavities. Large-scale magnetic flux accumulates on the SMBH, becomes dynamically important, and leads to a magnetically arrested disk state. The SMBH accretes at 2% of the Bondi rate ($\dot{M}\simeq 2.4\times {10}^{-3}{M}_{\odot }\,{\mathrm{yr}}^{-1}$ for M87*) and launches twin jets at η = 150% efficiency. These jets are powerful enough (P jets ≃ 2 × 10 44 erg s -1 ) to escape along the SMBH spin axis and end the short-lived intermittent jet state, whose transient nature can account for the rarity of XRGs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Multizone Modeling of Black Hole Accretion and Feedback in 3D GRMHD: Bridging Vast Spatial and Temporal Scales

Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. Elaborating on H. Cho et al., we describe and test a "multizone" technique, which is designed to tackle this difficult problem in three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. While short-timescale variability should be interpreted with caution, the method is demonstrated to be well-suited for finding dynamical steady states over a wide dynamic range. We simulate accretion on a nonspinning SMBH ($a\ast$ = 0) using initial conditions and the external galactic potential from a large-scale galaxy simulation and achieve a steady state over eight decades in radius. As found in H. Cho et al., the density scales with radius as ρ ∝ r –1 inside the Bondi radius R B , which is located at R B = 2 × 10 5 r g (≈60 pc for M87), where r g is the gravitational radius of the SMBH; the plasma-β is ~ unity, indicating an extended magnetically arrested state; the mass accretion rate $\dot{M}$ is ≈1% of the analytical Bondi accretion rate ${\dot{M}}_{{\rm{B}}};$ and there is continuous energy feedback out to ≈100R B (or beyond > kpc) at a rate $\approx 0.02\dot{M}{c}^{2}$. Surprisingly, no ordered rotation in the external medium survives as the magnetized gas flows to smaller radii, and the final steady solution is very similar to when the exterior has no rotation. Using the multizone method, we simulate GRMHD accretion over a wide range of Bondi radii, R B ~ 10 2 –10 7 r g , and find that $\dot{M}/{\dot{M}}_{{B}}\approx {({R}_{{B}}/6\,{r}_{g})}^{-0.5}$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

SILCC – VII. Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities

ABSTRACT We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities Σgas = 10, 30, 50, and 100 $\mathrm{M_\odot \, pc^{-2}}$. The silcc project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionizing UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between Σgas and the star formation rate surface density ΣSFR. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher Σgas and the mass loading factors of the CR-supported outflows are of order unity independent of ΣSFR. Similar to observations, vertical velocity dispersions of the warm ionized medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as $\sigma _\mathrm{z} \propto \Sigma _\mathrm{SFR}^a$, with a ∼ 0.20. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor ∼2.2 higher than that of the CNM, in agreement with local observations. For $\Sigma _\mathrm{SFR} \gtrsim 1.5 \times 10^{-2}\, \mathrm{M}_\odot \, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$ the WIM motions become supersonic.

Astronomy & Astrophysics↗
DOE OSTI2024

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

3D Numerical Simulation of Kink-driven Rayleigh–Taylor Instability Leading to Fast Magnetic Reconnection

Fast magnetic reconnection involving non-MHD microscale physics is thought to underlie both solar eruptions and laboratory plasma current disruptions. While there is extensive research on both the MHD macroscale physics and the non-MHD microscale physics, the process by which large-scale MHD couples to the microscale physics is not well understood. An MHD instability cascade from a kink to a secondary Rayleigh–Taylor instability in the Caltech astrophysical jet laboratory experiment provides new insights into this coupling and motivates a 3D numerical simulation of this transition from large to small scale. A critical finding from the simulation is that the axial magnetic field inside the current-carrying dense plasma must exceed the field outside. In addition, the simulation verifies a theoretical prediction and experimental observation that, depending on the strength of the effective gravity produced by the primary kink instability, the secondary instability can be Rayleigh–Taylor or mini-kink. Finally, it is shown that the kink-driven Rayleigh–Taylor instability generates a localized electric field sufficiently strong to accelerate electrons to very high energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing with 3D Adaptive Mesh Refinement and Local Adaptive Time Stepping

General relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a GPU-accelerated GRMHD code H-AMR with multifaceted optimizations that, collectively, accelerate computation by 2–5 orders of magnitude for a wide range of applications. First, it introduces a spherical grid with 3D adaptive mesh refinement that operates in each of the three dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-resolution computational performance. Second, we demonstrate that local adaptive time stepping on a logarithmic spherical-polar grid accelerates computation by a factor of ≲10 compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of ~10 9 zone cycles per second per node on 5400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We illustrate H-AMR's computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of 13,440 × 4608 × 8092 cells and a total of ≲22 billion cells and ~0.65 × 10 8 time steps, it is among the largest astrophysical simulations ever performed. We find that frame dragging by the black hole tears up the disk into two independently precessing subdisks. The innermost subdisk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Minidisk Accretion onto Spinning Black Hole Binaries: Quasi-periodicities and Outflows

We perform a full 3D general relativistic magnetohydrodynamical (GRMHD) simulation of an equal-mass, spinning, binary black hole approaching merger, surrounded by a circumbinary disk and with a minidisk around each black hole. For this purpose, we evolve the ideal GRMHD equations on top of an approximated spacetime for the binary that is valid in every position of space, including the black hole horizons, during the inspiral regime. We use relaxed initial data for the circumbinary disk from a previous long-term simulation, where the accretion is dominated by a m = 1 overdensity called the lump. We compare our new spinning simulation with a previous non-spinning run, studying how spin influences the minidisk properties. We analyze the accretion from the inner edge of the lump to the black hole, focusing on the angular momentum budget of the fluid around the minidisks. We find that minidisks in the spinning case have more mass over a cycle than the non-spinning case. However, in both cases we find that most of the mass received by the black holes is delivered by the direct plunging of material from the lump. We also analyze the morphology and variability of the electromagnetic fluxes, and we find they share the same periodicities of the accretion rate. In the spinning case, we find that the outflows are stronger than the non-spinning case. Our results will be useful to understand and produce realistic synthetic light curves and spectra, which can be used in future observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

On the Geometry of the Near-core Magnetic Field in Massive Stars

It is well-known that the cores of massive stars sustain a stellar dynamo with a complex magnetic field configuration. However, the same cannot be said for the field's strength and geometry at the convective–radiative boundary, which are crucial when performing asteroseismic inference. In this Letter, we present 3D magnetohydrodynamic (MHD) simulations of a 7 M ⊙ mid-main-sequence star, with particular attention given to the convective–radiative boundary in the near-core region. Our simulations reveal that the toroidal magnetic field is significantly stronger than the poloidal field in this region, contrary to recent assumptions. Moreover, the rotational shear layer, also important for asteroseismic inference, is specifically confined within the extent of the Brunt–Väisälä frequency peak. These results, which are based on the inferred properties of HD 43317, have widespread implications for asteroseismic studies of rotation, mixing, and magnetism in stars. While we expect our results to be broadly applicable across stars with similar Brunt–Väisälä frequency profiles and stellar masses, we also expect the MHD parameters (e.g., Re m ) and the initial stellar rotation rate to impact the geometry of the field and differential rotation at the convective–radiative interface.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2025

Self-similar Cosmic-Ray Transport in High-resolution Magnetohydrodynamic Turbulence

We study the propagation of cosmic rays (CRs) through a simulation of magnetohydrodynamic (MHD) turbulence at an unprecedented resolution of 10,240 3 . We drive turbulence that is subsonic and super-Alfvénic, characterized by δB rms /B 0 = 2. The high resolution enables an extended inertial range such that the Alfvén scale l A , where δB(l A ) ≈ B 0 , is well resolved. This allows us to properly capture how the cascade transitions from large amplitudes on large scales to small amplitudes on small scales. We find that sharp bends in the magnetic field are key mediators of particle transport even on small scales via resonant curvature scattering. We further find that particle scattering in the turbulence shows strong hints of self-similarity: (1) the diffusion has a weak energy dependence over almost two decades in particle energy, and (2) the particles’ random walk exhibits a broad power-law distribution of collision times such that the diffusion is dominated by the rarest, long-distance excursions. Our results suggest that large-amplitude MHD turbulence can provide efficient scattering over a wide range of CR energies and may help explain many CR observations above a ∼TeV: the flattening of the B/C spectrum, the hardening of CR primary spectra, and the weak dependence of arrival anisotropy on CR energy.

Kempski, Philipp [Princeton Univ., NJ (United Stat↗
DOE OSTI2021

Exosphere Modeling of Proxima b: A Case Study of Photochemical Escape with a Venus-like Atmosphere

Exoplanets orbiting M dwarfs within habitable zones are exposed to stellar environments more extreme than that terrestrial planets experience in our solar system, which can significantly impact the atmospheres of the exoplanets and affect their habitability and sustainability. This study provides the first prediction of hot oxygen corona structure and the associated photochemical loss from a 1 bar CO 2 -dominated atmosphere of a Venus-like rocky exoplanet, where dissociative recombination of O 2 + ions is assumed to be the major source reaction for the escape of neutral O atoms and formation of the hot O corona (or exospheres) as on Mars and Venus. We employ a 3D Monte Carlo code to simulate the exosphere of Proxima Centauri b (PCb) based on the ionosphere simulated by a 3D magnetohydrodynamic model. Our simulation results show that variability of the stellar wind dynamic pressure over one orbital period of PCb does not affect the overall spatial structure of the hot O corona but contributes to the change in the global hot O escape rate that varies by an order of magnitude. The escape increases dramatically when the planet possesses its intrinsic magnetic fields as the ionosphere becomes more extended with the presence of a global magnetic field. Furthermore, the extended hot O corona may lead to a more extended H exosphere through collisions between thermal H and hot O, which exemplifies the importance of considering nonthermal populations in exospheres to interpret future observations.

79 ASTRONOMY AND ASTROPHYSICS↗