Com o surgimento das novas tecnologias aplicadas aos equipamentos e arranjos geométricos das fontes renováveis de energia, o sistema de aterramento passou a operar com correntes de várias frequências, diferentes da componente fundamental. Estas correntes são provenientes das novas cargas não lineares conectadas e das condições de falta, tornando mais difícil a equipotencialização do sistema de aterramento, devido, principalmente, ao seu tamanho e sua geometria. Neste cenário, esta pesquisa tem como objetivo o estudo do comportamento de um sistema de aterramento operando na presença de correntes com distorção harmônica. Para a modelagem do sistema de aterramento teste, foi utilizada a teoria de linhas de transmissão, considerando parâmetros distribuídos, pela qual o condutor foi modelado em porções ao longo de seu comprimento. Para a modelagem e análise do sistema de aterramento foram utilizados o programa ATP (Alternative Transients Program) e uma ferramenta computacional comercial, com o intuito de validar os valores da resistência da malha de aterramento observados. Pelos resultados encontrados, têm-se argumentos para alertar os projetos futuros de sistemas de aterramento para a necessidade e uma melhor análise das situações operacionais passíveis de ocorrência no sistema elétrico que venham a apresentar uma composição harmônica diferente e não desprezível em relação ao componente fundamental. / With the rise of new technology applied to devices and geometric arrangements related to renewable energy sources, grounding systems have started to operate with currents of multiple frequencies, deviating from its fundamental component. These currents come from the connection of new nonlinear loads and fault conditions, which makes the equipotential bonding harder for the grounding system, mainly due to its size and geometry. In this scenario, this research aims to study the behavior of a grounding system operating in the presence of harmonically distorted currents. The transmission line theory was used to model the testing grounding system, based on distributed parameters, through which the conductor was modeled in increments along its length. For modeling and studying the grounding system, the ATP (Alternative Transients Program) and a commercial software tool were used in order to validate the resistance values observed in the grounding grid. Based on the results found, there are arguments to alert future grounding system projects about the need for a better assessment of operational scenarios likely to occur in the electrical system, which may display a different harmonic composition that is not negligible in relation to the fundamental component.

O estudo do transporte de partículas tem uma importância fundamental em diversas áreas de pesquisas como, por exemplo, na formação de gelo em uma aeronave pois pode afetar a sua sustentação e estabilidade. Tamanha é a preocupação com a segurança de vôo que diversos estudos têm sido realizados, resultando em códigos computacionais como o LEWICE nos Estados Unidos, TRAJICE no Reino Unido, ONERA na França e CANICE no Canadá. No Brasil, um dos estudo é feito pela EMBRAER em parceria com algumas instituições. O objetivo deste trabalho é desenvolver um algoritmo que possa ser empregado na trajetória das partículas, utilizando uma formulação euleriana que elimina a dificuldade da semeadura de partículas específica da formulação lagrangiana na determinação da fração de volume da partícula. O método empregado é dos volumes finitos em malhas não-estruturadas cuja principal chave está na reconstrução mínimos quadrados de alta ordem com restrição nos contornos. O desenvolvimento do trabalho engloba 3 etapas: definição da geometria e geração das malhas; utilização de um solver para o tratamento do escoamento do ar e obtenção do campo de velocidade; implementação e utilização do esquema numérico com reconstrução mínimos quadrados de alta ordem para simular o cálculo da fração de volume com imposição de condições limites apropriadas no contorno do corpo. Os resultados dos testes realizados mostram que o esquema numérico com reconstrução mínimos quadrados pode ser empregado na resolução de equações que apresentam uma região de descontinuidade, como é o caso da região de sombra, reduzindo a largura da banda de difusão numérica e overshoots. / The particle transport study has a fundamental importance in diverse research area like in the icing accretion on an aircraft because that can affect its sustentation and stability. The concern is so big that many researches have been carried through, resulting in computational codes like the LEWICE in the United States, TRAJICE in the United Kingdom, ONERA in France and CANICE in Canada. In Brazil, one of the study has been made by the EMBRAER with some institutes. The goal of this work is to develop an algorithm that can be used in the particles trajectory study, using an Eulerian method that eliminates the difficulty particle sowing, particular of the Lagrangian method, in the determination of the droplet fraction volume. This is made by the finite volume method on unstructured meshes whose main key is the high order reconstruction with restriction on the boundary. The development of the work involves 3 stages: geometry definition and mesh generation; using code for the treatment of the air flow and obtained flow velocity; use of the high order numerical scheme least square reconstruction to simulate the droplet fraction volume result with imposition of appropriate limit conditions in the body contour. The realized simulations shown that Least Square method can be used in problem resolution that present descontinuos region like is shadow region reducing numerical diffusion and overshoots.

Apresenta-se a evolução dos sitemas de comunicação, com ênfase especial nos sistemas com tecnologia óptica. Discute-se a necessidade contínua do aumento de capacidade de tais sistemas de comunicação, e a consequente repercussão sobre os futuros sistemas ópticos. Em vista da necessidade do aumento de capacidade dos futuros sistemas de comunicação óptica, apresentam-se em seguida duas propostas recentes da literatura, sendo uma referente à realização de um gerador de pulsos ópticos estreitos, e a outra referente à implementação de um extrator de relógio realizado com técnicas ópticas. Apresenta-se um breve resumo da teoria das malhas de sincronismo de fase, ou PLLs, mostrando como as duas propostas discutidas realizam funções típicas desses sistemas. Ressalta-se a necessidade dos PLLs de sistemas ópticos possuirem ganhos de malha (parâmetro K) elevados. Após a caracterização dos requisitos necessários para PLLs de futuros sistemas ópticos, e após um resumo de alguns conceitos necessários da teoria de redes e filtros elétricos, apresentam-se dois tipos de projeto de filtros para aqueles PLLs. As duas formas de projeto tem como objetivo viabilizar o uso de filtros de tipo e ordem arbitrários. Um tipo de projeto visa a realização da função de transferência do PLL com a curva de resposta igual à de um filtro escolhido. O outro tipo de projeto parte da realização do filtro de loop do PLL com as características de um tipo de filtro escolhido, e define métodos para ajustar a função de transferência resultante para esse PLL. Apresentam-se algoritmos para o cálculo de parâmetros importantes nos dois procedimentos. Após a discussão dos dois pontos de vista de projeto, apresentam-se exemplos de realização de PLLs de acordo com as técnicas apresentadas. Para cada exemplo, mostram-se as curvas de resposta em frequência tanto do PLL como do correspondente filtro de loop, bem como o lugar das raízes e a resposta de captura do PLL obtido. O processo de captura foi estudado por simulações que procuram reproduzir o mais fielmente possível as condições reais de implementação, sem entretanto considerar efeitos de ruído. Finalmente, mencionam-se brevemente possíveis linhas de pesquisa futuras, sendo o foco principal o uso de filtros com pólos e zeros finitos. / The evolution of communication systems is discussed, with emphasis on optical technology. Special consideration is given to the continuous need for increasing the capacity of such systems, and the impact over future optical communication. In view of the great demands imposed over the capacity of future optical systems, an overview is presented of two recent proposals found in the literature, one of such proposals being the implementation of a generator of short optical pulses, and the other being a clock extractor device realized through the use of optical techniques. A brief review is made of phase-locked loop (or PLL) theory, to show how the discussed proposals could be used to realize tipical functions found in these systems. The very high loop gains (the so-called parameter K) that must be used in PLLs of optical communication systems are emphasized. After discussion of the necessary characteristics for PLLs of future optical systems, and also after a review of some concepts of the theory of electrical networks and filters, two design procedures for filters to be used in such PLLs are presented. Both designs have the goal of allowing the use of loop filters with any type and order. The first type of design has the objective to realize a PLL transfer function that has a frequency response identical to the response of a chosen type of filter. The other design starts with a chosen type of filter for a PLL loop filter, arriving to an suitable PLL transfer function. Some algorithms for determination of important design parameters are also presented. After the discussion of the two types of design, some examples of PLLs obtained by such methods are presented. For each example, frequency response curves are presented for the PLL and the respective loop filter, as well as the root locus and the capture response for the PLL so obtained. The capture process was studied through the use of simulations with parameters intended to approximate real implementation conditions, although noise effects are not considered. Finally, some possible research lines are discussed, whose main focus is on filters with finite poles and zeros.

Made available in DSpace on 2015-03-04T18:50:26Z (GMT). No. of bitstreams: 1 TeseEduardoCoppoliFinal30ago2011.pdf: 3181243 bytes, checksum: 92f18a357bb128aeb133a3854958d1e9 (MD5) Previous issue date: 2011-10-27 / Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior / Métodos sem malha constituem uma classe de métodos numéricos usualmente aplicados na solução de equações diferenciais parciais (PDE). A principal característica destes métodos é que os mesmos não necessitam de uma malha como a usada nos métodos de elementos finitos. Esta característica pode ser bastante útil quando se esta modelando estruturas móveis, situação onde geralmente se necessita de um processo de reconstrução da malha. Este trabalho apresenta soluções para alguns dos principais aspectos encontrados na modelagem de dispositivos eletromagnéticos através do método sem malha denominado Element-Free Galerkin Method. Problemas como descontinuidade de materiais, domínios periódicos, imposição de condições de contorno essenciais, acoplamento entre as equações de campo e de circuito e modelagem de movimento são tratados nesta tese. Problemas de eletromagnetismo com materiais descontínuos constituem uma situação muito comum e são a causa de alguns problemas para os métodos sem malha. Geralmente esta característica leva a oscilações espúrias nas curvas do campo elétrico e magnético. Para se resolver este tipo de problema é proposta a utilização de funções de peso especiais na implementação do método. O uso de funções de peso especiais também permitira que a imposição das vii condições de contorno essenciais seja feita de maneira similar ao Método de Elementos Finitos. Este tipo de função propiciara também a simplificação da formulação variacional e do código computacional. Uma nova abordagem para o tratamento de problemas com fronteiras periódicas é proposta, fazendo uso de funções de peso especiais e na manipulação do sistema de equações. Este trabalho também apresenta um método que permite o acoplamento entre equações de circuito e do campo magnético no domínio do tempo para o EFGM. Para tal, uma maquina de indução trifásica é usada como exemplo, sendo apresentadas suas equações acopladas entre campo magnético e circuito magnético e circuito elétrico. Como etapa final do trabalho uma abordagem de modelagem do movimento desta maquina elétrica é proposta.

O processo de injeção de combustíveis nas câmaras de combustão é de vital importância, sob diversos aspectos, para o funcionamento ideal de cada tipo de motor. O elemento responsável pela introdução de propelentes líquidos e gasosos na câmara de combustão de foguetes e pela transformação de massas líquidas em "sprays" é o injetor que converte a energia potencial dos propelentes em energia cinética, pela queda de pressão em seu interior, formando um jato ou folha líquida que se desintegra então em gotas[

Presented here is the effort concerning the application of Computational Fluid Dynamics (CFD) as a swirl injector analysis tool. A literature review of the major inviscid swirl models is provided along with the viscous correction procedure. The bi-propellant atomizer design process is described and the features of previous designs are also detailed. The study was undertaken by focusing on the 5-[kN]-thrust rocket engine currently in development. Theoretical predictions of the discharge coefficient, spray angle and liquid film thickness were obtained for both the inner and outer swirlers of the core injection elements. The governing equations are solved based on the laminar volume of fluid (VOF) interface capturing method. Results from cold flow experiments and particle image velocimetry (PIV) are compared to the predictions of the swirl models and the numerical results. The laminar VOF model was able to predict the spray angle with reasonable confidence, however, a deviation of 25 % was observed in the mass flow rate and discharge coefficient. Although the laminar VOF model has proven inadequate, it constitutes a good starting point in the procedure needed to assess swirl injector performance.

Filtered two-fluid models used to perform large scale simulations of gas-solid fluidized flows of industrial risers require closures for filtered parameters such as filtered and residual stresses, and interphase interaction forces mainly effective drag. Closure models for those filtered parameters may be derived by averaging over results of highly resolved simulations with microscopic two-fluid modeling. This work is a contribution in that context. Recent models for filtered parameters have been written as functions of filter size, filtered solid volume fraction, and filtered slip velocity. A recent study showed that macro-scale variables like domain average solid volume fraction and gas Reynolds number also significantly affect the filtered parameters. In the current work, in addition to these filtered and macro-scale variables, the effects of two new variables over the filtered parameters are investigated: filtered solid kinetic energy and sub-grid gas turbulence. It is shown that the filtered solid kinetic energy should be accounted for in the concerning correlations, thereby improving accuracy. Regarding gas turbulence, literature shows it has no significant effects on the motion of high Stokes particles. Extending on literature, this work investigates the sub-grid gas turbulence effects on meso-scale structures formed of high Stokes particles. Results showed that sub-grid gas turbulence has no significant effects on the meso-scale structures and corresponding filtered parameters. The open source code MFIX was used for all simulations. Ranges of dilute concentration of solid and gas Reynolds number typical of riser flow regimes were considered. A modified two-fluid model with microscopic formulation was used. The sub-grid gas turbulence was generated by means of a forcing function procedure which was implemented in the physical space, over the gravitational term of the momentum conservation equation of the gas phase. First, numerical simulations of the gas phase alone were performed, accounting for literature available data, in order to set a turbulent gas field and calibrate the turbulence intensity. Then, the forcing function was introduced in to the two-fluid model and various gas-solid flows were simulated. While the current results show the necessity of accounting for additional variables in the filtered parameter correlation, they also make it clear the necessity of further developments that are required in the search for better accuracy. / Modelos filtrados de dois-fluidos usados em simulações de grandes escalas de escoamentos fluidizados de gás-sólido de risers industriais exigem fechamentos para parâmetros filtrados tais como as tensões filtradas e residuais, e forças interativas interfases, principalmente arrasto efetivo. Modelos de fechamento para estes parâmetros filtrados podem ser gerados a partir de procedimentos de media aplicados sobre resultados de simulações altamente resolvidas com modelagem microscópica de dois-fluidos. Este trabalho é uma contribuição neste contexto. Modelos de fechamento recentes para parâmetros filtrados tem sido formulados em função de tamanho de filtro, fração volumétrica de sólido filtrada, e velocidade de deslizamento filtrada. Estudo recente mostrou que variáveis de macro-escalas como fração volumétrica de sólido e número de Reynolds de gás médios no domínio também afetam significativamente os parâmetros filtrados. No presente trabalho, além dessas variáveis filtradas e de macro-escala, os efeitos de duas novas variáveis sobre os parâmetros filtrados são investigados: energia cinética filtrada do sólido e turbulência submalha do gás. Em relação à energia cinética filtrada do sólido, mostra-se que a sua consideração refina as correlações em questão, contribuindo assim para melhor acuracidade. Com relação à turbulência do gás, a literatura mostra que não tem efeitos significativos no movimento de particulados de elevados números de Stokes. Acrescentando à literatura, este trabalho investiga os efeitos da turbulência sub-malha do gás sobre estruturas de meso-escala formados de particulados de elevados números de Stokes. Os resultados mostraram que a turbulência sub-malha do gás não tem efeitos significativos sobre estruturas de meso-escalas e parâmetros filtrados correspondentes. O código aberto MFIX foi usado para todas as simulações. Faixas de concentração diluída de sólido e número de Reynolds típicos de escoamentos em risers foram considerados. Um modelo modificado de dois fluidos com formulação microscópica foi utilizado. A turbulência sub-malha do gás foi gerada por meio de um procedimento de \'forcing function\' que foi implementado no espaço físico, sobre o termo fonte gravitacional da equação de momentum da fase gás. Primeiramente, simulações numéricas da fase gás foram realizadas separadamente, levando-se em conta dados disponíveis na literatura, a fim de gerar um campo de gás turbulento e calibrar a intensidade de turbulência. Posteriormente, a \'forcing function\' foi introduzida no modelo de dois-fluidos e vários escoamentos de gás-sólido foram simulados. Enquanto os resultados obtidos mostram a necessidade de consideração de variáveis adicionais para correlação de parâmetros filtrados, também deixam claro a necessidade de desenvolvimentos mais aprofundados na busca de melhor acuracidade.

O Método dos Elementos Finitos Generalizados, MEFG, compartilha importantes características dos métodos sem malha. As funções de aproximação do MEFG, atreladas aos pontos nodais, são enriquecidas de modo análogo ao refinamento p realizado no Método das Nuvens hp. Por outro lado, por empregar uma malha de elementos para construir as funções partição da unidade, ele também pode ser entendido como uma forma não convencional do Método dos Elementos Finitos. Neste trabalho, ambas as interpretações são consideradas. Os métodos sem malha, particularmente o Método de Galerkin Livre de Elementos e o Método das Nuvens hp, são introduzidos com o propósito de estabelecer os conceitos fundamentais para a descrição do MEFG. Na seqüência, apresentam-se aplicações numéricas em análise linear e evidenciam-se características que tornam o MEFG interessante para a simulação da propagação de descontinuidades. Após discutir os modelos de dano adotados para representar o comportamento não-linear do material, são introduzidos exemplos de aplicação, inicialmente do Método das Nuvens hp e depois do MEFG, na análise de estruturas de concreto. Os resultados obtidos servem de argumento para a implementação de um procedimento p-adaptativo, particularmente com o MEFG. Propõe-se, então a adaptação do Método dos Resíduos em Elementos Equilibrados à formulação do MEFG. Com vistas ao seu emprego em problemas não-lineares, algumas modificações são introduzidas à formulação do estimador. Mostra-se que a medida obtida para representar o erro, apesar de fundamentada em diversas hipóteses nem sempre possíveis de serem satisfeitas, ainda assim viabiliza a análise não-linear p-adaptativa. Ao final, são enumeradas propostas para a aplicação do MEFG em problemas caracterizados pela propagação de defeitos / The Generalized Finite Element Method, GFEM, shares several features with the so called meshless methods. The approximation functions used in the GFEM are associated with nodal points like in meshless methods. In addition, the enrichment of the approximation spaces can be done in the same fashion as in the meshless hp-Cloud method. On the other hand, the partition of unity used in the GFEM is provided by Lagrangian finite element shape functions. Therefore, this method can also be understood as a variation of the Finite Element Method. Indeed, both interpretations of the GFEM are valid and give unique insights into the method. The meshless character of the GFEM justified the investigation of meshless methods in this work. Among them, the Element Free Galerkin Method and the hp-Cloud Method are described aiming to introduce key concepts of the GFEM formulation. Following that, several linear problems are solved using these three methods. Such linear analysis demonstrates several features of the GFEM and its suitability to simulate propagating discontinuities. Next, damage models employed to model the nonlinear behavior of concrete structures are discussed and numerical analysis using the hp-Cloud Method and the GFEM are presented. The results motivate the implementation of a p-adaptive procedure tailored to the GFEM. The technique adopted is the Equilibrated Element Residual Method. The estimator is modified to take into account nonlinear peculiarities of the problems considered. The hypotheses assumed in the definition of the error measure are sometimes violated. Nonetheless, it is shown that the proposed error indicator is effective for the class of p-adaptive nonlinear analysis investigated. Finally, several suggestions are enumerated considering future applications of the GFEM, specially for the simulation of damage and crack propagation

Nas últimas três décadas os phase locked loops (PLLs) totalmente digitais têm recebido muita atenção devido, principalmente, às vantagens que eles oferecem em comparação aos PLLs analógicos. Essas vantagens incluem melhor desempenho, maior velocidade e confiabilidade, tamanho reduzido e menor custo. Os PLLs também são amplamente utilizados em sistemas de comunicações e em outras aplicações digitais. A presente dissertação é uma contribuição no campo dos PLLs digitais adaptativos e otimizados para a sua implementação em hardware. É feito uma análise de suas características dinâmicas e proposta uma nova estrutura de PLL digital capaz de melhorar a resposta da malha em termos de tempo de aquisição e largura de banda. A Malha Síncrona Digital \"Tanlock\" com Estimação de Frequência e Ganho Adaptativo para Convergência Rápida, como é chamada, foi desenvolvida a partir da malha digital \"Tanlock\", utilizando-se teoremas de ponto fixo e mapas contrativos para determinar as condições de ganho que garantam convergência rápida e melhor utilização da largura de banda. Resultados das simulações são comparados com os obtidos teoricamente para avaliar o desempenho da malha proposta. / In the last three decades, fully-digital Phase-Locked-Loops (PLLs) systems have received a lot of attention due to its advantages in comparison with analog PLLs. These advantages include improved transient response, reliability and also reduced size and cost. The PLLs are widely used in communications systems and many other digital applications. This dissertation is a contribution to the field of digital adaptive PLLs optimized to hardware implementation. Here, a new PLL structure is presented; the Frequency Sensing Adaptive TDTL is an improvement to the classic Time-Delay Tanlock structure, alowing fast convergence to the synchronous states, using fixed-point theorems and contractive maps to determine the gain conditions which ensure the rapid convergence and also providing wider bandwidth. The results of simulations are compared with those obtained theoretically in order to assess the loop performance.

Malhas dinâmicas são comumente utilizadas em problemas de simulação sobre dominios cuja geometria varia com o tempo. Sempre que o domínio onde a malha está definida é alterado, as molas são acionadas movimentando os vértices para que estes se conformem com a nova descrição do domínio. Os tipos de molas mais utilizadas são: as longitudinais, as torcionais e as semi-torcionais. Nesta tese uma nova mola é proposta, a mola altura, que além de evitar sobreposição de elementos, é conceitualmente simples e fácel de ser implementada. Outra contribuição desse trabalho é o mecanismo de vértices ativos, que permite economia de processamento durante a resolução da malha dinâmica. Quando a fronteira do domínio sofre grandes alterações, o processo dinâmico pode não ter êxito na correção da malha. Para contornar esse problema, a fronteira deve ser alterada aos poucos. Uma nova estratégia para realizar grandes deformações em pequenos passos é introduzida nesta tese. Em algumas aplicações, o movimento da fronteira da malha pode comprometer células da própria fronteira. A correção da fronteira e um processo delicado, já que em muitos casos ele implica em alterar a descrição do domínio. Um novo método para efetuar a correção da fronteira é apresentado neste trabalho. Ele é baseado em malhas dinâmicas e utiliza um novo conceito de molas, as molas conservativas. Todas as contribuições citadas acima tiveram aplicação prática na industria aeronáutica, sendo utilizadas na implementação de uma metodologia inovadora para acoplar um simulador de escoamento de fluidos tridimensional com uma ferramenta de projeto inverso de aerofólios que roda em um contexto bidimensional. O outro assunto abordado e o remalhamento de triangulações superficiais. Foi proposto um novo método, chamado ANTS (Anisotropic Triangulations on Surfaces) que produz triangulações anisotrópicas de qualidade sobre superfícies descrevendo objetos com geometria complexa. O método ANTS é caracterizado por efetuar o remalhamento diretamente na triangulação inicial, isto é, ele não faz uso de qualquer tipo de parametrização, seja global ou local. O processo de remalhamento é feito por meio de quatro operadores: inserção, remoção e movimento de vértices e alternância de arestas. Os operadores de inserção e remoção de vértices possibilitam controlar a densidade de vértices no domínio, permitindo que nós sejam inseridos em regiões com densidade baixa ou eliminados onde a densidade é alta. A qualidade dos triângulos é controlada por meio dos operadores de movimento de vértices e de alternância (flipping) de arestas. O operador de movimento é utilizado no núcleo do processo de remalhamento. Para evitar que o remalhamento danifique a forma original da superfície, as quinas e os córneres são detectados no inicio do processo e preservados durante o remalhamento. A densidade de vértices sobre o domínio é controlada por uma função de espalhamento. Tal função pode ser passada como entrada para o ANTS ou calculada pelo próprio método. O ANTS foi aplicado com êxito em diversos exemplos gerando malhas de boa qualidade / This thesis intends to make a contribution on the field of dynamic meshes. Dynamic meshes are commonly used in the simulation of problems on domains whose geometry varies in time. Virtual springs are placed in the mesh to rearrange its vertices whenever the domain is changed. The most commonly used types of springs are: longitudinal, torsional and semi-torsional. In this thesis a new type of spring is introduced, the height spring, that is conceptually simple but produces good results. Another contribution of this thesis is the active vertices mechanism, that can improve the CPU processing time of the dynamic mesh. When the mesh domain undergoes large deformations, the proposed dynamic mesh algorithm may fail in correcting the mesh. A solution to this problem is perform large deformations in smal steps. A new strategy for this purpose is presented. Sometimes the motion of the mesh boundary can damage cells on the boundary itself. This is a trick problem to solve since the correction of boundary might change the domain geometry. A new method to correct the boundary cells is also presented in this study. The method is based on the dynamic mesh concept and uses a new type of spring, the conservative spring. All the mentioned contributions had been applied in the aeronautics industry. The techniques developed here has been used to implement an innovative methodology to couple a three-dimensional fluid dynamic solver with a two-dimensional inverse design tool for airfoils. This thesis also deals with remeshing. It is presented the ANTS, a practical method for remeshing anisotropic triangulations on surfaces of complex geometry. The method is capable of performing refinement and coarsening during the same process using the well-known remeshing operators: vertex motion, vertex deletion (by collapsing edges), vertex insertion, and edge flipping. An interesting feature is that vertex motion is used in the core of the process instead of in a post-processing smoothing step. The ANTS uses the input mesh as the geometrical description and works directly on the surface mesh without using any other auxiliary structure (besides the input mesh itself) to preserve the geometrical shape. Moreover, neither global nor local parameterization are applied. Sharp edges and points are identified at the beginning and kept during the process in order to preserve ridges and details. The method has been successfully applied to several examples producing high quality meshes