Correction of Curvilinear Grids for Numerical Solution of Groundwater Flow Problems

A new variational method is proposed for grid generation in numerical solution of groundwater flow problems requiring imposition of constraints on the space step. The proposed procedure utilizes the penalty function method for the solution of an optimization problem with inequality constraints. The efficiency of the method is illustrated by the generation of solution-dependent grids, which allow such constraints to be satisfied for complex-geometry domains.