We present a multiâfidelity blackâbox optimization approach for integrated design and control (IDC) of constrained nonlinear systems in the presence of uncertainty. The IDC framework is becoming increasingly important for the systematic design of nextâgeneration (flexible) manufacturing and energy systems. However, identifying optimal solutions to realistic IDC problems is intractable when (i) the dynamics occur on much shorter timescales than the system lifetime, (ii) the uncertainties are described by continuous random variables with high variance, and (iii) operational decisions involve a mixture of discrete and continuous variables. Instead of aggressively simplifying the problem to improve tractability, we develop a simulationâbased optimization procedure using highâquality decision rules that map information that can be measured online to optimal control actions. In particular, we rely on the Bayesian optimization (BO) framework that has been shown to perform very well on noisy and expensiveâtoâevaluate objective functions. We also discuss how BO can be extended to take advantage of computationally cheaper lowâfidelity approximations to the highâfidelity IDC cost function. Three major lowâfidelity approximation strategies are described in this work, which are related to the simplification of the system simulator, decision rule solution method, and time grid. Lastly, we demonstrate the advantages of multiâfidelity BO on the design of a solarâpowered building heating/cooling system (with battery and grid support) under uncertain weather and demand conditions with hourly variation over a yearâlong planning horizon.