[pt] Dados multidimensionais, oriundos tanto de medições quanto de simulações, são tipicamente apresentados como campos escalares, vetoriais ou tensoriais. A análise da topologia desses campos, em particular das suas singularidades, ajuda a compreender de forma sucinta a estrutura do dado. Para essa análise, é preciso uma fase de edição para separar o ruído de aquisição ou de representação numérica. O presente trabalho propõe técnicas para edição a partir da topologia, no contexto de campos de cruzes. Esses campos aparecem naturalmente a partir de malhas quadrangulares, que evidenciam a estrutura global do campo. Inicialmente é proposto um operador para realizar refinamentos preservando a estrutura quadrangular da malha e controlando o efeito sobre os vértices singulares. Esse operador permite também construir modelos de forma incremental, aumentando a complexidade a cada passo. No caso de campos de cruzes gerais, é analisado o efeito global de suavizações gaussianas repetidas sobre o campo tensorial associado. Isso permite aplicar técnicas de espaços de escala no caso de campos de cruzes. Os pontos singulares geralmente se cancelam ao longo da escala, desde que se use uma caracterização de pontos singulares consistente com a topologia, como proposto nesse trabalho. Para aumentar o grau de controle sobre as singularidades, é proposta uma técnica para a manipulação individual das singularidades. Essa técnica também permite que os pares de singularidades sejam cancelados de forma local. Utilizando estratégias similares à persistência de campos escalares, é possível simplificar progressivamente o campo de cruzes e atingir configurações além do alcance dos espaços de escala gaussiano. / [en] Measurement and simulation of multidimensional data is typically presented as scalar, vector or tensor fields. The topological analysis of those fields, in particular of their singularities, provides a succinct understanding for the structure of the data. Such analysis requires some editing in order to separate the noise due to limitations in the acquisition and numerical representation. This work proposes techniques for such topology-aware editing, in the context of cross fields. Those fields naturally appear from quadrangular meshes, which summarize the global structure of the field. It is first proposed an operator to refine the mesh preserving its quadrangular structure and controlling the effect on singular vertices. This operator further allows constructing models incrementally, increasing complexity step by step. For general cross fields, the global effect of repeated Gaussian smoothing on the associated tensor field is analyzed. This allows applying scale-space techniques to cross fields. The singular points generally cancel along the scale if using a characterization of singular points consistent with the topology, as the one proposed in this work. To better control the singularities, a technique for individually handling singularities is proposed. It allows canceling pairs of singularities in a local way. Using strategies similar to scalar fields persistence, it is possible to progressively simplify the cross field, reaching configurations beyond Gaussian scale-spaces.

TEIXEIRA, Daniel Nascimento. Uma técnica de decomposição a priori para geração paralela de malhas bidimensionais. 2014. 94 f. : Dissertação (mestrado) - Universidade Federal do Ceará, Centro de Ciências, Departamento de Computação, Fortaleza-CE, 2014. / Submitted by guaracy araujo (guaraa3355@gmail.com) on 2016-06-15T19:57:36Z No. of bitstreams: 1 2014_dis_dnteixeira.pdf: 17919971 bytes, checksum: 092ad12b33cf64a31552e6a839a5a5bc (MD5) / Approved for entry into archive by guaracy araujo (guaraa3355@gmail.com) on 2016-06-15T19:58:41Z (GMT) No. of bitstreams: 1 2014_dis_dnteixeira.pdf: 17919971 bytes, checksum: 092ad12b33cf64a31552e6a839a5a5bc (MD5) / Made available in DSpace on 2016-06-15T19:58:41Z (GMT). No. of bitstreams: 1 2014_dis_dnteixeira.pdf: 17919971 bytes, checksum: 092ad12b33cf64a31552e6a839a5a5bc (MD5) Previous issue date: 2014 / This work describes a technique of two-dimensional domain decomposition for parallel mesh generation. This technique works for both distributed and shared memory and has the freedom to use any data structure that manages rectangular regions parallel to the axes to decompose the domain given as input, such as a quaternary tree (quadtree) or a binary space decomposition (bsp), for example. Any process of mesh generation that respects the prerequisites established can be used in the subdomains created, for instance, Delaunay or Advancing Front, among others. This technique is called a priori because the mesh on the interface of the subdomains is generated prior to the their internal meshes. The load estimation for each sub-domain in this work is performed with the aid of a refined quadtree, whose level of refinement guides the creation of edges that are defined from the bounderies of only inner cells. This way of estimate load produces results that accurately represent the number of elements to be generated in each subdomain. That contributes to a good partitioning of the domain, making the mesh generation in parallel be significantly faster than the serial generation. Furthermore, the quality of the generated mesh in parallel is qualitatively equivalent to that generated serially within acceptable limits. / Este trabalho descreve uma técnica de decomposição de domínios bidimensionais para geração em paralelo de malhas. Esta técnica funciona tanto para memória distribuída quanto compartilhada, além de permitir que se utilize qualquer estrutura de dados que gere regiões quadrangulares paralelas aos eixos para decompor o domínio dado como entrada. Pode se utilizar por exemplo, uma árvore quaternária (quadtree) ou uma partição binária do espaço (bsp). Além disso, qualquer processo de geração de malha que respeite os pré-requisitos estabelecidos pode ser empregado nos subdomínios criados, como as técnicas de Delaunay ou Avanço de Fronteira, dentre outras. A técnica proposta é dita a priori porque a malha de interface entre os subdomínios é gerada antes das suas malhas internas. A estimativa de carga de processamento associada a cada subdomínio é feita nesse trabalho com a ajuda de uma quadtree refinada, cujo nível de refinamento orienta a criação das arestas que são definidas a partir da discretização das fronteiras das células internas. Essa maneira de estimar carga produz resultados que representam, com boa precisão, o número de elementos a serem gerados em cada subdomínio. Isso contribui para um bom particionamento do domínio, fazendo com que a geração de malha em paralelo seja significativamente mais rápida do que a geração serial. Além disso, a qualidade da malha gerada em paralelo é qualitativamente equivalente àquela gerada serialmente, dentro de limites aceitáveis.

This work describes a method based on three stages for reconstructing a model from a given set of scanned meshes obtained from 3D scanners. Meshes scanned from different scanner s view points have their representation in local coordinate systems. Therefore, for final model reconstruction, an alignment of the meshes is required. The most popular algorithm for cloud data registration is the ICP algorithm. However, ICP requires an initial estimate of mesh alignment, which is, many times, done manually. To automate this process, this work uses a surface representation called spin-images to identify overlap areas between the meshes and to estimate their alignment. After this initial registration, the alignment is refined by the ICP algorithm, and finally the model is reconstructed using a method called VRIP. / Fundação de Amparo a Pesquisa do Estado de Alagoas / Este trabalho descreve um método baseado em três etapas para reconstrução de modelos a partir de malhas capturadas de scanners 3D. Malhas obtidas a partir de diferentes pontos de visão de um scanner têm sua representação em sistemas de coordenadas local. Portanto, para a reconstrução final do modelo, é necessário realizar um alinhamento dessas malhas, ou registro. O algoritmo mais famoso para realizar registro de nuvens de pontos é o algoritmo ICP. Porém, um dos requisitos desse algoritmo é uma estimativa inicial do alinhamento das malhas, que muitas vezes é feita manualmente. Para automatizar esse processo, este trabalho utiliza descritores spin-image para identificar regiões de sobreposição entre as malhas e estimar seus alinhamentos. Após este registro inicial, o alinhamento é refinado através do algoritmo ICP, e finalmente o modelo é reconstruído usando uma técnica chamada VRIP.

This work has the goal of explore the triangulated mesh visualization problem in order that geometry and texture information could be computed adaptively according to the view position. As only the essential information from the mesh is going to persist, the main applications involved are linked with real-time mesh visualization which concerns a large number of scientific areas. To adapt the mesh geometry, its resolution is being locally changed through the insertion/removal of vertices, edges and faces, view-dependent considering the number of polygons to preserve smooth forms near to the original mesh. In the texture adaptation, the texture map association of the mesh is being adapted to the resolution changes minimizing or eliminating the texture distortion resulting from the collapse of vertices, edges and faces / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Este trabalho trata do problema de visualização de malhas triangulares de modo que a geometria e a textura das mesmas possam ser calculadas adaptativamente de acordo com a posição do orientador. Como apenas as informações necessárias de malha tendem a permanecer, as principais aplicações envolvidas estão relacionadas à visualização de malhas em tempo real, as quais têm sido de grande interesse de pesquisa em diversas áreas científicas. Na adaptação da geometria da malha, sua resolução vai sendo localmente alterada através de inserção/remoção de vértices, arestas e faces, de acordo com a posição do observador, de modo que o número de polígonos preserve uma forma visual suave e bem próxima da malha original. Na adaptação da textura, o mapeamento de textura associado `a malha vai sendo adaptado `as mudança as de sua resolução de modo a reduzir ou eliminar as distorções de textura resultantes do colapso de vértices, arestas e faces

[pt] A exploração segura de reservatórios de petróleo necessita uma boa modelagem numérica dos objetos geológicos da sub superfície, que inclui entre outras etapas: interpretação sísmica e geração de malha. Esta tese apresenta um estudo nessas duas áreas. O primeiro estudo é uma contribuição para interpretação de dados sísmicos, que se baseia na detecção automática de falhas sísmicas usando redes neurais profundas. Em particular, usamos Redes Neurais Convolucionais (RNCs) diretamente sobre mapas de amplitude sísmica, com a particularidade de usar dados sintéticos para treinar a rede com o objetivo final de classificar dados reais. Num segundo estudo, propomos um novo algoritmo para geração de malhas bidimensionais de quadrilaterais para estudos geomecânicos, baseado numa abordagem inovadora do método de quadtree: definimos novos padrões de subdivisão para adaptar a malha de maneira eficiente a qualquer geometria de entrada. As malhas obtidas podem ser usadas para simulações com o Método de Elementos Finitos (MEF). / [en] Safe oil exploration requires good numerical modeling of the subsurface geobodies, which includes among other steps: seismic interpretation and mesh generation. This thesis presents a study in these two areas. The first study is a contribution to data interpretation, examining the possibilities of automatic seismic fault detection using deep learning methods. In particular, we use Convolutional Neural Networks (CNNs) on seismic amplitude maps, with the particularity to use synthetic data for training with the goal to classify real data. In the second study, we propose a new two-dimensional all-quadrilateral meshing algorithm for geomechanical domains, based on an innovative quadtree approach: we define new subdivision patterns to efficiently adapt the mesh to any input geometry. The resulting mesh is suited for Finite Element Method (FEM) simulations.

Esta dissertação tem por objetivo a implementação de um método adaptativo anisotrópico de malhas e seu acoplamento a um programa de simulação de escoamentos de fluidos compressíveis numa faixa ampla de número de Mach utilizando o Método dos Elementos Finitos (Bubnov-Galerkin). Mais especificamente, implementa-se um código para simulação de escoamentos compressíveis utilizando o método das linhas ou direções características (CBS), o qual possui boa performance para uma amplo espectro de números de Mach. A discretização espacial é feita utilizando-se elementos finitos triangulares lineares e tetraédricos lineares para problemas bidimensionais e tridimensionais, respectivamente. A discretização temporal é feita por uma expansão em séries de Taylor. Implementa-se uma adaptação anisotrópica com capacidade de refinamento e desrefinamento do domínio computacional para se obter soluções precisas a baixos custos computacionais, onde a estimativa de erro é conduzida por uma métrica Riemanniana, a qual permite o tratamento contínuo do erro através de um tensor de erro anisotrópico. O código é modificado para possibilitar o uso de arquiteturas de memória compartilhada a fim de aumentar a velocidade de processamento da simulação. Por fim, diversos estudos, tanto para problemas bidimensionais quanto tridimensionais, são realizados, validando a metodologia e realizando estudos adicionais em problemas mais complexos. / The objectives of this dissertation are the implementation of an anisotropic mesh adaptation procedure and its coupling with a compressible flow simulation code with large Mach range using the Finite Element Method (Bubnov-Galerkin). Specifically, a computer code using the Characteristic Based Split method called CBS, which has a good performance for a wide range of Mach numbers, is implemented. The spatial discretization is performed using linear triangular and linear tetrahedral finite elements for two and three dimensional cases, respectively. The temporal discretization is done with Taylor series. An anisotropic mesh adaption technique, based on a Riemann metric to evaluate the errors (which allows continuous error analysis through an anisotropic error tensor evaluation), including mesh refinement and mesh coarsening procedures is also implemented looking for accurate and computationally cheaper solutions. The code is able to use shared memory architecture in order to reduce computer process time. Finally, two-dimensional and three-dimensional studies are performed in order to validate the proposed methodology and additional complex tests are also analysed.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / O estudo do fluxo de água e do transporte escalar em reservatórios hidrelétricos é importante para a determinação da qualidade da água durante as fases iniciais do enchimento e durante a vida útil do reservatório. Neste contexto, um código de elementos finitos paralelo 2D foi implementado para resolver as equações de Navier-Stokes para fluido incompressível acopladas a transporte escalar, utilizando o modelo de programação de troca de mensagens, a fim de realizar simulações em um ambiente de cluster de computadores. A discretização espacial é baseada no elemento MINI, que satisfaz as condições de Babuska-Brezzi (BB), que permite uma formulação mista estável. Todas as estruturas de dados distribuídos necessárias nas diferentes fases do código, como pré-processamento, solução e pós-processamento, foram implementadas usando a biblioteca PETSc. Os sistemas lineares resultantes foram resolvidos usando o método da projeção discreto com fatoração LU por blocos. Para aumentar o desempenho paralelo na solução dos sistemas lineares, foi empregado o método de condensação estática para resolver a velocidade intermediária nos vértices e no centróide do elemento MINI separadamente. Os resultados de desempenho do método de condensação estática com a abordagem da solução do sistema completo foram comparados. Os testes mostraram que o método de condensação estática apresenta melhor desempenho para grandes problemas, às custas de maior uso de memória. O desempenho de outras partes do código também são apresentados. / The study of the water flow and scalar transport in water reservoirs is important for the determination of the water quality during the initial stages of the reservoir filling and during the life of the reservoir. For this scope, a parallel 2D finite element code for solving the incompressible Navier-Stokes equations coupled with scalar transport was implemented using the message-passing programming model, in order to perform simulations of hidropower water reservoirs in a computer cluster environment. The spatial discretization is based on the MINI element that satisfies the Babuska-Brezzi (BB) condition, which provides sufficient conditions for a stable mixed formulation. All the distributed data structures needed in the different stages of the code, such as preprocessing, solving and post processing, were implemented using the PETSc library. The resulting linear systems were solved using the projection method implemented by an approximate block LU factorization. In order to increase the parallel performance in the solution of the linear systems, we employ the static condensation method for solving the intermediate velocity at the vertex and centroid nodes separately. We compare performance results of the static condensation method with the approach of solving the complete system. In our tests the static condensation method shows better performance for large problems, at the cost of an increased memory usage. Performance results for other intensive parts of the code in a computer cluster are also presented.

[pt] A visualização volumétrica é uma importante técnica para a exploração de dados tridimensionais complexos, como, por exemplo, o resultado de análises numéricas usando o método dos elementos finitos. A aplicação eficiente dessa técnica a malhas não-estruturadas tem sido uma importante área de pesquisa nos últimos anos. Há dois métodos básicos para a visualização dos dados volumétricos: extração de superfícies e renderização direta de volumes. Na primeira, iso-superfícies de um campo escalar são extraídas explicitamente. Na segunda, que é a utilizada neste trabalho, dados escalares são classificados a partir de uma função de transferência, que mapeia valores do campo escalar em cor e opacidade, para serem visualizados. Com a evolução das placas gráficas (GPU) dos computadores pessoais, foram desenvolvidas novas técnicas para visualização volumétrica interativa de malhas não-estruturadas. Os novos algoritmos tiram proveito da aceleração e da possibilidade de programação dessas placas, cujo poder de processamento cresce a um ritmo superior ao dos processadores convencionais (CPU). Este trabalho avalia e compara dois algoritmos para visualização volumétrica de malhas não-estruturadas, baseados em GPU: projeção de células independente do observador e traçado de raios. Adicionalmente, são propostas duas adaptações dos algoritmos estudados. Para o algoritmo de projeção de células, propõe-se uma estruturação dos dados na GPU para eliminar o alto custo de transferência de dados para a placa gráfica. Para o algoritmo de traçado de raios, propõe-se fazer a integração da função de transferência na GPU, melhorando a qualidade da imagem final obtida e permitindo a alteração da função de transferência de maneira interativa. / [en] Volume visualization is an important technique for the exploration of threedimensional complex data sets, such as the results of numerical analysis using the finite elements method. The efficient application of this technique to unstructured meshes has been an important area of research in the past few years. There are two basic methods to visualize volumetric data: surface extraction and direct volume rendering. In the first, the iso-surfaces of the scalar field are explicitly extracted. In the second, which is the one used in this work, scalar data are classified by a transfer function, which maps the scalar values to color and opacity, to be visualized. With the evolution of personal computer graphics cards (GPU), new techniques for volume visualization have been developed. The new algorithms take advantage of modern programmable graphics cards, whose processing power increases at a faster rate than the one observed in conventional processors (CPU). This work evaluates and compares two GPU- based algorithms for volume visualization of unstructured meshes: view- independent cell projection (VICP) and ray-tracing. In addition, two adaptations of the studied algorithms are proposed. For the cell projection algorithm, we propose a GPU data structure in order to eliminate the high costs of the CPU to GPU data transfer. For the raytracing algorithm, we propose to integrate the transfer function in the GPU, which increases the quality of the generated image and allows to interactively change the transfer function.

[pt] A geração de malhas geomecânicas de reservatórios ainda é uma tarefa tediosa que consome muito tempo. Para acelerar este processo, soluções que reconstroem analiticamente a geometria do reservatório têm sido propostas, mas essas soluções não são as mais adequadas para modelagem de objetos naturais. Este trabalho propõe uma modelagem discreta para a geometria do reservatório, onde os vértices da malha são posicionados por meio da solução de um problema de otimização quadrático e convexo. O problema de otimização é modelado de forma a garantir que as malhas geomecânicas de saída sejam suaves e que ao mesmo tempo respeitem as restrições do reservatório e dos horizontes presentes. Além disso, a metodologia proposta permite uma implementação eficiente, paralelizável e de baixo consumo de memória. Casos de teste com milhões de variáveis são apresentados para validar essa abordagem. Finalmente, a metodologia proposta neste trabalho para malhas de geomecânica pode ser naturalmente estendida para a modelagem estrutural de sub-superfícies na interpretação sísmica e de restauração geológica. / [en] Geomechanical mesh generation of complex reservoirs remains a tedious task prone to errors. Recently proposed solutions based on analytical reconstruction of the sub-surfaces are not capable to represent all the geometric details of natural objects. This work proposes a discrete model where the mesh vertices are positioned based on a convex quadratic optimization process. The optimization problem seeks to guarantee smooth meshes that conform with prescribed constraints. The resulting mesh therefore respects, as far as possible, the finite volume mesh of the reservoir pay zone and the existing horizons. Finally, the proposed methodology for Geomechanical meshes can be easily extend to model sub-surfaces present in the structural interpretation and geological restauration.

Made available in DSpace on 2014-06-11T19:23:39Z (GMT). No. of bitstreams: 0 Previous issue date: 2005-08-29Bitstream added on 2014-06-13T18:07:17Z : No. of bitstreams: 1 oliveira_kpm_me_ilha.pdf: 4592597 bytes, checksum: 54231cd87b58b957eaaa00ebf38e513b (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Este trabalho trata da elaboração de um aplicativo computacional em Visual Basic capaz de gerar malhas estruturadas e não estruturadas sobre domínios bidimensionais multiplamente conexos. Esta geração deverá ocorrer de modo bastante automático com pouca intervenção do usuário, a qual será efetuada através de uma interface gráfica amigável. Para armazenamento das malhas definiu-se uma estrutura de dados de fácil compatibilidade com aplicativos computacionais baseados no método dos elementos finitos para solução de problemas do tipo convectivo-difusivo. Os tipos de células (elementos finitos) que foram implementadas são: células triangulares lineares e células quadrilaterais quadráticas. Adicionalmente gerou-se malhas bidimensionais para solução de problemas clássicos do tipo convectivo-difusivo, utilizando-se códigos de elementos finitos já desenvolvidos por pesquisadores do grupo de pesquisa. / The mesh generation is needed in many applications of numerical methods such as finite difference, finite volume and finite element methods. In this work the algebraic method has been applied to generate 2D structured and unstructured mesh of quadrilateral and triangular elements by using Visual Basic. Both linear and quadratic elements can be generated. The connectivity, the nodes coordinates and contour nodes can be saved in an automatic way for a posterior use in, for example, a solver of finite element methods. A friendly interface has been developed for easy usage by users. Some tests have been done in applications of convective-diffusive fluid flows problems using solvers previously constructed, based on finite elements methods to demonstrate the capabilities of the mesh generator.