An official website of the United States government
Abstract We compare hybrid (kinetic proton, fluid electron) and particle-in-cell (kinetic proton, kinetic electron) simulations of the solar wind termination shock with parameters similar to those observed by Voyager 2 during its crossing. The steady-state results show excellent agreement between the downstream variations in the density, plasma velocity, and magnetic field. The quasi-perpendicular shock accelerates interstellar pickup ions to a maximum energy limited by the size of the computational domain, with somewhat higher fluxes and maximal energies observed in the particle-in-cell simulation, likely due to differences in the cross-shock electric field arising from electron kinetic-scale effects. The higher fluxes may help address recent discrepancies noted between observations and large-scale hybrid simulations.
more »
« less
This content will become publicly available on January 20, 2027
A “Faux-Shock” Method for Hybrid Simulations of Astrophysical Shocks
Abstract We demonstrate a novel setup for hybrid particle-in-cell simulations designed to isolate the physics of the shock precursor over long time periods for significantly lower computational cost than previous methods. This is achieved using a “faux-shock” or shock-like boundary condition on one edge of our simulation domain such that particles that interact with the boundary either pass through it or are reflected off of it with a change in momentum that mimics scattering in the downstream. We show that our faux-shock setup reproduces the same fluid quantities and phase spaces as traditional shock simulations, including those which could otherwise only be done in 3D, with higher particle resolution and for reduced computational cost. While the method involves an assumed boundary condition, it nonetheless captures the essential physics of interest, establishing it as a reliable and efficient tool for future self-consistent studies of instabilities driven by cosmic rays in a shock upstream medium.
more »
« less
- PAR ID:
- 10669309
- Publisher / Repository:
- ApJ
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 997
- Issue:
- 2
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 144
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract Shock waves from underwater and air explosions are significant threats to surface and underwater vehicles and structures. Recent studies on the mechanical and thermal properties of various phase-separated elastomers indicate the possibility of applying these materials as a coating to mitigate shock-induced structural failures. To demonstrate this approach and investigate its efficacy, this paper presents a fluid-structure coupled computational model capable of predicting the dynamic response of air-backed bilayer (i.e. elastomer coating – metal substrate) structures submerged in water to hydrostatic and underwater explosion loads. The model couples a three-dimensional multiphase finite volume computational fluid dynamics model with a nonlinear finite element computational solid dynamics model using the FIVER (FInite Volume method with Exact multi-material Riemann solvers) method. The kinematic boundary condition at the fluid-structure interface is enforced using an embedded boundary method that is capable of handling large structural deformation and topological changes. The dynamic interface condition is enforced by formulating and solving local, one-dimensional fluid-solid Riemann problems, which is well-suited for transferring shock and impulsive loads. The capability of this computational model is demonstrated through a numerical investigation of hydrostatic and shock-induced collapse of aluminum tubes with polyurea coating on its inner surface. The thickness of the structure is resolved explicitly by the finite element mesh. The nonlinear material behavior of polyurea is accounted for using a hyper-viscoelastic constitutive model featuring a modified Mooney-Rivlin equation and a stress relaxation function in the form of prony series. Three numerical experiments are conducted to simulate and compare the collapse of the structure in different loading conditions, including a constant pressure, a fluid environment initially in hydrostatic equilibrium, and a two-phase fluid flow created by a near-field underwater explosion.more » « less
-
Abstract Understanding the conditions that enable particle acceleration at nonrelativistic collisionless shocks is essential to unveil the origin of cosmic rays. We employ 2D and 3D hybrid simulations (with kinetic ions and fluid electrons) to explore particle acceleration and magnetic field amplification in nonrelativistic perpendicular shocks, focusing on the role of shock drift acceleration and its dependence on the shock Mach number. We perform an analysis of the ion injection process and demonstrate why efficient acceleration is only observed in 3D. In particular, we show that ion injection critically depends on the “porosity” of the magnetic turbulence in the downstream region near the shock, a property describing how easily the postshock region allows particles to traverse it and return upstream without being trapped. This effect can only be properly captured in 3D. Additionally, we explore the impact of numerical resolution on ion energization, highlighting how resolving small-scale turbulence—on scales below the thermal ion gyroradius—is essential for accurately modeling particle injection. Overall, our results emphasize the necessity of high-resolution 3D simulations to capture the fundamental microphysics driving particle acceleration at perpendicular shocks.more » « less
-
ABSTRACT Spider pulsars are compact binary systems composed of a millisecond pulsar and a low-mass companion. Their X-ray emission – modulated on the orbital period – is interpreted as synchrotron radiation from high-energy electrons accelerated at the intrabinary shock. We perform global two-dimensional particle-in-cell simulations of the intrabinary shock, assuming that the shock wraps around the companion star. When the pulsar spin axis is nearly aligned with the orbital angular momentum, we find that the magnetic energy of the relativistic pulsar wind – composed of magnetic stripes of alternating field polarity – efficiently converts to particle energy at the intrabinary shock, via shock-driven reconnection. The highest energy particles accelerated by reconnection can stream ahead of the shock and be further accelerated by the upstream motional electric field. In the downstream, further energization is governed by stochastic interactions with the plasmoids/magnetic islands generated by reconnection. We also extend our earlier work by performing simulations that have a larger (and more realistic) companion size and a more strongly magnetized pulsar wind. We confirm that our first-principles synchrotron spectra and light curves are in good agreement with X-ray observations.more » « less
-
Abstract Fast radio bursts (FRBs) are extremely energetic radio transients, where some are generated in magnetar magnetospheres and winds. Despite a growing number of observations, their emission mechanisms remain elusive. It has recently been proposed that Alfvénic perturbations can convert into superluminalO-modes at magnetized shocks and propagate downstream as a radio signal. We validate this superluminal wave activation mechanism using pair-plasma theory and particle-in-cell simulations. Theory predicts two different downstream modes: nonpropagating Alfvénic perturbations and propagating superluminalO-modes. Superluminal wave activation occurs if the frequency of upstream perturbations in the shock frame exceeds the downstream plasma frequency. 1D particle-in-cell simulations confirm wavenumber and frequency jumps across the shock for upstream perturbations with frequencies well above the plasma frequency. Our simulations model both monochromatic upstream waves and broadband spectra with the downstream plasma frequency acting like a high-pass filter for superluminalO-modes. We discuss implications for FRB generation in relativistic magnetized winds.more » « less
