Flow topology optimization (TopOpt) based on Darcy’s source term is widely used in the field of TopOpt. It has a high degree of freedom, making it suitable for conceptual aerodynamic design. Two problems of TopOpt are addressed in this paper to apply the TopOpt method to high-Reynolds-number turbulent flow that is often encountered in aerodynamic design. First, a strategy for setting Darcy’s source term is proposed based on the relationship between the magnitude of the source term and some characteristic variables of the flow (length scale, freestream velocity, and fluid viscosity). Second, we construct two modified turbulence models, a modified Launder–Sharma k − ϵ (LSKE) model and a modified shear stress transport (SST) model, that consider the influence of Darcy’s source term on turbulence and the wall-distance field. The TopOpt of a low-drag profile in turbulent flow is studied using the modified LSKE model. It is demonstrated by comparing velocity profiles that the model can reflect the influence of solids on turbulence at Reynolds numbers as high as one million. The TopOpt of a rotor-like geometry, which is of great importance in aerodynamic design, is conducted using the modified SST model. In all the cases considered, the drag, the total pressure loss, and the energy dissipation are significantly reduced by TopOpt, indicating the proposed model’s ability to handle the TopOpt of turbulent flow.