Large-eddy simulation (LES) of turbulent convection is discussed in various versions (mixing-length theory, modal theory and spectral theory) in respect to the application to stellar convective zones. For the model construction, the non-local mixing-length formalism is suitable. However, for the determination of basic flow patterns and of mixing-length, the quasi-linear and nonlinear modal theories are useful. The eddy diffusivities are essential in these theories, and the nonlinear treatment of convection consistent with turbulent diffusivities (of effective Reynolds number of about 20 and Prandtl number of 0.4) offers a simple method of constructing stellar models without the use of the mixing-length.