zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: Renaming K to kappa and kappa to kappaRad in solid thermo and

Browse Source 
updating tutorials
This commit is contained in:
sergio
2012-05-08 11:02:39 +01:00
parent 060a2b0f42
commit fa4fd968a8
87 changed files with 430 additions and 417 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousSolid/setPorousRegionSolidFields.H
View File

@ -9,8 +9,8 @@
tmp<volScalarField> tcp = thermo.Cp();
const volScalarField& cp = tcp();
tmp<volScalarField> tkappa = thermo.K();
//tmp<volSymmTensorField> tkappa = thermo.directionalK()*betav;
tmp<volScalarField> tkappa = thermo.kappa();
//tmp<volSymmTensorField> tkappa = thermo.directionalKappa()*betav;
const volScalarField& kappa = tkappa();
//const volSymmTensorField& K = tK();

2
applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/solid/setRegionSolidFields.H
View File

@ -7,7 +7,7 @@
tmp<volScalarField> tcp = thermo.Cp();
const volScalarField& cp = tcp();
tmp<volScalarField> tkappa = thermo.K();
tmp<volScalarField> tkappa = thermo.kappa();
//tmp<volSymmTensorField> tkappa = thermo.directionalkappa();
const volScalarField& kappa = tkappa();

4
applications/solvers/heatTransfer/chtMultiRegionFoam/porousSolid/setPorousRegionSolidFields.H
View File

@ -9,8 +9,8 @@
tmp<volScalarField> tcp = thermo.Cp();
const volScalarField& cp = tcp();
tmp<volScalarField> tkappa = thermo.K();
//tmp<volSymmTensorField> tkappa = thermo.directionalK()*betav;
tmp<volScalarField> tkappa = thermo.kappa();
//tmp<volSymmTensorField> tkappa = thermo.directionalKappa()*betav;
const volScalarField& kappa = tkappa();
//const volSymmTensorField& K = tK();

4
applications/solvers/heatTransfer/chtMultiRegionFoam/solid/setRegionSolidFields.H
View File

@ -7,9 +7,9 @@
tmp<volScalarField> tcp = thermo.Cp();
const volScalarField& cp = tcp();
tmp<volScalarField> tkappa = thermo.K();
tmp<volScalarField> tkappa = thermo.kappa();
const volScalarField& kappa = tkappa();
//tmp<volSymmTensorField> tkappa = thermo.directionalK();
//tmp<volSymmTensorField> tkappa = thermo.directionalKappa();
//const volSymmTensorField& kappa = tkappa();
volScalarField& T = thermo.T();

2
applications/utilities/postProcessing/wall/solidWallHeatFlux/solidWallHeatFlux.C
View File

@ -73,7 +73,7 @@ int main(int argc, char *argv[])
forAll(wallHeatFlux.boundaryField(), patchi)
{
wallHeatFlux.boundaryField()[patchi] =
thermo().K(patchi)
thermo().kappa(patchi)
*T.boundaryField()[patchi].snGrad();
if (isA<wallFvPatch>(mesh.boundary()[patchi]))

8
src/regionModels/pyrolysisModels/noPyrolysis/noPyrolysis.C
View File

@ -151,15 +151,15 @@ const tmp<volScalarField> noPyrolysis::Cp() const
}
const volScalarField& noPyrolysis::kappa() const
const volScalarField& noPyrolysis::kappaRad() const
{
return (solidThermo_->kappa());
return (solidThermo_->kappaRad());
}
const volScalarField& noPyrolysis::K() const
const volScalarField& noPyrolysis::kappa() const
{
return (solidThermo_->K());
return (solidThermo_->kappa());
}

4
src/regionModels/pyrolysisModels/noPyrolysis/noPyrolysis.H
View File

@ -124,10 +124,10 @@ public:
virtual const tmp<volScalarField> Cp() const;
//- Return the region absorptivity [1/m]
virtual const volScalarField& kappa() const;
virtual const volScalarField& kappaRad() const;
//- Return the region thermal conductivity [W/m/k]
virtual const volScalarField& K() const;
virtual const volScalarField& kappa() const;
//- Return the total gas mass flux to primary region [kg/m2/s]
virtual const surfaceScalarField& phiGas() const;

4
src/regionModels/pyrolysisModels/pyrolysisModel/pyrolysisModel.H
View File

@ -218,10 +218,10 @@ public:
virtual const tmp<volScalarField> Cp() const = 0;
//- Return the region absorptivity [1/m]
virtual const volScalarField& kappa() const = 0;
virtual const volScalarField& kappaRad() const = 0;
//- Return the region thermal conductivity [W/m/k]
virtual const volScalarField& K() const = 0;
virtual const volScalarField& kappa() const = 0;
//- Return the total gas mass flux to primary region [kg/m2/s]
virtual const surfaceScalarField& phiGas() const = 0;

18
src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C
View File

@ -126,7 +126,7 @@ void reactingOneDim::updateQr()
const label cellI = cells[k];
const point& Cf1 = regionMesh().cellCentres()[cellI];
const scalar delta = mag(Cf1 - Cf0);
kappaInt += kappa_[cellI]*delta;
kappaInt += kappaRad_[cellI]*delta;
Qr_[cellI] = Qr0*exp(-kappaInt);
Cf0 = Cf1;
}
@ -294,7 +294,7 @@ void reactingOneDim::solveEnergy()
fvScalarMatrix TEqn
(
fvm::ddt(rhoCp, T_)
- fvm::laplacian(K_, T_)
- fvm::laplacian(kappa_, T_)
==
chemistrySh_
+ fvc::div(phiQr)
@ -347,8 +347,8 @@ reactingOneDim::reactingOneDim(const word& modelType, const fvMesh& mesh)
pyrolysisModel(modelType, mesh),
solidChemistry_(solidChemistryModel::New(regionMesh())),
solidThermo_(solidChemistry_->solidThermo()),
kappaRad_(solidThermo_.kappaRad()),
kappa_(solidThermo_.kappa()),
K_(solidThermo_.K()),
rho_(solidThermo_.rho()),
Ys_(solidThermo_.composition().Y()),
T_(solidThermo_.T()),
@ -449,8 +449,8 @@ reactingOneDim::reactingOneDim
pyrolysisModel(modelType, mesh, dict),
solidChemistry_(solidChemistryModel::New(regionMesh())),
solidThermo_(solidChemistry_->solidThermo()),
kappaRad_(solidThermo_.kappaRad()),
kappa_(solidThermo_.kappa()),
K_(solidThermo_.K()),
rho_(solidThermo_.rho()),
Ys_(solidThermo_.composition().Y()),
T_(solidThermo_.T()),
@ -584,7 +584,7 @@ scalar reactingOneDim::solidRegionDiffNo() const
surfaceScalarField KrhoCpbyDelta
(
regionMesh().surfaceInterpolation::deltaCoeffs()
* fvc::interpolate(K_)
* fvc::interpolate(kappa_)
/ fvc::interpolate(Cp()*rho_)
);
@ -619,15 +619,15 @@ const tmp<volScalarField> reactingOneDim::Cp() const
}
const volScalarField& reactingOneDim::kappa() const
const volScalarField& reactingOneDim::kappaRad() const
{
return kappa_;
return kappaRad_;
}
const volScalarField& reactingOneDim::K() const
const volScalarField& reactingOneDim::kappa() const
{
return K_;
return kappa_;
}

8
src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.H
View File

@ -84,10 +84,10 @@ protected:
// Reference to solid thermo properties
//- Absorption coefficient [1/m]
const volScalarField& kappa_;
const volScalarField& kappaRad_;
//- Thermal conductivity [W/m/K]
const volScalarField& K_;
const volScalarField& kappa_;
//- Density [kg/m3]
volScalarField& rho_;
@ -230,10 +230,10 @@ public:
virtual const tmp<volScalarField> Cp() const;
//- Return the region absorptivity [1/m]
virtual const volScalarField& kappa() const;
virtual const volScalarField& kappaRad() const;
//- Return the region thermal conductivity [W/m/k]
virtual const volScalarField& K() const;
virtual const volScalarField& kappa() const;
//- Return the total gas mass flux to primary region [kg/m2/s]
virtual const surfaceScalarField& phiGas() const;

6
src/regionModels/thermoBaffleModels/derivedFvPatchFields/temperatureThermoBaffle/temperatureThermoBaffleFvPatchScalarField.H
View File

@ -35,7 +35,7 @@ Description
// Coupled BC.
neighbourFieldName T;
K basicThermo;
kappa basicThermo;
KName none;
@ -57,10 +57,10 @@ Description
//- thermo properties
rho rho [1 -3 0 0 0 0 0] 80;
Cp Cp [0 2 -2 -1 0 0 0] 15;
K K [1 1 -3 -1 0 0 0] 0.01;
kappa kappa [1 1 -3 -1 0 0 0] 0.01;
//- radiative properties
kappa kappa [0 -1 0 0 0 0 0] 0;
kappaRad kappaRad [0 -1 0 0 0 0 0] 0;
sigmaS sigmaS [0 -1 0 0 0 0 0] 0;
emissivity emissivity [0 0 0 0 0 0 0] 1;

8
src/regionModels/thermoBaffleModels/noThermo/noThermo.C
View File

@ -95,9 +95,9 @@ const tmp<volScalarField> noThermo::Cp() const
);
}
const volScalarField& noThermo::kappa() const
const volScalarField& noThermo::kappaRad() const
{
FatalErrorIn("const volScalarField& noThermo::kappa() const")
FatalErrorIn("const volScalarField& noThermo::kappaRad() const")
<< "kappa field not available for " << type() << abort(FatalError);
return volScalarField::null();
}
@ -111,9 +111,9 @@ const volScalarField& noThermo::rho() const
}
const volScalarField& noThermo::K() const
const volScalarField& noThermo::kappa() const
{
FatalErrorIn("const volScalarField& noThermo::K() const")
FatalErrorIn("const volScalarField& noThermo::kappa() const")
<< "K field not available for " << type() << abort(FatalError);
return volScalarField::null();
}

4
src/regionModels/thermoBaffleModels/noThermo/noThermo.H
View File

@ -105,7 +105,7 @@ public:
virtual const tmp<volScalarField> Cp() const;
//- Return solid absortivity [1/m]
virtual const volScalarField& kappa() const;
virtual const volScalarField& kappaRad() const;
//- Return the film mean temperature [K]
virtual const volScalarField& T() const;
@ -114,7 +114,7 @@ public:
virtual const volScalarField& rho() const;
//- Return thermal conductivity [W/m/K]
virtual const volScalarField& K() const;
virtual const volScalarField& kappa() const;
// Evolution

16
src/regionModels/thermoBaffleModels/thermoBaffle2D/thermoBaffle2D.C
View File

@ -96,7 +96,7 @@ void thermoBaffle2D::solveEnergy()
volScalarField& Q = tQ();
volScalarField rhoCp("rhoCp", thermo_->rho()*thermo_->Cp()());
volScalarField K("K", thermo_->K());
volScalarField kappa("kappa", thermo_->kappa());
//If region is one-dimension variable thickness
@ -119,7 +119,7 @@ void thermoBaffle2D::solveEnergy()
rhoCp[cellId] *= delta_.value()/thickness_[localFaceI];
K[cellId] *= delta_.value()/thickness_[localFaceI];
kappa[cellId] *= delta_.value()/thickness_[localFaceI];
}
}
}
@ -131,7 +131,7 @@ void thermoBaffle2D::solveEnergy()
fvScalarMatrix TEqn
(
fvm::ddt(rhoCp, T_)
- fvm::laplacian(K, T_)
- fvm::laplacian(kappa, T_)
==
Q
);
@ -316,9 +316,9 @@ const tmp<volScalarField> thermoBaffle2D::Cp() const
}
const volScalarField& thermoBaffle2D::kappa() const
const volScalarField& thermoBaffle2D::kappaRad() const
{
return thermo_->kappa();
return thermo_->kappaRad();
}
@ -328,9 +328,9 @@ const volScalarField& thermoBaffle2D::rho() const
}
const volScalarField& thermoBaffle2D::K() const
const volScalarField& thermoBaffle2D::kappa() const
{
return thermo_->K();
return thermo_->kappa();
}
@ -362,7 +362,7 @@ void thermoBaffle2D::info() const
(
mag(regionMesh().Sf().boundaryField()[patchI])
* pT.snGrad()
* thermo_->K(patchI)
* thermo_->kappa(patchI)
) << endl;
}
}

4
src/regionModels/thermoBaffleModels/thermoBaffle2D/thermoBaffle2D.H
View File

@ -155,7 +155,7 @@ public:
virtual const tmp<volScalarField> Cp() const;
//- Return solid absortivity [1/m]
virtual const volScalarField& kappa() const;
virtual const volScalarField& kappaRad() const;
//- Return the film mean temperature [K]
virtual const volScalarField& T() const;
@ -164,7 +164,7 @@ public:
virtual const volScalarField& rho() const;
//- Return thermal conductivity [W/m/K]
virtual const volScalarField& K() const;
virtual const volScalarField& kappa() const;
// Helper functions

4
src/regionModels/thermoBaffleModels/thermoBaffleModel/thermoBaffleModel.H
View File

@ -211,10 +211,10 @@ public:
virtual const tmp<volScalarField> Cp() const = 0;
//- Return the region absorptivity [1/m]
virtual const volScalarField& kappa() const = 0;
virtual const volScalarField& kappaRad() const = 0;
//- Return the region thermal conductivity [W/m/k]
virtual const volScalarField& K() const = 0;
virtual const volScalarField& kappa() const = 0;
// Evolution

12
src/thermophysicalModels/basicSolidThermo/basicSolidThermo/basicSolidThermo.C
View File

@ -78,11 +78,11 @@ Foam::basicSolidThermo::basicSolidThermo(const fvMesh& mesh)
mesh,
dimMass/dimVolume
),
kappa_
kappaRad_
(
IOobject
(
"kappa",
"kappaRad",
mesh.time().timeName(),
mesh,
IOobject::NO_READ,
@ -165,11 +165,11 @@ Foam::basicSolidThermo::basicSolidThermo
mesh,
dimMass/dimVolume
),
kappa_
kappaRad_
(
IOobject
(
"kappa",
"kappaRad",
mesh.time().timeName(),
mesh,
IOobject::NO_READ,
@ -239,9 +239,9 @@ Foam::volScalarField& Foam::basicSolidThermo::rho()
}
const Foam::volScalarField& Foam::basicSolidThermo::kappa() const
const Foam::volScalarField& Foam::basicSolidThermo::kappaRad() const
{
return kappa_;
return kappaRad_;
}

12
src/thermophysicalModels/basicSolidThermo/basicSolidThermo/basicSolidThermo.H
View File

@ -73,7 +73,7 @@ protected:
volScalarField rho_;
//- Absorption coefficient [1/m]
volScalarField kappa_;
volScalarField kappaRad_;
//- Scatter coeffcient [1/m]
volScalarField sigmaS_;
@ -157,7 +157,7 @@ public:
virtual volScalarField& rho();
//- Absorption coefficient [1/m]
virtual const volScalarField& kappa() const;
virtual const volScalarField& kappaRad() const;
//- Emissivity coefficient
virtual const volScalarField& sigmaS() const;
@ -169,10 +169,10 @@ public:
// Derived thermal properties
//- Thermal conductivity [W/m/K]
virtual tmp<volScalarField> K() const = 0;
virtual tmp<volScalarField> kappa() const = 0;
//- Thermal conductivity [W/m/K]
virtual tmp<volSymmTensorField> directionalK() const = 0;
virtual tmp<volSymmTensorField> directionalKappa() const = 0;
//- Specific heat capacity [J/kg/K]
virtual tmp<volScalarField> Cp() const = 0;
@ -187,10 +187,10 @@ public:
// Per patch calculation
//- Thermal conductivity [W//m/K]
virtual tmp<scalarField> K(const label patchI) const = 0;
virtual tmp<scalarField> kappa(const label patchI) const = 0;
//- Thermal conductivity [W//m/K]
virtual tmp<symmTensorField> directionalK(const label) const = 0;
virtual tmp<symmTensorField> directionalKappa(const label) const = 0;
//- Specific heat capacity [J/kg/K)]
virtual tmp<scalarField> Cp(const label patchI) const = 0;

70
src/thermophysicalModels/basicSolidThermo/constSolidThermo/constSolidThermo.C
View File

@ -46,36 +46,36 @@ Foam::constSolidThermo::constSolidThermo
:
basicSolidThermo(mesh, dict),
dict_(dict.subDict(typeName + "Coeffs")),
constK_(dimensionedScalar(dict_.lookup("K"))),
K_
constKappa_(dimensionedScalar(dict_.lookup("kappa"))),
kappa_
(
IOobject
(
"K",
"kappa",
mesh.time().timeName(),
mesh,
IOobject::NO_READ,
IOobject::NO_WRITE
),
mesh,
constK_
constKappa_
),
constRho_(dimensionedScalar(dict_.lookup("rho"))),
constCp_(dimensionedScalar(dict_.lookup("Cp"))),
constHf_(dimensionedScalar(dict_.lookup("Hf"))),
constEmissivity_(dimensionedScalar(dict_.lookup("emissivity"))),
constKappa_(dimensionedScalar(dict_.lookup("kappa"))),
constKappaRad_(dimensionedScalar(dict_.lookup("kappaRad"))),
constSigmaS_(dimensionedScalar(dict_.lookup("sigmaS")))
{
read();
K_ = constK_;
kappa_ = constKappa_;
rho_ = constRho_;
emissivity_ = constEmissivity_;
kappa_ = constKappa_;
kappaRad_ = constKappaRad_;
sigmaS_ = constSigmaS_;
}
@ -85,36 +85,36 @@ Foam::constSolidThermo::constSolidThermo(const fvMesh& mesh)
:
basicSolidThermo(mesh),
dict_(subDict(typeName + "Coeffs")),
constK_(dimensionedScalar(dict_.lookup("K"))),
K_
constKappa_(dimensionedScalar(dict_.lookup("kappa"))),
kappa_
(
IOobject
(
"K",
"kappa",
mesh.time().timeName(),
mesh,
IOobject::NO_READ,
IOobject::NO_WRITE
),
mesh,
constK_
constKappa_
),
constRho_(dimensionedScalar(dict_.lookup("rho"))),
constCp_(dimensionedScalar(dict_.lookup("Cp"))),
constHf_(dimensionedScalar(dict_.lookup("Hf"))),
constEmissivity_(dimensionedScalar(dict_.lookup("emissivity"))),
constKappa_(dimensionedScalar(dict_.lookup("kappa"))),
constKappaRad_(dimensionedScalar(dict_.lookup("kappaRad"))),
constSigmaS_(dimensionedScalar(dict_.lookup("sigmaS")))
{
read();
K_ = constK_;
kappa_ = constKappa_;
rho_ = constRho_;
emissivity_ = constEmissivity_;
kappa_ = constKappa_;
kappaRad_ = constKappaRad_;
sigmaS_ = constSigmaS_;
}
@ -132,26 +132,26 @@ void Foam::constSolidThermo::correct()
{}
Foam::tmp<Foam::volScalarField> Foam::constSolidThermo::K() const
Foam::tmp<Foam::volScalarField> Foam::constSolidThermo::kappa() const
{
return K_;
return kappa_;
}
Foam::tmp<Foam::volSymmTensorField> Foam::constSolidThermo::directionalK() const
Foam::tmp<Foam::volSymmTensorField> Foam::constSolidThermo::directionalKappa() const
{
dimensionedSymmTensor t
(
constK_.name(),
constK_.dimensions(),
constKappa_.name(),
constKappa_.dimensions(),
symmTensor
(
constK_.value(),
constKappa_.value(),
0.0,
0.0,
constK_.value(),
constKappa_.value(),
0.0,
constK_.value()
constKappa_.value()
)
);
return tmp<volSymmTensorField>
@ -160,7 +160,7 @@ Foam::tmp<Foam::volSymmTensorField> Foam::constSolidThermo::directionalK() const
(
IOobject
(
"K",
"kappa",
mesh_.time().timeName(),
mesh_,
IOobject::NO_READ,
@ -215,28 +215,28 @@ Foam::tmp<Foam::volScalarField> Foam::constSolidThermo::Hf() const
}
Foam::tmp<Foam::scalarField> Foam::constSolidThermo::K
Foam::tmp<Foam::scalarField> Foam::constSolidThermo::kappa
(
const label patchI
) const
{
return (K_.boundaryField()[patchI]);
return (kappa_.boundaryField()[patchI]);
}
Foam::tmp<Foam::symmTensorField> Foam::constSolidThermo::directionalK
Foam::tmp<Foam::symmTensorField> Foam::constSolidThermo::directionalKappa
(
const label patchI
) const
{
symmTensor t
(
constK_.value(),
constKappa_.value(),
0.0,
0.0,
constK_.value(),
constKappa_.value(),
0.0,
constK_.value()
constKappa_.value()
);
return tmp<symmTensorField>
(
@ -291,19 +291,19 @@ bool Foam::constSolidThermo::read(const dictionary& dict)
{
constRho_ = dimensionedScalar(dict.lookup("rho"));
constCp_ = dimensionedScalar(dict.lookup("Cp"));
constK_ = dimensionedScalar(dict.lookup("K"));
constKappa_ = dimensionedScalar(dict.lookup("kappa"));
constHf_ = dimensionedScalar(dict.lookup("Hf"));
constEmissivity_ = dimensionedScalar(dict.lookup("emissivity"));
constKappa_ = dimensionedScalar(dict_.lookup("kappa"));
constKappaRad_ = dimensionedScalar(dict_.lookup("kappaRad"));
constSigmaS_ = dimensionedScalar(dict_.lookup("sigmaS"));
Info<< "Constructed constSolidThermo with" << nl
<< " rho : " << constRho_ << nl
<< " Cp : " << constCp_ << nl
<< " K : " << constK_ << nl
<< " kappa : " << constKappa_ << nl
<< " Hf : " << constHf_ << nl
<< " emissivity : " << constEmissivity_ << nl
<< " kappa : " << constKappa_ << nl
<< " kappaRad : " << constKappaRad_ << nl
<< " sigmaS : " << constSigmaS_ << nl
<< endl;
@ -316,9 +316,9 @@ bool Foam::constSolidThermo::writeData(Ostream& os) const
bool ok = basicSolidThermo::writeData(os);
os.writeKeyword("rho") << constRho_ << token::END_STATEMENT << nl;
os.writeKeyword("Cp") << constCp_ << token::END_STATEMENT << nl;
os.writeKeyword("K") << constK_ << token::END_STATEMENT << nl;
os.writeKeyword("Hf") << constHf_ << token::END_STATEMENT << nl;
os.writeKeyword("kappa") << constKappa_ << token::END_STATEMENT << nl;
os.writeKeyword("Hf") << constHf_ << token::END_STATEMENT << nl;
os.writeKeyword("kappaRad") << constKappaRad_ << token::END_STATEMENT << nl;
os.writeKeyword("sigmaS") << constSigmaS_ << token::END_STATEMENT << nl;
os.writeKeyword("emissivity") << constEmissivity_ << token::END_STATEMENT
<< nl;

14
src/thermophysicalModels/basicSolidThermo/constSolidThermo/constSolidThermo.H
View File

@ -58,10 +58,10 @@ private:
dictionary dict_;
//- Constant thermal conductivity [W/(m.K)]
dimensionedScalar constK_;
dimensionedScalar constKappa_;
//- Thermal conductivity field[W/(m.K)]
volScalarField K_;
volScalarField kappa_;
//- Density [kg/m3]
dimensionedScalar constRho_;
@ -76,7 +76,7 @@ private:
dimensionedScalar constEmissivity_;
//- Absorptivity [1/m]
dimensionedScalar constKappa_;
dimensionedScalar constKappaRad_;
//- Scatter [1/m]
dimensionedScalar constSigmaS_;
@ -110,10 +110,10 @@ public:
// Derived thermal properties
//- Thermal conductivity [W/m/K]
virtual tmp<volScalarField> K() const;
virtual tmp<volScalarField> kappa() const;
//- Thermal conductivity [W/m/K]
virtual tmp<volSymmTensorField> directionalK() const;
virtual tmp<volSymmTensorField> directionalKappa() const;
//- Specific heat capacity [J/(kg.K)]
virtual tmp<volScalarField> Cp() const;
@ -125,10 +125,10 @@ public:
// Per patch calculation
//- Thermal conductivity [W//m/K]
virtual tmp<scalarField> K(const label patchI) const;
virtual tmp<scalarField> kappa(const label patchI) const;
//- Thermal conductivity [W//m/K]
virtual tmp<symmTensorField>directionalK(const label) const;
virtual tmp<symmTensorField>directionalKappa(const label) const;
//- Specific heat capacity [J/kg/K)]
virtual tmp<scalarField> Cp(const label patchI) const;

40
src/thermophysicalModels/basicSolidThermo/directionalKSolidThermo/directionalKSolidThermo.C
View File

@ -59,11 +59,11 @@ Foam::directionalKSolidThermo::directionalKSolidThermo
)
:
interpolatedSolidThermo(mesh, typeName + "Coeffs", dict),
directionalK_
directionalKappa_
(
IOobject
(
"K",
"kappa",
mesh_.time().timeName(),
mesh_,
IOobject::NO_READ,
@ -93,11 +93,11 @@ Foam::directionalKSolidThermo::directionalKSolidThermo
Foam::directionalKSolidThermo::directionalKSolidThermo(const fvMesh& mesh)
:
interpolatedSolidThermo(mesh, typeName + "Coeffs"),
directionalK_
directionalKappa_
(
IOobject
(
"K",
"kappa",
mesh_.time().timeName(),
mesh_,
IOobject::NO_READ,
@ -134,7 +134,8 @@ Foam::directionalKSolidThermo::~directionalKSolidThermo()
void Foam::directionalKSolidThermo::init()
{
KValues_ = Field<vector>(subDict(typeName + "Coeffs").lookup("KValues"));
kappaValues_ =
Field<vector>(subDict(typeName + "Coeffs").lookup("kappaValues"));
// Determine transforms for cell centres
forAll(mesh_.C(), cellI)
@ -362,51 +363,51 @@ void Foam::directionalKSolidThermo::correct()
Foam::tmp<Foam::volSymmTensorField>
Foam::directionalKSolidThermo::directionalK() const
Foam::directionalKSolidThermo::directionalKappa() const
{
return directionalK_;
return directionalKappa_;
}
void Foam::directionalKSolidThermo::calculate()
{
// Correct directionalK
Field<vector> localK
Field<vector> localKappa
(
interpolateXY
(
T_.internalField(),
TValues_,
KValues_
kappaValues_
)
);
// Transform into global coordinate system
transformField
(
directionalK_.internalField(),
directionalKappa_.internalField(),
ccTransforms_.internalField(),
localK
localKappa
);
forAll(directionalK_.boundaryField(), patchI)
forAll(directionalKappa_.boundaryField(), patchI)
{
directionalK_.boundaryField()[patchI] == this->directionalK(patchI)();
directionalKappa_.boundaryField()[patchI] == this->directionalKappa(patchI)();
}
}
Foam::tmp<Foam::symmTensorField> Foam::directionalKSolidThermo::directionalK
Foam::tmp<Foam::symmTensorField> Foam::directionalKSolidThermo::directionalKappa
(
const label patchI
) const
{
const fvPatchScalarField& patchT = T_.boundaryField()[patchI];
Field<vector> localK(interpolateXY(patchT, TValues_, KValues_));
Field<vector> localKappa(interpolateXY(patchT, TValues_, kappaValues_));
tmp<symmTensorField> tglobalK(new symmTensorField(localK.size()));
transformField(tglobalK(), ccTransforms_.boundaryField()[patchI], localK);
tmp<symmTensorField> tglobalK(new symmTensorField(localKappa.size()));
transformField(tglobalK(), ccTransforms_.boundaryField()[patchI], localKappa);
return tglobalK;
}
@ -421,7 +422,7 @@ bool Foam::directionalKSolidThermo::read()
bool Foam::directionalKSolidThermo::read(const dictionary& dict)
{
coordSys_ = coordinateSystem(dict, mesh_);
KValues_ = Field<vector>(subDict(typeName + "Coeffs").lookup("KValues"));
kappaValues_ = Field<vector>(subDict(typeName + "Coeffs").lookup("kappaValues"));
return true;
}
@ -429,7 +430,8 @@ bool Foam::directionalKSolidThermo::read(const dictionary& dict)
bool Foam::directionalKSolidThermo::writeData(Ostream& os) const
{
bool ok = interpolatedSolidThermo::writeData(os);
os.writeKeyword("KValues") << KValues_ << token::END_STATEMENT << nl;
os.writeKeyword("kappaValues")
<< kappaValues_ << token::END_STATEMENT << nl;
return ok && os.good();
}

16
src/thermophysicalModels/basicSolidThermo/directionalKSolidThermo/directionalKSolidThermo.H
View File

@ -54,10 +54,10 @@ class directionalKSolidThermo
// Private data
//- Thermal conductivity [W/(m.K)]
volSymmTensorField directionalK_;
volSymmTensorField directionalKappa_;
//- Thermal conductivity vector
Field<vector> KValues_;
Field<vector> kappaValues_;
//- Coordinate system used for the directional properties
coordinateSystem coordSys_;
@ -117,26 +117,26 @@ public:
//- Access functions
//- Thermal conductivity [W/m/K]
virtual tmp<volScalarField> K() const
virtual tmp<volScalarField> kappa() const
{
notImplemented("directionalKSolidThermo::K() const");
notImplemented("directionalKSolidThermo::kappa() const");
return tmp<volScalarField>(NULL);
}
//- Thermal conductivity [W/m/K]
virtual tmp<volSymmTensorField> directionalK() const;
virtual tmp<volSymmTensorField> directionalKappa() const;
// Per patch calculation
//- Thermal conductivity [W//m/K]
virtual tmp<scalarField> K(const label patchI) const
virtual tmp<scalarField> kappa(const label patchI) const
{
notImplemented("directionalKSolidThermo::K(const label) const");
notImplemented("directionalKSolidThermo::kappa(const label) const");
return tmp<scalarField>(NULL);
}
//- Thermal conductivity [W//m/K]
virtual tmp<symmTensorField> directionalK(const label) const;
virtual tmp<symmTensorField> directionalKappa(const label) const;
// I-O

6
src/thermophysicalModels/basicSolidThermo/interpolatedSolidThermo/interpolateSolid/interpolateSolid.C
View File

@ -39,7 +39,7 @@ Foam::interpolateSolid::interpolateSolid(const dictionary& dict)
<< " cp : " << cpValues_ << nl
<< " Hf : " << HfValues_ << nl
<< " emissivity : " << emissivityValues_ << nl
<< " kappa : " << kappaValues_ << nl
<< " kappaRad : " << kappaRadValues_ << nl
<< " sigmaS : " << sigmaSValues_ << nl
<< endl;
@ -88,7 +88,7 @@ bool Foam::interpolateSolid::writeData(Ostream& os) const
os.writeKeyword("HfValues") << HfValues_ << token::END_STATEMENT << nl;
os.writeKeyword("emissivityValues") << emissivityValues_ <<
token::END_STATEMENT << nl;
os.writeKeyword("kappaValues") << kappaValues_
os.writeKeyword("kappaRadValues") << kappaRadValues_
<< token::END_STATEMENT << nl;
os.writeKeyword("sigmaSValues") << sigmaSValues_
<< token::END_STATEMENT << nl;
@ -103,7 +103,7 @@ bool Foam::interpolateSolid::read(const dictionary& dict)
TValues_ = Field<scalar>(dict.lookup("TValues"));
rhoValues_ = Field<scalar>(dict.lookup("rhoValues"));
cpValues_ = Field<scalar>(dict.lookup("cpValues"));
kappaValues_ = Field<scalar>(dict.lookup("kappaValues"));
kappaRadValues_ = Field<scalar>(dict.lookup("kappaRadValues"));
sigmaSValues_ = Field<scalar>(dict.lookup("sigmaSValues"));
HfValues_ = Field<scalar>(dict.lookup("HfValues"));
emissivityValues_ = Field<scalar>(dict.lookup("emissivityValues"));

2
src/thermophysicalModels/basicSolidThermo/interpolatedSolidThermo/interpolateSolid/interpolateSolid.H
View File

@ -63,7 +63,7 @@ protected:
Field<scalar> emissivityValues_;
Field<scalar> kappaValues_;
Field<scalar> kappaRadValues_;
Field<scalar> sigmaSValues_;

12
src/thermophysicalModels/basicSolidThermo/interpolatedSolidThermo/interpolatedSolidThermo.C
View File

@ -69,7 +69,7 @@ Foam::tmp<Foam::scalarField> Foam::interpolatedSolidThermo::emissivity
}
Foam::tmp<Foam::scalarField> Foam::interpolatedSolidThermo::kappa
Foam::tmp<Foam::scalarField> Foam::interpolatedSolidThermo::kappaRad
(
const label patchI
) const
@ -82,7 +82,7 @@ Foam::tmp<Foam::scalarField> Foam::interpolatedSolidThermo::kappa
(
T_.boundaryField()[patchI],
TValues_,
kappaValues_
kappaRadValues_
)
)
);
@ -178,16 +178,16 @@ void Foam::interpolatedSolidThermo::calculate()
// Correct absorptivity
kappa_.internalField() = interpolateXY
kappaRad_.internalField() = interpolateXY
(
T_.internalField(),
TValues_,
kappaValues_
kappaRadValues_
);
forAll(kappa_.boundaryField(), patchI)
forAll(kappaRad_.boundaryField(), patchI)
{
kappa_.boundaryField()[patchI] == this->kappa(patchI)();
kappaRad_.boundaryField()[patchI] == this->kappaRad(patchI)();
}

2
src/thermophysicalModels/basicSolidThermo/interpolatedSolidThermo/interpolatedSolidThermo.H
View File

@ -69,7 +69,7 @@ private:
tmp<scalarField> sigmaS(const label) const;
//- Absorption coefficient [1/m]
tmp<scalarField> kappa(const label) const;
tmp<scalarField> kappaRad(const label) const;
//- Emissivity []
tmp<scalarField> emissivity(const label) const;

40
src/thermophysicalModels/basicSolidThermo/isotropicKSolidThermo/isotropicKSolidThermo.C
View File

@ -57,11 +57,11 @@ Foam::isotropicKSolidThermo::isotropicKSolidThermo
)
:
interpolatedSolidThermo(mesh, typeName + "Coeffs", dict),
K_
kappa_
(
IOobject
(
"K",
"kappa",
mesh.time().timeName(),
mesh,
IOobject::NO_READ,
@ -70,7 +70,10 @@ Foam::isotropicKSolidThermo::isotropicKSolidThermo
mesh,
dimEnergy/dimTime/(dimLength*dimTemperature)
),
KValues_ (Field<scalar>(subDict(typeName + "Coeffs").lookup("KValues")))
kappaValues_
(
Field<scalar>(subDict(typeName + "Coeffs").lookup("kappaValues"))
)
{
correct();
}
@ -79,11 +82,11 @@ Foam::isotropicKSolidThermo::isotropicKSolidThermo
Foam::isotropicKSolidThermo::isotropicKSolidThermo(const fvMesh& mesh)
:
interpolatedSolidThermo(mesh, typeName + "Coeffs"),
K_
kappa_
(
IOobject
(
"K",
"kappa",
mesh.time().timeName(),
mesh,
IOobject::NO_READ,
@ -92,7 +95,10 @@ Foam::isotropicKSolidThermo::isotropicKSolidThermo(const fvMesh& mesh)
mesh,
dimEnergy/dimTime/(dimLength*dimTemperature)
),
KValues_ (Field<scalar>(subDict(typeName + "Coeffs").lookup("KValues")))
kappaValues_
(
Field<scalar>(subDict(typeName + "Coeffs").lookup("kappaValues"))
)
{
correct();
}
@ -101,20 +107,20 @@ Foam::isotropicKSolidThermo::isotropicKSolidThermo(const fvMesh& mesh)
void Foam::isotropicKSolidThermo::correct()
{
// Correct K
K_.internalField() = interpolateXY
kappa_.internalField() = interpolateXY
(
T_.internalField(),
TValues_,
KValues_
kappaValues_
);
forAll(K_.boundaryField(), patchI)
forAll(kappa_.boundaryField(), patchI)
{
K_.boundaryField()[patchI] == interpolateXY
kappa_.boundaryField()[patchI] == interpolateXY
(
T_.boundaryField()[patchI],
TValues_,
KValues_
kappaValues_
);
}
@ -122,31 +128,31 @@ void Foam::isotropicKSolidThermo::correct()
}
Foam::tmp<Foam::volScalarField> Foam::isotropicKSolidThermo::K() const
Foam::tmp<Foam::volScalarField> Foam::isotropicKSolidThermo::kappa() const
{
return K_;
return kappa_;
}
Foam::tmp<Foam::scalarField> Foam::isotropicKSolidThermo::K
Foam::tmp<Foam::scalarField> Foam::isotropicKSolidThermo::kappa
(
const label patchI
) const
{
return K_.boundaryField()[patchI];
return kappa_.boundaryField()[patchI];
}
bool Foam::isotropicKSolidThermo::read()
{
KValues_ = Field<scalar>(subDict(typeName + "Coeffs").lookup("KValues"));
kappaValues_ = Field<scalar>(subDict(typeName + "Coeffs").lookup("kappaValues"));
return true;
}
bool Foam::isotropicKSolidThermo::writeData(Ostream& os) const
{
os.writeKeyword("KValues") << KValues_ << token::END_STATEMENT << nl;
os.writeKeyword("kappaValues") << kappaValues_ << token::END_STATEMENT << nl;
bool ok = interpolatedSolidThermo::writeData(os);
return ok && os.good();

16
src/thermophysicalModels/basicSolidThermo/isotropicKSolidThermo/isotropicKSolidThermo.H
View File

@ -53,10 +53,10 @@ class isotropicKSolidThermo
// Private data
//- Thermal conductivity [W/(m.K)]
volScalarField K_;
volScalarField kappa_;
//- Thermal conductivity vector
Field<scalar> KValues_;
Field<scalar> kappaValues_;
public:
@ -87,26 +87,26 @@ public:
// Access functions
//- Constant access to K
virtual tmp<volScalarField> K() const;
virtual tmp<volScalarField> kappa() const;
//- Thermal conductivity [W/m/K]
virtual tmp<volSymmTensorField> directionalK() const
virtual tmp<volSymmTensorField> directionalKappa() const
{
notImplemented("directionalKSolidThermo::directionalK() const");
notImplemented("directionalKSolidThermo::directionalKappa() const");
return tmp<volSymmTensorField>(NULL);
}
// Per patch calculation
//- Thermal conductivity [W//m/K]
virtual tmp<scalarField> K(const label patchI) const;
virtual tmp<scalarField> kappa(const label patchI) const;
//- Thermal conductivity [W//m/K]
virtual tmp<symmTensorField> directionalK(const label) const
virtual tmp<symmTensorField> directionalKappa(const label) const
{
notImplemented
(
"directionalKSolidThermo::directionalK(const label) const"
"directionalKSolidThermo::directionalKappa(const label) const"
);
return tmp<symmTensorField>(NULL);
}

4
src/thermophysicalModels/basicSolidThermo/solidMixtureThermo/mixtures/basicSolidMixture/basicSolidMixture.H
View File

@ -114,13 +114,13 @@ public:
virtual scalar rho(scalar T, label celli) const = 0;
//- Absorption coefficient
virtual scalar kappa(scalar T, label celli) const = 0;
virtual scalar kappaRad(scalar T, label celli) const = 0;
//- Scatter coefficient
virtual scalar sigmaS(scalar T, label celli) const = 0;
//- Thermal conductivity
virtual scalar K(scalar T, label celli) const = 0;
virtual scalar kappa(scalar T, label celli) const = 0;
//- Emissivity coefficient
virtual scalar emissivity(scalar T, label celli) const = 0;

8
src/thermophysicalModels/basicSolidThermo/solidMixtureThermo/mixtures/multiComponentSolidMixture/multiComponentSolidMixture.C
View File

@ -159,7 +159,7 @@ Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::h
template<class ThermoSolidType>
Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::kappa
Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::kappaRad
(
scalar T, label celli
) const
@ -167,7 +167,7 @@ Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::kappa
scalar tmp = 0.0;
forAll(solidData_, i)
{
tmp += solidData_[i].kappa(T)*X(i, celli, T);
tmp += solidData_[i].kappaRad(T)*X(i, celli, T);
}
return tmp;
}
@ -189,7 +189,7 @@ Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::sigmaS
template<class ThermoSolidType>
Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::K
Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::kappa
(
scalar T, label celli
) const
@ -197,7 +197,7 @@ Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::K
scalar tmp = 0.0;
forAll(solidData_, i)
{
tmp += solidData_[i].K(T)*X(i, celli, T);
tmp += solidData_[i].kappa(T)*X(i, celli, T);
}
return tmp;
}

4
src/thermophysicalModels/basicSolidThermo/solidMixtureThermo/mixtures/multiComponentSolidMixture/multiComponentSolidMixture.H
View File

@ -99,13 +99,13 @@ public:
virtual scalar rho(scalar T, label celli) const;
//- Absorption coefficient
virtual scalar kappa(scalar T, label celli) const;
virtual scalar kappaRad(scalar T, label celli) const;
//- Scatter coefficient
virtual scalar sigmaS(scalar T, label celli) const;
//- Thermal conductivity
virtual scalar K(scalar T, label celli) const;
virtual scalar kappa(scalar T, label celli) const;
//- Emissivity coefficient
virtual scalar emissivity(scalar T, label celli) const;

42
src/thermophysicalModels/basicSolidThermo/solidMixtureThermo/solidMixtureThermo/solidMixtureThermo.C
View File

@ -34,8 +34,8 @@ void Foam::solidMixtureThermo<MixtureType>::calculate()
{
scalarField& rhoCells = rho_.internalField();
scalarField& KCells = K_.internalField();
scalarField& kappaCells = kappa_.internalField();
scalarField& kappaRadCells = kappaRad_.internalField();
scalarField& sigmaSCells = sigmaS_.internalField();
scalarField& emissivityCells = emissivity_.internalField();
@ -44,9 +44,9 @@ void Foam::solidMixtureThermo<MixtureType>::calculate()
forAll(iT, celli)
{
rhoCells[celli] = MixtureType::rho(iT[celli], celli);
kappaCells[celli] = MixtureType::kappa(iT[celli], celli);
kappaRadCells[celli] = MixtureType::kappaRad(iT[celli], celli);
sigmaSCells[celli] = MixtureType::sigmaS(iT[celli], celli);
KCells[celli] = MixtureType::K(iT[celli], celli);
kappaCells[celli] = MixtureType::kappa(iT[celli], celli);
emissivityCells[celli] = MixtureType::emissivity(iT[celli], celli);
}
@ -55,8 +55,8 @@ void Foam::solidMixtureThermo<MixtureType>::calculate()
forAll(bT, patchI)
{
rho_.boundaryField()[patchI] == this->rho(patchI)();
K_.boundaryField()[patchI] == this->K(patchI)();
kappa_.boundaryField()[patchI] == this->kappa(patchI)();
kappaRad_.boundaryField()[patchI] == this->kappaRad(patchI)();
sigmaS_.boundaryField()[patchI] == this->sigmaS(patchI)();
emissivity_.boundaryField()[patchI] == this->emissivity(patchI)();
}
@ -109,7 +109,7 @@ Foam::tmp<Foam::scalarField> Foam::solidMixtureThermo<MixtureType>::sigmaS
template<class MixtureType>
Foam::tmp<Foam::scalarField> Foam::solidMixtureThermo<MixtureType>::kappa
Foam::tmp<Foam::scalarField> Foam::solidMixtureThermo<MixtureType>::kappaRad
(
const label patchI
) const
@ -118,16 +118,16 @@ Foam::tmp<Foam::scalarField> Foam::solidMixtureThermo<MixtureType>::kappa
const polyPatch& pp = mesh_.boundaryMesh()[patchI];
const labelUList& cells = pp.faceCells();
tmp<scalarField> tKappa(new scalarField(patchT.size()));
scalarField& kappa = tKappa();
tmp<scalarField> tKappaRad(new scalarField(patchT.size()));
scalarField& kappaRad = tKappaRad();
forAll(patchT, celli)
{
kappa[celli] =
MixtureType::kappa(patchT[celli], cells[celli]);
kappaRad[celli] =
MixtureType::kappaRad(patchT[celli], cells[celli]);
}
return tKappa;
return tKappaRad;
}
@ -163,11 +163,11 @@ Foam::solidMixtureThermo<MixtureType>::solidMixtureThermo
:
basicSolidThermo(mesh),
MixtureType(*this, mesh),
K_
kappa_
(
IOobject
(
"K",
"kappa",
mesh.time().timeName(),
mesh,
IOobject::NO_READ,
@ -190,11 +190,11 @@ Foam::solidMixtureThermo<MixtureType>::solidMixtureThermo
:
basicSolidThermo(mesh, dict),
MixtureType(*this, mesh),
K_
kappa_
(
IOobject
(
"K",
"kappa",
mesh.time().timeName(),
mesh,
IOobject::NO_READ,
@ -225,9 +225,9 @@ void Foam::solidMixtureThermo<MixtureType>::correct()
template<class MixtureType>
Foam::tmp<Foam::volScalarField> Foam::solidMixtureThermo<MixtureType>::K() const
Foam::tmp<Foam::volScalarField> Foam::solidMixtureThermo<MixtureType>::kappa() const
{
return K_;
return kappa_;
}
@ -352,7 +352,7 @@ Foam::solidMixtureThermo<MixtureType>::Hf() const
template<class MixtureType>
Foam::tmp<Foam::scalarField> Foam::solidMixtureThermo<MixtureType>::K
Foam::tmp<Foam::scalarField> Foam::solidMixtureThermo<MixtureType>::kappa
(
const label patchI
) const
@ -361,15 +361,15 @@ Foam::tmp<Foam::scalarField> Foam::solidMixtureThermo<MixtureType>::K
const polyPatch& pp = mesh_.boundaryMesh()[patchI];
const labelUList& cells = pp.faceCells();
tmp<scalarField> tK(new scalarField(patchT.size()));
scalarField& K = tK();
tmp<scalarField> tkappa(new scalarField(patchT.size()));
scalarField& kappa = tkappa();
forAll(patchT, celli)
{
K[celli] = MixtureType::K(patchT[celli], cells[celli]);
kappa[celli] = MixtureType::kappa(patchT[celli], cells[celli]);
}
return tK;
return tkappa;
}

16
src/thermophysicalModels/basicSolidThermo/solidMixtureThermo/solidMixtureThermo/solidMixtureThermo.H
View File

@ -57,7 +57,7 @@ protected:
// Protected data
//- Thermal conductivity [W/m/K]
volScalarField K_;
volScalarField kappa_;
private:
@ -74,7 +74,7 @@ private:
tmp<scalarField> sigmaS(const label patchI) const;
tmp<scalarField> kappa(const label patchI) const;
tmp<scalarField> kappaRad(const label patchI) const;
tmp<scalarField> emissivity(const label patchI) const;
@ -118,12 +118,12 @@ public:
// Derived thermal properties
//- Thermal conductivity [W/m/K]
virtual tmp<volScalarField> K() const;
virtual tmp<volScalarField> kappa() const;
//- Thermal conductivity [W/m/K]
virtual tmp<volSymmTensorField> directionalK() const
virtual tmp<volSymmTensorField> directionalKappa() const
{
notImplemented("solidMixtureThermo::directionalK() const");
notImplemented("solidMixtureThermo::directionalKappa() const");
return tmp<volSymmTensorField>(NULL);
}
@ -140,14 +140,14 @@ public:
// Per patch calculation
//- Thermal conductivity [W/(m.K)]
virtual tmp<scalarField> K(const label patchI) const;
virtual tmp<scalarField> kappa(const label patchI) const;
//- Thermal conductivity [W/(m.K)]
virtual tmp<symmTensorField> directionalK(const label) const
virtual tmp<symmTensorField> directionalKappa(const label) const
{
notImplemented
(
"solidMixtureThermo::directionalK(const label) const"
"solidMixtureThermo::directionalKappa(const label) const"
);
return tmp<symmTensorField>(NULL);
}

4
src/thermophysicalModels/solid/radiation/const/constSolidRad.C
View File

@ -37,7 +37,7 @@ template<class thermo>
constSolidRad<thermo>::constSolidRad(const dictionary& dict)
:
thermo(dict),
kappa_(readScalar(dict.subDict("radiation").lookup("kappa"))),
kappaRad_(readScalar(dict.subDict("radiation").lookup("kappaRad"))),
sigmaS_(readScalar(dict.subDict("radiation").lookup("sigmaS"))),
emissivity_(readScalar(dict.subDict("radiation").lookup("emissivity")))
{}
@ -49,7 +49,7 @@ template<class thermo>
Ostream& operator<<(Ostream& os, const constSolidRad<thermo>& pg)
{
os << static_cast<const thermo&>(pg);
os << tab << pg.kappa_ << tab << pg.sigmaS_ << tab << pg.emissivity_;
os << tab << pg.kappaRad_ << tab << pg.sigmaS_ << tab << pg.emissivity_;
os.check("Ostream& operator<<(Ostream& os, const constSolidRad& st)");
return os;

6
src/thermophysicalModels/solid/radiation/const/constSolidRad.H
View File

@ -62,7 +62,7 @@ class constSolidRad
// Private data
//- Absorption coefficient
scalar kappa_;
scalar kappaRad_;
//- Scattering coefficient
scalar sigmaS_;
@ -75,7 +75,7 @@ class constSolidRad
inline constSolidRad
(
const thermo& t,
const scalar kappa,
const scalar kappaRad,
const scalar sigmaS,
const scalar emissivity
);
@ -101,7 +101,7 @@ public:
// Member functions
//- Return absorption coefficient [1/m]
inline scalar kappa(scalar T) const;
inline scalar kappaRad(scalar T) const;
//- Return scattering [1/m]
inline scalar sigmaS(scalar T) const;

10
src/thermophysicalModels/solid/radiation/const/constSolidRadI.H
View File

@ -36,13 +36,13 @@ template<class thermo>
inline constSolidRad<thermo>::constSolidRad
(
const thermo& t,
const scalar kappa,
const scalar kappaRad,
const scalar sigmaS,
const scalar emissivity
)
:
thermo(t),
kappa_(kappa),
kappaRad_(kappaRad),
sigmaS_(sigmaS),
emissivity_(emissivity)
{}
@ -59,7 +59,7 @@ inline constSolidRad<thermo>::constSolidRad
)
:
thermo(name, pg),
kappa_(pg.kappa_),
kappaRad_(pg.kappaRad_),
sigmaS_(pg.sigmaS_),
emissivity_(pg.emissivity_)
{}
@ -82,9 +82,9 @@ inline autoPtr<constSolidRad<thermo> > constSolidRad<thermo>::New(Istream& is)
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
template<class thermo>
inline scalar constSolidRad<thermo>::kappa(scalar) const
inline scalar constSolidRad<thermo>::kappaRad(scalar) const
{
return kappa_;
return kappaRad_;
}

4
src/thermophysicalModels/solid/transport/const/constSolidTransport.C
View File

@ -32,7 +32,7 @@ template<class thermo>
Foam::constSolidTransport<thermo>::constSolidTransport(const dictionary& dict)
:
thermo(dict),
K_(readScalar(dict.subDict("transport").lookup("K")))
kappa_(readScalar(dict.subDict("transport").lookup("kappa")))
{}
@ -46,7 +46,7 @@ Foam::Ostream& Foam::operator<<
)
{
operator<<(os, static_cast<const thermo&>(ct));
os << tab << ct.K_;
os << tab << ct.kappa_;
os.check("Ostream& operator<<(Ostream& os,const constSolidTransport& ct)");

6
src/thermophysicalModels/solid/transport/const/constSolidTransport.H
View File

@ -65,7 +65,7 @@ class constSolidTransport
// Private data
//- Constant thermal conductivity.
scalar K_;
scalar kappa_;
// Private Member Functions
@ -74,7 +74,7 @@ class constSolidTransport
inline constSolidTransport
(
const thermo& t,
const scalar k
const scalar kappa
);
@ -93,7 +93,7 @@ public:
// Member functions
//- Thermal conductivity [W/mK]
inline scalar K(const scalar T) const;
inline scalar kappa(const scalar T) const;
// Ostream Operator

10
src/thermophysicalModels/solid/transport/const/constSolidTransportI.H
View File

@ -29,11 +29,11 @@ template<class thermo>
inline Foam::constSolidTransport<thermo>::constSolidTransport
(
const thermo& t,
const scalar k
const scalar kappa
)
:
thermo(t),
K_(k)
kappa_(kappa)
{}
@ -45,16 +45,16 @@ inline Foam::constSolidTransport<thermo>::constSolidTransport
)
:
thermo(name, ct),
K_(ct.K_)
kappa_(ct.kappa_)
{}
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
template<class thermo>
inline Foam::scalar Foam::constSolidTransport<thermo>::K(const scalar T) const
inline Foam::scalar Foam::constSolidTransport<thermo>::kappa(const scalar T) const
{
return K_;
return kappa_;
}

6
src/thermophysicalModels/solid/transport/exponential/exponentialSolidTransport.C
View File

@ -35,12 +35,12 @@ Foam::exponentialSolidTransport<thermo>::exponentialSolidTransport
)
:
thermo(dict),
K0_(0.0),
kappa0_(0.0),
n0_(0.0),
Tref_(0.0)
{
const dictionary& subDict = dict.subDict("transport");
K0_ = readScalar(subDict.lookup("K0"));
kappa0_ = readScalar(subDict.lookup("kappa0"));
n0_ = readScalar(subDict.lookup("n0"));
Tref_ = readScalar(subDict.lookup("Tref"));
}
@ -55,7 +55,7 @@ Foam::Ostream& Foam::operator<<
)
{
operator<<(os, static_cast<const thermo&>(et));
os << tab << et.K0_ << tab << et.n0_ << tab << et.Tref_;
os << tab << et.kappa0_ << tab << et.n0_ << tab << et.Tref_;
os.check
(

6
src/thermophysicalModels/solid/transport/exponential/exponentialSolidTransport.H
View File

@ -64,7 +64,7 @@ class exponentialSolidTransport
// Private data
//- Constant thermal coefficient.
scalar K0_;
scalar kappa0_;
//- Exponent coefficient
scalar n0_;
@ -79,7 +79,7 @@ class exponentialSolidTransport
inline exponentialSolidTransport
(
const thermo& t,
const scalar K0,
const scalar kappa0,
const scalar n0,
const scalar Tref
);
@ -104,7 +104,7 @@ public:
// Member functions
//- Thermal conductivity [W/mK]
inline scalar K(const scalar T) const;
inline scalar kappa(const scalar T) const;
// Ostream Operator

10
src/thermophysicalModels/solid/transport/exponential/exponentialSolidTransportI.H
View File

@ -29,13 +29,13 @@ template<class thermo>
inline Foam::exponentialSolidTransport<thermo>::exponentialSolidTransport
(
const thermo& t,
const scalar K0,
const scalar kappa0,
const scalar n0,
const scalar Tref
)
:
thermo(t),
K0_(K0),
kappa0_(kappa0),
n0_(n0),
Tref_(Tref)
{}
@ -49,7 +49,7 @@ inline Foam::exponentialSolidTransport<thermo>::exponentialSolidTransport
)
:
thermo(name, ct),
K0_(ct.K0_),
kappa0_(ct.kappa0_),
n0_(ct.n0_),
Tref_(ct.Tref_)
{}
@ -58,12 +58,12 @@ inline Foam::exponentialSolidTransport<thermo>::exponentialSolidTransport
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
template<class thermo>
inline Foam::scalar Foam::exponentialSolidTransport<thermo>::K
inline Foam::scalar Foam::exponentialSolidTransport<thermo>::kappa
(
const scalar T
) const
{
return (K0_*pow(T/Tref_, n0_));
return (kappa0_*pow(T/Tref_, n0_));
}

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

@ -86,7 +86,7 @@ externalWallHeatFluxTemperatureFvPatchScalarField
)
:
mixedFvPatchScalarField(ptf, p, iF, mapper),
temperatureCoupledBase(patch(), ptf.KMethod(), ptf.KName()),
temperatureCoupledBase(patch(), ptf.KMethod(), ptf.kappaName()),
oldMode_(unknown),
q_(ptf.q_, mapper),
h_(ptf.h_, mapper),
@ -181,7 +181,7 @@ externalWallHeatFluxTemperatureFvPatchScalarField
)
:
mixedFvPatchScalarField(tppsf, iF),
temperatureCoupledBase(patch(), tppsf.KMethod(), tppsf.KName()),
temperatureCoupledBase(patch(), tppsf.KMethod(), tppsf.kappaName()),
oldMode_(tppsf.oldMode_),
q_(tppsf.q_),
h_(tppsf.h_),
@ -199,7 +199,7 @@ void Foam::externalWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
}
scalarField q(size(), 0.0);
scalarField KDelta(K(*this)*patch().deltaCoeffs());
scalarField KDelta(kappa(*this)*patch().deltaCoeffs());
if (oldMode_ == fixedHeatFlux)
{
@ -223,7 +223,7 @@ void Foam::externalWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
{
if (q[i] > 0) //in
{
this->refGrad()[i] = q[i]/K(*this)()[i];
this->refGrad()[i] = q[i]/kappa(*this)()[i];
this->refValue()[i] = 0.0;
this->valueFraction()[i] = 0.0;
}
@ -239,7 +239,7 @@ void Foam::externalWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
if (debug)
{
scalar Q = gSum(K(*this)*patch().magSf()*snGrad());
scalar Q = gSum(kappa(*this)*patch().magSf()*snGrad());
Info<< patch().boundaryMesh().mesh().name() << ':'
<< patch().name() << ':'

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

@ -34,12 +34,12 @@ Description
myWallPatch
{
type externalWallHeatFluxTemperature;
K solidThermo; // solidThermo or lookup
kappa solidThermo; // solidThermo or lookup
q uniform 1000; // Heat flux / [W/m2]
Ta uniform 300.0; // Tambient temperature /[K]
h uniform 10.0; // Heat transfer coeff /[W/Km2]
value uniform 300.0; // Initial temperature / [K]
KName none;
kappaName none;
}
Note:

42
src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/temperatureCoupledBase/temperatureCoupledBase.C
View File

@ -58,12 +58,12 @@ Foam::temperatureCoupledBase::temperatureCoupledBase
(
const fvPatch& patch,
const word& calculationType,
const word& KName
const word& kappaName
)
:
patch_(patch),
method_(KMethodTypeNames_[calculationType]),
KName_(KName)
kappaName_(kappaName)
{}
@ -74,14 +74,14 @@ Foam::temperatureCoupledBase::temperatureCoupledBase
)
:
patch_(patch),
method_(KMethodTypeNames_.read(dict.lookup("K"))),
KName_(dict.lookup("KName"))
method_(KMethodTypeNames_.read(dict.lookup("kappa"))),
kappaName_(dict.lookup("kappaName"))
{}
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
Foam::tmp<Foam::scalarField> Foam::temperatureCoupledBase::K
Foam::tmp<Foam::scalarField> Foam::temperatureCoupledBase::kappa
(
const scalarField& Tp
) const
@ -111,7 +111,7 @@ Foam::tmp<Foam::scalarField> Foam::temperatureCoupledBase::K
(
"solidThermophysicalProperties"
);
return thermo.K(patch_.index());
return thermo.kappa(patch_.index());
}
break;
@ -124,23 +124,23 @@ Foam::tmp<Foam::scalarField> Foam::temperatureCoupledBase::K
(
"solidThermophysicalProperties"
);
return n & thermo.directionalK(patch_.index()) & n;
return n & thermo.directionalKappa(patch_.index()) & n;
}
break;
case LOOKUP:
{
if (mesh.objectRegistry::foundObject<volScalarField>(KName_))
if (mesh.objectRegistry::foundObject<volScalarField>(kappaName_))
{
return patch_.lookupPatchField<volScalarField, scalar>(KName_);
return patch_.lookupPatchField<volScalarField, scalar>(kappaName_);
}
else if
(
mesh.objectRegistry::foundObject<volSymmTensorField>(KName_)
mesh.objectRegistry::foundObject<volSymmTensorField>(kappaName_)
)
{
const symmTensorField& KWall =
patch_.lookupPatchField<volSymmTensorField, scalar>(KName_);
patch_.lookupPatchField<volSymmTensorField, scalar>(kappaName_);
const vectorField n(patch_.nf());
@ -149,12 +149,12 @@ Foam::tmp<Foam::scalarField> Foam::temperatureCoupledBase::K
else
{
FatalErrorIn("temperatureCoupledBase::K() const")
<< "Did not find field " << KName_
<< "Did not find field " << kappaName_
<< " on mesh " << mesh.name() << " patch " << patch_.name()
<< endl
<< "Please set 'K' to one of " << KMethodTypeNames_.toc()
<< " and 'KName' to the name of the volScalar"
<< " or volSymmTensor field (if K=lookup)"
<< "Please set 'kappa' to one of " << KMethodTypeNames_.toc()
<< " and 'kappaName' to the name of the volScalar"
<< " or volSymmTensor field (if kappa=lookup)"
<< exit(FatalError);
return scalarField(0);
@ -163,11 +163,11 @@ Foam::tmp<Foam::scalarField> Foam::temperatureCoupledBase::K
default:
{
FatalErrorIn("temperatureCoupledBase::K() const")
FatalErrorIn("temperatureCoupledBase::kappa() const")
<< "Unimplemented method " << method_ << endl
<< "Please set 'K' to one of " << KMethodTypeNames_.toc()
<< " and 'KName' to the name of the volScalar"
<< " or volSymmTensor field (if K=lookup)"
<< "Please set 'kappa' to one of " << KMethodTypeNames_.toc()
<< " and 'kappaName' to the name of the volScalar"
<< " or volSymmTensor field (if kappa=lookup)"
<< exit(FatalError);
}
break;
@ -178,9 +178,9 @@ Foam::tmp<Foam::scalarField> Foam::temperatureCoupledBase::K
void Foam::temperatureCoupledBase::write(Ostream& os) const
{
os.writeKeyword("K") << KMethodTypeNames_[method_]
os.writeKeyword("kappa") << KMethodTypeNames_[method_]
<< token::END_STATEMENT << nl;
os.writeKeyword("KName") << KName_ << token::END_STATEMENT << nl;
os.writeKeyword("kappaName") << kappaName_ << token::END_STATEMENT << nl;
}

19
src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/temperatureCoupledBase/temperatureCoupledBase.H
View File

@ -27,13 +27,14 @@ Class
Description
Common functions for use in temperature coupled boundaries. For now only
K() : heat conduction at patch. Gets supplied how to lookup/calculate K:
kappa() : heat conduction at patch. Gets supplied how to lookup/calculate
kappa:
- 'lookup' : lookup volScalarField (or volSymmTensorField) with name
- 'basicThermo' : use basicThermo and default compressible::turbulenceModel
to calculate K
- 'solidThermo' : use basicSolidThermo K()
- 'directionalSolidThermo' directionalK()
to calculate kappa
- 'solidThermo' : use basicSolidThermo kappa()
- 'directionalSolidThermo' directionalKappa()
SourceFiles
temperatureCoupledBase.C
@ -59,7 +60,7 @@ namespace Foam
class temperatureCoupledBase
{
public:
//- Type of supplied K
//- Type of supplied Kappa
enum KMethodType
{
BASICTHERMO,
@ -81,7 +82,7 @@ private:
const KMethodType method_;
//- Name of thermal conductivity field (if looked up from database)
const word KName_;
const word kappaName_;
public:
@ -113,13 +114,13 @@ public:
}
//- Name of thermal conductivity field
const word& KName() const
const word& kappaName() const
{
return KName_;
return kappaName_;
}
//- Given patch temperature calculate corresponding K field
tmp<scalarField> K(const scalarField& Tp) const;
tmp<scalarField> kappa(const scalarField& Tp) const;
//- Write
void write(Ostream&) const;

2
src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/temperatureThermoBaffle1D/temperatureThermoBaffle1DFvPatchScalarField.C
View File

@ -292,7 +292,7 @@ void temperatureThermoBaffle1DFvPatchScalarField<solidType>::updateCoeffs()
// Create fields for solid properties
forAll(KDeltaw, i)
{
KDeltaw[i] = solid_().K((Tp[i] + nbrTp[i])/2.0)/thickness_[i];
KDeltaw[i] = solid_().kappa((Tp[i] + nbrTp[i])/2.0)/thickness_[i];
}
const scalarField q

4
src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/temperatureThermoBaffle1D/temperatureThermoBaffle1DFvPatchScalarField.H
View File

@ -114,9 +114,9 @@ class temperatureThermoBaffle1DFvPatchScalarField
}
scalar K(const scalar T) const
scalar kappa(const scalar T) const
{
return solidPtr_().K(T);
return solidPtr_().kappa(T);
}
};

10
src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.C
View File

@ -93,7 +93,7 @@ turbulentHeatFluxTemperatureFvPatchScalarField
)
:
fixedGradientFvPatchScalarField(ptf, p, iF, mapper),
temperatureCoupledBase(patch(), ptf.KMethod(), ptf.KName()),
temperatureCoupledBase(patch(), ptf.KMethod(), ptf.kappaName()),
heatSource_(ptf.heatSource_),
q_(ptf.q_, mapper)
{}
@ -124,7 +124,7 @@ turbulentHeatFluxTemperatureFvPatchScalarField
)
:
fixedGradientFvPatchScalarField(thftpsf),
temperatureCoupledBase(patch(), thftpsf.KMethod(), thftpsf.KName()),
temperatureCoupledBase(patch(), thftpsf.KMethod(), thftpsf.kappaName()),
heatSource_(thftpsf.heatSource_),
q_(thftpsf.q_)
{}
@ -138,7 +138,7 @@ turbulentHeatFluxTemperatureFvPatchScalarField
)
:
fixedGradientFvPatchScalarField(thftpsf, iF),
temperatureCoupledBase(patch(), thftpsf.KMethod(), thftpsf.KName()),
temperatureCoupledBase(patch(), thftpsf.KMethod(), thftpsf.kappaName()),
heatSource_(thftpsf.heatSource_),
q_(thftpsf.q_)
{}
@ -188,12 +188,12 @@ void turbulentHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
case hsPower:
{
const scalar Ap = gSum(patch().magSf());
gradient() = q_/(Ap*K(Tp));
gradient() = q_/(Ap*kappa(Tp));
break;
}
case hsFlux:
{
gradient() = q_/K(Tp);
gradient() = q_/kappa(Tp);
break;
}
default:

2
src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.H
View File

@ -36,7 +36,7 @@ Description
type compressible::turbulentHeatFluxTemperature;
heatSource flux; // power [W]; flux [W/m2]
q uniform 10; // heat power or flux
K basicThermo; // calculate K by alphaEff*thermo.Cp
kappa basicThermo; // calculate kappa=alphaEff*thermo.Cp
value uniform 300; // initial temperature value
}

12
src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/turbulentTemperatureCoupledBaffleMixed/turbulentTemperatureCoupledBaffleMixedFvPatchScalarField.C
View File

@ -68,7 +68,7 @@ turbulentTemperatureCoupledBaffleMixedFvPatchScalarField
)
:
mixedFvPatchScalarField(ptf, p, iF, mapper),
temperatureCoupledBase(patch(), ptf.KMethod(), ptf.KName()),
temperatureCoupledBase(patch(), ptf.KMethod(), ptf.kappaName()),
neighbourFieldName_(ptf.neighbourFieldName_)
{}
@ -131,7 +131,7 @@ turbulentTemperatureCoupledBaffleMixedFvPatchScalarField
)
:
mixedFvPatchScalarField(wtcsf, iF),
temperatureCoupledBase(patch(), wtcsf.KMethod(), wtcsf.KName()),
temperatureCoupledBase(patch(), wtcsf.KMethod(), wtcsf.kappaName()),
neighbourFieldName_(wtcsf.neighbourFieldName_)
{}
@ -186,11 +186,11 @@ void turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs()
scalarField nbrIntFld(nbrField.patchInternalField());
distMap.distribute(nbrIntFld);
// Swap to obtain full local values of neighbour K*delta
scalarField nbrKDelta(nbrField.K(nbrField)*nbrPatch.deltaCoeffs());
// Swap to obtain full local values of neighbour kappa*delta
scalarField nbrKDelta(nbrField.kappa(nbrField)*nbrPatch.deltaCoeffs());
distMap.distribute(nbrKDelta);
tmp<scalarField> myKDelta = K(*this)*patch().deltaCoeffs();
tmp<scalarField> myKDelta = kappa(*this)*patch().deltaCoeffs();
// Both sides agree on
@ -219,7 +219,7 @@ void turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs()
if (debug)
{
scalar Q = gSum(K(*this)*patch().magSf()*snGrad());
scalar Q = gSum(kappa(*this)*patch().magSf()*snGrad());
Info<< patch().boundaryMesh().mesh().name() << ':'
<< patch().name() << ':'

14
src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/turbulentTemperatureCoupledBaffleMixed/turbulentTemperatureCoupledBaffleMixedFvPatchScalarField.H
View File

@ -41,18 +41,20 @@ Description
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K lookup;
KName K;
kappa lookup;
kappaName kappa;
value uniform 300;
}
Needs to be on underlying mapped(Wall)FvPatch.
Note: K : heat conduction at patch. Gets supplied how to lookup/calculate K:
Note: kappa : heat conduction at patch. Gets supplied how to lookup/calculate
kappa:
- 'lookup' : lookup volScalarField (or volSymmTensorField) with name
- 'basicThermo' : use basicThermo and compressible::RASmodel to calculate K
- 'solidThermo' : use basicSolidThermo K()
- 'directionalSolidThermo' directionalK()
- 'basicThermo' : use basicThermo and compressible::RASmodel to calculate
kappa
- 'solidThermo' : use basicSolidThermo kappa()
- 'directionalSolidThermo' directionalKappa()
Note: runs in parallel with arbitrary decomposition. Uses mapped
functionality to calculate exchange.

8
src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/turbulentTemperatureRadCoupledMixed/turbulentTemperatureRadCoupledMixedFvPatchScalarField.C
View File

@ -67,7 +67,7 @@ turbulentTemperatureRadCoupledMixedFvPatchScalarField
)
:
mixedFvPatchScalarField(psf, p, iF, mapper),
temperatureCoupledBase(patch(), psf.KMethod(), psf.KName()),
temperatureCoupledBase(patch(), psf.KMethod(), psf.kappaName()),
TnbrName_(psf.TnbrName_),
QrNbrName_(psf.QrNbrName_),
QrName_(psf.QrName_)
@ -134,7 +134,7 @@ turbulentTemperatureRadCoupledMixedFvPatchScalarField
)
:
mixedFvPatchScalarField(psf, iF),
temperatureCoupledBase(patch(), psf.KMethod(), psf.KName()),
temperatureCoupledBase(patch(), psf.KMethod(), psf.kappaName()),
TnbrName_(psf.TnbrName_),
QrNbrName_(psf.QrNbrName_),
QrName_(psf.QrName_)
@ -182,11 +182,11 @@ void turbulentTemperatureRadCoupledMixedFvPatchScalarField::updateCoeffs()
// Swap to obtain full local values of neighbour K*delta
scalarField KDeltaNbr(nbrField.K(nbrField)*nbrPatch.deltaCoeffs());
scalarField KDeltaNbr(nbrField.kappa(nbrField)*nbrPatch.deltaCoeffs());
distMap.distribute(KDeltaNbr);
scalarField KDelta(K(*this)*patch().deltaCoeffs());
scalarField KDelta(kappa(*this)*patch().deltaCoeffs());
scalarField Qr(Tp.size(), 0.0);
if (QrName_ != "none")

13
src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/turbulentTemperatureRadCoupledMixed/turbulentTemperatureRadCoupledMixedFvPatchScalarField.H
View File

@ -35,8 +35,8 @@ Description
{
type compressible::turbulentTemperatureRadCoupledMixed;
TNbr T; // name of T field on neighbour region
K lookup;
KName K;
kappa lookup;
KappaName kappa;
QrNbr Qr; // or none. Name of Qr field on neighbour region
Qr Qr; // or none. Name of Qr field on local region
value uniform 300;
@ -44,12 +44,13 @@ Description
Needs to be on underlying mapped(Wall)FvPatch.
Note: K : heat conduction at patch. Gets supplied how to lookup/calculate
K:
Note: kappa : heat conduction at patch. Gets supplied how to
lookup/calculate
kappa:
- 'lookup' : lookup volScalarField (or volSymmTensorField) with name
- 'basicThermo' : use basicThermo and compressible::RASmodel to calculate K
- 'solidThermo' : use basicSolidThermo K()
- 'directionalSolidThermo' directionalK()
- 'solidThermo' : use basicSolidThermo kappa()
- 'directionalSolidThermo' directionalKappa()
Note: runs in parallel with arbitrary decomposition. Uses mapped
functionality to calculate exchange.

4
tutorials/combustion/fireFoam/les/oppositeBurningPanels/0/T
View File

@ -50,10 +50,10 @@ boundaryField
{
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
K basicThermo;
kappa basicThermo;
QrNbr none;
Qr Qr;
KName none;
kappaName none;
value $internalField;
}
}

4
tutorials/combustion/fireFoam/les/oppositeBurningPanels/0/panelRegion/T
View File

@ -41,8 +41,8 @@ boundaryField
{
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
QrNbr Qr;
Qr none;
value uniform 298.15;

8
tutorials/combustion/fireFoam/les/oppositeBurningPanels/constant/panelRegion/solidThermophysicalProperties
View File

@ -26,12 +26,12 @@ vCoeffs
{
transport
{
K 0.135;
kappa 0.135;
}
radiation
{
sigmaS 0.0;
kappa 0.0; //opaque
kappaRad 0.0; //opaque
emissivity 0.17;
}
thermodynamics
@ -49,12 +49,12 @@ charCoeffs
{
transport
{
K 0.4;
kappa 0.4;
}
radiation
{
sigmaS 0.0;
kappa 0.0;
kappaRad 0.0;
emissivity 0.85;
}
thermodynamics

4
tutorials/heatTransfer/buoyantSimpleFoam/circuitBoardCooling/system/baffleRegion/changeDictionaryDict
View File

@ -24,8 +24,8 @@ dictionaryReplacement
{
type compressible::temperatureThermoBaffle;
neighbourFieldName T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
value uniform 300;
}

4
tutorials/heatTransfer/buoyantSimpleFoam/circuitBoardCooling/system/changeDictionaryDict
View File

@ -94,12 +94,12 @@ dictionaryReplacement
Qs uniform 100; // heat flux [W/m2]
transport
{
K 1.0;
kappa 1.0;
}
radiation
{
sigmaS 0;
kappa 0;
kappaRad 0;
emissivity 0;
}
thermodynamics

4
tutorials/heatTransfer/buoyantSimpleFoam/circuitBoardCooling/system/changeDictionaryDict.baffle
View File

@ -94,12 +94,12 @@ dictionaryReplacement
Qs uniform 300; // heat flux [W/m2]
transport
{
K 1.0;
kappa 1.0;
}
radiation
{
sigmaS 0;
kappa 0;
kappaRad 0;
emissivity 0;
}
thermodynamics

6
tutorials/heatTransfer/buoyantSimpleFoam/circuitBoardCooling/system/changeDictionaryDict.baffleRegion
View File

@ -95,7 +95,7 @@ dictionaryReplacement
// Coupled BC.
neighbourFieldName T;
K basicThermo;
kappa basicThermo;
KName none;
@ -117,10 +117,10 @@ dictionaryReplacement
//- thermo properties
rho rho [1 -3 0 0 0 0 0] 80;
Cp Cp [0 2 -2 -1 0 0 0] 15;
K K [1 1 -3 -1 0 0 0] 0.01;
kappa kappa [1 1 -3 -1 0 0 0] 0.01;
//- radiative properties
kappa kappa [0 -1 0 0 0 0 0] 0;
kappaRad kappaRad [0 -1 0 0 0 0 0] 0;
sigmaS sigmaS [0 -1 0 0 0 0 0] 0;
emissivity emissivity [0 0 0 0 0 0 0] 1;

2
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/Allrun
View File

@ -37,7 +37,7 @@ do
done
# Run
#runParallel `getApplication` 4
runParallel `getApplication` 4
# Reconstruct
#for i in bottomAir topAir heater leftSolid rightSolid

21
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/constant/heater/solidThermophysicalProperties
View File

@ -23,10 +23,10 @@ constSolidThermoCoeffs
//- thermo properties
rho rho [1 -3 0 0 0 0 0] 8000;
Cp Cp [0 2 -2 -1 0 0 0] 450;
K K [1 1 -3 -1 0 0 0] 80;
kappa kappa [1 1 -3 -1 0 0 0] 80;
//- radiation properties
kappa kappa [0 -1 0 0 0 0 0] 0;
kappaRad kappaRad [0 -1 0 0 0 0 0] 0;
sigmaS sigmaS [0 -1 0 0 0 0 0] 0;
emissivity emissivity [0 0 0 0 0 0 0] 0;
@ -41,11 +41,11 @@ isotropicKSolidThermoCoeffs
TValues (100 1000);
rhoValues (1000 1700);
cpValues (1700 1700);
KValues (80 40);
kappaValues (80 40);
//- radiation properties
emissivityValues (1 1);
kappaValues (0 0);
kappaRadValues (0 0);
sigmaSValues (0 0);
//- chemical properties
@ -63,7 +63,7 @@ directionalKSolidThermoCoeffs
TValues (100 1000);
rhoValues (1000 1700);
cpValues (1700 1700);
KValues ((10 10 10) (40 40 40));
kappaValues ((10 10 10) (40 40 40));
coordinateSystem
{
@ -81,7 +81,7 @@ directionalKSolidThermoCoeffs
//- radiation properties
emissivityValues (1 1);
kappaValues (1 1);
kappaRadValues (1 1);
sigmaSValues (1 1);
}
@ -94,14 +94,15 @@ pmmaCoeffs
{
transport
{
K0 0.152;
kappa0 0.152;
n0 1.3;
Tref 300; //K0*pow(T/Tref_, n0))
}
radiation
{
sigmaS 0.0;
kappa 0.0;
kappaRad 0.0;
emissivity 0.1;
}
thermodynamics
@ -121,7 +122,7 @@ charCoeffs
{
transport
{
K0 0.4;
kappa0 0.4;
n0 1.3;
Tref 300; //K0*pow(T/Tref_, n0))
}
@ -129,7 +130,7 @@ charCoeffs
radiation
{
sigmaS 0.0;
kappa 0.0;
kappaRad 0.0;
emissivity 0.0;
}

4
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/bottomAir/changeDictionaryDict
View File

@ -57,8 +57,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K basicThermo;
KName none;
kappa basicThermo;
kappaName none;
value uniform 300;
}
}

4
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/heater/changeDictionaryDict
View File

@ -47,8 +47,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
value uniform 300;
}
minY

4
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/leftSolid/changeDictionaryDict
View File

@ -43,8 +43,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
value uniform 300;
}
}

4
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/rightSolid/changeDictionaryDict
View File

@ -43,8 +43,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
value uniform 300;
}
}

4
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/topAir/changeDictionaryDict
View File

@ -68,8 +68,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K basicThermo;
KName none;
kappa basicThermo;
kappaName none;
value uniform 300;
}
}

4
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionLiquidHeater/system/bottomWater/changeDictionaryDict
View File

@ -70,8 +70,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K basicThermo;
KName none;
kappa basicThermo;
kappaName none;
value uniform 300;
}
}

4
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionLiquidHeater/system/bottomWater/changeDictionaryDict.save
View File

@ -57,8 +57,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K basicThermo;
KName none;
kappa basicThermo;
kappaName none;
value uniform 300;
}
}

4
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionLiquidHeater/system/heater/changeDictionaryDict
View File

@ -48,8 +48,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
value uniform 300;
}
minY

4
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionLiquidHeater/system/leftSolid/changeDictionaryDict
View File

@ -43,8 +43,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
value uniform 300;
}
}

4
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionLiquidHeater/system/rightSolid/changeDictionaryDict
View File

@ -43,8 +43,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
value uniform 300;
}
}

4
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionLiquidHeater/system/topAir/changeDictionaryDict
View File

@ -69,8 +69,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K basicThermo;
KName none;
kappa basicThermo;
kappaName none;
value uniform 300;
}
}

20
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/constant/heater/solidThermophysicalProperties
View File

@ -25,10 +25,10 @@ constSolidThermoCoeffs
//- thermo properties
rho rho [1 -3 0 0 0 0 0] 8000;
Cp Cp [0 2 -2 -1 0 0 0] 450;
K K [1 1 -3 -1 0 0 0] 80;
kappa kappa [1 1 -3 -1 0 0 0] 80;
//- radiation properties
kappa kappa [0 -1 0 0 0 0 0] 0;
kappaRad kappaRad [0 -1 0 0 0 0 0] 0;
sigmaS sigmaS [0 -1 0 0 0 0 0] 0;
emissivity emissivity [0 0 0 0 0 0 0] 1;
@ -43,11 +43,11 @@ isotropicKSolidThermoCoeffs
TValues (100 1000);
rhoValues (1000 1700);
cpValues (1700 1700);
KValues (80 40);
kappaValues (80 40);
//- radiation properties
emissivityValues (1 1);
kappaValues (0 0);
kappaRadValues (0 0);
sigmaSValues (0 0);
//- chemical properties
@ -65,7 +65,7 @@ directionalKSolidThermoCoeffs
TValues (100 1000);
rhoValues (1700 1700);
cpValues (1700 1700);
KValues ((10 10 10) (40 40 40));
kappaValues ((10 10 10) (40 40 40));
coordinateSystem
{
@ -83,7 +83,7 @@ directionalKSolidThermoCoeffs
//- radiation properties
emissivityValues (1 1);
kappaValues (1 1);
kappaRadValues (1 1);
sigmaSValues (1 1);
}
@ -96,14 +96,14 @@ pmmaCoeffs
{
transport
{
K0 0.152;
kappa0 0.152;
n0 1.3;
Tref 300; //K0*pow(T/Tref_, n0))
}
radiation
{
sigmaS 0.0;
kappa 0.0;
kappaRad 0.0;
emissivity 0.1;
}
thermodynamics
@ -123,14 +123,14 @@ charCoeffs
{
transport
{
K0 0.4;
kappa0 0.4;
n0 1.3;
Tref 300; //K0*pow(T/Tref_, n0))
}
radiation
{
sigmaS 0.0;
kappa 0.0;
kappaRad 0.0;
emissivity 0.0;
}
thermodynamics

4
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/bottomAir/changeDictionaryDict
View File

@ -57,8 +57,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K basicThermo;
KName none;
kappa basicThermo;
kappaName none;
value uniform 300;
}
}

4
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/heater/changeDictionaryDict
View File

@ -47,8 +47,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
value uniform 300;
}
minY

4
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/leftSolid/changeDictionaryDict
View File

@ -43,8 +43,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
value uniform 300;
}
}

4
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/rightSolid/changeDictionaryDict
View File

@ -43,8 +43,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
value uniform 300;
}
}

4
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/topAir/changeDictionaryDict
View File

@ -68,8 +68,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K basicThermo;
KName none;
kappa basicThermo;
kappaName none;
value uniform 300;
}
}

20
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/constant/heater/solidThermophysicalProperties
View File

@ -25,10 +25,10 @@ constSolidThermoCoeffs
//- thermo properties
rho rho [1 -3 0 0 0 0 0] 8000;
Cp Cp [0 2 -2 -1 0 0 0] 450;
K K [1 1 -3 -1 0 0 0] 80;
kappa kappa [1 1 -3 -1 0 0 0] 80;
//- radiation properties
kappa kappa [0 -1 0 0 0 0 0] 0;
kappaRad kappaRad [0 -1 0 0 0 0 0] 0;
sigmaS sigmaS [0 -1 0 0 0 0 0] 0;
emissivity emissivity [0 0 0 0 0 0 0] 1;
@ -43,11 +43,11 @@ isotropicKSolidThermoCoeffs
TValues (100 1000);
rhoValues (1000 1700);
cpValues (1700 1700);
KValues (80 40);
kappaValues (80 40);
//- radiation properties
emissivityValues (1 1);
kappaValues (0 0);
kappaRadValues (0 0);
sigmaSValues (0 0);
//- chemical properties
@ -65,7 +65,7 @@ directionalKSolidThermoCoeffs
TValues (100 1000);
rhoValues (1700 1700);
cpValues (1700 1700);
KValues ((10 10 10) (40 40 40));
kappaValues ((10 10 10) (40 40 40));
coordinateSystem
{
@ -83,7 +83,7 @@ directionalKSolidThermoCoeffs
//- radiation properties
emissivityValues (1 1);
kappaValues (1 1);
kappaRadValues (1 1);
sigmaSValues (1 1);
}
@ -96,7 +96,7 @@ pmmaCoeffs
{
transport
{
K0 0.152;
kappa0 0.152;
n0 1.3;
Tref 300; //K0*pow(T/Tref_, n0))
}
@ -104,7 +104,7 @@ pmmaCoeffs
radiation
{
sigmaS 0.0;
kappa 0.0;
kappaRad 0.0;
emissivity 0.1;
}
@ -126,7 +126,7 @@ charCoeffs
{
transport
{
K0 0.4;
kappa0 0.4;
n0 1.3;
Tref 300; //K0*pow(T/Tref_, n0))
}
@ -134,7 +134,7 @@ charCoeffs
radiation
{
sigmaS 0.0;
kappa 0.0;
kappaRad 0.0;
emissivity 0.0;
}

20
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/constant/heater/solidThermophysicalProperties
View File

@ -25,10 +25,10 @@ constSolidThermoCoeffs
//- thermo properties
rho rho [1 -3 0 0 0 0 0] 8000;
Cp Cp [0 2 -2 -1 0 0 0] 450;
K K [1 1 -3 -1 0 0 0] 80;
kappa kappa [1 1 -3 -1 0 0 0] 80;
//- radiation properties
kappa kappa [0 -1 0 0 0 0 0] 0;
kappaRad kappaRad [0 -1 0 0 0 0 0] 0;
sigmaS sigmaS [0 -1 0 0 0 0 0] 0;
emissivity emissivity [0 0 0 0 0 0 0] 0.5;
@ -43,11 +43,11 @@ isotropicKSolidThermoCoeffs
TValues (100 1000);
rhoValues (1000 1700);
cpValues (1700 1700);
KValues (80 40);
kappaValues (80 40);
//- radiation properties
emissivityValues (1 1);
kappaValues (0 0);
kappaRadValues (0 0);
sigmaSValues (0 0);
//- chemical properties
@ -65,7 +65,7 @@ directionalKSolidThermoCoeffs
TValues (100 1000);
rhoValues (1700 1700);
cpValues (1700 1700);
KValues ((10 10 10) (40 40 40));
kappaValues ((10 10 10) (40 40 40));
coordinateSystem
{
@ -83,7 +83,7 @@ directionalKSolidThermoCoeffs
//- radiation properties
emissivityValues (1 1);
kappaValues (1 1);
kappaRadValues (1 1);
sigmaSValues (1 1);
}
@ -96,7 +96,7 @@ pmmaCoeffs
{
transport
{
K0 0.152;
kappa0 0.152;
n0 1.3;
Tref 300; //K0*pow(T/Tref_, n0))
}
@ -104,7 +104,7 @@ pmmaCoeffs
radiation
{
sigmaS 0.0;
kappa 0.0;
kappaRad 0.0;
emissivity 0.1;
}
@ -126,7 +126,7 @@ charCoeffs
{
transport
{
K0 0.4;
kappa0 0.4;
n0 1.3;
Tref 300; //K0*pow(T/Tref_, n0))
}
@ -134,7 +134,7 @@ charCoeffs
radiation
{
sigmaS 0.0;
kappa 0.0;
kappaRad 0.0;
emissivity 0.0;
}

4
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/system/bottomAir/changeDictionaryDict
View File

@ -55,10 +55,10 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
K basicThermo;
kappa basicThermo;
QrNbr none;
Qr Qr;
KName none;
kappaName none;
value uniform 300;
}
}

12
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/system/heater/changeDictionaryDict
View File

@ -48,10 +48,10 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
K solidThermo;
kappa solidThermo;
QrNbr Qr;
Qr none;
KName none;
kappaName none;
value uniform 300;
}
@ -59,10 +59,10 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
K solidThermo;
kappa solidThermo;
QrNbr Qr;
Qr none;
KName none;
kappaName none;
value uniform 300;
}
@ -70,8 +70,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
value uniform 300;
}
minY

12
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/system/leftSolid/changeDictionaryDict
View File

@ -44,10 +44,10 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
K solidThermo;
kappa solidThermo;
QrNbr Qr;
Qr none;
KName none;
kappaName none;
value uniform 300;
}
@ -55,10 +55,10 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
K solidThermo;
kappa solidThermo;
QrNbr Qr;
Qr none;
KName none;
kappaName none;
value uniform 300;
}
@ -66,8 +66,8 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K solidThermo;
KName none;
kappa solidThermo;
kappaName none;
value uniform 300;
}
}

8
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/system/rightSolid/changeDictionaryDict
View File

@ -44,10 +44,10 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
K solidThermo;
kappa solidThermo;
QrNbr Qr;
Qr none;
KName none;
kappaName none;
value uniform 300;
}
@ -55,10 +55,10 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
K solidThermo;
kappa solidThermo;
QrNbr Qr;
Qr none;
KName none;
kappaName none;
value uniform 300;
}

4
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/system/topAir/changeDictionaryDict
View File

@ -66,10 +66,10 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
K basicThermo;
kappa basicThermo;
QrNbr none;
Qr Qr;
KName none;
kappaName none;
value uniform 300;
}
}