<p>This paper presents the shape optimization of an idealized arterial bypass graft under uncertainties. The underlying blood flow problem is numerically solved using computational fluid dynamics (CFD) simulations, able to account for mechanical hemolysis and non-Newtonian viscosity properties of blood. The employed hemolysis and non-Newtonian models utilize a set of parameters that are considered uncertain in this work. To this end, the optimization problem is rendered robust and targets the minimization of the first two statistical moments of a hemolysis index. The propagation of the uncertainties to the hemolysis index is done through the first-order second-moment method (FOSM). The necessary derivatives are computed by the adjoint method and finite differences. Several steady-state robust shape optimization simulations of the idealized bypass graft are presented. Selected optimized shapes are further assessed by means of additional fluid–structure interaction (FSI) simulations under unsteady conditions.</p>