Merge pull request #126 from ramin1728/main

binarySystemOfParticles updated V1.0

tutorials/sphereGranFlow/binarySystemOfParticles/caseSetup/interaction

@@ -6,16 +6,16 @@ objectName  interaction;
 objectType  dicrionary;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
-materials     (prop1);          // properties of material
+materials                     (prop1);           // properties of material
 
-densities     (1000.0);         // density of materials [kg/m3]
+densities                    (1000.0);           // density of materials [kg/m3]
 
 contactListType     sortedContactList;
 
 contactSearch
 {
   
-   method          NBS;          // method for broad search 
+   method          NBS;                          // method for broad search 
     
    updateInterval   10;
 
@@ -34,18 +34,18 @@ model
 
    // Property (solid-solid Properties)
 
-   Yeff     (1.0e6);          // Young modulus [Pa]
+   Yeff     (1.0e6);                             // Young modulus [Pa]
 
-   Geff     (0.8e6);          // Shear modulus [Pa]
+   Geff     (0.8e6);                             // Shear modulus [Pa]
 
-   nu       (0.25);           // Poisson's ratio [-]
+   nu       (0.25);                              // Poisson's ratio [-]
 
-   en       (0.7);            // coefficient of normal restitution
+   en       (0.7);                               // coefficient of normal restitution
 
-   et       (1.0);            // coefficient of tangential restitution 
+   et       (1.0);                               // coefficient of tangential restitution 
 
-   mu       (0.3);            // dynamic friction 
+   mu       (0.3);                               // dynamic friction 
 
-   mur      (0.1);            // rolling friction     
+   mur      (0.1);                               // rolling friction     
 }
 

tutorials/sphereGranFlow/binarySystemOfParticles/caseSetup/particleInsertion

@@ -6,6 +6,6 @@ objectName 	particleInsertion;
 objectType 	dicrionary;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
-collisionCheck     No;          // is checked          -> yes or no
+collisionCheck     No;                           // is checked          -> yes or no
 
-active             No;          // is insertion active -> yes or no
+active             No;                           // is insertion active -> yes or no

tutorials/sphereGranFlow/binarySystemOfParticles/caseSetup/shapes

@@ -8,6 +8,6 @@ fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
 names 	  (smallSphere largeSphere); 	         // names of shapes 
 
-diameters (0.003 0.005);	                     // diameter of shapes 
+diameters             (0.003 0.005);	         // diameter of shapes 
 
-materials (prop1 prop1);	                     // material names for shapes 
+materials             (prop1 prop1);	         // material names for shapes 

tutorials/sphereGranFlow/binarySystemOfParticles/settings/domainDict

@@ -6,7 +6,7 @@ objectName 	domainDict;
 objectType 	dictionary;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
-globalBox          // Simulation domain: every particles that goes outside this domain will be deleted
+globalBox                                        // Simulation domain: every particles that goes outside this domain will be deleted
 {
 	min     (-0.12 -0.12 0);
 
@@ -21,43 +21,48 @@ decomposition
 boundaries
 {
 	// Determines how often (how many iterations) do you want to 
+
 	// rebuild the list of particles in the neighbor list 
+
 	// of all boundaries in the simulation domain 
+
 	neighborListUpdateInterval     50; 
 
 	// Determines how often do you want to update the new changes in the boundary
+
 	updateInterval                 10;
 
 	// The distance from the boundary plane within which particles are marked to be in the boundary list 
+
 	neighborLength              0.004;
 
 	left
 	{
-		type     exit;	        // other options: periodict, reflective 
+		type     exit;	                         // other options: periodict, reflective 
 	}
 
 	right
 	{
-		type     exit;          // other options: periodict, reflective 
+		type     exit;                           // other options: periodict, reflective 
 	}
 
 	bottom
 	{
-		type     exit;          // other options: periodict, reflective 
+		type     exit;                           // other options: periodict, reflective 
 	}
 
 	top
 	{
-		type     exit;          // other options: periodict, reflective 
+		type     exit;                           // other options: periodict, reflective 
 	}
 
 	rear
 	{
-		type     exit;          // other options: periodict, reflective 
+		type     exit;                           // other options: periodict, reflective 
 	}
 
 	front
 	{
-		type     exit;          // other options: periodict, reflective 
+		type     exit;                           // other options: periodict, reflective 
 	}
 }

tutorials/sphereGranFlow/binarySystemOfParticles/settings/geometryDict

@@ -6,17 +6,17 @@ objectName 	geometryDict;
 objectType 	dictionary;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
-motionModel     rotatingAxis;          // motion model can be rotatingAxis or stationary or vibrating
+motionModel     rotatingAxis;                    // motion model can be rotatingAxis or stationary or vibrating
 
-rotatingAxisInfo          // information for rotatingAxis motion model 
+rotatingAxisInfo                                 // information for rotatingAxis motion model 
 {
 	rotAxis 
 	{
-		p1 (0.0 0.0 0.0);	       // first point for the axis of rotation 
+		p1 (0.0 0.0 0.0);	                     // first point for the axis of rotation 
 
-		p2 (0.0 0.0 1.0);	       // second point for the axis of rotation
+		p2 (0.0 0.0 1.0);	                     // second point for the axis of rotation
 		
-		omega      1.214; 		   // rotation speed (rad/s)
+		omega      1.214; 		                 // rotation speed (rad/s)
 	}
 }
 
@@ -24,21 +24,21 @@ surfaces
 {
 	cylinder
 	{
-		type     cylinderWall;          // other options: cuboidWall and planeWall
+		type     cylinderWall;                   // other options: cuboidWall and planeWall
 
-		p1 (0.0 0.0 0.0);          // begin point of cylinder axis
+		p1 (0.0 0.0 0.0);                        // begin point of cylinder axis
 
-		p2 (0.0 0.0 0.1);          // end point of cylinder axis
+		p2 (0.0 0.0 0.1);                        // end point of cylinder axis
 
-		radius1     0.12;		   // radius at p1
+		radius1     0.12;		                 // radius at p1
 
-		radius2     0.12;		   // radius at p2
+		radius2     0.12;		                 // radius at p2
 
-		resolution    24;          // number of divisions
+		resolution    24;                        // number of divisions
 
-		material   prop1;          // material name of this wall
+		material   prop1;                        // material name of this wall
 
-		motion   rotAxis;		   // motion component name 
+		motion   rotAxis;		                 // motion component name 
 	}
 
 	/*
@@ -47,19 +47,19 @@ surfaces
 
 	wall1
 	{
-		type     planeWall;			 // other options: cuboidWall and  cylinderWall
+		type       planeWall;			         // other options: cuboidWall and  cylinderWall
 
-		p1 (-0.12 -0.12 0.0);	       // first point of the wall
+		p1 (-0.12 -0.12 0.0);	                 // first point of the wall
 
-		p2  (0.12 -0.12 0.0);          // second point of the wall
+		p2  (0.12 -0.12 0.0);                    // second point of the wall
 
-		p3   (0.12 0.12 0.0);          // third point of the wall
+		p3   (0.12 0.12 0.0);                    // third point of the wall
 
-		p4  (-0.12 0.12 0.0);          // fourth point of the wall
+		p4  (-0.12 0.12 0.0);                    // fourth point of the wall
 
-		material       prop1;          // material name of the wall 
+		material       prop1;                    // material name of the wall 
 
-		motion       rotAxis;		   // motion component name 
+		motion       rotAxis;		             // motion component name 
 	}
 
 	/*
@@ -68,18 +68,18 @@ surfaces
 
 	wall2
 	{
-		type     planeWall;          // other options: cuboidWall and  cylinderWall 
+		type       planeWall;                    // other options: cuboidWall and  cylinderWall 
 
-		p1 (-0.12 -0.12 0.1);          // first point of the wall
+		p1 (-0.12 -0.12 0.1);                    // first point of the wall
 
-		p2  (0.12 -0.12 0.1);          // second point of the wall
+		p2  (0.12 -0.12 0.1);                    // second point of the wall
 
-		p3   (0.12 0.12 0.1);          // third point of the wall
+		p3   (0.12 0.12 0.1);                    // third point of the wall
 
-		p4  (-0.12 0.12 0.1);          // fourth point of the wall
+		p4  (-0.12 0.12 0.1);                    // fourth point of the wall
 
-		material       prop1;          // material name of the wall 
+		material       prop1;                    // material name of the wall 
 
-		motion       rotAxis;          // motion component name 
+		motion       rotAxis;                    // motion component name 
 	}
 }

tutorials/sphereGranFlow/binarySystemOfParticles/settings/particlesDict

@@ -9,65 +9,66 @@ fileFormat  ASCII;
 setFields
 {
 	/*
-	    	Default value for fields defined for particles:
+	    Default value for fields defined for particles:
+
 		These fields should always be defined for simulations with spherical particles
 	*/
 
 	defaultValue 
 	{
-		velocity 		realx3 	(0 0 0); 			       // linear velocity (m/s)
+		velocity 		     realx3 (0 0 0);     // linear velocity (m/s)
 
-		acceleration 	realx3 	(0 0 0); 			       // linear acceleration (m/s2)
+		acceleration 	     realx3 (0 0 0); 	 // linear acceleration (m/s2)
 
-		rVelocity 		realx3 	(0 0 0); 			       // rotational velocity (rad/s)
+		rVelocity 		     realx3 (0 0 0); 	 // rotational velocity (rad/s)
 
-		shapeName 		word	smallSphere; 		       // name of the particle shape 
+		shapeName 		word    smallSphere; 	 // name of the particle shape 
 	}
 
 	selectors
 	{
 		shapeAssigne
 		{
-			selector 	stridedRange; 				       // other options: box, cylinder, sphere, randomPoints 
+			selector 	stridedRange; 	         // other options: box, cylinder, sphere, randomPoints 
  
 			stridedRangeInfo
 			{
-				begin 	0;			 				       // begin index of points
+				begin 	    0;			         // begin index of points
 
-				end 	30000;		 				       // end index of points 
+				end 	30000;		 			 // end index of points 
 
-				stride  3;			 				       // stride for selector 
+				stride      3;			 		 // stride for selector 
 			}
 
-			fieldValue  							       // fields that the selector is applied to 
+			fieldValue  					     // fields that the selector is applied to 
 			{
-				shapeName 	word 	largeSphere;           // sets shapeName of the selected points to largeSphere
+				shapeName 	word    largeSphere; // sets shapeName of the selected points to largeSphere
 			}
 		}
 	}
 }
 
-positionParticles 				                           // positions particles 
+positionParticles 				                 // positions particles 
 {
-	method ordered;                                        // other options: random and empty
+	method ordered;                              // other options: random and empty
 
 	orderedInfo
 	{		
-		diameter  0.005; 	                               // diameter of particles
+		diameter  0.005; 	                     // diameter of particles
 
-		numPoints 30000; 	                               // number of particles in the simulation
+		numPoints 30000; 	                     // number of particles in the simulation
 
-		axisOrder (z x y);                                 // axis order for filling the space with particles
+		axisOrder (z x y);                       // axis order for filling the space with particles
 	}
 	
-	regionType cylinder;                                   // other options: box and sphere  
+	regionType cylinder;                         // other options: box and sphere  
 
-	cylinderInfo  				                           // cylinder information for positioning particles 
+	cylinderInfo  				                 // cylinder information for positioning particles 
 	{
-		p1 (0.0 0.0 0.003); 	                           // begin point of cylinder axis
+		p1 (0.0 0.0 0.003); 	                 // begin point of cylinder axis
 
-		p2 (0.0 0.0 0.097);		                           // end point of cylinder axis
+		p2 (0.0 0.0 0.097);		                 // end point of cylinder axis
 
-		radius 		0.117;		                           // radius of cylinder 
+		radius 		0.117;		                 // radius of cylinder 
 	}
 }

tutorials/sphereGranFlow/binarySystemOfParticles/settings/settingsDict

@@ -8,28 +8,31 @@ fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
 run 	binarySystemofParticles;
 
-dt       0.00001; 	       // time step for integration (seconds)
+dt       0.00001; 	                             // time step for integration (seconds)
 
-startTime    0.0; 		   // start time for simulation 
+startTime    0.0; 		                         // start time for simulation 
 
-endTime     10.0; 		   // end time for simulation 
+endTime     10.0; 		                         // end time for simulation 
 
-saveInterval 0.1; 		   // time interval for saving the simulation
+saveInterval 0.1; 		                         // time interval for saving the simulation
 
-timePrecision  6;		   // maximum number of digits for time folder 
+timePrecision  6;		                         // maximum number of digits for time folder 
 
-g 	  (0 -9.8 0);          // gravity vector (m/s2) 
+g 	  (0 -9.8 0);                                // gravity vector (m/s2) 
 
 // save necessary (i.e., required) data on disk
-includeObjects (diameter);                                                 
+
+includeObjects (diameter);  
+
 // exclude unnecessary data from saving on disk
+
 excludeObjects (rVelocity.dy1 pStructPosition.dy1 pStructVelocity.dy1); 
 
-integrationMethod 		AdamsBashforth2; 	      // integration method 
+integrationMethod 		AdamsBashforth2; 	     // integration method 
 
-writeFormat 			          ascii;          // data writting format (ascii or binary)
+writeFormat 			          ascii;         // data writting format (ascii or binary)
 
-timersReport 			            Yes;  	      // report timers
+timersReport 			            Yes;  	     // report timers
 
-timersReportInterval   	           0.01;	      // time interval for reporting timers
+timersReportInterval   	           0.01;	     // time interval for reporting timers
  

tutorials/sphereGranFlow/layeredSiloFilling/caseSetup/interaction

@@ -6,11 +6,11 @@ objectName interaction;
 objectType dicrionary;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
-materials   (lightMat heavyMat wallMat);          // a list of materials names
+materials   (lightMat heavyMat wallMat);         // a list of materials names
 
-densities            (1000 1500.0 2500);          // density of materials [kg/m3]
+densities            (1000 1500.0 2500);         // density of materials [kg/m3]
 
-contactListType   sortedContactList;
+contactListType       sortedContactList;
 
 contactSearch
 {
@@ -29,39 +29,40 @@ model
 {
    contactForceModel nonLinearLimited;
 
-   rollingFrictionModel normal;
+   rollingFrictionModel        normal;
 
    /*
+
    Property (lightMat-lightMat   lightMat-heavyMat    lightMat-wallMat
                                  heavyMat-heavyMat    heavyMat-wallMat
-                                                      wallMat-wallMat );
+                                                      wallMat-wallMat  );
    */
 
-   Yeff (1.0e6 1.0e6 1.0e6            // Young modulus [Pa]
+   Yeff (1.0e6 1.0e6 1.0e6                       // Young modulus [Pa]
                1.0e6 1.0e6
                      1.0e6);
    
-   Geff (0.8e6 0.8e6 0.8e6            // Shear modulus [Pa]
+   Geff (0.8e6 0.8e6 0.8e6                       // Shear modulus [Pa]
                0.8e6 0.8e6
                      0.8e6);
 
-   nu    (0.25 0.25  0.25             // Poisson's ratio [-]
+   nu    (0.25 0.25  0.25                        // Poisson's ratio [-]
                0.25  0.25
                      0.25);
 
-   en (0.97  0.97    0.85             // coefficient of normal restitution
+   en (0.97  0.97    0.85                        // coefficient of normal restitution
              0.97    0.85
                      1.00);
 
-   et (1.0   1.0  1.0                 // coefficient of tangential restitution 
+   et (1.0   1.0  1.0                            // coefficient of tangential restitution 
              1.0  1.0
                   1.0);
 
-   mu (0.65   0.65 0.35               // dynamic friction 
+   mu (0.65   0.65 0.35                          // dynamic friction 
               0.65 0.35
                    0.35);
 
-   mur (0.1  0.1 0.1                  // rolling friction 
+   mur (0.1  0.1 0.1                             // rolling friction 
              0.1 0.1
                  0.1);       
 }

tutorials/sphereGranFlow/layeredSiloFilling/caseSetup/particleInsertion

@@ -6,9 +6,9 @@ objectName particleInsertion;
 objectType dicrionary;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
-active           Yes;          // is insertion active -> yes or no
+active           Yes;                            // is insertion active -> yes or no
 
-checkForCollision No;          // is checked          -> yes or no
+checkForCollision No;                            // is checked          -> yes or no
 
 /*
    one layers of particles are packed 

tutorials/sphereGranFlow/layeredSiloFilling/settings/domainDict

@@ -6,7 +6,7 @@ objectName 	domainDict;
 objectType 	dictionary;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
-globalBox                    // Simulation domain: every particles that goes outside this domain will be deleted
+globalBox                                        // Simulation domain: every particles that goes outside this domain will be deleted
 {
 	min (-0.11 -0.11 -0.11);
 	
@@ -16,44 +16,49 @@ globalBox                    // Simulation domain: every particles that goes out
 boundaries
 {
 	// Determines how often (how many iterations) do you want to 
+
 	// rebuild the list of particles in the neighbor list 
+
 	// of all boundaries in the simulation domain 
+
 	neighborListUpdateInterval 30;
 	
 	// Determines how often do you want to update the new changes in the boundary
+
 	updateInterval             10;
 
 	// The distance from the boundary plane within which particles are marked to be in the boundary list 
+
 	neighborLength          0.004;
 
 	left
 	{
-		type exit;	         // other options: periodict, reflective 
+		type exit;	                             // other options: periodict, reflective 
 	}
 
 	right
 	{
-		type exit;           // other options: periodict, reflective 
+		type exit;                               // other options: periodict, reflective 
 	}
 
 	bottom
 	{
-		type exit;           // other options: periodict, reflective 
+		type exit;                               // other options: periodict, reflective 
 	}
 
 	top
 	{
-		type exit;           // other options: periodict, reflective 
+		type exit;                               // other options: periodict, reflective 
 	}
 
 	rear
 	{
-		type exit;           // other options: periodict, reflective 
+		type exit;                               // other options: periodict, reflective 
 	}
 
 	front
 	{
-		type exit;           // other options: periodict, reflective 
+		type exit;                               // other options: periodict, reflective 
 	}
 }    
 

tutorials/sphereGranFlow/layeredSiloFilling/settings/geometryDict

@@ -6,7 +6,7 @@ objectName 	geometryDict;
 objectType 	dictionary;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
-motionModel stationary;          // motion model can be rotatingAxis or stationary or vibrating
+motionModel            stationary;               // motion model can be rotatingAxis or stationary or vibrating
 
 stationaryInfo
 {
@@ -17,36 +17,36 @@ surfaces
 {
 	cylinderShell
 	{
-		type 		 cylinderWall;  	    // other options: cuboidWall and planeWall
+		type 		 cylinderWall;  	         // other options: cuboidWall and planeWall
 
-		p1 			(0.0 0.0 0.0);          // begin point of cylinder axis
+		p1 			(0.0 0.0 0.0);               // begin point of cylinder axis
 
-		p2 			(0.0 0.0 0.4);          // end point of cylinder axis
+		p2 			(0.0 0.0 0.4);               // end point of cylinder axis
 
-		radius1 	          0.1;		    // radius at p1
+		radius1 	          0.1;		         // radius at p1
  
-		radius2 	          0.1;		    // radius at p2
+		radius2 	          0.1;		         // radius at p2
 
-		resolution 	           36;      	// number of divisions
+		resolution 	           36;      	     // number of divisions
 
-		material 	      wallMat;   	    // material name of this wall
+		material 	      wallMat;   	         // material name of this wall
 	}
 
 	coneShell
 	{
-		type 		 cylinderWall;  	    // other options: cuboidWall and planeWall
+		type 		 cylinderWall;  	         // other options: cuboidWall and planeWall
 
-		p1 		   (0.0 0.0 -0.1);          // begin point of cylinder axis
+		p1 		   (0.0 0.0 -0.1);               // begin point of cylinder axis
 
-		p2 			(0.0 0.0 0.0);          // end point of cylinder axis
+		p2 			(0.0 0.0 0.0);               // end point of cylinder axis
 
-		radius1 	         0.02;		    // radius at p1
+		radius1 	         0.02;		         // radius at p1
 
-		radius2 	          0.1;		    // radius at p2
+		radius2 	          0.1;		         // radius at p2
+  
+		resolution 	           36;      	     // number of divisions
 
-		resolution 	           36;      	// number of divisions
-
-		material 	      wallMat;   	    // material name of this wall
+		material 	      wallMat;   	         // material name of this wall
 	}
 
 	/*
@@ -55,17 +55,17 @@ surfaces
 
 	exitGate
 	{
-		type planeWall;                     // other options: cuboidWall and cylinderWall
+		type        planeWall;                   // other options: cuboidWall and cylinderWall
 
-		p1 (-0.02 -0.02 -0.1);              // first point of the wall
+		p1 (-0.02 -0.02 -0.1);                   // first point of the wall
 
-		p2 ( 0.02 -0.02 -0.1);              // second point of the wall
+		p2 ( 0.02 -0.02 -0.1);                   // second point of the wall
 
-		p3 ( 0.02  0.02 -0.1);              // third point of the wall
+		p3 ( 0.02  0.02 -0.1);                   // third point of the wall
 
-		p4 (-0.02  0.02 -0.1);              // fourth point of the wall
+		p4 (-0.02  0.02 -0.1);                   // fourth point of the wall
 		 
-		material wallMat;                   // material name of the wall 
+		material      wallMat;                   // material name of the wall 
 	}
 }
 

tutorials/sphereGranFlow/layeredSiloFilling/settings/particlesDict

@@ -9,37 +9,39 @@ fileFormat  ASCII;
 setFields
 {
 	/*
-	    	Default value for fields defined for particles
+	    Default value for fields defined for particles
+
 		These fields should always be defined for simulations with 
+
 		spherical particles.
 	*/
 
 	defaultValue 
 	{
-		velocity 		realx3 	    (0 0 0); 		  // linear velocity (m/s)
+		velocity 		realx3 	    (0 0 0);     // linear velocity (m/s)
 
-		acceleration 	realx3 	    (0 0 0); 		  // linear acceleration (m/s2)
+		acceleration 	realx3 	    (0 0 0);     // linear acceleration (m/s2)
 
-		rVelocity 		realx3 	    (0 0 0);  		  // rotational velocity (rad/s)
+		rVelocity 		realx3 	    (0 0 0);  	 // rotational velocity (rad/s)
 
-		shapeName 		word	lightSphere; 	      // name of the particle shape 
+		shapeName 		word	lightSphere; 	 // name of the particle shape 
 	}
 
 	selectors
 	{}
 }
 
-positionParticles                                          // positions particles 
+positionParticles                                // positions particles 
 {
-	method                empty;          // other options: ordered and random 
+	method                empty;                 // other options: ordered and random 
 
-	regionType              box;          // other options: cylinder and sphere  
+	regionType              box;                 // other options: cylinder and sphere  
 
-	boxInfo                               // box region for positioning particles 
+	boxInfo                                      // box region for positioning particles 
 	{
-		min (-0.08 -0.08 0.015);          // lower corner point of the box 
+		min (-0.08 -0.08 0.015);                 // lower corner point of the box 
 
-		max ( 0.08  0.08 0.098);	      // upper corner point of the box 
+		max ( 0.08  0.08 0.098);	             // upper corner point of the box 
 	} 
 }
 

tutorials/sphereGranFlow/layeredSiloFilling/settings/settingsDict

@@ -8,22 +8,26 @@ fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
 run   layerdSiloFilling;
 
-dt 				0.00001; 	      // time step for integration (s)
+dt 				0.00001; 	                     // time step for integration (s)
 
-startTime 		    0.0; 		  // start time for simulation 
+startTime 		    0.0; 		                 // start time for simulation 
 
-endTime 		    5.0;	      // end time for simulation 
+endTime 		    5.0;	                     // end time for simulation 
 
-saveInterval 	   0.05; 		  // time interval for saving the simulation
+saveInterval 	   0.05; 		                 // time interval for saving the simulation
 
-timePrecision         6;	      // maximum number of digits for time folder 
+timePrecision         6;	                     // maximum number of digits for time folder 
 
-g 			 (0 0 -9.8);          // gravity vector (m/s2) 
+g 			 (0 0 -9.8);                         // gravity vector (m/s2) 
 
 // save data objects that are not automatically saved on disk.
+
 // overrides the default behavior 
-includeObjects (diameter);                                                 
+
+includeObjects (diameter);    
+
 // exclude unnecessary data from saving on disk
+
 excludeObjects (rVelocity.dy1 pStructPosition.dy1 pStructVelocity.dy1); 
 
 integrationMethod 		AdamsBashforth2; 	      // integration method