É notório que a degradação ambiental vem ao longo do tempo posicionando-se como um dos principais problemas do mundo moderno. Dentre as várias questões do interesse ambiental podemos destacar a ascensão da bolha de metano, em reservatórios hidrelétricos, desde o sedimento anóxico no fundo do reservatório até a interface água atmosfera. Neste contexto, a presente tese vêm propor uma nova modelagem matemática para a ascensão da bolha axissimétrica em fluidos newtonianos/nãonewtonianos e mostrar resultados numéricos simulados. Desta forma, o estado da arte estaria elevado a posição de permitir, via Matemática e Simulação Numérica-Computacional, a análise do transporte de metano em reservatórios de hidrelétricas através da bolha / It is well-known that environmental degradation has come along the time positioning as one of the main problems from modern world. Among several questions of the environmental interest may emphasize the methane bubble rise in hydroelectric reservoirs from the anoxic sediment in the bottom of reservoir until water interface atmosphere. In this context, the current thesis has come to propose a new mathematic modeling to the rise of the axisymmetric bubble in the newtonian/non-newtonian fluids and display numerical results simulated. Therefore, the state of art would be ascended to a position to permit, via Mathematics and Computing-Numeric Simulation, the analysis of transport of methane in hydroelectric reservoirs through of bubble

No contexto do método MAC e baseado em esquemas de diferenças finitas, este trabalho apresenta três estudos: i) uma análise de estabilidade, ii) o desenvolvimento de técnicas implícitas e, iii) a construção de métodos de projeção para escoamentos com superfície livre. Na análise de estabilidade, o principal resultado mostra que o método de Crank-Nicolson torna-se condicionalmente estável quando aplicado para uma malha deslocada com a discretiza ção explícita das condições de contorno do tipo Dirichlet. Entretanto, o mesmo método com condições de contorno implícitas é incondicionalmente estável. Para obter métodos mais estáveis, formulações implícitas são desenvolvidas para a equação da pressão na superfície livre, derivada da condição de tensão normal. Esta estratégia resulta no acoplamento dos campos de velocidade e pressão, o que exige a introdução de novos métodos de projeção. Os métodos de projeção assim desenvolvidos resultam em novas metodologias para escoamentos com superfície livre que são apropriados para o tratamento de problemas com baixo número de Reynolds. Além disso, mostra-se que os métodos propostos podem ser aplicados para fluidos viscoelásticos. Novas estratégias são derivadas para obter métodos de projeção de segunda ordem de precisão para escoamentos com superfícies livres. Além dos resultados teóricos sobre a estabilidade de esquemas numéricos, técnicas implícitas e métodos de projeção, testes computacionais são realizados e comparados para consolidação da teoria apresentada. Os resultados numéricos são obtidos no sistema FREEFLOW. A eficiência e robustez das técnicas desenvolvidas neste trabalho são demonstradas na solução de problemas tridimensionais complexos com superfície livre e baixo número de Reynolds, incluindo os problemas do jato oscilante e do inchamento do extrudado / In the context of the MAC method and based on finite difference schemes, this work presents three studies: i) a stability analysis, ii) the development of implicit techniques, and iii) the construction of projection methods for free surface flows. In the stability analysis, the main result shows a precise stability restriction on the Crank-Nicolson method when one uses a staggered grid with Dirichlet explicit boundary conditions. However, the same method with implicit boundary conditions becomes unconditionally stable. In order to obtain more stable methods, implicit formulations are applied for the pressure equation at the free surface, which is derived from the normal stress condition. This approach results in a coupling of the velocity and pressure fields; hence new projection methods for free surface flows need to be developed. The developed projection methods result in new methodologies for low Reynolds number free surface flows. It is also shown that the proposed methods can be applied for viscoelastic fluids. New strategies are derived for obtaining second-order accurate projection methods for free surface flows. In addition to the theoretical results on the stability of numerical schemes, implicit techniques and projection methods, computational tests are carried out and the results compared to consolidate the theory. The numerical results are obtained by the FREEFLOW system. The eficiency and robustness of the techniques in this work are demonstrated by solving complex tridimensional problems involving free surface and low Reynolds numbers, including the jet buckling and the extrudate swell problems

Escoamentos viscoelásticos não estacionários com superfícies livres são comuns em muitos processos industriais e diversas técnicas numéricas têm sido empregadas para reproduzir computacionalmente estes processos. A maioria dos modelos empregados utiliza equações diferenciais na definição do tensor de tensões. Porém, para alguns grupos de fluidos complexos, por exemplo, fluidos de Boger, os modelos integrais mostram-se mais capacitados em fornecer uma boa aproximação para os comportamentos não lineares desses fluidos. Este trabalho trata da solução numérica do modelo constitutivo integral KBKZ-PSM para escoamentos transientes bidimensionais com superfícies livres. O método numérico proposto é uma técnica numérica que utiliza diferenças finitas para simular escoamentos com superfícies livres na presença de paredes sólidas. As principais características do método numérico proposto são: solução das equações de conservação de quantidade de movimento e massa utilizando um método semi-implícito; a condição de contorno na superfície livre é acoplada à equação de Poisson, o que garante conservação de massa; a discretização do tempo t é realizada por uma nova técnica numérica; o tensor de Finger é calculado pelo método dos campos de deformação e avançado no tempo pelo método de Euler modificado. Essa nova técnica é verificada em escoamentos cisalhantes e elongacionais. Adicionalmente, uma solução analítica desenvolvida para escoamentos em canais bidimensionais é empregada para verificar e analisar a convergência do método proposto. Com relação a escoamentos com superfícies livres, a convergência é verificada por meio de refinamento de malha nas simulações de um jato incidindo sobre placa rígida e no problema do inchamento do extrudado. Finalmente, o método é aplicado para investigar os problemas jet buckling e inchamento do extrudado de fluidos KBKZ-PSM. / Unsteady viscoelastic free surface flows are common in many industrial processes and a variety of numerical techniques have been employed to simulate these flows. The majority of constitutive models employed are based on differential equations to define the extra stress tensor. However, for some complex fluids, for instance, Boger fluids, integral models are more adequate to approximate the nonlinear behaviour of these fluids. This work deals with the numerical solution of the integral constitutive model KBKZ-PSM for two-dimensional unsteady free surface flows. The proposed numerical method is a numerical technique that employs finite differences to simulate moving free surface flows that interact with solid walls. The main features of the method are: numerical solution of the momentum and mass equations by an implicit method; the pressure condition on the free surface is implicitly coupled with the Poisson equation for obtaining the pressure field from mass conservation; a novel scheme for defining the past times t is employed; the Finger tensor is calculated by the deformation fields method and is advanced in time by the modified Euler method. This new technique is verified by solving shear and uniaxial elongational flows. Moreover, an analytic solution for channel flows is obtained that is used in the verification and convergence analysis of the proposed methodology. For free surface flows, the assessment of convergence lies on the mesh refinement on the simulation of a jet impinging on a flat surface and the extrudade swell problem. Finally, the new method is applied to investigate the jet buckling phenomenon and extrudate swell of KBKZ-PSM fluids.

[ES] La forma más exacta de conocer el desplazamiento de los neutrones a través de un medio material se consigue resolviendo la Ecuación del Transporte Neutrónico. Tres diferentes aproximaciones de esta ecuación se han investigado en esta tesis: Ecuación del transporte neutrónico resuelta por el método de Ordenadas Discretas, Ecuación de la Difusión y Ecuación de Armónicos Esféricos Simplificados. Para resolver estás ecuaciones se estudian diferentes esquemas del Método de Diferencias Finitas. La solución a estas ecuaciones describe la población de neutrones y las reacciones ocasionadas dentro de un reactor nuclear. A su vez, estas variables están relacionadas con el flujo y la potencia, parámetros fundamentales para el Análisis de Seguridad Nuclear. La tesis introduce la definición de las ecuaciones mencionadas y en particular se detallan para el estado estacionario. Se plantea el Método Modal como solución a los problemas de autovalores definidos por dichas ecuaciones. Primero se desarrollan varios algoritmos para la resolución del estado estacionario de la Ecuación del Transporte de Neutrones con el Método de Ordenadas Discretas para la discretización angular y el Método de Diferencias Finitas para la discretización espacial. Se ha implementado una formulación capaz de resolver el problema de autovalores para cualquier número de grupos energéticos con upscattering y anisotropía. Varias cuadraturas utilizadas por este método en su resolución angular han sido estudiadas e implementadas para cualquier orden de aproximación de Ordenadas Discretas. Además, otra formulación se desarrolla para la solución del problema fuente de la ecuación del transporte neutrónico. A continuación, se lleva a cabo un algoritmo que permite resolver la Ecuación de la Difusión de Neutrones con dos variantes del método de diferencias Finitas, una centrada en celda y otra en vértice o nodo. Se utiliza también el Método Modal calculando cualquier número de autovalores para varios grupos de energía y con upscattering. También se implementan los dos esquemas del Método de Diferencias Finitas anteriormente mencionados en el desarrollo de diferentes algoritmos para resolver las Ecuaciones de Armónicos Esféricos Simplificados. Además, se ha realizado un análisis de diferentes aproximaciones de las condiciones de contorno. Finalmente, se han realizado cálculos de la constante de multiplicación, los modos subcríticos, el flujo neutrónico y la potencia para diferentes tipos de reactores nucleares. Estas variables resultan esenciales en Análisis de Seguridad Nuclear. Además, se han realizado diferentes estudios de sensibilidad de parámetros como tamaño de malla, orden utilizado en cuadraturas o tipo de cuadraturas. / [CA] La forma més exacta de conèixer el desplaçament dels neutrons a través d'un mitjà material s'aconsegueix resolent l'Equació del Transport Neutrònic. Tres diferents aproximacions d'esta equació s'han investigat en aquesta tesi: Equació del Transport Neutrònic resolta pel mètode d'Ordenades Discretes, Equació de la Difusió i Equació d'Ármonics Esfèrics Simplificats. Per a resoldre estes equacions s'estudien diferents esquemes del Mètode de Diferències Finites. La solució a estes equacions descriu la població de neutrons i les reaccions ocasionades dins d'un reactor nuclear. Al seu torn, estes variables estan relacionades amb el flux i la potència, paràmetres fonamentals per a l'Anàlisi de Seguretat Nuclear. La tesi introduïx la definició de les equacions mencionades i en particular es detallen per a l'estat estacionari. Es planteja el Mètode Modal com a solució als problemes d'autovalors definits per les dites equacions. Primer es desenvolupen diversos algoritmes per a la resolució de l'estat estacionari de l'Equació del Transport de Neutrons amb el Mètode d'Ordenades Discretes per a la discretiztació angular i el Mètode de Diferències Finites per a la discretització espacial. S'ha implementat una formulació capaç de resoldre el problema d'autovalors per a qualsevol nombre de grups energètics amb upscattering i anisotropia. Diverses quadratures utilitzades per este mètode en la seua resolució angular han sigut estudiades i implementades per a qualsevol orde d'aproximació d'Ordenades Discretes. A més, una altra formulació es desenvolupa per a la solució del problema font de l'Equació del Transport Neutrònic. A continuació, es du a terme un algoritme que permet resoldre l'Equació de la Difusió de Neutrons amb dos variants del mètode de Diferències Finites, una centrada en cel·la i una altra en vèrtex o node. S'utilitza també el Mètode Modal calculant qualsevol nombre d'autovalors per a diversos grups d'energia i amb upscattering. També s'implementen els dos esquemes del Mètode de Diferències Finites anteriorment mencionats en el desenvolupament de diferents algoritmes per a resoldre les Equacions d'Harmònics Esfèrics Simplificats. A més, s'ha realitzat una anàlisi de diferents aproximacions de les condicions de contorn. Finalment, s'han realitzat càlculs de la constant de multiplicació, els modes subcrítics, el flux neutrònic i la potència per a diferents tipus de reactors nuclears. Estes variables resulten essencials en Anàlisi de Seguretat Nuclear. A més, s'han realitzat diferents estudis de sensibilitat de paràmetres com la grandària de malla, orde utilitzat en quadratures o tipus de quadratures. / [EN] The most accurate way to know the movement of the neutrons through matter is achieved by solving the Neutron Transport Equation. Three different approaches to solve this equation have been investigated in this thesis: Discrete Ordinates Neutron Transport Equation, Neutron Diffusion Equation and Simplified Spherical Harmonics Equations. In order to solve the equations, different schemes of the Finite Differences Method were studied. The solution of these equations describes the population of neutrons and the occurred reactions inside a nuclear system. These variables are related with the flux and power, fundamental parameters for the Nuclear Safety Analysis. The thesis introduces the definition of the mentioned equations. In particular, they are detailed for the steady state case. The Modal Method is proposed as a solution to the eigenvalue problems determined by the equations. First, several algorithms for the solution of the steady state of the Neutron Transport Equation with the Discrete Ordinates Method for the angular discretization and Finite Difference Method for spatial discretization are developed. A formulation able to solve eigenvalue problems for any number of energy groups, with scattering and anisotropy has been developed. Several quadratures used by this method for the angular discretization have been studied and implemented for any order of approach of the discrete ordinates. Furthermore, an adapted formulation has been developed as a solution of the source problem for the Neutron Transport Equation. Next, an algorithm is carried out that allows to solve the Neutron Diffusion Equation with two variants of the Finite Difference Method, one with cell centered scheme and another edge entered. The Modal method is also used for calculating any number of eigenvalues for several energy groups and upscattering. Both Finite Difference schemes mentioned before are also implemented to solve the Simplified Spherical Harmonics Equations. Moreover, an analysis of different approaches of the boundary conditions is performed. Finally, calculations of the multiplication factor, subcritical modes, neutron flux and the power for different nuclear reactors were carried out. These variables result essential in Nuclear Safety Analysis. In addition, several sensitivity studies of parameters like mesh size, quadrature order or quadrature type were performed. / Me gustaría dar las gracias al Ministerio de Economía, Industria y Competitividad y a la Agencia Estatal de Investigación de España por la concesión de mi contrato predoctoral de formación de personal investigador con referencia BES-2016-076782. La ayuda económica proporcionada por este contrato fue esencial para el desarrollo de esta tesis, así como para el financiamiento de una estancia. / Morato Rafet, S. (2020). Contributions to solve the Multi-group Neutron Transport equation with different Angular Approaches [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/159271

[pt] Para avaliar projetos de investimento florestal, a Teoria de Opções Reais é utilizada com o propósito de incorporar questões relacionadas a incertezas e flexibilidades gerenciais. O objetivo desta dissertação é desenvolver uma modelagem para valoração de um projeto florestal, no qual a variação do estoque de árvores se aproxima do crescimento real de uma floresta, utilizando uma adaptação da equação logística. Desta forma, para um determinado nível de saturação, o estoque se estabiliza. Adicionalmente, o trabalho busca quantificar os benefícios econômicos de uma política ótima de produção a partir do corte de árvores. Para a análise proposta, os resultados são obtidos através do método de diferenças finitas explícitas. O problema apresenta três variáveis independentes - estoque, tempo e preço, sendo este modelado como um Movimento Geométrico Browniano - e duas variáveis dependentes - a taxa de corte e o valor da opção de investimento. É apresentada uma aplicação para uma floresta de eucalipto e os resultados são comparados considerando outras alternativas para a evolução do estoque de árvores, bem como para a decisão da quantidade a ser cortada. Para esta comparação, avalia-se o caso em que o estoque segue um modelo estocástico e o caso em que a taxa de corte é fixa. Os resultados mostram que é vantajoso adotar uma política ótima de corte, corroborando resultados obtidos em trabalhos anteriores de opções reais na área florestal. Além disso, o maior valor da opção é obtido quando o estoque é modelado pela equação logística para o crescimento. / [en] In order to incorporate uncertainty and managerial flexibilities in forestry investment projects, financial literature uses the real options approach. The aim of this work is to develop a model for valuation of a forestry project, in which the inventory growth of trees follows a logistic equation based in the estimated real growth of a forest. Thus, for a given level of saturation, the inventory stabilizes. In addition, this dissertation aims to quantify the economic benefits of an optimal production policy from cutting trees. For the proposal analysis, the results are obtained by the finite difference method in the explicit form. The problem has three independent variables – inventory, time and price, which is modeled as a Geometric Brownian Motion – and two dependent variables – the cutting rate and the value of the investment option. An application for an eucalyptus forest is presented and the results are compared considering other alternatives for the evolution of the inventory and for the decision on the cutting amount. For this purpose, two different assumptions are made, considering a stochastic model for the inventory and a fixed shear rate. The results show that it is advantageous to adopt a optimal cutting policy, confirming results obtained in previous studies of real options in forestry. Moreover, the greatest option value is obtained when the inventory is modeled by logistic equation for growth.

O objetivo desse trabalho é desenvolver métodos numéricos para simular escoamentos de cristais líquidos nemáticos governados pelas equações dinâmicas de Ericksen-Leslie. São apresentados dois métodos numéricos para a simulação de escoamentos de cristais líquidos nemáticos. O primeiro método foi desenvolvido para simular escoamentos tridimensionais de cristais líquidos nemáticos sob efeito de forte campo magnético enquanto que o segundo método foi desenvolvido para a simulação de escoamentos bidimensionais. Utilizando a notação de Einstein, as equações dinâmicas de Ericksen-Leslie são apresentadas. Empregando variáveis primitivas e coordenadas cartesianas, as equações governantes para escoamentos de cristais líquidos nemáticos são derivadas e as formulações matemáticas para a obtenção dos métodos numéricos são apresentadas. As equações descrevendo os métodos numéricos são resolvidas por um método numérico baseado na metodologia GENSMAC3D para o caso tridimensional enquanto que o método bidimensional é baseado na metodologia GENSMAC (GENeralized-Simplified-Marker-And-Cell). Em ambos os métodos, a técnica de diferenças finitas em uma malha deslocada é utilizada. As equações que descrevem as técnicas numéricas desenvolvidas foram incorporadas aos ambientes de simulação Freeflow2D e Freeflow3D. As condições de contorno para cada tipo de contorno são descritas em detalhes. A solução analítica apresentada por Stewart para o escoamento entre duas placas paralelas é utilizada para a validação do método numérico tridimensional. Empregando as hipóteses de escoamento desenvolvido e que o ângulo de orientação do diretor é pequeno, uma solução analítica para o escoamento em um canal bidimensional é encontrada. O método numérico bidimensional é então validado utilizando a solução analítica obtida. Utilizando refinamento de malha, resultados de convergência dos métodos numéricos são apresentados. Os métodos numéricos desenvolvidos nesse trabalho são aplicados para a simulação dos seguintes problemas: escoamento de um cristal líquido nemático em um canal tridimensional; investigação numérica do escoamento em L-canais e escoamento através de uma contração 4:1 e de uma expansão planar 4:1 / The aim of this work is to develop numerical methods capable of simulating nematic liquid crystal flows described by the dynamic Ericksen-Leslie equations. Two numerical techniques have been proposed: a numerical method for solving three-dimensional flows of nematic liquid crystals under strong magnetic fields and a method for simulating flows of nematic liquid crystals in two-dimensional complex geometries. Both methods employ Cartesian coordinates using primitive variables of pressure and velocity. These techniques are based on the GENSMAC (GENeralized-Simplified-Marker-And-Cell) methodology and a detailed description of the equations involved is presented. The resulting governing equations are solved by the finite difference method on a staggered grid. The three-dimensional technique was applied to solve fully developed flow between two paralel plates for which an analytic solution exists. By using this analytic solution, validation and convergence results of the developed numerical technique were obtained. To validate the two-dimensional method developed herein, an analytic solution for steady state flow in a 2D-channel was found which was used to obtain validation and convergence results. The method was then applied to simulate the flow in L-shaped channels, flow through a 4:1 contraction and flow in a 4:1 expansion

Uma grande limitação dos métodos de diferenças finitas é que eles estão restritos a malhas e domínios retangulares. Para descrever escoamentos em domínios complexos, como, por exemplo, problemas com superfícies livres, faz-se necessário o uso de técnicas acessórias. O método de interfaces imersas é uma dessas técnicas. Nesse trabalho, primeiramente foi desenvolvido um método de projeção, totalmente livre de pressão, para as equações de Navier-Stokes com variáveis primitivas em malha deslocada. Esse método é baseado em diferenças finitas compactas, possuindo segunda ordem temporal e quarta ordem espacial. Esse método foi combinado com o método de interface imersa de Linnick e Fasel [2] para resolver numericamente as equações de Stokes com quarta ordem de precisão. A verificação do código foi feita por meio do método das soluções manufaturadas e da comparação com resultados de outros autores em problemas clássicos da literatura. / A great limitation of finite differences methods is that they are restricted to retangular meshes and domains. In order to describe flows in complex domains, e.g. free surface problems, it is necessary to use accessory techniques. The immersed interface method is one of such techniques. In the present work, firstly, a projection method was developed, which is completely pressure-free, for the Navier-Stokes equations with primitive variables in a staggered mesh. This method is based on compact finite differences, with temporal second-order precision and spatial foruth-order precision. This method was combined with the immersed interface method from Linnick e Fasel [2] in order to numerically solve the Stokes equations with fourth-order precision. The verification of the code was performed with the manufactured solutions method and by comparing results with other authors for some classical problems in the literature.

Made available in DSpace on 2016-08-10T10:40:19Z (GMT). No. of bitstreams: 1 Noeli Antonia Pimentel Vaz.pdf: 27939265 bytes, checksum: a285e763b07345fb2801fce0e9ed4cfe (MD5) Previous issue date: 2013-05-15 / The generation of a pattern of navigation in Virtual Learning Environments (VLE) closer to the characteristics of the student is a question relevant of research. This paper takes as its starting point a Hybrid Intelligent Tutoring System (HITS) already proposed in the literature that uses the Artificial Neural Networks MLP (Multi Layer Perceptron) and rules of faculty experts, and adds a methodology using Partial Differential Equations (PDE) calculated by the method of finite differences, in order to improve the navigation decision to lead the student in a personalized way, adding to their biggest gains learning process. The original HITS developed for web platform has been adapted using technologies that enable the evolution system and the integration of the system settings for navigation mid-level, system without adaptive capacity, and the smart navigation were conducted to verify the effectiveness of the proposed system. Data obtained were evaluated by statistical tests with 5% level of significance, which proved the efficiency of the proposed model. / A geração de um padrão de navegação em Ambientes Virtuais de Aprendizagem (AVA) mais próximo das características do estudante é uma questão relevante de pesquisa. Esta dissertação toma como ponto de partida um Sistema Tutor Inteligente (STI) Híbrido já proposto na literatura que utiliza Redes Neurais Artificiais do tipo MLP (Multi Layer Perceptron) e regras de docentes especialistas; e agrega uma metodologia utilizando Equações Diferenciais Parciais (EDP) calculadas pelo método das diferenças finitas, com o objetivo de melhorar a decisão de navegação, para conduzir o estudante de forma personalizada, agregando melhores ganhos ao seu processo de aprendizado. O STI original, desenvolvido para plataforma web, foi adaptado utilizando tecnologias que possibilitem a evolução do sistema e a integração a Ambientes Virtuais de Aprendizagem, como o Moodle, por exemplo. Experimentos com o sistema em configurações para a navegação nível médio (sistema sem capacidade adaptativa) e a navegação inteligente foram realizados para verificação da efetividade do sistema proposto. Os dados obtidos foram avaliados por meio de testes estatísticos com nível de significância de 5%, que comprovam a eficiência do modelo proposto.

Este trabalho tem como objetivo desenvolver um método numérico para simular escoamentos incompressíveis, isotérmicos, confinados ou com superfícies livres, de fuidos viscoelásticos governados pelos modelos integrais de Maxwell e K-BKZ (Kaye-Bernstein, Kearsley e Zapas). A técnica numérica apresentada é uma extensão do método GENSMAC (Tomé McKee - J. Comp. Phys., (110), pp 171--186, 1994 ) para a solução das equações de conservação, juntamente com as equações constitutivas integrais de Maxwell e K-BKZ. As equações governantes são resolvidas pelo método de diferenças finitas em uma malha deslocada. O tensor de Finger, B_t\'(t) é calculado com base nas idéias do método de campos de deformação (Peters et al. - J. Non-Newtonian Fluid Mech. (89), de maneira que não há a necessidade de seguir a trajetória da partícula de fuido para descrever a história de deformação da partícula. Uma abordagem diferente para a discretização do instante passado é utilizada e o tensor de Finger e o tensor das tensões são calculados utilizando um método de segunda ordem. A validação do método numérico descrito nesse trabalho foi feita utilizando o escoamento em um canal bidimensional e a solução numérica obtida para a velocidade e para as componentes de tensão com o modelo de Maxwell foram comparadas com as respectivas soluções analíticas no estado estacionário, mostrando excelente concordância. Os resultados numéricos para a simulação do escoamento em uma contração planar 4 : 1 mostraram bons resultados, tanto qualitativos quanto quantitativos, quando comparados com os resultados experimentais de Quinzani et al. ( J. Non-Newtonian Fluid Mech. (52), pp 1?36, 1994 ). Além disso, utilizando os modelos Maxwel e K-BKZ, o escoamento em uma contração planar 4 : 1 foi simulado para vários números de Weissenberg e os resultados obtidos estão de acordo com os encontrados na literatura. Resultados numéricos de escoamentos com superfícies livres modelados pelas equações integrais de Maxwell e K-BKZ são apresentados. Em particular, a simulação numérica do jato oscilante para diferentes números de Weissenberg e diferentes números de Reynolds é apresentada. / The aim of this work is to develop a numerical technique for simulating incompressible, isothermal, free surface (also con¯ned) viscoelastic flows of fuids governed by the integral models of Maxwell and K-BKZ (Kaye-Bernstein, Kearsley and Zapas). The numerical technique described herein is an extension of the GENSMAC method (Tome and McKee, J. Comput. Phys., 110, pp. 171-186, 1994) to the solution of the momentuum and mass conservation equations together with the integral constitutive Maxwell and K-BKZ equations. The governing equations are solved by the finite difference method on a staggered grid using a Marker-and-Cell approach. The fluid is represented by marker particles on the fluid surface only. This provides the visualization and location of the fluid free surface so that the free surface stress conditions can be applied. The Finger tensor Bt0(t) is computed using the ideias of the deformation fields method (Peters et al. J. Non-Newtonian Fluid Mech., 89, pp. 209-228, 2001) so that it is not necessary to track a fluid particle in order to calculate its deformation history. However, in this work modifcations to the deformation fields method are introduced: the past time is discretized using a different formula, the Finger tensor Bt0(x; t) is obtained by a second order method and the stress tensor ? (x; t) is computed by a second order quadrature formula. The numerical method presented in this work is validated by simulating the flow of a Maxwell fluid in a two-dimensional channel and the numerical solutions of the velocity and the stress components are compared with the respective analytic solutions providing a good agreement. Further, the flow through a 4:1 planar contraction of a specific fuid studied experimentally by Quinzani et al. (J. Non-Newtonian Fluid Mech., 52, pp. 1-36, 1994) was simulated and the numerical results were compared qualitatively and quantitatively with the experimental results and very good agreement was obtained. The Maxwell and the K-BKZ models were applied to simulate the 4:1 planar contraction problem using various Weissenberg numbers and the numerical results were in agreement with those published in the literature. Finally, numerical results of free surface flows using the Maxwell and K-BKZ integral constitutive equations are presented. In particular, the numerical simulation of jet buckling using several Weissenberg numbers and various Reynolds numbers are presented

Para o projeto de estruturas com perfis de aço formados a frio, é fundamental a compreensão dos fenômenos da instabilidade local e global, uma vez que estes apresentam alta esbeltez e baixa rigidez à torção. A determinação do carregamento crítico e a identificação do modo de instabilidade contribuem para o entendimento do comportamento dessas estruturas. Este trabalho avalia três metodologias para a análise linear de estabilidade de perfis de aço formados a frio isolados, com o objetivo de determinar os carregamentos críticos elásticos de bifurcação e os modos de instabilidade associados. Estritamente, analisa-se perfis de seção U enrijecido e Z enrijecido isolados, de diversos comprimentos e diferentes condições de vinculação e carregamento. Determinam-se os carregamentos críticos elásticos de bifurcação e os modos de instabilidade globais e locais por meio de: (i) análise com o Método das Faixas Finitas (MFF), através do uso do programa computacional CUFSM; (ii) análise com elementos finitos de barra baseados na Teoria Generalizada de Vigas (MEF-GBT), via uso do programa GBTUL; e (iii) análise com elementos finitos de casca (MEF-cascas) por meio do uso do programa ABAQUS. Algumas restrições e ressalvas com relação ao uso do MFF são apresentadas, assim como limitações da Teoria Generalizada de Viga e precauções a serem tomadas nos modelos de cascas. Analisa-se também a influência do grau de discretização da seção transversal. No entanto, não é feita avaliação em relação aos procedimentos normativos e tampouco análises não lineares, considerando as imperfeições geométricas iniciais, tensões residuais e o comportamento elastoplástico do material. / For the design of cold formed steel members, it is essential to understand the effects of local and global instability, since these members typically have a high slenderness and low torsion stiffness. The determination of critical loads and the associated buckling modes contribute to understand the behavior of these members. This work performs a evaluation of three methods for linear stability analysis of isolated cold-formed steel members in order to determine the elastic critical loads and the corresponding buckling modes. Specifically, Ue and Ze shape members were studied with various length, different boundary conditions and loads. The elastic critical loads and buckling modes are determined by means of: (i) analysis with the Finite Strip Method (FSM), by the computer program CUFSM, (ii) beam finite element analysis based on the Generalized Beam Theory (FEM-GBT), by GBTUL program, and (iii) Finite Element Method with shell analysis using ABAQUS program. Some restrictions and warnings regarding the use of the FSM are presented, as well as limitations of the Generalized Beam Theory and precautions to be taken in the shell models. It is also analyzed the influence of the degree of discretization of the cross section. In the present study, no evaluation was made with respect to normative procedures neither nonlinear analyses considering the initial geometric imperfections, residual stresses and elastoplastic behavior of the material.