Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aeqi_v1_0.tar.gz(25 Kbytes)
Manuscript Title: Migpore, a code package for the estimation of migration of radioactive species in a porous medium
Authors: Soubhadra Sen, N. Mohankumar
Program title: Migpore
Catalogue identifier: AEQI_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 185(2014)302
Programming language: Fortran 77.
Computer: PC under Linux or Windows.
Operating system: Ubuntu 10.04 (Kernel version 2.6.32-47-generic), Windows-XP.
Keywords: Radioactivity migration, diffusion, porous matrix, parallel fractures and repository.
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 migpor code package needs approximately 50 minutes to run the test inputs. For other cases the run time will depend on the species, the distance and time step.