A bi-projection method for Bingham type flows

We propose and study a new numerical scheme to compute the isothermal and unsteady flow of an incompressible viscoplastic Bingham medium. The main difficulty, for both theoretical and numerical approaches, is due to the non-differentiability of the plastic part of the stress tensor in regions where the rate-of-strain tensor vanishes. This is handled by reformulating the definition of the plastic stress tensor in terms of a projection. A new time scheme, based on the classical incremental projection method for the Newtonian Navier-Stokes equations, is proposed. The plastic tensor is treated implicitly in the first sub-step of the projection scheme and is computed by using a fixed point procedure. A pseudo-time relaxation is added into the Bingham projection whose effect is to ensure a geometric convergence of the fixed point algorithm. This is a key feature of the bi-projection scheme which provides a fast and accurate computation of the plastic tensor. Stability and error analyses of the numerical scheme are provided. The error induced by the pseudo-time relaxation term is controlled by a prescribed numerical parameter so that a first-order estimate of the time error is derived for the velocity field. A second-order cell-centred finite volume scheme on staggered grids is applied for the spatial discretisation. The scheme is assessed against previously published benchmark results for both Newtonian and Bingham flows in a two-dimensional lid-driven cavity for Reynolds number equals 1000. Moreover, the proposed numerical scheme is able to reproduce the fundamental property of cessation in finite time of a viscoplastic medium in the absence of any energy source term in the equations. For a fixed value (100) of the Bingham number, various numerical simulations for a range of Reynolds numbers up to 200000 were performed with the bi-projection scheme on a grid with 1024 2 mesh points. The effect of this (physical) parameter on the flow behaviour is discussed.

[1]  R. Glowinski,et al.  Numerical Analysis of Variational Inequalities , 1981 .

[2]  R. Glowinski,et al.  OPERATOR-SPLITTING METHODS FOR THE SIMULATION OF BINGHAM VISCO-PLASTIC FLOW , 2002 .

[3]  Jie Shen.  Remarks on the pressure error estimates for the projection methods , 1994 .

[4]  Christophe Ancey,et al.  Plasticity and geophysical flows: A review , 2007 .

[5]  W. Hager.  Review: R. Glowinski, J. L. Lions and R. Trémolières, Numerical analysis of variational inequalities , 1983 .

[6]  Simone Scacchi,et al.  Parallel multilevel solvers for the cardiac electro-mechanical coupling , 2015 .

[7]  T. Papanastasiou.  Flows of Materials with Yield , 1987 .

[8]  Georgios C. Georgiou,et al.  Cessation of Couette and Poiseuille flows of a Bingham plastic and finite stopping times , 2005 .

[9]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[10]  Jie Shen.  On error estimates of projection methods for Navier-Stokes equations: first-order schemes , 1992 .

[11]  Andreas G. Boudouvis,et al.  Flows of viscoplastic materials: Models and computations , 1997 .

[12]  Hantaek Bae.  Navier-Stokes equations , 1992 .

[13]  Georgios C. Georgiou,et al.  Performance of the finite volume method in solving regularised Bingham flows: inertia effects in the lid-driven cavity flow , 2014, ArXiv.

[14]  Michel Bercovier,et al.  A finite-element method for incompressible non-Newtonian flows , 1980 .

[15]  R. Temam,et al.  Navier-Stokes equations: theory and numerical analysis: R. Teman North-Holland, Amsterdam and New York. 1977. 454 pp. US $45.00 , 1978 .

[16]  D. Vola,et al.  Laminar unsteady flows of Bingham fluids: a numerical strategy and some benchmark results , 2003 .

[17]  R. Byron Bird,et al.  The Rheology and Flow of Viscoplastic Materials , 1983 .

[18]  Evan Mitsoulis,et al.  Flow of Bingham plastics in a lid-driven square cavity , 2001 .

[19]  O. Botella,et al.  BENCHMARK SPECTRAL RESULTS ON THE LID-DRIVEN CAVITY FLOW , 1998 .

[20]  Roland Glowinski,et al.  On the numerical simulation of Bingham visco-plastic flow: Old and new results , 2007 .

[22]  F. J. Sánchez.  Application of a First-Order Operator Splitting Method to Bingham Fluid Flow Simulation , 1998 .

[23]  Jie Shen,et al.  An overview of projection methods for incompressible flows , 2006 .

[24]  L. Muravleva.  Uzawa-like methods for numerical modeling of unsteady viscoplastic Bingham medium flows , 2015 .

[25]  J. Lions,et al.  Les inéquations en mécanique et en physique , 1973 .

[26]  M. Olshanskii.  Analysis of semi-staggered finite-difference method with application to Bingham flows , 2009 .

[27]  M. Fortin,et al.  Augmented Lagrangian methods : applications to the numerical solution of boundary-value problems , 1983 .

[28]  R. Glowinski,et al.  Qualitative properties and approximation of solutions of Bingham flows: On the stabilization for large time and the geometry of the support , 2010 .

[29]  Yongmin Zhang,et al.  Error estimates for the numerical approximation of time-dependent flow of Bingham fluid in cylindrical pipes by the regularization method , 2003, Numerische Mathematik.

[30]  M. Olshanskii,et al.  Two finite-difference schemes for calculation of Bingham fluid flows in a cavity , 2008 .

[31]  Pierre Saramito,et al.  An adaptive finite element method for Bingham fluid flows around a cylinder , 2003 .