Submitted by RODRIGO KOITI ISHIZAKA (rodrigo_ishizaka@hotmail.com) on 2018-04-09T20:10:27Z No. of bitstreams: 1 Dissertação Final PDF.pdf: 10461391 bytes, checksum: a78974d2db3d166540e03d4cee7def37 (MD5) / Rejected by ALESSANDRA KUBA OSHIRO ASSUNÇÃO (alessandra@fct.unesp.br), reason: Solicitamos que realize correções na submissão seguindo as orientações abaixo: - Número de página das referências no sumário está incorreto (consta p.86 quando é p.87); - No sumário consta Apêndices, porém no trabalho não. É finalizado nas referências. Agradecemos a compreensão. on 2018-04-10T12:44:45Z (GMT) / Submitted by RODRIGO KOITI ISHIZAKA (rodrigo_ishizaka@hotmail.com) on 2018-04-10T16:09:35Z No. of bitstreams: 1 Dissertação Final PDF.pdf: 10461372 bytes, checksum: 679a6e28fa80deb9079156b9b6b70c36 (MD5) / Approved for entry into archive by ALESSANDRA KUBA OSHIRO ASSUNÇÃO (alessandra@fct.unesp.br) on 2018-04-10T20:02:48Z (GMT) No. of bitstreams: 1 ishizaka_rk_me_prud.pdf: 10461372 bytes, checksum: 679a6e28fa80deb9079156b9b6b70c36 (MD5) / Made available in DSpace on 2018-04-10T20:02:48Z (GMT). No. of bitstreams: 1 ishizaka_rk_me_prud.pdf: 10461372 bytes, checksum: 679a6e28fa80deb9079156b9b6b70c36 (MD5) Previous issue date: 2018-02-28 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / O presente trabalho de mestrado consiste em apresentar uma modelagem para escoamentos isotérmicos viscoelásticos em que a viscosidade varia de acordo com a pressão e estudos numéricos para alguns problemas bidimensionais, como escoamentos entre placas paralelas e expansão planar 1:4. Para este trabalho é utilizado o modelo viscoelástico Oldroyd-B e uma modelagem linear da viscosidade com relação a pressão. O método numérico desenvolvido é baseado no método da projeção para desacoplar velocidade e pressão nas equações de Navier-Stokes e depois calcula os tensores pela equação constitutiva com a informação da viscosidade variando com a pressão. As equações são discretizadas em uma malha deslocada pelo método de diferenças finitas. Os resultados numéricos são obtidos das simulações de escoamentos em um canal e expansão planar 1:4 bidimensional, cujo foco é destacar algumas diferenças entre o modelo Oldroyd-B Standard e o modelo Oldroyd-B sob influência da viscosidade que varia linearmente com a pressão. Este trabalho de mestrado também tem por objetivo resolver problemas, com fluidos viscoelásticos, em que haja regiões com variações da pressão e estudar estes resultados numéricos comparando-os com o modelo Oldroyd-B Standard. / The present master's work consists in presenting a model for viscoelastic isothermal flows with pressure-dependent viscosity and numerical studies for some two-dimensional problems, such as channel flow and planar expansion 1:4. In this work will be considered the Oldroyd-B viscoelastic model and a linear modeling viscosity with a linear relationship between viscosity and pressure. The numerical method developed is based on the projection method to decouple velocity and pressure in the Navier-Stokes equations and then calculates the tensor from the constitutive equation taking account the pressure-dependent viscosity. The equations are discretized on a staggered grid by using the finite diference method. The numerical results are obtained for two problems: two-dimensional channel and planar expansion 1:4 flows, whose focus is to highlight some diferences between the standard Oldroyd-B model and the Oldroyd-B model under the influence of viscosity that varies linearly with pressure. This master's work aims to solve problems of viscoelastic fluids in which there exist regions where pressure varies and to study these numerical results comparing them with the standard Oldroyd-B model.

[pt] Filmes finos de líquidos estão presentes em uma gama de aplicações industriais, como processos de atomização e revestimento de substrato. O processo de quebra pode ser divido em duas etapas: o estágio de ruptura, e o estágio de retração. O primeiro, movido pelas forças de van der Waals, ocorre quando uma pequena perturbação cresce e provoca o aparecimento de um pequeno furo no filme. O segundo, movido por forças capilares, provoca o crescimento desse furo levando à desintegração do filme de líquido. A estabilidade de uma cortina de líquido depende das características da perturbação, da espessura do filme e das propriedades do fluido. Análises experimentais mostraram que uma cortina super fina pode ser obtida pela utilização de fluidos viscoelásticos. Os mecanismos físicos associados à esta estabilidade, contudo, não são totalmente compreendidos. Este trabalho apresenta um estudo numérico e teórico dos efeitos das propriedades viscoelásticas na estabilidade de uma cortina de fluido, englobando ambos os estágio do processo. As análises numéricas foram desenvolvidas através da expansão assintótica das variáveis do escoamento com aplicação de um esquema de integração no tempo totalmente implícito. A partir da análise teórica da dinâmica de ruptura foi possível obter um critério de estabilidade linear para perturbações planares e axissimétricas em fluidos Newtonianos e não-Newtonianos. O tempo de ruptura e a velocidade de retração do filme foram calculados numericamente como função das propriedades viscoelásticas do líquido. Resultados mostraram que as forças elásticas atuam de forma a dificultar o processo de quebra e retração. Análises da evolução da espessura mostraram que as propriedades reológicas do fluído também interferem no formato que o filme de fluido assume durante o processo de retração. Para regimes de baixa viscosidade, as forças elásticas atuaram evitando a formação de ondas capilares observadas em fluidos Newtonianos. / [en] Thin free liquid sheets are ubiquitous in many industrial processes, such as atomization and curtain coating. Liquid sheets are susceptible to instabilities at the interface, which can grow, triggering a breakup process. This process can be divided into two different stages: the rupture stage and retraction. The first, driven by van der Waals force, occurs when a small instability grows until it pinches-off the sheet. The second, driven by capillary forces, induces the growth of the hole caused by the pinch-off, leading to the full disintegration of the liquid sheet. The stability of a liquid sheet depends on disturbance characteristics, sheet thickness, and fluid properties. Experimental analyses have shown that thinner stable liquid curtain can be obtained with viscoelastic liquids. The underlyning physical mechanisms associated with increased stability are, however, not fully understood. This work presents a theoretical and numerical analysis of the effect of viscoelasticity on the stability of a thin liquid sheet during both stages of the breakup process. We first analyze the rupture dynamics, deriving linear stability criteria for both planar and axisymmetric perturbations of Newtonian and Oldroyd-B liquids. The time evolution of planar and axisymmetric perturbations in an Oldroyd-B liquid sheet is evaluated using the asymptotic expansion of the flow variables and a fully-implicit time integration scheme. The rupture time and retraction velocity are calculated as a function of the viscoelastic properties. The results show that the liquid rheological behavior does not influence the linear stability criterion. Nevertheless, it has a strong effect on the growth rate of the disturbance and retraction velocity, increasing, thus, the breakup time. The results show that elastic forces act to hinder the rupture and retraction stages. Analysis of the temporal evolution of the thickness profile reveals that liquid rheological behavior also affects the shape of the liquid sheet. For low viscosity regime, the elastic forces damp the capillary waves that arise during the retraction of Newtonian sheets.