A compressible 3D finite volume approach for the simulation of unsteady viscoelastic cavitating flows. (May 2022)

• Record Type: Journal Article
• Authors: Lang, Christian; Boolakee, Oliver; Schmidt, Steffen J.; Adams, Nikolaus A.
• Is Part Of: International journal of multiphase flow. Volume 150(2022)
• Publication Date: 2022-05
• Subjects: Viscoelasticity -- Cavitation -- Compressible flow -- Objective time derivative -- Stress tensor symmetry boundary condition -- Bubble dynamics; Multiphase flow -- Periodicals; Écoulement polyphasique -- Périodiques; Multiphase flow; Periodicals; 620.1064
• DOI: 10.1016/j.ijmultiphaseflow.2022.103981
• Languages: English
• ISSNs: 0301-9322

Abstract

Abstract: We present a fully compressible, density-based finite volume solver for the simulation of 3D cavitating flows in viscoelastic Maxwell-/Oldroyd-like fluids. The upper-convected Maxwell model, the Oldroyd-B model and the linear and exponential simplified Phan-Thien Tanner models are implemented as viscoelastic constitutive equations in conservative formulation, and we identified the Truesdell rate as appropriate objective time derivative for compressible flows. Cavitation is modeled by a single-fluid homogeneous mixture equilibrium approach considering condensation and evaporation assuming volume averaged mixture quantities. The corresponding simplified quasilinear system is analyzed, and the wave speeds are calculated in order to adapt the employed four step Runge–Kutta explicit time stepping concerning the viscoelastic transport equations. We introduce a novel ghost cell boundary condition for the viscoelastic stress tensor. The approach is tested against (semi-)analytical unsteady and steady-state references and shows very good agreement. 3D simulations of the spherical vapor bubble collapse are performed for all implemented viscoelastic models and show a distinct influence compared to the Newtonian case. For the upper-convected Maxwell fluid a variation of the relaxation time exhibits its perspicuous influence on the dynamics of the collapse. Highlights: Compressible density based finite volume solver for viscoelastic cavitating flows. UCM, Oldroyd-B, simplified LPTT and EPTT model and compressible objective rate. Single-fluid homogeneous mixture equilibrium cavitation model. Novel symmetry ghost cell boundary condition for the viscoelastic stress tensor. Viscoelasticity significantly changes the spherical vapor bubble collapse.