Numerical Simulations of HH 555

We present 3D gasdynamic simulations of the Herbig Haro object HH 555. HH 555 is a bipolar jet emerging from the tip of an elephant trunk entering the Pelican Nebula from the adjacent molecular cloud. Both beams of HH 555 are curved away from the center of the H II region. This indicates that they are being deflected by a side-wind probably coming from a star located inside the nebula or by the expansion of the nebula itself. HH 555 is most likely an irradiated jet emerging from a highly embedded protostar, which has not yet been detected. In our simulations we vary the incident photon flux, which in one of our models is equal to the flux coming from a star 1 pc away emitting 5x10^48 ionizing (i. e., with energies above the H Lyman limit) photons per second. An external, plane-parallel flow (a ``side-wind'') is coming from the same direction as the photoionizing flux. We have made four simulations, decreasing the photon flux by a factor of 10 in each simulation. We discuss the properties of the flow and we compute Halpha emission maps (integrated along lines of sight). We show that the level of the incident photon flux has an important influence on the shape and visibility of the jet. If the flux is very high, it causes a strong evaporation of the neutral clump, producing a photoevaporated wind traveling in the direction opposite to the incident flow. The interaction of the two flows creates a double shock ``working surface'' around the clump protecting it and the jet from the external flow. The jet only starts to curve when it penetrates through the working surface.Comment: 14 pages, 4 figures, accepted by Ap

Traveling foreshocks and transient foreshock phenomena

We use the multispacecraft capabilities of the Cluster and Time History of Events and Macroscale Interactions during Substorms (THEMIS) missions to show that two types of foreshock may be detected in spacecraft data. One is the global foreshock that appears upstream of the Earth's quasi-parallel bow shock under steady or variable interplanetary magnetic field. Another type is a traveling foreshock that is bounded by two rotational discontinuities in the interplanetary magnetic field and propagates along the bow shock. Foreshock compressional boundaries are found at the edges of both types of foreshock. We show that isolated foreshock cavities are a subset of the traveling foreshocks that form when two bounding rotational discontinuities are so close that the ultralow-frequency waves do not develop in the region between them. We also report observations of a spontaneous hot flow anomaly inside a traveling foreshock. This means that other phenomena, such as foreshock cavitons, may also exist inside this type of foreshock. In the second part of this work we present statistical properties of phenomena related to the foreshock, namely, foreshock cavities, cavitons, spontaneous hot flow anomalies, and foreshock compressional boundaries. We show that spontaneous hot flow anomalies are the most depleted transient structures in terms of the B field and plasma density inside them and that the foreshock compressional boundaries and foreshock cavities are closely related structures

Transmission of foreshock waves through Earth’s bow shock

The Earth's magnetosphere and its bow shock, which is formed by the interaction of the supersonic solar wind with the terrestrial magnetic field, constitute a rich natural laboratory enabling in situ investigations of universal plasma processes. Under suitable interplanetary magnetic field conditions, a foreshock with intense wave activity forms upstream of the bow shock. So-called 30 s waves, named after their typical period at Earth, are the dominant wave mode in the foreshock and play an important role in modulating the shape of the shock front and affect particle reflection at the shock. These waves are also observed inside the magnetosphere and down to the Earth's surface, but how they are transmitted through the bow shock remains unknown. By combining state-of-the-art global numerical simulations and spacecraft observations, we demonstrate that the interaction of foreshock waves with the shock generates earthward-propagating, fast-mode waves, which reach the magnetosphere. These findings give crucial insight into the interaction of waves with collisionless shocks in general and their impact on the downstream medium.Peer reviewe

Backstreaming ions at a high Mach number interplanetary shock: Solar Orbiter measurements during the nominal mission phase

Solar Orbiter, a mission developed by the European Space Agency, explores in situ plasma across the inner heliosphere while providing remote-sensing observations of the Sun. Our study examines particle observations for the 30 October 2021 shock. The particles provide clear evidence of ion reflection up to several minutes upstream of the shock. Additionally, the magnetic and electric field observations contain complex electromagnetic structures near the shock, and we aim to investigate how they are connected to ion dynamics. The main goal of this study is to advance our understanding of the complex coupling between particles and the shock structure in high Mach number regimes of interplanetary shocks. We used observations of magnetic and electric fields, probe-spacecraft potential, and thermal and energetic particles to characterize the structure of the shock front and particle dynamics. Furthermore, ion velocity distribution functions were used to study reflected ions and their coupling to the shock. To determine shock parameters and study waves, we used several methods, including cold plasma theory, singular-value decomposition, minimum variance analysis, and shock Rankine-Hugoniot relations. To support the analysis and interpretation of the experimental data, test-particle analysis, and hybrid particle in-cell simulations were used. The ion velocity distribution functions show clear evidence of particle reflection in the form of backstreaming ions several minutes upstream. The shock structure has complex features at the ramp and whistler precursors. The backstreaming ions may be modulated by the complex shock structure, and the whistler waves are likely driven by gyrating ions in the foot. Supra-thermal ions up to 20 keV were observed, but shock-accelerated particles with energies above this were not

EMISSION LINE RATIOS FROM VARIABLE VELOCITY JET MODELS

RESUMEN En este artículo, presentamos una red de simulaciones numéricas axisimétricas de jets con velocidad de eyección variable. En estos modelos, suponemos que los jets son eyectados con una velocidad que varía sinusoidalmente, y con una densidad constante. La red de modelos entonces cubre un intervalo de diferentes amplitudes y períodos de la variabilidad de velocidad de eyección. Como las simulaciones incluyen un tratamiento de ionización fuera de equilibrio del gas, podemos hacer predicciones de la emisión en un conjunto de diferentes líneas espectrales. Así, obtenemos razones de las líneas de emisión para los nudos sucesivos a lo largo de los jets (que coresponden a las "superficies de trabajo internas" formadas como resultado de la variabilidad de la velocidad de eyección), las que se pueden comparar con observaciones de las cadenas de nudos a lo largo de jets HH. ABSTRACT In this paper we present a grid of axisymmetric numerical simulations of variable ejection velocity jets. In these models we assume that the jets are ejected with a sinusoidally varying ejection velocity and a time-independent ejection density. The grid of models then spans a range of different velocity variability amplitudes and periods. Because the simulations include a treatment of the non-equilibrium ionization state of the gas, we are able to make predictions of the emission in a set of different emission lines. In this way, we obtain predicted emission line ratios for the successive knots along the jets (which correspond to the "internal working surfaces" formed as a result of the ejection velocity variability), which can be compared directly with observations of the chains of knots along HH jets

SUPRATHERMAL ELECTRON STRAHL WIDTHS IN THE PRESENCE OF NARROW-BAND WHISTLER WAVES IN THE SOLAR WIND

International audienc

High resolution simulations of a variable HH jet

Context.In many papers, the flows in Herbig-Haro (HH) jets have been modeled as collimated outflows with a time-dependent ejection. In particular, a supersonic variability of the ejection velocity leads to the production of “internal working surfaces” which (for appropriate forms of the time-variability) can produce emitting knots that resemble the chains of knots observed along HH jets. Aims.In this paper, we present axisymmetric simulations of an “internal working surface” in a radiative jet (produced by an ejection velocity variability). We concentrate on a given parameter set (i.e., on a jet with a constante ejection density, and a sinusoidal velocity variability with a 20 yr period and a 40 km s-1 half-amplitude), and carry out a study of the behaviour of the solution for increasing numerical resolutions. Methods.In our simulations, we solve the gasdynamic equations together with a 17-species atomic/ionic network, and we are therefore able to compute emission coefficients for different emission lines. Results.We compute 3 adaptive grid simulations, with 20, 163 and 1310 grid points (at the highest grid resolution) across the initial jet radius. From these simulations we see that successively more complex structures are obtained for increasing numerical resolutions. Such an effect is seen in the stratifications of the flow variables as well as in the predicted emission line intensity maps. Conclusions.We find that while the detailed structure of an internal working surface depends on resolution, the predicted emission line luminosities (integrated over the volume of the working surface) are surprisingly stable. This is definitely good news for the future computation of predictions from radiative jet models for carrying out comparisons with observations of HH objects