A New Method to Simulate Free Surface Flows for Viscoelastic Fluid

Free surface flows arise in a variety of engineering applications. To predict the dynamic characteristics of such problems, specific numerical methods are required to accurately capture the shape of free surface. This paper proposed a new method which combined the Arbitrary Lagrangian-Eulerian (ALE) technique with the Finite Volume Method (FVM) to simulate the time-dependent viscoelastic free surface flows. Based on an open source CFD toolbox called OpenFOAM, we designed an ALE-FVM free surface simulation platform. In the meantime, the die-swell flow had been investigated with our proposed platform to make a further analysis of free surface phenomenon. The results validated the correctness and effectiveness of the proposed method for free surface simulation in both Newtonian fluid and viscoelastic fluid.

[1]  Hrvoje Jasak,et al.  OpenFOAM: Open source CFD in research and industry , 2009 .

[2]  M. F. Tomé,et al.  A finite difference technique for simulating unsteady viscoelastic free surface flows , 2002 .

[3]  H. R. Tamaddon-Jahromi,et al.  Simulation of viscoelastic and viscoelastoplastic die-swell flows , 2013 .

[4]  N. Phan-Thien,et al.  Viscoelastic finite volume method , 1999 .

[5]  S. Osher,et al.  Algorithms Based on Hamilton-Jacobi Formulations , 1988 .

[6]  G. Georgiou,et al.  Converged solutions of the Newtonian extrudate-swell problem , 1999 .

[7]  C. W. Hirt,et al.  An Arbitrary Lagrangian-Eulerian Computing Method for All Flow Speeds , 1997 .

[8]  H. Jasak,et al.  A moving mesh finite volume interface tracking method for surface tension dominated interfacial fluid flow , 2012 .

[9]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[10]  Caicheng Lu,et al.  Incomplete LU preconditioning for large scale dense complex linear systems from electromagnetic wave scattering problems , 2003 .

[11]  H. Jasak,et al.  Automatic mesh motion for the unstructured finite volume method , 2007 .

[12]  G. Georgiou,et al.  A study of various factors affecting Newtonian extrudate swell , 2012 .

[13]  V. G. Ferreira,et al.  The MAC method , 2008 .

[14]  T. Phillips,et al.  Spectral element predictions of die-swell for Oldroyd-B fluids , 2011 .

[15]  Michel Fortin,et al.  A finite element procedure for viscoelastic flows , 1987 .

[16]  Xiao-Wei Guo,et al.  Non-equilibrium steady states of entangled polymer mixtures under shear flow , 2015 .

[17]  S. Muzaferija.  Computation of free-surface flows using the finite-volume method and moving grids , 1997 .

[18]  M. F. Webster,et al.  Recovery and stress-splitting schemes for viscoelastic flows , 1998 .

[19]  M. A. Ajiz,et al.  A robust incomplete Choleski‐conjugate gradient algorithm , 1984 .

[20]  R. Keunings,et al.  Finite-element Analysis of Die Swell of a Highly Elastic Fluid , 1982 .

[21]  Shun Zou,et al.  Interface instabilities and chaotic rheological responses in binary polymer mixtures under shear flow , 2014 .

[22]  Hrvoje Jasak,et al.  Viscoelastic fluid analysis in internal and in free surface flows using the software OpenFOAM , 2010, Comput. Chem. Eng..

[23]  S. Claus.  Numerical simulation of complex viscoelastic flows using discontinuous galerkin spectral/hp element methods , 2013 .

[24]  J. Clermont,et al.  Numerical simulation of extrudate swell for Oldroyd-B fluids using the stream-tube analysis and a streamline approximation , 1993 .