Modelling sediment resuspension in industrial tanks using SPH

Remise en suspension modélisation en citernes industriels qui utilisent des SPH

University of Manchester
email : georgios.fourtakas@postgrad.manchester.ac.uk, benedict.rogers@manchester.ac.uk, dominique.laurence@manchester.ac.uk

Abstract

The present paper uses the weakly compressible smoothed particle hydrodynamics (SPH) to model numerically the sediment suspension and the shear layer induced by rapid flows resulting in scouring of the sediment with two yield criteria based on the Mohr-Coulomb and Drucker-Prager formulations. To the best of the authors’ knowledge this is the first time that the two yield criteria have been rigorously assessed in SPH where they are both expressed in terms of Coulomb-type parameters. The sediment is treated as a pseudo-Newtonian fluid where two different conditions distinguish between the state of the sediment – either still or suspended- according to the yield criteria. A variable viscosity is used to depict accurately the suspension layer. The models have been implemented using the DualSPHysics code, an SPH solver capable of dealing with high deformations and free surface flows. Three validation cases have been used to assess the validity of the yield criteria. The still sediment/water case was used to investigate the stability of the interface and the interaction with the boundaries that produced acceptable results for both criteria. The global relative error in the still case was used to assess the yield behaviour of the solid phase with both models performing satisfactorily. It was shown that the Mohr-Coulomb criterion was more conservative in particle movement. The second validation case was the tangential flow between two co-rotating co-axial cylinders using an analytical solution to confirm the rheological behaviour of the viscous sediment. With this case it was shown that the variable viscosity of the Drucker-Prager is lower than the Mohr-Coulomb model making the sediment less viscous. Finally, the erosional dam break was used to visualize the scour induced by the fluid. Following on from the results of the previous test cases in terms of variable viscosity and interface interaction the scour predicted with the Drucker-Prager model was larger and closer to the experimental results with a larger and deeper shear layer capturing the dune and scour. The Mohr-Coulomb yield criterion appears to underpredict the scouring process in addition to a small lag in the toe front of the solid phase that may be accounted to the smaller shear layer of the sediment. A sensitivity analysis on the Coulomb parameters and a convergence study must be performed in order to assess any deviations from the results discussed herein.

Résumé

Ce papier présente une discrétisation particulaire Lagrangienne (SPH) avec hypothèse faiblement compressible de particules (SPH smoothed particle hydrodynamics) étendue à la modélisation de sédiments en suspension ainsi que le cisaillement dans la couche de dépôt induit par des écoulements rapides entrainant le récurage du sédiment, et compare deux critères Newtonien-plastiques de Mohr-Coulomb et de Drucker-Prager. A notre connaissance ceci est la première fois que les deux critères de rendement sont rigoureusement évalués dans la méthode SPH où les deux états du sédiment, pseudo fluide-newtonien ou statique, sont modélisés. Une viscosité variable est employée pour représenter précisément la couche en suspension. Les modèles ont été mis en oeuvre le code de DualSPHysics, un solver SPH capable de traiter des déformations fortes et des surfaces libres. Trois cas de validation permettent d’évaluer la validité des critères de fluidification. Le cas statique de réservoir sédiment/eau pour étudier la stabilité de l’interface et l’interaction avec les frontières a produit des résultats satisfaisant pour les deux critères en termes d’erreur relative globale. On montre que le critère de Mohr-Coulomb était plus conservateur dans le mouvement de particules. Le deuxième cas de validation est l’écoulement entre deux cylindres coaxiaux co-tournants, qui a une solution analytique. Il a été employé pour confirmer le comportement rhéologique du sédiment visqueux. Dans ce cas on observe que la viscosité variable du modèle Drucker-Prager est plus faible que dans le modèle de Mohr-Coulomb rendant le sédiment moins visqueux. Enfin la rupture de barrage avec érosion a été employée pour visualiser le lessivage induit par le fluide. En cohérence avec les tests précédents, en termes de viscosité variable et d’interaction à l’interface de lessivage, le modèle de Drucker-Prager était dynamique et plus près des résultats expérimentaux avec une plus grande et plus profonde couche de cisaillement capturant la formation dunaire et le lessivage. Le critère de Mohr-Coulomb semble sous-estimer le processus de récurage et présente un léger retard à l’avant du front de produits lessivés qui peut être du au plus faible cisaillement du sédiment. Une analyse de sensibilité sur les paramètres de Coulomb et une étude de convergence doivent être réalisées afin d’évaluer les déviations des résultats discutés ci-dessus.