39 namespace functionObjects
49 Foam::functionObjects::thermoCoupleProbes::thermoCoupleProbes
54 const bool loadFromFiles,
60 UName_(
dict.getOrDefault<
word>(
"U",
"U")),
61 radiationFieldName_(
dict.
get<
word>(
"radiationField")),
75 probeDict.readEntry(
"Tc",
Ttc_);
111 if (radiationFieldName_ !=
"none")
120 rhoc = this->sample(thermo_.rho()());
121 kappac = this->sample(thermo_.kappa()());
122 muc = this->sample(thermo_.mu()());
123 Cpc = this->sample(thermo_.Cp()());
150 derivatives(
x,
y, dfdt);
152 const label
n = nEqns();
154 for (label i=0; i<
n; i++)
156 for (label j=0; j<
n; j++)
168 (void) prepare(ACTION_WRITE);
170 const auto& Tfield = thermo_.T();
171 writeValues(Tfield.name(), Ttc_, Tfield.time().timeOutputValue());
174 probeDict.add(
"Tc", Ttc_);
175 setProperty(typeName, probeDict);
187 scalar dt = mesh_.time().deltaTValue();
188 scalar t = mesh_.time().value();
189 odeSolver_->solve(t, t+dt, Ttc_, dt);
201 dict.readEntry(
"rho", rho_);
202 dict.readEntry(
"Cp", Cp_);
203 dict.readEntry(
"d", d_);
204 dict.readEntry(
"epsilon", epsilon_);
Different types of constants.
List< ReturnType > get(const UPtrList< T > &list, const AccessOp &aop)
List of values generated by applying the access operation to each list item.
dimensionedScalar Pr("Pr", dimless, laminarTransport)
virtual bool write()
Sample and write.
Abstract base-class for fluid and solid thermodynamic properties.
Abstract base class for the systems of ordinary differential equations.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
A list of keyword definitions, which are a keyword followed by a number of values (eg...
virtual bool read(const dictionary &)
Read.
const dimensionedScalar G
Newtonian constant of gravitation.
bool getDict(const word &entryName, dictionary &dict) const
Set dictionary, return true if set.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
autoPtr< ODESolver > odeSolver_
ODESolver.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
bool empty() const noexcept
True if List is empty (ie, size() is zero)
const word dictName("faMeshDefinition")
Set of locations to sample.
dimensionedScalar sqrt(const dimensionedScalar &ds)
virtual const volScalarField & T() const
Temperature [K].
const fluidThermo & thermo_
Fluid thermo reference.
dimensionedScalar sigma("sigma", dimMass/sqr(dimTime), transportProperties)
Class to control time during OpenFOAM simulations that is also the top-level objectRegistry.
Macros for easy insertion into run-time selection tables.
virtual label nEqns() const
Number of ODE's to solve.
virtual void derivatives(const scalar x, const scalarField &y, scalarField &dydx) const
Calculate the derivatives in dydx.
virtual bool execute()
Execute. Evaluates the ODESolver.
static const char *const typeName
Typename for Field.
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for expressions::valueTypeCode::INVALID.
void readFields(const typename GeoFieldType::Mesh &mesh, const IOobjectList &objects, const NameMatchPredicate &selectedFields, DynamicList< regIOobject *> &storedObjects)
Read the selected GeometricFields of the templated type and store on the objectRegistry.
const wordList area
Standard area field types (scalar, vector, tensor, etc)
A class for handling words, derived from Foam::string.
Fundamental fluid thermodynamic properties.
constexpr scalar pi(M_PI)
const wordList volume
Standard volume field types (scalar, vector, tensor, etc)
defineTypeNameAndDebug(combustionModel, 0)
scalarField Re(const UList< complex > &cmplx)
Extract real component.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
scalarField Ttc_
Cached thermocouple temperature.
dimensionedScalar pow4(const dimensionedScalar &ds)
virtual bool read(const dictionary &)
Read the probes.
tmp< Field< Type > > sample(const VolumeField< Type > &) const
Sample a volume field at all locations.
virtual void jacobian(const scalar t, const scalarField &y, scalarField &dfdt, scalarSquareMatrix &dfdy) const
Calculate the Jacobian of the system.
static autoPtr< ODESolver > New(const ODESystem &ode, const dictionary &dict)
Select null constructed.
addToRunTimeSelectionTable(functionObject, pointHistory, dictionary)
static constexpr const zero Zero
Global zero (0)