zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

wallHeatFlux: Updated to use the thermophysicalTransportModel::q() function

Browse Source 
providing consistency between the reported heat fluxes and those used in the
energy equation.
This commit is contained in:
Henry Weller
2020-09-22 15:30:28 +01:00
parent 77ee78b47e
commit 16da54acda
3 changed files with 15 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files

24
src/functionObjects/field/wallHeatFlux/wallHeatFlux.C
View File

@ -27,7 +27,7 @@ License
#include "thermophysicalTransportModel.H"
#include "solidThermo.H"
#include "surfaceInterpolate.H"
#include "fvcSnGrad.H"
#include "fvcGrad.H"
#include "wallPolyPatch.H"
#include "addToRunTimeSelectionTable.H"
@ -59,11 +59,7 @@ void Foam::functionObjects::wallHeatFlux::writeFileHeader(const label i)
Foam::tmp<Foam::volScalarField>
Foam::functionObjects::wallHeatFlux::calcWallHeatFlux
(
const volScalarField& alpha,
const volScalarField& he
)
Foam::functionObjects::wallHeatFlux::calcWallHeatFlux(const volVectorField& q)
{
tmp<volScalarField> twallHeatFlux
(
@ -78,14 +74,16 @@ Foam::functionObjects::wallHeatFlux::calcWallHeatFlux
volScalarField::Boundary& wallHeatFluxBf =
twallHeatFlux.ref().boundaryFieldRef();
const volScalarField::Boundary& heBf = he.boundaryField();
const volScalarField::Boundary& alphaBf = alpha.boundaryField();
const volVectorField::Boundary& qBf = q.boundaryField();
forAllConstIter(labelHashSet, patchSet_, iter)
{
const label patchi = iter.key();
wallHeatFluxBf[patchi] = alphaBf[patchi]*heBf[patchi].snGrad();
const vectorField& Sfp = mesh_.Sf().boundaryField()[patchi];
const scalarField& magSfp = mesh_.magSf().boundaryField()[patchi];
wallHeatFluxBf[patchi] = (-Sfp/magSfp) & qBf[patchi];
}
if (foundObject<volScalarField>("qr"))
@ -211,7 +209,7 @@ bool Foam::functionObjects::wallHeatFlux::execute()
return store
(
name,
calcWallHeatFlux(ttm.alphaEff(), ttm.thermo().he())
calcWallHeatFlux(ttm.q())
);
}
else if (foundObject<solidThermo>(solidThermo::dictName))
@ -219,7 +217,11 @@ bool Foam::functionObjects::wallHeatFlux::execute()
const solidThermo& thermo =
lookupObject<solidThermo>(solidThermo::dictName);
return store(name, calcWallHeatFlux(thermo.alpha(), thermo.he()));
return store
(
name,
calcWallHeatFlux(-thermo.alpha()*fvc::grad(thermo.he()))
);
}
else
{

6
src/functionObjects/field/wallHeatFlux/wallHeatFlux.H
View File

@ -108,11 +108,7 @@ protected:
virtual void writeFileHeader(const label i);
//- Calculate the heat-flux
tmp<volScalarField> calcWallHeatFlux
(
const volScalarField& alpha,
const volScalarField& he
);
tmp<volScalarField> calcWallHeatFlux(const volVectorField& q);
public:

3
src/functionObjects/field/wallShearStress/wallShearStress.C
View File

@ -82,9 +82,8 @@ Foam::functionObjects::wallShearStress::calcShearStress
const vectorField& Sfp = mesh_.Sf().boundaryField()[patchi];
const scalarField& magSfp = mesh_.magSf().boundaryField()[patchi];
const symmTensorField& taup = tau.boundaryField()[patchi];
wallShearStressBf[patchi] = (-Sfp/magSfp) & taup;
wallShearStressBf[patchi] = (-Sfp/magSfp) & tau.boundaryField()[patchi];
}
return twallShearStress;