Aplicação das fórmulas de Vincenty nos cálculos das correções dos efeitos do relevo na gravidade e na altura geoidal.

O presente trabalho apresenta a aplicação das fórmulas de Vincenty nos cálculos das correções do terreno e do efeito indireto, que desempenham papel relevante na construção de cartas geoidais. Implementa-se um programa de processamento que realiza a integração numérica sobre o modelo digital do terreno, discretizado em células triangulares de Delaunay. O sistema foi desenvolvido com a linguagem de programação FORTRAN, para a execução de intensos algoritmos numéricos usando compiladores livres e robustos. Para o cálculo do efeito indireto, considera-se a redução gravimétrica efetuada com base no segundo método de condensação de Helmert, face ao pequeno valor de efeito indireto no cálculo do geóide, em função da mudança que este produz no potencial da gravidade devido ao deslocamento da massa topográfica. Utiliza-se, o sistema geodésico SIRGAS 2000 como sistema de referência para o cômputo das correções. Simplificando o exame dos resultados alcançados, distingue-se o processamento e desenvolvimento do trabalho em etapas como a escolha de ferramentas geodésicas para máxima precisão dos resultados, elaboração de subrotinas e comparação de resultados com cálculos anteriores. Os resultados encontrados foram de geração sadia e satisfatória e podem ser perfeitamente empregados no cálculo do geóide em qualquer área do globo.

Aplicação da geocomputação a estudos de susceptibilidade de deslizamentos de terra

Esta tese propôs uma metodologia para detecção de áreas susceptíveis a deslizamentos de terra a partir de imagens aéreas, culminando no desenvolvimento de uma ferramenta computacional, denominada SASD/T, para testar a metodologia. Para justificar esta pesquisa, um levantamento sobre os desastres naturais da história brasileira relacionada a deslizamentos de terra e as metodologias utilizadas para a detecção e análise de áreas susceptíveis a deslizamentos de terra foi realizado. Estudos preliminares de visualização 3D e conceitos relacionados ao mapeamento 3D foram realizados. Estereoscopia foi implementada para visualizar tridimensionalmente a região selecionada. As altitudes foram encontradas através de paralaxe, a partir dos pontos homólogos encontrados pelo algoritmo SIFT. Os experimentos foram realizados com imagens da cidade de Nova Friburgo. O experimento inicial mostrou que o resultado obtido utilizando SIFT em conjunto com o filtro proposto, foi bastante significativo ao ser comparado com os resultados de Fernandes (2008) e Carmo (2010), devido ao número de pontos homólogos encontrados e da superfície gerada. Para detectar os locais susceptíveis a deslizamentos, informações como altitude, declividade, orientação e curvatura foram extraídas dos pares estéreos e, em conjunto com as variáveis inseridas pelo usuário, forneceram uma análise de quão uma determinada área é susceptível a deslizamentos. A metodologia proposta pode ser estendida para a avaliação e previsão de riscos de deslizamento de terra de qualquer outra região, uma vez que permite a interação com o usuário, de modo que este especifique as características, os itens e as ponderações necessárias à análise em questão.

A contribuição da geomática na geografia da saúde sob uma abordagem da Teoria Geral de Sistemas.

O conjunto de técnicas de Geomática tem sido um poderoso aliado no apoio à tomada de decisões e manipulação de dados espaciais. Nas últimas décadas tornou-se quase indispensável na gestão de recursos humanos e recursos naturais. Na área da Saúde, a Geomática tem viabilizado estudos sobre: análise da distribuição de pacientes; variações na ocorrência de epidemias; monitoramento de vetores; avaliação em tempo real de situações de emergência ou catastróficas, entre outros. A associação da Medicina com a Geografia é antiga. A associação entre a saúde e a espacialização das doenças deve considerar a estruturação física do espaço geográfico a partir da caracterização de cada parte componente do sistema como preconiza a perspectiva sistêmica na Teoria Geral de Sistema, pois a população é parte integrante do espaço geográfico no qual ocorrem determinados agravos à saúde. Nestes pontos recai a relevância da presente pesquisa. O estudo foi desenvolvido na área da Bacia Hidrográfica de Jacarepaguá no município do Rio de Janeiro. Com aproximadamente 295 km2 é composta por 19 bairros e tem uma população aproximada de 682000 habitantes. O estudo pretendeu estruturar uma metodologia para inserção das técnicas de Geomática na prática de análises de dados de saúde, apoiada na Teoria Geral de Sistemas, visando contribuir no processo de vigilância ambiental em saúde permitindo realizar um prognóstico de situações de risco de determinada população. A base cartográfica que apoiou o estudo foi construída através de cartas topográficas da DSG na escala 1:50000. Das cartas foram retiradas informações de planimetria relevantes ao estudo e as informações de altimetria que foram utilizadas na elaboração do modelo digital de terreno (MDT). Com o MDT foi gerado o mapa de declividade por grau que permitiu criar divisões em 5 classes para posterior geração de mapas temáticos com os dados de saúde. Os dados de saúde foram retirados do Sistema de Informações de Notificação de Agravos (SINAN) disponibilizado pelo Ministério da Saúde. A doença analisada foi a leptospirose. Os casos foram localizados pontualmente a partir dos endereços de residência. A hipótese aventada neste estudo era a existência de um alto grau de correlação entre as partes componentes representadas pela estruturação física do espaço e a ocorrência de leptospirose. Com os resultados comprovou-se a hipótese nula, pois somente ocorrem doenças em áreas associadas ao Fundo Chato da Baixada de Jacarepaguá e nos Declives Suaves dos Morros Isolados, também localizados no Subsistema Baixada. Isto demonstra que em termos desta doença existe uma ausência de população nesta componente espacial em áreas de maiores declividades.

Caracterização geomorfométrica multiescala de modelos digitais de terreno.

O estudo das formas de relevo através da geomorfometria ciência que quantifica as superfícies topográficas demanda que múltiplas escalas sejam consideradas simultaneamente. Esse fato soma-se à continuidade espacial inerente à topografia, a qual raramente apresenta limites naturais abruptos, para conceder um caráter vago às formas de relevo. Essa dupla vagueza, no entanto, é raramente abordada de forma explícita em trabalhos científicos que, usualmente, valem-se de limiares e definições subjetivas e arbitrárias para descreverem o relevo. Dentro desse contexto, o trabalho apresenta um método inovador de caracterização geomorfométrica multiescala de modelos digitais de terreno (MDTs). A abordagem proposta baseia-se no método de extração de feições morfométricas de Wood (1996), mas o modifica em um número de maneiras: (i) expandindo o número de classes identificadas; (ii) transformando-o em um sistema fuzzy, cujos conjuntos fuzzy são parametrizados automaticamente e; (iii) limitando localmente a escala máxima de análise de maneira não-supervisionada. Como resultado, obtém-se um mapa de feições fundamentais e um mapa de escalas fundamentais que, juntos, sintetizam a estrutura multiescala das superfícies. Além disso, são produzidos mapas de pertinências fuzzy e de índice de confusão para cada escala analisada, assim como versões multiescala dos mesmos. Para avaliar a transferibilidade e o caráter não-supervisionado do método, foram analisados cinco MDTs oriundos de bases de dados distintas, com diferentes resoluções e extensões espaciais, compreendendo regiões continentais, do fundo dos oceanos e do planeta Marte. Os resultados obtidos foram avaliados em relação à utilização de escalas e parametrizações fixas, atestando a capacidade do método de empreender caracterizações geomorfométricas mais completas do que abordagens convencionais. Duas aplicações foram ainda propostas: a parametrização geomorfométrica multiescala e o desenvolvimento de uma assinatura morfométrica multiescala, demonstrando claros caminhos para a continuidade da pesquisa. Como conclusão geral, não obstante a algumas limitações apontadas, considerou-se que o trabalho apresentado alcançou seu objetivo de prover a caracterização geomorfométrica multiescala de modelos digitais de terreno.

Simulação de elevação do nível do mar na praia do Leme - RJ.

Os crescentes alertas sobre as mudanças climáticas e suas consequências vêm preocupando a sociedade de modo geral. Geram dúvidas sobre as reais modificações que podem ser ocasionadas, principalmente com relação aos efeitos da elevação do nível médio dos oceanos e seus efeitos nas regiões costeiras. Por vezes, as informações veiculadas contribuem para uma percepção limitada da extensão espaço-temporal dos fenômenos. Em paralelo, a democratização de mapas 2D e 3D está cada vez mais ampla, não sendo mais uma ferramenta direcionada apenas aos profissionais ligados à área da Cartografia e demais Geociências. Este é um dos enormes benefícios proporcionados pela Cartografia Digital. Porém, não teria tanto alcance sem a associação das informações a um poderoso veículo de divulgação como a Web. A Internet é um dos meios de comunicação mais democráticos do mundo e uma ferramenta de grande alcance na sociedade. O uso dessas ferramentas pode ser ainda mais explorado para a construção do conhecimento, permitindo aos cidadãos entenderem mais facilmente um determinado fenômeno por meio da sua visualização/simulação em tempo real. Isto melhora a formação de opinião e o posicionamento da sociedade sobre o caso em questão. Com base nessa argumentação, esta dissertação tem o objetivo de simular cenários de elevações do nível médio do mar e disponibilizá-los em formato de mapas interativos na Web utilizando um Modelo Digital de Elevação suficientemente acurado. Para isto, a área de estudo escolhida foi a Praia do Leme, situada no município do Rio de Janeiro, RJ. Os processos para atingir os objetivos envolveram etapas de coleta de dados de campo, confecção de um banco de dados integrado a dados pré-existentes, interpolações, simulação de cenários da elevação do nível do mar e implementação de uma página Web. Os cenários investigados foram obtidos por simulações de elevações do nível do mar a cada 0,5m até 4,0m acima da cota atual. Os resultados obtidos mostram que a Praia do Leme sofrerá modificações geomorfológicas, com perda de até 18,0m da praia, para uma elevação do nível médio do mar de 2,0m. Entretanto, a orla litorânea sofrerá impactos significativos a partir de uma elevação do nível do mar de 3,0m, em que o avanço planimétrico do mar pode ficar entre 76 e 90m, o que atingiria o calçadão e a região do emissário de águas pluviais do Leme, respectivamente. Isto resultaria em um volume adicional de água na praia de, aproximadamente, 106.000m3 e uma perda de até 41.400m3 de areia. Somente a partir de 3,5m de elevação, o mar atingiria a Avenida Atlântica, sendo que, desta altitude em diante, todo o bairro do Leme seria gradualmente inundado pelo avanço do mar.

Cenários de mudanças climáticas usando modelagem dinâmica na Bacia do Alto Taquari.

A região Centro-Oeste do Brasil tornou-se nos últimos 40 anos grande produtora de grãos e carne bovina. As condições edafoclimáticas, o sistema de manejo do solo e o descumprimento de leis ambientais trouxeram conseqüências drásticas à região como o agravamento do processo hídrico erosivo, principalmente na Bacia do Alto Taquari (BAT). Cerca de 90% da BAT localiza-se na porção norte do estado de Mato Grosso do Sul (MS), porém os efeitos do transporte de sedimentos e volume de água são refletidos a jusante dos rios, na Bacia do Pantanal. Utilizando-se pressupostos do Painel Intergovernamental de Mudanças Climáticas (IPCC) foram estabelecidos cenários de mudanças climáticas na Bacia do Alto Taquari, visando identificar áreas com maior vulnerabilidade ao processo erosivo em função de pressões de uso da terra. Usando a modelagem dinâmica no TerraME (Environment Modeling) foram gerados cenários topopluviais até 2100, considerando-se para a temperatura do ar média anual um aumento de 1C, em cenário otimista e, em pessimista, elevações térmicas de 3C. Para a precipitação pluvial média anual um cenário foi com aumento de 15% e outro com reduções de 15%. Os dados foram espacializados no ArcGis 9.2 e exportados para o TerraView 3.2, criando-se espaços celulares e integrando-se com as informações do modelo digital do terreno do Shuttle Radar Topography Mission (SRTM) para geração dos mapas topoclimáticos e simulações de cenários no TerraMe. Os resultados apontam que 85% da área da BAT nas condições atuais as temperaturas médias variam entre 23,6 a 25,7C. As simulações térmicas no cenário otimista indicam que em 40 anos as temperaturas tendem a superar o maior limite térmico médio nas áreas ao longo do rio Taquari, no sentido Oeste-Leste. Esses valores evidenciam elevações nas taxas evapotranspiratórias de matas ciliares, indicando reduções na vazão do Taquari. Em cenário pessimista essas temperaturas antecipam sua ocorrência, em um prazo de 20 anos. Os cenários com acréscimo de 15% na precipitação pluvial mostram aumentos no volume de água precipitada na parte norte da Bacia, região mais vulnerável aos problemas de erosão hídrica. Cenários do regime térmico-hídrico apontam áreas mais sensíveis às mudanças climáticas na parte oeste da BAT e impactos ambientais também na Bacia do Pantanal. Conclui-se que o TerraME é indicado para gerar cenários de mudanças climáticas em bacias hidrográficas.

Evolução geológica do Complexo Quirino, Terreno Paraíba do Sul, Setor Central da Faixa Ribeira, com base em dados isotópicos de Sm-Nd e Sr

A área estudada está inserida na Faixa Ribeira, Segmento Central da Província Mantiqueira (Almeida et al., 1973, 1977, 1981), que representa um cinturão de dobramentos e empurrões gerado no Neo-proterozóico/Cambriano, durante a Orogênese Brasiliana, na borda sul/sudeste do Cráton do São Francisco (Almeida, 1971, 1977; Cordani et al., 1967, 1973; Cordani & Brito Neves, 1982; Teixeira & Figueiredo, 1991). Neste contexto, o Complexo Quirino é o embasamento retrabalhado do Terreno Paraíba do Sul (Heilbron et al., 2004). O Complexo Quirino é formado por extensos corpos de ortognaisses foliados a homogêneos, leuco a mesocráticos, de granulometria média à grossa, composicionalmente variando entre granitóides tonalíticos/granodioríticos a graníticos, e apresentando enclaves de rochas ultramáficas, máficas e cálcio-silicáticas (ricas em tremolita). Os ortognaisses tonalíticos/granodioríticos apresentam porfiroblastos de plagioclásio e a hornblenda como máfico principal, contrastando com os de composição granítica que apresentam porfiroblastos de K-feldspato e biotita predominante. Como acessórios aparecem zircão, titanita, apatita e epidoto. Também estão associados a estes ortognaisses, granitóides neoproterozóicos que formam corpos individualizados ou lentes anatéticas no conjunto paleoproterozóico. Estes são compostos predominantemente por biotita gnaisse e hornblenda-biotita gnaisse. A análise litogeoquímicas dos ortognaisses do Complexo Quirino demonstrou a existência de duas séries magmáticas distintas. A primeira pertencente à série cálcio-alcalina de alto-K apresenta uma composição mais expandida granítica-adamelítica/granodioritica/tonalítica e é correlacionável aos bt-ortognaisses e alguns hb-bt-ortognaisses. Os ortognaisses da série médio-K apresentam composição predominantemente tonalítica, sendo correlacionáveis à maioria dos hornblenda-biotita gnaisses. Enclaves lenticulares de metapiroxeníticos e anfibolíticos ocorrem em muitos afloramentos. Também ocorrem granitóides neoproterozóicos de composição graníticas a quartzo-monzoníticas O estudo isotópico de Sm-Nd e Sr demonstrou que os ortognaisses da série cálcio-alcalina de alto-K e aqueles da série cálcio-alcalina de médio-K possuem idades modelo TDM variando entre paleoproterozóicas a arqueanas, consistentes com dados U-Pb em zircão publicados na literatura. A série cálcio-alcalina de alto-K é mais antiga (2308 9,2 Ma a 2185 8 Ma) do que a série calcio-alcalina de médio-K (2169 3 a 2136 14 Ma) e a existência de zircões herdados com idades mínimas de 2846 Ma e 2981 Ma para série de médio-K e 3388 16 para série de alto-K. Os granitóides brasilianos possuem idades de cristalização neoproterozóica correlacionada a Orogênese Brasiliana (602 a 627 Ma) (Viana, 2008; Valladares et al., 2002)./Com base nos dados de Sr e Sm-Nd foi possível caracterizar 4 grupos distintos. Os grupos 1 e 2 são formados por rochas de idade paleoproterozóica (2,1 a 2,3 Ga) com idades modelo TDM variando de 2,9 e 3,4 Ga, εNd entre -8,1 e -5,8 e 87Sr/86Sr(t) = 0,694707 (Grupo 1) e TDM variando de 2,5 a 2,7 Ga, εNd entre -5,8 e -3,1 e 87Sr/86Sr(t) = 0,680824 (Grupo 2), formados no paleoproterozóico com contribuição de uma crosta arqueana. O grupo 3 é formado por rochas juvenis de idade paleoproterozóica, com idades de cristalização variando entre 2,0 e 2,2 Ga e com idades modelo TDM variando de 2,1 a 2,2 Ga e εNd entre + 1,5 e + 1,2. O grupo 4 é formado durante o neoproterozóico (645 Ma) por rochas possivelmente de idade paleoproterozóico com idades modelo TDM igual a 1,7 Ga e εNd igual a -8,3.

Informações georreferenciadas no suporte a decisão em arqueologia brasileira: proposta para uma carta arqueológica digital.

O desenvolvimento do projeto consistiu na modelagem da informação para gestão em arqueologia numa base de dados convergente de informação, caracterizada por um ambiente cartográfico digital, integrando mapas temáticos e base de dados alfanumérica num modelo protótipo e interativo para ajudar a integrar as informações na gestão em arqueologia, permitindo o registro, armazenamento e visualização da informação arqueológica correlacionada a mapas temáticos.

Os efeitos da cidade digital: o direito à privacidade e à liberdade de informação

O presente trabalho analisa a formação da cidade digital nas relações sociais, ressaltando os efeitos da garantia do direito à privacidade no ambiente dos navegantes de sites e redes sociais, em função das repercussões jurídicas do vazamento de informações da vida pessoal dos usuários da rede, e do tratamento dos dados coletados pelos prestadores de serviço. Através do ciberespaço formam-se comunidades virtuais que ultrapassam a necessidade de localidade e sociabilidade, criando um isolamento social e abandonando as interações face a face em ambientes reais, originando uma sociabilidade baseada no individualismo. Avaliamos os novos padrões de interação que se originam nesta nova formatação de coletividade informacional e suas repercussões no âmbito do direito. Em uma perspectiva mais detalhada, esse estudo indica quais as hipóteses de responsabilidade civil dos provedores na Internet em decorrência de atos ilícitos cometidos por terceiros e as alternativas de um sistema de tutela da privacidade à proteção de dados, face à lesão no ambiente informacional. O levantamento das possíveis situações de responsabilização civil foi efetuado através da análise da jurisprudência e da doutrina dominante, ressaltando os aspectos fáticos que caracterizam sua formatação. Esse modelo se impõe, através de uma relação hierárquica a uma multiplicidade de indivíduos, criando um encarceramento perfeito através do exercício do biopoder. Tais papéis são reforçados por uma cultura consumista e a sociedade do espetáculo, que transforma o indivíduo em mercadoria levantando perfis de usuários conectados em rede, propiciando uma categorização dos consumidores. Nesse contexto, apresentamos os riscos de uma sociedade de vigilância que se apresenta factível como um produto das relações de mercado, que possibilita dispor livremente de um conjunto crescente de informações. Esta constante vigilância invade todos os espaços, custodiando nosso comportamento independente do tempo, com uma implacável memória no âmbito das comunicações eletrônicas, tornando nosso passado eternamente visível e fazendo surgir situações constrangedoras a nos assombrar.

Brasil rumo aos padrões internacionais de contabilidade para o setor público: uma análise sob a ótica do Financial Management Reform Process Model de Lüder.

O termo Sistema Contábil define toda a modelagem pela qual são tratados os dados internos e externos, de forma estruturada, para atender às demandas atribuídas à profissão contábil (tanto nos sistemas privados como nos governamentais), englobando fatores intrínsecos e extrínsecos da contabilidade e seus inter-relacionamentos. As variações desses fatores levam ao desenvolvimento de sistemas nacionalmente específicos. Nas classificações internacionais desses sistemas, são identificadas influências como a dos Estados Unidos, a do Reino Unido e a da Europa, ou ainda uma orientação micro ou macro, também definidas como modelo anglo-saxão e modelo continental. Observa-se, ainda, uma relação entre os níveis de desenvolvimento contábil públicos e privados. As principais causas determinantes da diversidade contábil se relacionam com as variáveis do modelo proposto por Lüder para se avaliar reformas governamentais financeiras (o FMR) e esse fato, por conseguinte, determinou a escolha do referido modelo para realizar as análises deste estudo. A pesquisa foi realizada com procedimentos metodológicos de natureza teórica com caráter exploratório e descritivo, através da qual se pode observar que a IFAC tem sido considerada responsável pelos principais esforços no processo de harmonização contábil do setor público e que boa parte dos países vem aceitando a IFAC como órgão emissor de padrões contábeis internacionais, ainda que os países que adotam integralmente esses padrões sem adaptações, sejam exceções. Os padrões de contabilidade governamental emitidos pela IFAC (IPSAS), na verdade, têm servido como um padrão de qualidade, e não como uma norma a ser adotada integralmente. A análise das principais características nacionais determinantes da diversidade contábil frente à adoção ou adaptação dos padrões internacionais da IFAC para o setor público, sob a ótica do FMR, indica que, no Brasil, o sistema contábil do setor público é mais propício a um processo de convergência aos padrões internacionais de Contabilidade da IFAC.

Uma técnica híbrida para geração de rotas em espaço geográfico com obstáculos.

Este trabalho está inserido no campo da Geomática e se concentra, mais especificamente, no estudo de métodos para exploração e seleção de rotas em espaços geográficos sem delimitação prévia de vias trafegáveis. As atividades que poderiam se beneficiar de estudos desse tipo estão inseridas em áreas da engenharia, logística e robótica. Buscou-se, com as pesquisas realizadas nesse trabalho, elaborar um modelo computacional capaz de consultar as informações de um terreno, explorar uma grande quantidade de rotas viáveis e selecionar aquelas rotas que oferecessem as melhores condições de trajetória entre dois pontos de um mapa. Foi construído um sistema a partir do modelo computacional proposto para validar sua eficiência e aplicabilidade em diferentes casos de estudo. Para que esse sistema fosse construído, foram combinados conceitos de sistemas baseados em agentes, lógica nebulosa e planejamento de rotas em robótica. As informações de um terreno foram organizadas, consumidas e apresentadas pelo sistema criado, utilizando mapas digitais. Todas as funcionalidades do sistema foram construídas por meio de software livre. Como resultado, esse trabalho de pesquisa disponibiliza um sistema eficiente para o estudo, o planejamento ou a simulação de rotas sobre mapas digitais, a partir de um módulo de inferência nebuloso aplicado à classificação de rotas e um módulo de exploração de rotas baseado em agentes autônomos. A perspectiva para futuras aplicações utilizando o modelo computacional apresentado nesse trabalho é bastante abrangente. Acredita-se que, a partir dos resultados alcançados, esse sistema possa ajudar a reduzir custos e automatizar equipamentos em diversas atividades humanas.

Uma Metodologia de Estudo de Simulação Tridimensional de Escoamento Turbulento Estratificado no Reservatório de Plantas Hidrelétricas.

Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada.