songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

ENH: changing externalWallHeatFluxTemperature BC to work using mixed type.

Browse Source
This commit is contained in:
sergio 2014-04-30 12:00:17 +01:00 committed by Andrew Heather
parent 60c761d0ca
commit 222e403794
2 changed files with 17 additions and 24 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

39
src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/externalWallHeatFluxTemperature/externalWallHeatFluxTemperatureFvPatchScalarField.C
View File

@ -72,9 +72,9 @@ externalWallHeatFluxTemperatureFvPatchScalarField
thicknessLayers_(),
kappaLayers_()
{
this->refValue() = 0.0;
this->refGrad() = 0.0;
this->valueFraction() = 1.0;
refValue() = 0.0;
refGrad() = 0.0;
valueFraction() = 1.0;
}
@ -242,17 +242,13 @@ void Foam::externalWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
return;
}
scalarField q(size(), 0.0);
const scalarField Tc(patchInternalField());
const scalarField Tp(*this);
const scalarField KWall(kappa(Tp));
const scalarField KDelta(KWall*patch().deltaCoeffs());
scalarField hp(patch().size(), 0.0);
switch (mode_)
{
case fixedHeatFlux:
{
q = q_;
break;
}
case fixedHeatTransferCoeff:
@ -269,7 +265,7 @@ void Foam::externalWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
}
}
}
q = (Ta_ - Tp)/(1.0/h_ + totalSolidRes);
hp = 1.0/(1.0/h_ + totalSolidRes);
break;
}
default:
@ -283,27 +279,24 @@ void Foam::externalWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
}
}
forAll(*this, i)
if (mode_ == fixedHeatFlux)
{
if (q[i] > 0) //in
{
this->refGrad()[i] = q[i]/KWall[i];
this->refValue()[i] = 0.0;
this->valueFraction()[i] = 0.0;
}
else //out
{
this->refGrad()[i] = 0.0;
this->refValue()[i] = q[i]/KDelta[i] + Tc[i];
this->valueFraction()[i] = 1.0;
}
refGrad() = q_/kappa(Tp);
refValue() = 0.0;
valueFraction() = 0.0;
}
else if (mode_ == fixedHeatTransferCoeff)
{
refGrad() = 0.0;
refValue() = Ta_;
valueFraction() = hp/(hp + kappa(Tp)*patch().deltaCoeffs());
}
mixedFvPatchScalarField::updateCoeffs();
if (debug)
{
scalar Q = gSum(KWall*patch().magSf()*snGrad());
scalar Q = gSum(kappa(Tp)*patch().magSf()*snGrad());
Info<< patch().boundaryMesh().mesh().name() << ':'
<< patch().name() << ':'

2
src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/externalWallHeatFluxTemperature/externalWallHeatFluxTemperatureFvPatchScalarField.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License