41 int main(
int argc,
char *argv[])
43 timeSelector::addOptions();
52 runTime.setTime(timeDirs[timeI], timeI);
53 Info<<
"Time = " << runTime.timeName() <<
endl;
56 #include "createFields.H" 70 const surfaceScalarField::Boundary& patchHeatFlux =
73 const volScalarField::Boundary& patchRadHeatFlux =
76 const surfaceScalarField::Boundary& magSf =
77 mesh.magSf().boundaryField();
79 Info<<
"\nWall heat fluxes [W]" <<
endl;
88 <<
" convective: " << convFlux << endl
89 <<
" radiative: " << radFlux << endl
90 <<
" total: " << convFlux + radFlux <<
endl;
107 volScalarField::Boundary& wallHeatFluxBf =
115 wallHeatFlux.write();
133 heatFlux.dimensions(),
138 volScalarField::Boundary& totalWallHeatFluxBf =
139 totalWallHeatFlux.boundaryFieldRef();
143 totalWallHeatFluxBf[
patchi] =
147 totalWallHeatFlux.write();
autoPtr< compressible::turbulenceModel > turbulence
List< instant > instantList
List of instants.
#define forAll(list, i)
Loop across all elements in list.
Ostream & endl(Ostream &os)
Add newline and flush stream.
tmp< surfaceScalarField > interpolate(const RhoType &rho)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
psiReactionThermo & thermo
Type gSum(const FieldField< Field, Type > &f)
const Boundary & boundaryField() const
Return const-reference to the boundary field.
const dimensionedScalar h
Planck constant.
Boundary & boundaryFieldRef()
Return a reference to the boundary field.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
volScalarField Qr(IOobject("Qr", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::NO_WRITE), mesh, dimensionedScalar("Qr", dimMass/pow3(dimTime), 0.0))
Foam::argList args(argc, argv)
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)