SEECCM III 3 rd South-East European Conference on Computational Mechanics- an ECCOMAS and IACM Special Interest Conference M. Papadrakakis, M. Kojic, I. Tuncer (eds.) Kos Island, Greece, 12–14 June 2013 SOLUTION OF VISCOPLASTIC FLOWS WITH THE FINITE VOLUME / MULTIGRID METHOD Alexandros Syrakos 1 , Georgios C. Georgiou 2 , and Andreas N. Alexandrou 3 1 Oceanography Centre, University of Cyprus PO Box 20537, 1678 Nicosia, Cyprus e-mail: syrakos.alexandros@ucy.ac.cy 2 Department of Mathematics and Statistics, University of Cyprus PO Box 20537, 1678 Nicosia, Cyprus e-mail: georgios@ucy.ac.cy 3 Department of Mechanical and Manufacturing Engineering, University of Cyprus PO Box 20537, 1678 Nicosia, Cyprus e-mail: andalexa@ucy.ac.cy Keywords: Finite Volume Method, Viscoplastic Flows, Papanastasiou Regularisation, Multi- grid, Lid-driven Cavity Abstract. We investigate the performance of the finite volume method in solving viscoplastic flows. The square lid-driven cavity flow of a Bingham plastic is chosen as the test case and the constitutive equation is regularised as proposed by Papanastasiou [J. Rheol. 31 (1987) 385- 404]. The numerical results obtained for Bingham numbers up to 10 and Reynolds numbers up to 5000 compare favourably with reported results obtained through other methods. The effects of the Reynolds and Bingham numbers are also investigated. It is shown that the convergence rate of the standard SIMPLE pressure-correction algorithm, which is used to solve the alge- braic equation system that is produced by the finite volume discretisation, severely deteriorates as the Bingham number increases, with a corresponding increase in the non-linearity of the equations. Using the SIMPLE algorithm in a multigrid context [Syrakos & Goulas, Int. J. Nu- mer. Methods Fluids 52 (2006) 1215-1245] dramatically improves convergence, although the multigrid convergence rates are not as high as those for Newtonian flows. 70