Abstract
Three dimensional computation process by the Boundary Integral Equation Method has been developed for flows in an unconfined aquifer. An availability of the computation process presented here was applied to a duct flow and two dimensional groundwater flow through a dam, and was confirmed by the comparison of the theoretical values with computed ones concerning the potential, flux and groundwater surface profile. Finally, a transient three dimensional flow pattern arround an underground dam was shown as a typical computation example