Este trabalho apresenta uma metodologia para simulação e análise de uma combinação veicular de transporte de carga do tipo bitrem. Inicia-se com o desenvolvimento e calibração de um modelo numérico da suspensão do veículo, que é um item fundamental para a análise dinâmica do semi-reboque. Dessa forma garante-se que esse subsistema importante esteja coerente com a realidade facilitando a calibração do modelo dinâmico do veículo completo. Na segunda etapa, definem-se os pavimentos com base em medições de rugosidade realizadas em estradas nacionais e realiza-se uma simulação de tráfego com objetivo de obtenção dos carregamentos na estrutura do veículo analisado. Através da passagem do veículo real e do modelo numérico por um obstáculo conhecido realiza-se uma calibração. Com os carregamentos obtidos realiza-se uma análise via o método dos elementos finitos. Em virtude da complexidade da estrutura utiliza-se uma metodologia global-local para obtenção dos históricos de tensão nos pontos críticos da estrutura em cada pavimento, que são comparados com medições experimentais através de extensiometria em pavimentos similares. Através dos históricos de tensões realiza-se uma análise de confiabilidade quanto à probabilidade de falha por sobrecarga em cada pavimento. Já a análise de fadiga nesses pontos é realizada através de duas metodologias, uma estatística e a segunda clássica. Na metodologia estatística, foi empregado um modelo de acúmulo de dano baseado na freqüência média da excitação. A metodologia clássica utilizada é a de Wöhler e consiste na contagem de ciclos através do método de Rainflow, correção da tensão média por Goodman e cálculo de dano por Palmgren-Miner. As duas metodologias forneceram indicações similares sobre quais pontos são mais criticamente solicitados. / This work presents the development and application of a methodology for simulation and analysis of a bi-train type vehicular composition. Firstly, the numerical model of the suspension subsystem is built and calibrated using experimental data. This enables an easier calibration of the composition model in the dynamic case. In a second phase, several road tracks are defined through measurements accomplished in typical Brazilian roads and to be used in the dynamic simulations of the vehicle. The objective of these simulations is to evaluate the loads acting on the semi-trailer structure. A comparison with experimental results is carried out for the case of the vehicle passing on a known obstacle. Because of the complexity of the actual structure, a global-local finite element analysis is necessary in order to obtain the stress histories at critical points caused by the road roughness. These results are compared with experimental measurements through strain gages in similar roads. A reliability analysis is carried out to evaluate the failure probability due to overloads in each type of road. The fatigue analysis at critical points is accomplished through two methodologies, one statistic based and other based on a classic method. In the statistical methodology, a model of damage accumulation was used, which is based on the medium frequency of the loading. The classic methodology is the well known Wöhler approach, along with Rainflow cycle counting, Goodman medium stress correction and Palmgren-Miner damage prediction. Both methodologies provided similar results about the points more critically loaded.

Neste trabalho constrói-se uma representação analítica para a soluçãao da equação de cinética espacial de nêutrons em geometria cilíndrica. O domínio estudado e unidimensional e homogêneo. As equações foram resolvidas, primeiramente, para os casos com um grupo de energia e um grupo de precursores de nêutrons atrasados e, posteriormente, com dois grupos de energia e seis de precursores de nêutrons atrasados. A originalidade do trabalho consiste em inserir uma dependência temporal nas seções de choque de absorção e fissão na equação de cinética. A ideia principal para a representação das soluções reside na obtenção da solução das equações com seções de choque e fissão constantes, utilizando resultados conhecidos para a resolução de sistemas de equações diferenciais lineares de primeira ordem. Em seguida, aplicando a ideia de decomposição, e formado um sistema recursivo de maneira que os termos de correções são considerados como termos fontes, os quais são compostos pelas soluções geradas nas etapas anteriores. Para finalizar, apresentam-se simulações numéricas comparando o método do trabalho e a ferramenta numérica ODE15S (MATLAB) nos casos em que as seções de choque variam com o tempo de forma linear, senoidal e degrau. / In this work, we build an analytical representation for the solution of the neutron space kinetics equations in cylinder coordinates. The domain studied is one-dimensional and homogeneous. The equations were rst solved for the cases with one energy group and one group of delayed neutron precursors and further with two energy groups and six groups of delayed neutron precursors. The originality of the work is to insert a time dependence in the absorption and ssion cross sections in the space kinetics equations. The main idea for the representation of solutions lies in obtaining the solution of the equations with constant absorption/ ssion cross sections, using known results for the solution of system of rst-order linear diffe- rential equations. Then, applying the idea of decomposition, a recursive system is formed so that the terms of corrections are considered as source terms, which are composed of the solutions generated in the previous steps. Finally, numerical simu- lations are presented comparing the method of the present work and the numerical tool ODE15S (MATLAB) in cases where the cross sections are linear, sine and step functions of time.

Carbon-dots (C-dots, pontos de carbono) ou Carbon Quantum Dots (CQDs, pontos quânticos de carbono) consistem em uma interessante classe de nanopartículas de carbono, que vêm recebendo grande atenção devido à sua combinação de propriedades como luminescência, excelente biocompatibilidade e dispersabilidade em água. C-dots são nanopartículas quase esféricas, formados por fragmentos de grafeno empilhados, com diâmetros entre 1 à 10 nm, contendo grupos funcionais na superfície. Embora exista um significativo número de estudos experimentais na literatura, ainda há uma ampla discussão sobre os detalhes estruturais destas nanopartículas e seus efeitos nas propriedades óticas apresentadas. Neste trabalho, utilizamos simulações de dinâmica molecular para explorar detalhes da estrutura atômica dos C-dots, propondo estruturas que apresentem características estruturais, morfológicas e de composição conforme observados em experimentos. Para tanto, criamos inicialmente modelos dessas nanopartículas baseados no empilhamento de fragmentos de grafeno funcionalizados nas bordas. Observou-se que após a relaxação estrutural, os fragmentos foram reorientados, levando a um alinhamento turbostrático entre camadas, com espaçamento próximo aos característicos do grafite Em seguida, calculamos a energia de adesão das camadas dos C-dots, para avaliar a estabilidade desses aglomerados, e os resultados mostram que a presença de grupos funcionais oxigenados contribui para uma melhor adesão entre as camadas. Ainda, criamos modelos semelhantes para essas nanopartículas, incluindo átomos próximos as bordas dos nanoflakes utilizando um método hibrido de MC/MD e os resultados mostram reconstrução da superfície devido a formação espontânea de ligações covalentes nas bordas com a presença de pentágonos, hexágonos e defeitos. Usamos também simulações de dinâmica molecular reativa para analisar a formação destas partículas, partindo de fragmentos resultantes da quebra de moléculas orgânicas, conforme observado experimentalmente. Estas simulações mostraram que a estrutura dos C-dots pode ser um pouco mais complexa do que as atualmente propostas, contendo outros tipos de interação entre os nanoflakes de grafeno, que podem afetar significativamente as propriedades do material. / Carbon-dots (C-dots) or Carbon Quantum Dots consist of an interesting class of carbon nanoparticles, which have received great attention due to their combination of properties such as luminescence, excellent biocompatibility and dispersibility in water. C-dots are almost spherical nanoparticles, formed by stacked graphene fragments, with diameters between 1 and 10 nm, containing functional groups on the surface. Although there is a significant number of experimental studies in the literature, there is still a wide discussion about the structural details of these nanoparticles and their effects on the presented optical properties. In this work, we use molecular dynamics simulations to explore details of the atomic structure of C-dots, which are proposed on basis of experimental characterization data available from the literature. We initially created models of these nanoparticles based on the stacking of functionalized graphene fragments at the edges. It was observed that after the structural relaxation, the fragments were reoriented, leading to a turbostratic alignment between layers, with spacing close to the characteristics of the graphite Then, we calculate the adhesion energy of the C-dots layers, to evaluate the stability of these clusters, and the results show that the presence of oxygenated functional groups contributes to a better adhesion between the layers. Furthermore, we have created similar models for these nanoparticles, including atoms near the edges of nanoflakes using a hybrid MC / MD method and the results show the spontaneous surface reconstruction due to the formation of covalent bonds at the edge with the presence of pentagons, hexagons and defects. We also used simulations of reactive molecular dynamics to analyze the initial steps for the formation of these particles, starting from fragments resulting from the breakdown of organic molecules, as observed experimentally. These simulations have shown that the structure of the C-dots may be what more complex than those currently proposed, containing other types of interaction between the graphene fragments, which can significantly affect the properties of the material.

A combustão possuí inúmeras aplicações comerciais e industriais. Consequentemente, o entendimento dos processos de combustão possuí elevado valor comercial em diferentes setores como aviação, transportes, produção de energia, entre outros. Diversos equipamentos utilizam uma chama para realizar a conversão de energia através da queima. Exemplos desses dispositivos são os fornos industriais e as turbinas a gás. Garantir que a chama permaneça acesa por tempo indeterminado é um problema comum na etapa de desenvolvimento desses queimadores. Por este motivo, muitas pesquisas são realizadas para entender os mecanismos envolvidos na estabilização de uma chama. Em particular, o estudo dos fenômenos termoacústicos tem ganhado destaque nos últimos anos, pois verificou-se que mesmo pequenas oscilações acústicas, isto é oscilações de pressão, podem interferir na estabilidade de uma chama e reduzir a eficiência do queimador. Alguns efeitos da interação da chama com o campo acústico ainda não estão totalmente compreendidos e os modelos utilizados para descrever esses sistemas são ainda raros ou restritos. Considerando essa situação, o presente trabalho busca combinar os modelos existentes e propor um equacionamento que possibilite analisar a acústica de chamas de maneira rápida e suficientemente precisa. O número de equações resolvidas numericamente é reduzido através do desacoplamento da parte química do problema e sua solução é feita previamente de forma analítica. O conjunto de equações obtido nesse processo é implementado em uma linguagem de programação e são realizados 3 testes para comparar a precisão dos dados obtidos com os dados presentes na literatura. Verifica-se que é possível obter boa concordância com os dados experimentais mesmo após algumas simplificações e o desacoplamento das equações químicas. Como o número de equações resolvidas numericamente ´e menor, o custo computacional do problema também é reduzido. Acredita-se que o conjunto de equações apresentado possa ser melhorado continuamente e que possa ser aplicado para finalidades comerciais, de pesquisa e de desenvolvimento na ´area de queimadores de alto desempenho. / The combustion has numerous commercial and industrial applications. Consequently, the understanding of the process of combustion has high commercial value in different sectors such as aviation, transport, energy production, among others. Several devices use a flame to perform the conversion of energy by burning. Examples of these devices are industrial burners and gas turbines. A common problem in the design of these burners is ensure that the flame stays on indefinitely. For this reason, research is being carried out to understand the mechanisms involved in stabilizing a flame. In particular, the study of thermoacoustic phenomena has gained prominence in recent years because it was observed that even small acoustic oscillations, i.e., pressure oscillations may interfere the stability of the flame and reduce the flare efficiency. Some effects of the interaction of the flame with the acoustic field are not yet fully understood and the models used to describe these systems are still rare or restricted. Given this situation, this work seeks to combine existing models and propose a set of equations that allows analysing the acoustics of flames in quickly and sufficiently precise manner. The number of equations numerically solved is reduced by decoupling the chemical part of the problem, whose solution is previously obtained analytically. The set of equations obtained in this process is implemented in a programming language and are performed three tests to compare the accuracy of the data with the data in the literature. Good agreement with the experimental data is obtained even after some simplifications and the decoupling of chemical equations. Because the number of equations solved numerically is smaller, the computational cost of the problem is also reduced. It is believed that the presented set of equations can be continually improved and can be applied for commercial applications, research and development of high-performance burners.

A interação fluido estrutura faz parte de uma ampla gama de estudos da atualidade. Em função das interações dinâmicas entre fluido e estrutura, baseadas tanto nas propriedades do fluido e do escoamento quanto nas características estruturais do corpo submerso, poderá ocorrer vibração induzida por vórtices (VIV) e/ou a vibração induzida pelo escoamento (FIV). Neste trabalho a abordagem desse problema é do tipo numérica, mediante o uso do código Incompact3d, acoplado com um modelo dinâmico de oscilador simples. O esquema do Incompact3d permite a solução das equações de Navier-Stokes mediante Simulação Numérica Direta (DNS). O caso específico para o qual aplica-se o modelo numérico é o caso de dois cilindros circulares lisos idênticos alinhados com o escoamento. A validação do código é feita para números de Reynolds (Re) entre 100 e 300, espaçamento entre centro dos cilindros (S=D) entre 1; 5 e 10; 0 e velocidades reduzidas (U0) entre 2 e 14. Os cenários simulados são: cilindros de montante e de jusante fixos (MF JF), cilindros de montante e de jusante livres na direção transversal ao escoamento (MLy JLy) e cilindros de montante e de jusante livres nas duas direções (MLxy JLxy). Devido à carência de estudos comparativos entres estes três cenários, para números de Reynolds da ordem de 102 e espaçamento crítico (S=D 3; 5, que é função do Re), o presente trabalho estuda as estruturas formadas no escoamento, além das vibrações induzidas, a partir de simulações 2D para Re = 200 e 300, S=D = 3; 5 e 2 < U0 < 14, que corresponde à região onde acontecem os fenômenos tanto de ressonância quanto de galloping. Encontra-se que para o cenário MLyJLy a máximas oscilações ocorrem para 5 U0 8 e são maiores para Re = 300. Em termos de amplitudes de vibração, os cenários MLxy JLxy tiveram os maiores deslocamentos relativos, em relação aos MLy JLy. / The fluid-structure interaction makes part of a great number of current studies. As a function of dynamic interactions between fluid and structure, based on fluid and flow properties as on structural characteristics of the submerged body, could occur vortex induced vibration (VIV) and/or flow induced vibration (FIV). This work has a numerical approach to the problem, using Incompact3d code coupled with a simple oscillator dynamic model. The Incompact3d numerical scheme allows to solve the Navier-Stokes equations via Direct Numerical Simulation (DNS ). The coupled numerical code was applied to case of two circular smooth identical cylinders in cross-flow (tandem). The code validation is made for Reynolds numbers between 100 and 300, spacing between center-to-center cylinder (S=D) between 1; 5 and 10; 0 and reduced velocities (U0) between 2 and 14. The conditions simulated are: both cylinders fixed (MF JF), both cylinders free in transversal streamwise direction (MLy JLy) and both cylinders free to oscillate parallel and transverse to the streamwise direction (MLxy JLxy). Due to the absence of comparative studies between these settings, for Reynolds number O(102) and critical spacing (S=D 3; 5, which is a function of Re), the present work studies the flow structures, moreover induced vibrations, via 2D simulations for Re = 200 and Re = 300, S=D = 3; 5 and 2 < U0 < 14, that correspond to the region where resonance and galloping phenomena take place. It is found for MLy JLy escenary that the greatest oscillations occur for 5 U0 8 and they are greater than the correspond to Re = 300. All the MLxy JLxy settings have greater oscillations than the MLy JLy settings.

Um algoritmo, com velocidade e pressão como variáveis primárias e com condição de Neumann para a pressão, visando analisar escoamentos incompressíveis e aplicado ao estudo do problema da cavidade, é apresentado neste trabalho. São utilizados métodos em diferenças central, up-wind e semi-lagrangeano e formulados matricialmente. O campo da pressão é atualizado de maneira direta, através da equação de Poisson, e inicializado por mínimos quadrados. Simulações são feitas com variados números de Reynolds, obtendo-se o deslocamento do vórtice central e a aparição de vórtices secundários e terciários. / An algorithm with velocity and pressure as primary unknowns and with Neumann pressure boundary condition, to analyse incompressible ows is presented in this work, and it is applied to study the driven cavity ow. Central, up-wind and semi-lagrangean nite di erences schemes, formulated in matrix terms, are employed. The pressure eld is updated directly from Poisson's equation and initialized by least-squares. Simulations with several Reynolds numbers are made, obtaining the displacement of the central vortex as well as the development of secondary and terciary vortices.

Tese (doutorado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Mecânica, 2012. / Submitted by Alaíde Gonçalves dos Santos (alaide@unb.br) on 2013-05-22T15:17:01Z No. of bitstreams: 1 2012_AnnaPauladeSousaParenteRodrigues.PDF: 6383616 bytes, checksum: 9b60bc89be69d6f650b784ecb59ea9d3 (MD5) / Approved for entry into archive by Guimaraes Jacqueline(jacqueline.guimaraes@bce.unb.br) on 2013-05-22T15:27:53Z (GMT) No. of bitstreams: 1 2012_AnnaPauladeSousaParenteRodrigues.PDF: 6383616 bytes, checksum: 9b60bc89be69d6f650b784ecb59ea9d3 (MD5) / Made available in DSpace on 2013-05-22T15:27:53Z (GMT). No. of bitstreams: 1 2012_AnnaPauladeSousaParenteRodrigues.PDF: 6383616 bytes, checksum: 9b60bc89be69d6f650b784ecb59ea9d3 (MD5) / Este trabalho apresenta uma metodologia de otimização para geometria de pás de turbinas hidráulicas axiais. O principal objetivo é determinar a geometria de pá que melhor se adaptará às condições de escoamento e funcionamento impostas, melhorando assim o desempenho dessa turbina. A motivação deste trabalho, decorre da necessidade cada vez maior do "melhor" projeto, em um curto período de tempo, fazendo com que esforços para a concepção de novas metodologias de otimização e ciclos de projeto, sejam desenvolvidos. A otimização se dará através da junção das técnicas de otimização heurísticas: Algoritmos Genéticos (AG) e Otimização por Enxame de Partículas (PSO- Particle Swarm Optimization). O principal objetivo dessa junção é a busca por um método de otimização mais robusto e eficiente na determinação da melhor geometria dessas pás. Dessa forma, o desempenho das pés é considerada como função objetivo. Métodos de dinâmica dos fluidos computacional (CFD- Computational Fluid Dynamics) são utilizados como ferramenta para o cálculo desse desempenho. Essa metodologia é aplicada para uma turbina hidrocinética e uma turbina do tipo bulbo. Os resultados provenientes da otimização da turbina hidrocinética são comparados, para diversos pontos de operação, com valores obtidos através de metodologias de otimizações anteriores. O que proporcionou a comprovação do melhor desempenho do algoritmo híbrido AG-PSO, sendo este posteriormente aplicado à otimização de forma das pás de uma turbina bulbo. Os resultados apresentam um ganho na potência (turbina hidrocinética) e rendimento (turbina bulbo) de saída do rotor otimizado via AG-PSO, além de uma melhora no tempo de processamento demonstrando assim a viabilidade do uso dessa metodologia para a geração de geometrias ótimas de pás de rotores axiais. _______________________________________________________________________________________ ABSTRACT / This thesis presents a methodology as optimization method for the blades of axial turbine runners. The aim is to determine the optimal shape for the blade, considering a specific flow condition and operation point, improving this turbine performance. The motivation behind this work is based on search of the "best"design, increasingly, in a short time, resulting in efforts to new optimization methods and design cycles. The optimization method combines two heuristic optimization techniques: Genetics algorithms (GA) and particle swarm optimization (PSO). The addition of these two techniques in one algorithm (GA-PSO) allows an optimization method more robust and efficient for the purpose of optimal runner's shape. The objective function to be optimized is the performance generated by the rotor. To computation of the performance by the turbine blades, a commercial CFD (Computational Fluid Dynamics) code is used. This methodology is applied to a hidrokinetic turbine and a bulb turbine. The performances of the hidrokinetic turbine, for different operating conditions, of the optimized rotor with GA-PSO were compared to results of previous optimization methods. Proving that the hybrid GA-PSO algorithm has the best performance, being applied to optimization of an bulb runner shape. The results showed a gain in power output (hidrokinetic turbine) and efficiency (bulb turbine), beyond a improvement in the the processing time, which demonstrated the viability of the use of GA-PSO algorithm as a good design procedure for shape optimization of axial hydraulic turbines.

Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Elétrica, 2009. / Submitted by Larissa Ferreira dos Angelos (ferreirangelos@gmail.com) on 2010-04-27T19:13:15Z No. of bitstreams: 1 2009_JoaoPauloLeite.pdf: 2236316 bytes, checksum: 0d7b1e65f81c58bbc15f5459be01fefc (MD5) / Approved for entry into archive by Lucila Saraiva(lucilasaraiva1@gmail.com) on 2010-04-27T22:06:53Z (GMT) No. of bitstreams: 1 2009_JoaoPauloLeite.pdf: 2236316 bytes, checksum: 0d7b1e65f81c58bbc15f5459be01fefc (MD5) / Made available in DSpace on 2010-04-27T22:06:53Z (GMT). No. of bitstreams: 1 2009_JoaoPauloLeite.pdf: 2236316 bytes, checksum: 0d7b1e65f81c58bbc15f5459be01fefc (MD5) Previous issue date: 2009-08-03 / Técnicas de comunicação que utilizam múltiplas antenas (MIMO, multiple-input, multiple-output) são uma importante área de pesquisa, que têm atraído considerável atenção. Em especial, os sistemas MIMO multiusuário são capazes de combinar o aumento de capacidade proporcionado por sistemas MIMO e, quando informação sobre o estado de canal (CSI, channel state information) está disponível no transmissor, as propriedades espaciais do canal podem ser exploradas por meio da pré-codificação, de forma que os diversos usuários podem ser multiplexados espacialmente. Entretanto, a natureza variante no tempo do canal rádio móvel não é levada em consideração ao analisar o desempenho dos esquemas de précodificação. Este trabalho apresenta como primeira contribuição uma proposta de pré-codificador para o enlace direto de sistemas MIMO multiusuário que utiliza o conceito de perturbação contínua. Os resultados de simulação mostram que a técnica proposta apresenta desempenho superior ao das técnicas lineares atualmente relatadas na literatura. A segunda contribuição diz respeito à influência da informação desatualizada sobre o estado do canal no desempenho das técnicas de pré-codificação. É proposto um preditor adaptativo de canal que utiliza o algoritmo set-membership affine projection (SM-AP) como forma de compensar a desatualização da CSI causada pela variação temporal do canal de comunicação. Os resultados mostram que o preditor SM-AP, além de exibir baixo custo computacional, apresenta desempenho superior aos algoritmos clássicos NLMS (normalized least mean square) e RLS (recursive least squares. _________________________________________________________________________________________ ABSTRACT / Multiple-input, multiple-output communication techniques are an important area of research, and have attracted considerable attention. In particular, multiuser MIMO systems are capable of offering the higher link capacity of MIMO systems and, when channel state information (CSI) is available at the transmitter, the spatial properties of the channel can be exploited by precoding, so that multiple users can be spatially multiplexed. However, the time-varying nature of the mobile-radio channel is not taken into account when analyzing the performance of the precoding schemes. The first contribution of this work is the proposal of a precoding technique based on the concept of continuous perturbation. Simulation results show that the proposed technique outperforms the linear precoding techniques currently reported in the literature. The second contribution concerns the influence of the outdated CSI in the performance of the precoding techniques. An adaptive channel predictor based on the set-membership affine projection (SM-AP) filtering is proposed as a way to compensate the outdated CSI caused by the temporal variations of the communication channel. The results show that the SM-AP predictor has low computational cost and it has superior performance when compared to the classical algorithms NLMS (normalized least mean square) and RLS (recursive least squares).

Este trabalho apresenta um estudo sobre a viabilidade de utilização do protocolo conservativo CMB para sincronização de simulações distribuídas em diferentes plataformas MIMD com memória distribuída, utilizando granulosidade grossa com poucos processos paralelos. A técnica utilizada para a realização das análises consiste em efetuar coleta de dados na execução da simulação de um grande número de modelos. Para a simulação desses modelos foi desenvolvida, e é apresentada neste trabalho, a extensão funcional ParSMPL, que implementa o protocolo de sincronização CMB. Pode-se subdividir os resultados obtidos com a análise dos dados das simulações em enfoques que apresentam contribuições distintas. O primeiro enfoque avalia como a plataforma adotada e o modelo a ser simulado influenciam no speedup alcançado. Neste caso, define-se quando uma aplicação tem tendência a ser eficiente com o paradigma da simulação distribuída utilizando o protocolo CMB. O segundo enfoque refere-se à necessidade do usuário conhecer a melhor forma de se utilizar a simulação distribuída. A partir das análises realizadas foi estabelecido um conjunto de diretrizes que visam auxiliar o processo de desenvolvimento de uma simulação distribuída conservativa. Com a ajuda dessas diretrizes e a disponibilidade da extensão ParSMPL espera-se auxiliar o usuário a desenvolver aplicações eficientes utilizando o paradigma da simulação distribuída conservativa, sem exigir desse usuário o conhecimento das características e particularidades do protocolo de sincronização utilizado / This thesis approaches a study about the viability of using the CMB conservative protocol for distributed simulation synchronization on different distributed memory M vID platforms and considering coarse granularity and few parailel processes. The technique used to analyze the results comprises data acquisition during simulation execution for a large number of models. The simulation of these models is performed on a special-built distributed simulation environment (ParSMPL) developed and presented in this work, being a CMB synchronization protocol implementation. The results obtained in this work are split according to different views taking to distinct sets of contributions. The first view evaluates the influence of the model and the execution platform on the speedup reached. In this case it is defined when an application can reach efficiency through the adoption of the distributed simulation paradigm using the CMB protocol. The second view refers to the need of the users knowing the best way to make use of distributed simulation. Thus, following the analysis performed in this thesis, it is established a set of procedures to help in the distributed simulation development process adopting the conservative approach. By following the procedures proposed and using the ParSMPL an user can hope precious help in the development of efficient conservative distributed sirnulation programs, without the need of knowing the features and particularities of the CMB protocol

A qualidade de produtos obtidos pelo processo de co-extrusão depende fortemente da estabilidade do escoamento multicamada nas cavidades internas da co-extrusora. Problemas de instabilidades interfaciais levam a distribuições de espessura de camadas não uniforme e, em situações mais drásticas, ao colapso de camadas. A possibilidade de ocorrência deste tipo de defeito pode limitar a escolha dos materiais a serem utilizados e a espessura mínima de cada camada do produto, bem como restringir a janela de processamento para uma determinada combinação de materiais. Na literatura existem diversos trabalhos que analisam a influência das propriedades reológicas dos materiais utilizados na formação destas instabilidades interfaciais. Um dos critérios de analise é o do sinal de TNSD (Total Normal Stress Difference), o qual serve para analisar a possibilidade de ocorrência de instabilidades interfaciais do tipo onda. A utilização deste critério requer o conhecimento das propriedades reológicas elongacionais e de cisalhamento dos materiais envolvidos, bem como a representação do estado de tensões viscoelásticas em cada uma das camadas durante o escoamento. Tendo como foco central a análise da ocorrência de instabilidades interfaciais em co-extrusão, neste trabalho foi realizada uma série de estudos relacionados ao desenvolvimento e aplicação de uma metodologia para realização desta análise em nosso grupo de pesquisa. Estes estudos envolveram os seguintes tópicos: i) adaptação da técnica Fiber-Windup para obtenção de dados de viscosidade elongacional no reômetro rotacional ARES (Advanced Rheometric Expansion System); ii) estimação de parâmetros dos modelos constitutivos utilizados; iii) solução numérica para o estado de tensões em escoamentos multicamada de co-extrusão; iv) análise experimental da ocorrência de instabilidades interfaciais em sistemas bicamadas envolvendo materiais poliestirênicos. A adaptação da técnica Fiber-Windup mostrou reprodutibilidade tendo desempenho comparável ao de outras técnicas. Também foi proposto um método de correção dos dados experimentais, visando melhor utilização destes na estimação de parâmetros de modelos constitutivos. Para a estimação dos parâmetros não-lineares do modelo constitutivo de Leonov multimodal foi utilizado um procedimento de estimação por etapas e com método estocástico de otimização, obtendo-se ajustes satisfatórios dos dados experimentais. Nos experimentos de co-extrusão foram observadas instabilidades interfaciais do tipo onda, e a ocorrência destas analisadas pelo critério do sinal de TNSD, variável obtida da solução do campo de tensões através de simulações computacionais das condições experimentais. O formato das curvas da quantidade TNSD pareceu refletir o comportamento do estado de tensões normais do escoamento, bem como os efeitos das condições de processamento neste. Porém, o critério de sinal de TNSD pareceu não ser aplicável ao caso estudado, o que, no entanto, pode ser conseqüência do método numérico utilizado para obter o campo de tensões. / The quality of thermoplastic products obtained by co-extrusion process strongly depends of multilayer flow stability in the internal co-extrusion cavities. Interfacial instabilities troubles lead up to thickness non-uniformity layers and, in more drastic situations, to collapsed layers. The occurrence possibility of this type of defect can limit the choice of materials to be used, such as restrict processing window to certain materials matching. In the literature there are a lot of papers analyzing the rheological properties influence of chosen materials in interfacial instabilities occurrence. One analyze criterion is the TNSD-signal (Total Normal Stress Difference), which is used to predict the occurrence of wave type interfacial instabilities. To using this criterion, one needs to know elongational and shear rheological properties of the used materials, such as the viscoelastic stress state presentation in each layers flow. Having as main focus the interfacial instabilities occurrence analysis in co-extrusion, in this work was done a series of researches related to development and application of a methodology to doing this analyze in our work group. This researches involved the topics as follow: i) adaptation of the Fiber Windup technique to obtain elongational viscosity data in ARES (Advanced Rheometric Expansion System) rotational rheometer; ii) parameter estimation of the constitutive models used; iii) numerical solution to the stress state in multilayer co-extrusion flows; iv) experimental analyzes of interfacial instabilities occurrence in bi-layers co-extrusion systems employing polystyrene. The Fiber Windup technique adaptation shown reproducibility, having performance comparable to other techniques. Also was proposed a correction method to experimental data, which aim for their better using in parameter estimation of constitutive models. To the estimation of non-linear parameters of the multimodal Leonov constitutive model, was used a step estimation procedure and optimization stochastic method, with which was obtained satisfactory fitting of experimental data. In the co-extrusion experiments were observed the wave type interfacial instabilities, and their occurrence analyzed by the TNSD-signal criterion, variable obtained from the stress field solution by computational simulations of the experimental conditions. The TNSD curves shape seems to show the normal stress state behavior of flow, as well as the processing conditions effects in this. However, the TNSD-signal criterion seems not be applicable in the studied case, what, yet, can be consequence of the numerical method used to obtain the stress field.