997 resultados para Método dos elementos finitos - Programas de computador
Resumo:
Este trabalho tem por objetivo estudar a transferência de calor em tubos capilares cilíndricos utilizados na técnica de separação de moléculas denominada Eletroforese Capilar. Esta técnica é usada, por exemplo, na análise de biomoléculas e no sequenciamento de DNA, onde o controle da temperatura está diretamente ligado ao desempenho destes métodos e à qualidade dos resultados. Para empregar esta técnica, tensões elétricas da ordem de 20 kV são aplicadas entre as extremidades dos tubos capilares, que possuem normalmente 50 cm de comprimento, 350 µm de diâmetro externo e 50 µm de diâmetro interno, preenchidos por uma solução aquosa. Tais tensões geram uma corrente elétrica na solução, provocando aquecimento distribuído por Efeito Joule. Os tubos capilares são construídos em quartzo amorfo e protegidos por uma camada de material polimérico (poliimida). Para implementar o controle da temperatura, os tubos capilares são colocados em contato com um fluido de resfriamento. Num primeiro momento, os estudos são realizados por simulação numérica, empregando o Método dos Volumes Finitos em rotinas escritas em FORTRAN. São simulados casos onde os tubos são recobertos por camadas cilíndricas de materiais com uma condutividade térmica relativamente boa, com o objetivo de aumentar a superfície de troca de calor com o fluido de resfriamento. Como resultado, obtêm-se curvas da temperatura no centro dos tubos capilares em função do coeficiente de transferência de calor por convecção. Um caso de interesse é quando os tubos capilares são posicionados excentricamente ao recobrimento cilíndrico Num segundo momento, é utilizado o software de simulação numérica ANSYS CFX®, onde é simulado o resfriamento dos mesmos tubos capilares expostos a um escoamento transversal de ar a 15°C. Neste caso, também são aplicados os recobrimentos cilíndricos e, além disso, opta-se por simular o resfriamento de um arranjo de vários tubos (sistema multicapilar) dispostos entre placas de vidro, no formato de um sanduíche. Como resultados mais importantes salientam-se: a) o aumento do raio do recobrimento resulta no aumento da transferência de calor, fazendo com que a temperatura no núcleo do capilar fique estacionada em valores baixos que não comprometem as separações/análises; b) chegou-se a um valor de raio crítico da ordem de 10 mm para a condição de operação mais típicas na área da Eletroforese Capilar; c) as montagens com o tubo capilar concêntrico e excêntrico ao recobrimento não apresentam diferenças significativas no perfil de temperatura da solução tampão; e finalmente d) observa-se que o uso de duas placas de material dielétrico com os capilares posicionados em forma de sanduíche entre elas permite uma eficiente dissipação do calor gerado na solução tampão.
Análise de escoamentos incompressíveis utilizando simulação de grandes escalas e adaptação de malhas
Resumo:
No presente estudo, são apresentadas soluções numéricas de problemas de Engenharia, na área de Dinâmica dos Fluidos Computacional, envolvendo fluidos viscosos, em escoamentos incompressíveis, isotérmicos e não isotérmicos, em regime laminar e turbulento, podendo envolver transporte de massa. Os principais objetivos deste trabalho são a formulação e a aplicação de uma estratégia de adaptação automática de malhas e a inclusão de modelos de viscosidade turbulenta, integrados com um algoritmo utilizado para simular escoamentos de fluidos viscosos bi e tridimensionais, no contexto de malhas não estruturadas. O estudo é dirigido no sentido de aumentar o conhecimento a respeito das estruturas de escoamentos turbulentos e de estudar os efeitos físicos no transporte de quantidades escalares propiciando, através de técnicas de adaptação automática de malhas, a obtenção de soluções numéricas precisas a um custo computacional otimizado. As equações de conservação de massa, de balanço de quantidade de movimento e de quantidade escalar filtradas são utilizadas para simular as grandes escalas de escoamentos turbulentos e, para representar as escalas submalha, são utilizados dois modelos de viscosidade turbulenta: o modelo de Smagorinsky clássico e o modelo dinâmico. Para obter soluções numéricas com precisão, é desenvolvida e implementada uma estratégia de adaptação automática de malhas, a qual é realizada simultaneamente e interativamente com a obtenção da solução. O estudo do comportamento da solução numérica é fundamentado em indicadores de erro, com o propósito de mapear as regiões onde certos fenômenos físicos do escoamento ocorrem com maior intensidade e de aplicar nestas regiões um esquema de adaptação de malhas. A adaptação é constituída por processos de refinamento/desrefinamento e por um processo de suavização laplaciana. Os procedimentos para a implementação dos modelos de viscosidade turbulenta e a estratégia de adaptação automática de malhas são incorporados ao código computacional de elementos finitos tridimensionais, o qual utiliza elementos tetraédricos lineares. Aplicações de escoamentos de fluidos viscosos, incompressíveis, isotérmicos e não isotérmicos em regime laminar e turbulento são simuladas e os resultados são apresentados e comparados com os obtidos numérica ou experimentalmente por outros autores.
Resumo:
A presente dissertação está dividida em 6 (seis) capítulos, um dos quais é esta própria introdução, com a qual se pretendeu, além de mostrar a importância do problema que será analisado, também realçar o seu papel no panorama de Processamento de Dados brasileiro de hoje. O capítulo seguinte apresenta o trabalho alvo deste estudo, ou seja, programação de computadores, situando-o em relação ao todo da função de processamento de dados na em presa e caracterizando-se também os aspectos organizacionais dessa função além de realçar-se os aspectos básicos do ciclo de vida dos sistemas automatizados. Em seguida, no terceiro capítulo, são enumerados e analisados os vários fatores que atuam sobre a tarefa de programação de computadores os quais encontram-se separados em três categorias; Fatores ambientais, Variáveis Inerentes a Cada Programa e Experiência do Programador. No quarto capítulo são consideradas vá- rias técnicas de Medida do Trabalho, examinando-se cada uma sob a ótica de uma possível aplicação com programadores de computador. A metodologia proposta para se estabelecer padrões de tempo e produtividade para esses profissionais está apresentada no quinto capítulo. Resumidamente ela consiste em se escolher e observar um grupo de programas de computador, quantificando simultaneamente suas variáveis peculiares conforme conceituadas no terceiro capítulo e o esforço dispendido em cada um. Através do uso da análise de regressão múltipla são então obtidas as equações que darão os tempos normais necessários para se produzir cada tipo de programa em função das variáveis relevantes em cada tipo. Finalmente no último capítulo são mostrados os principais benefícios obtidos com a aplicação da metodologia sugerindo-se possibilidades para se expandir o estudo
Resumo:
A introdução de um novo material para a execução de estruturas, implica um amplo estudo sobre o mesmo, para que seja aceite como um material de utilização, ou solução, corrente na engenharia civil. Os varões plásticos reforçados com fibra de vidro, GFRP,apresentam-se com diversas vantagens para uma solução de substituição das tradicionais armaduras de aço no betão armado. Um Estado da Arte foi apurado para determinar o ponto de situação do conhecimento e dos estudos realizados por vários autores, desde o início das pesquisas até aos mais recentes estudos disponíveis. Neste seguimento, o presente estudo incide na análise do comportamento à flexão de vigas contínuas de betão armado com varões de GFRP, submetidas a um carregamento distribuído. O trabalho realiza-se com base em 18 modelos de elementos finitos para uma análise paramétrica, sendo os parâmetros estudados, o vão, a percentagem de armadura de flexão e a esbelteza da secção. Para validação do modelo foram, modeladas 7 vigas reais ensaiadas no âmbito de trabalhos de investigação de dois autores Os resultados das modelações das análises paramétricas foram comparados com uma análise elástica e uma análise plástica, utilizando a metodologia do ACI 440.1R-06, para a carga última de dimensionamento. Outras comparações foram feitas através dos limites de deformação do EN 1992-1-1, 2004. Analisa-se a distribuição de esforços, não elásticos, devido à perda de rigidez da viga. Os modelos das vigas mais esbeltas e com menor percentagem de armadura mostraram maior concordância com a análise realizada. As deformações comparadas com o limite L/250 mostram concordância de resultados.
Resumo:
O dimensionamento de uma estrutura de betão armado depende em grande parte do tipo de metodologia de cálculo utilizado, sendo a análise não linear a que melhor reflete o comportamento da estrutura. Neste tipo de análise é preciso ter em conta alguns aspetos relativos à ductilidade como, a abertura de fendas, a formação de rótulas plásticas e as relações não lineares do material. Entende-se por ductilidade a capacidade de deformação da estrutura sem perda significativa de resistência após atingida a sua capacidade elástica. Esta pode ser quantificada através de diversos indicadores. No caso das vigas é corrente o uso do índice de ductilidade de flecha, que pode ser útil no pré-dimensionamento estrutural pois pode permitir aos projetistas pré-dimensionar vigas segundo um critério de ductilidade. Assim realizou-se uma pesquisa de maneira a identificar quais os fatores mais importantes na ductilidade de uma viga e através do programa cálculo de elementos finitos DIANA 9.4 simularam-se 162 vigas variando a percentagem de armadura mecânica total, o comprimento do vão, a altura útil da secção, a redistribuição de momentos e a classe de betão. Observou-se que o aumento da percentagem de armadura mecânica total origina uma diminuição no índice de ductilidade de cerca de 25% a cada incremento de 0.1 de percentagem de armadura mecânica total, a esbelteza provoca um aumento de 5% a cada 5 unidades no betão C30 e de 2 a 3% no betão C70 e a redistribuição de momentos provocou um aumento na ductilidade. Os resultados mostram que os fatores que mais influenciam na ductilidade das vigas são a percentagem mecânica de armadura total e a classe de resistência do betão, no entanto existem outras que deverão sofrer maior aprofundamento em trabalhos futuros como a utilização de armadura à compressão e os níveis de confinamento do betão.
Resumo:
O presente trabalho pretende contribuir para o aprofundamento do conhecimento geológico e geotécnico das formações vulcânicas da Ilha da Madeira, com especial destaque para os parâmetros geomecânicos obtidos através de retroanálise, tendo por base as leituras em obra do Túnel do Cortado (Faial, Santana; Madeira). Assim, o trabalho incide inicialmente sobre os túneis rodoviários na Ilha, com principal atenção nas características das suas secções transversais e emboquilhamentos. Segue-se um levantamento de todos os túneis construídos entre os anos de 1990 e 2013. Uma vez terminada esta apresentação, é feita uma caracterização geológica-geotécnica dos maciços da região, com especial atenção para as zonas geotécnicas afetas ao túnel em questão. São ainda abordadas metodologias de classificações geomecânicas, de retroanálise e modelação computacional. Segue-se um processo de retroanálise através do qual foram gerados vários modelos por forma a efetuar uma análise bidimensional das secções do túnel, onde se procurou estabelecer uma relação entre as leituras de obra e os parâmetros de resistência e deformabilidade do maciço em questão. Por fim, os resultados referentes aos parâmetros de resistência obtidos através dos modelos são analisados, procedendo-se à sua comparação com outros estudos existentes para a Ilha.
Resumo:
Oil wells subjected to cyclic steam injection present important challenges for the development of well cementing systems, mainly due to tensile stresses caused by thermal gradients during its useful life. Cement sheath failures in wells using conventional high compressive strength systems lead to the use of cement systems that are more flexible and/or ductile, with emphasis on Portland cement systems with latex addition. Recent research efforts have presented geopolymeric systems as alternatives. These cementing systems are based on alkaline activation of amorphous aluminosilicates such as metakaolin or fly ash and display advantageous properties such as high compressive strength, fast setting and thermal stability. Basic geopolymeric formulations can be found in the literature, which meet basic oil industry specifications such as rheology, compressive strength and thickening time. In this work, new geopolymeric formulations were developed, based on metakaolin, potassium silicate, potassium hydroxide, silica fume and mineral fiber, using the state of the art in chemical composition, mixture modeling and additivation to optimize the most relevant properties for oil well cementing. Starting from molar ratios considered ideal in the literature (SiO2/Al2O3 = 3.8 e K2O/Al2O3 = 1.0), a study of dry mixtures was performed,based on the compressive packing model, resulting in an optimal volume of 6% for the added solid material. This material (silica fume and mineral fiber) works both as an additional silica source (in the case of silica fume) and as mechanical reinforcement, especially in the case of mineral fiber, which incremented the tensile strength. The first triaxial mechanical study of this class of materials was performed. For comparison, a mechanical study of conventional latex-based cementing systems was also carried out. Regardless of differences in the failure mode (brittle for geopolymers, ductile for latex-based systems), the superior uniaxial compressive strength (37 MPa for the geopolymeric slurry P5 versus 18 MPa for the conventional slurry P2), similar triaxial behavior (friction angle 21° for P5 and P2) and lower stifness (in the elastic region 5.1 GPa for P5 versus 6.8 GPa for P2) of the geopolymeric systems allowed them to withstand a similar amount of mechanical energy (155 kJ/m3 for P5 versus 208 kJ/m3 for P2), noting that geopolymers work in the elastic regime, without the microcracking present in the case of latex-based systems. Therefore, the geopolymers studied on this work must be designed for application in the elastic region to avoid brittle failure. Finally, the tensile strength of geopolymers is originally poor (1.3 MPa for the geopolymeric slurry P3) due to its brittle structure. However, after additivation with mineral fiber, the tensile strength became equivalent to that of latex-based systems (2.3 MPa for P5 and 2.1 MPa for P2). The technical viability of conventional and proposed formulations was evaluated for the whole well life, including stresses due to cyclic steam injection. This analysis was performed using finite element-based simulation software. It was verified that conventional slurries are viable up to 204ºF (400ºC) and geopolymeric slurries are viable above 500ºF (260ºC)
Resumo:
The oil production in mature areas can be improved by advanced recovery techniques. In special, steam injection reduces the viscosity of heavy oils, thus improving its flow to surrounding wells. On the other hand, the usually high temperatures and pressures involved in the process may lead to cement cracking, negatively affecting both the mechanical stability and zonal isolation provided by the cement sheath of the well. The addition of plastic materials to the cement is an alternative to prevent this scenario. Composite slurries consisting of Portland cement and a natural biopolymer were studied. Samples containing different contents of biopolymer dispersed in a Portland cement matrix were prepared and evaluated by mechanical and rheological tests in order to assess their behavior according to API (American Petroleum Institute) guidelines. FEM was also applied to map the stress distribution encountered by the cement at bottom bole. The slurries were prepared according to a factorial experiment plan by varying three parameters, i.e., cement age, contents of biopolymer and water-to-cement ratio. The results revealed that the addition of the biopolymer reduced the volume of free water and the setting time of the slurry. In addition, tensile strength, compressive strength and toughness improved by 30% comparing hardened composites to plain Portland slurries. FEM results suggested that the stresses developed at bottomhole may be 10 to 100 times higher than the strength of the cement as evaluated in the lab by unconfined mechanical testing. An alternative approach is proposed to adapt the testing methodology used to evaluate the mechanical behavior of oilwell cement slurries by simulating the confined conditions encountered at bottornhole
Resumo:
The research and development of wind turbine blades are essential to keep pace with worldwide growth in the renewable energy sector. Although currently blades are typically produced using glass fiber reinforced composite materials, the tendency for larger size blades, particularly for offshore applications, has increased the interest on carbon fiber reinforced composites because of the potential for increased stiffness and weight reduction. In this study a model of blade designed for large generators (5 MW) was studied on a small scale. A numerical simulation was performed to determine the aerodynamic loading using a Computational Fluid Dynamics (CFD) software. Two blades were then designed and manufactured using epoxy matrix composites: one reinforced with glass fibers and the other with carbon fibers. For the structural calculations, maximum stress failure criterion was adopted. The blades were manufactured by Vacuum Assisted Resin Transfer Molding (VARTM), typical for this type of component. A weight comparison of the two blades was performed and the weight of the carbon fiber blade was approximately 45% of the weight of the fiberglass reinforced blade. Static bending tests were carried out on the blades for various percentages of the design load and deflections measurements were compared with the values obtained from finite element simulations. A good agreement was observed between the measured and calculated deflections. In summary, the results of this study confirm that the low density combined with high mechanical properties of carbon fibers are particularly attractive for the production of large size wind turbine blades
Resumo:
The aim of this work is the numerical simulation of the mechanical performance of concrete affected by Alkali-Aggregate Reaction or RAA, reported by Stanton in 1940. The RAA has aroused attention in the context of Civil Engineering from the early 80, when they were reported consequences of his swelling effect in concrete structures, including cracking, failure and loss of serviceability. Despite the availability of experimental results the problem formulation still lacks refinement so that your solution remains doubtful. The numerical simulation is important resource for the assessment of damages in structures caused by the reaction, and their recoveries The tasks of support of this work were performed by means of the finite element approach, about orthotropic non-linear formulation, and, thermodynamic model of deformation by RAA. The results obtained revealed that the swelling effect of RAA induced decline of the mechanical performance of concrete by decreasing the margin of safety prior to the material failure. They showed that the temperature influences, exclusively, the kinetics of the reaction, so that the failure was the more precocious the higher the temperature of the solid mass of concrete
Resumo:
The processing of materials through plasma has been growing enough in the last times in several technological applications, more specifically in surfaces treatment. That growth is due, mainly, to the great applicability of plasmas as energy source, where it assumes behavior thermal, chemical and/or physical. On the other hand, the multiplicity of simultaneous physical effects (thermal, chemical and physical interactions) present in plasmas increases the complexity for understanding their interaction with solids. In that sense, as an initial step for the development of that subject, the present work treats of the computational simulation of the heating and cooling processes of steel and copper samples immersed in a plasma atmosphere, by considering two experimental geometric configurations: hollow and plane cathode. In order to reach such goal, three computational models were developed in Fortran 90 language: an one-dimensional transient model (1D, t), a two-dimensional transient model (2D, t) and a two-dimensional transient model (2D, t) which take into account the presence of a sample holder in the experimental assembly. The models were developed based on the finite volume method and, for the two-dimensional configurations, the effect of hollow cathode on the sample was considered as a lateral external heat source. The main results obtained with the three computational models, as temperature distribution and thermal gradients in the samples and in the holder, were compared with those developed by the Laboratory of Plasma, LabPlasma/UFRN, and with experiments available in the literature. The behavior showed indicates the validity of the developed codes and illustrate the need of the use of such computational tool in that process type, due to the great easiness of obtaining thermal information of interest
Resumo:
Multiphase flows in ducts can adopt several morphologies depending on the mass fluxes and the fluids properties. Annular flow is one of the most frequently encountered flow patterns in industrial applications. For gas liquid systems, it consists of a liquid film flowing adjacent to the wall and a gas core flowing in the center of the duct. This work presents a numerical study of this flow pattern in gas liquid systems in vertical ducts. For this, a solution algorithm was developed and implemented in FORTRAN 90 to numerically solve the governing transport equations. The mass and momentum conservation equations are solved simultaneously from the wall to the center of the duct, using the Finite Volumes Technique. Momentum conservation in the gas liquid interface is enforced using an equivalent effective viscosity, which also allows for the solution of both velocity fields in a single system of equations. In this way, the velocity distributions across the gas core and the liquid film are obtained iteratively, together with the global pressure gradient and the liquid film thickness. Convergence criteria are based upon satisfaction of mass balance within the liquid film and the gas core. For system closure, two different approaches are presented for the calculation of the radial turbulent viscosity distribution within the liquid film and the gas core. The first one combines a k- Ɛ one-equation model and a low Reynolds k-Ɛ model. The second one uses a low Reynolds k- Ɛ model to compute the eddy viscosity profile from the center of the duct right to the wall. Appropriate interfacial values for k e Ɛ are proposed, based on concepts and ideas previously used, with success, in stratified gas liquid flow. The proposed approaches are compared with an algebraic model found in the literature, specifically devised for annular gas liquid flow, using available experimental results. This also serves as a validation of the solution algorithm
Resumo:
Annular flow is the prevailing pattern in transport and energy conversion systems and therefore, one of the most important patterns in multiphase flow in ducts. The correct prediction of the pressure gradient and heat transfer coefficient is essential for optimizing the system s capacity. The objective of this work is to develop and implement a numerical algorithm capable of predicting hydrodynamic and thermal characteristics for upflow, vertical, annular flow. The numerical algorithm is then complemented with the physical modeling of phenomena that occurs in this flow pattern. These are, turbulence, entrainment and deposition and phase change. For the development of the numerical model, axial diffusion of heat and momentum is neglected. In this way the time-averaged equations are solved in their parabolic form obtaining the velocity and temperature profiles for each axial step at a time, together with the global parameters, namely, pressure gradient, mean film thickness and heat transfer coefficient, as well as their variation in the axial direction. The model is validated for the following conditions: fully-developed laminar flow with no entrainment; fully developed laminar flow with heat transfer, fully-developed turbulent flow with entrained drops, developing turbulent annular flow with entrained drops, and turbulent flow with heat transfer and phase change
Resumo:
A critical problem in mature gas wells is the liquid loading. As the reservoir pressure decreases, gas superficial velocities decreases and the drag exerted on the liquid phase may become insufficient to bring all the liquid to the surface. Liquid starts to drain downward, flooding the well and increasing the backpressure which decreases the gas superficial velocity and so on. A popular method to remedy this problem is the Plunger Lift. This method consists of dropping the "plunger"to the bottom of the tubing well with the main production valve closed. When the plunger reaches the well bottom the production valve is opened and the plunger carry the liquid to the surface. However, models presented in literature for predicting the behavior in plunger lift are simplistic, in many cases static (not considering the transient effects). Therefore work presents the development and validation of a numerical algorithm to solve one-dimensional compressible in gas wells using the Finite Volume Method and PRIME techniques for treating coupling of pressure and velocity fields. The code will be then used to develop a dynamic model for the plunger lift which includes the transient compressible flow within the well
Resumo:
Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)