Development, implementation and application of electronic structural descriptors to the analysis of the chemical bonding, aromaticity and chemical reactivity

En la literatura sobre mecànica quàntica és freqüent trobar descriptors basats en la densitat de parells o la densitat electrònica, amb un èxit divers segons les aplicacions que atenyin. Per tal de que tingui sentit químic un descriptor ha de donar la definició d'un àtom en una molècula, o ésser capaç d'identificar regions de l'espai molecular associades amb algun concepte químic (com pot ser un parell solitari o zona d'enllaç, entre d'altres). En aquesta línia, s'han proposat diversos esquemes de partició: la teoria d'àtoms en molècules (AIM), la funció de localització electrònica (ELF), les cel·les de Voroni, els àtoms de Hirshfeld, els àtoms difusos, etc. L'objectiu d'aquesta tesi és explorar descriptors de la densitat basats en particions de l'espai molecular del tipus AIM, ELF o àtoms difusos, analitzar els descriptors existents amb diferents nivells de teoria, proposar nous descriptors d'aromaticitat, així com estudiar l'habilitat de totes aquestes eines per discernir entre diferents mecanismes de reacció.

Algorithm for generating defective graphene sheets

An algorithm is presented for the generation of molecular models of defective graphene fragments, containing a majority of 6-membered rings with a small number of 5- and 7-membered rings as defects. The structures are generated from an initial random array of points in 2D space, which are then subject to Delaunay triangulation. The dual of the triangulation forms a Voronoi tessellation of polygons with a range of ring sizes. An iterative cycle of refinement, involving deletion and addition of points followed by further triangulation, is performed until the user-defined criteria for the number of defects are met. The array of points and connectivities are then converted to a molecular structure and subject to geometry optimization using a standard molecular modeling package to generate final atomic coordinates. On the basis of molecular mechanics with minimization, this automated method can generate structures, which conform to user-supplied criteria and avoid the potential bias associated with the manual building of structures. One application of the algorithm is the generation of structures for the evaluation of the reactivity of different defect sites. Ab initio electronic structure calculations on a representative structure indicate preferential fluorination close to 5-ring defects.

Computational modes and grid imprinting on five quasi-uniform spherical C-grids

Currently, most operational forecasting models use latitude-longitude grids, whose convergence of meridians towards the poles limits parallel scaling. Quasi-uniform grids might avoid this limitation. Thuburn et al, JCP, 2009 and Ringler et al, JCP, 2010 have developed a method for arbitrarily-structured, orthogonal C-grids (TRiSK), which has many of the desirable properties of the C-grid on latitude-longitude grids but which works on a variety of quasi-uniform grids. Here, five quasi-uniform, orthogonal grids of the sphere are investigated using TRiSK to solve the shallow-water equations. We demonstrate some of the advantages and disadvantages of the hexagonal and triangular icosahedra, a Voronoi-ised cubed sphere, a Voronoi-ised skipped latitude-longitude grid and a grid of kites in comparison to a full latitude-longitude grid. We will show that the hexagonal-icosahedron gives the most accurate results (for least computational cost). All of the grids suffer from spurious computational modes; this is especially true of the kite grid, despite it having exactly twice as many velocity degrees of freedom as height degrees of freedom. However, the computational modes are easiest to control on the hexagonal icosahedron since they consist of vorticity oscillations on the dual grid which can be controlled using a diffusive advection scheme for potential vorticity.

Hydrological cycle in the Danube basin in present-day and XXII century simulations by IPCCAR4 global climate models

We present an intercomparison and verification analysis of 20 GCMs (Global Circulation Models) included in the 4th IPCC assessment report regarding their representation of the hydrological cycle on the Danube river basin for 1961–2000 and for the 2161–2200 SRESA1B scenario runs. The basin-scale properties of the hydrological cycle are computed by spatially integrating the precipitation, evaporation, and runoff fields using the Voronoi-Thiessen tessellation formalism. The span of the model- simulated mean annual water balances is of the same order of magnitude of the observed Danube discharge of the Delta; the true value is within the range simulated by the models. Some land components seem to have deficiencies since there are cases of violation of water conservation when annual means are considered. The overall performance and the degree of agreement of the GCMs are comparable to those of the RCMs (Regional Climate Models) analyzed in a previous work, in spite of the much higher resolution and common nesting of the RCMs. The reanalyses are shown to feature several inconsistencies and cannot be used as a verification benchmark for the hydrological cycle in the Danubian region. In the scenario runs, for basically all models the water balance decreases, whereas its interannual variability increases. Changes in the strength of the hydrological cycle are not consistent among models: it is confirmed that capturing the impact of climate change on the hydrological cycle is not an easy task over land areas. Moreover, in several cases we find that qualitatively different behaviors emerge among the models: the ensemble mean does not represent any sort of average model, and often it falls between the models’ clusters.

Modelling the reorientation of sea-ice faults as the wind changes direction

A discrete-element model of sea ice is used to study how a 90° change in wind direction alters the pattern of faults generated through mechanical failure of the ice. The sea-ice domain is 400km in size and consists of polygonal floes obtained through a Voronoi tessellation. Initially the floes are frozen together through viscous–elastic joints that can break under sufficient compressive, tensile and shear deformation. A constant wind-stress gradient is applied until the initially frozen ice pack is broken into roughly diamond-shaped aggregates, with crack angles determined by wing-crack formation. Then partial refreezing of the cracks delineating the aggregates is modelled through reduction of their length by a particular fraction, the ice pack deformation is neglected and the wind stress is rotated by 90°. New cracks form, delineating aggregates with a different orientation. Our results show the new crack orientation depends on the refrozen fraction of the initial faults: as this fraction increases, the new cracks gradually rotate to the new wind direction, reaching 90° for fully refrozen faults. Such reorientation is determined by a competition between new cracks forming at a preferential angle determined by the wing-crack theory and at old cracks oriented at a less favourable angle but having higher stresses due to shorter contacts across the joints

Mesh adaptation on the sphere using optimal transport and the numerical solution of a Monge-Ampère type equation

An equation of Monge-Ampère type has, for the first time, been solved numerically on the surface of the sphere in order to generate optimally transported (OT) meshes, equidistributed with respect to a monitor function. Optimal transport generates meshes that keep the same connectivity as the original mesh, making them suitable for r-adaptive simulations, in which the equations of motion can be solved in a moving frame of reference in order to avoid mapping the solution between old and new meshes and to avoid load balancing problems on parallel computers. The semi-implicit solution of the Monge-Ampère type equation involves a new linearisation of the Hessian term, and exponential maps are used to map from old to new meshes on the sphere. The determinant of the Hessian is evaluated as the change in volume between old and new mesh cells, rather than using numerical approximations to the gradients. OT meshes are generated to compare with centroidal Voronoi tesselations on the sphere and are found to have advantages and disadvantages; OT equidistribution is more accurate, the number of iterations to convergence is independent of the mesh size, face skewness is reduced and the connectivity does not change. However anisotropy is higher and the OT meshes are non-orthogonal. It is shown that optimal transport on the sphere leads to meshes that do not tangle. However, tangling can be introduced by numerical errors in calculating the gradient of the mesh potential. Methods for alleviating this problem are explored. Finally, OT meshes are generated using observed precipitation as a monitor function, in order to demonstrate the potential power of the technique.

Determinação dos parâmetros de fluidização do carvão de charqueadas a temperatura ambiente

Foi construído um equipamento que possibilitou a determinação dos parâmetros de fluidização de amostras de granulometrias simples e de misturas binárias, de carvão de Charqueadas com granulonetrias entre 9 e 28 malhas (mesh number), à temperatura ambiente, tendo-se obtido uma correlação para o cálculo da velocidade mínima de fluidização em função do diâmetro das partículas.

Simulação numérica do processo de fusão e solidificação da água no interior de cilindros

Durante a realização de projetos de sistemas de refrigeração assistidos por sistemas de armazenagem de energia, na forma latente, é necessário que se saiba de antemão as taxas de transferência de energia dos mesmos. Na literatura soluções para este problema são escassos e não abordam todas as geometrias utilizadas para o armazenamento de energia latente, inexistindo correlaçõesmatemáticas que permitam aos projetistas prever o comportamento da fusão e da solidificação do material de mudança de fase (MMF) e as taxas de transferência de energia térmica. Na solução destes problemas são usualmente adotadas hipóteses simplificativas para amenizar o grau de complexidade da descrição matemática do problema Uma destas hipóteses constitui-se em desconsiderar o processo advectivo na transferência de calor no MMF, em sua fase líquida. Outra hipótese adotada, quando o MMF é a água pura, consiste em utilizar um comportamento linear para a massa específica ou um comportamento médio, consistindo em um afastamento da realidade. Este trabalho descreve esforços no sentido de aproximar estas soluções cada vez mais da realidade fisica do problema. São apresentadas simulações para a mudança de fase da água pura levando em conta a inversão da massa específica, em geometria polar, na forma transiente, utilizando o método de Volumes Finitos para solução das equações governantes. As soluções apresentadas, otimizadas quanto às malhas espacial e temporal, envolvem condições de contorno de primeira e terceira espécies Como resultado das soluções são apresentados o Nusselt local, ao longo da parede do cilindro além do Nusselt médio. Correlações matemáticas para o número de Nusselt médio durante a fusão e para o volume de MMF fundido, são propostos e apresentados. A evolução geométrica da fronteira de mudança de fase também é avaliada para as duas condições de contorno e diversos números de Rayleigh. Comparações entre o volume fundido considerando apenas o processo difusivo e o volume fundido utilizando a convecção são realizadas. Para a verificação da independ~encia de malha aplica-se o índice de convergência de malha (GCI) . Algumas comparações qualitativas com soluções experimentais disponíveis na literatura são apresentadas.

Imunorreatividade à FMRF-amida no sistema nervoso central e na musculatura pediosa de Megalobulimus oblongus

O mapeamento da imunorreatividade ao neuropeptídio FMRF-amida no SNC e na musculatura pediosa de Megalobulimus oblongus foi realizado com o objetivo de dar continuidade aos estudos que vem sendo realizados nesta espécie, com o interesse em estabelecê-lo como modelo experimental em estudos neurobiológicos. O neuropeptídio FMRF-amida foi o primeiro peptídio nativo de invertebrados a ser purificado e seqüenciado, a partir do extrato de gânglios cardioativos do molusco Macrocallista nimbosa. Ele atuaria sobre os receptores metabotrópicos, induzindo a abertura dos canais iônicos de K+ tipo S (sensíveis à serotonina (5HT)), hiperpolarizando a membrana celular e assim, aumentando o limiar para o desencadeamento do potencial de ação. A distribuição da imunorreatividade à FMRF-amida em M. oblongus foi investigada empregando-se a técnica imunohistoquímica com anticorpo não-marcado de Sternberger (1979). Todos os gânglios nervosos centrais apresentaram neurônios e fibras FLI (FMRF-amide-like immunoreactivity), assim como suas comissuras, conetivos e a maioria dos nervos emergentes destes gânglios. Nos cortes da musculatura pediosa observou-se também uma intensa imunorreatividade a FMRF-amida. O maior número de neurônios FLI foi encontrado nos gânglios cerebrais com seus tamanhos variando de pequeno a grande. O corpo dorsal apresentou grande quantidade de fibras viii imunomarcadas. No mesocérebro foram detectados alguns agrupamentos FLI junto ao neuropilo e à comissura cerebral e neurônios dispersos na sua porção ântero-medial. O pró-cérebro mostrou-se repleto de neurônios globulosos FLI. Vários agrupamentos foram observados nos lobos pedal, pleural e comissural do pós-cérebro. Nos gânglios pedais houve intensa imunomarcação em agrupamentos neuronais localizados principalmente nas regiões posterior e lateral dos gânglios. Na região anterior foi observado um par de neurônios gigantes FLI. Uma distribuição homogênea de neurônios FLI foi observada por toda a extensão dos gânglios pleurais, com a exceção de apenas um neurônio grande localizado no gânglio pleural esquerdo. A diferença no tamanho dos dois gânglios parietais foi constatada também em M. oblongus. No hemigânglio esquerdo, o menor, foi identificado um neurônio gigante (300 µm). O gânglio direito apresentou um número maior de neurônios FLI, de tamanhos variados. O gânglio visceral foi a estrutura que mais apresentou neurônios grandes e gigantes por toda a sua extensão, com imunorreatividades variáveis. O par de gânglios bucais possuía agrupamentos neuronais FLI médios e grandes. Os neurônios gigantes não se mostraram imunorreativos, ou foram apenas levemente imunomarcados no trofospôngio. Na musculatura pediosa foram observadas fibras nervosas de diferentes calibres nas regiões dorsal, medial e ventral e sobre as células musculares. Identificaram-se gânglios nervosos nos pontos de intersecção dos ramos nervosos e malhas nervosas que formavam plexos nervosos nas regiões ventral e dorsal da musculatura pediosa de M. oblongus.

Biologia alimentar da taxocenose de peixes do Rio Silveira (cabeceira do Rio Pelotas), São José dos Ausentes, Rio Grande do Sul, Brasil

A região de São José dos Ausentes (planalto sul-riograndense), onde se situa o rio Silveira pertencente a cabeceira do rio Pelotas (bacia do rio Uruguai), caracteriza-se por uma ictiofauna composta de espécies endêmicas e várias espécies novas ainda não descritas. O conhecimento referente a estrutura trófica da ictiofauna desta região ainda é inexistente, motivo pelo qual este trabalho teve com objetivo caracterizar a biologia alimentar das espécies de peixes encontradas no rio Silveira. O período de amostragem iniciou em julho de 2000 e concluiu-se em julho 2001. Os indivíduos foram capturados mensalmente através de redes de espera (malhas 1,5 cm entrenós) e puçás, em ambientes de corredeiras e remansos. Foram encontradas seis espécies de peixes neste trecho do rio Silveira: Astyanax sp. n. 4 (Characidae), Eurycheilichthys pantherinus (Loricariidae), Hemipsilichthys sp. n. I (Loricariidae), Cichlasoma facetum (Cichlidae), Rhamdia sp. (Pimelodidae) e Bryconamericus sp. n. 1 (Characidae). O conteúdo estomacal foi identificado para cada uma destas espécies, sendo somente analisado através do método de freqüência de ocorrência, composição percentual e pelo cálculo do índice de importância alimentar para as três espécies mais freqüentes e abundantes. Dentre os itens encontrados destacam-se para Astyanax sp. n. 4 matéria vegetal e Ephemeroptera; para Eurycheilichthys pantherinus larvas de Diptera Simulidae e Ephemeroptera; e em Hemipsilichthys sp. n. I detritos, onde se verificou a presença de diatomáceas e clorofítas. O cálculo do coeficiente de sobreposição alimentar não apresentou valor significativo para estas três espécies. As demais espécies apresentaram os seguintes itens em sua dieta, Cichlasoma facetum larvas de Tricoptera, Ephemeroptera, Gastropoda e insetos alóctones; Rhamdia sp. escamas, Gastropoda, larvas de Lepidoptera, Tricoptera, Ephemeroptera, sementes, matéria vegetal, sedimento; e Bryconamericus sp. n. 1 larvas de Diptera (Simulidae, Psychodidae, Chironomidae), larvas de Tricoptera, Ephemeroptera, Diatomáceas, Clorofítas, Diptera adulto, Hymenoptera, Coleoptera. Esta dieta sugere os seguintes hábitos alimentares para Astyanax sp. n. 4 onívoro, Eurycheilichthys pantherinus insetívoro, Hemipsilichthys sp. n. I detritívoro, Cichlasoma facetum insetívoro, Rhamdia sp. onívoro, e Bryconamericus sp. n. 1 insetívoro.

Desenvolvimento de um método de volumes finitos de alta ordem para a simulação de escoamentos de fluidos viscoelásticos

O uso da mecânica de fluidos computacional no estudo de processos envolvendo o escoamento de fluidos poliméricos está cada vez mais presente nas indústrias de transformação de polímeros. Um código computacional voltado a esta função, para que possa ser aplicado com sucesso, deve levar a predições mais próximas possível da realidade (modelagem), de uma forma relativamente rápida e eficiente (simulação). Em relação à etapa de modelagem, o ponto chave é a seleção de uma equação constitutiva que represente bem as características reológicas do fluido, dentre as diversas opções existentes. Para a etapa de simulação, ou seja, a resolução numérica das equações do modelo, existem diversas metodologias encontradas na literatura, cada qual com suas vantagens e desvantagens. Neste tópico se enquadra o trabalho em questão, que propõe uma nova metodologia para a resolução das equações governantes do escoamento de fluidos viscoelásticos. Esta se baseia no método dos volumes finitos, usando o arranjo co-localizado para as variáveis do problema, e na utilização de aproximações de alta ordem para os fluxos médios lineares e não-lineares e para outros termos não lineares que surgem da discretização das equações constitutivas. Nesta metodologia, trabalha-se com os valores médios das variáveis nos volumes durante todo o processo de resolução, sendo que os valores pontuais são obtidos ao final do procedimento via deconvolução. A solução do sistema de equações não lineares, resultante da discretização das equações, é feita de forma simultânea, usando o método de Newton São mostrados então, resultados da aplicação da metodologia proposta em problemas envolvendo escoamentos de fluidos newtonianos e fluidos viscoelásticos. Para descrever o comportamento reológico destes últimos, são usadas duas equações constitutivas, que são o modelo de Oldroyd-B e o modelo de Phan-Thien-Tanner Simplificado. Por estes resultados pode-se ver que a metodologia é muito promissora, apresentando algumas vantagens frente às metodologias convencionais em volumes finitos. A implementação atual da metodologia desenvolvida está restrita a malhas uniformes e, consequentemente, soluções para problemas com geometrias complexas, que necessitam de refinamento localizado da malha, foram obtidas somente para baixos números de Weissenberg, devido a limitação do custo computacional. Esta restrição pode ser contornada, tornando o seu uso competitivo.

A cidade colonial ibero-americana : a malha urbana

Esta dissertação estuda a cidade colonial ibero-americana, a partir de seu traçado, de sua configuração espacial e como produto da milenar tradição urbana do ocidente, em suas variáveis erudita e popular. Foram analisadas as circunstâncias políticas, econômicas e socioculturais que condicionaram os três séculos do período colonial na Ibero-américa e que influenciaram, de um ou de outro modo, o arranjo espacial das cidades. O trabalho identifica os elementos da arquitetura grega, romana, medieval cristã, muçulmana, renascentista, pré-colombiana e barroca que foram naturalmente selecionados, sintetizados e re-elaborados em sua implementação na cidade ibero-americana. Uma tipologia de malhas urbanas é proposta e, a partir da observação e redesenho de plantas urbanas do período colonial, o trabalho analisa e classifica 21 assentamentos produzidos pela colonização espanhola e portuguesa A análise mostra que a cidade colonial ibero-americana constitui de fato um tipo especifico dentro da categoria maior de cidade tradicional, anterior ao movimento moderno. Por ser uma cidade nova, tem implícita na sua gênese a atividade de planejamento. O traçado em malha é o instrumento regulador essencial. Em seu processo de adaptação às determinantes locais, na busca de uma ordem espacial, o traçado em malha passa por diferentes graus de deformação geométrica, o que condiciona a forma e o posicionamento das partes, ou seja, dos elementos da arquitetura urbana - a praça, a rua, o quarteirão, os edifícios singulares e a estrutura predial de tipos recorrentes de edificação -, gerando assim grande diversidade e riqueza de situações espaciais. O trabalho compara as cidades coloniais espanhola e a portuguesa e consta o predomínio das similaridades pelo fato de terem a malha como denominador comum. As diferenças mais relevantes ficam por conta das implantações, das adaptações ao contexto e de outras circunstâncias específicas, não constituindo fator determinante serem elas espanholas ou portuguesas.

Modelagem e simulação de reatores de polimerização em leito fluidizado

O presente trabalho tem como objetivo o desenvolvimento de um modelo dinâmico genérico para reatores de polimerização em leito fluidizado com capacidade de predizer produção, condições de operação e propriedades do produto, como o índice de fluidez e densidade. Este modelo genérico é então aplicado ao caso específico da polimerização de eteno em reatores do processo UNIPOL, utilizando catalisador de óxido de cromo. De modo a obter um modelo validado, os parâmetros cinéticos são estimados utilizando experimentos realizados em uma planta piloto UNIPOL. Também é apresentado um modelo simplificado em relação ao modelo proposto inicialmente, cujo objetivo é remover as interações dos parâmetros das malhas de controle com os parâmetros cinéticos a serem estimados, aproveitando medidas de planta não utilizadas pelo outro modelo. Junto com a parte de estimação de parâmetros é apresentada a metodologia da análise de sensibilidade para o planejamento de experimentos e estimação de parâmetros.

Otimização de pré-formas e matrizes em problemas bidimensionais de forjamento

Este trabalho apresenta uma sistemática para realizar a otimização numérica de pré-formas e de matrizes em problemas de forjamento axissimétricos e em estado plano de deformações. Para este fim, desenvolveu-se um código computacional composto basicamente de três módulos: módulo de pré-processamento, módulo de análise e módulo de otimização. Cada um destes foi elaborado acrescentando rotinas em programas comerciais ou acadêmicos disponíveis no GMAp e no CEMACOM. Um programa gerenciador foi desenvolvido para controlar os módulos citados no processo de otimização. A abordagem proposta apresenta uma nova função objetivo a minimizar, a qual está baseada em uma operação booleana XOR (exclusive or) sobre os dois polígonos planos que representam a geometria desejada para o componente e a obtida na simulação, respectivamente. Esta abordagem visa eliminar possíveis problemas geométricos associados com as funções objetivo comumente utilizadas em pesquisas correlatas. O trabalho emprega análise de sensibilidade numérica, via método das diferenças finitas. As dificuldades associadas a esta técnica são estudadas e dois pontos são identificados como limitadores da abordagem para problemas de conformação mecânica (grandes deformações elastoplásticas com contato friccional): baixa eficiência e contaminação dos gradientes na presença de remalhamentos. Um novo procedimento de diferenças finitas é desenvolvido, o qual elimina as dificuldades citadas, possibilitando a sua aplicação em problemas quaisquer, com características competitivas com as da abordagem analítica Malhas não estruturadas são tratadas mediante suavizações Laplacianas, mantendo as suas topologias. No caso de otimização de pré-formas, o contorno do componente a otimizar é parametrizado por B-Splines cujos pontos de controle são adotados como variáveis de projeto. Por outro lado, no caso de otimização de matrizes, a parametrização é realizada em termos de segmentos de reta e arcos de circunferências. As variáveis de projeto adotadas são, então, as coordenadas das extremidades das retas, os raios e centros dos arcos, etc. A sistemática é fechada pela aplicação dos algoritmos de programação matemática de Krister Svanberg (Método das Assíntotas Móveis Globalmente Convergente) e de Klaus Schittkowski (Programação Quadrática Sequencial – NLPQLP). Resultados numéricos são apresentados mostrando a evolução das implementações adotadas e o ganho de eficiência obtido.

Modelamento paramétrico e geração de malha em superfícies para aplicações em engenharia

Este trabalho apresenta um conjunto de técnicas para a modelagem paramétrica e geração de malhas de superfícies para uso em sistemas de análise e simulações numéricas pelo Método dos Elementos Finitos. Foram desenvolvidos algoritmos para a geração paramétrica de superfícies, para a determinação das curvas de interseções entre superfícies e para a geração de malhas em superfícies curvas e recortadas. Foram implementas linhas e curvas paramétricas básicas, a partir das quais são geradas superfícies paramétricas de vários tipos que proporcionam uma grande flexibilidade de modelamento geométrico. Curvas e superfícies são geradas e manipuladas de forma interativa. São apresentadas técnicas que simplificam a implementação de linhas e superfícies paramétricas. Foi desenvolvido um algoritmo para determinar as curvas de interseção entre superfícies paramétricas, que são utilizadas como linhas de recorte (trimming lines) para obter geometrias complexas e compostas de várias superfícies. O algoritmo desenvolvido emprega técnicas de subdivisão adaptativa, por quadtrees, em função da curvatura local das superfícies. Primeiramente, obtém-se uma aproximação das curvas de interseção no espaço 3D, através da aproximação por triângulos. Os resultados iniciais são refinados e projetados sobre as duas superfícies envolvidas com algoritmos que permitem obter grande precisão. As curvas de interseção finais são mapeadas nos espaços paramétricos das duas superfícies, porém com uma parametrização única, o que facilita a junção com superfícies adjacentes Um algoritmo de geração de malha foi desenvolvido para gerar malhas triangulares de qualidade sobre as superfícies curvas e recortadas. O algoritmo utiliza um processo de subdivisão adaptativa por quadtrees, similar ao utilizado no algoritmo de interseção, para definir tamanhos de elementos em função da curvatura local. Em seguida, aplica-se um algoritmo tipo advancing front para gerar a malha sobre a superfície. Os algoritmos foram implementados e testados em um ambiente gráfico interativo especialmente desenvolvido para este trabalho. São apresentados vários exemplos que comprovam a eficiência das técnicas e algoritmos propostos, incluindo exemplos de matrizes de conformação mecânica para uso com código de análise METAFOR, análise de sensibilidade para otimização de pré-formas e de modelagem de superfícies compostas recortadas com geração de malhas de qualidade, para uso em análise por Elementos Finitos ou como contorno para geração de elementos tridimensionais.