Two uncertainty quantification (UQ) techniques, latin-hypercube sampling (LHS) and polynomial chaos expansion (PCE), have been used in an initial UQ study to calculate the effect of boundary condition uncertainty on Large-eddy spray simulations. Liquid and vapor penetration as well as multidimensional liquid and vapor data were used as response variables. The Morris one-at-a-time (MOAT) screening method was used to identify the most important boundary conditions. The LHS and PCE methods both predict the same level of variability in the response variables, which was much larger than the corresponding experimental uncertainty. Nested grids were used in conjunction with the PCE method to examine the effects of subsets of boundary condition variables. Numerical modeling parameters had a much larger effect on the resulting spray predictions; the uncertainty in spray penetration or multidimensional spray contours from physically derived boundary conditions was close to the uncertainty of the measurements.