createFields.H
Go to the documentation of this file.
1  Info<< "Reading field psi\n" << endl;
3  (
4  IOobject
5  (
6  "psi",
7  runTime.timeName(),
8  mesh,
9  IOobject::MUST_READ,
10  IOobject::AUTO_WRITE
11  ),
12  mesh
13  );
14 
15  Info<< "Reading transportProperties\n" << endl;
16 
17  IOdictionary transportProperties
18  (
19  IOobject
20  (
21  "transportProperties",
22  runTime.constant(),
23  mesh,
24  IOobject::MUST_READ_IF_MODIFIED,
25  IOobject::NO_WRITE
26  )
27  );
28 
29  List<magnet> magnets(transportProperties.lookup("magnets"));
30 
32  (
33  IOobject
34  (
35  "murf",
36  runTime.timeName(),
37  mesh
38  ),
39  mesh,
40  1
41  );
42 
44  (
45  IOobject
46  (
47  "Mrf",
48  runTime.timeName(),
49  mesh
50  ),
51  mesh,
52  dimensionedScalar("Mr", dimensionSet(0, 1, 0, 0, 0, 1, 0), 0)
53  );
54 
56  {
57  label magnetZonei = mesh.faceZones().findZoneID(magnets[i].name());
58 
59  if (magnetZonei == -1)
60  {
61  FatalIOErrorIn(args.executable().c_str(), transportProperties)
62  << "Cannot find faceZone for magnet " << magnets[i].name()
63  << exit(FatalIOError);
64  }
65 
66  const labelList& faces = mesh.faceZones()[magnetZonei];
67 
68  const scalar muri = magnets[i].mur();
69  const scalar Mri = magnets[i].Mr().value();
70  const vector& orientationi = magnets[i].orientation();
71 
72  const surfaceVectorField& Sf = mesh.Sf();
73 
74  forAll(faces, i)
75  {
76  label facei = faces[i];
77  murf[facei] = muri;
78  Mrf[facei] = Mri*(orientationi & Sf[facei]);
79  }
80  }
forAll(Y, i)
Definition: createFields.H:99
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
errorManipArg< error, int > exit(error &err, const int errNo=1)
Definition: errorManip.H:124
surfaceScalarField murf(IOobject("murf", runTime.timeName(), mesh), mesh, 1)
GeometricField< vector, fvsPatchField, surfaceMesh > surfaceVectorField
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:253
Vector< scalar > vector
A scalar version of the templated Vector.
Definition: vector.H:49
const word & executable() const
Name of executable without the path.
Definition: argListI.H:30
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
dynamicFvMesh & mesh
List< label > labelList
A List of labels.
Definition: labelList.H:56
#define FatalIOErrorIn(functionName, ios)
Report an error message using Foam::FatalIOError.
Definition: error.H:325
word name(const complex &)
Return a string representation of a complex.
Definition: complex.C:47
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
surfaceScalarField Mrf(IOobject("Mrf", runTime.timeName(), mesh), mesh, dimensionedScalar("Mr", dimensionSet(0, 1, 0, 0, 0, 1, 0), 0))
messageStream Info
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Foam::argList args(argc, argv)
Info<< "Reading field psi\n"<< endl;volScalarField psi(IOobject("psi", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading transportProperties\n"<< endl;IOdictionary transportProperties(IOobject("transportProperties", runTime.constant(), mesh, IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE));List< magnet > magnets(transportProperties.lookup("magnets"))
const volScalarField & psi
Definition: createFields.H:30
IOerror FatalIOError