Search NASA⌕ Search

SEARCH · Search NASA

Results for “Space simulators”

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

The Feasibility of Directly Imaging Nearby Cold Jovian Planets with MIRI/JWST

The upcoming launch of the James Webb Space Telescope (JWST) will dramatically increase our understanding of exoplanets, particularly through direct imaging. Microlensing and radial velocity surveys indicate that some M dwarfs host long-period giant planets. Some of these planets will likely be just a few parsecs away and a few astronomical units from their host stars, a parameter space that cannot be probed by existing high-contrast imagers. We studied whether the coronagraphs on the Mid-infrared Instrument on JWST can detect Jovian-type planets around nearby M dwarfs. For a sample of 27 very nearby M dwarfs, we simulated a sample of Saturn–Jupiter-mass planets with three atmospheric configurations and three orbital separations, observed in three different filters. We found that the f1550c 15.5 μm filter is best suited for detecting Jupiter-like planets. Jupiter-like planets with patchy cloud cover, 2 au from their star, are detectable at 15.5 μm around 14 stars in our sample, while Jupiters with clearer atmospheres are detectable around all stars in the sample. Saturns were most detectable at 10.65 and 11.4 μm (f1065c and f1140c filters), but only with cloud-free atmospheres and within 3 pc (six stars). Surveying all 27 stars would take <170 hr of JWST integration time, or just a few hours for a shorter survey of the most favorable targets. There is one potentially detectable known planet in our sample: GJ 832 b. Observations aimed at detecting this planet should occur in 2024–2026, when the planet is maximally separated from the star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Modeling a ring magnet in ALEGRA

We show here that Sandia's ALEGRA software can be used to model a permanent magnet in 2D and 3D, with accuracy matching that of the open-source commercial software FEMM. This is done by conducting simulations and experimental measurements for a commercial-grade N42 neodymium alloy ring magnet with a measured magnetic field strength of approximately 0.4 T in its immediate vicinity. Transient simulations using ALEGRA and static simulations using FEMM are conducted. Comparisons are made between simulations and measurements, and amongst the simulations, for sample locations in the steady-state magnetic field. The comparisons show that all models capture the data to within 7%. The FEMM and ALEGRA results agree to within approximately 2%. The most accurate solutions in ALEGRA are obtained using quadrilateral or hexahedral elements. In the case where iron shielding disks are included in the magnetized space, ALEGRA simulations are considerably more expensive because of the increased magnetic diffusion time, but FEMM and ALEGRA results are still in agreement. The magnetic field data are portable to other software interfaces using the Exodus file format.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An Alfvenic reconnecting plasmoid thruster

A new concept for the generation of thrust for space propulsion is introduced. Energetic thrust is generated in the form of plasmoids (confined plasma in closed magnetic loops) when magnetic helicity (linked magnetic field lines) is injected into an annular channel. Using a novel configuration of static electric and magnetic fields, the concept utilizes a current-sheet instability to spontaneously and continuously create plasmoids via magnetic reconnection. Here, the generated low-temperature plasma is simulated in a global annular geometry using the extended magnetohydrodynamic model. Because the system-size plasmoid is an Alfvenic outflow from the reconnection site, its thrust is proportional to the square of the magnetic field strength and does not ideally depend on the mass of the ion species of the plasma. Exhaust velocities in the range of 20 to $500\ \textrm {km}\ \textrm {s}^{-1}$ , controllable by the coil currents, are observed in the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Atmospheric Escape From TOI-700 d: Venus versus Earth Analogs

The recent discovery of an Earth-sized planet (TOI-700 d) in the habitable zone of an early-type M-dwarf by the Transiting Exoplanet Survey Satellite constitutes an important advance. In this Letter, we assess the feasibility of this planet to retain an atmosphere – one of the chief ingredients for surface habitability – over long timescales by employing state-of-the-art magnetohydrodynamic models to simulate the stellar wind and the associated rates of atmospheric escape. We take two major factors into consideration, namely, the planetary atmospheric composition and magnetic field. In all cases, we determine that the atmospheric ion escape rates are potentially a few orders of magnitude higher than the inner Solar system planets, but TOI-700 d is nevertheless capable of retaining a 1 bar atmosphere over gigayear timescales for certain regions of the parameter space. The simulations show that the unmagnetized TOI-700 d with a 1 bar Earth-like atmosphere could be stripped away rather quickly (< 1 gigayear), while the unmagnetized TOI-700 d with a 1 bar CO 2 -dominated atmosphere could persist for many billions of years; we find that the magnetized Earth-like case falls in between these two scenarios. We also discuss the prospects for detecting radio emission of the planet (thereby constraining its magnetic field) and discerning the presence of an atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Impact of bar resonances in the velocity–space distribution of the solar neighbourhood stars in a self-consistent N -body Galactic disc simulation

ABSTRACT The velocity–space distribution of the solar neighbourhood stars shows complex substructures. Most of the previous studies use static potentials to investigate their origins. Instead we use a self-consistent N-body model of the Milky Way, whose potential is asymmetric and evolves with time. In this paper, we quantitatively evaluate the similarities of the velocity–space distributions in the N-body model and that of the solar neighbourhood, using Kullback–Leibler divergence (KLD). The KLD analysis shows the time evolution and spatial variation of the velocity–space distribution. The KLD fluctuates with time, which indicates the velocity–space distribution at a fixed position is not always similar to that of the solar neighbourhood. Some positions show velocity–space distributions with small KLDs (high similarities) more frequently than others. One of them locates at $(R,\phi)=(8.2\,\,\rm{\mathrm{kpc}}, 30^\circ)$, where R and ϕ are the distance from the galactic centre and the angle with respect to the bar’s major axis, respectively. The detection frequency is higher in the inter-arm regions than in the arm regions. In the velocity maps with small KLDs, we identify the velocity–space substructures, which consist of particles trapped in bar resonances. The bar resonances have significant impact on the stellar velocity–space distribution even though the galactic potential is not static.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Unified many-worlds browsing of arbitrary physics-based animations

Manually tuning physics-based animation parameters to explore a simulation outcome space or achieve desired motion outcomes can be notoriously tedious. This problem has motivated many sophisticated and specialized optimization-based methods for fine-grained (keyframe) control, each of which are typically limited to specific animation phenomena, usually complicated, and, unfortunately, not widely used. In this paper, we propose Unified Many-Worlds Browsing (UMWB), a practical method for sample-level control and exploration of physics-based animations. Our approach supports browsing of large simulation ensembles of arbitrary animation phenomena by using a unified volumetric WORLDPACK representation based on spatiotemporally compressed voxel data associated with geometric occupancy and other low-fidelity animation state. Beyond memory reduction, the WORLDPACK representation also enables unified query support for interactive browsing: it provides fast evaluation of approximate spatiotemporal queries, such as occupancy tests that find ensemble samples ("worlds") where material is either IN or NOT IN a user-specified spacetime region. WORLDPACKS also support real-time hardware-accelerated voxel rendering by exploiting the spatially hierarchical and temporal RLE raster data structure. Our UMWB implementation supports interactive browsing (and offline refinement) of ensembles containing thousands of simulation samples, and fast spatiotemporal queries and ranking. We show UMWB results using a wide variety of physics-based animation phenomena---not just JELL-O ® .

Computer Science↗
DOE OSTI2025

Digital Twin Framework for PIP-II Linac: AI-Driven Multi-Scale Modeling from Ion Source to 800 MeV

The PIP-II linac will enable >1.2 MW beam power for DUNE, requiring unprecedented operational reliability across its warm front-end (RFQ, MEBT) and five distinct SRF sections operating at 162.5/325/650 MHz. We present a comprehensive digital twin framework uniquely combining a fully differentiable fast beam transport code with neural network surrogates trained on high-fidelity PIC simulations, capturing space charge and nonlinear dynamics beyond traditional envelope codes while achieving 10⁴ speedup at <1% accuracy. End-to-end differentiability enables gradient-based optimization across 500+ parameters simultaneously previously impossible with conventional tools while the model incorporates static/dynamic errors and serves as a virtual commissioning platform for diverse hardware integration. The framework facilitates reinforcement learning for pulsed/CW mode transitions, predictive maintenance through anomaly detection, and autonomous tuning algorithm development with real-time execution capability. Validation against physics simulations shows excellent agreement for the front-end, with initial results demonstrating potential for 30% commissioning time reduction and proactive fault mitigation, providing a scalable blueprint for operating next-generation high-intensity accelerators.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗
DOE OSTI2025

Digital Twin Framework for PIP-II Linac: AI-Driven Multi-Scale Modeling from Ion Source to 800 MeV

The PIP-II linac will enable >1.2 MW beam power for DUNE, requiring unprecedented operational reliability across its warm front-end (RFQ, MEBT) and five distinct SRF sections operating at 162.5/325/650 MHz. We present a comprehensive digital twin framework uniquely combining a fully differentiable fast beam transport code with neural network surrogates trained on high-fidelity PIC simulations, capturing space charge and nonlinear dynamics beyond traditional envelope codes while achieving 10⁴× speedup at <1% accuracy. End-to-end differentiability enables gradient-based optimization across 500+ parameters simultaneously—previously impossible with conventional tools—while the model incorporates static/dynamic errors and serves as a virtual commissioning platform for diverse hardware integration. The framework facilitates reinforcement learning for pulsed/CW mode transitions, predictive maintenance through anomaly detection, and autonomous tuning algorithm development with real-time execution capability. Validation against physics simulations shows excellent agreement for the front-end, with initial results demonstrating potential for 30% commissioning time reduction and proactive fault mitigation, providing a scalable blueprint for operating next-generation high-intensity accelerators.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗
DOE OSTI2020

Modeling Stochastic Variability in Multiband Time-series Data

In preparation for the era of time-domain astronomy with upcoming large-scale surveys, we propose a state-space representation of a multivariate damped random walk process as a tool to analyze irregularly-spaced multifilter light curves with heteroscedastic measurement errors. We adopt a computationally efficient and scalable Kalman filtering approach to evaluate the likelihood function, leading to maximum O(k 3 n) complexity, where k is the number of available bands and n is the number of unique observation times across the k bands. This is a significant computational advantage over a commonly used univariate Gaussian process that can stack up all multiband light curves in one vector with maximum O(k 3 n 3 ) complexity. Using such efficient likelihood computation, we provide both maximum likelihood estimates and Bayesian posterior samples of the model parameters. Three numerical illustrations are presented: (i) analyzing simulated five-band light curves for a comparison with independent single-band fits; (ii) analyzing five-band light curves of a quasar obtained from the Sloan Digital Sky Survey Stripe 82 to estimate short-term variability and timescale; (iii) analyzing gravitationally lensed g- and r-band light curves of Q0957+561 to infer the time delay. Two R packages, Rdrw and timedelay, are publicly available to fit the proposed models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Updated Parameters and a New Transmission Spectrum of HD 97658b

Recent years have seen increasing interest in the characterization of sub-Neptune-sized planets because of their prevalence in the Galaxy, contrasted with their absence in our solar system. HD 97658 is one of the brightest stars hosting a planet of this kind, and we present the transmission spectrum of this planet by combining four Hubble Space Telescope transits, 12 Spitzer/IRAC transits, and eight MOST transits of this system. Our transmission spectrum has a higher signal-to-noise ratio than those from previous works, and the result suggests that the slight increase in transit depth from wavelength 1.1–1.7 μm reported in previous works on the transmission spectrum of this planet is likely systematic. Nonetheless, our atmospheric modeling results are inconclusive, as no model provides an excellent match to our data. Nonetheless, we find that atmospheres with high C/O ratios (C/O ≳ 0.8) and metallicities of ≳100× solar metallicity are favored. We combine the mid-transit times from all of the new Spitzer and MOST observations and obtain an updated orbital period of P = 9.489295 ± 0.000005, with a best-fit transit time center at T {sub 0} = 2456361.80690 ± 0.00038 (BJD). No transit timing variations are found in this system. We also present new measurements of the stellar rotation period (34 ± 2 days) and stellar activity cycle (9.6 yr) of the host star HD 97658. Finally, we calculate and rank the Transmission Spectroscopy Metric of all confirmed planets cooler than 1000 K and with sizes between 1 R {sub ⊕} and 4 R {sub ⊕}. We find that at least a third of small planets cooler than 1000 K can be well characterized using James Webb Space Telescope, and of those, HD 97658b is ranked fifth, meaning that it remains a high-priority target for atmospheric characterization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Pluto's Ultraviolet Spectrum, Surface Reflectance, and Airglow Emissions

During the New Horizons spacecraft’s encounter with Pluto, the Alice ultraviolet spectrograph conducted a series of observations that detected emissions from both the interplanetary medium (IPM) and Pluto. In the direction of Pluto, the IPM was found to be 133.4 ± 0.6 R at Lyα, 0.24 ± 0.02 R at Lyβ, and <0.10 R at He i 584 Å. We analyzed 3900 s of data obtained shortly before closest approach to Pluto and detect airglow emissions from H i, N i, N ii, N{sub 2}, and CO above the disk of Pluto. We find Pluto’s brightness at Lyα to be 29.3 ± 1.9 R, in good agreement with preencounter estimates. The detection of the N ii multiplet at 1085 Å marks the first direct detection of ions in Pluto’s atmosphere. We do not detect any emissions from noble gases and place a 3σ upper limit of 0.14 R on the brightness of the Ar i 1048 Å line. We compare preencounter model predictions and predictions from our own airglow model, based on atmospheric profiles derived from the solar occultation observed by New Horizons, to the observed brightness of Pluto’s airglow. Although completely opaque at Lyα, Pluto’s atmosphere is optically thin at wavelengths longer than 1425 Å. Consequently, a significant amount of solar far-UV light reaches the surface, where it can participate in space weathering processes. From the brightness of sunlight reflected from Pluto, we find the surface has a reflectance factor (I/F) of 17% between 1400 and 1850 Å. We also report the first detection of a C{sub 3} hydrocarbon molecule, methylacetylene, in absorption, at a column density of ∼5 × 10{sup 15} cm{sup −2}, corresponding to a column-integrated mixing ratio of 1.6 × 10{sup −6}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Analysis of strategies to meet ASHRAE S241 infectious aerosol control targets by space type and region using EnergyPlus™ simulations

Professional organizations, such as ASHRAE, have recently proposed new voluntary standards for controlling infectious aerosols. Specifically, ASHRAE Standard 241 (S241) defines equivalent clean air targets for a range of space types that can be met through combinations of mitigation measures. This paper seeks to inform the selection of measures by space type and climate zone through building simulations that quantify the impacts of increased outdoor air ventilation; increasing filtration and/or adding germicidal ultraviolet (GUV) in the central heating, ventilation, and air-conditioning (HVAC) systems; or using portable air cleaners (PACs), upper room GUV, or whole room GUV approaches to increase equivalent clean air delivery. The measures are assessed individually and in practical combinations for their ability to meet S241 in four space types (offices, classrooms, dining areas, and healthcare waiting rooms) using prototype buildings modeled using EnergyPlus. The mitigation measures are compared holistically against baseline building operations using metrics for equivalent clean air, energy, and comfort. The results in this study show that single measures can meet S241 for offices with minimal impacts on energy and comfort, while either upper room GUV systems or combinations of measures such as MERV 13 HVAC filtration with PACs are needed to meet S241 for classrooms. The dining area and healthcare waiting room targets cannot be met in this study when assuming design occupancy and MS2 as the challenge agent, even when combining multiple measures together. This paper provides valuable considerations when designing measures to meet S241 for a range of spaces and scenarios.

GUV↗
DOE OSTI2020

The Multiplanet System TOI-421: A Warm Neptune and a Super Puffy Mini-Neptune Transiting a G9 V Star in a Visual Binary

We report the discovery of a warm Neptune and a hot sub-Neptune transiting TOI-421 (BD-14 1137, TIC 94986319), a bright (V = 9.9) G9 dwarf star in a visual binary system observed by the Transiting Exoplanet Survey Satellite (TESS) space mission in Sectors 5 and 6. We performed ground-based follow-up observations—comprised of Las Cumbres Observatory Global Telescope transit photometry, NIRC2 adaptive optics imaging, and FIbre-fed Echellé Spectrograph, CORALIE, High Accuracy Radial velocity Planet Searcher, High Resolution Échelle Spectrometer, and Planet Finder Spectrograph high-precision Doppler measurements—and confirmed the planetary nature of the 16 day transiting candidate announced by the TESS team. We discovered an additional radial velocity signal with a period of five days induced by the presence of a second planet in the system, which we also found to transit its host star. We found that the inner mini-Neptune, TOI-421 b, has an orbital period of P{sub b} = 5.19672 ± 0.00049 days, a mass of M{sub b} = 7.17 ± 0.66 M{sub ⊕}, and a radius of R{sub b} = 2.68{sub −0.18}{sup +0.19} R{sub ⊕}, whereas the outer warm Neptune, TOI-421 c, has a period of P{sub c} = 16.06819 ± 0.00035 days, a mass of M{sub c} = 16.42{sub −1.04}{sup +1.06} M{sub ⊕}, a radius of R{sub c} = 5.09{sub −0.15}{sup +0.16} R{sub ⊕}, and a density of ρ{sub c} = 0.685{sub −0.072}{sup +0.080} g cm{sup −3}. With its characteristics, the outer planet (ρ{sub c} = 0.685{sub −0.072}{sup +0.080} g cm{sup −3}) is placed in the intriguing class of the super-puffy mini-Neptunes. TOI-421 b and TOI-421 c are found to be well-suited for atmospheric characterization. Our atmospheric simulations predict significant Lyα transit absorption, due to strong hydrogen escape in both planets, as well as the presence of detectable CH{sub 4} in the atmosphere of TOI-421 c if equilibrium chemistry is assumed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Astraea: Predicting Long Rotation Periods with 27 Day Light Curves

The rotation periods of planet-hosting stars can be used for modeling and mitigating the impact of magnetic activity in radial velocity measurements and can help constrain the high-energy flux environment and space weather of planetary systems. Millions of stars and thousands of planet hosts are observed with the Transiting Exoplanet Survey Satellite (TESS). However, most will only be observed for 27 contiguous days in a year, making it difficult to measure rotation periods with traditional methods. This is especially problematic for field M dwarfs, which are ideal candidates for exoplanet searches, but which tend to have periods in excess of the 27 day observing baseline. We present a new tool, Astraea, for predicting long rotation periods from short-duration light curves combined with stellar parameters from Gaia DR2. Using Astraea, we can predict the rotation periods from Kepler 4 yr light curves with 13% uncertainty overall (and a 9% uncertainty for periods >30 days). By training on 27 day Kepler light-curve segments, Astraea can predict rotation periods up to 150 days with 9% uncertainty (5% for periods >30 days). After training this tool on these 27 day Kepler light-curve segments, we applied Astraea to real TESS data. For the 195 stars that were observed by both Kepler and TESS, we were able to predict the rotation periods with 55% uncertainty despite the wild differences in systematics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Search for FeH in Hot-Jupiter Atmospheres with High-dispersion Spectroscopy

Most of the molecules detected thus far in exoplanet atmospheres, such as water and CO, are present for a large range of pressures and temperatures. In contrast, metal hydrides exist in much more specific regimes of parameter space, and so can be used as probes of atmospheric conditions. Iron hydride (FeH) is a dominant source of opacity in low-mass stars and brown dwarfs, and evidence for its existence in exoplanets has recently been observed at low resolution. We performed a systematic search of archival CARMENES near-infrared data for signatures of FeH during transits of 12 exoplanets. These planets span a large range of equilibrium temperatures (600 ≲ T eq ≲ 4000 K) and surface gravities (2.5 ≲ logg ≲ 3.5). We did not find a statistically significant FeH signal in any of the atmospheres, but obtained potential low-confidence signals (signal-to-noise ratio ~ 3) in two planets, WASP-33b and MASCARA-2b. Previous modeling of exoplanet atmospheres indicate that the highest volume mixing ratios (VMRs) of 10 -7 to 10 -9 are expected for temperatures between 1800 and 3000 K and log g≳3. The two planets for which we find low-confidence signals are in the regime where strong FeH absorption is expected. We performed injection and recovery tests for each planet and determined that FeH would be detected in every planet for VMRs &NestedGreaterGreater;10 -6 , and could be detected in some planets for VMRs as low as 10 -9.5 . Additional observations are necessary to conclusively detect FeH and assess its role in the temperature structures of hot-Jupiter atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Lpp positions peptidoglycan at the AcrA-TolC interface in the AcrAB-TolC multidrug efflux pump

The multidrug efflux pumps of Gram-negative bacteria are a class of complexes that span the periplasm, coupling both the inner and outer membranes to expel toxic molecules. The best-characterized example of these tripartite pumps is the AcrAB-TolC complex of Escherichia coli. However, how the complex interacts with the peptidoglycan (PG) cell wall, which is anchored to the outer membrane (OM) by Braun’s lipoprotein (Lpp), is still largely unknown. In this work, we present molecular dynamics simulations of a complete, atomistic model of the AcrAB-TolC complex with the inner membrane, OM, and PG layers all present. We find that the PG localizes to the junction of AcrA and TolC, in agreement with recent cryo-tomography data. Free-energy calculations reveal that the positioning of PG is determined by the length and conformation of multiple Lpp copies anchoring it to the OM. The distance between the PG and OM measured in cryo-electron microscopy images of wild-type E. coli also agrees with the simulation-derived spacing. Finally, sequence analysis of AcrA suggests a conserved role for interactions with PG in the assembly and stabilization of efflux pumps, one that may extend to other trans-envelope complexes as well.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Anisotropic Heating and Parallel Heat Flux in Electron-only Magnetic Reconnection with Intense Guide Fields

Electron-only reconnection (E-REC) is a process recently observed in the Earth’s magnetosheath, where magnetic reconnection occurs at electron kinetic scales, and ions do not couple to the reconnection process. Electron-only reconnection is likely to have a significant impact on the energy conversion and dissipation of turbulence cascades at kinetic scales in some settings. This paper investigates E-REC under different intensities of strong guide fields (the ratio between the guide field and the in-plane asymptotic field strength is 5, 10 and 20, respectively) via two-dimensional fully kinetic particle-in-cell simulations, focusing on electron heating. The simulations are initialized with a force-free current sheet equilibrium under various intensities of strong guide fields. Similarly to previous experimental studies, electron temperature anisotropy along separatrices is observed, which is found to be mainly caused by the variations of parallel temperature. Both regions of anisotropy and parallel temperature increase/decrease along separatrices become thinner with increasing guide fields. Besides, we find a transition from a quadrupolar to a hexapolar (six-polar) to an octopolar (eight-polar) structure in temperature anisotropy and parallel temperature as the guide field intensifies. Non-Maxwellian electron velocity distribution functions (EVDFs) at different locations in the three simulations are observed. Our results show that parallel electron velocity varies notably with different guide field intensities and finite parallel electron heat flux density is observed. The three simulations exhibit features of the Chew–Goldberger–Low theory, with the level of consistency increasing as the guide field strength increases. This explains the electron parallel temperature variations and the shape of the EVDFs observed along the separatrices. This work may provide insights into the understanding of electron heating and parallel heat flux density in E-REC observed in the turbulent magnetosheath.

79 ASTRONOMY AND ASTROPHYSICS↗