1 Info<<
"Reading thermal properties\n" <<
endl;
3 IOdictionary thermalProperties
10 IOobject::MUST_READ_IF_MODIFIED,
15 Switch thermalStress(thermalProperties.lookup(
"thermalStress"));
61 autoPtr<volScalarField> CPtr;
63 const dictionary& CDict(thermalProperties.subDict(
"C"));
64 word CType(CDict.lookup(
"type"));
65 if (CType ==
"uniform")
67 scalar CValue(
readScalar(CDict.lookup(
"value")));
77 dimensionSet(0, 2, -2 , -1, 0),
84 else if (CType ==
"field")
86 CIO.readOpt() = IOobject::MUST_READ;
100 <<
"Valid type entries are uniform or field for C" 115 const dictionary& kDict(thermalProperties.subDict(
"k"));
116 word kType(kDict.lookup(
"type"));
117 if (kType ==
"uniform")
119 scalar rhoKValue(
readScalar(kDict.lookup(
"value")));
129 dimensionSet(1, 1, -3 , -1, 0),
136 else if (kType ==
"field")
138 rhoKIO.readOpt() = IOobject::MUST_READ;
152 <<
"Valid type entries are uniform or field for K" 158 autoPtr<volScalarField> alphaPtr;
170 const dictionary& alphaDict(thermalProperties.subDict(
"alpha"));
171 word alphaType(alphaDict.lookup(
"type"));
173 if (alphaType ==
"uniform")
175 scalar alphaValue(
readScalar(alphaDict.lookup(
"value")));
190 else if (alphaType ==
"field")
192 alphaIO.readOpt() = IOobject::MUST_READ;
206 <<
"Valid type entries are uniform or field for alpha" 212 Info<<
"Normalising k : k/rho\n" <<
endl;
215 Info<<
"Calculating thermal coefficients\n" <<
endl;
217 threeKalpha = threeK*
alpha;
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
dimensionedSphericalTensor inv(const dimensionedSphericalTensor &dt)
Ostream & endl(Ostream &os)
Add newline and flush stream.
label k
Boltzmann constant.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Info<< "Reading thermal properties\"<< endl;IOdictionary thermalProperties(IOobject("thermalProperties", runTime.constant(), mesh, IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE));Switch thermalStress(thermalProperties.lookup("thermalStress"));volScalarField threeKalpha(IOobject("threeKalpha", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mesh, dimensionedScalar(dimensionSet(0, 2, -2, -1, 0), 0));volScalarField DT(IOobject("DT", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mesh, dimensionedScalar(dimensionSet(0, 2, -1, 0, 0), 0));autoPtr< volScalarField > rhoKPtr
const dimensionSet dimTemperature(0, 0, 0, 1, 0, 0, 0)
errorManip< error > abort(error &err)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
volScalarField alpha(IOobject("alpha", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))