Migpore, a code package for the estimation of migration of radioactive species in a porous medium

In the event of an accidental leakage of high level radioactive waste buried deep in repositories surrounded by rock, the build up of the concentration of the radioactive species within the rock needs to be assessed. Towards this, we follow the model of Chen and Li and provide a numerical code to solve the relevant partial differential equations using a compact finite difference scheme.Program summaryProgram title: MigporeCatalogue identifier: AEQI_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEQI_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 5396No. of bytes in distributed program, including test data, etc.: 26483Distribution format: tar.gzProgramming language: Fortran 77.Computer: PC under Linux or Windows.Operating system: Ubuntu 10.04 (Kernel version 2.6.32-47-generic), Windows-XP.Classification: 4.3, 4.12, 13.External routines:We included the following Lapack subroutines in migpor.for program.: DCOPY, DGBTF2, DGBTRF, DGBTRS, DGEMM, DGEMV, DGER, DLASWP, DCAL, DSWAP, DTBSV, DTRSM, IDAMAX, IEEECK, ILAENV, IPARMQ, LSAME, XERBLA.Nature of problem:We consider the case of an accidental leakage of high level radioactive waste into the surrounding rock (a porous medium).Solution method:We follow the model of Chen and Li to estimate the buildup of radioactivity due to migration from the source. This amounts to the solution of two coupled partial differential equations. The numerical algorithm that we provide solves these equations by a higher order compact finite difference scheme.Running time:The Migpore code package needs approximately 50 min to run the test inputs. For other cases the run time will depend on the species, the distance and time step.