We introduce the Lévy fractional Brownian field family to model the turbulent wavefront phase. This generalized model allows us to overcome the limitations found in a previous approach [Perez et al., J. Opt. Soc. Am. A 21, 1962 (2004)]. More precisely, our new model provides stationary phase increments over the full inertial range. Thus it successfully extends classical results to non-Kolmogorov turbulence without any approximation.