Scholarlink

Abstract <p>Numerical modeling of galactic cosmic rays (GCRs) penetration through the heliosphere to the vicinity of the Sun is considered. Galactic cosmic rays are charged particles with energies exceeding 10 MeV/nucl., originating from far beyond the boundaries of our Solar System. As they penetrate through the heliosphere—the region of space filled by the solar wind—they interact strongly with the interplanetary magnetic field. In this paper, we present numerical approaches to solving the so-called Parker transport equation for the isotropic velocity distribution function of GCRs. This equation includes a convective term, anisotropic diffusion, adiabatic cooling, and drifts. Additionally, the diffusion coefficient is spatially and energy-dependent, varying by several orders of magnitude. Our numerical approaches are based on the finite-difference method (Crank–Nicolson scheme) and the stochastic differential equations (SDE) method. The numerical methods were validated against a known analytical solution under simplified conditions. For the general problem formulation, which involves anisotropic diffusion and the Parker spiral interplanetary magnetic field configuration, we used the most efficient and flexible SDE method and compared the numerical results with the data from the works of Kota and Jokipii [<CitationRef CitationID="CR1">1</CitationRef>] and Burger [<CitationRef CitationID="CR2">2</CitationRef>]. Special attention was devoted to incorporating drift along the heliospheric current sheet in the model.</p>

错误：搜索内容不能为空，请输入英文关键词
错误：关键词超出字数限制，请精简
高级检索

Numerical Modeling of Galactic Cosmic Ray Modulation in the Heliosphere

  • 期刊论文

  • 发表日期：

    2026年1月1日

查看全文
  • D. A. Shestakov,
  • V. V. Izmodenov

摘要

Abstract

Numerical modeling of galactic cosmic rays (GCRs) penetration through the heliosphere to the vicinity of the Sun is considered. Galactic cosmic rays are charged particles with energies exceeding 10 MeV/nucl., originating from far beyond the boundaries of our Solar System. As they penetrate through the heliosphere—the region of space filled by the solar wind—they interact strongly with the interplanetary magnetic field. In this paper, we present numerical approaches to solving the so-called Parker transport equation for the isotropic velocity distribution function of GCRs. This equation includes a convective term, anisotropic diffusion, adiabatic cooling, and drifts. Additionally, the diffusion coefficient is spatially and energy-dependent, varying by several orders of magnitude. Our numerical approaches are based on the finite-difference method (Crank–Nicolson scheme) and the stochastic differential equations (SDE) method. The numerical methods were validated against a known analytical solution under simplified conditions. For the general problem formulation, which involves anisotropic diffusion and the Parker spiral interplanetary magnetic field configuration, we used the most efficient and flexible SDE method and compared the numerical results with the data from the works of Kota and Jokipii [1] and Burger [2]. Special attention was devoted to incorporating drift along the heliospheric current sheet in the model.