FreeFOAM The Cross-Platform CFD Toolkit
createFields.H
Go to the documentation of this file.
1  Info<< "Reading physicalProperties\n" << endl;
2 
3  IOdictionary physicalProperties
4  (
5  IOobject
6  (
7  "physicalProperties",
8  runTime.constant(),
9  mesh,
10  IOobject::MUST_READ,
11  IOobject::NO_WRITE
12  )
13  );
14 
15  dimensionedScalar epsilon0
16  (
17  physicalProperties.lookup("epsilon0")
18  );
19 
21  (
22  physicalProperties.lookup("k")
23  );
24 
25 
26  Info<< "Reading field phi\n" << endl;
28  (
29  IOobject
30  (
31  "phi",
32  runTime.timeName(),
33  mesh,
34  IOobject::MUST_READ,
35  IOobject::AUTO_WRITE
36  ),
37  mesh
38  );
39 
40 
41  Info<< "Reading field rho\n" << endl;
43  (
44  IOobject
45  (
46  "rho",
47  runTime.timeName(),
48  mesh,
49  IOobject::MUST_READ,
50  IOobject::AUTO_WRITE
51  ),
52  mesh
53  );
54 
55 
56  Info<< "Calculating field rhoFlux\n" << endl;
57  surfaceScalarField rhoFlux
58  (
59  IOobject
60  (
61  "rhoFlux",
62  runTime.timeName(),
63  mesh,
64  IOobject::NO_READ,
65  IOobject::NO_WRITE
66  ),
67  -k*mesh.magSf()*fvc::snGrad(phi)
68  );
69 
70 // ************************ vim: set sw=4 sts=4 et: ************************ //