File tree 1 file changed +16
-14
lines changed
turek-hron-fsi3/fluid-openfoam/0
1 file changed +16
-14
lines changed Original file line number Diff line number Diff line change @@ -16,7 +16,7 @@ boundaryField
16
16
cylinder
17
17
{
18
18
type noSlip;
19
- }
19
+ }
20
20
21
21
flap
22
22
{
@@ -43,24 +43,26 @@ boundaryField
43
43
44
44
code
45
45
#{
46
- const vectorField& Cf = patch().Cf();
47
- vectorField& field = *this;
48
-
49
- const scalar pi= 3.14159265358979;
50
- const scalar Umean= 2.0;
46
+ const vectorField& Cf = patch().Cf();
47
+ vectorField& field = *this;
48
+
49
+ const scalar pi = 3.14159265358979;
50
+ const scalar Umean = 2.0;
51
51
const scalar t = this->db().time().value();
52
- if(t<2){
53
- forAll(Cf,faceI)
52
+ if (t < 2)
53
+ {
54
+ forAll(Cf, faceI)
54
55
{
55
- const scalar y= Cf[faceI][1];
56
- field[faceI]= vector((1.5* Umean* 4.0/ 0.1681)*y* (0.41-y)* ((1- cos(pi/ 2 *t))/ 2),0, 0);
56
+ const scalar y = Cf[faceI][1];
57
+ field[faceI] = vector((1.5 * Umean * 4.0 / 0.1681) * y * (0.41 - y) * ((1 - cos(pi / 2 * t)) / 2), 0, 0);
57
58
}
58
- }else
59
+ }
60
+ else
59
61
{
60
- forAll(Cf,faceI)
62
+ forAll(Cf, faceI)
61
63
{
62
- const scalar y= Cf[faceI][1];
63
- field[faceI]= vector((1.5* Umean* 4.0/ 0.1681)*y* (0.41- y),0, 0);
64
+ const scalar y = Cf[faceI][1];
65
+ field[faceI] = vector((1.5 * Umean * 4.0 / 0.1681) * y * (0.41 - y), 0, 0);
64
66
}
65
67
}
66
68
#};
You can’t perform that action at this time.
0 commit comments