Reasonable indentation for parabolic inlet profile

precice · Mar 26, 2024 · fc05b0e · fc05b0e
1 parent 401c1bb
commit fc05b0e
Showing 1 changed file with 16 additions and 14 deletions.
diff --git a/turek-hron-fsi3/fluid-openfoam/0/U b/turek-hron-fsi3/fluid-openfoam/0/U
@@ -16,7 +16,7 @@ boundaryField
     cylinder
     {
         type            noSlip;
-    }    
+    }
 
     flap
     {
@@ -43,24 +43,26 @@ boundaryField
 
         code
         #{
-            const       vectorField& Cf = patch().Cf();
-            vectorField&    field = *this;
-            
-            const scalar pi=3.14159265358979;
-            const scalar Umean=2.0;            
+            const vectorField& Cf = patch().Cf();
+            vectorField& field = *this;
+
+            const scalar pi = 3.14159265358979;
+            const scalar Umean = 2.0;
             const scalar t = this->db().time().value();
-            if(t<2){
-                forAll(Cf,faceI)
+            if (t < 2)
+            {
+                forAll(Cf, faceI)
                 {
-                    const scalar y=Cf[faceI][1];
-                    field[faceI]=vector((1.5*Umean*4.0/0.1681)*y*(0.41-y)*((1-cos(pi/2 *t))/2),0,0);
+                    const scalar y = Cf[faceI][1];
+                    field[faceI] = vector((1.5 * Umean * 4.0 / 0.1681) * y * (0.41 - y) * ((1 - cos(pi / 2 * t)) / 2), 0, 0);
                 }
-            }else
+            }
+            else
             {
-                forAll(Cf,faceI)
+                forAll(Cf, faceI)
                 {
-                    const scalar y=Cf[faceI][1];
-                    field[faceI]=vector((1.5*Umean*4.0/0.1681)*y*(0.41-y),0,0);
+                    const scalar y = Cf[faceI][1];
+                    field[faceI] = vector((1.5 * Umean * 4.0 / 0.1681) * y * (0.41 - y), 0, 0);
                 }
             }
         #};