Sem resumo / Sem abstract

O projeto enfocou a aplicação da difração de raios-X (DRX) na caracterização e quantificação dos compostos cristalinos do clínquer de cimento Portland, através do método de Rietveld, constituindo-se em contribuição pioneira sobre o tema em âmbito nacional. Foram utilizadas 40 amostras de clínquer provenientes de cinco diferentes unidades fabris, visando ampla representatividade do material de estudo. O clínquer de cimento Portland é o material sinterizado e peletizado, resultante da calcinação de uma mistura adequada de calcário e argila e, eventualmente, de componentes corretivos. Os compostos metaestáveis do clínquer Portland podem ser subdivididos em três grupos distintos: os silicatos cálcicos (C3S e C2S), a fase intersticial (C4AF, C3A, C12A7), e o grupo dos componentes menos freqüentes como o periclásio, a cal livre e os sulfatos. As proporções destes compostos são parâmetros importantes no controle de processo industrial de clinquerização. O método de Rietveld tem por base a simulação de todo o perfil difratométrico a partir de parâmetros estruturais das fases constituintes, permitindo refinar parâmetros de natureza instrumental e cristalográfica. A comparação do difratograma calculado com o observado e redução das diferenças através do método de mínimos quadrados permitem a obtenção de resultados quantitativos. DRX-Rietveld apresentou-se como uma técnica de quantificação de elevada reprodutibilidade com vantagens de cunho técnico e logístico com relação aos dois métodos correntemente utilizados no Brasil (microscopia e cálculo potencial de Bogue). Técnicas analíticas adicionais permitiram comparar resultados quantitativos obtidos por DRX-Rietveld e também correlacionar características dos compostos com o seu perfil difratométrico. A microscopia óptica foi a técnica de maior importância para comparações tanto qualitativas como quantitativas. A microscopia eletrônica (MEV-EDS) permitiu a aferição de fases não identificadas por microscopia óptica, como o C12A7 e sulfatos. A técnica de dissolução seletiva, aplicada para a concentração da fase intersticial, deu suporte para a aferição quantitativa de teores de polimorfos do C3A. A técnica mostrou resultados coerentes com a microscopia e o cálculo potencial de Bogue, sendo que constituintes não quantificados nestes puderam ser introduzidos (C12A7 e sulfatos), bem como possibilitou a distinção entre polimorfos de um mesmo composto (C3A). A redução no tempo de análise e a diminuição da subjetividade das análises, face às metodologias usuais, constituem fatores importantes da técnica visando atender os interesses da indústria do cimento. / The project has focused on the application of X-ray diffraction (XRD) on the characterization and quantification of the Portland cement clinker crystalline compounds using the Rietveld method. The present research represents a pioneer scientific contribution on the theme in Brazil. Overall forty clinker samples from five distinct kiln lines were collected for analysis aiming to get a broad representativeness of various cement process parameters. Portland cement clinker is the sintered and pelletized product from calcination of an adequate mix of limestone and clay and minor corrective materials. The metastable Portland clinker compounds are subdivided into three main groups: calcium silicates (C3S and C2S), matrix (C4AF, C3A, C12A7) and minor components as periclase, free lime and sulfates. The proportioning of these phases are important parameters to the industry clinkering process. The Rietveld XRD method is based on the simulation of the whole diffraction spectrum from the components structural data, allowing for refining instrumental and crystallographic parameters. By comparing the calculated and actual diffractograms and minimizing differences mathematically through a least squares method quantitative values are obtained. The Rietveld XRD has shown to be a high reproducible quantification technique, with technical and logistics advantages in comparison to the more usual microscopy and Bogue potential calculation. Additional analytical techniques have given reference data to compare quantitative results obtained from Rietveld XRD and to correlate characteristics of the compounds with their diffractogram profile. Optical microscopy was the most relevant technique for comparison both qualitative and quantitatively. Scanning electronic microscopy - energy dispersive system has allowed recognition of phases that could not otherwise be identified by optical microscopy, like C12A7 and sulfates. Selective dissolution of silicates carried out in order to concentrate matrix compounds sustained the quantitative results of C3A polimorphs found by Rietveld XRD. Rietveld has shown coherent results with both microscopy and Bogue potential calculation, but additionnally made it possible quantifying other compounds like C12A7 and sulfates, as well as distinguishing C3A polimorphs. The significant time saving and subjectivity minimization it provides makes up key-factors for the cement industry needs.

Com a popularização de equipamentos tradicionais de captura de imagens, como câmeras digitais, e o avanço tecnológico dos dispositivos não invasivos, como tomografia e ressonância, cresce também a necessidade e consequente uso de métodos numéricos para simulação de fenômenos físicos em domínios definidos por imagens. Um dos pré-requisitos para a aplicação de tais métodos numéricos consiste na discretização do domínio em questão, num processo denominado geração de malhas. Embora diversos métodos de geração de malha tenham sido propostos para discretizar domínios definidos por primitivas geométricas, pouco tem sido feito no sentido de gerar uma decomposição diretamente a partir de imagens. Neste trabalho, apresentamos uma abordagem de geração de malhas de qualidade a partir de domínios definidos por imagens com textura. Mais especificamente, a pesquisa descrita nesta tese contribui com a melhoria do algoritmo Imesh, ao sanar três de suas principais limitações: tratamento de imagens com texturas; controle do nível de refinamento da malha e suporte a outros tipos de elementos. Estas contribuições flexibilizam o processo de geração da malha, e ainda ampliam o domínio de aplicações do algoritmo Imesh, à medida que são considerados domínios definidos por imagens com textura e o uso de métodos numéricos para elementos não simpliciais torna-se possível. O algoritmo de melhoria da malha gerada utiliza uma abordagem inovadora de remalhamento baseada em templates e guiada por retalhos de Bézier / With the spreading of traditional image capturing devices, such as digital cameras, and the technological advancement of more specific imaging devices such as CT and MRI, also increased the need and the following use of numerical methods for simulation of physical phenomena in domains defined by images. One of the prerequisites for the application of such numerical methods is the discretization of the corresponding domain, in a process called mesh generation. Although several methods of mesh generation have been proposed to discretize domains defined by geometric primitives, little has been done to generate a decomposition directly from images. We present an approach to generate quality meshes from domains defined by images with texture. More specifically, the research described in this thesis contributes to the improvement of the Imesh algorithm, removing three of its main limitations: treatment textured images, control of the level of the mesh refinement and support for other types of non-simplicial elements. These contributions provide flexibility to the mesh generation process, and extend the range of applications of Imesh by both handling textured images and considering the use of numerical methods for non-simplicial elements. The mesh quality improvement algorithm uses a new approach based on mesh templates and it is guided by Bezier patches

[pt] Imprescindível para quaisquer aplicações de visão computacional ou realidade aumentada, a calibração de câmera é o processo no qual se obtém os parâmetros intrínsecos e extrínsecos da câmera, tais como distância focal, ponto principal e valores que mensuram a distorção ótica da lente. Atualmente o método mais utilizado para calibrar uma câmera envolve o uso de imagens de um padrão planar em diferentes perspectivas, a partir das quais se extrai pontos de controle para montar um sistema de equações lineares cuja solução representa os parâmetros da câmera, que são otimizados com base no erro de reprojeção 2D. Neste trabalho, foi escolhido o padrão de calibração aneliforme por oferecer maior precisão na detecção dos pontos de controle. Ao aplicarmos técnicas como transformação frontal-paralela, refinamento iterativo dos pontos de controle e segmentação adaptativa de elipses, nossa abordagem apresentou melhoria no resultado do processo de calibração. Além disso, propomos estender o modelo de otimização ao redefinir a função objetivo, considerando não somente o erro de reprojeção 2D, mas também o erro de colinearidade 2D. / [en] Crucial for any computer vision or augmented reality application, the camera calibration is the process in which one gets the intrinsics and the extrinsics parameters of a camera, such as focal length, principal point and distortions values. Nowadays, the most used method to deploy the calibration comprises the use of images of a planar pattern in different perspectives, in order to extract control points to set up a system of linear equations whose solution represents the camera parameters, followed by an optimization based on the 2D reprojection error. In this work, the ring calibration pattern was chosen because it offers higher accuracy on the detection of control points. Upon application of techniques such as fronto-parallel transformation, iterative refinement of the control points and adaptative segmentation of ellipses, our approach has reached improvements in the result of the calibration process. Furthermore, we proposed extend the optimization model by modifying the objective function, regarding not only the 2D reprojection error but also the 2D collinearity error.

[pt] A otimização topológica tem como objetivo encontrar a distribuição mais eficiente de material (ótima topologia) em uma determinada região, satisfazendo as restrições de projeto estabelecidas pelo usuário. Na abordagem tradicional atribui-se uma variável de projeto, constante, denominada densidade, para cada elemento finito da malha. Dessa forma, a qualidade da representação dos novos contornos da estrutura depende do nível de discretização da malha: quanto maior a quantidade de elementos, mais bem definida será a topologia da estrutura otimizada. No entanto, a utilização de malhas super-refinadas implica em um elevado custo computacional, principalmente na etapa de solução numérica das equações de equilíbrio pelo método dos elementos finitos. Este trabalho propõe uma nova estratégia computacional para o refinamento adaptativo local de malhas utilizando elementos finitos poligonais em domínios bidimensionais arbitrários. A ideia consiste em realizar um refinamento da malha nas regiões de concentração de material, sobretudo nos contornos internos e externos, e um desrefinamento nas regiões de baixa concentração de material, como por exemplo, nos furos internos. Desta forma, é possível obter topologias ótimas, com alta resolução e relativamente baixo custo computacional. Exemplos representativos são apresentados para demonstrar a robustez e a eficiência da metodologia proposta por meio de comparações com resultados obtidos com malhas super-refinadas e mantidas constantes durante todo o processo de otimização topológica. / [en] Topology optimization aims to find the most efficient distribution of material (optimal topology) in a given domain, subjected to design constraints defined by the user. The quality of the new boundary representation depends on the level of mesh refinement: the greater the number of elements in the mesh, the better will be the representation of the optimized structure. However, the use of super refined meshes implies in a high computational cost, especially regarding the numerical solution of the linear systems of equations that arise from the finite element method. This work proposes a new computational strategy for adaptive local mesh refinement using polygonal finite elements in arbitrary two-dimensional domains. The idea is to perform a mesh refinement in regions of material concentration, mostly in inner and outer boundaries, and a mesh derefinement in regions of low material concentration such as the internal holes. Thus, it is possible to obtain optimal topologies with high resolution and relatively low computational cost. Representative examples are presented to demonstrate the robustness and efficiency of the proposed methodology by comparing the results obtained herein with the ones from the literature where super refined meshes are held constant throughout all topology optimization process.

A simulação computacional de escoamentos incompressíveis multifásicos tem avançado continuamente e é uma área extremamente importante em Dinâmica de Fluidos Computacional (DFC) por suas várias aplicações na indústria, em medicina e em biologia, apenas para citar alguns exemplos. Apresentamos modelos matemáticos e métodos numéricos tendo em vista simulações computacionais de fluidos bifásicos newtonianos e viscoelásticos (não newtonianos), em seus regimes transiente e estacionário de escoamento. Os ingredientes principais requeridos são o Modelo de Um Fluido e o Método da Fronteira Imersa em malhas adaptativas, usados em conjunto com os métodos da Projeção de Chorin-Temam e de Uzawa. Tais metodologias são obtidas a partir de equações a derivadas parciais simples as quais, naturalmente, são resolvidas em malhas adaptativas empregando métodos multinível-multigrid. Em certas ocasiões, entretanto, para escoamentos modelados pelas equações de Navier-Stokes (e.g. em problemas onde temos altos saltos de massa específica), tem-se problemas de convergência no escopo destes métodos. Além disso, no caso de escoamentos estacionários, resolver as equações de Stokes em sua forma discreta por tais métodos não é uma tarefa fácil. Verificamos que zeros na diagonal do sistema linear resultante impedem que métodos de relaxação usuais sejam empregados. As dificuldades mencionadas acima motivaram-nos a pesquisar por, a propor e a desenvolver alternativas à metodologia multinível-multigrid. No presente trabalho, propomos métodos para obter explicitamente as matrizes que representam os sistemas lineares oriundos da discretização daquelas equações a derivadas parciais simples que são a base dos métodos de Projeção e de Uzawa. Ter em mãos estas representações matriciais é vantajoso pois com elas podemos caracterizar tais sistemas lineares em termos das propriedades de seus raios espectrais, suas definições e simetria. Muito pouco (ou nada) se sabe efetivamente sobre estes sistemas lineares associados a discretizações em malhas compostas bloco-estruturadas. É importante salientarmos que, além disso, ganhamos acesso ao uso de bibliotecas numéricas externas, como o PETSc, com seus pré-condicionadores e métodos numéricos, seriais e paralelos, para resolver sistemas lineares. Infraestrutura para nossos desenvolvimentos foi propiciada pelo código denominado ``AMR2D\'\', um código doméstico para problemas em DFC que vem sendo cuidado ao longo dos anos pelos grupos de pesquisa em DFC do IME-USP e da FEMEC-UFU. Estendemos este código, adicionando módulos para escoamentos viscoelásticos e para escoamentos estacionários modelados pelas equações de Stokes. Além disso, melhoramos de maneira notável as rotinas de cálculo de valores fantasmas. Tais melhorias permitiram a implementação do Método dos Gradientes Bi-Conjugados, baseada em visitas retalho-a-retalho e varreduras da estrutura hierárquica nível-a-nível, essencial à implementação do Método de Uzawa. / Numerical simulation of incompressible multiphase flows has continuously of advanced and is an extremely important area in Computational Fluid Dynamics (CFD) because its several applications in industry, in medicine, and in biology, just to mention a few of them. We present mathematical models and numerical methods having in sight the computational simulation of two-phase Newtonian and viscoelastic fluids (non-Newtonian fluids), in the transient and stationary flow regimes. The main ingredients required are the One-fluid Model and the Immersed Boundary Method on dynamic, adaptive meshes, in concert with Chorin-Temam Projection and the Uzawa methods. These methodologies are built from simple linear partial differential equations which, most naturally, are solved on adaptive grids employing mutilevel-multigrid methods. On certain occasions, however, for transient flows modeled by the Navier-Stokes equations (e.g. in problems where we have high density jumps), one has convergence problems within the scope of these methods. Also, in the case of stationary flows, solving the discrete Stokes equations by those methods represents no straight forward task. It turns out that zeros in the diagonal of the resulting linear systems coming from the discrete equations prevent the usual relaxation methods from being used. Those difficulties, mentioned above, motivated us to search for, to propose, and to develop alternatives to the multilevel-multigrid methodology. In the present work, we propose methods to explicitly obtain the matrices that represent the linear systems arising from the discretization of those simple linear partial differential equations which form the basis of the Projection and Uzawa methods. Possessing these matrix representations is on our advantage to perform a characterization of those linear systems in terms of their spectral, definition, and symmetry properties. Very little is known about those for adaptive mesh discretizations. We highlight also that we gain access to the use of external numerical libraries, such as PETSc, with their preconditioners and numerical methods, both in serial and parallel versions, to solve linear systems. Infrastructure for our developments was offered by the code named ``AMR2D\'\' - an in-house CFD code, nurtured through the years by IME-USP and FEMEC-UFU CFD research groups. We were able to extend that code by adding a viscoelastic and a stationary Stokes solver modules, and improving remarkably the patchwise-based algorithm for computing ghost values. Those improvements proved to be essential to allow for the implementation of a patchwise Bi-Conjugate Gradient Method which ``powers\'\' Uzawa Method.

Neste trabalho é apresentado um método de refinamento para especificações de sistemas embarcados, baseado na linguagem de especificação gráfica UML-RT e na linguagem de especificação formal CSP-OZ. A linguagem UML-RT é utilizada para descrever a arquitetura de sistemas de tempo real distribuídos e esses mapeados para uma especificação formal através de CSP-OZ. A linguagem de especificação formal CSP-OZ é a combinação da linguagem orientada a objetos Object-Z e a algebra de processos CSP, que descreve o comportamento de processos concorrentes. O método de refinamento proposto é baseado na integração de dois métodos: o de bi-simulação, para refinar a parte comportamental da especificação descrita por CSP; e o de equivalência de especificações, para refinar as estruturas de dados descritas por Object-Z, permitindo assim que características de orientação a objetos possam ser utilizadas. Com o método proposto é possível refinar especificações e, conseqüentemente, verificá-las com sua implementação. O desenvolvimento desse método é rigoroso, incluindo a definição formal para um metamodelo da UML-RT. Um exemplo detalhado é apresentado no final deste trabalho. / In this work, a method of refinement of embedded systems specifications based on the graphical specification language UML-RT and the formal specification CSP-OZ is introduced. The UML-RT is used to model real time distributed architecture systems and these are mapped onto formal specifications using CSP-OZ. The CSP-OZ formal specification language is a combination of the state-based object oriented language Object-Z and the CSP process algebra that describes behavioral models of concurrent processes. The rationale of the proposed refinement method is twofold, the use of bisimulation to refine the behavioral part and the specification matching algorithm to refine the state-based part, supporting object-oriented characteristics. Using this result, an equivalence between the specification-matching algorithm and simulation rules is showed. Using the proposed method it is possible to refine CSP-OZ specifications and verify them against their implementations. The development of the proposed refinement method is rigorous, including a formal definition for a UML-RT metamodel. A detailed study case is given at the end of this work.

A simulação computacional de escoamentos incompressíveis multifásicos tem avançado continuamente e é uma área extremamente importante em Dinâmica de Fluidos Computacional (DFC) por suas várias aplicações na indústria, em medicina e em biologia, apenas para citar alguns exemplos. Apresentamos modelos matemáticos e métodos numéricos tendo em vista simulações computacionais de fluidos bifásicos newtonianos e viscoelásticos (não newtonianos), em seus regimes transiente e estacionário de escoamento. Os ingredientes principais requeridos são o Modelo de Um Fluido e o Método da Fronteira Imersa em malhas adaptativas, usados em conjunto com os métodos da Projeção de Chorin-Temam e de Uzawa. Tais metodologias são obtidas a partir de equações a derivadas parciais simples as quais, naturalmente, são resolvidas em malhas adaptativas empregando métodos multinível-multigrid. Em certas ocasiões, entretanto, para escoamentos modelados pelas equações de Navier-Stokes (e.g. em problemas onde temos altos saltos de massa específica), tem-se problemas de convergência no escopo destes métodos. Além disso, no caso de escoamentos estacionários, resolver as equações de Stokes em sua forma discreta por tais métodos não é uma tarefa fácil. Verificamos que zeros na diagonal do sistema linear resultante impedem que métodos de relaxação usuais sejam empregados. As dificuldades mencionadas acima motivaram-nos a pesquisar por, a propor e a desenvolver alternativas à metodologia multinível-multigrid. No presente trabalho, propomos métodos para obter explicitamente as matrizes que representam os sistemas lineares oriundos da discretização daquelas equações a derivadas parciais simples que são a base dos métodos de Projeção e de Uzawa. Ter em mãos estas representações matriciais é vantajoso pois com elas podemos caracterizar tais sistemas lineares em termos das propriedades de seus raios espectrais, suas definições e simetria. Muito pouco (ou nada) se sabe efetivamente sobre estes sistemas lineares associados a discretizações em malhas compostas bloco-estruturadas. É importante salientarmos que, além disso, ganhamos acesso ao uso de bibliotecas numéricas externas, como o PETSc, com seus pré-condicionadores e métodos numéricos, seriais e paralelos, para resolver sistemas lineares. Infraestrutura para nossos desenvolvimentos foi propiciada pelo código denominado ``AMR2D\'\', um código doméstico para problemas em DFC que vem sendo cuidado ao longo dos anos pelos grupos de pesquisa em DFC do IME-USP e da FEMEC-UFU. Estendemos este código, adicionando módulos para escoamentos viscoelásticos e para escoamentos estacionários modelados pelas equações de Stokes. Além disso, melhoramos de maneira notável as rotinas de cálculo de valores fantasmas. Tais melhorias permitiram a implementação do Método dos Gradientes Bi-Conjugados, baseada em visitas retalho-a-retalho e varreduras da estrutura hierárquica nível-a-nível, essencial à implementação do Método de Uzawa. / Numerical simulation of incompressible multiphase flows has continuously of advanced and is an extremely important area in Computational Fluid Dynamics (CFD) because its several applications in industry, in medicine, and in biology, just to mention a few of them. We present mathematical models and numerical methods having in sight the computational simulation of two-phase Newtonian and viscoelastic fluids (non-Newtonian fluids), in the transient and stationary flow regimes. The main ingredients required are the One-fluid Model and the Immersed Boundary Method on dynamic, adaptive meshes, in concert with Chorin-Temam Projection and the Uzawa methods. These methodologies are built from simple linear partial differential equations which, most naturally, are solved on adaptive grids employing mutilevel-multigrid methods. On certain occasions, however, for transient flows modeled by the Navier-Stokes equations (e.g. in problems where we have high density jumps), one has convergence problems within the scope of these methods. Also, in the case of stationary flows, solving the discrete Stokes equations by those methods represents no straight forward task. It turns out that zeros in the diagonal of the resulting linear systems coming from the discrete equations prevent the usual relaxation methods from being used. Those difficulties, mentioned above, motivated us to search for, to propose, and to develop alternatives to the multilevel-multigrid methodology. In the present work, we propose methods to explicitly obtain the matrices that represent the linear systems arising from the discretization of those simple linear partial differential equations which form the basis of the Projection and Uzawa methods. Possessing these matrix representations is on our advantage to perform a characterization of those linear systems in terms of their spectral, definition, and symmetry properties. Very little is known about those for adaptive mesh discretizations. We highlight also that we gain access to the use of external numerical libraries, such as PETSc, with their preconditioners and numerical methods, both in serial and parallel versions, to solve linear systems. Infrastructure for our developments was offered by the code named ``AMR2D\'\' - an in-house CFD code, nurtured through the years by IME-USP and FEMEC-UFU CFD research groups. We were able to extend that code by adding a viscoelastic and a stationary Stokes solver modules, and improving remarkably the patchwise-based algorithm for computing ghost values. Those improvements proved to be essential to allow for the implementation of a patchwise Bi-Conjugate Gradient Method which ``powers\'\' Uzawa Method.

Este é o primeiro trabalho que apresenta simulações tridimensionais completamente adaptativas de um modelo de campo de fase para um fluido incompressível com densidade de massa constante e viscosidade variável, conhecido como Modelo H. Solucionando numericamente as equações desse modelo em malhas refinadas localmente com a técnica AMR, simulamos computacionalmente escoamentos bifásicos tridimensionais. Os modelos de campo de fase oferecem uma aproximação física sistemática para investigar fenômenos que envolvem sistemas multifásicos complexos, tais como fluidos com camadas de mistura, a separação de fases sob forças de cisalhamento e a evolução de micro-estruturas durante processos de solidificação. Como as interfaces são substituídas por delgadas regiões de transição (interfaces difusivas), as simulações de campo de fase requerem muita resolução nessas regiões para capturar corretamente a física do problema em estudo. Porém essa não é uma tarefa fácil de ser executada numericamente. As equações que caracterizam o modelo de campo de fase contêm derivadas de ordem elevada e intrincados termos não lineares, o que exige uma estratégia numérica eficiente capaz de fornecer precisão tanto no tempo quanto no espaço, especialmente em três dimensões. Para obter a resolução exigida no tempo, usamos uma discretização semi-implícita de segunda ordem para solucionar as equações acopladas de Cahn-Hilliard e Navier-Stokes (Modelo H). Para resolver adequadamente as escalas físicas relevantes no espaço, utilizamos malhas refinadas localmente que se adaptam dinamicamente para recobrir as regiões de interesse do escoamento, como por exemplo, as vizinhanças das interfaces do fluido. Demonstramos a eficiência e a robustez de nossa metodologia com simulações que incluem a separação dos componentes de uma mistura bifásica, a deformação de gotas sob cisalhamento e as instabilidades de Kelvin-Helmholtz. / This is the first work that introduces 3D fully adaptive simulations for a phase field model of an incompressible fluid with matched densities and variable viscosity, known as Model H. Solving numerically the equations of this model in meshes locally refined with AMR technique, we simulate computationally tridimensional two-phase flows. Phase field models offer a systematic physical approach to investigate complex multiphase systems phenomena such as fluid mixing layers, phase separation under shear and microstructure evolution during solidification processes. As interfaces are replaced by thin transition regions (diffuse interfaces), phase field simulations need great resolution in these regions to capture correctly the physics of the studied problem. However, this is not an easy task to do numerically. Phase field model equations have high order derivatives and intricate nonlinear terms, which require an efficient numerical strategy that can achieve accuracy both in time and in space, especially in three dimensions. To obtain the required resolution in time, we employ a semi-implicit second order discretization scheme to solve the coupled Cahn-Hilliard/Navier-Stokes equations (Model H). To resolve adequatly the relevant physical scales in space, we use locally refined meshes which adapt dynamically to cover special flow regions, e.g., the vicinity of the fluid interfaces. We demonstrate the efficiency and robustness of our methodology with simulations that include spinodal decomposition, the deformation of drops under shear and Kelvin-Helmholtz instabilities.

Neste trabalho é apresentado um método de refinamento para especificações de sistemas embarcados, baseado na linguagem de especificação gráfica UML-RT e na linguagem de especificação formal CSP-OZ. A linguagem UML-RT é utilizada para descrever a arquitetura de sistemas de tempo real distribuídos e esses mapeados para uma especificação formal através de CSP-OZ. A linguagem de especificação formal CSP-OZ é a combinação da linguagem orientada a objetos Object-Z e a algebra de processos CSP, que descreve o comportamento de processos concorrentes. O método de refinamento proposto é baseado na integração de dois métodos: o de bi-simulação, para refinar a parte comportamental da especificação descrita por CSP; e o de equivalência de especificações, para refinar as estruturas de dados descritas por Object-Z, permitindo assim que características de orientação a objetos possam ser utilizadas. Com o método proposto é possível refinar especificações e, conseqüentemente, verificá-las com sua implementação. O desenvolvimento desse método é rigoroso, incluindo a definição formal para um metamodelo da UML-RT. Um exemplo detalhado é apresentado no final deste trabalho. / In this work, a method of refinement of embedded systems specifications based on the graphical specification language UML-RT and the formal specification CSP-OZ is introduced. The UML-RT is used to model real time distributed architecture systems and these are mapped onto formal specifications using CSP-OZ. The CSP-OZ formal specification language is a combination of the state-based object oriented language Object-Z and the CSP process algebra that describes behavioral models of concurrent processes. The rationale of the proposed refinement method is twofold, the use of bisimulation to refine the behavioral part and the specification matching algorithm to refine the state-based part, supporting object-oriented characteristics. Using this result, an equivalence between the specification-matching algorithm and simulation rules is showed. Using the proposed method it is possible to refine CSP-OZ specifications and verify them against their implementations. The development of the proposed refinement method is rigorous, including a formal definition for a UML-RT metamodel. A detailed study case is given at the end of this work.