A numerical procedure for optimizing the shape of three-dimensional sedimentation tank is presented to maximize its sedimentation efficiency. The response surface based optimization is used as an optimization technique with Reynolds-averaged Navier-Stokes analysis for multi-phase flow. Standard $k-{epsilon}$ model is used as a turbulence closure. Three design variables such as, tank height to center feed wall diameter ratio, blockage ratio of center feed wall and angle of distributor are chosen as design variables. Sedimentation efficiency is defined as an objective function. Full-factorial method is used to determine the training points as a means of design of experiment. Sensitivity of each design variable on the objective function has been evaluated. And, optimal values of the design variables have been obtained.