A theoretical study is presented of three-dimensional turbulent flow provoked in a boundary layer by an array of low-profile vortex generators (VGs) on the surface. The typical VG sits in the logarithmic region of the incident boundary layer, and the turbulence model used seems representative in this region. The governing equations yield a forward-marching three-dimensional vortex-type system, which is solved computationally and analytically for spanwise periodic VG arrays. Streamwise vortex patterns of various strengths are produced downstream, owing to three-dimensional distortion of the original logarithmic profile and to the turbulent stresses present. Predictions are given for certain basic VG shapes, e.g. triangular, with various spanwise spacings, and the predictions are found to agree favourably overall with recent experiments. In addition, the analytical formulae obtained prove useful in suggesting designs for favourable VG distributions, based on three factors : close spanwise packing, increased VG length, and suitably non-smooth spanwise shaping.