FreeFOAM The Cross-Platform CFD Toolkit
createFields.H
Go to the documentation of this file.
1 Info<< "Reading combustion properties\n" << endl;
2 
3 IOdictionary engineGeometry
4 (
5  IOobject
6  (
7  "engineGeometry",
8  runTime.constant(),
9  mesh,
10  IOobject::MUST_READ,
11  IOobject::NO_WRITE
12  )
13 );
14 
15 vector swirlAxis
16 (
17  engineGeometry.lookup("swirlAxis")
18 );
19 
20 vector swirlCenter
21 (
22  engineGeometry.lookup("swirlCenter")
23 );
24 
25 dimensionedScalar swirlRPMRatio
26 (
27  engineGeometry.lookup("swirlRPMRatio")
28 );
29 
30 dimensionedScalar swirlProfile
31 (
32  engineGeometry.lookup("swirlProfile")
33 );
34 
36 (
37  engineGeometry.lookup("bore")
38 );
39 
41 (
42  engineGeometry.lookup("rpm")
43 );
44 
45 
46 Info<< "Reading field U\n" << endl;
48 (
49  IOobject
50  (
51  "U",
52  runTime.timeName(),
53  mesh,
54  IOobject::MUST_READ,
55  IOobject::AUTO_WRITE
56  ),
57  mesh
58 );
59 
60 vector zT = swirlAxis;
61 vector yT = vector(0, zT.z(), -zT.y());
62 vector xT = vector(zT.y()*zT.y() + zT.z()*zT.z(), -zT.x()*zT.y(), -zT.x()*zT.z());
63 
64 // if swirl is around (1, 0, 0) we have to find another transformation
65 if (mag(yT) < SMALL)
66 {
67  yT = vector(zT.y(), -zT.x(), 0);
68  xT = vector(-zT.x()*zT.z(), -zT.y()*zT.z(), zT.x()*zT.x() + zT.y()*zT.y());
69 }
70 
71 //swirlAxis doesn't have to be of unit length.
72 xT /= mag(xT);
73 yT /= mag(yT);
74 zT /= mag(zT);
75 
76 // ************************ vim: set sw=4 sts=4 et: ************************ //