Flexible slender body fluid interaction: Vector-based discrete element method with Eulerian smoothed particle hydrodynamics

Abstract In this paper Eulerian weakly compressible SPH (EWCSPH) is coupled to the Vector-based discrete element model (V-model) to produce a robust mesh-less SPH-V model for fluid-structure interaction (FSI) applications including thin highly flexible bodies. The V-model, originally proposed as a model for granular structures, is extended herein to model non-linear deformation of a continuum representation of elastic structures. Non-linear static validation tests show that the extended V-model provides an accurate representation for problems involving large deformation. Undamped free vibration tests show the V-model is also accurate for transient problems. The SPH-V model with the immersed boundary method and SPH of Eulerian weakly compressible form is applied to established FSI tests for validation. These are oscillation of a flexible plate in quiescent fluids and flow past a cylinder with a rigidly connected flexible flat plate. The complex flow through flexible valve leaflets representative of flow through deep leg veins is demonstrated. Results thus show the SPH-V model is able to accurately capture complex flow features and also non-linear structural mechanics in fully coupled FSI problems.

[1]  X Zheng,et al.  A coupled sharp-interface immersed boundary-finite-element method for flow-structure interaction with application to human phonation. , 2010, Journal of biomechanical engineering.

[2]  Xiangmin Jiao,et al.  Fluid–structure interactions of the mitral valve and left heart: Comprehensive strategies, past, present and future , 2010, International journal for numerical methods in engineering.

[3]  JiSeok Lee,et al.  Fluid–structure interaction analysis on a flexible plate normal to a free stream at low Reynolds numbers , 2012 .

[4]  Alistair Revell,et al.  Eulerian weakly compressible smoothed particle hydrodynamics (SPH) with the immersed boundary method for thin slender bodies , 2019, Journal of Fluids and Structures.

[5]  A. Colagrossi,et al.  δ-SPH model for simulating violent impact flows , 2011 .

[6]  Michael E. Plesha,et al.  Discrete element method for modelling solid and particulate materials , 2007 .

[7]  F. Rosendaal,et al.  The Effect of Aging on Venous Valves , 2010, Arteriosclerosis, thrombosis, and vascular biology.

[8]  J. Morris,et al.  Modeling Low Reynolds Number Incompressible Flows Using SPH , 1997 .

[9]  C. Peskin.  Flow patterns around heart valves: A numerical method , 1972 .

[10]  T. Wick,et al.  Finite elements for fluid–structure interaction in ALE and fully Eulerian coordinates , 2010 .

[11]  Hitoshi Gotoh,et al.  Development of a fully Lagrangian MPS-based coupled method for simulation of fluid-structure interaction problems , 2014 .

[12]  Moubin Liu,et al.  Smoothed particle hydrodynamics and element bending group modeling of flexible fibers interacting with viscous fluids. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  Boyce E. Griffith.  Erratum: Immersed boundary model of aortic heart valve dynamics with physiological driving and loading conditions [International Journal for Numerical Methods in Biomedical Engineering 28,3 (2012) 317-345] DOI: 10.1002/cnm.1445 , 2013 .

[14]  Dominique Laurence,et al.  DNS and LES of 3-D wall-bounded turbulence using Smoothed Particle Hydrodynamics , 2015 .

[15]  J. Vierendeels,et al.  Performance of partitioned procedures in fluid-structure interaction , 2010 .

[16]  David Beeman,et al.  Some Multistep Methods for Use in Molecular Dynamics Calculations , 1976 .

[17]  K. Namkoong,et al.  Computation of dynamic fluid-structure interaction in two-dimensional laminar flows using combined formulation , 2005 .

[18]  D. Young,et al.  Tables of characteristic functions representing normal modes of vibration of a beam , 1949 .

[19]  Holger Wendland,et al.  Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree , 1995, Adv. Comput. Math..

[20]  S. Turek,et al.  Proposal for Numerical Benchmarking of Fluid-Structure Interaction between an Elastic Object and Laminar Incompressible Flow , 2006 .

[21]  Mihai Basa,et al.  Permeable and non‐reflecting boundary conditions in SPH , 2009 .

[22]  J. Monaghan,et al.  Shock simulation by the particle method SPH , 1983 .

[23]  Barbara M. Johnston,et al.  Non-Newtonian blood flow in human right coronary arteries: steady state simulations. , 2004, Journal of biomechanics.

[24]  H. Kolsky,et al.  Dynamics of vibrations , 1965 .

[25]  V. Kuzkin,et al.  Vector-based model of elastic bonds for simulation of granular solids. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[26]  Eli J Weinberg,et al.  Transient, Three-dimensional, Multiscale Simulations of the Human Aortic Valve , 2007, Cardiovascular engineering.

[27]  Annalisa Quaini,et al.  Extended ALE Method for fluid-structure interaction problems with large structural displacements , 2017, J. Comput. Phys..

[28]  Adrian R.G. Harwood,et al.  LUMA: A many-core, Fluid-Structure Interaction solver based on the Lattice-Boltzmann Method , 2018, SoftwareX.

[29]  Ernst Rank,et al.  Computation of fluid-structure interaction on lightweight structures , 2001 .

[30]  A. Colagrossi,et al.  Numerical simulation of interfacial flows by smoothed particle hydrodynamics , 2003 .

[31]  Abouzied Nasar,et al.  Eulerian and Lagrangian smoothed particle hydrodynamics as models for the interaction of fluids and flexible structures in biomedical flows , 2016 .

[32]  Jan Vierendeels,et al.  Fluid-Structure Interaction Simulation of Prosthetic Aortic Valves: Comparison between Immersed Boundary and Arbitrary Lagrangian-Eulerian Techniques for the Mesh Representation , 2016, PloS one.

[33]  D. V. Griffiths,et al.  Modelling of elastic continua using a grillage of structural elements based on discrete element concepts , 2001 .

[34]  Donald Stephen Dugdale,et al.  Elasticity for engineers , 1971 .

[35]  P. G. Jayathilake,et al.  Three-dimensional numerical simulations of human pulmonary cilia in the periciliary liquid layer by the immersed boundary method , 2012 .

[36]  Fedor Lurie,et al.  Mechanism of venous valve closure and role of the valve in circulation: a new concept. , 2003, Journal of vascular surgery.

[37]  Yucang Wang,et al.  A new algorithm to model the dynamics of 3-D bonded rigid bodies with rotations , 2009 .

[38]  K.A. Hoo,et al.  Numerical studies on the hemodynamics in the human vein and venous valve , 2008, 2008 American Control Conference.

[39]  A. Yu. Panchenko,et al.  Application of the pair torque interaction potential to simulate the elastic behavior of SLMoS2 , 2016 .

[40]  B. Rogers,et al.  Smoothed particle hydrodynamics method applied to pulsatile flow inside a rigid two‐dimensional model of left heart cavity , 2012, International journal for numerical methods in biomedical engineering.