32 template<
class CloudType>
36 return cloudCopyPtr_();
40 template<
class CloudType>
41 inline const typename CloudType::particleType::constantProperties&
48 template<
class CloudType>
49 inline typename CloudType::particleType::constantProperties&
56 template<
class CloudType>
63 template<
class CloudType>
70 template<
class CloudType>
77 template<
class CloudType>
81 return heatTransferModel_;
85 template<
class CloudType>
93 template<
class CloudType>
100 template<
class CloudType>
107 <<
"Radiation field requested, but radiation model not active" 115 template<
class CloudType>
122 <<
"Radiation field requested, but radiation model not active" 130 template<
class CloudType>
137 <<
"Radiation field requested, but radiation model not active" 145 template<
class CloudType>
152 <<
"Radiation field requested, but radiation model not active" 160 template<
class CloudType>
167 <<
"Radiation field requested, but radiation model not active" 171 return radAreaPT4_();
175 template<
class CloudType>
182 <<
"Radiation field requested, but radiation model not active" 186 return radAreaPT4_();
190 template<
class CloudType>
198 template<
class CloudType>
206 template<
class CloudType>
214 template<
class CloudType>
222 template<
class CloudType>
228 Info<<
"hsTrans min/max = " <<
min(hsTrans()).value() <<
", " 229 <<
max(hsTrans()).value() <<
nl 230 <<
"hsCoeff min/max = " <<
min(hsCoeff()).value() <<
", " 231 <<
max(hsCoeff()).value() <<
endl;
236 if (this->
solution().semiImplicit(
"h"))
240 Vdt(this->
mesh().V()*this->db().time().deltaT());
245 + hsCoeff()/(Cp*Vdt)*hs;
252 fvm.
source() = -hsTrans()/(this->db().time().deltaT());
262 template<
class CloudType>
271 this->
name() +
":radiation:Ep",
286 const scalar dt = this->db().time().deltaTValue();
288 const scalar
epsilon = constProps_.epsilon0();
289 const scalarField& sumAreaPT4 = radAreaPT4_->field();
298 template<
class CloudType>
307 this->
name() +
":radiation:ap",
322 const scalar dt = this->db().time().deltaTValue();
324 const scalar
epsilon = constProps_.epsilon0();
327 ap = sumAreaP*epsilon/V/dt;
334 template<
class CloudType>
344 this->
name() +
":radiation:sigmap",
359 const scalar dt = this->db().time().deltaTValue();
361 const scalar
epsilon = constProps_.epsilon0();
362 const scalar
f = constProps_.f0();
365 sigmap *= sumAreaP*(1.0 -
f)*(1.0 - epsilon)/V/dt;
372 template<
class CloudType>
398 template<
class CloudType>
fvMatrix< scalar > fvScalarMatrix
tmp< volScalarField > sigmap() const
Return tmp equivalent particulate scattering factor.
Templated heat transfer model class.
DimensionedField< scalar, volMesh > & hsTrans()
Sensible enthalpy transfer [J/kg].
scalar Tmin() const
Minimum temperature.
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
const parcelType::constantProperties & constProps() const
Return the constant properties.
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
Ostream & endl(Ostream &os)
Add newline and flush stream.
const dimensionedScalar sigma
Stefan-Boltzmann constant: default SI units: [W/m2/K4].
const SLGThermo & thermo() const
Return const access to thermo package.
DimensionedField< scalar, volMesh > & radAreaPT4()
Radiation sum of parcel projected area*temperature^4 [m2K4].
const Type & value() const
Return const reference to value.
Templated base class for thermodynamic cloud.
DimensionedField< scalar, volMesh > & radAreaP()
Radiation sum of parcel projected areas [m2].
DimensionedField< scalar, volMesh > & radT4()
Radiation sum of parcel temperature^4 [K4].
tmp< fvMatrix< Type > > SuSp(const DimensionedField< scalar, volMesh > &, const GeometricField< Type, fvPatchField, volMesh > &)
scalar Tmax() const
Maximum temperature.
A special matrix type and solver, designed for finite volume solutions of scalar equations. Face addressing is used to make all matrix assembly and solution loops vectorise.
Thermo package for (S)olids (L)iquids and (G)ases Takes reference to thermo package, and provides:
forAllConstIter(PtrDictionary< phaseModel >, mixture.phases(), phase)
errorManip< error > abort(error &err)
const HeatTransferModel< ThermoCloud< CloudType > > & heatTransfer() const
Return reference to heat transfer model.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
bool radiation() const
Radiation flag.
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
void reduce(const List< UPstream::commsStruct > &comms, T &Value, const BinaryOp &bop, const int tag, const label comm)
word name(const complex &)
Return a string representation of a complex.
const dimensionSet dimEnergy
const volScalarField & T() const
Return const access to the carrier temperature field.
dimensionedScalar pow3(const dimensionedScalar &ds)
const ThermoCloud & cloudCopy() const
Return a reference to the cloud copy.
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
Top level model for Integration schemes.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
const volScalarField & p() const
Return const access to the carrier prressure field.
tmp< fvScalarMatrix > Sh(volScalarField &hs) const
Return sensible enthalpy source term [J/kg/m3/s].
Selector class for relaxation factors, solver type and solution.
const dimensionSet dimMass(1, 0, 0, 0, 0, 0, 0)
A class for managing temporary objects.
T & ref() const
Return non-const reference or generate a fatal error.
DimensionedField< scalar, volMesh > & hsCoeff()
Return coefficient for carrier phase hs equation.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
tmp< volScalarField > ap() const
Return tmp equivalent particulate absorption.
const scalarIntegrationScheme & TIntegrator() const
Return reference to velocity integration.
tmp< volScalarField > Ep() const
Return tmp equivalent particulate emission.