BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Horst Fichtner

Horst Fichtner appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRphysics.space-phphysics.plasm-phastro-ph.HEastro-ph.EPphysics.flu-dynastro-ph.IMcond-mat.stat-mechmath.NAphysics.comp-ph

Catalog footprint

What is connected

25works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

On the properties of 0.11 keV to 344 MeV ion spectra in the inner heliosheath using regularized $κ$-distributions

The shape of the ion energy spectra plays a critical role toward determining the ion energetics, the acceleration mechanisms and the possible sources of different plasma and suprathermal ion populations. The determination of the exact shape of the total particle spectrum, provide the necessary means to address the inner heliosheath (IHS) dynamics. Apart from various modeling efforts, a direct fit to the measured ion spectra for an extended energy range of $\sim$0.11 to 344 MeV has not been performed to date. We use an extended set of combined 0.11-55 keV ENA measurements from the Interstellar Boundary Explorer (IBEX-Lo and IBEX-Hi) and Cassini/Ion and Neutral Camera (INCA), converted to protons, together with $\sim$28 keV to 344 MeV ion measurements from the Low Energy Charged Particle (LECP) and Cosmic Ray Subsystem (CRS) experiments on Voyager 2, over the declining phase of Solar Cyle 23 (SC23) and ascending phase Solar Cylce 24 (SC24) (2009-2016) to study the characteristics of the particle energy spectrum. We fit the 0.11 keV to 344 MeV composite spectra with a set of regularized isotropic $κ$-distribution functions (RKD) allowing the determination of the macroscopic physical properties. We demonstrate that the 2009-2012 spectrum that corresponds to the declining phase of SC23 is well fitted by three different RKDs, while the 2013-2016 spectrum, associated with the rise of SC24, can only be approximated with six different RKDs. Our results are generally consistent with shock accelerated particles that undergo additional acceleration inside the IHS. We identify a low energy transmitted population of particles, a suprathermal reflected population and a very high energy component that is modulated by GCRs. The 2013-2016 time period is most likely associated with a mixture of particles from SC23 and SC24, which is reflected by the need to employ six RDKs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Revisiting the Ulysses electron data with a triple fit of velocity distributions

Given their uniqueness, the Ulysses data can still provide us with valuable new clues about the properties of plasma populations in the solar wind, and, especially, about their variations with heliographic coordinates. We revisit the electron data reported by by the SWOOPS instrument on-board of the Ulysses spacecraft between 1990 to early 2008. These observations reveal velocity distributions out of thermal equilibrium, with anisotropies (e.g., parallel drifts or/and different temperatures, parallel and perpendicular to the background magnetic field), and quasi-thermal and suprathermal populations with different properties. We apply a 2D nonlinear least square fitting procedure, using the Levenberg-Marquardt algorithm, to simultaneously fit the velocity electron data (up to a few keV) with a triple model combining three distinct populations: the more central quasi-thermal core and suprathermal halo, and a second suprathermal population consisting mainly of the electron strahl (or beaming population with a major field-aligned drift). The recently introduced $κ$-cookbook is used to describe each component with the following anisotropic distribution functions (recipes): Maxwellian, regularized kappa-, and generalized kappa-distributions. The temperature anisotropy quantified by the best fits is considered as a case study of the main parameters characterizing electron populations. By comparison to the core, both suprathermal populations exhibit higher temperature anisotropies, which slightly increase with the energy of electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Towards a realistic evaluation of transport coefficients in non-equilibrium space plasmas

Recent studies have outlined the interest for the evaluation of transport coefficients in space plasmas, where the observed velocity distributions of plasma particles are conditioned not only by the binary collisions, e.g., at low energies, but also by the energisation of particles from their interaction with wave turbulence and fluctuations, generating the suprathermal Kappa-distributed populations. This paper provides a first estimate of the main transport coefficients based on regularised Kappa distributions (RKDs), which, unlike standard Kappa distributions (SKDs), enable macroscopic parameterisation without mathematical divergences or physical inconsistencies. All transport coefficients derived here, i.e., the diffusion and mobility coefficients, electric conductivity, thermoelectric coefficient and thermal conductivity, are finite and well defined for all values of $κ> 0$. Moreover, for low values of $κ$ (i.e., below the SKD poles), the transport coefficients can be orders of magnitudes higher than the corresponding Maxwellian limits, meaning that significant underestimations can be made if suprathermal electrons are ignored.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Turbulence in the outer heliosphere

The solar wind (SW) and local interstellar medium (LISM) are turbulent media. Their interaction is governed by complex physical processes and creates heliospheric regions with significantly different properties in terms of particle populations, bulk flow and turbulence. Our knowledge of the solar wind turbulence \nature and dynamics mostly relies on near-Earth and near-Sun observations, and has been increasingly improving in recent years due to the availability of a wealth of space missions, including multi-spacecraft missions. In contrast, the properties of turbulence in the outer heliosphere are still not completely understood. In situ observations by Voyager and New Horizons, and remote neutral atom measurements by IBEX strongly suggest that turbulence is one of the critical processes acting at the heliospheric interface. It is intimately connected to charge exchange processes responsible for the production of suprathermal ions and energetic neutral atoms. This paper reviews the observational evidence of turbulence in the distant SW and in the LISM, advances in modeling efforts, and open challenges.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A new class of discontinuous solar wind solutions

A new class of one-dimensional solar wind models is developed within the general polytropic, single-fluid hydrodynamic framework. The particular case of quasi-adiabatic radial expansion with a localized heating source is considered. We consider analytical solutions with continuous Mach number over the entire radial domain while allowing for jumps in the flow velocity, density, and temperature, provided that there exists an external source of energy in the vicinity of the critical point which supports such jumps in physical quantities. This is substantially distinct from both the standard Parker solar wind model and the original nozzle solutions, where such discontinuous solutions are not permissible. We obtain novel sample analytic solutions of the governing equations corresponding to both slow and fast wind.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Linear Dispersion Theory of Parallel Electromagnetic Modes for Regularized Kappa-Distributions

The velocity particle distributions measured in-situ in space plasmas deviate from Maxwellian (thermal) equilibrium, showing enhanced suprathermal tails which are well described by the standard Kappa-distribution (SKD). Despite its successful application, the SKD is frequently disputed due to a series of unphysical implications like diverging velocity moments, preventing a macroscopic description of the plasma. The regularized Kappa-distribution (RKD) has been introduced to overcome these limitations, but the dispersion properties of RKD-plasmas are not explored yet. In the present paper we compute the wavenumber dispersion of the frequency and damping or growth rates for the electromagnetic modes in plasmas characterized by the RKD. This task is accomplished by using the grid-based kinetic dispersion solver LEOPARD developed for arbitrary gyrotropic distributions [P. Astfalk and F. Jenko, J. Geophys. Res. 122, 89 (2017)]. By reproducing previous results obtained for the SKD and Maxwellian, we validate the functionality of the code. Furthermore, we apply the isotropic as well as the anisotropic RKDs to investigate stable electromagnetic electron-cyclotron (EMEC) and ion-cyclotron (EMIC) modes as well as temperature-anisotropy-driven instabilities, both for the case $T_\perp / T_\parallel > 1$ (EMEC and EMIC instabilities) and for the case $T_\perp / T_\parallel < 1$ (proton and electron firehose instabilities), where $\parallel$ and $\perp$ denote directions parallel and perpendicular to the local time-averaged magnetic field. Provided that the cutoff parameter $α$ is small enough, the results show that the RKDs reproduce the dispersion curves of the SKD plasmas at both qualitative and quantitative levels. For higher values, however, physically significant deviation occurs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the Diversity of M-Star Astrospheres and the Role of Galactic Cosmic Rays Within

With upcoming missions such as the James Webb Space Telescope (JWST), the European Extremely Large Telescope (ELT), and the Atmospheric Remote-sensing Infrared Exoplanet Large-survey (ARIEL), we soon will be on the verge of detecting and characterizing Earth-like exoplanetary atmospheres for the first time. These planets are most likely to be found around smaller and cooler K- and M-type stars. However, recent observations showed that their radiation environment might be much harsher than that of the Sun. Thus, the exoplanets are most likely exposed to an enhanced stellar radiation environment, which could affect their habitability, for example, in the form of a hazardous flux of energetic particles. Knowing the stellar radiation field, and being able to model the radiation exposure on the surface of a planet is crucial to assess its habitability. In this study, we present 3D magnetohydrodynamic (MHD)-based model efforts investigating M-stars, focusing on V374 Peg, Proxima Centauri, and LHS 1140, chosen because of their diverse astrospheric quantities. We show that V374 Peg has a much larger astrosphere (ASP) than our Sun, while Proxima Centauri and LHS 1140 most likely have ASPs comparable or even much smaller than the heliosphere, respectively. Based on a 1D transport model, for the first time, we provide numerical estimates of the modulation of Galactic cosmic rays (GCRs) within the three ASPs. We show that the impact of GCRs on the Earth-like exoplanets Proxima Centauri b and LHS 1140 b cannot be neglected in the context of exoplanetary habitability.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The $κ$-cookbook: a novel generalizing approach to unify $κ$-like distributions for plasma particle modeling

In the literature different so-called $κ$-distribution functions are discussed to fit and model the velocity (or energy) distributions of solar wind species, pickup ions or magnetospheric particles. Here we introduce a generalized (isotropic) $κ$-distribution as a "cookbook", which admits as special cases, or "recipes", all the other known versions of $κ$-models. A detailed analysis of the generalized distribution function is performed, providing general analytical expressions for the velocity moments, Debye length, and entropy, and pointing out a series of general requirements that plasma distribution functions should satisfy. From a contrasting analysis of the recipes found in the literature, we show that all of them lead to almost the same macroscopic parameters with a small standard deviation between them. However, one of these recipes called the regularized $κ$-distribution provides a functional alternative for macroscopic parameterization without any constraint for the power-law exponent $κ$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An empirical model of Energetic Neutral Atom imaging of the heliosphere and its implications for future heliospheric missions at great heliocentric distances

Several concepts for heliospheric missions operating at heliocentric distances far beyond Earth orbit are currently investigated by the scientific community. The mission concept of the Interstellar Probe (McNutt et al. 2018), e.g., aims at reaching a distance of 1000 au away from the Sun within this century. This would allow the coming generation to obtain a global view of our heliosphere from an outside vantage point by measuring the Energetic Neutral Atoms (ENAs) originating from the various plasma regions. It would also allow for direct sampling of unperturbed interstellar medium, and for many observation opportunities beyond heliospheric science, such as visits to Kuiper Belt Objects, a comprehensive view on the interplanetary dust populations, and infrared astronomy free from the foreground emission of the Zodiacal cloud. In this study, we present a simple empirical model of ENAs from the heliosphere and derive basic requirements for ENA instrumentation onboard a spacecraft at great heliocentric distances. We consider the full energy range of heliospheric ENAs from 10 eV to 100 keV because each part of the energy spectrum has its own merits for heliospheric science. To cover the full ENA energy range, two or three different ENA instruments are needed. Thanks to parallax observations, some insights about the nature of the IBEX Ribbon and the dimensions of the heliosphere can already be gained by ENA imaging from a few au heliocentric distance. To directly reveal the global shape of the heliosphere, measurements from outside the heliosphere are, of course, the best option.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A generalized two-component model of solar wind turbulence and ab initio diffusion mean free paths and drift lengthscales of cosmic rays

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Generalized multi-polytropic Rankine-Hugoniot relations and the entropy condition

The study aims at a derivation of generalized \RH relations, especially that for the entropy, for the case of different upstream/downstream polytropic indices and their implications. We discuss the solar/stellar wind interaction with the interstellar medium for different polytropic indices and concentrate on the case when the polytropic index changes across hydrodynamical shocks. We use first a numerical mono-fluid approach with constant polytropic index in the entire integration region to show the influence of the polytropic index on the thickness of the helio-/astrosheath and on the compression ratio. Second, the Rankine-Hugonoit relations for a polytropic index changing across a shock are derived analytically, particularly including a new form of the entropy condition. In application to the/an helio-/astrosphere, we find that the size of the helio-/astrosheath as function of the polytropic index decreases in a mono-fluid model for indices less than $γ=5/3$ and increases for higher ones and vice versa for the compression ratio. Furthermore, we demonstrate that changing polytropic indices across a shock are physically allowed only for sufficiently high Mach numbers and that in the hypersonic limit the compression ratio depends only on the downstream polytropic index, while the ratios of the temperature and pressure as well as the entropy difference depend on both, the upstream and downstream polytropic indices.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Toward more realistic analytic models of the heliotail: Incorporating magnetic flattening via distortion flows

Both physical arguments and simulations of the global heliosphere indicate that the tailward heliopause is flattened considerably in the direction perpendicular to both the incoming flow and the large-scale interstellar magnetic field. Despite this fact, all of the existing global analytical models of the outer heliosheath's magnetic field assume a circular cross section of the heliotail. To eliminate this inconsistency, we introduce a mathematical procedure by which any analytically or numerically given magnetic field can be deformed in such a way that the cross sections along the heliotail axis attain freely prescribed, spatially dependent values for their total area and aspect ratio. The distorting transformation of this method honors both the solenoidality condition and the stationary induction equation with respect to an accompanying flow field, provided that both constraints were already satisfied for the original magnetic and flow fields prior to the transformation. In order to obtain realistic values for the above parameters, we present the first quantitative analysis of the heliotail's overall distortion as seen in state-of-the-art three-dimensional hybrid MHD-kinetic simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An exact analytical solution for the interstellar magnetic field in the vicinity of the heliosphere

An analytical representation of the interstellar magnetic field in the vicinity of the heliosphere is derived. The three-dimensional field structure close to the heliopause is calculated as a solution of the induction equation under the assumption that it is frozen into a prescribed plasma flow resembling the characteristic interaction of the solar wind with the local interstellar medium. The usefulness of this analytical solution as an approximation to self-consistent magnetic field configurations obtained numerically from the full MHD equations is illustrated by quantitative comparisons.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmic rays in astrospheres

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Implementing turbulence transport in the CRONOS framework and application to the propagation of CMEs

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énic. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the Geometry of the IBEX Ribbon

The Energetic Neutral Atom (ENA) full-sky maps obtained with the Interstellar Boundary Explorer (IBEX) show an unexpected bright narrow band of increased intensity. This so-called ENA ribbon results from charge exchange of interstellar neutral atoms with protons in the outer heliosphere or beyond. Amongst other hypotheses it has been argued that this ribbon may be related to a neutral density enhancement, or H-wave, in the local interstellar medium. Here we quantitatively demonstrate, on the basis of an analytical model of the principal large-scale heliospheric structure, that this scenario for the ribbon formation leads to results that are fully consistent with the observed location of the ribbon in the full-sky maps at all energies detected with high-energy sensor IBEX-Hi.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Parameter estimation of superdiffusive motion of energetic particles upstream of heliospheric shocks

In-situ spacecraft observations recently suggested that the transport of energetic particles accelerated at heliospheric shocks can be anomalous, i.e. the mean square displacement can grow non-linearly in time. In particular, a new analysis technique has permitted the study of particle transport properties from energetic particle time profiles upstream of interplanetary shocks. Indeed, the time/spatial power laws of the differential intensity upstream of several shocks are indicative of superdiffusion. A complete determination of the key parameters of superdiffusive transport comprises the power-law index, the superdiffusion coefficient, the related transition scale at which the energetic particle profiles turn to decay as power laws, and the energy spectral index of the shock accelerated particles. Assuming large-scale spatial homogeneity of the background plasma, the power-law behaviour can been derived from both a (microscopic) propagator formalism and a (macroscopic) fractional transport equation. We compare the two approaches and find a relation between the diffusion coefficients used in the two formalisms. Based on the assumption of superdiffusive transport, we quantitatively derive these parameters by studying energetic particle profiles observed by the Ulysses and Voyager 2 spacecraft upstream of shocks in the heliosphere, for which a superdiffusive particle transport has previously been observed. Further, we have jointly studied the electron energy spectra, comparing the values of the spectral indices observed with those predicted by the standard diffusive shock acceleration theory and by a model based on superdiffusive transport. For a number of interplanetary shocks and for the solar wind termination shock, for the first time we obtain the anomalous diffusion constants and the scale at which the probability of particle free paths changes to a power-law...

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Shock structures of astrospheres

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmic Ray Transport in Heliospheric Magnetic Structures: I. Modeling Background Solar Wind Using the CRONOS MHD Code

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The return of the bow shock

Recently it has been discussed whether a bow shock ahead of the heliospheric stagnation region does exist or not. This discussion was triggered by measurements indicating that the Alfvén speed and that of fast magnetosonic waves are higher than the flow speed of the local interstellar medium (LISM) relative to the heliosphere and resulted in the conclusion that there might exist either a bow wave or a slow magnetosonic shock. We demonstrate here that including the He$^{+}$ component of the LISM yields both an Alfvén and fast magnetosonic wave speed lower than the LISM flow speed. Consequently, the scenario of a bow shock in front of the heliosphere as modelled in numerous simulations of the interaction of the solar wind with the LISM remains valid.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The space-fractional diffusion-advection equation: Analytical solutions and critical assessment of numerical solutions

The present work provides a critical assessment of numerical solutions of the space-fractional diffusion-advection equation, which is of high significance for applications in various natural sciences. In view of the fact that, in contrast to the case of normal (Gaussian) diffusion, no standard methods and corresponding numerical codes for anomalous diffusion problems have been established yet, it is of importance to critically assess the accuracy and performance of existing approaches. Three numerical methods, namely a finite-difference method, the so-called matrix transfer technique, and a Monte-Carlo method based on the solution of stochastic differential equations, are analyzed and compared by applying them to three selected test problems for which analytical or semi-analytical solutions were known or are newly derived. The accuracy and performance differences are critically discussed with the result that the use of stochastic differential equations appears to be advantageous.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Generalized Diffusion Tensor for Fully Anisotropic Diffusion of Energetic Particles in the Heliospheric Magnetic Field

The spatial diffusion of cosmic rays in turbulent magnetic fields can, in the most general case, be fully anisotropic, i.e. one has to distinguish three diffusion axes in a local, field-aligned frame. We reexamine the transformation for the diffusion tensor from this local to a global frame, in which the Parker transport equation for energetic particles is usually formulated and solved. Particularly, we generalize the transformation formulas to allow for an explicit choice of two principal local perpendicular diffusion axes. This generalization includes the 'traditional' diffusion tensor in the special case of isotropic perpendicular diffusion. For the local frame, we motivate the choice of the Frenet-Serret trihedron which is related to the intrinsic magnetic field geometry. We directly compare the old and the new tensor elements for two heliospheric magnetic field configurations, namely the hybrid Fisk and the Parker field. Subsequently, we examine the significance of the different formulations for the diffusion tensor in a standard 3D model for the modulation of galactic protons. For this we utilize a numerical code to evaluate a system of stochastic differential equations equivalent to the Parker transport equation and present the resulting modulated spectra. The computed differential fluxes based on the new tensor formulation deviate from those obtained with the 'traditional' one (only valid for isotropic perpendicular diffusion) by up to 60% for energies below a few hundred MeV depending on heliocentric distance.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A novel code for numerical 3-D MHD studies of CME expansion

A recent third-order, essentially non-oscillatory central scheme to advance the equations of single-fluid magnetohydrodynamics (MHD) in time has been implemented into a new numerical code. This code operates on a 3-D Cartesian, non-staggered grid, and is able to handle shock-like gradients without producing spurious oscillations. To demonstrate the suitability of our code for the simulation of coronal mass ejections (CMEs) and similar heliospheric transients, we present selected results from test cases and perform studies of the solar wind expansion during phases of minimum solar activity. We can demonstrate convergence of the system into a stable Parker-like steady state for both hydrodynamic and MHD winds. The model is subsequently applied to expansion studies of CME-like plasma bubbles, and their evolution is monitored until a stationary state similar to the initial one is achieved. In spite of the model's (current) simplicity, we can confirm the CME's nearly self-similar evolution close to the Sun, thus highlighting the importance of detailed modelling especially at small heliospheric radii. Additionally, alternative methods to implement boundary conditions at the coronal base, as well as strategies to ensure a solenoidal magnetic field, are discussed and evaluated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Anisotropic diffusion of galactic cosmic ray protons and their steady-state azimuthal distribution

Galactic transport models for cosmic rays involve the diffusive motion of these particles in the interstellar medium. Due to the large-scale structured galactic magnetic field this diffusion is anisotropic with respect to the local field direction. We included this transport effect along with continuous loss processes in a quantitative model of galactic propagation for cosmic ray protons which is based on stochastic differential equations. We calculated energy spectra at different positions along the Sun's galactic orbit and compared them to the isotropic diffusion case. The results show that a larger amplitude of variation as well as different spectral shapes are obtained in the introduced anisotropic diffusion scenario and emphasize the need for accurate galactic magnetic field models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the low-frequency boundary of Sun-generated MHD turbulence in the slow solar wind

New aspects of the slow solar wind turbulent heating and acceleration are investigated. A physical meaning of the lower boundary of the Alfvén wave turbulent spectra in the solar atmosphere and the solar wind is studied and the significance of this natural parameter is demonstrated. Via an analytical and quantitative treatment of the problem we show that a truncation of the wave spectra from the lower frequency side, which is a consequence of the solar magnetic field structure and its cyclic changes, results in a significant reduction of the heat production and acceleration rates. An appropriate analysis is presented regarding the link of the considered problem with existing observational data and slow solar wind initiation scenarios.

0 citations0 reviewsNo code linkedOpen access