Cellular automaton simulations of cement hydration and microstructure development

Cellular automaton algorithms, which operate on a starting digital image of a water-cement suspension, are described. The algorithms simulate the microstructure development process due to hydration reactions that occurs between cement and water. This paper describes the evolution of the cement model from a simple model, which treated the cement particles as single-phase materials, with a greatly simplified hydration chemistry, into a model which has many more chemical species and includes numerous reactions which eventually convert the viscous water-cement suspension into a rigid porous solid. Methods are presented for generating two- and three-dimensional images representing suspension initial conditions; these are derived from both micrographs of real cements and computer-based algorithms. The 2D initial images are based on the processing of backscattered electron and x-ray images of real cement suspensions. The 3D images employ either spheres to represent cement particles, or more realistic randomly shaped particles via an algorithm which smooths and thresholds a 3D lattice whose sites are initially populated with random white noise. A convenient measure of the point at which the initial paste turns into a solid material is the percolation threshold of the solids. Consideration of these models has already led to the prediction and subsequent experimental observation of a sharply defined onset of shear wave propagation, from ultrasonic measurements through hydrating cement slurries. The amount of hydration needed to reach the percolation threshold can be determined in the present simulations, and our results are compared with time of shear wave onset in actual cement slurries. Variants of the basic model provide insight into both early-time behaviour that is of primary interest to oil well cementing and the later-time microstructural properties that are of interest in the construction industry.

[1]  W. Daniel Hillis,et al.  The connection machine: A computer architecture based on cellular automata , 1984 .

[2]  E. Garboczi,et al.  Computer simulation of the diffusivity of cement-based materials , 1992 .

[3]  P. Coveney,et al.  Simple chemical clock reactions: application to cement hydration , 1993 .

[4]  D. R. Janecky,et al.  A lattice gas automata model for heterogeneous chemical reactions at mineral surfaces and in pore networks , 1991 .

[5]  K. G. Eggert,et al.  Lattice gas automata for flow through porous media , 1991 .

[6]  Daniel H. Rothman,et al.  Cellular‐automaton fluids: A model for flow in porous media , 1988 .

[7]  Edward J. Garboczi,et al.  Percolation of phases in a three-dimensional cement paste microstructural model , 1991 .

[8]  E. Garboczi,et al.  The elastic moduli of a sheet containing circular holes , 1992 .

[9]  I. Odler,et al.  On the origin of Portland cement setting , 1992 .

[10]  P. Brown Early Hydration of Tetracalcium Aluminoferrite in Gypsum and Lime‐Gypsum Solutions , 1987 .

[11]  Bruce M. Boghosian,et al.  A cellular automata simulation of two-phase flow on the CM-2 Connection Machine computer , 1988, Proceedings Supercomputing Vol.II: Science and Applications.

[12]  E. Garboczi,et al.  Modelling the leaching of calcium hydroxide from cement paste: effects on pore space percolation and diffusivity , 1992 .

[13]  Jayanth R. Banavar,et al.  Image‐based models of porous media: Application to Vycor glass and carbonate rocks , 1991 .

[14]  Colin M. Sayers,et al.  Ultrasonic propagation through hydrating cements , 1993 .

[15]  Edward J. Garboczi,et al.  Digital simulation of the aggregate–cement paste interfacial zone in concrete , 1991 .

[16]  S. Wolfram Cellular automaton fluids 1: Basic theory , 1986 .

[17]  R. Bogue,et al.  Calculation of the Compounds in Portland Cement , 1929 .

[18]  J. F. Young,et al.  Phase Composition of Hydrated DSP Cement Pastes , 1993 .

[19]  Dab,et al.  Lattice-gas automata for coupled reaction-diffusion equations. , 1991, Physical review letters.

[20]  E. Garboczi,et al.  Cellular automaton algorithm for surface mass transport due to curvature gradients simulations of sintering , 1992 .