An effective method for performing the thermal optimization of stationary and rotating MCM disks with an unconfined round-jet impingement under space limitation constraint has been successfully developed. The design variables of stationary and rotating MCM disks with an unconfined round-jet impingement include: the ratio of jet separation distance to nozzle diameter (H/d), steady-state Grashof number (Grs), jet Reynolds number (Rej), rotational Reynolds number (Rer). The total experimental cases for stationary and rotating MCM disks are statistically designed by the Central Composite Design (CCD) method. In addition, a sensitivity analysis, the so-called ANOVA, for the design factors has been performed. In the stationary MCM disk with an unconfined round-jet impingement, the contribution percentage of jet Reynolds number on the thermal performance is 95.86%. The effect of jet Reynolds numbers on chip temperature distribution is more significant than that of the H/d ratio and steady-state Grashof number. In rotating MCM disk with an unconfined round-jet impingement, the effect of jet Reynolds number, which has the contribution percentage of 91.81%, dominates the thermal performance. Furthermore, the comparisons between the predictions by using the quadratic Response Surface Methodology (RSM) and the experimental data are made. The maximum deviations for transient stagnation Nusselt number and transient average Nusselt number for the cases of stationary MCM disk are 10.05% and 11.82%, respectively; and 9.41% and 12.44% for the cases of rotating MCM disk, respectively. Finally, with the Sequential Quadratic Programming (SQP) technique, a series of thermal optimal designs under space limitation constraint H/d≤12 has been efficiently performed. Comparisons between the numerical optimization results and the experimental data are made with a satisfactory agreement.