// Construct from components
dsmcMixedDiffuseSpecularWallFieldPatch::dsmcMixedDiffuseSpecularWallFieldPatch
(
Time& t,
const polyMesh& mesh,
dsmcCloud& cloud,
const dictionary& dict
)
:
dsmcPatchBoundary(t, mesh, cloud, dict),
propsDict_(dict.subDict(typeName + "Properties")),
boundaryT_
(
volScalarField
(
IOobject
(
"boundaryT",
mesh_.time().timeName(),
mesh_,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh_
)
),
boundaryU_
(
volVectorField
(
IOobject
(
"boundaryU",
mesh_.time().timeName(),
mesh_,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh_
)
),
diffuseFraction_(readScalar(propsDict_.lookup("diffuseFraction")))
{
writeInTimeDir_ = false;
writeInCase_ = false;
measurePropertiesAtWall_ = true;
}
int main(int argc, char *argv[])
{
#include "setRootCase.H"
#include "createTime.H"
#include "createMesh.H"
#include "createGravity.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Info<< nl << "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
<< nl << endl;
Info<< "End\n" << endl;
Info<< "Reading field thetaf\n" << endl;
surfaceScalarField thetaf
(
IOobject
(
"thetaf",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
surfaceScalarField bf = surfaceScalarField("bf", thetaf * g.unitFaceNormal());
bf.write();
volVectorField b = volVectorField("b", fvc::reconstruct(bf * mesh.magSf()));
b.write();
volScalarField theta = volScalarField("theta", b & g.unit());
theta.write();
return 0;
}
