Working karman vortex street 2d seq

BasicSolver/2D-seq/karman.par

@@ -19,15 +19,15 @@ re    5050.0	   # Reynolds number
 
 u_init    1.0      # initial value for velocity in x-direction
 v_init    0.0      # initial value for velocity in y-direction
-p_init    1.0      # initial value for pressure
+p_init    0.0      # initial value for pressure
 
 # Geometry Data:
 # -------------
 
 xlength    30.0     # domain size in x-direction
-ylength    4.0	   # domain size in y-direction
+ylength    8.0	   # domain size in y-direction
-imax       200      # number of interior cells in x-direction
+imax       400      # number of interior cells in x-direction
-jmax       50	   # number of interior cells in y-direction
+jmax       200	   # number of interior cells in y-direction
 
 # Time Data:
 # ---------
@@ -54,9 +54,9 @@ injectTimePeriod    1.0
 writeTimePeriod     0.5
 
 x1                  0.0
-y1                  1.8
+y1                  3.8
 x2                  0.0
-y2                  2.1
+y2                  4.1
 
 # Obstacle Geometry Data:
 # -----------------------
@@ -64,7 +64,7 @@ y2                  2.1
 
 shape               2 
 xCenter             5.0
-yCenter             2.0
+yCenter             4.0
 xRectLength         2.0
 yRectLength         1.0
 circleRadius        1.0

BasicSolver/2D-seq/src/particletracing.c

@@ -103,7 +103,7 @@ void advanceParticles(ParticleTracer* particletracer, double* restrict u, double
                 particletracer->particlePool[i].flag = false;
             }
             int i_new = new_x/dx, j_new = new_y/dy;
-            if(S(i_new, j_new) == LOCAL)
+            if(S(i_new, j_new) != NONE)
             {
                 particletracer->particlePool[i].flag = false;
             }

BasicSolver/2D-seq/vis_files/karman_animate.plot

@@ -2,9 +2,9 @@ unset border; unset tics; unset key;
 set term gif animate delay 10
 set output "trace.gif"
 set xrange [0:30]
-set yrange [0:4]
+set yrange [0:8]
 set size ratio -1
-set object 1 circle front at 5.0,2.0 size 1.0 fillcolor rgb "black" lw 2
+set object 1 circle front at 5.0,4.0 size 1.0 fillcolor rgb "black" lw 2
 
 
 do for [ts=0:500] {