2D MPI enhanced solver

EnhancedSolver/2D-seq/Makefile

@@ -70,6 +70,18 @@ format:
 	done
 	@echo "Done"
 
+vis:
+	$(info ===>  GENERATE VISUALIZATION)
+	@gnuplot -e "filename='pressure.dat'" ./surface.plot
+	@gnuplot -e "filename='velocity.dat'" ./vector.plot
+
+vis_clean:
+	$(info ===>  CLEAN VISUALIZATION)
+	@rm -f *.dat
+	@rm -f *.png
+	@rm -f ./vis_files/*.dat
+	@rm -f ./vis_files/*.gif
+
 $(BUILD_DIR):
 	@mkdir $(BUILD_DIR)
 

EnhancedSolver/2D-seq/backstep.par

@@ -15,7 +15,7 @@ bcRight    3			#
 gx    0.0      # Body forces (e.g. gravity)
 gy    0.0      #
 
-re    36500.0	   # Reynolds number
+re    20000.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
@@ -45,11 +45,18 @@ rho      0.52
 omg      1.8       # relaxation parameter for SOR iteration
 gamma    0.9       # upwind differencing factor gamma
 
+# Multigrid data:
+# ---------
+
+levels        3         # Multigrid levels
+presmooth     5         # Pre-smoothning iterations
+postsmooth    5         # Post-smoothning iterations
+
 # Particle Tracing Data:
 # -----------------------
 
 numberOfParticles   200
-startTime           100
+startTime           10
 injectTimePeriod    1.0
 writeTimePeriod     0.5
 

EnhancedSolver/2D-seq/canal.par

@@ -44,7 +44,13 @@ eps      0.0001   # stopping tolerance for pressure iteration
 rho      0.52
 omg      1.8       # relaxation parameter for SOR iteration
 gamma    0.9       # upwind differencing factor gamma
-levels   5         # Multigrid levels
+
+# Multigrid data:
+# ---------
+
+levels        3         # Multigrid levels
+presmooth     5         # Pre-smoothning iterations
+postsmooth    5         # Post-smoothning iterations
 
 # Particle Tracing Data:
 # -----------------------

EnhancedSolver/2D-seq/config.mk

@@ -1,5 +1,5 @@
 # Supported: GCC, CLANG, ICC
-TAG ?= CLANG
+TAG ?= ICC
 ENABLE_OPENMP ?= false
 # Supported: sor, rb, mg
 SOLVER ?= mg

EnhancedSolver/2D-seq/dcavity.par

@@ -44,7 +44,13 @@ eps      0.001		# stopping tolerance for pressure iteration
 rho      0.5
 omg      1.8		# relaxation parameter for SOR iteration
 gamma    0.9		# upwind differencing factor gamma
-levels   5         # Multigrid levels
+
+# Multigrid data:
+# ---------
+
+levels        3         # Multigrid levels
+presmooth     10        # Pre-smoothning iterations
+postsmooth    5         # Post-smoothning iterations
 
 # Particle Tracing Data:
 # -----------------------

EnhancedSolver/2D-seq/karman.par

@@ -44,13 +44,19 @@ eps      0.001   # stopping tolerance for pressure iteration
 rho      0.52
 omg      1.75      # relaxation parameter for SOR iteration
 gamma    0.9       # upwind differencing factor gamma
-levels   5         # Multigrid levels
+
+# Multigrid data:
+# ---------
+
+levels        3         # Multigrid levels
+presmooth     5         # Pre-smoothning iterations
+postsmooth    5         # Post-smoothning iterations
 
 # Particle Tracing Data:
 # -----------------------
 
 numberOfParticles   200
-startTime           201
+startTime           50
 injectTimePeriod    1.0
 writeTimePeriod     0.5
 

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

@@ -98,7 +98,7 @@ static void advanceParticles(
             if (!gridIsFluid(p->grid, newI, newJ)) {
                 p->particlePool[i].flag = false;
                 p->removedParticles++;
-                printf("Forbidden movement of particle into obstacle!\n");
+                // printf("Forbidden movement of particle into obstacle!\n");
             }
         }
     }