validated code

clean.sh

+rm output/grid* output/output.dat src/*.o src/*.mod src/pafiX
+rm -r results
+mkdir results
+rm *.e*

input/input.dat

+pafiX input data
+----------------
+dimensions in x,y,z-direction (m,m,m) =
+0.24,0.04,0.04 
+number of cells in x,y,z-direction (-,-,-) =
+60,60,60
+start time-step (-), time-steps to compute (-) =
+1,10000
+interval of writing result files, restart files (-,-) =
+10,10000
+CFL number (-) =
+0.3
+initial bulk flow (m/s), pressure gradient (N/m**3) =
+3.65,-5.8
+fluid density (kg/m**3), fluid kinematic viscosity (m**2/s) =
+1.2,1.46E-5
+particle number density (-/m**3), seet at nt (-) =
+1E8,10000
+particle radius (m), material density (kg/m**3) =
+50E-6,1200
+particle charge (C) =
+1.E-13
+gravity vector in x,y,z-dir. (m/s**2,m/s**2,m/s**2) =
+9.81,0,0

monitor.plt

+#set term postscript enhanced eps 10 color linewidth 1 'Times-Roman'
+#set term dumb
+#set output "post.eps"
+#set xrange [0:]
+set multiplot layout 3,2 title 'pafiX -- by H. Grosshans'
+f = 'output/monitor.dat'
+
+set key box opaque width 2 t l
+set xlabel offset 0,0.5
+set ylabel offset 2
+set linetype 1 ps 0.5 lc 3
+set linetype 2 ps 0.5 lc 4
+set xlabel 's'
+
+set ylabel 'm/s'
+p f u 2:4 t 'u_{cl}', f u 2:5 t 'u_{av}'
+
+set ylabel '-'
+p f u 2:8 t 'Re_c',   f u 2:9 t 'Re_b'
+
+set ylabel 'm/s'
+p f u 2:6 t 'v_{rms}', f u 2:7 t 'w_{rms}'
+
+set ylabel '-'
+p f u 2:($6/$11) t 'v_{rms}/u_{tau}', f u 2:($7/$11) t 'w_{rms}/u_{tau}'
+
+#set ylabel '1/s'
+#p f u 2:14 t 'gradu'
+
+set ylabel 'N/m**2'
+p f u 2:10 t 'tau_w'
+
+set ylabel '-'
+p f u 2:12 t 'Re_{tau}'
+
+pause 20
+reread

run.sh

+mpiexec -np 4 src/pafiX