This page provides information on research conducted in the Earth and Environmental Sciences Area of Berkeley Lab on diffusion through porous and fractured media. Currently, it includes the following:
- A PDF file of the report “A set of semianalytical solutions for parameter estimation in diffusion cell experiments” (by G. Moridis, Report No. LBNL-41857). A much shortened version of this report was published under the title “Semianalytical solutions for parameter estimation in diffusion cell experiments” in Water Resources Research, 35(6), 1729-1740, 1999.
- The set of FORTRAN 77 codes of the semianalytical solutions discussed in the report, as well as representative data sets. The codes provided here are:
- DIFF_TH.f: For the analysis of data from diffusion cell experiments. Provides output of boundary concentrations and mass balances over a desired time. Uses the De Hoog algorithm for the numerical inversion.
- DIFF_TS.f: For the analysis of data from diffusion cell experiments. Provides output of boundary concentrations and mass balances over a desired time. Uses the Stehfest algorithm for the numerical inversion.
- ttt: A sample data file, to be used by either DIFF_TH.f or DIFF_TS.f.
- DIFF_XH.f: For the study of concentration distribution in samples or barrier performance. Provides output of boundary concentrations and mass balances over a the sample/barrier thickness at a desired time. Uses the De Hoog algorithm for the numerical inversion.
- DIFF_XS.f: For the study of concentration distribution in samples or barrier performance. Provides output of boundary concentrations and mass balances over a the sample/barrier thickness at a desired time. Uses the Stehfest algorithm for the numerical inversion.
- xxx: A sample data file, to be used by either DIFF_XH.f or DIFF_XS.f.
Compilation and execution is straightforward. Supporting information can be found in the input data files, as well as in the codes. You are welcome to copy, change and use these files as you see fit. Please let me know by e-mail if you do so. I would be delighted to receive feedback from users of the code concerning performance and potential problems.
For more information please contact:
Dr. George Moridis
Hydrocarbon Resources Program Head
Geological Senior Scientist