The complex turbulent flow regimes encountered in many thermal-fluid engineering applications have proven resistant to the effective application of systematic design because of the computational expense of model evaluation and the inherent variability of turbulent systems. In this thesis the integration of the Proper Orthogonal Decomposition (POD) for reduced order modeling of turbulent convection with the application of robust design principles is proposed as a practical design approach. The POD has been used successfully to create low dimensional steady state flow models within a prescribed range of parameters. The underlying foundation of robust design is to determine superior solutions to design problems by minimizing the effects of variation on system performance, without eliminating their causes. The integration of these constructs utilizing the compromise Decision Support Problem (DSP) results in an efficient, effective robust design approach for complex turbulent convective systems. The efficacy of the approach is illustrated through application to the configuration of data center server cabinets. Data centers are computing infrastructures that house large quantities of data processing equipment. The data processing equipment is stored in 2 m high enclosures known as cabinets. The demand for increased computational performance has led to very high power density cabinet design, with a single cabinet dissipating up to 20 kW. The computer servers are cooled by turbulent convection and have unsteady heat generation and cooling air flows, yielding substantial inherent variability, yet require some of the most stringent operational requirements of any engineering system. Through variation of the power load distribution and flow parameters, such as the rate of cooling air supplied, thermally efficient configurations that are insensitive to variations in operating conditions are determined. This robust design approach is applied to three common data center server cabinet designs, in increasing levels of modeling detail and complexity. Results of the application of this approach to the example problems studied show that the resulting thermally efficient configurations are capable of dissipating up to a 50% greater heat load and 15% decrease in the temperature variability using the same cooling infrastructure. These results are validated rigorously, including comparison of detailed CFD simulations with experimentally gathered temperature data of a mock server cabinet. Finally, with the approach validated, augmentations to the approach are considered for multi-scale design, extending approaches domain of applicability.