Simulation of pollutant dispersion in the planetary boundary layer, solving the langevin equation by decomposition method

Neste trabalho é apresentado uma solução analítica para a equação de Langevin tridimensional e estocastica aplicada na dispersão de poluentes na atmosfera considerando as seguintes funções densidade probabilidade (PDF): Gaussiana, Bi-Gaussiana e Gram- Charlier. A solução é obtida usando o Método de Decomposição Adomian (ADM), que é um método para resolver equações diferenciais não-lineares sem linearização. O método de decomposição consiste na expansão da solução em série de funções e o termo não-linear em série de polinômios definidos por Adomian. Substituindo estas expansões na equação µa ser resolvida, é construído um sistema linear recursivo, que é então resolvido de maneira analítica. Também é apresentado um estudo de estabilidade baseado na teoria de Lyapunov, bem como é introduzido um novo índice estatístico para a validação do modelo. Os resultados obtidos por esta metodologia são comparados com os dados do experimento de Copenhagen, bem como com os resultados obtidos a partir de outros modelos Lagrangeanos: Ito, ILS e a solução analítica. Na comparação com os dados experimentais obtidos pelo modelo proposto e o método ILS foram o que apresentaram os melhores resultados.This work presents an analytical solution of the three-dimensional stochastic Langevin equation and applied to the dispersion of pollutants in the atmosphere considering the fol- lowing probability density functions (PDF): Gaussian, Bi-Gaussian and Gram-Charlier. A solution is obtained using the Adomian Decomposition Method (ADM), which is a method for solving non-linear di®erential equations without the use of linearization. The decompo- sition method consists in expanding the solution in a series of fuction and the non-linear term in a series de¯ned by Adomian polynomials. Upon substitution of these expansions in the equation to be solved, one may built a linear recursive system which is then solved analytically. Further, a study of convergence stability based on Lyapunov theory is presented and a new statistical index for model validation is introduced. The results obtained by this method are compared with the experimental data from Copenhagen, as well as the results obtained from other Lagrangian models: Ito, ILS and an analytical solution. In comparison with the experimental data the proposed model and the ILS method showed the best results