Style: Changed formatting

Author: jibranhaider

applications/solvers/solidFoam/Make/options

@@ -4,8 +4,6 @@ EXE_INC = \
     -I$(PWD)/../../../src/models/lnInclude \
     -I$(LIB_SRC)/finiteVolume/lnInclude \
     -I$(LIB_SRC)/meshTools/lnInclude \
-    
-
 
 EXE_LIBS = \
     -L$(FOAM_USER_LIBBIN) \
@@ -14,4 +12,4 @@ EXE_LIBS = \
     -lModels \
     -L$(FOAM_LIBBIN) \
     -lfiniteVolume \
-    -lmeshTools
+    -lmeshTools

applications/solvers/solidFoam/compile

@@ -4,4 +4,4 @@ printf "\n\n........................................\n"
 printf "Compiling solidFoam solver\n"
 printf "........................................\n\n"
 
-wclean && wmake
+wclean && wmake

applications/solvers/solidFoam/createFields.H

@@ -123,7 +123,6 @@
     volScalarField Up_time (model.Up_nonLinear());
 
 
-
 //-----------------------------------------------------------------------------------
 // Gradients
 //-----------------------------------------------------------------------------------
@@ -161,7 +160,6 @@
     grad.gradient(P, PxGrad, PyGrad, PzGrad);
 
 
-
 //-----------------------------------------------------------------------------------
 // Reconstruction
 //-----------------------------------------------------------------------------------
@@ -201,7 +199,6 @@
 // Riemann solver
 //-----------------------------------------------------------------------------------
 
-
     // Contact linear momentum
     surfaceVectorField lmC
     (
@@ -240,23 +237,19 @@
     // Constrained class
     interpolationSchemes interpolate(mesh); 
 
-
     // Cell averaged linear momentum
     volVectorField lmR (interpolate.surfaceToCell(lmC));                                     
 
-
     // Local gradient of cell averaged linear momentum
     volTensorField lmRgrad (grad.localGradient(lmR, lmC, AvCinv));  
 
-
     // Constrained fluxes
     interpolate.volToPoint(lmR, lmRgrad, lmN);      
     #include "strongBCs.H"
     lmN.correctBoundaryConditions();
     interpolate.pointToSurface(lmN, lmC);       
 
 
-
 //-----------------------------------------------------------------------------------
 // Angular momentum
 //-----------------------------------------------------------------------------------

applications/solvers/solidFoam/gEqns.H

@@ -6,18 +6,13 @@ if (angularMomentumConservation == "yes")
     AM.AMconservation( rhsLm, rhsLm1, rhsAm, RKstage);
 }
 
-
 // Update coordinates
 x += deltaT * (lm/rho);
 xF += deltaT * (lmC/rho);                       
 xN += deltaT * (lmN/rho);
 
-
 // Update linear momentum
 lm += deltaT * rhsLm;   
 
-
 // Update deformation gradient tensor
-F += deltaT * fvc::surfaceIntegrate( (lmC/rho) * Sf );  
-
-
+F += deltaT * fvc::surfaceIntegrate( (lmC/rho) * Sf );  

applications/solvers/solidFoam/meshData.H

@@ -5,15 +5,13 @@
     // Point Mesh (pMesh)
         pointMesh pMesh(mesh);
 
-
     // Material face area
         const surfaceScalarField& magSf (mesh.magSf());
 
-
     // Material face area normal vector
         const surfaceVectorField& Sf (mesh.Sf());
 
-        
+
 //----------------------------------------------------------------------
 // Mesh coordinates
 //----------------------------------------------------------------------
@@ -40,7 +38,6 @@
         xN.primitiveFieldRef() = mesh.points();
 
 
-
     // Material nodal coordinates
         pointVectorField xN_0
         (
@@ -57,7 +54,6 @@
         surfaceVectorField xF = mesh.Cf();
 
 
-
 //----------------------------------------------------------------------
 // Normals
 //----------------------------------------------------------------------
@@ -70,8 +66,6 @@
         );          
 
 
-
-
 // EDGE TO POINT CONNECTIVITY (edges)
 //-----------------------------------
     const edgeList& edges = mesh.edges();
@@ -105,8 +99,6 @@
     }
 
 
-
-
 // BOUNDARY PATCHES
 //-----------------
 
@@ -117,7 +109,6 @@
     const label& symmetricZpatchID = bm.findPatchID("symmetricZ");  
 
 
-
 // A INVERSE FOR CALCULATION OF GRADIENT (Ainv)
 //---------------------------------------------
 
@@ -129,12 +120,8 @@
     );
 
 
-    
-
 // Ainv FOR CALCULATION OF GRADIENT FROM CONTACT VALUES (AvCinv)
 //--------------------------------------------------------------
 
     // Initialise AvCinv    
-    volTensorField AvCinv(Ainv);
-
-
+    volTensorField AvCinv(Ainv);

applications/solvers/solidFoam/riemannSolver.H

@@ -1,8 +1,4 @@
-//-------------------------
-// ACOUSTIC RIEMANN SOLVER 
-//------------------------- 
-
-
+// Acoustic Riemann solver
 lmC = 0.5*(lm_M+lm_P) + ( 0.5 * ((nCn/fvc::interpolate(Up)) + (iMnCn/fvc::interpolate(Us))) & (t_P-t_M) );  
 tC  = 0.5*(t_M+t_P) + ( 0.5 * ((fvc::interpolate(Up)*nCn) + (fvc::interpolate(Us)*iMnCn)) & (lm_P-lm_M) );
 
@@ -22,6 +18,7 @@ if( Pstream::parRun() )
 
     forAll(mesh.boundary(), patchID)
     {
+        // Riemann solver for inter-processor boundaries        
         if (mesh.boundary()[patchID].coupled())             
         {           
             vectorField lm_nei = lm.boundaryField()[patchID].patchNeighbourField();
@@ -78,8 +75,7 @@ if( Pstream::parRun() )
 }
 
 
-
-
+// Compute boundary values
 lm_b.correctBoundaryConditions();
 t_b.correctBoundaryConditions();
 

applications/solvers/solidFoam/solidFoam.C

@@ -96,4 +96,4 @@ int main(int argc, char *argv[])
     return 0;
 }
 
-// ************************************************************************* //
+// ************************************************************************* //

applications/solvers/solidFoam/strongBCs.H

@@ -64,4 +64,4 @@
                 lmN[nodeID] = iMnCn_ & lmN[nodeID];
             }
         }
-    }
+    }

applications/solvers/solidFoam/updateVariables.H

@@ -1,4 +1,3 @@
-
 //-----------------------------------------------------------------------------------
 // Deformation dependent properties
 //-----------------------------------------------------------------------------------           
@@ -8,13 +7,11 @@
     nCn = n*n;
     iMnCn = (tensor::I)-(nCn);
 
-
     // Strain measures
     Finv = inv(F);       
     H = det(F)*Finv.T();    
     J = det(F); 
 
-
     // Wave speeds
     Up_time = model.Up_nonLinear();     
 
@@ -25,7 +22,6 @@
     P = model.piola();
 
 
-
 //-----------------------------------------------------------------------------------
 // Linear reconstruction
 //-----------------------------------------------------------------------------------           
@@ -34,23 +30,20 @@
     grad.gradient(lm, lmGrad);      
     grad.gradient(P, PxGrad, PyGrad, PzGrad);
 
-
     // Reconstruction   
     grad.reconstruct( lm, lmGrad, lm_M, lm_P ); 
     grad.reconstruct( P, PxGrad, PyGrad, PzGrad, P_M, P_P) ;
     t_M = P_M & N;
     t_P = P_P & N;
 
 
-
 //-----------------------------------------------------------------------------------
 // Riemann solver
 //-----------------------------------------------------------------------------------
 
     #include "riemannSolver.H"              
 
 
-
 //-----------------------------------------------------------------------------------
 // Constrained approach
 //-----------------------------------------------------------------------------------
@@ -62,10 +55,5 @@
     #include "strongBCs.H"
     lmN.correctBoundaryConditions();    
 
-
     // Constrained fluxes
-    interpolate.pointToSurface(lmN, lmC);
-
-
-
-
+    interpolate.pointToSurface(lmN, lmC);

applications/utilities/initialConditions/initialConditions.C

@@ -107,4 +107,4 @@ int main(int argc, char *argv[])
     return 0;
 }
 
-// ************************************************************************* //
+// ************************************************************************* //