We describe the magnetohydrodynamics (MHD) code CRONOS, which has been used in astrophysics and space physics studies in recent years. CRONOS has been designed to be easily adaptable to the problem at hand, where the user can expand or exchange core modules or add new functionality to the code. This modularity comes about through its implementation using a C++ class structure. The core components of the code include solvers for both hydrodynamical (HD) and MHD problems. These problems are solved on different rectangular grids, which currently support Cartesian, spherical, and cylindrical coordinates. CRONOS uses a finite-volume description with different approximate Riemann solvers that can be chosen at runtime. Here, we describe the implementation of the code with a view toward its ongoing development. We illustrate the code's potential through several (M)HD test problems and some astrophysical applications.

This paper summarizes the results obtained by the team "Heliosheath Processes and the Structure of the Heliopause: Modeling Energetic Particles, Cosmic Rays, and Magnetic Fields" supported by the International Space Science Institute in Bern, Switzerland. We focus on the physical processes occurring in the outer heliosphere, especially at its boundary called the heliopause (HP), and in the LISM. The importance of magnetic field, charge exchange between atoms and ions, and solar cycle on the heliopause topology and observed heliocentric distances to different heliospheric discontinuities are discussed. It is shown that time-dependent boundary conditions are necessary to describe the heliospheric asymmetries detected by the Voyager spacecraft. We also discuss the structure of the HP, especially due to its instability and magnetic reconnection. It is demonstrated that the Rayleigh-Taylor instability of the nose of the HP creates consecutive layers of the interstellar and heliospheric plasma which are magnetically connected to different sources. This may be a possible explanation of abrupt changes in the galactic and anomalous cosmic ray fluxes observed by Voyager 1 when it was crossing the HP structure for a period of about one month in the summer of 2012. This paper also discusses the plausibility of fitting simulation results to a number of observational data sets obtained by in situ and remote measurements. The distribution of magnetic field in the vicinity of the HP is discussed in the context of Voyager measurements. We discuss the transport of energetic particles in the inner and outer heliosheath, concentrating on the anisotropic spatial diffusion diffusion tensor and the pitch-angle dependence of perpendicular diffusion and demonstrate that the latter can explain the observed pitch-angle anisotropies of both the anomalous and galactic cosmic rays in the outer heliosheath.

We extend a two-component model for the evolution of fluctuations in the solar wind plasma so that it is fully three-dimensional (3D) and also coupled self-consistently to the large-scale magnetohydrodynamic (MHD) equations describing the background solar wind. The two classes of fluctuations considered are a high-frequency parallel-propagating wave-like piece and a low-frequency quasi-two-dimensional component. For both components, the nonlinear dynamics is dominanted by quasi-perpendicular spectral cascades of energy. Driving of the fluctuations, by, for example, velocity shear and pickup ions, is included. Numerical solutions to the new model are obtained using the Cronos framework, and validated against previous simpler models. Comparing results from the new model with spacecraft measurements, we find improved agreement relative to earlier models that employ prescribed background solar wind fields. Finally, the new results for the wave-like and quasi-two-dimensional fluctuations are used to calculate ab initio diffusion mean free paths and drift lengthscales for the transport of cosmic rays in the turbulent solar wind.

The interaction between a supersonic stellar wind and a (super-)sonic interstellar wind has recently been viewed with new interest. We here first give an overview of the modeling, which includes the heliosphere as an example of a special astrosphere. Then we concentrate on the shock structures of fluid models, especially of hydrodynamic (HD) models. More involved models taking into account radiation transfer and magnetic fields are briefly sketched. Even the relatively simple HD models show a rich shock structure, which might be observable in some objects. We employ a single fluid model to study these complex shock structures, and compare the results obtained including heating and cooling with results obtained without these effects. Furthermore, we show that in the hypersonic case valuable information of the shock structure can be obtained from the Rankine-Hugoniot equations. We solved the Euler equations for the single fluid case and also for a case including cooling and heating. We also discuss the analytical Rankine-Hugoniot relations and their relevance to observations. We show that the only obtainable length scale is the termination shock distance. Moreover, the so-called thin shell approximation is usually not valid. We present the shock structure in the model that includes heating and cooling, which differs remarkably from that of a single fluid scenario in the region of the shocked interstellar medium. We find that the heating and cooling is mainly important in this region and is negligible in the regions dominated by the stellar wind beyond an inner boundary.

We present the implementation of turbulence transport equations in addition to the Reynolds-averaged MHD equations within the Cronos framework. The model is validated by comparisons with earlier findings before it is extended to be applicable to regions in the solar wind that are not highly super-Alfv\'enic. We find that the respective additional terms result in absolute normalized cross-helicity to decline more slowly, while a proper implementation of the mixing terms can even lead to increased cross-helicities in the inner heliosphere. The model extension allows to place the inner boundary of the simulations closer to the Sun, where we choose its location at 0.1 AU for future application to the Wang-Sheeley-Arge model. Here, we concentrate on effects on the turbulence evolution for transient events by injecting a coronal mass ejection (CME). We find that the steep gradients and shocks associated with these structures result in enhanced turbulence levels and reduced cross-helicity. Our results can now be used straightforwardly for studying the transport of charged energetic particles, where the elements of the diffusion tensor can now benefit from the self-consistently computed solar wind turbulence. Furthermore, we find that there is no strong back-reaction of the turbulence on the large-scale flow so that CME studies concentrating on the latter need not be extended to include turbulence transport effects.

Cosmic rays passing through large astrospheres can be efficiently cooled inside these "cavities" in the interstellar medium. Moreover, the energy spectra of these energetic particles are already modulated in front of the astrospherical bow shocks. We study the cosmic ray flux in and around lambda Cephei as an example for an astrosphere. The large-scale plasma flow is modeled hydrodynamically with radiative cooling. We studied the cosmic ray flux in a stellar wind cavity using a transport model based on stochastic differential equations. The required parameters, most importantly, the elements of the diffusion tensor, are based on the heliospheric parameters. The magnetic field required for the diffusion coefficients is calculated kinematically. We discuss the transport in an astrospheric scenario with varying parameters for the transport coefficients. We show that large stellar wind cavities can act as sinks for the galactic cosmic ray flux and thus can give rise to small-scale anisotropies in the direction to the observer. Small-scale cosmic ray anisotropies can naturally be explained by the modulation of cosmic ray spectra in huge stellar wind cavities.

The transport of energetic particles such as Cosmic Rays is governed by the properties of the plasma being traversed. While these properties are rather poorly known for galactic and interstellar plasmas due to the lack of in situ measurements, the heliospheric plasma environment has been probed by spacecraft for decades and provides a unique opportunity for testing transport theories. Of particular interest for the 3D heliospheric transport of energetic particles are structures such as corotating interaction regions (CIRs), which, due to strongly enhanced magnetic field strengths, turbulence, and associated shocks, can act as diffusion barriers on the one hand, but also as accelerators of low energy CRs on the other hand as well. In a two-fold series of papers we investigate these effects by modeling inner-heliospheric solar wind conditions with a numerical magnetohydrodynamic (MHD) setup (this paper), which will serve as an input to a transport code employing a stochastic differential equation (SDE) approach (second paper). In this first paper we present results from 3D MHD simulations with our code CRONOS: for validation purposes we use analytic boundary conditions and compare with similar work by Pizzo. For a more realistic modeling of solar wind conditions, boundary conditions derived from synoptic magnetograms via the Wang-Sheeley-Arge (WSA) model are utilized, where the potential field modeling is performed with a finite-difference approach (FDIPS) in contrast to the traditional spherical harmonics expansion often utilized in the WSA model. Our results are validated by comparing with multi-spacecraft data for ecliptical (STEREO-A/B) and out-of-ecliptic (Ulysses) regions.

Maps of the radial magnetic field at a heliocentric distance of ten solar radii are used as boundary conditions in the MHD code CRONOS to simulate a 3D inner-heliospheric solar wind emanating from the rotating Sun out to 1 AU. The input data for the magnetic field are the result of solar surface flux transport modelling using observational data of sunspot groups coupled with a current sheet source surface model. Amongst several advancements, this allows for higher angular resolution than that of comparable observational data from synoptic magnetograms. The required initial conditions for the other MHD quantities are obtained following an empirical approach using an inverse relation between flux tube expansion and radial solar wind speed. The computations are performed for representative solar minimum and maximum conditions, and the corresponding state of the solar wind up to the Earths orbit is obtained. After a successful comparison of the latter with observational data, they can be used to drive outer-heliospheric models.