Numerical methods for solution of the stochastic differential equations equivalent to the non-stationary Parkers transport equation

Wawrzyńczak, A.; Modzelewska, Renata; Kluczek, Mateusz

doi:10.1088/1742-6596/633/1/012058

Cited by 4 publications

(1 citation statement)

References 8 publications

(13 reference statements)

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Value of the distribution function is assumed as the statistical mean among the values of a local interstellar function for all simulated pseudoparticles. The more details on the solution of the set SDEs ( 2 ) can be found in [ 16 , 17 , 21 ].…”

Section: Applied Methodologymentioning

confidence: 99%

Approximate Bayesian Computation for Estimating Parameters of Data-Consistent Forbush Decrease Model

Wawrzyńczak

Kopka

2018

Entropy

View full text Add to dashboard Cite

Realistic modeling of complex physical phenomena is always quite a challenging task. The main problem usually concerns the uncertainties surrounding model input parameters, especially when not all information about a modeled phenomenon is known. In such cases, Approximate Bayesian Computation (ABC) methodology may be helpful. The ABC is based on a comparison of the model output data with the experimental data, to estimate the best set of input parameters of the particular model. In this paper, we present a framework applying the ABC methodology to estimate the parameters of the model of Forbush decrease (Fd) of the galactic cosmic ray intensity. The Fd is modeled by the numerical solution of the Fokker–Planck equation in five-dimensional space (three spatial variables, the time and particles energy). The most problematic in Fd modeling is the lack of detailed knowledge about the spatial and temporal profiles of the parameters responsible for the creation of the Fd. Among these parameters, the diffusion coefficient plays a central role. We employ the ABC Sequential Monte Carlo algorithm, scanning the space of the diffusion coefficient parameters within the region of the heliosphere where the Fd is created. Assessment of the correctness of the proposed parameters is done by comparing the model output data with the experimental data of the galactic cosmic ray intensity. The particular attention is put on the rigidity dependence of the rigidity spectrum exponent. The proposed framework is adopted to create the model of the Fd observed by the neutron monitors and ground muon telescope in November 2004.

show abstract

Section: Applied Methodologymentioning

confidence: 99%