The Drag Functional(Hydro dynamical Force acting on the boundary)is chosen as objective functional for shape optimization of Navier-Stokes boundary.Since conjugate gradient methods to compute optimization must do numerical differential for 3D Stress tensor and Gateaux derivative of solutions of Navier-Stokes equation with respect to the shape of boundary.Thus is a difficult and no efficiently problem.Our contributions are that all computation for conjugate gradient method for this kind of optimization do not need numerical differentiation for stress tensor and Gateaux derivative of solutions of Navier-Stokes equation with respect to the shape of boundary it is only to solve two dimensional boundary layer equations.This method is based on “A Finite Element Splitting Methods for Rotating Navier-Stokes Equations”.A well-known example of 3D rotating Navier-Stokes equations is to model the channel flows in a turbo machinery.Since the boundary geometry is very complex,it is difficulty to perform efficient numerical simulation.The method can be applied to the boundary geometry to establish a dimensional splitting method including one dimensional problem and a series of two dimensional problems on a two dimensional manifold which is called 2D-3D Navier-Stokes equations.This method does not need to generate 3D mesh and a bi-parallel algorithm will be developed.