diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/Ksl b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/Ksl
index da68e074..097284a2 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/Ksl
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/Ksl
@@ -24,27 +24,27 @@ boundaryField
     {
         type            cyclic;
     }
-    
+
     top
     {
         type            cyclic;
     }
-    
+
     lr1
     {
         type            cyclic;
     }
-    
+
     lr2
     {
         type            cyclic;
     }
-    
+
     vh1
     {
         type            cyclic;
     }
-    
+
     vh2
     {
         type            cyclic;
diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/T b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/T
index 9783870d..a39c499b 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/T
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/T
@@ -1,9 +1,9 @@
 /*--------------------------------*- C++ -*----------------------------------*\
-| =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.6                                   |
-|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
-|    \\/     M anipulation  |                                                 |
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  6
+     \\/     M anipulation  |
 \*---------------------------------------------------------------------------*/
 FoamFile
 {
@@ -24,27 +24,27 @@ boundaryField
     {
         type            cyclic;
     }
-    
+
     top
     {
         type            cyclic;
     }
-    
+
     lr1
     {
         type            cyclic;
     }
-    
+
     lr2
     {
         type            cyclic;
     }
-    
+
     vh1
     {
         type            cyclic;
     }
-    
+
     vh2
     {
         type            cyclic;
diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/U b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/U
index 885d559d..7e80556e 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/U
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/U
@@ -24,27 +24,27 @@ boundaryField
     {
         type            cyclic;
     }
-    
+
     top
     {
         type            cyclic;
     }
-    
+
     lr1
     {
         type            cyclic;
     }
-    
+
     lr2
     {
         type            cyclic;
     }
-    
+
     vh1
     {
         type            cyclic;
     }
-    
+
     vh2
     {
         type            cyclic;
diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/Us b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/Us
index 584bbc9a..39613c01 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/Us
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/Us
@@ -24,27 +24,27 @@ boundaryField
     {
         type            cyclic;
     }
-    
+
     top
     {
         type            cyclic;
     }
-    
+
     lr1
     {
         type            cyclic;
     }
-    
+
     lr2
     {
         type            cyclic;
     }
-    
+
     vh1
     {
         type            cyclic;
     }
-    
+
     vh2
     {
         type            cyclic;
diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/p b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/p
index 6ab594d3..e11f415d 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/p
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/p
@@ -24,27 +24,27 @@ boundaryField
     {
         type            cyclic;
     }
-    
+
     top
     {
         type            cyclic;
     }
-    
+
     lr1
     {
         type            cyclic;
     }
-    
+
     lr2
     {
         type            cyclic;
     }
-    
+
     vh1
     {
         type            cyclic;
     }
-    
+
     vh2
     {
         type            cyclic;
diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/rho b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/rho
index 466c36c4..45681188 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/rho
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/rho
@@ -24,27 +24,27 @@ boundaryField
     {
         type            cyclic;
     }
-    
+
     top
     {
         type            cyclic;
     }
-    
+
     lr1
     {
         type            cyclic;
     }
-    
+
     lr2
     {
         type            cyclic;
     }
-    
+
     vh1
     {
         type            cyclic;
     }
-    
+
     vh2
     {
         type            cyclic;
diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/sSmoothField b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/sSmoothField
index 895f49f0..b63b83b6 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/sSmoothField
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/sSmoothField
@@ -1,9 +1,9 @@
 /*--------------------------------*- C++ -*----------------------------------*\
-| =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.6                                   |
-|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
-|    \\/     M anipulation  |                                                 |
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  6
+     \\/     M anipulation  |
 \*---------------------------------------------------------------------------*/
 FoamFile
 {
@@ -25,27 +25,27 @@ boundaryField
     {
         type            cyclic;
     }
-    
+
     top
     {
         type            cyclic;
     }
-    
+
     lr1
     {
         type            cyclic;
     }
-    
+
     lr2
     {
         type            cyclic;
     }
-    
+
     vh1
     {
         type            cyclic;
     }
-    
+
     vh2
     {
         type            cyclic;
diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/vSmoothField b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/vSmoothField
index 56fb9e07..1f1754fe 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/vSmoothField
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/vSmoothField
@@ -1,9 +1,9 @@
 /*--------------------------------*- C++ -*----------------------------------*\
-| =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.6                                   |
-|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
-|    \\/     M anipulation  |                                                 |
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  6
+     \\/     M anipulation  |
 \*---------------------------------------------------------------------------*/
 FoamFile
 {
@@ -25,27 +25,27 @@ boundaryField
     {
         type            cyclic;
     }
-    
+
     top
     {
         type            cyclic;
     }
-    
+
     lr1
     {
         type            cyclic;
     }
-    
+
     lr2
     {
         type            cyclic;
     }
-    
+
     vh1
     {
         type            cyclic;
     }
-    
+
     vh2
     {
         type            cyclic;
diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/voidfraction b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/voidfraction
index f797ffe9..3ce47938 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/voidfraction
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/0/voidfraction
@@ -24,27 +24,27 @@ boundaryField
     {
         type            cyclic;
     }
-    
+
     top
     {
         type            cyclic;
     }
-    
+
     lr1
     {
         type            cyclic;
     }
-    
+
     lr2
     {
         type            cyclic;
     }
-    
+
     vh1
     {
         type            cyclic;
     }
-    
+
     vh2
     {
         type            cyclic;
diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/system/blockMeshDict b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/system/blockMeshDict
index 862008f6..a03b7919 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/CFD/system/blockMeshDict
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/CFD/system/blockMeshDict
@@ -1,9 +1,9 @@
 /*--------------------------------*- C++ -*----------------------------------*\
-| =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  dev                                   |
-|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
-|    \\/     M anipulation  |                                                 |
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  6
+     \\/     M anipulation  |
 \*---------------------------------------------------------------------------*/
 FoamFile
 {
@@ -16,7 +16,7 @@ FoamFile
 
 convertToMeters 1.0;
 
-vertices 
+vertices
 (
     (-0.004 -0.004 -0.016)
     ( 0.004 -0.004 -0.016)
@@ -30,7 +30,7 @@ vertices
 
 blocks
 (
-    hex (0 1 2 3 4 5 6 7) (32 32 128) 
+    hex (0 1 2 3 4 5 6 7) (32 32 128)
     simpleGrading (1 1 1)
 );
 
@@ -42,38 +42,38 @@ boundary
 (
     bottom
     {
-	type cyclic;
-	neighbourPatch  top;
-        faces	((0 3 2 1));
+        type cyclic;
+        neighbourPatch  top;
+        faces ((0 3 2 1));
     }
     top
     {
-	type cyclic;
-	neighbourPatch  bottom;
-        faces	((4 5 6 7));
+        type cyclic;
+        neighbourPatch  bottom;
+        faces ((4 5 6 7));
     }
     lr1
     {
-	type cyclic;
-	neighbourPatch  lr2;
-        faces	((0 3 7 4));
+        type cyclic;
+        neighbourPatch  lr2;
+        faces ((0 3 7 4));
     }
     lr2
     {
-	type cyclic;
-	neighbourPatch  lr1;
-	faces ((1 5 6 2));
+        type cyclic;
+        neighbourPatch  lr1;
+        faces ((1 5 6 2));
     }
     vh1
     {
-	type cyclic;
-	neighbourPatch  vh2;
+        type cyclic;
+        neighbourPatch  vh2;
         faces ((0 1 4 5));
     }
     vh2
     {
-	type cyclic;
-	neighbourPatch  vh1;
+        type cyclic;
+        neighbourPatch  vh1;
         faces ((2 6 7 3));
     }
 );
diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/DEM/in.liggghts_init b/tutorials/cfdemSolverPimple/PeriodicBox/DEM/in.liggghts_init
index 2dce10d5..f3c7ef4f 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/DEM/in.liggghts_init
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/DEM/in.liggghts_init
@@ -1,10 +1,10 @@
 ################################################################################
-# brief: Periodic box - creating particle bed                        		#
+# brief: Periodic box - creating particle bed                                  #
 #                                                                              #
 # Periodic box simulations of Geldart A particles (d=0.075mm, rho=1500kg/m3)   #
-# Cf. Radl, S. and Sundaresan, S., A drag model for filtered Euler–Lagrange 	#
-#     simulations of clustered gas–particle suspensions. 			#
-#     Chemical engineering science, (2014).  					#
+# Cf. Radl, S. and Sundaresan, S., A drag model for filtered Euler–Lagrange    #
+#     simulations of clustered gas–particle suspensions.                       #
+#     Chem. Eng. Sci., (2014)                                                  #
 # mass loading = 0.05                                                          #
 # authors: Behrad Esgandari                                                    #
 # date: Aug 2023                                                               #
@@ -50,10 +50,10 @@ neigh_modify    delay 0
 
 # define the material properties required for granular pair styles
 # type 1 = particles with friction
-fix m1 all property/global youngsModulus peratomtype 1e6 
-fix m2 all property/global poissonsRatio peratomtype 0.42 
-fix m3 all property/global coefficientRestitution peratomtypepair 1 0.9  
-fix m4 all property/global coefficientFriction peratomtypepair 1 0.1 
+fix m1 all property/global youngsModulus peratomtype 1e6
+fix m2 all property/global poissonsRatio peratomtype 0.42
+fix m3 all property/global coefficientRestitution peratomtypepair 1 0.9
+fix m4 all property/global coefficientFriction peratomtypepair 1 0.1
 
 # specify contact model to use
 pair_style  gran model hertz tangential incremental_history limitForce on
diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/DEM/in.liggghts_run b/tutorials/cfdemSolverPimple/PeriodicBox/DEM/in.liggghts_run
index a7a9c06d..8bdc0ba8 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/DEM/in.liggghts_run
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/DEM/in.liggghts_run
@@ -1,10 +1,10 @@
 ################################################################################
-# brief: Periodic box - run                       				#
+# brief: Periodic box - run                                                    #
 #                                                                              #
 # Periodic box simulations of Geldart A particles (d=0.075mm, rho=1500kg/m3)   #
-# Cf. Radl, S. and Sundaresan, S., A drag model for filtered Euler–Lagrange 	#
-#     simulations of clustered gas–particle suspensions. 			#
-#     Chemical engineering science, (2014).  					#
+# Cf. Radl, S. and Sundaresan, S., A drag model for filtered Euler–Lagrange    #
+#     simulations of clustered gas–particle suspensions.                       #
+#     Chem. Eng. Sci., (2014)                                                  #
 # mass loading = 0.05                                                          #
 # authors: Behrad Esgandari                                                    #
 # date: Aug 2023                                                               #
@@ -48,10 +48,10 @@ neigh_modify every 5 delay 0 check no
 
 # define the material properties required for granular pair styles
 # type 1 = particles with friction
-fix m1 all property/global youngsModulus peratomtype 1e6 
-fix m2 all property/global poissonsRatio peratomtype 0.42 
-fix m3 all property/global coefficientRestitution peratomtypepair 1 0.9  
-fix m4 all property/global coefficientFriction peratomtypepair 1 0.1 
+fix m1 all property/global youngsModulus peratomtype 1e6
+fix m2 all property/global poissonsRatio peratomtype 0.42
+fix m3 all property/global coefficientRestitution peratomtypepair 1 0.9
+fix m4 all property/global coefficientFriction peratomtypepair 1 0.1
 
 # specify contact model to use
 pair_style  gran model hertz tangential incremental_history limitForce on
@@ -84,7 +84,7 @@ compute_modify  thermo_temp dynamic yes
 # output particle properties (same frequency as CFD output)
 dump dmp all custom/vtk 10000 ../DEM/post/liggghts_run_*.vtk id type x y z vx vy vz &
      fx fy fz omegax omegay omegaz radius
-     
+
 dump dmp2 all custom 10000 ../DEM/post/liggghts_run_*.run id type x y z vx vy vz &
      fx fy fz omegax omegay omegaz radius
 
diff --git a/tutorials/cfdemSolverPimple/PeriodicBox/parDEMrun.sh b/tutorials/cfdemSolverPimple/PeriodicBox/parDEMrun.sh
index abbd57a6..a9448871 100755
--- a/tutorials/cfdemSolverPimple/PeriodicBox/parDEMrun.sh
+++ b/tutorials/cfdemSolverPimple/PeriodicBox/parDEMrun.sh
@@ -1,10 +1,9 @@
 #!/bin/bash
-
-#===================================================================#
+#------------------------------------------------------------------------------
 # DEM run script for periodic box testcase
 # init periodic box
 # Behrad Esgandari - August 2023
-#===================================================================#
+#------------------------------------------------------------------------------
 
 #- source CFDEM env vars
 . ~/.bashrc
@@ -13,7 +12,7 @@
 source $CFDEM_PROJECT_DIR/etc/functions.sh
 
 echo "starting DEM run in parallel..."
-#--------------------------------------------------------------------------------#
+#------------------------------------------------------------------------------
 #- define variables
 casePath="$(dirname "$(readlink -f ${BASH_SOURCE[0]})")"
 logpath="$casePath"
@@ -23,7 +22,7 @@ solverName="in.liggghts_init"
 nrProcs=32
 machineFileName="none"
 debugMode="off"
-#--------------------------------------------------------------------------------#
+#------------------------------------------------------------------------------
 
 #- call function to run DEM case
 parDEMrun $logpath $logfileName $casePath $headerText $solverName $nrProcs $machineFileName $debugMode