SPH simulations of thixo-viscoplastic fluid flow past a cylinder

• Published in: Journal of non-Newtonian fluid mechanics 2022-10, Vol.308, p.104891, Article 104891
• Main Authors: Rossi, E., Garcia de Beristain, I., Vazquez-Quesada, A., López-Aguilar, J.E., Ellero, M.
• Format: Article
• Language: English
• Subjects: SPH; Thixotropy; Viscoplasticity
• ISSN: 0377-0257; 1873-2631
• DOI: 10.1016/j.jnnfm.2022.104891

Thixotropic materials are complex fluids that display time-dependent viscosity and/or yield-stress response upon the application of a fixed deformation, while recovering their original structured-state when the deformation is discontinued. Thixotropic effects are presents in many different systems and applications, ranging from food products, such as ketchup, to metals, such as molten aluminum. In this work we present a first attempt to simulate the rheological properties of thixo-viscoplastic flows using a Smoothed Particle Hydrodynamic (SPH) method. The study set up is a 2D flow around a circular cylinder as well as a simple shear flow between parallel plates to validate our numerical results. SPH solutions are compared with simulations performed using the open-source Finite Volume Method solver RheoTool, based on OpenFOAM. The viscoplastic model used in this work is the Papanastasiou model combined with a recently developed microstructural one, in order to include thixotropy. In this thixo-viscoplastic framework, we analyze the flow properties in terms of yield-fronts, streamlines and structure-parameter fields at different Bingham and Thixotropy numbers, through microstructural thixotropic and yield-stress parameters variation. Obtained results show an important novelty: an asymmetry in the thixo-viscoplastic flow around the cylinder. •First attempt to simulate thixo-viscoplastic flows using SPH.•2D flow around a circular cylinder as well as a simple shear flow are analyzed.•Papanastasiou model is combined with a microstructural one to include thixotropy.•Flow properties are analysed using different Bingham and Thixotropy numbers.•Main finding is an asymmetry in the thixo-viscoplastic flow around a cylinder.