Added square tube mixed wet case

example/SquareTube/MixedWet.sh

@@ -0,0 +1,8 @@
+#!/bin/bash
+
+LBPM_DIR=../../tests
+
+python SquareTube.py
+mpirun -np 1 $LBPM_DIR/lbpm_serial_decomp input.db
+mpirun -np 4 $LBPM_DIR/lbpm_segmented_pp input.db
+mpirun -np 4 $LBPM_DIR/lbpm_color_simulator input.db

example/SquareTube/SquareTube.py

@@ -0,0 +1,43 @@
+import numpy
+
+nx=96
+ny=24
+nz=24
+N=nx*ny*nz
+
+mesh=(nx,ny,nz)
+data=numpy.ones(mesh,dtype=numpy.int8)
+
+LabelTop=-1
+LabelBottom=-2
+LabelLeft=-3
+LabelRight=-4
+
+#print(data)
+print("Writing parllel plates")
+print("Mesh size: "+repr(mesh))
+print("Top label: "+repr(LabelTop))
+print("Bottom label: "+repr(LabelBottom))
+print("Left label: "+repr(LabelLeft))
+print("Right label: "+repr(LabelRight))
+
+# assign a bubble in the middle
+for x in range(12,72):
+  for y in range(0,ny):
+    for z in range(0,nz):
+       data[x,y,z]=2
+
+# solid walls with different component labels
+for x in range(0,nx):
+  for y in range(0,ny):
+    data[x,y,0]=LabelBottom
+    data[x,y,nz-1]=LabelTop
+
+
+for x in range(0,nx):
+  for z in range(0,nz):
+    data[x,0,z]=LabelLeft
+    data[x,ny-1,z]=LabelRight
+
+
+data.tofile("SquareTube.raw")

example/SquareTube/input.db

@@ -0,0 +1,45 @@
+Color {
+    tauA   = 0.7;  
+    tauB   = 0.7;  
+    rhoA   = 1.0;  
+    rhoB   = 1.0;  
+    alpha = 1e-3;
+    beta  = 0.95;
+    F = 0, 0, 0
+    Restart = false
+    pBC = 0
+    din = 1.0
+    dout = 1.0
+    timestepMax = 3000
+    interval = 1000
+    tol = 1e-5;
+    das = 0.1
+    dbs = 0.9
+    flux = 0.0
+    ComponentLabels = -1, -2, -3, -4
+    ComponentAffinity = -1.0, -0.5, -0.8, -0.2;
+}
+
+Domain {
+    Filename = "SquareTube.raw"
+    nproc = 1, 1, 4     // Number of processors (Npx,Npy,Npz)
+    n = 24, 24, 24      // Size of local domain (Nx,Ny,Nz)
+    N = 24, 24, 96      // size of the input image
+    n_spheres = 1       // Number of spheres
+    L = 1, 1, 1         // Length of domain (x,y,z)
+    BC = 0              // Boundary condition type
+    ReadValues = -2, -1, 1, 2
+    WriteValues = -2, -1, 1, 2
+}
+
+Analysis {
+    blobid_interval = 1000      // Frequency to perform blob identification
+    analysis_interval = 1000    // Frequency to perform analysis
+    restart_interval = 1000    // Frequency to write restart data
+    visualization_interval = 1000        // Frequency to write visualization data
+    restart_file = "Restart"    // Filename to use for restart file (will append rank)
+    N_threads    = 4            // Number of threads to use
+    load_balance = "independent" // Load balance method to use: "none", "default", "independent"
+}
+
+