Previous research on self-similar mixing caused by Rayleigh–Taylor (RT) instability is summarized and a recent series of high resolution large eddy simulations is described. Mesh sizes of approximately 2000 ×1000 × 1000 are used to investigate the properties of high Reynolds number self-similar RT mixing at a range of density ratios from 1.5 : 1 to 20 : 1. In some cases, mixing evolves from ‘small random perturbations’. In other cases, random long wavelength perturbations (
k
−3
spectrum) are added to give self-similar mixing at an enhanced rate, more typical of that observed in experiments. The properties of the turbulent mixing zone (volume fraction distributions, turbulence kinetic energy, molecular mixing parameter, etc.) are related to the RT growth rate parameter,
α
. Comparisons are made with experimental data on the internal structure and the asymmetry of the mixing zone (spike distance/bubble distance). The main purpose of this series of simulations is to provide data for calibration of engineering models (e.g. Reynolds-averaged Navier–Stokes models). It is argued that the influence of initial conditions is likely to be significant in most applications and the implications of this for engineering modelling are discussed.