56 int main(
int argc,
char *argv[])
67 scalar b1 =
j1(swirlProfile).
value();
68 scalar b2 = 2.0*b1/swirlProfile.value() -
j0(swirlProfile).
value();
70 scalar omega = 0.125*(Vphi*bore*swirlProfile/b2).value();
72 scalar cylinderRadius = 0.5*bore.value();
80 if (r <= cylinderRadius)
82 scalar
b =
j1(swirlProfile*r/cylinderRadius).
value();
83 scalar vEff = omega*
b;
85 U[celli] = ((vEff/r)*(c & yT))*xT + (-(vEff/r)*(c & xT))*yT;
86 Umax =
max(Umax,
mag(
U[celli]));