From 7419b0bf98f4dec3dc45a0f6ef804bcadbcb099c Mon Sep 17 00:00:00 2001 From: Andrew Heather <> Date: Mon, 13 Dec 2021 14:37:13 +0000 Subject: [PATCH] ENH: heatTransferCoeff function object - use ::he(p,T) instead of ::he() --- .../heatTransferCoeffModel/heatTransferCoeffModel.C | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) diff --git a/src/functionObjects/field/heatTransferCoeff/heatTransferCoeffModels/heatTransferCoeffModel/heatTransferCoeffModel.C b/src/functionObjects/field/heatTransferCoeff/heatTransferCoeffModels/heatTransferCoeffModel/heatTransferCoeffModel.C index 7b2ffd04b7..92e5965493 100644 --- a/src/functionObjects/field/heatTransferCoeff/heatTransferCoeffModels/heatTransferCoeffModel/heatTransferCoeffModel.C +++ b/src/functionObjects/field/heatTransferCoeff/heatTransferCoeffModels/heatTransferCoeffModel/heatTransferCoeffModel.C @@ -63,7 +63,8 @@ Foam::heatTransferCoeffModel::q() const const auto& turb = mesh_.lookupObject(cmpTurbModel::propertiesName); - const volScalarField& he = turb.transport().he(); + // Note: calling he(p,T) instead of he() + const volScalarField he(turb.transport().he(turb.transport().p(), T)); const volScalarField::Boundary& hebf = he.boundaryField(); const volScalarField alphaEff(turb.alphaEff()); @@ -79,7 +80,8 @@ Foam::heatTransferCoeffModel::q() const const auto& thermo = mesh_.lookupObject(fluidThermo::dictName); - const volScalarField& he = thermo.he(); + // Note: calling he(p,T) instead of he() + const volScalarField he(thermo.he(thermo.p(), T)); const volScalarField::Boundary& hebf = he.boundaryField(); const volScalarField& alpha(thermo.alpha()); @@ -93,7 +95,8 @@ Foam::heatTransferCoeffModel::q() const else { FatalErrorInFunction - << "Unable to find a valid thermo model to evaluate q" << nl + << "Unable to find a valid thermo model to evaluate q. " << nl + << "Database contents are: " << mesh_.objectRegistry::sortedToc() << exit(FatalError); }