Shape optimization of a Printed circuit heat exchanger (PCHE) has been performed by using three-dimensional Reynolds-Averaged Navier-Stokes (3-D RANS) analysis and surrogate modeling techniques. The objective function is defined as a linear combination of effectiveness of the PCHE term and pressure drop in the cold channels of the PCHE. The cold channel angle and the ellipse aspect ratio of the cold channel are used as design variables for the optimization. Design points are selected through Latin-hypercube sampling. The optimal point is determined through surrogate-based optimization method which uses 3-D RANS analyses at design points. The results of three types of surrogate model are compared each other. The results of the optimizations indicate improved performance in friction loss but low performance in effectiveness than the reference shape.