A new efficient method of generating solution-adaptive grids for cascade problems is introduced. This scheme uses Laplace equations, which are then transformed by using stretching functions to the final computational domain so that the generated grid can be clustered in the desired regions. Thus, the resulting generating equations are linear and uncoupled. To adapt the Mid to solutions, the control functions are chosen to depend upon the curvature and the gradient of the solution at each grid point and the grid spacing is controlled by these values. The scheme is then tested on model problems in two dimensions. Lastly, this solution-adaptive scheme is employed for transonic flow calculations.