... response surface models in CFD-based shape optimization. For each numerical optimization objective considered herein, we observe that polynomial response surfaces of orders higher than 2 (the default choice in the literature) not only improve the RSM fit but also produce better optimization results...