We propose that the high energy Cosmic Ray particles around the spectral turn-down commonly called the {\it knee} and up to the upturn, commonly called the {\it ankle}, mostly come from Blue Super Giant star explosions. At the upturn, i.e. the {\it ankle}, Cosmic Rays probably switch to another source class, most likely extragalactic sources. To show this we recently compiled a set of Radio Supernova data to list the magnetic field, shock speed and radius scale (Biermann et al. 2018) \cite{Biermann18}. Using particle acceleration theory at shocks, those numbers can be transformed into characteristic {\it knee} and {\it ankle} energies. Without adjusting any free parameters both of these observed energies are directly indicated by the supernova data. We now proceed to the next step in the argument, and use the Supernova Remnant data of the starburst galaxy M82. Assuming that they are Blue Supergiant star explosions, the shock will race to their outer edge with a magnetic field that follows ${B (r) \, r \, \sim \, const.}$. We argue that the shock runs through the entire magnetic plasma wind region at full speed all the way out to the wind-shell, which is of order parsec scale. The speed is observed to be $\sim \, 0.1 \, c$ at about ${10^{16} \, {\rm cm}}$ radius in the plasma wind. This demonstrates how Blue Supergiant star explosions can provide the Cosmic Ray particles across the {\it knee} and up to the {\it ankle} energy range. The data from the CREAM (Cosmic Ray Energetics and Mass Experiment) mission will test this cosmic ray concept which is reasonably well grounded in two independent radio supernova data sets. The next step in developing our understanding is to obtain accurate Cosmic Ray data near to the {\it knee}, and use unstable isotopes of Cosmic Ray nuclei at high energy to probe the "piston" driving the explosion. We plan to combine these data with the physics of the budding black hole which is probably forming in each of these stars to learn more.

The Particle-In-Cell (PIC) method has been developed in order to investigate microscopic phenomena, and with the advances of computing power, newly developed codes have been used for several fields such as astrophysical, magnetospheric, and solar plasmas. PIC applications have grown extensively with large computing powers available on supercomputers such as Pleiades and Blue Waters in the US. For astrophysical plasma research PIC methods have been utilized for several topics such as reconnection, pulsar dynamics, non-relativistic shocks, relativistic shocks, relativistic jets, etc. PIC simulations of relativistic jets have been reviewed with the emphasis on the physics involved in the simulations. This review summarizes PIC simulations, starting with the Weibel instability in slab models of jets, and then focuses on global jet evolution in helical magnetic field geometry. In particular we address kinetic Kelvin-Helmholtz instabilities and mushroom instabilities.

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron–proton (e^- - p⁺) and electron–positron (e^±) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the e^- - p⁺ jet simulation recollimation-like instability occurs and jet electrons are strongly perturbed. In the e^± jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.