Numerical calculation of thermal convection in porous media in two dimensions with applications to nuclear repository design

A computer code HP2D was developed to calculate the effect. Transient and steady state heat and fluid flow in porous media are calculated. The code allows for anisotropy of hydraulic conductivities, temperature dependent fluid density and viscosity. A heat source and a variety of temperature and flow boundary conditions can be easily employed. The code solves two coupled governing equations. The time dependent heat transport is solved by an alternating direction implicit technique; and the steady state flow of water is solved by a Gauss-Seidel Poisson solver. A sample calculation of thermal convection about a repository is given.