Bubble Benchmark


Reference data

The benchmark data is provided in the following ascii text files. File names use the notation c#g#l# with c# indicating the test case, g# the group and l# the refinement level:

c1Test Case 1
c2Test Case 2
g1TU Dortmund (TP2D)
g2EPFL Lausanne (FreeLIFE)
g3Uni Magdebug (MooNMD)
l#Grid Refinement Level (higher number means denser grid)

FileAction
All in one zip fileDownload
c1g1l4.txtDownload
c1g1l5.txtDownload
c1g1l6.txtDownload
c1g1l7.txtDownload
c1g2l1.txtDownload
c1g2l2.txtDownload
c1g2l3.txtDownload
c1g3l1.txtDownload
c1g3l2.txtDownload
c1g3l3.txtDownload
c1g3l4.txtDownload
c2g1l4.txtDownload
c2g1l5.txtDownload
c2g1l6.txtDownload
c2g1l7.txtDownload
c2g1l8.txtDownload
c2g2l1.txtDownload
c2g2l2.txtDownload
c2g2l3.txtDownload
c2g3l2.txtDownload
c2g3l3.txtDownload
c2g3l4.txtDownload

The data is given in ascii format with columns corresponding to the following quantities.

Benchmark file format
Column Quantity
1 Time
2 Bubble mass or area
3 Circularity
4 Center of mass (y-coordinate)
5 Rise velocity

The benchmark curves can the be plotted with standard postprocessing tools like Matlab or Gnuplot:

Gnuplot code

plot "c1g1l4.txt" using 1:2 with \
lines linestyle 1 title 'Bubble mass/area'; set xlabel 'Time'
plot "c1g1l4.txt" using 1:3 with \
lines linestyle 1 title 'Circularity'; set xlabel 'Time'
plot "c1g1l4.txt" using 1:4 with \
lines linestyle 1 title 'Center of mass'; set xlabel 'Time'
plot "c1g1l4.txt" using 1:5 with \
lines linestyle 1 title 'Rise velocity'; set xlabel 'Time'

Matlab code

load c1g1l4.txt, plot(c1g1l4(:,1),c1g1l4(:,2)), ...
title('Bubble mass/area'), xlabel('Time')
load c1g1l4.txt, plot(c1g1l4(:,1),c1g1l4(:,3)), ...
title('Circularity'), xlabel('Time')
load c1g1l4.txt, plot(c1g1l4(:,1),c1g1l4(:,4)), ...
title('Center of mass'), xlabel('Time')
load c1g1l4.txt, plot(c1g1l4(:,1),c1g1l4(:,5)), ...
title('Rise velocity'), xlabel('Time')

The data for the bubble shapes can be found in the files that end with a s character, for example c1g1l4s.txt. The shapes are stored as line segments with x-coordinates in the first column and y-coordinates in the second in consecutive line pairs (that is the first line end point is stored at odd numbered lines and the other end points on the following lines). The shapes can for example be plotted in Gnuplot or Matlab with the following code:

Gnuplot code

plot "c1g1l4s.txt" with lines linestyle 1; \
unset key; set size ratio -1

Matlab code

load c1g1l4s.txt
data = c1g1l4s;
hold on
nSegments = size(data,1)/2;
for i=1:nSegments
  plot(data(2*(i-1)+1:2*i,1),data(2*(i-1)+1:2*i,2))
end
axis equal

Downloads

FileAction
All in one zip fileDownload
c1g1l4s.txtDownload
c1g1l5s.txtDownload
c1g1l6s.txtDownload
c1g1l7s.txtDownload
c1g2l1s.txtDownload
c1g2l2s.txtDownload
c1g2l3s.txtDownload
c1g3l4s.txtDownload
c2g1l4s.txtDownload
c2g1l5s.txtDownload
c2g1l6s.txtDownload
c2g1l7s.txtDownload
c2g1l8s.txtDownload
c2g2l1s.txtDownload
c2g2l2s.txtDownload
c2g2l3s.txtDownload
c2g3l4s.txtDownload


Funding

Funded by:
  • German Research foundation (DFG); Grant Number: To143/9
  • Paketantrag PAK178; Grant Numbers: Tu102/27-1, Ku1530/5-1
  • Sonderforschungsbereich SFB708; Grant Number: TP B7
  • SFB TR R30; Grant Number: TPC3
  • Swiss National Science Foundation; Grant Number: 112166

Back to top