A Reynolds-averaged Navier-Stokes solver, a laminar boundary-layer code and two e N -database methods for the prediction of transition due to Tollmien-Schlichting and crossflow instabilities were coupled to be applied to 3-D high-lift aircraft configurations which are of industrial relevance. The first application of the coupled system to a wing-body configuration with a three-element wing consisting of slat, main wing, and flap is described and documented in this paper. The prediction of the laminar-turbulent transition lines was done in a fully automatic manner. It will be shown that complex aircraft configurations can be handled without a priori knowledge of the transition characteristics of the specific flow problem. The computational results are compared with experimental data. Nomenclature b = semispan c = local chord length c p = pressure coefficient k cyc = number of RANS cycles for the transition location iteration, which represents the interval between two calls of the transition prediction module n T = global number of transition points Tu = turbulence intensity x = longitudinal coordinate of the configuration in the global coordinate system of the RANS solver x T = longitudinal coordinate value of the transition point z = spanwise coordinate being perpendicular to the longitudinal axis = angle of attack = value of the intermittency F = flap deflection angle S = slat deflection angle = nondimensional spanwise coordinate Subscripts CF = crossflow crit = critical j = counter of the transition points TS = Tollmien-Schlichting 1 = freestream value