The proposed numerical procedure for static open boundary problems combines
fast adaptive multigrid methods with spatial mappings of the infinite regi
on to a finite computational domain. A rigorous analysis of the finite elem
ent solution error for ballooning and for spatial mappings is presented. Th
e method is applied to geophysical problems of de logging.