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