Aerodynamic shape design using stochastic optimization methods, such as simulated annealing method to optimize objective functions evaluated by modern state-of-the-art computational uid dynamics solvers, normally requires enormous computation time to search for the global optimal design. Aerodynamic shape optimization of internal ow systems is studied using Euler/Navier-Stokes solvers and parallel simulated annealing algorithm, which is implemented on parallel computing platforms. A variety of inverse and direct design of internal ow systems are carried out to examine the ef ciency and speedup of the parallel simulated annealing algorithms. The results demonstrate that parallel simulated annealing can be a feasible global optimizer for aerodynamic shape design resulting in considerable reductions in wall-clock time on multiple processors.