Este trabalho apresenta simulações numéricas de escoamento de fluidos através de placas de orifício, elementos primários de medição de vazão, em comparação aos resultados esperados por normas internacionais. O coeficiente de descarga usado para a determinação da vazão de um escoamento em uma tubulação foi obtido numericamente através da aplicação do modelo de turbulência de baixo Reynolds proposto por Launder-Sharma, resolvido através de um código aberto de Dinâmica de Fluidos Computacional. Diferentes configurações de placa de orifício foram analisadas, tais como: diâmetro do orifício de restrição, número de Reynolds, posição de instalação, ângulo de chanfro. A principal abordagem deste trabalho é apresentar uma metodologia para determinação do coeficiente de descarga para a montagem invertida da placa de orifício e apresentar os parâmetros que influenciam na determinação da vazão nesta condição. Foi observado que os melhores resultados foram obtidos para os maiores diâmetros de placa de orifício, salientando que todos os resultados apresentaram desvios inferiores a 2% em relação à norma ISO-5167 (2003). Na montagem invertida da placa de orifício foi observada uma redução na diferença de pressão entre as faces da placa de orifício em comparação à montagem correta, resultando em um maior valor necessário para o coeficiente de descarga. De acordo com os resultados deste trabalho a determinação do coeficiente de descarga para a montagem invertida da placa de orifício depende do número de Reynolds, ângulo de chanfro e da dimensão do orifício de restrição. Este trabalho apresenta ainda resultados para a determinação da vazão em condição de montagem invertida da placa, que podem ser abordados em trabalhos futuros para estimativa do coeficiente de descarga em condições adversas de montagem.

This work is part of the research and development, at the Institute for Advanced Studies (IEAv), of the first Brazilian hypersonic vehicle prototype, the 14-X airplane. As this vehicle will be propelled by scramjet (supersonic combustion ramjet) engines, this work presents detailed two-dimensional CFD analyses of the airflow in the intake system of such engines based on the 14-XB scramjet geometry and the expected flight conditions. The main objective is to study the airflow in the intake of the 14-XB at nominal flight condition and also for some off-design flight conditions and geometry using numerical methods and models available in the Fluent code. Off-design values of the vehicle velocity, angle of attack and altitude as well as of the angle of the inlet compression ramp and the number of inlet compression ramps were chosen to show how these changes impact the overall intake airflow. In this study are presented results for the airflow in the entire intake system and of specific flow variables at the engine combustor entrance, as well as calculation results of some intake performance parameters. Both, wall temperature and free stream flow turbulence effects on the intake airflow have also been analyzed. Investigation of viscous flow modeling and of the effects of temperature-dependent air properties has also been performed. Inviscid flow calculations have been performed to serve as a comparison basis for the viscous flow effects and as preliminary information of the airflow. A model validation analysis of the k-kl-? and Transition SST transition models has shown that both models can calculate BL and shock wave interactions (SWBLI) quite well, although, the k-kl-? is better to calculate the separation region whereas the Transition SST is superior to predict the reattachment point. Wall temperature has shown to affect quite significantly SWBLI while viscous flow modeling has shown to have an important impact on the intake airflow with some degradation of the intake system performance.

Na usinagem, visando a processos produtivos com menor geração de resíduos, maior eficiência e mais sustentáveis, são desenvolvidas novas técnicas, como a usinagem com mínima quantidade de lubrificantes (MQL) ou mesmo a seco. A inserção de novas tecnologias também implica investimento financeiro e alterações no consumo de energia. O presente trabalho tem por objetivo avaliar as dimensões técnico-econômica, social e ambiental do fresamento de Ferro Fundido Vermicular (CGI - Cast Grey Iron) utilizando as técnicas de Lubri-refrigeração a seco, MQL e inundação. O ambiente de estudo foi eleito conforme seu nível de criticidade, pois o CGI possui alta taxa de grafita e o fluido de corte de base vegetal de Pinhão-manso (Jatropha curcas) contém alta taxa de matéria orgânica em sua composição. Como critério técnico-econômico estudou-se o consumo energético da máquina-ferramenta e de equipamentos auxiliares específicos, além da quantificação do volume de fluido de corte gasto e seu gerenciamento. Na dimensão ambiental, estudou-se a qualidade do ar (fungos ambientais) e os resíduos gerados. Por fim, na dimensão social, foram estudados indicadores relacionados à saúde do operador. Avaliando a quantidade de óleo presente em suspensão (névoa) no ambiente de trabalho e a poluição sonora gerada com a aplicação de cada técnica. Levando em consideração os aspectos técnico-econômicos, ambientais e sociais analisados, concluiu-se que a técnica de Lubri-refrigeração a seco é a mais indicada para o fresamento de CGI. Os custos envolvidos provaram ser menores para essa condição, que produz resíduos de tratamento mais facilitado. Os resultados também indicam menores problemas de saúde ocupacional. Adicionalmente, conclui-se que o ambiente de trabalho estudado não está apto a trabalhar com o processo de Lubri-refrigeração por MQL, pois ainda faltam projetos de máquinas a serem melhorados.

Nas últimas décadas, as melhorias na tecnologia de extit{hardware} principalmente no que diz respeito a unidade de processamento central, memória e capacidade de armazenamento, diversas ferramentas numéricas aplicadas na engenharia foram desenvolvidas com o intuito de se obter respostas rápidas e mais baratas do que ensaios, experimentos e testes laboratoriais. A Dinâmica dos Fluidos Computacional (DFC) faz uma aproximação numérica como um meio de simular escoamentos. Um dos grandes desafios em soluções numéricas na DFC é a garantia da consistência, estabilidade e convergência garantindo com isso, a confiabilidade do código numérico. As turbomáquinas operam em uma ampla faixa de velocidades no campo do escoamento, podendo ter regiões com baixos números de Mach. Nesse ponto, surge um problema, pois códigos da DFC para o cálculo de escoamento compressível geralmente são programados seguindo formulações matemáticas baseadas no vetor de variáveis conservadas, podendo ter problemas numéricos e apresentar soluções inconsistentes devido a pequena variação da densidade do fluido. Chorin através de alguns estudos constatou que os códigos de DFC possuem uma diminuição de sua convergência e estabilidade quando o número de Mach do escoamento torna-se menor do que o limite de compressibilidade. Com o objetivo de solucionar este problema, Chorin criou a técnica de compressibilidade artificial. Choi e Merkle; Turkel e outros pesquisadores modificaram a técnica de compressibilidade artificial e desenvolveram a técnica de pré-condicionamento para diversas aplicações para garantir a convergência, estabilidade e melhorias na soluções numéricas dos códigos DFC para uma ampla faixa de velocidades do escoamento. No presente trabalho, foi realizada a implementação do método de pré-condicionamento de Turkel, Vatsa e Radaspiel em um programa de DFC desenvolvido por Tomita para calcular escoamentos compressíveis, internos e externos. A fim de melhorar a robustez do código numérico o pré-condicionador foi modificado seguindo o proposto no trabalho de Damorfal e Siu. O pré-condicionamento foi implementado primeiramente para as Equações de Euler e depois para escoamentos viscosos e turbulentos. Diversas experiências numéricas foram realizadas e resultados interessantes utilizando o código de DFC compressível e o código de DFC compressível com o pré-condicionador original e modificado foram obtidos. Todos estes casos estão apresentados e comentados no presente trabalho. A implementação do pré-condicionador em um código de DFC desenvolvido para o cálculo de escoamento tipicamente compressível resultou em uma ferramenta numérica com maior potencial para o cálculo do escoamento numa ampla faixa de número de Mach.

En el presente trabajo se ha realizado el diseño mecánico de una planta piloto de extracción por fluido supercrítico que utiliza dióxido de carbono (CO2) en estado supercrítico como solvente. La planta piloto está diseñada para poder procesar 5 litros de biomateria por cada corrida y trabajar con una presión de hasta 300 bar y una temperatura de 100C con un flujo de CO2 de 1 kg/min, los parámetros de estado han sido determinados a partir de un análisis de la biomateria disponible en el Perú con potencial para beneficiarse de la extracción por fluido supercrítico, en el cual destaca el Sacha Inchi, con alto contenido de Omega-3 como componente minoritario, como el material de referencia debido a que sus parámetros de extracción son los más altos del grupo encontrado (Pext = 300 bar y Text = 60 C), de modo que la planta sea capaz de procesar productos cuyos parámetros de extracción se encuentren por debajo del rango mencionado. El uso de CO2 supercrítico como solvente está justificado debido a la relativa facilidad para alcanzar su condición de fluido supercrítico (Pcrit = 73:8 bar y Tcrit = 31 C) respecto de otros fluidos cuyas presiones y temperaturas críticas son mucho más altas lo cual impacta directamente en el costo de la planta y en el caso de la temperatura posiblemente en la integridad de los extractos obtenidos. Por otra parte, el CO2 es un gas relativamente barato y no genera residuos dañinos ni para el producto deseado ni para el medio ambiente. El presente trabajo está enfocado en el diseño y la selección de los componentes principales de la planta piloto, a partir de lo cual se ha realizado la cotización de la planta obteniendo un costo de US$ 60812.27 incluido IGV. Cabe indicar que el apartado del control y automatización de la planta se ha planteado únicamente a nivel de ingeniería básica por lo cual su desarollo está previsto como trabajo a futuro. / Tesis

Este trabalho aborda o escoamento ao redor de cilindros flexíveis com grande razão de aspecto, representando os riser utilizados na indústria petrolífera. O escoamento foi dividido em seções bidimensionais ao longo do comprimento do riser. O método dos vórtices discretos foi utilizado para estimativa dos coeficientes hidrodinâmicos nas seções correspondentes e o acoplamento se dá através da dinâmica estrutural. Esta dinâmica é resolvida pelo método dos elementos finitos implementado no código Anflex (Mourelle et al. (2001)). Processamento paralelo é utilizado para acelerar o desempenho do método numérico. Um esquema master-slave (\"mestre-escravo\") utilizando MPI (Message Passing Interface1) é utilizado para que seja explorado o paralelismo da modelagem. As seções hidrodinâmicas são igualmente divididas ao longo dos nós de um cluster de computadores utilizado nos cálculos. Cada nó do cluster (núcleo de processamento) resolve o escoamento nas seções hidrodinâmicas necessárias. As forças hidrodinâmicas atuando nas seções correspondentes são enviadas ao processo mestre que também é responsável pela resolução da dinâmica estrutural. Uma das contribuições do presente trabalho é a possibilidade de análise de risers em regime pós-crítico, através da imposição do ponto de separação da camada-limite baseado em resultados experimentais, e da modelagem de casos com risers em tandem. O Método dos Vórtices Discretos (MVD) está detalhado de forma a facilitar a implementação computacional, levando em consideração a possibilidade de estudo de interferência entre esteiras de múltiplos corpos. Como primeira aproximação, foram simulados, casos de cilindros rígidos fixos e em base elástica para os quais foram examinadas as amplitudes máximas de oscilação e os coeficientes de arrasto e sustentação. Posteriormente, o MVD foi integrado ao Anflex dando origem ao AnflexCFD, com capacidade de processamento paralelo da parte correspondente ao CFD. Os resultados obtidos numericamente foram comparados a resultados experimentais realizados em laboratório. A metodologia de utilização do MVD para cálculo dos coeficientes hidrodinâmicos é baseada na comparação com resultados experimentais obtidos em laboratório para cilindros rígidos fixos e em base elástica para os quais se observa uma aderência satisfatória com relação aos resultados obtidos numericamente. O resultado do presente trabalho foi a criação de uma ferramenta alternativa no projeto de risers de grande razão de aspecto, como uma forma de avaliação do comportamento dinâmico da estrutura submetida a diferentes perfis de corrente e, consequentemente, a estimativa do tempo de sua vida útil. A característica puramente Lagrangiana do MVD possibilita a simulação de centenas de ciclos de oscilação da estrutura sendo que para os métodos Eulerianos tradicionais de CFD, como o método dos volumes finitos e o método dos elementos finitos, a simulação necessita de um tempo computacional inviável para projeto. A aderência do modelo adotado em relação aos resultados experimentais, aliada a eficiência do MVD, é um indicativo da vantagem da utilização de métodos Lagrangianos na modelagem do escoamento ao redor de cilindros flexíveis com grande razão de aspecto, em relação aos métodos Eulerianos tradicionais. A metodologia apresentada neste trabalho pode ser aplicada utilizando outros métodos de CFD, lembrando que isto pode implicar num aumento significativo do tempo computacional necessário para uma análise. / In this work the dynamic response of a high aspect ratio flexible cylinder due to vortex shedding is numerically investigated. The model is divided in two-dimensional sections along the riser length. The discrete vortex method (DVM) is employed for the assessment of the hydrodynamic forces acting on these two-dimensional sections. The hydrodynamic sections are solved independently, and the coupling among the sections is taken into account by the solution of the structure in the time domain by the finite element method implemented in Anflex code (Mourelle et al. (2001)). Parallel processing is employed to improve the performance of the method. A master-slave approach via MPI (Message Passing Interface) is used to exploit the parallelism of the present code. The riser sections are equally divided among the nodes of the cluster. Each node solves the hydrodynamic sections assigned to it. The forces acting on the sections are then passed to the master processor, which is responsible for the calculation of the displacement of the whole structure. One of the main contributions of the present work is the possibility of simulating the flow around flexible cylinders in the pos-critical regime and around bundle of risers.

Uma das mais recentes abordagens propostas na literatura para tratar o problema do alto número de Weissenberg (We) é a Formulação Log-Conformação (FLC). Nesta formulação, a equação constitutiva viscoelástica utilizada é reescrita em termos de uma variável Ψ, que é o logaritmo do tensor conformação. Apesar do potencial de aplicação da FLC, pouca atenção tem sido dirigida para análise da acurácia da solução obtida para o campo de tensões quando se utiliza esta formulação. Assim, o objetivo do presente trabalho foi estudar a acurácia da solução obtida pela FLC na análise de escoamentos de fluidos viscoelásticos usando duas geometrias padrão de estudo: placas paralelas e cavidade quadrada com tampa móvel. Primeiramente, a FLC foi implementada no pacote de CFD OpenFOAM. Em seguida foram verificados os limites do número de Weissenberg na formulação numérica padrão (Welim,P), onde para a geometria de placas paralelas foi encontrado Welim,P = 0,3 e para a geometria da cavidade quadrada com tampa móvel foi encontrado Welim,P = 0,8. Depois o código implementado foi aplicado em ambas as geometrias, comparando-se a solução obtida pela FLC com aquela da formulação padrão na faixa de We < Welim,P. Os resultados obtidos na geometria de placas paralelas apresentaram boa concordância com a solução padrão e solução analítica. Para a geometria da cavidade quadrada com tampa móvel, que não possui solução analítica, boa concordância dos resultados também foi observada em comparação com a solução padrão. Posteriormente foram comparados os resultados obtidos pela FLC na faixa de We > Welim,P. Na geometria de placas paralelas, além da boa concordância com a solução analítica, obteve-se convergência em todos os casos estudados neste trabalho, com o maior número de Weissenberg utilizado sendo igual a 8 Os resultados da geometria da cavidade quadrada com tampa móvel também apresentaram boa concordância em comparação com dados da literatura, porém a convergência foi obtida até para We = 2. Com respeito à comparação das formulações numéricas com a solução analítica, feita apenas na geometria de placas paralelas, foi observado um erro máximo de 7,57% na solução padrão e de 12,33% na FLC. Em relação à análise da qualidade das tensões usando os resíduos da equação constitutiva viscoelástica como critério de acurácia, foi verificado nas duas geometrias que os valores de tensão obtidos usando a FLC são menos acurados que aqueles obtidos pela formulação explícita no tensor das tensões nos casos em que esta última converge. Também foi observado que a acurácia diminui com o aumento do We. Esse efeito pôde ser melhor notado na geometria de placas paralelas. Uma razão para a perda de acurácia da tensão provavelmente ocorre devido à natureza matemática da transformação algébrica inversa de Ψxx para τxx. O novo solver implementado neste trabalho apresentou convergência e soluções corretas para as duas geometrias, logo foi implementado corretamente. Ele também potencializa o solver de partida viscoelastiFluidFoam ao estender simulações para uma faixa maior do número de Weissenberg. / A recent approach proposed in the literature to deal with the High Weissenberg Number Problem is the Log-Conformation formulation (LCF). In this formulation the viscoelastic constitutive equation is rewritten in terms of the logarithm of the conformation tensor Ψ. Despite the great potential application of the LCF, little attention has been given in the literature to the accuracy of the obtained stress fields. The purpose of this work was to study the solution obtained by LCF in the analysis of viscoelastic flows using two benchmark geometries: parallel plates and lid driven cavity. Firstly, the LCF was implemented in the OpenFOAM CFD package. Then, the limits of Weissenberg number for the standard numerical formulation (Welim,P) were verified, obtaining Welim,P = 0.3 for the parallel plates and Welim,P = 0.8 for the lid driven cavity. When comparing the solution obtained by the LCF with that of the standard formulation in a range of We < Welim,P, the results obtained for the parallel plates geometry showed good agreement with the standard solution and the analytical solution. For the lid driven cavity geometry, for which there is not analytical solution, good agreement with the standard solution was also observed. For We > Welim,P in the parallel plates geometry, in addition to the good agreement with the analytical solution, it was possible to obtain convergence in all the cases studied in this work, with the largest number of Weissenberg used being equal to 8 The results of the lid driven cavity geometry also presented good agreement in comparison with literature data, but convergence was obtained up to We = 2. With respect to the comparison of the numerical formulations with the analytical solution for the parallel plates geometry, a maximum error of 7.57% was observed in the standard solution and of 12.33% in the LCF. When using the residues of the viscoelastic constitutive equation as a criterion of accuracy, it was verified that for the two geometries the stress values obtained using the LCF were less accurate than those obtained by the explicit formulation in the stress tensor. It has also been observed that accuracy decreases with increasing of We. One reason for the loss of stress accuracy probably occurs because of the mathematical nature of the inverse algebraic transformation from Ψxx to τxx. The new solver implemented in this work presented convergence and correct solutions for the two geometries, so it was implemented correctly. It also potentiates the viscoelastiFluidFoam starting solver by extending simulations to a larger range of Weissenberg number.

Publicación a texto completo no autorizada por el autor / El documento digital no refiere asesor / Presenta las herramientas necesarias para evaluar hidráulicamente el nuevo sistema de bombeo ubicado en la Planta Yanacocha Norte (PYN). Dentro de los objetivos específicos se destaca los siguientes: el dimensionamiento de la capacidad del tanque (norma API 650), la configuración de la geometría del tanque, definir las características de las bombas: caudal, ADT de la instalación (TDH), tipo, potencia y NPSH disponible, seleccionar el sistema de tuberías (material, diámetro y espesor) y verificar que las velocidades del flujo en las tuberías de los sistemas se ajusten a los límites de velocidad establecidos en los criterios de diseño de tuberías de MYSRL, y calcular las presiones máximas y mínimas en las tuberías de impulsión que se ocasionan a causa del fenómeno de golpe de ariete el cual a su vez se ha producido por un corte intempestivo del suministro eléctrico, seguido a esto se utilizará para el alivio de presión una válvula anticipadora de onda. / Trabajo de suficiencia profesional

Projecto realizado na CeNTI - Centro de Nanotecnologia e Materiais Técnicos, Funcionais e Inteligentes, orient. Doutor Tiago Sotto Mayor / Tese de Mestrado Integrado. Engenharia Química. Universidade do Porto. Faculdade de Engenharia. 2010

Tese de doutoramento. Engenharia Química e Biológica. Faculdade de Engenharia. Universidade do Porto. 2011