In this paper, a mathematical programming technique for numerical simulatio
n of unconfined flow through porous media is presented. Considering the ori
ginal free boundary problem as a shape optimization problem, we perform bou
ndary elements discretization. Taking the state variable and free boundary
variable as independent variables, we treat the discretized problem as nonl
inear mathematical program and apply interior point algorithm to solve it.
This simple, accurate and computationaly efficient technique can be easily
applied to 2D real size problems and extended to 3D problems. Numerical res
ults for an illustrative 2D test problem of an earth dam are discussed. (C)
2001 Elsevier Science Ltd. All rights reserved.