exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Bugfix in fixedHeatFluxTemperature BC

Browse source 
1. Reading in gradient from dictionary if present,
2. No need to distinguish between RAS and LES models: use turbulenceModel
   instead.
This commit is contained in:
Vuko Vukcevic 2018-06-11 14:23:48 +02:00
parent 27e5ca1f6b
commit 3d6293352a
1 changed files with 15 additions and 22 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

37
src/turbulenceModels/incompressible/LES/derivedFvPatchFields/fixedHeatFluxTemperature/fixedHeatFluxTemperatureFvPatchScalarField.C
View file

@ -48,9 +48,7 @@ fixedHeatFluxTemperatureFvPatchScalarField
Prt_(0.0),
rhoRef_(0.0),
c_(0.0)
{
this->gradient() = 0.0;
}
{}
fixedHeatFluxTemperatureFvPatchScalarField::
@ -68,8 +66,11 @@ fixedHeatFluxTemperatureFvPatchScalarField
rhoRef_(0.0), // Initialized in constructor body
c_(0.0) // Initialized in constructor body
{
// Set dummy gradient
this->gradient() = 0.0;
// Read the gradient entry from the dictionary
if (dict.found("gradient"))
{
this->gradient() = scalarField("gradient", dict, p.size());
}
// Read the value entry from the dictionary
if (dict.found("value"))
@ -275,28 +276,20 @@ void fixedHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
// Get this patch index
const label patchID = this->patch().index();
if (mesh.foundObject<incompressible::RASModel>("RASProperties"))
if (mesh.foundObject<incompressible::turbulenceModel>("turbulenceModel"))
{
const incompressible::RASModel& ras =
mesh.lookupObject<incompressible::RASModel>("RASProperties");
const incompressible::turbulenceModel& turb =
mesh.lookupObject
<
incompressible::turbulenceModel
>("turbulenceModel");
// Calculate effective kappa at the patch
const scalarField kappaEffp =
ras.nu().boundaryField()[patchID]/Pr_
+ ras.nut()().boundaryField()[patchID]/Prt_;
this->gradient() = heatFlux_/(kappaEffp*rhoRef_*c_);
}
else if (mesh.foundObject<incompressible::LESModel>("LESProperties"))
{
const incompressible::LESModel& les =
mesh.lookupObject<incompressible::LESModel>("LESProperties");
// Calculate effective kappa at the patch
const scalarField kappaEffp =
les.nu().boundaryField()[patchID]/Pr_
+ les.nut()().boundaryField()[patchID]/Prt_;
turb.nu().boundaryField()[patchID]/Pr_
+ turb.nut()().boundaryField()[patchID]/Prt_;
// Calculate gradient at the boundary
this->gradient() = heatFlux_/(kappaEffp*rhoRef_*c_);
}
else