We develop a novel calibration approach to address the problem of predictive k − ε RANS simulations of jet-incrossflow. Our approach is based on the hypothesis that predictive k − ε parameters can be obtained by estimating them from a strongly vortical flow, specifically, flow over a square cylinder. In this study, we estimate three k − ε parameters, C µ , C ε2 and C ε1 by fitting 2D RANS simulations to experimental data. We use polynomial surrogates of 2D RANS for this purpose. We conduct an ensemble of 2D RANS runs using samples of (C µ ,C ε2 ,C ε1 ) and regress Reynolds stresses to the samples using a simple polynomial. We then use this surrogate of the 2D RANS model to infer a joint distribution for the k − ε parameters by solving a Bayesian inverse problem, conditioned on the experimental data. The calibrated (C µ ,C ε2 ,C ε1 ) distribution is used to seed an ensemble of 3D jet-in-crossflow simulations. We compare the ensemble's predictions of the flowfield, at two planes, to PIV measurements and estimate the predictive skill of the calibrated 3D RANS model. We also compare it against 3D RANS predictions using the nominal (uncalibrated) values of (C µ ,C ε2 ,C ε1 ), and find that calibration delivers a significant improvement to the predictive skill of the 3D RANS model. We repeat the calibration using surrogate models based on kriging and find that the calibration, based on these more accurate models, is not much better that those obtained with simple polynomial surrogates. We discuss the reasons for this rather surprising outcome.