Milne, a routine for the numerical solution of Milne's problem

Published: 01-11-2010| Version 1 | DOI: 10.17632/5wvgtcmytx.1
Ajay Rawat,
N. Mohankumar


Abstract The routine Milne provides accurate numerical values for the classical Milne's problem of neutron transport for the planar one speed and isotropic scattering case. The solution is based on the Case eigen-function formalism. The relevant X functions are evaluated accurately by the Double Exponential quadrature. The calculated quantities are the extrapolation distance and the scalar and the angular fluxes. Also, the H function needed in astrophysical calculations is evaluated as a byproduct. Title of program: Milne Catalogue Id: AEGS_v1_0 Nature of problem The X functions are integral expressions. The convergence of these regular and Cauchy Principal Value integrals are impaired by the singularities of the integrand in the complex plane. The DE quadrature scheme tackles these singularities in a robust manner compared to the standard Gauss quadrature. Versions of this program held in the CPC repository in Mendeley Data AEGS_v1_0; Milne; 10.1016/j.cpc.2010.07.023 This program has been imported from the CPC Program Library held at Queen's University Belfast (1969-2018)