Desenvolvimento de um método numérico para simular escoamentos viscoelásticos axissimétricos com superfícies livres: Modelo PTT

Neste trabalho é apresentado um método numérico para simular escoamentos viscoelásticos axissimétricos com superfícies livres modelados pela equação constitutiva PTT (Phan-Thien-Tanner). As equações governantes para escoamentos axissimétricos transientes incompressíveis de um fluido PTT são resolvidas utilizando o método de diferenças finitas numa malha deslocada. As derivadas temporais da equação do momento e da equação constitutiva são integradas pelo método de Euler explícito. Os termos convectivos são aproximados pelo método de alta ordem CUBISTA (Convergent and Universally Bounded Interpolation Scheme for the Treament of Advection), os termos difusivos e as derivadas espaciais são aproximados por diferenças centrais. O fluido é modelado utilizando a técnica Marker-and-Cell (MAC) o que permite visualizar e localizar a superfície livre do fluido. O método numérico descrito neste trabalho foi implementado na plataforma de alto desempenho denominada FREEFLOW-AXI e verificado comparando-se os resultados numéricos obtidos da simulação do escoamento totalmente desenvolvido em um tubo com a solução analítica deste problema. Foram obtidos alguns resultados numéricos para escoamentos axissimétricos com superfícies livres utilizando um fluido PTT, tais como, o impacto de uma gota esférica de fluido viscoelástico numa superfície rígida (Impacting Drop), o impacto de uma gota esférica numa superfície líquida (Splashing Drop) e o problema do inchamento do extrudado (Dieswell). Nestas simulações, o foco principal foi estudar numericamente os efeitos dos parâmetros do modelo PTT e nestes escoamentosIn this work is presented a numerical method for simulating axisymmetric viscoelastic free surface flows modeled by the constitutive equation PTT (Phan-Thien-Tanner). The governing equations for incompressible transient axisymmetric flows of a PTT fluid are solved using the finite difference method in a staggered grid. The time derivatives on momentum and constitutive equations are integrated by forward Euler method. The convective terms are approximated by the high order method CUBISTA (Convergent and Universally Bounded Interpolation Scheme for the Treament of Advection), the diffusive terms and the spatial derivatives are approximated by central differences. The fluid is modeled using the strategy “Marker-and-Cell” (MAC) which allows to view and locate the free surface of the fluid. The numerical method described in this work was implemented on the high performance platform called FREEFLOW-AXI and verified by comparing the numerical results of fully developed tube flow with the analytic solution of this problem. Some numerical results for axisymmetric free surface flows using a PTT fluid were obtained, such as, impacting drop, splashing drop and the dieswell problems. In these simulations, the numerical studies about the influence that the parameters of the PTT model and exert on the flows were the main focusFundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP