Utilização da transformada de características invariante a escala (SIFT) na automatização da obtenção de pontos do Sistema de Imagens Tridimensional Híbrido (SITH)

Esta dissertação apresenta um aperfeiçoamento para o Sistema de Imagens Tridimensional Híbrido (SITH) que é utilizado para obtenção de uma superfície tridimensional do relevo de uma determinada região a partir de dois aerofotogramas consecutivos da mesma. A fotogrametria é a ciência e tecnologia utilizada para obter informações confiáveis a partir de imagens adquiridas por sensores. O aperfeiçoamento do SITH consistirá na automatização da obtenção dos pontos através da técnica de Transformada de Características Invariantes a Escala (SIFT - Scale Invariant Feature Transform) dos pares de imagens estereoscópicas obtidos por câmeras aéreas métricas, e na utilização de técnicas de interpolação por splines cúbicos para suavização das superfícies tridimensionais obtidas pelo mesmo, proporcionando uma visualização mais clara dos detalhes da área estudada e auxiliando em prevenções contra deslizamentos em locais de risco a partir de um planejamento urbano adequado. Os resultados computacionais mostram que a incorporação destes métodos ao programa SITH apresentaram bons resultados.

Avaliação tridimensional da via aérea superior de pacientes com apneia obstrutiva do sono que utilizaram aparelho de avanço mandibular

Os objetivos deste estudo foram avaliar o efeito do aparelho de avanço mandibular (Twin block - TB) no volume das vias aéreas superiores, por meio de tomografia computadorizada cone beam (CBCT); analisar, por meio da polissonografia, as mudanças no índice de apneia e hipopneia (IAH) e índice de apneia por hora de sono (IA), saturação de oxi-hemoglobina e eficiência do sono; e correlacionar o volume na CBCT e as polissonografias. Dezesseis pacientes portadores de apneia obstrutiva do sono, idade média de 47,06 anos, participaram deste estudo prospectivo, com acompanhamento médio de 7 meses. Foram feitas polissonografias iniciais (T1) e de acompanhamento (T2) com o TB em posição, e CBCT sem e com TB em posição. A segmentação e obtenção dos volumes das vias aéreas superiores foram realizadas e utilizados os testes t de Student pareado, de Wilcoxon e o índice de correlação de Spearman, com 5% de significância. Os resultados das polissonografias mostraram diferenças estatisticamente significativas entre T1 e T2 apenas para IAH (p<0,05). Houve aumento do volume da via aérea superior com TB quando comparado com o volume sem TB (p<0,05). Foi possível estabelecer-se correlação positiva entre volume da via aérea superior sem TB e IAH e IA em T1 (p<0,05), mas não houve correlação entre o volume da via aérea com TB e índices polissonográficos em T2. Pode-se concluir que, houve aumento de volume da via aérea superior com o TB e houve redução do IAH em T2 porém, sem correlação entre estes dados.

Resultados da expansão de maxila cirurgicamente assistida sem disjunção ptérigomaxilar: uma avaliação tridimensional

A expansão rápida da maxila cirurgicamente assistida tornou-se amplamente utilizada e muito aceitável no tratamento da deficiência maxilar de pacientes adolescentes e adultos. Diversas técnicas cirúrgicas foram propostas ao longo dos anos com o objetivo de solucionar este problema de forma eficiente, com estabilidade dos resultados e baixa morbidade. Controvérsias em relação ao procedimento cirúrgico persistem, principalmente relacionadas a quais osteotomias devem ser realizadas para se obter bons resultados. O objetivo deste trabalho foi avaliar os resultados da expansão ortocirúrgica da maxila realizando osteotomias nas paredes laterais da maxila e na sutura palatina mediana. Foram selecionados dezessete pacientes adultos portadores de deficiência transversa maxilar, com média de idade de 24 anos e 8 meses; todos foram submetidos a exames de tomografia computadorizada convencional e moldagens maxilares previamente ao procedimento cirúrgico e após três meses, no mínimo, do término de ativação e estabilização do aparelho expansor. As medidas do pós-cirúrgico foram confrontadas com as do pré-cirúrgico e os resultados foram comparados e analisados estatisticamente. Foi obtida a expansão desejada clinicamente em todos os pacientes. No entanto, a quantidade de expansão na região de molares foi estatisticamente maior nas áreas referentes aos dentes, enquanto que os resultados obtidos referentes aos caninos se mostraram similares nas três regiões maxilares avaliadas. Quando comparadas às regiões de caninos e molares entre si, a expansão intercaninos foi maior na altura dos forames palatinos e o inverso ocorreu nas regiões de processo alveolar e dentária, nas quais a expansão intermolar foi maior.

Estudo do efeito das cargas de nutrientes no reservatório APM-Manso: uma abordagem através da modelagem numérica tridimensional.

Os gestores de recursos hídricos estão encarregados da gestão de longo prazo, a regulação e a proteção dos recursos hídricos. No entanto, reconhece-se que a estes gestores devem levar em conta a multiplicidade de usos dos recursos hídricos que são apresentadas pelas partes interessadas, tais como agricultores, fornecedores de água e grupos de ambientalistas. Dada a complexidade do sistema hidrológico, o desenvolvimento e a utilização de modelos matemáticos são muitas das vezes necessários. Neste contexto a modelagem ambiental é frequentemente realizada para avaliar os impactos da degradação do ecossistema devido à ação humana. Essa aplicação orientada a investigações proporciona um importante meio pelo qual os cientistas podem interagir e influenciar nas políticas a nível local, regional, nacional e internacional. No Mato Grosso, durante a implantação da hidroelétrica de Aproveitamento Múltiplo de Manso foram adotadas medidas de mitigação dos impactos socioeconômicos causados. Essas medidas geram uma tendência de aumento populacional associado a uma mudança das características sócio-econômicas, para toda a região do entorno do Reservatório, o que agrava o problema de poluição por nutrientes e denota que existe uma necessidade proeminente de estudos do impacto que estas cargas causariam no ecossistema do reservatório. Utilizando o modelo hidrodinâmico e termodinâmico tridimensional ELCOM, associado ao modelo biogeoquímico Caedym, este trabalho tem a finalidade de criar uma modelagem representativa das cargas dos principais nutrientes causadores da eutrofização cultural, sendo eles: a amônia (NH4), o nitrato (NO3) e o Ortofosfato (PO4), com a finalidade de estudar os efeitos das dinâmicas espaciais e temporais destas cargas no estado trófico deste reservatório no em torno dos pontos de lançamento de esgoto e na sua totalidade.

Variação temporal de uma comunidade fitoplanctônica do reservatório de APM-Manso através de modelagem ecológica tridimensional.

A ecologia de reservatórios, que são ecossistemas complexos, dinâmicos e artificiais, vem assumindo destaque no Brasil. O objetivo deste trabalho foi avaliar a viabilidade da aplicação, no reservatório de APM-Manso, de um modelo ecológico tridimensional em estudos sobre a dinâmica fitoplanctônica, simulando a variação temporal do fitoplâncton para cenários distintos de carga de nutrientes. O modelo CAEDYM foi acoplado ao ELCOM e simulação foi realizada em duas etapas: uma hidrodinâmica e outra ecológica. Escolheu-se para as simulações o período de cinco meses, a partir de 1 de setembro de 2005. Foram construídos dois cenários de simulação, o primeiro contendo os valores reais de carga de nutrientes dos principais rios contribuintes medidos em campo, e o segundo com redução na carga nutricional destes rios, simulando um possível processo de substituição de áreas florestadas por áreas de pastagem na bacia do rio Manso. A comunidade fitoplanctônica simulada apresentou rápidas respostas à disponibilidade nutricional do ambiente, e os resultados obtidos corroboraram com diversas teorias sobre as estratégias adaptativas e sobre as dinâmicas algais. Dentre as classes simuladas, Bacillariophyceae e Cryptophyceae se mostraram mais sensíveis às reduções de carga, enquanto Chrolophyceae e Cyanophyceae, apesar de terem suas biomassas reduzidas, sofreram menos com o impacto, sugerindo estarem mais adaptadas à limitação de nutrientes. Os picos chuvosos influenciaram positivamente as taxas de crescimento das Bacillariophyceae apenas no Cenário 1, uma vez que a limitação por nutrientes foi mais decisiva para esta classe no Cenário 2. Observou-se em ambas as simulações uma tendência de substituição na dominância de Cyanophyceae por Chlorophyceae.

Caracterização de imagens de microtomografia de raios X baseada em descritores de textura

A microtomografia computadorizada (computed microtomography - μCT) permite uma análise não destrutiva de amostras, além de possibilitar sua reutilização. A μCT permite também a reconstrução de objetos tridimensionais a partir de suas seções transversais que são obtidas interceptando a amostra através de planos paralelos. Equipamentos de μCT oferecem ao usuário diversas opções de configurações que alteram a qualidade das imagens obtidas afetando, dessa forma, o resultado esperado. Nesta tese foi realizada a caracterização e análise de imagens de μCT geradas pelo microtomógrafo SkyScan1174 Compact Micro-CT. A base desta caracterização é o processamento de imagens. Foram aplicadas técnicas de realce (brilho, saturação, equalização do histograma e filtro de mediana) nas imagens originais gerando novas imagens e em seguida a quantificação de ambos os conjuntos, utilizando descritores de textura (probabilidade máxima, momento de diferença, momento inverso de diferença, entropia e uniformidade). Os resultados mostram que, comparadas às originais, as imagens que passaram por técnicas de realce apresentaram melhoras quando gerados seus modelos tridimensionais.

Foto-assemblage

FOTO-ASSEMBLAGE consiste em nomenclatura sugestionada para definir os trabalhos que tenho produzido a partir da junção de fotografias digitais. As elaborações e fundamentações desses trabalhos representam também o cerne das pesquisas que resultaram na presente dissertação. Em princípio, o termo foto-assemblage haveria de referir-se a questões técnicas ou formais dessa prática. Contudo, ao desenvolver as pesquisas alguns procedimentos acabaram por determinar certas nuances que revelaram aspectos comuns também em seus e conteúdos. Como resultado de construções artísticas juntando fotografias desde 2009, cheguei às composições sintéticas aqui apresentadas, construídas a partir de duas fotografias. Aventei o nome foto-assemblage por observar nas imagens resultantes ressalvas que as distinguiriam de certas convenções atribuídas à ideia de fotografia. Ao mesmo tempo, as referidas imagens proporiam um possível desdobramento ao entendimento de assemblage enquanto técnica artística. Ainda que não seja uma regra, fotografias revelam imagens de momentos. Em sua relação com a compreensão humana de tempo ou espaço, fotografias quase sempre contêm instâncias mínimas. Fotografias, contudo, podem ser também compreendidas como uma contração de um percurso de tempo. Toda imagem fotográfica pode ser assimilada como resultante de determinados acontecimentos anteriores e mesmo tida como elemento gerador de conseqüências futuras. Seguindo esse entendimento, o que proponho com a foto-assemblage é que essa lide com segmentos de tempo ou de espaço contidos numa mesma imagem. Essas fotografias originárias ganhariam uma nova atribuição, sendo retiradas de seu contexto original, serviriam de balizas do percurso de tempo ou espaço suprimido e subjetivado entre elas. Poeticamente, eventos ocorridos entre as fotografias originárias estariam contidos nas imagens produzidas. O termo assemblage foi incorporado às artes a partir de 1953, por Jean Dubuffet, para descrever trabalhos que seriam algo mais do que simples colagem. A ideia de assemblage se baseia no princípio de que todo e qualquer material ou objeto colhido de nosso mundo cotidiano pode ser incorporado a uma obra de arte, criando um novo conjunto, sem que perca seu sentido original. Esse objeto é arrancado de seu uso habitual e inserido num novo contexto, tecendo laços de relação com os demais elementos, construindo narrativas num novo ambiente, o da obra. Na ideia da foto-assemblage, entretanto, é sugerido uso das imagens fotográficas originárias não como objetos que estariam em um mundo cotidiano, mas sim como imagem na concepção do que seria uma entidade mental. Adoto como que uma visão mágica onde as imagens originárias e básicas estariam numa outra dimensão, num plano bidimensional, não manipulável por nós habitantes da tridimensionalidade. Nesse ambiente imaginário ou não, as fotografias são assentadas consolidando a foto-assemblage. Quando a foto-assemblage se concretiza, se corporifica numa mídia, sendo impressa para uma contemplação, ai então, passaria a integrar nosso mundo tridimensional. O resultado poderia ser admitido como um híbrido, uma terceira coisa, a partir de duas que já não se dissociam mais no ensejo de uma compreensão estética. Ao final da dissertação, apresento experiências práticas que resultaram em quatro séries de imagens em foto-assemblage. Cada série enfatiza aspectos peculiares do que denomino paisagem expandida, representando percursos de tempo, espaço ou trajetos entre o mundo concreto e mundos do inconsciente.

A influência da espessura óssea na medição da altura da borda alveolar em imagens de tomografia computadorizada de feixe cônico e espiral

Imagens de tomografia computadorizada (TC) permitem a visualização, sem distorções ou sobreposições, do complexo maxilo-facial, principalmente do osso alveolar. Estudos demonstraram boa reprodutibilidade e precisão da mensuração da altura da borda alveolar, todavia a influência da espessura óssea ainda é pouco descrita. Através da comparação com a mensuração direta, o objetivo deste estudo foi avaliar a precisão, reprodutibilidade e a influência da espessura óssea, na mensuração da altura da borda alveolar em imagens volumétricas e imagens bidimensionais multiplanares em TC de feixe cônico (TCFC) e em TC espiral (TCE). Utilizando 10 mandíbulas secas de humanos, 57 dentes anteriores foram tomografados em equipamentos iCAT (Imaging Science International, Hatfield, PA, EUA) e Brilliance 64 canais (Philips Eletronics, Eindhoven, Holanda), ambos utilizando voxels de 0,25 mm. Através de imagens volumétricas (3D) e imagens bidimensionais (2D) de cortes multiplanares, foi comparada a mensuração da altura da borda alveolar dessas imagens com a mensuração direta nas mandíbulas, feita por vestibular e lingual, por três avaliadores, com o auxílio de um paquímetro, totalizando 114 bordas alveolares medidas. Alta reprodutibilidade intra-avaliador (0,999 a 0,902) e interavaliador (0,998 e 0,868) foi observada através do índice de correlação intraclasse (ICC). Observou-se alta correlação entre a mensuração direta e indireta da altura da borda alveolar em imagens 2D, sendo r=0,923** e 0,916**, e em imagens 3D, com r=0,929** e 0,954*, em TCFC e TCE, respectivamente. Imagens 2D superestimam a altura da borda alveolar em 0,32 e 0,49 mm e imagens 3D em 0,34 e 0,30 mm, em TCFC e TCE respectivamente. Quando o osso alveolar apresenta espessura de no mínimo 0,6 mm a média da diferença entre medidas diretas e indiretas é de 0,16 e 0,28 mm em imagens 2D e de 0,12 e 0,03 mm em imagens 3D para TCFC e TCE respectivamente, sendo que 95% do limite de concordância varia de -0,46 a 0,79 mm e -0,32 a 0,88 mm em imagens 2D, e de -0,64 a 0,67 mm e -0,57 a 0,62 mm em imagens 3D, para TCFC e TCE respectivamente. Quando o osso alveolar é mais fino do que 0,6 mm a TC é imprecisa, pois 95% do limite de concordância variou de -1,74 a 5,42 mm e -1,64 a 5,42 mm em imagens 2D, e de -3,70 a 4,28 mm e -3,49 a 4,25 mm em imagens 3D, para TCFC e TCE respectivamente. Conclui-se que a mensuração da altura da borda alveolar através de imagens tomográficas apresenta alta reprodutibilidade, sendo que quando a borda alveolar apresenta pelo menos 0,6 mm, a precisão da mensuração é alta, todavia quando esta espessura é menor do que 0,6 mm a técnica é imprecisa.

Avaliação tridimensional da sínfise mandibular em crânios secos e sua associação com diferentes classificações de tipos faciais

O objetivo deste trabalho foi desenvolver um estudo morfológico quantitativo e qualitativo da região da sínfise mandibular (SM), através da construção de modelos tridimensionais (3D) e avaliar o seu grau de associação com diferentes classificações de padrões faciais. Foram avaliados 61 crânios secos humanos de adultos jovens com oclusão normal, com idade entre 18 e 45 anos e dentadura completa. Tomografias computadorizadas de feixe cônico (TCFC) de todos os crânios foram obtidas de forma padronizada. O padrão facial foi determinado por método antropométrico e cefalométrico. Utilizando o critério antropométrico, tomando como referência o índice facial (IF), o padrão facial foi classificado em: euriprósopo (≤84,9), mesoprósopo (85,0 - 89,9) e leptoprósopo (≥90,0). Pelo critério cefalométrico, o plano mandibular (FMA) determinou o padrão facial em curto (≤21,0), médio (21,1 - 29,0) e longo (≥29,1); e o índice de altura facial (IAF) classificou a face em hipodivergente (≥0,750), normal (0,749 - 0,650) e hiperdivergente (≤0,649). A construção de modelos 3D, representativos da região da SM, foi realizada com o auxílio do software ITK-SNAP. Os dentes presentes nesta região, incisivos, caninos e pré-molares inferiores, foram separados do modelo por técnica de segmentação semi-automática, seguida de refinamento manual. Em seguida, foram obtidos modelos 3D somente com o tecido ósseo, possibilitando a mensuraçãodo volume ósseo em mm3 (VOL) e da densidade radiográfica, pela média de intensidade dos voxels (Mvox). No programa Geomagic Studio 10 foi feita uma superposição anatômica dos modelos 3D em bestfit para estabelecer um plano de corte padronizado na linha média. Para cada sínfise foi medida a altura (Alt), a largura (Larg) e calculado o índice de proporção entre altura e largura (PAL). A avaliação da presença de defeitos alveolares foi feita diretamente na mandíbula,obtendo-se a média de todas as alturas ósseas alveolares (AltOss) e a média da dimensão das deiscências presentes (Medef). O índice de correlação intra-classe (ICC) com valores entre 0,923 a 0,994,indicou alta reprodutibilidade e confiabilidade das variáveis medidas. As diferenças entre os grupos, determinados pelas classificações do padrão facial (IF, FMA e IAF), foram avaliadas através da análise de variância (oneway ANOVA) seguida do teste post-hoc de Tukey. O grau de associação entre o padrão facial e as variáveis Vol, Mvox, PAL, Alt, Larg, AltOss e Medef foi avaliado pelo coeficiente de correlação de Pearson com um teste t para r. Os resultados indicaram ausência de diferença ou associação entre o volume, densidade radiográfica e presença de defeitos alveolares da SM e o padrão facial quando determinado pelo IF, FMA e IAF. Verificou-se tendência de SM mais longas nos indivíduos com face alongada, porém a largura não mostrou associação com o padrão facial. Estes resultados sugerem que as classificações utilizadas para determinar o padrão facial não representam satisfatoriamente o caráter 3D da face humana e não estão associadas com a morfologia da SM.

Contribuição ao estudo do efeito tridimensional de instalação e de grupo em estacas cravadas em areia.

O efeito de instalação de estacas cravadas em areia, que promove um aumento da compacidade do solo e de seus parâmetros de resistência e deformabilidade, é analisado nesta dissertação em conjunto com o efeito de grupo. O procedimento estabelecido objetiva quantificar a melhoria após a instalação de um grupo de estacas e foi desenvolvido nesta pesquisa a partir da contribuição de Alves (1998). Para a utilização deste procedimento, foram propostas correlações do módulo de Young do solo em função da resistência à penetração NSPT (ou da resistência à penetração normalizada N60), a partir de um banco de dados obtido da literatura, visando não apenas as aplicações implementadas neste estudo, como também nortear futuras aplicações práticas. Na análise de casos documentados da literatura, foi utilizado o programa Plaxis 3D Foundation. Inicialmente foi procedida a calibração do programa através da reprodução de um caso de instrumentação de um grupo de estacas na argila rija de Londres, muito bem documentado por Cooke et al. (1980). Os resultados das simulações, com parâmetros do solo obtidos da literatura, se ajustaram muito bem aos resultados experimentais, seja nos valores de recalque, seja nos diagramas de transferência de carga. Posteriormente, a melhoria do solo por efeito da instalação foi estimada, através do procedimento estabelecido nesta dissertação, baseado em Alves (1998), a partir das características iniciais do solo antes da instalação de um grupo de estacas em modelo reduzido em solo arenoso. Obtidos os parâmetros melhorados previstos para o solo após a instalação, estes foram aplicados a um grupo de 9 estacas, para diferentes espaçamentos relativos, como dados de entrada no programa Plaxis. A comparação dos resultados experimentais com aqueles simulados numericamente sinaliza para o potencial tanto do programa Plaxis, na reprodução do comportamento do conjunto, como do modelo de Alves (1998), que norteou o procedimento proposto nesta pesquisa. O fator de escala, ressaltado na dissertação, bem como outros aspectos relevantes observados nas análises realizadas, são propostos como temas futuros a serem investigados na continuidade desta linha de pesquisa.

Comparação tridimensional dos efeitos dentários de duas mecânicas para fechamento de espaço: estudo preliminar

As más oclusões de Classe I com biprotrusão dentária são caracterizadas pela projeção dos incisivos superiores e inferiores e podem ser corrigidas com extrações de primeiros pré-molares, permitindo reposicionar os incisivos, usando como ancoragem os dentes posteriores. Esse trabalho teve como objetivo comparar, por meio de superposição de modelos tridimensionais, a movimentação de molares e de incisivos com duas técnicas para fechamento de espaços. Foram selecionados, tomografados (T0) e distribuídos, aleatoriamente em dois grupos, seis pacientes com indicação de exodontia de quatro primeiros pré-molares. No grupo 1 (n=3) foi feito o fechamento dos espaços das extrações em uma única etapa (fechamento em massa) e no grupo 2 (n=3), em duas etapas (distalização de caninos seguida da retração de incisivos). Após o completo fechamento dos espaços foi solicitada nova tomografia (T1). Os modelos tridimensionais em T0 e T1 para os dois grupos foram superpostos para descrever e quantificar a movimentação ocorrida nos primeiros molares. A retração de incisivos foi avaliada de maneira bidimensional. Os molares superiores e inferiores do grupo 1 sofreram movimentação mesial de translação, enquanto os do grupo 2, apresentaram deslocamento de coroa maior do que de ápice radicular, demonstrando movimento de inclinação mesial de coroa. Ao se comparar o movimento mesial dos molares entre os dois grupos, observou-se que no grupo 1 houve maior movimentação mesial da raiz mésio-vestibular dos molares superiores (p=0,009). No entanto, em ambos os grupos os molares superiores giraram sobre a raiz palatina. No arco inferior verificou-se maior movimento mesial de coroa para o grupo 2 (p=0,015). A análise da movimentação dos incisivos apresentou variações semelhantes para as medidas lineares e angulares dos incisivos superiores nos dois grupos. Já, para as medidas angulares dos incisivos inferiores, observou-se que o grupo 1 apresentou menor controle de torque. O tempo de tratamento de todos os pacientes da amostra foi superior a 30 meses, exceto um paciente que apresentou tempo de tratamento inferior a 24 meses. Pôde-se concluir que a perda de ancoragem ocorreu por movimento de translação no fechamento de espaços em massa, e por inclinação mesial no fechamento de espaços em duas etapas, para ambos os arcos. Parece haver maior movimentação mesial radicular nos molares superiores no fechamento em massa e maior movimento mesial de coroa dos molares inferiores para o fechamento em duas etapas. Sugere-se que não há diferença na retração dos incisivos superiores e maior retroinclinação dos incisivos inferiores no fechamento de espaços em duas etapas.

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.