A robust two-phase flow Large Eddy Simulation (LES) algorithm has been developed and applied to predict the primary breakup of an axisymmetric water jet injected into a surrounding coaxial air flow. The high liquid/gas density and viscosity ratios are known to represent a significant challenge in numerical modelling of the primary breakup process. In the current LES methodology, an extrapolated liquid velocity field was used to minimise discretisation errors, whilst maintaining sharp treatment of fluid properties across the interface. The proposed numerical approach showed excellent robustness and high accuracy in predicting coaxial liquid jet primary breakup. Since strong turbulence structures will develop inside the injector at high Reynolds numbers and affect the subsequent primary breakup, the Rescaling and Recycling Method ((RM)-M-2) was implemented to facilitate generation of appropriate unsteady LES inlet conditions for both phases. The influence of inflowing liquid and gas turbulent structures on the initial interface instability was investigated. It is shown that liquid turbulent eddies play the dominant role in the initial development of liquid jet surface disturbance and distortion for the flow conditions considered. When turbulent inflows were specified by the (RM)-M-2 technique, the predicted core breakup lengths at different air/water velocities agreed closely with experimental data. (C) 2013 The Authors. Published by Elsevier Ltd. All rights reserved.