For calculations in complex geometries boundary-fitted non-orthogonal curvlinear grids is usually used.