Simulations of quasi-static foam flow through a diverging-converging channel
Simon Cox and I. Tudur Davies
Department of Mathematics, Aberystwyth University, Aberystwyth SY23 3BZ, UK
and Coleg Cymraeg Cenedlaethol
published as
S.J. Cox and I.T. Davies (2016) Simulations of quasi-static foam flow through a diverging-converging channel Korea-Aust. Rheol. J. 28: 181--186.
Preprint available here as wedge_flow_3mono.pdf.
=====================
Simulations were conducted in Surface Evolver, Version 2.71c.
See http://facstaff.susqu.edu/brakke/evolver/evolver.html
We saved a dump file at each time step and from this extracted the
information given below.
For each of the two simulations (monodisperse, polydisperse) there are four files containing
information required to reconstruct the fields of stress, texture, topological changes, and
velocity. To generate the figures in the paper we analysed this data using bespoke C programs
and plotted the result using gnuplot, Version 4.4 (http://www.gnuplot.info/).
The format of these files, with XXXX replaced by mono or poly, is:
stress_725_XXXX_1.dat
The stress is an integral of the surface tension over all edges, i.e. the orientation of the
tangent. The files contain one block of data for each time step, separated by two blank lines. Each
block is prefaced with a comment #YYYY giving the number of lines in each block. Each line
corresponds to one edge of the foam and has six pieces of information: (edge id, x coordinate of
midpoint, y coordinate of midpoint, xx component of stress, xy component of stress, yy component of
stress).
t1data_725_XXXX_1.dat
Three columns containing (time step, x position of T1, y position of T1)
texture_725_XXXX_1.dat
To calculate strain, we complement the velocity data below with this file containing, for each edge,
the id of the bubble on each side of it. The files contain one block of data for each time step,
separated by two blank lines. Each block is prefaced with a comment #YYYY giving the number of lines
in each block. Each line corresponds to one edge of the foam and gives the two bubble ids.
vels_725_XXXX_1.dat
The files contain one block of data for each time step, separated by two blank lines. Each block is
prefaced with a comment #YYYY giving the number of lines (corresponding to the number of bubbles,
always 725) in each block. Each line corresponds to one bubble and gives: (x coordinate of bubble
centre, x coordinate of bubble centre, bubble pressure).
======================
SJC
September 2016