La dynamique de la stratosphère au printemps et en été reste à ce jour largement inexplorée. Or dans les contextes actuels du recouvrement de la couche d’ozone et de l’augmentation des émissions de gaz à effet de serre, une amélioration de la compréhension des processus dynamiques contrôlant la stratosphère s’avère nécessaire, afin de mieux appréhender l’évolution du climat dans le futur. Des observations satellitaires récentes du printemps/été arctique ont montré l’existence de phénomènes de transport irréversibles depuis les régions tropicales vers la région arctique. Cependant, les mécanismes associés à ces évènements restent mal connus. Ce travail de thèse consiste en l’analyse dynamique et climatologique de ces phénomènes, afin d’évaluer les mécanismes responsables de leur développement et de leur fréquence d’apparition. Une attention particulière est donnée aux rares évènements, où la signature de l’intrusion persiste dans une anomalie anticyclonique jusqu’en été, soit plusieurs mois après son établissement en région polaire. Les données des instruments satellitaires MLS/Aura et MIPAS/ENVISAT, de l’instrument ballon SPIRALE, et le modèle d’advection MIMOSA ont permis d’identifier, caractériser et quantifier ces évènements. L’analyse des conditions dynamiques a été faite à partir des données météorologiques réanalysées de L’ECMWF. Enfin, le développement d’un algorithme de détection systématique de ces intrusions a permis d’en établir une climatologie entre les années 1980 et 2011. Parmi les résultats majeurs de cette étude, il apparaît que la fréquence de ces évènements, contrôlés par l’activité ondulatoire, a fortement augmenté depuis les années 2000. Nous montrons aussi que leur développement au printemps est fortement lié à l’évolution dynamique de la stratosphère durant l’hiver et au régime de circulation intertropicale. / The stratosphere dynamics remains largely unexplored in summer and spring. In the context of the ozone layer recovery and the increasing of greenhouse gases emissions, efforts must be provided to improve our knowledge of the dynamical processes driving the stratosphere. Such improvements would lead to better future climate trends estimates. Recently, spring and summer satellites observations revealed occurrences of irreversible air masses transport from the tropics to the Arctic region. However, the associated mechanisms are poorly understood. The present work consists of dynamical and climatological analyses of these events in order to identify their causes and their occurrence frequency. In particular, we focused on the sporadic events, where the intrusion signal persists several months in the polar region, trapped within an anticyclonic anomaly before disappearing in summer. We used MLS/Aura and MIPAS/ENVISAT satellites data, the SPIRALE balloon borne data and the results of the advection model MIMOSA to identify, characterize and quantify these events. Stratospheric dynamical conditions are investigated using the reanalyses data of the ECMWF. Finally, a systematic algorithm to detect low-latitude intrusions has been developed and applied on MIMOSA results to perform a climatology between 1980 and 2011 The results suggest that the frequency of these events, driven by the planetary wave activity, is increasing since the 21st century. Furthermore, their occurrence in spring appears to depend on the stratospheric dynamical evolution during winter and on the tropical region dynamical regime.

Tendo em vista o aumento considerável da poltúção do ar provocado em grande parte pela industrialização e o aumento da emissão de poluentes resultantes da queima de combustíveis fósseis por veículos automotores, o presente trabalho tem como objetivo melhorar a previsão e o entendimento da dispersão turbulenta atmosférica. Para tanto, apresenta-se, pela primeira vez, uma representação analít ica para a equação de advecção-difusão-reação tridimensional transiente, com perfil de vento e coeficientes de difusão tmbulenta dependentes da altura, que modelam a dispersão de poluentes na atmosfera. A solução da equação é obtida pela combinação do método GILTT ( Generalized Integral Laplace Transform Technique) com o método da Decomposição de Adomian modificado. Consideram-se dois casos para a aplicação do modelo: no primeiro modela-se a dispersão de um poluente secundário formado por uma reação fotoquímica e no segundo caso, utiliza-se o modelo para determinar o campo de concentração de um poluente que sofre perdas e ganhos devido a influência da radiação solar. Para poder realizar essas análises propôs-se uma parametrização para o termo de reação fotoquímica. São apresentados os resultados numéricos e estatísticos, comparandose com os dados da campanha experimental da Usina Termelétrica de Candiota e com os dados de medições realizadas pela Fundação Estadual de Proteção Ambiental Henrique Luiz Roessler (FEPAM). / In view of the considerable increase of air pollution caused largely by industrialization and the increase of emission pollutants resulting from burning of fossil fuels by motor vehicles, the present work aims to improve the prediction and understanding of atmospheric turbu- lent dispersion. Therefore, is presented, for the rst time, an analytical representation to the transient three-dimensional advection-diffusion-reaction equation, with wind pro le and turbulent diffusion coefficients dependent of height, modeling the dispersion of pollutants in the atmosphere. The solution of the equation is obtained by combining of the GILTT method (Generalized Integral Laplace Transform Technique) with the modi ed Adomian Decomposition method. It is considered two cases for the application of the model: in the rst is modeled the dispersion of a secondary pollutant formed by a photochemical reaction, and in the second case the model is used to determine the concentration eld of a pollutant that suffers losses and gains due to the in uence of solar radiation. To realise these analisis a parameterization for the photochemical reaction term is proposed. Numerical and statistical results are presented, comparing with the experimental campaign data of the thermoelectric plant of Candiota and with data from measurements performed by the \Funda c~ao Estadual de Prote c~ao Ambiental Henrique Luiz Roessler" (FEPAM).

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Uma análise utilizando a série de Taylor é apresentada para se estimar a priori os erros envolvidos na solução numérica da equação de advecção unidimensional com termo fonte, através do Método dos Volumes Finitos em uma malha do tipo uniforme e uma malha não uniforme. Também faz-se um estudo a posteriori para verificar a magnitude do erro de discretização e corroborar os resultados obtidos através da análise a priori. Por meio da técnica de solução manufaturada tem-se uma solução analítica para o problema, a qual facilita a análise dos resultados numéricos encontrados, e estuda-se ainda a influência das funções de interpolação UDS e CDS e do parâmetro u na solução numérica. / An analysis based on Taylor series is presented for estimating a priori the errors involved in the numerical solution of advection equation one-dimensional with source term, using the Finite Volume Method in a mesh uniform and a nonuniform mesh. Also is accomplished a study to determine the magnitude of discretization error and corroborate the results obtained on analyzing a priori. By using the technique of solution manufactured is produced an analytical solution for the problem, which facilitates analysis of the numeric results, and was also studied the influence functions of interpolation UDS and CDS and of parameter u in the numerical solution.

Tendo em vista o aumento considerável da poltúção do ar provocado em grande parte pela industrialização e o aumento da emissão de poluentes resultantes da queima de combustíveis fósseis por veículos automotores, o presente trabalho tem como objetivo melhorar a previsão e o entendimento da dispersão turbulenta atmosférica. Para tanto, apresenta-se, pela primeira vez, uma representação analít ica para a equação de advecção-difusão-reação tridimensional transiente, com perfil de vento e coeficientes de difusão tmbulenta dependentes da altura, que modelam a dispersão de poluentes na atmosfera. A solução da equação é obtida pela combinação do método GILTT ( Generalized Integral Laplace Transform Technique) com o método da Decomposição de Adomian modificado. Consideram-se dois casos para a aplicação do modelo: no primeiro modela-se a dispersão de um poluente secundário formado por uma reação fotoquímica e no segundo caso, utiliza-se o modelo para determinar o campo de concentração de um poluente que sofre perdas e ganhos devido a influência da radiação solar. Para poder realizar essas análises propôs-se uma parametrização para o termo de reação fotoquímica. São apresentados os resultados numéricos e estatísticos, comparandose com os dados da campanha experimental da Usina Termelétrica de Candiota e com os dados de medições realizadas pela Fundação Estadual de Proteção Ambiental Henrique Luiz Roessler (FEPAM). / In view of the considerable increase of air pollution caused largely by industrialization and the increase of emission pollutants resulting from burning of fossil fuels by motor vehicles, the present work aims to improve the prediction and understanding of atmospheric turbu- lent dispersion. Therefore, is presented, for the rst time, an analytical representation to the transient three-dimensional advection-diffusion-reaction equation, with wind pro le and turbulent diffusion coefficients dependent of height, modeling the dispersion of pollutants in the atmosphere. The solution of the equation is obtained by combining of the GILTT method (Generalized Integral Laplace Transform Technique) with the modi ed Adomian Decomposition method. It is considered two cases for the application of the model: in the rst is modeled the dispersion of a secondary pollutant formed by a photochemical reaction, and in the second case the model is used to determine the concentration eld of a pollutant that suffers losses and gains due to the in uence of solar radiation. To realise these analisis a parameterization for the photochemical reaction term is proposed. Numerical and statistical results are presented, comparing with the experimental campaign data of the thermoelectric plant of Candiota and with data from measurements performed by the \Funda c~ao Estadual de Prote c~ao Ambiental Henrique Luiz Roessler" (FEPAM).

FundaÃÃo de Amparo à Pesquisa do Estado do Cearà / O presente trabalho apresenta um estudo de resistÃncia à injeÃÃo de carga orgÃnica na zona estuarina de um importante manancial da RegiÃo Metropolitana de Fortaleza, no caso o Rio CocÃ. Para tanto, foi desenvolvido um programa em C++, com a finalidade de simular o escoamento e a dispersÃo de poluentes em um trecho especÃfico do rio, que recebe influÃncia direta da marÃ. Para a simulaÃÃo do escoamento, foi utilizada uma forma das equaÃÃes de Saint Vennant, com as condiÃÃes de contorno apropriadas. Na resoluÃÃo numÃrica foi aplicado um esquema explÃcito de diferenÃas finitas, no caso o de McCormack. Uma vez calculados os perfis temporais do escoamento, a equaÃÃo de AdvecÃÃo-DispersÃo (A-D) foi utilizada, usando um fator de decaimento de 1 ordem, para simular o transporte de massa de DBO (Demanda BioquÃmica de OxigÃnio). O programa desenvolvido aplicou um mÃtodo pseudoespectral com passo de tempo explÃcito. Para a soluÃÃo dos sistemas gerados foi utilizado o algoritmo de Gauss sem condensaÃÃo pivotal. Os resultados encontrados para o escoamento atestam o perfeito ajuste do esquema de McCormack a situaÃÃes de estuÃrio, onde o movimento de marà faz com que o fluxo esteja uma hora em um sentido e em outra hora no sentido oposto. Jà os resultados de A-D demonstraram a sensibilidade dos mÃtodos espectrais Ãs condiÃÃes de contorno impostas. Ao se resolver a equaÃÃo para o trecho em estudo do Rio CocÃ, houve acÃmulo de massa no estuÃrio, indicando que alguma condiÃÃo de contorno de escoamento nÃo era condizente com a realidade. Aplicando-se a resoluÃÃo para um estuÃrio hipotÃtico, dez vezes maior que o considerado, o acÃmulo de massa desapareceu, e resultados ficaram de acordo com a fÃsica do problema. Este fato ratificou a hipÃtese que a condiÃÃo de contorno da interface com o rio (aproximadamente no KM 20) nÃo poderia ser estacionÃria, resultado que està de acordo com o encontrado por Villela (1988). / The present work represents a study of the resistance of the estuarine zone of the Cocà river in Fortaleza, Cearà - Brazil, to the injection of organic pollutant. In order to accomplish this goal a mathematical model was developed and a C++ code was written to solve the model equations. The flow and pollutant dispersion were simulated for a particular section of the Cocà River under tidal influence. The Saint-Vennant equations were solved using boundary conditions appropriated to that situation. The McCormackâs explicit approach was applied to implement a numerical solution of the equations. The mass transport was simulated by an Advection-Dispersion equation using a first order decay. A pseudo-spectral method with explicit time step was applied and the solution of the system of equations was performed by Gauss algorithm without pivot condensation. The results of the McCormackâs approach were reliable to estimate estuary conditions where the flow was influenced by the tidal movement. The results of Advection â Dispersion (A-D) equations have shown that spectral methods were sensible to the boundary conditions of the problem. On the other hands, the results of A-D equations presented a mass retention in the estuarine zone. That fact lets one to believe that the defined boundary conditions to the problem were not appropriated to that situation. Thus, another boundary conditions were defined by using a hypothetical new estuarine zone ten times bigger than the Cocà River estuarine zone. As a result of that change, the mass retention cease to exist, and the problem became physically correct. The Villelaâs (1988) hypothesis on nonstationary behavior of the estuarine zone â Cocà River interface was confirmed.

Desde o importante trabalho de Barth e Frederickson (1990), um certo número de pesquisadores têm estudado o método de Volumes Finitos de alta-ordem k-exato, por exemplo o grupo do Prof. Ollivier-Gooch: Ollivier-Gooch e van Altena (2002), Nejat (2007), Michalak (2009), etc. Outras discretizações espaciais de alta-ordem bastante populares são o método Galerkin Descontínuo e o método de Diferença Espectral; processos iterativos que involucram estes esquemas tem sido acelerados, nos últimos anos, por métodos p-multigrid. Porém, esta aceleração não tem sido aplicada no contexto do método de Volumes Finitos de alta-ordem, pelo menos para conhecimento do autor desta tese. Por isso, o objetivo desta pesquisa é adaptar o p-multigrid desenvolvido por Liang et al. (2009b) no contexto da Diferença Espectral, para o ambiente dos Volumes Finitos estudado pelo Prof. Ollivier-Gooch. A pesquisa começa implementando o solver VF-RK, de Volumes Finitos com avanço Runge-Kutta, para resolver as equações de advecção-difusão e de Euler aplicados a problemas estacionários, por exemplo, o escoamento transônico ao redor do NACA 0012. Depois, estuda-se o método p-multigrid no contexto da Diferença Espectral; o p-multigrid acelera o processo iterativo comutando níveis polinomiais de alta e de baixa-ordem. Após esse estudo, a adaptação ao âmbito dos Volumes Finitos é realizada resultando num p-multigrid relativamente mais simples porque, em contraposição com o p-multigrid para Diferença Espectral, não precisa de operadores de restrição e prolongação para a comunicação entre diferentes níveis polinomiais. A pesquisa conclui com uma comparação com o método de Volumes Finitos de 4a ordem sem p-multigrid (solver VF-RK). Nesse sentido, implementa-se o solver pMG, baseado no p-multigrid proposto, para resolver os problemas estacionários considerados na primeira parte do trabalho; o smoother do p-multigrid é o esquema Runge-Kutta do código VF-RK, e cada problema estacionário é resolvido utilizando diferentes Vciclos procurando sempre soluções de 4a ordem. Os resultados indicam que o método p-multigrid proposto é mais eficiente que o método de Volumes Finitos de 4a ordem sem p-multigrid, isto é, os dois métodos oferecem a mesma precisão mas o primeiro pode levar menos de 50% do tempo de CPU do segundo. / Since Barth and Frederickson\'s important work (Barth e Frederickson, 1990), a number of researchers have studied high-order k-exact Finite Volume method, for example Prof. Ollivier-Gooch\'s group: Ollivier-Gooch e van Altena (2002), Nejat (2007), Michalak (2009), etc. Other quite popular high-order spatial discretizations are the Discontinuous Galerkin methods and the Spectral Difference methods; the iterative processes involving these schemes have been accelerated in recent years by p-multigrid methods. However, this acceleration has not been applied in the context of the high-order Finite Volume method, at least for the knowledge of the author of this thesis. Therefore, the objective of this research is to adapt the p-multigrid developed by Liang et al. (2009b) in the context of Spectral Difference methods, to the environment of Finite Volume studied by Prof. Ollivier-Gooch. This research begins by implementing the solver VF-RK, Finite Volume solver with Runge-Kutta advance, to compute the advection-diffusion equation and Euler equations applied to steady state problems, for example, the transonic flow around NACA 0012. Then, it is studied the p-multigrid method in the context of Spectral Difference schemes; p-multigrid accelerates the iterative process by switching polynomial levels of high- and low-order. After this study, the adaptation to the context of the Finite Volume scheme is performed resulting in a relatively simple p-multigrid because, in contrast to the p-multigrid for Spectral Difference schemes, it doesn\'t need restriction and prolongation operators for communication between different polynomial levels. The research concludes with a comparison with 4th order Finite Volume method without p-multigrid (solver VF-RK). Accordingly, the solver pMG, based on the proposed p-multigrid, is implemented to resolve the steady state problems considered in the first part of the work; the p-multigrid smoother is the Runge-Kutta scheme from VF-RK code, and each steady state problem is solved using different Vcycles, looking for 4th order solutions ever. The results indicate that the proposed p-multigrid method is more efficient than the 4th order Finite Volume method without p-multigrid: the two methods give the same accuracy but the first one can take less than 50% of second one\'s CPU time.

Este trabalho apresenta um estudo sobre a solução numérica da equação geral de advecção-difusão usando uma metodologia numérica conservativa. Para a discretização espacial, é usado o Método de Volumes Finitos devido à natureza conservativa da equação em questão. O método é configurado de modo a ter suas variáveis centradas em centro de célula e, para as variáveis, como a velocidade, centradas nas faces um método de interpolação de segunda ordem é utilizado para um ajuste numérico ao centro. Embora a implementação computacional tenha sido feita de forma paramétrica de maneira a acomodar outros esquemas numéricos, a discretização temporal dá ênfase ao Método de Crank-Nicolson. Tal método numérico, sendo ele implícito, dá origem a um sistema linear de equações que, aqui, é resolvido empregando-se o Método Multigrid-Multinível. A corretude do código implementado é verificada a partir de testes por soluções manufaturadas, de modo a checar se a ordem de convergência prevista em teoria é alcançada pelos métodos numéricos. Um jato laminar é simulado, com o acoplamento entre a equação de Navier-Stokes e a equação geral de advecção-difusão, em um domínio computacional tridimensional. O jato é uma forma de vericar se o algoritmo de geração de malhas adaptativas funciona corretamente. O módulo produzido neste trabalho é baseado no código computacional AMR3D-P desenvolvido pelos grupos de pesquisa do IME-USP e o MFLab/FEMEC-UFU (Laboratório de Dinâmica de Fluidos da Universidade Federal de Uberlândia). A linguagem FORTRAN é utilizada para o desenvolvimento da metodologia numérica e as simulações foram executadas nos computadores do LabMAP(Laboratório da Matemática Aplicada do IME-USP) e do MFLab/FEMEC-UFU. / This work presents a study about the numerical solution of variable coecients advectiondi usion equation, or simply, general advection-diusion equation using a conservative numerical methodology. The Finite Volume Method is choosen as discretisation of the spatial domain because the conservative nature of the focused equation. This method is set up to have the scalar variable in a cell centered scheme and the vector quantities, such velocity, are face centered and they need a second order interpolation to get adjusted to the cell center. The computational code is parametric, in which, any implicit temporal discretisation can be choosen, but the emphasis relies on Crank-Nicolson method, a well-known second order method. The implicit nature of aforementioned method gives a linear system of equations which is solved here by the Multilevel-Multigrid method. The correctness of the computational code is checked by manufactured solution method used to inspect if the theoretical order of convergence is attained by the numerical methods. A laminar jet is simulated, coupling the Navier-Stokes equation and the general advection-diusion equation in a 3D computational domain. The jet is a good way to check the corectness of adaptative mesh renement algorithm. The module designed here is based in a previous implemented code AMR3D-P designed by IME-USP and MFLab/FEMEC-UFU (Fluid Dynamics Laboratory, Federal University of Uberlândia). The programming language used is FORTRAN and the simulations were run in LabMAP(Applied Mathematics Laboratoy at IME-USP) and MFLab/FEMEC-UFU computers.

Ce travail a porté sur la variabilité de la salinité de surface (SSS) de l'océan dans le gyre subtropical Nord Atlantique. J'ai étudié la variabilité saisonnière de la SSS en lien avec les flux d'eau douce échangés avec l'atmosphère et la circulation océanique à méso échelle, au cours de plus de deux ans, d'août 2012 à décembre 2014. Les produits issus de la mission satellitaire soil moisture and ocean salinity (SMOS) corrigés de biais systématiques aux grandes échelles ont été testés et utilisés pour restituer la variabilité méso-échelle de SSS. Nous avons de surcroit utilisé les simulations numériques à haute résolution PSY2V4R2-R4 de Mercator. Les champs issus de SMOS et des simulations ont été comparés aux données in situ de bouées dérivantes et de thermosalinographes recueillies pendant l'expérience SPURS, avec des résultats satisfaisants, en particulier en hiver, et des écarts-type de différences typiques de l'ordre de 0.15 pss. Le flux d’eau douce échangé avec l’atmosphère est le terme dominant dans le bilan saisonnier de la SSS. Ce sont des termes associés à la dynamique océanique qui le compensent partiellement. En particulier, l’entrainement des eaux sous-jacentes contribue fortement en début d’hiver. Il agit d’ordinaire à réduire la SSS, à l’exception de la région au sud du maximum de SSS, où c’est au contraire une augmentation qu’il induit. L’advection est une seconde contribution importante à la variabilité de la SSS. Elle transfert ainsi vers le nord les eaux ‘salinisées’ plus au sud dans la région du maximum de perte d’eau douce vers l’atmosphère. La contribution d'advection est fortement dépend du type de données utilisées et leur résolution spatiale. / The focus of this work is on sea surface salinity (SSS) variability in the North Atlantic subtropical gyre. We study seasonal SSS variability and its link to the atmospheric freshwater flux at the ocean surface and to ocean dynamics at meso-scales for the period August 2012 – December 2014. The products from the soil moisture and ocean salinity (SMOS) satellite mission corrected from large scale systematic errors are tested and used to retrieve meso-scale salinity features. Furthermore, the PSY2V4R2-R4 simulation produced by Mercator with a high spatial resolution is also used. The comparison of corrected SMOS SSS data and Mercator simulation with drifter's in situ and TSG measurements from the SPURS experiment shows a reasonable agreement with RMS differences on the order of 0.15 pss.The freshwater seasonal flux is the leading term in the SSS seasonal budget. To balance its effect the ocean dynamics strongly contribute. The entrainment of deeper water is strong during the winter time. It usually acts to lower SSS, except in the south of the SSS–max region where it contributes to increase salinity. Advection is the second important component responsible for the SSS variability. It transfers further north the salty water from the evaporation maximum region. The contribution of advertion term is strongly dependent on the type of data used and their spatial resolution.

Orientadores: Simone Andrea Pozza, Varese Salvador Timóteo / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Tecnologia / Made available in DSpace on 2018-08-26T20:52:16Z (GMT). No. of bitstreams: 1 Pereira_MatheusFernando_M.pdf: 7089543 bytes, checksum: 5e3fe7245aa3630ebca8dda3ddd08eb3 (MD5) Previous issue date: 2014 / Resumo: Diversas técnicas vêm sendo apresentadas para resolução da equação da difusão, a qual é empregada para estimativas da concentração de poluentes em função do espaço e do tempo, levando-se conta fatores como fonte emissora, condições meteorológicas, características do meio e velocidade em que o poluente é carreado. Neste estudo, foi empregado um algoritmo de passo variável para a resolução da equação da difusão unidimensional e avaliação da influência do parâmetro de heterogeneidade do meio, da velocidade do fluxo e do coeficiente de dispersão na variação da concentração de poluentes em função do espaço e do tempo. As simulações foram realizadas utilizando as mesmas condições iniciais e de contorno adotadas em dois estudos abordados recentemente na literatura, e de acordo com os resultados, verificou-se que características como meios de menor heterogeneidade, baixa velocidade inicial do fluxo e baixo coeficiente de dispersão implicam em menores valores de concentração, facilitando a dispersão de poluentes. O método utilizado é caracterizado pela rápida convergência, simplicidade do código e baixo tempo computacional, podendo ser utilizado como base para resolução da equação da difusão bi e tridimensional / Abstract: Several techniques have been employed for solving advection-diffusion equation, which is used to estimate pollutants concentration as function of time and space, taking account factors such as emission source, meteorological conditions, medium characteristics and the velocity in which pollutant is adduced. In this study, we used an adaptive-step algorithm for solving one-dimensional advection-diffusion equation, and evaluating the influence of medium inhomogeneity parameter, flow velocity and dispersion coefficient in the pollutants concentration variation as function of space and time. Simulations were performed using the same initial and boundary conditions adopted by Kumar et al. (2010) and by Savovic and Djordjevich (2012), and according to the results, it was found that characteristics such as medium of less inhomogeneity, low initial flow velocity and low dispersion coefficient imply in lower concentration and facilitate pollutants dispersion. The method is characterized by rapid convergence, simplicity of the code and low computational time, and it can be used as a basis for solving the two and the three dimensional advection-diffusion equation / Mestrado / Tecnologia e Inovação / Mestre em Tecnologia

Zwei numerische Verfahren, das Upstream- und das MacCormack-Verfahren, werden im Hinblick auf ihre Eigenschaften bei der Lösung unterschiedlicher Advektionsprobleme verglichen. Das Letzere, welches eine Genauigkeit zweiter Ordnung liefert, wird in den Impulsadvektionsteil des numerischen Strömungsmodells MISKAM implementiert und anhand von Sensitivitätsstudien mit dem Upstream-Verfahren verglichen. Anschließend findet eine Modellevaluierung mit Hilfe von Datensätzen, welche am Hamburger Grenzschichtwindkanal erzeugt wurden, statt. Das Verfahren zeigt bei stabiler thermischer Modellschichtung signifikant veränderte Ergebnisse, deutlich zeigt sich die verringerte numerische Diffusion, vor allem im Bereich von Ecken und Kanten eines Hindernisses. / Two numerical schemes are compared concerning their numerical abilities when solving different advection problems. MacCormack\''s scheme which is of second order accuracy is implemented in the numerical flow model MISKAM in order to calculate the advection of momentum. It is compared to the upstream scheme with the help of sensitivity studies and with a model evaluation using wind tunnel data from the University of Hamburg. The scheme shows for stable thermical stratification significant differences near the edges of obstacles that result mainly from the reduced numerical diffusion which was the major problem when using the upstream scheme.