A Backwards-Tracking Lagrangian-Eulerian Method for Viscoelastic Two-Fluid Flows

A new Lagrangian–Eulerian method for the simulation of viscoelastic free surface flow is proposed. The approach is developed from a method in which the constitutive equation for viscoelastic stress is solved at Lagrangian nodes, which are convected by the flow, and interpolated to the Eulerian grid with radial basis functions. In the new method, a backwards-tracking methodology is employed, allowing for fixed locations for the Lagrangian nodes to be chosen a priori. The proposed method is also extended to the simulation of viscoelastic free surface flow with the volume of fluid method. No unstructured interpolation or node redistribution is required with the new approach. Furthermore, the total amount of Lagrangian nodes is significantly reduced when compared to the original Lagrangian–Eulerian method. Consequently, the method is more computationally efficient and robust. No additional stabilization technique, such as both-sides diffusion or reformulation of the constitutive equation, is necessary. A validation is performed with the analytic solution for transient and steady planar Poiseuille flow, with excellent results. Furthermore, the proposed method agrees well with numerical data from the literature for the viscoelastic die swell flow of an Oldroyd-B model. The capabilities to simulate viscoelastic free surface flow are also demonstrated through the simulation of a jet buckling case.

[1]  H. Marschall,et al.  Numerical simulation of viscoelastic two-phase flows using openFOAM® , 2011 .

[2]  Andreas Mark,et al.  A Lagrangian-Eulerian framework for simulation of transient viscoelastic fluid flow , 2019, Journal of Non-Newtonian Fluid Mechanics.

[3]  Transient 3D finite element method for predicting extrudate swell of domains containing sharp edges , 2019, Journal of Non-Newtonian Fluid Mechanics.

[4]  M. F. Tomé,et al.  A marker-and-cell approach to viscoelastic free surface flows using the PTT model , 2007 .

[5]  Manuel Laso,et al.  Numerical simulation of 3D viscoelastic flows with free surfaces , 2006, J. Comput. Phys..

[6]  F. Pinho,et al.  Numerical simulation of non-linear elastic flows with a general collocated finite-volume method , 1998 .

[7]  Johan S. Carlson,et al.  Process Simulation and Automatic Path Planning of Adhesive Joining , 2016 .

[8]  Henrik Koblitz Rasmussen,et al.  Simulation of transient viscoelastic flow with second order time integration , 1995 .

[9]  R. Keunings,et al.  The adaptive Lagrangian particle method for macroscopic and micro–macro computations of time-dependent viscoelastic flows , 1999 .

[10]  Patrick Patrick Anderson,et al.  Sintering of Two Viscoelastic Particles: A Computational Approach , 2017 .

[11]  M. F. Tomé,et al.  Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows , 2011 .

[12]  José Alberto Cuminato,et al.  An implicit technique for solving 3D low Reynolds number moving free surface flows , 2008, J. Comput. Phys..

[13]  Andreas Mark,et al.  Simulations of 3D bioprinting: predicting bioprintability of nanofibrillar inks , 2018, Biofabrication.

[14]  R. Keunings,et al.  The Lagrangian particle method for macroscopic and micro-macro viscoelastic flow computations , 1998 .

[15]  Vincent Legat,et al.  The backward-tracking Lagrangian particle method for transient viscoelastic flows , 2000 .

[16]  J. O. Cruickshank.  Low-Reynolds-number instabilities in stagnating jet flows , 1988, Journal of Fluid Mechanics.

[17]  V. G. Ferreira,et al.  A finite difference technique for solving the Oldroyd-B model for 3D-unsteady free surface flows , 2008 .

[18]  Erik Svenning,et al.  Simulation of a highly elastic structure interacting with a two-phase flow , 2012 .

[19]  Andreas Mark,et al.  Multicriteria Optimization of an Oven With a Novel ε-Constraint-Based Sandwiching Method , 2021 .

[20]  R. I. Issa,et al.  A Method for Capturing Sharp Fluid Interfaces on Arbitrary Meshes , 1999 .

[21]  M. Muradoglu,et al.  A front-tracking method for computational modeling of viscoelastic two-phase flow systems , 2015 .

[22]  Robert Rundqvist,et al.  Comparison Between Different Immersed Boundary Conditions for Simulation of Complex Fluid Flows , 2011 .

[23]  Erik Svenning,et al.  An immersed boundary method for simulation of flow with heat transfer , 2013 .

[24]  Andreas Mark,et al.  Computationally efficient viscoelastic flow simulation using a Lagrangian-Eulerian method and GPU-acceleration , 2020, Journal of Non-Newtonian Fluid Mechanics.

[25]  J. P. V. Doormaal,et al.  ENHANCEMENTS OF THE SIMPLE METHOD FOR PREDICTING INCOMPRESSIBLE FLUID FLOWS , 1984 .

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

[27]  M. F. Tomé,et al.  Numerical solution of the PTT constitutive equation for unsteady three-dimensional free surface flows , 2010 .

[28]  Dieter Bothe,et al.  An extended volume of fluid method and its application to single bubbles rising in a viscoelastic liquid , 2018, J. Comput. Phys..

[29]  R. Tanner.  A theory of die‐swell , 1970 .

[30]  Raanan Fattal,et al.  Flow of viscoelastic fluids past a cylinder at high Weissenberg number : stabilized simulations using matrix logarithms , 2005 .

[31]  Andreas Mark,et al.  Derivation and validation of a novel implicit second-order accurate immersed boundary method , 2008, J. Comput. Phys..

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

[33]  F. Baaijens,et al.  Viscoelastic flow past a confined cylinder of a polyisobutylene solution , 1995 .

[34]  Oliver G. Harlen,et al.  A split Lagrangian-Eulerian method for simulating transient viscoelastic flows , 1995 .

[35]  N. Phan-Thien,et al.  Numerical modelling of transient viscoelastic flows , 2004 .

[36]  Tommy Andersson,et al.  Multiobjective Optimization of a Heat-Sink Design Using the Sandwiching Algorithm and an Immersed Boundary Conjugate Heat Transfer Solver , 2018 .

[37]  F. Pinho,et al.  Benchmark solutions for the flow of Oldroyd-B and PTT fluids in planar contractions , 2003 .

[38]  Fernando T. Pinho,et al.  The flow of viscoelastic fluids past a cylinder : finite-volume high-resolution methods , 2001 .

[39]  R. Tanner.  A theory of die-swell revisited , 2005 .

[40]  Robert Bohlin,et al.  Optimisation of robotised sealing stations in paint shops by process simulation and automatic path planning , 2014, Int. J. Manuf. Res..

[41]  Johan S. Carlson,et al.  Math-Based Algorithms and Software for Virtual Product Realization Implemented in Automotive Paint Shops , 2017 .

[42]  J. Hattel,et al.  Numerical simulation of the planar extrudate swell of pseudoplastic and viscoelastic fluids with the streamfunction and the VOF methods , 2018 .

[43]  M. F. Tomé,et al.  Application of the log-conformation tensor to three-dimensional time-dependent free surface flows , 2012 .