variavel0=Amauri Pereira Oliveira - apdolive@usp.br Universidade de São Paulo Jacyra Soares - jacyra@usp.br Universidade de São Paulo Hugo Karam - hakaram@model.iag.usp.br Universidade de São Paulo Maxsuel Pereira - maxsuel@model.iag.usp.br Universidade de São Paulo Edson Marques Filho - emarques@model.iag.usp.br Universidade de São Paulo Abstract. This work describes the major available techniques to simulate the time and space evolution of the planetary boundary layer. For homogeneous and equilibrium conditions the structure of the planetary boundary layer can be diagnosed from the Monin-Obukhov, Free Convection, Local and Mixed Layer Similarity theories. For the other atmospheric conditions the planetary boundary layer can be numerically simulated using first and second order closure models and large eddy models. The closure models take into consideration the traditional statistical approach. Large eddy simulation models are based on the filtered equations of motion and require the statistical approach to estimate subgrid turbulence. Keywords. Numerical Modeling, Large eddy simulation, Convective boundary layer, Dispersion of inert pollutants, Planetary boundary layer.