Numerical pressure threshold method to simulate cement paste slump flow

In this study, we propose a novel pressure threshold method for the Bingham numerical model and apply the method to simulate cement paste slump flow. Calculation is divided into three steps. First, the cement paste is computed as a Newtonian fluid. An advection step and a non-advection step are computed in sequence. In the non-advection step, the fractional steps are implemented to calculate viscous, gravity, and pressure terms. Second, a pressure threshold judgment is used to modify the motion state of cement paste cells. Third, an iterative correction, which is an iterative process to determine the precise unyielded (rigid) and yielded (fluid) regions, is conducted for the Bingham computational model. The volume/surface integrated average-based multi-moment method scheme is used to compose a new finite volume formula for solving general fluid dynamic problems. The tangent of hyperbola for interface capturing approach is used to capture free boundaries in multi-fluid simulations. The computational modeling enables accurate simulation of the flow process of cement paste to analyze different mix designs and to evaluate the workability of cement paste. The proposed methodology is relevant to computational mechanics with applications in cement flow simulations.

[1]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[2]  Takashi Yabe,et al.  Cubic interpolated pseudo-particle method (CIP) for solving hyperbolic-type equations , 1985 .

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

[4]  Jingwei Zheng,et al.  GPU-based parallel algorithm for particle contact detection and its application in self-compacting concrete flow simulations , 2012 .

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

[6]  C. Ferraris,et al.  Workability of Self-Compacting Concrete , 2000 .

[7]  Feng Xiao,et al.  Unified formulation for compressible and incompressible flows by using multi-integrated moments I: one-dimensional inviscid compressible flow , 2004 .

[8]  Nicolas Roussel,et al.  “Fifty-cent rheometer” for yield stress measurements: From slump to spreading flow , 2005 .

[9]  Hans Muhlhaus,et al.  A Lagrangian integration point finite element method for large deformation modeling of viscoelastic geomaterials , 2003 .

[10]  R. Huilgol Variational inequalities in the flows of yield stress fluids including inertia: Theory and applications , 2002 .

[11]  I. Frigaard,et al.  On the usage of viscosity regularisation methods for visco-plastic fluid flow computation , 2005 .

[12]  Nicolas Roussel,et al.  Steady and transient flow behaviour of fresh cement pastes , 2005 .

[13]  Johan Silfwerbrand,et al.  Numerical simulation of fresh SCC flow: applications , 2011 .

[14]  Evan Mitsoulis,et al.  FLOWS OF VISCOPLASTIC MATERIALS: MODELS AND COMPUTATIONS , 2007 .

[15]  Hajime Okamura,et al.  Self-Compacting Concrete , 2000 .

[16]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

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

[18]  Georgios C. Georgiou,et al.  Solution of the square lid-driven cavity flow of a Bingham plastic using the finite volume method , 2013, ArXiv.

[19]  Frédéric Dufour,et al.  Computational modeling of concrete flow: General overview , 2007 .

[20]  A. Clappier A Correction Method for Use in Multidimensional Time-Splitting Advection Algorithms: Application to Two- and Three-Dimensional Transport , 1998 .

[21]  Atsushi Yashima,et al.  Estimating the impact force generated by granular flow on a rigid obstruction , 2009 .

[22]  D. Juric,et al.  A front-tracking method for the computations of multiphase flow , 2001 .

[23]  Rachid Bennacer,et al.  The mini-conical slump flow test: Analysis and numerical study , 2010 .

[24]  Bhushan Lal Karihaloo,et al.  Modelling the flow of self‐compacting concrete , 2011 .

[25]  Howard L. Schreyer,et al.  Axisymmetric form of the material point method with applications to upsetting and Taylor impact problems , 1996 .

[26]  Amir M. Halabian,et al.  Simulation of concrete flow in V-funnel test and the proper range of viscosity and yield stress for SCC , 2014 .

[27]  Knut Krenzer,et al.  Simulation of fresh concrete flow using Discrete Element Method (DEM): theory and applications , 2014 .

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

[29]  Feng Xiao,et al.  Unified formulation for compressible and incompressible flows by using multi-integrated moments II: Multi-dimensional version for compressible and incompressible flows , 2006, J. Comput. Phys..

[30]  Takashi Yabe,et al.  A universal solver for hyperbolic equations by cubic-polynomial interpolation I. One-dimensional solver , 1991 .

[31]  T. Yabe,et al.  An Exactly Conservative Semi-Lagrangian Scheme (CIP–CSL) in One Dimension , 2001 .

[32]  R. Glowinski,et al.  Numerical Methods for Nonlinear Variational Problems , 1985 .

[33]  A computational model for incompressible flow including surface tension , 1999 .

[34]  Feng Xiao,et al.  A simple algebraic interface capturing scheme using hyperbolic tangent function , 2005 .