A framework is presented for automated tuning of lumped parameter network boundary conditions in cardiovascular simulations. An automated Bayesian approach is employed, achieving excellent agreement with clinical targets. The method is demonstrated using patient specific data for six coronary bypass patients and one normal. Uncertainties are propagated to 3D simulation results to demonstrate success of the method in a multi-scale framework.