Numerical solution of the FENE-CR model in complex flows

Autor(es)

Paulo, G. S. ^[1] ; Oishi, C. M. ^[1] ; Tome, M. F. ^[2] ; Alves, M. A. ^[3] ; Pinho, F. T. ^[4]

Afiliação do(s) autor(es)

^[1] Univ Estadual Paulista, Dept Matemat & Comp, BR-19060900 Presidente Prudente - Brazil

^[2] Univ Sao Paulo, Dept Matemat Aplicada & Estat, Sao Carlos, SP - Brazil

^[3] Univ Porto, Fac Engn, Dept Engn Quim, Ctr Estudos Fenomenos Transporte, P-4200465 Oporto - Portugal

^[4] Univ Porto, Fac Engn, Dept Engn Mecan, Ctr Estudos Fenomenos Transporte, P-4200465 Oporto - Portugal

• Tipo de documento: Artigo Científico
• Fonte: Journal of Non-Newtonian Fluid Mechanics; v. 204, p. 50-61, FEB 2014.
• Citações Web of Science: 9

Resumo

A finite difference technique for solving the FENE-CR (Finite Extendable Non-linear Elastic - Chilcott and Rallison) closure constitutive model in complex flows has been developed and tested. The governing equations are solved using a Marker-and-Cell type method on a staggered grid. The momentum equation is integrated employing an implicit method while the FENE-CR constitutive equation is approximated by a second-order Runge-Kutta scheme. To demonstrate that the numerical technique can cope with complex flows governed by the FENE-CR model, three flow problems were analysed: the fully-developed channel flow, the 2D cross-slot flow and the impacting drop problem. The analytic solution for fully-developed channel flow of FENE-CR fluids with a solvent viscosity is also presented for validation purposes. This flow problem is used to verify the numerical method and to quantify its accuracy by comparing numerical results of fully-developed channel flow with the analytic solution. The second flow is employed to assess whether the numerical methodology is capable of capturing the purely-elastic instabilities predicted in the literature for 2D cross-slot confined flows. Additionally, the complex free surface flow corresponding to the filling of a 20 cross geometry has also been investigated. The last problem concerns the flow dynamics of a FENE-CR fluid drop impacting on a rigid surface, which allows the assessment of the capability of the model to deal with free surfaces. The effects of varying the Reynolds number, the Weissenberg number and the finite extensibility of the polymer molecules (L-2) on the resulting flow patterns are analysed. (C) 2013 Elsevier B.V. All rights reserved. (AU)

• Processo FAPESP: 13/07375-0 - CeMEAI - Centro de Ciências Matemáticas Aplicadas à Indústria
• Beneficiário: José Alberto Cuminato
• Linha de fomento: Auxílio à Pesquisa - Centros de Pesquisa, Inovação e Difusão - CEPIDs
• Processo FAPESP: 11/13930-0 - Desenvolvimento de um método numérico para simular escoamentos transientes viscoelásticos com superfícies livres
• Beneficiário: Gilcilene Sanchez de Paulo
• Linha de fomento: Auxílio à Pesquisa - Regular
• Processo FAPESP: 09/15892-9 - Estudo de métodos numéricos estáveis e acurados em escoamentos transientes: aperfeiçoamentos, implementações, aplicações com superfícies livres e modelos viscoelásticos
• Beneficiário: Cassio Machiaveli Oishi
• Linha de fomento: Auxílio à Pesquisa - Apoio a Jovens Pesquisadores