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

radiation: Corrected the name of the radiative heat flux from Qr to qr

Browse Source 
The standard naming convention for heat flux is "q" and this is used for the
conductive and convective heat fluxes is OpenFOAM.  The use of "Qr" for
radiative heat flux is an anomaly which causes confusion, particularly for
boundary conditions in which "Q" is used to denote power in Watts.  The name of
the radiative heat flux has now been corrected to "qr" and all models, boundary
conditions and tutorials updated.
This commit is contained in:
Henry Weller
2017-04-08 22:23:40 +01:00
parent 861b273e56
commit 5c62d81880
50 changed files with 283 additions and 284 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
applications/utilities/preProcessing/viewFactorsGen/viewFactorsGen.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -31,7 +31,7 @@ Description
Each view factor between the agglomerated faces i and j (Fij) is calculated
using a double integral of the sub-areas composing the agglomaration.
The patches involved in the view factor calculation are taken from the Qr
The patches involved in the view factor calculation are taken from the qr
volScalarField (radiative flux) when is greyDiffusiveRadiationViewFactor
otherwise they are not included.
@ -272,11 +272,11 @@ int main(int argc, char *argv[])
const label debug = viewFactorDict.lookupOrDefault<label>("debug", 0);
volScalarField Qr
volScalarField qr
(
IOobject
(
"Qr",
"qr",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
@ -338,17 +338,17 @@ int main(int argc, char *argv[])
labelList viewFactorsPatches(patches.size());
const volScalarField::Boundary& Qrb = Qr.boundaryField();
const volScalarField::Boundary& qrb = qr.boundaryField();
label count = 0;
forAll(Qrb, patchi)
forAll(qrb, patchi)
{
const polyPatch& pp = patches[patchi];
const fvPatchScalarField& QrpI = Qrb[patchi];
const fvPatchScalarField& qrpI = qrb[patchi];
if ((isA<fixedValueFvPatchScalarField>(QrpI)) && (pp.size() > 0))
if ((isA<fixedValueFvPatchScalarField>(qrpI)) && (pp.size() > 0))
{
viewFactorsPatches[count] = QrpI.patch().index();
viewFactorsPatches[count] = qrpI.patch().index();
nCoarseFaces += coarsePatches[patchi].size();
nFineFaces += patches[patchi].size();
count ++;

56
src/TurbulenceModels/compressible/turbulentFluidThermoModels/derivedFvPatchFields/thermalBaffle1D/thermalBaffle1DFvPatchScalarField.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -53,9 +53,9 @@ thermalBaffle1DFvPatchScalarField
Qs_(p.size()),
solidDict_(),
solidPtr_(nullptr),
QrPrevious_(p.size()),
QrRelaxation_(1),
QrName_("undefined-Qr")
qrPrevious_(p.size()),
qrRelaxation_(1),
qrName_("undefined-qr")
{}
@ -77,9 +77,9 @@ thermalBaffle1DFvPatchScalarField
Qs_(ptf.Qs_, mapper),
solidDict_(ptf.solidDict_),
solidPtr_(ptf.solidPtr_),
QrPrevious_(ptf.QrPrevious_, mapper),
QrRelaxation_(ptf.QrRelaxation_),
QrName_(ptf.QrName_)
qrPrevious_(ptf.qrPrevious_, mapper),
qrRelaxation_(ptf.qrRelaxation_),
qrName_(ptf.qrName_)
{}
@ -100,9 +100,9 @@ thermalBaffle1DFvPatchScalarField
Qs_(p.size(), 0),
solidDict_(dict),
solidPtr_(),
QrPrevious_(p.size(), 0.0),
QrRelaxation_(dict.lookupOrDefault<scalar>("relaxation", 1)),
QrName_(dict.lookupOrDefault<word>("Qr", "none"))
qrPrevious_(p.size(), 0.0),
qrRelaxation_(dict.lookupOrDefault<scalar>("relaxation", 1)),
qrName_(dict.lookupOrDefault<word>("qr", "none"))
{
fvPatchScalarField::operator=(scalarField("value", dict, p.size()));
@ -116,9 +116,9 @@ thermalBaffle1DFvPatchScalarField
Qs_ = scalarField("Qs", dict, p.size());
}
if (dict.found("QrPrevious"))
if (dict.found("qrPrevious"))
{
QrPrevious_ = scalarField("QrPrevious", dict, p.size());
qrPrevious_ = scalarField("qrPrevious", dict, p.size());
}
if (dict.found("refValue") && baffleActivated_)
@ -154,9 +154,9 @@ thermalBaffle1DFvPatchScalarField
Qs_(ptf.Qs_),
solidDict_(ptf.solidDict_),
solidPtr_(ptf.solidPtr_),
QrPrevious_(ptf.QrPrevious_),
QrRelaxation_(ptf.QrRelaxation_),
QrName_(ptf.QrName_)
qrPrevious_(ptf.qrPrevious_),
qrRelaxation_(ptf.qrRelaxation_),
qrName_(ptf.qrName_)
{}
@ -176,9 +176,9 @@ thermalBaffle1DFvPatchScalarField
Qs_(ptf.Qs_),
solidDict_(ptf.solidDict_),
solidPtr_(ptf.solidPtr_),
QrPrevious_(ptf.QrPrevious_),
QrRelaxation_(ptf.QrRelaxation_),
QrName_(ptf.QrName_)
qrPrevious_(ptf.qrPrevious_),
qrRelaxation_(ptf.qrRelaxation_),
qrName_(ptf.qrName_)
{}
@ -364,15 +364,15 @@ void thermalBaffle1DFvPatchScalarField<solidType>::updateCoeffs()
patch().template lookupPatchField<volScalarField, scalar>(TName_);
scalarField Qr(Tp.size(), 0.0);
scalarField qr(Tp.size(), 0.0);
if (QrName_ != "none")
if (qrName_ != "none")
{
Qr = patch().template lookupPatchField<volScalarField, scalar>
(QrName_);
qr = patch().template lookupPatchField<volScalarField, scalar>
(qrName_);
Qr = QrRelaxation_*Qr + (1.0 - QrRelaxation_)*QrPrevious_;
QrPrevious_ = Qr;
qr = qrRelaxation_*qr + (1.0 - qrRelaxation_)*qrPrevious_;
qrPrevious_ = qr;
}
tmp<scalarField> Ti = patchInternalField();
@ -393,7 +393,7 @@ void thermalBaffle1DFvPatchScalarField<solidType>::updateCoeffs()
scalarField KDeltaSolid(kappas/baffleThickness());
scalarField alpha(KDeltaSolid - Qr/Tp);
scalarField alpha(KDeltaSolid - qr/Tp);
valueFraction() = alpha/(alpha + myKDelta);
@ -435,9 +435,9 @@ void thermalBaffle1DFvPatchScalarField<solidType>::write(Ostream& os) const
solid().write(os);
}
QrPrevious_.writeEntry("QrPrevious", os);
os.writeKeyword("Qr")<< QrName_ << token::END_STATEMENT << nl;
os.writeKeyword("relaxation")<< QrRelaxation_
qrPrevious_.writeEntry("qrPrevious", os);
os.writeKeyword("qr")<< qrName_ << token::END_STATEMENT << nl;
os.writeKeyword("relaxation")<< qrRelaxation_
<< token::END_STATEMENT << nl;
}

18
src/TurbulenceModels/compressible/turbulentFluidThermoModels/derivedFvPatchFields/thermalBaffle1D/thermalBaffle1DFvPatchScalarField.H
View File

@ -31,8 +31,8 @@ Description
This BC solves a steady 1D thermal baffle.
The solid properties are specify as dictionary. Optionaly radiative heat
flux (Qr) can be incorporated into the balance. Some under-relaxation might
be needed on Qr. Baffle and solid properties need to be specified on the
flux (qr) can be incorporated into the balance. Some under-relaxation might
be needed on qr. Baffle and solid properties need to be specified on the
master side of the baffle.
Usage
@ -48,7 +48,7 @@ Usage
thickness uniform 0.005; // thickness [m]
Qs uniform 100; // heat flux [W/m2]
Qr none;
qr none;
relaxation 1;
// Solid thermo
@ -78,7 +78,7 @@ Usage
type compressible::thermalBaffle1D<hConstSolidThermoPhysics>;
samplePatch <masterPatchName>;
Qr none;
qr none;
relaxation 1;
}
\endverbatim
@ -132,14 +132,14 @@ class thermalBaffle1DFvPatchScalarField
//- Solid thermo
mutable autoPtr<solidType> solidPtr_;
//- Cache Qr for relaxation
scalarField QrPrevious_;
//- Cache qr for relaxation
scalarField qrPrevious_;
//- Relaxation for Qr
scalar QrRelaxation_;
//- Relaxation for qr
scalar qrRelaxation_;
//- Name of the radiative heat flux in local region
const word QrName_;
const word qrName_;
// Private members

38
src/TurbulenceModels/compressible/turbulentFluidThermoModels/derivedFvPatchFields/turbulentTemperatureRadCoupledMixed/turbulentTemperatureRadCoupledMixedFvPatchScalarField.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -48,8 +48,8 @@ turbulentTemperatureRadCoupledMixedFvPatchScalarField
mixedFvPatchScalarField(p, iF),
temperatureCoupledBase(patch(), "undefined", "undefined", "undefined-K"),
TnbrName_("undefined-Tnbr"),
QrNbrName_("undefined-QrNbr"),
QrName_("undefined-Qr"),
qrNbrName_("undefined-qrNbr"),
qrName_("undefined-qr"),
thicknessLayers_(0),
kappaLayers_(0),
contactRes_(0)
@ -72,8 +72,8 @@ turbulentTemperatureRadCoupledMixedFvPatchScalarField
mixedFvPatchScalarField(psf, p, iF, mapper),
temperatureCoupledBase(patch(), psf),
TnbrName_(psf.TnbrName_),
QrNbrName_(psf.QrNbrName_),
QrName_(psf.QrName_),
qrNbrName_(psf.qrNbrName_),
qrName_(psf.qrName_),
thicknessLayers_(psf.thicknessLayers_),
kappaLayers_(psf.kappaLayers_),
contactRes_(psf.contactRes_)
@ -91,8 +91,8 @@ turbulentTemperatureRadCoupledMixedFvPatchScalarField
mixedFvPatchScalarField(p, iF),
temperatureCoupledBase(patch(), dict),
TnbrName_(dict.lookupOrDefault<word>("Tnbr", "T")),
QrNbrName_(dict.lookupOrDefault<word>("QrNbr", "none")),
QrName_(dict.lookupOrDefault<word>("Qr", "none")),
qrNbrName_(dict.lookupOrDefault<word>("qrNbr", "none")),
qrName_(dict.lookupOrDefault<word>("qr", "none")),
thicknessLayers_(0),
kappaLayers_(0),
contactRes_(0.0)
@ -152,8 +152,8 @@ turbulentTemperatureRadCoupledMixedFvPatchScalarField
mixedFvPatchScalarField(psf, iF),
temperatureCoupledBase(patch(), psf),
TnbrName_(psf.TnbrName_),
QrNbrName_(psf.QrNbrName_),
QrName_(psf.QrName_),
qrNbrName_(psf.qrNbrName_),
qrName_(psf.qrName_),
thicknessLayers_(psf.thicknessLayers_),
kappaLayers_(psf.kappaLayers_),
contactRes_(psf.contactRes_)
@ -212,22 +212,22 @@ void turbulentTemperatureRadCoupledMixedFvPatchScalarField::updateCoeffs()
scalarField KDelta(kappa(Tp)*patch().deltaCoeffs());
scalarField Qr(Tp.size(), 0.0);
if (QrName_ != "none")
scalarField qr(Tp.size(), 0.0);
if (qrName_ != "none")
{
Qr = patch().lookupPatchField<volScalarField, scalar>(QrName_);
qr = patch().lookupPatchField<volScalarField, scalar>(qrName_);
}
scalarField QrNbr(Tp.size(), 0.0);
if (QrNbrName_ != "none")
scalarField qrNbr(Tp.size(), 0.0);
if (qrNbrName_ != "none")
{
QrNbr = nbrPatch.lookupPatchField<volScalarField, scalar>(QrNbrName_);
mpp.distribute(QrNbr);
qrNbr = nbrPatch.lookupPatchField<volScalarField, scalar>(qrNbrName_);
mpp.distribute(qrNbr);
}
valueFraction() = KDeltaNbr/(KDeltaNbr + KDelta);
refValue() = TcNbr;
refGrad() = (Qr + QrNbr)/kappa(Tp);
refGrad() = (qr + qrNbr)/kappa(Tp);
mixedFvPatchScalarField::updateCoeffs();
@ -261,8 +261,8 @@ void turbulentTemperatureRadCoupledMixedFvPatchScalarField::write
{
mixedFvPatchScalarField::write(os);
os.writeKeyword("Tnbr")<< TnbrName_ << token::END_STATEMENT << nl;
os.writeKeyword("QrNbr")<< QrNbrName_ << token::END_STATEMENT << nl;
os.writeKeyword("Qr")<< QrName_ << token::END_STATEMENT << nl;
os.writeKeyword("qrNbr")<< qrNbrName_ << token::END_STATEMENT << nl;
os.writeKeyword("qr")<< qrName_ << token::END_STATEMENT << nl;
thicknessLayers_.writeEntry("thicknessLayers", os);
kappaLayers_.writeEntry("kappaLayers", os);

14
src/TurbulenceModels/compressible/turbulentFluidThermoModels/derivedFvPatchFields/turbulentTemperatureRadCoupledMixed/turbulentTemperatureRadCoupledMixedFvPatchScalarField.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -38,8 +38,8 @@ Usage
\table
Property | Description | Required | Default value
Tnbr | name of the field | no | T
QrNbr | name of the radiative flux in the nbr region | no | none
Qr | name of the radiative flux in this region | no | none
qrNbr | name of the radiative flux in the nbr region | no | none
qr | name of the radiative flux in this region | no | none
thicknessLayers | list of thicknesses per layer [m] | no |
kappaLayers | list of thermal conductivites per layer [W/m/K] | no |
kappaMethod | inherited from temperatureCoupledBase | inherited |
@ -52,8 +52,8 @@ Usage
{
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
QrNbr Qr; // or none. Name of Qr field on neighbour region
Qr Qr; // or none. Name of Qr field on local region
qrNbr qr; // or none. Name of qr field on neighbour region
qr qr; // or none. Name of qr field on local region
thicknessLayers (0.1 0.2 0.3 0.4);
kappaLayers (1 2 3 4);
kappaMethod lookup;
@ -101,10 +101,10 @@ class turbulentTemperatureRadCoupledMixedFvPatchScalarField
const word TnbrName_;
//- Name of the radiative heat flux in the neighbout region
const word QrNbrName_;
const word qrNbrName_;
//- Name of the radiative heat flux in local region
const word QrName_;
const word qrName_;
//- Thickness of layers
scalarList thicknessLayers_;

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

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2016-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -78,12 +78,12 @@ void Foam::functionObjects::wallHeatFlux::calcHeatFlux
wallHeatFluxBf[patchi] = heatFluxBf[patchi];
}
if (foundObject<volScalarField>("Qr"))
if (foundObject<volScalarField>("qr"))
{
const volScalarField& Qr = lookupObject<volScalarField>("Qr");
const volScalarField& qr = lookupObject<volScalarField>("qr");
const volScalarField::Boundary& radHeatFluxBf =
Qr.boundaryField();
qr.boundaryField();
forAll(wallHeatFluxBf, patchi)
{

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

@ -58,7 +58,7 @@ void reactingOneDim::readReactingOneDimControls()
coeffs().lookup("minimumDelta") >> minimumDelta_;
coeffs().lookup("gasHSource") >> gasHSource_;
coeffs().lookup("QrHSource") >> QrHSource_;
coeffs().lookup("qrHSource") >> qrHSource_;
useChemistrySolvers_ =
coeffs().lookupOrDefault<bool>("useChemistrySolvers", true);
}
@ -92,41 +92,41 @@ bool reactingOneDim::read(const dictionary& dict)
}
void reactingOneDim::updateQr()
void reactingOneDim::updateqr()
{
// Update local Qr from coupled Qr field
Qr_ == dimensionedScalar("zero", Qr_.dimensions(), 0.0);
// Update local qr from coupled qr field
qr_ == dimensionedScalar("zero", qr_.dimensions(), 0.0);
// Retrieve field from coupled region using mapped boundary conditions
Qr_.correctBoundaryConditions();
qr_.correctBoundaryConditions();
volScalarField::Boundary& QrBf = Qr_.boundaryFieldRef();
volScalarField::Boundary& qrBf = qr_.boundaryFieldRef();
forAll(intCoupledPatchIDs_, i)
{
const label patchi = intCoupledPatchIDs_[i];
// Qr is positive going in the solid
// qr is positive going in the solid
// If the surface is emitting the radiative flux is set to zero
QrBf[patchi] = max(QrBf[patchi], scalar(0.0));
qrBf[patchi] = max(qrBf[patchi], scalar(0.0));
}
const vectorField& cellC = regionMesh().cellCentres();
tmp<volScalarField> kappa = kappaRad();
// Propagate Qr through 1-D regions
// Propagate qr through 1-D regions
label localPyrolysisFacei = 0;
forAll(intCoupledPatchIDs_, i)
{
const label patchi = intCoupledPatchIDs_[i];
const scalarField& Qrp = Qr_.boundaryField()[patchi];
const scalarField& qrp = qr_.boundaryField()[patchi];
const vectorField& Cf = regionMesh().Cf().boundaryField()[patchi];
forAll(Qrp, facei)
forAll(qrp, facei)
{
const scalar Qr0 = Qrp[facei];
const scalar qr0 = qrp[facei];
point Cf0 = Cf[facei];
const labelList& cells = boundaryFaceCells_[localPyrolysisFacei++];
scalar kappaInt = 0.0;
@ -136,7 +136,7 @@ void reactingOneDim::updateQr()
const point& Cf1 = cellC[celli];
const scalar delta = mag(Cf1 - Cf0);
kappaInt += kappa()[celli]*delta;
Qr_[celli] = Qr0*exp(-kappaInt);
qr_[celli] = qr0*exp(-kappaInt);
Cf0 = Cf1;
}
}
@ -202,9 +202,9 @@ void reactingOneDim::updatePhiGas()
void reactingOneDim::updateFields()
{
if (QrHSource_)
if (qrHSource_)
{
updateQr();
updateqr();
}
updatePhiGas();
@ -332,10 +332,10 @@ void reactingOneDim::solveEnergy()
hEqn += fvc::div(phiGas);
}
if (QrHSource_)
if (qrHSource_)
{
const surfaceScalarField phiQr(fvc::interpolate(Qr_)*nMagSf());
hEqn += fvc::div(phiQr);
const surfaceScalarField phiqr(fvc::interpolate(qr_)*nMagSf());
hEqn += fvc::div(phiqr);
}
if (regionMesh().moving())
@ -447,11 +447,11 @@ reactingOneDim::reactingOneDim
dimensionedScalar("zero", dimEnergy/dimTime/dimVolume, 0.0)
),
Qr_
qr_
(
IOobject
(
"Qr",
"qr",
time().timeName(),
regionMesh(),
IOobject::MUST_READ,
@ -465,7 +465,7 @@ reactingOneDim::reactingOneDim
totalGasMassFlux_(0.0),
totalHeatRR_(dimensionedScalar("zero", dimEnergy/dimTime, 0.0)),
gasHSource_(false),
QrHSource_(false),
qrHSource_(false),
useChemistrySolvers_(true)
{
if (active_)
@ -547,11 +547,11 @@ reactingOneDim::reactingOneDim
dimensionedScalar("zero", dimEnergy/dimTime/dimVolume, 0.0)
),
Qr_
qr_
(
IOobject
(
"Qr",
"qr",
time().timeName(),
regionMesh(),
IOobject::MUST_READ,
@ -565,7 +565,7 @@ reactingOneDim::reactingOneDim
totalGasMassFlux_(0.0),
totalHeatRR_(dimensionedScalar("zero", dimEnergy/dimTime, 0.0)),
gasHSource_(false),
QrHSource_(false),
qrHSource_(false),
useChemistrySolvers_(true)
{
if (active_)

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

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -125,10 +125,10 @@ protected:
//- Coupled region radiative heat flux [W/m2]
// Requires user to input mapping info for coupled patches
//volScalarField QrCoupled_;
//volScalarField qrCoupled_;
//- In depth radiative heat flux [W/m2]
volScalarField Qr_;
volScalarField qr_;
// Checks
@ -152,7 +152,7 @@ protected:
bool gasHSource_;
//- Add in depth radiation source term
bool QrHSource_;
bool qrHSource_;
//- Use chemistry solvers (ode or sequential)
bool useChemistrySolvers_;
@ -173,7 +173,7 @@ protected:
void updateMesh(const scalarField& mass0);
//- Update radiative flux in pyrolysis region
void updateQr();
void updateqr();
//- Update enthalpy flux for pyrolysis gases
void updatePhiGas();

32
src/regionModels/regionCoupling/derivedFvPatchFields/filmPyrolysisRadiativeCoupledMixed/filmPyrolysisRadiativeCoupledMixedFvPatchScalarField.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2013-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2013-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -110,7 +110,7 @@ filmPyrolysisRadiativeCoupledMixedFvPatchScalarField
filmRegionName_("surfaceFilmProperties"),
pyrolysisRegionName_("pyrolysisProperties"),
TnbrName_("undefined-Tnbr"),
QrName_("undefined-Qr"),
qrName_("undefined-qr"),
convectiveScaling_(1.0),
filmDeltaDry_(0.0),
filmDeltaWet_(0.0)
@ -135,7 +135,7 @@ filmPyrolysisRadiativeCoupledMixedFvPatchScalarField
filmRegionName_(psf.filmRegionName_),
pyrolysisRegionName_(psf.pyrolysisRegionName_),
TnbrName_(psf.TnbrName_),
QrName_(psf.QrName_),
qrName_(psf.qrName_),
convectiveScaling_(psf.convectiveScaling_),
filmDeltaDry_(psf.filmDeltaDry_),
filmDeltaWet_(psf.filmDeltaWet_)
@ -161,7 +161,7 @@ filmPyrolysisRadiativeCoupledMixedFvPatchScalarField
dict.lookupOrDefault<word>("pyrolysisRegion", "pyrolysisProperties")
),
TnbrName_(dict.lookup("Tnbr")),
QrName_(dict.lookup("Qr")),
qrName_(dict.lookup("qr")),
convectiveScaling_(dict.lookupOrDefault<scalar>("convectiveScaling", 1.0)),
filmDeltaDry_(readScalar(dict.lookup("filmDeltaDry"))),
filmDeltaWet_(readScalar(dict.lookup("filmDeltaWet")))
@ -207,7 +207,7 @@ filmPyrolysisRadiativeCoupledMixedFvPatchScalarField
filmRegionName_(psf.filmRegionName_),
pyrolysisRegionName_(psf.pyrolysisRegionName_),
TnbrName_(psf.TnbrName_),
QrName_(psf.QrName_),
qrName_(psf.qrName_),
convectiveScaling_(psf.convectiveScaling_),
filmDeltaDry_(psf.filmDeltaDry_),
filmDeltaWet_(psf.filmDeltaWet_)
@ -268,25 +268,25 @@ void filmPyrolysisRadiativeCoupledMixedFvPatchScalarField::updateCoeffs()
const pyrolysisModelType& pyrolysis = pyrModel();
const filmModelType& film = filmModel();
// Obtain Rad heat (Qr)
scalarField Qr(patch().size(), 0.0);
// Obtain Rad heat (qr)
scalarField qr(patch().size(), 0.0);
label coupledPatchi = -1;
if (pyrolysisRegionName_ == mesh.name())
{
coupledPatchi = patchi;
if (QrName_ != "none")
if (qrName_ != "none")
{
Qr = nbrPatch.lookupPatchField<volScalarField, scalar>(QrName_);
mpp.distribute(Qr);
qr = nbrPatch.lookupPatchField<volScalarField, scalar>(qrName_);
mpp.distribute(qr);
}
}
else if (pyrolysis.primaryMesh().name() == mesh.name())
{
coupledPatchi = nbrPatch.index();
if (QrName_ != "none")
if (qrName_ != "none")
{
Qr = patch().lookupPatchField<volScalarField, scalar>(QrName_);
qr = patch().lookupPatchField<volScalarField, scalar>(qrName_);
}
}
else
@ -344,7 +344,7 @@ void filmPyrolysisRadiativeCoupledMixedFvPatchScalarField::updateCoeffs()
scalarField qConv(ratio*htcwfilm*(Tfilm - Tp)*convectiveScaling_);
scalarField qRad((1.0 - ratio)*Qr);
scalarField qRad((1.0 - ratio)*qr);
scalarField alpha(KDeltaNbr - (qRad + qConv)/Tp);
@ -357,7 +357,7 @@ void filmPyrolysisRadiativeCoupledMixedFvPatchScalarField::updateCoeffs()
if (debug)
{
scalar Qc = gSum(qConv*patch().magSf());
scalar Qr = gSum(qRad*patch().magSf());
scalar qr = gSum(qRad*patch().magSf());
scalar Qt = gSum((qConv + qRad)*patch().magSf());
Info<< mesh.name() << ':'
@ -367,7 +367,7 @@ void filmPyrolysisRadiativeCoupledMixedFvPatchScalarField::updateCoeffs()
<< nbrPatch.name() << ':'
<< this->internalField().name() << " :" << nl
<< " convective heat[W] : " << Qc << nl
<< " radiative heat [W] : " << Qr << nl
<< " radiative heat [W] : " << qr << nl
<< " total heat [W] : " << Qt << nl
<< " wall temperature "
<< " min:" << gMin(*this)
@ -399,7 +399,7 @@ void filmPyrolysisRadiativeCoupledMixedFvPatchScalarField::write
pyrolysisRegionName_
);
os.writeKeyword("Tnbr")<< TnbrName_ << token::END_STATEMENT << nl;
os.writeKeyword("Qr")<< QrName_ << token::END_STATEMENT << nl;
os.writeKeyword("qr")<< qrName_ << token::END_STATEMENT << nl;
os.writeKeyword("convectiveScaling") << convectiveScaling_
<< token::END_STATEMENT << nl;
os.writeKeyword("filmDeltaDry") << filmDeltaDry_

10
src/regionModels/regionCoupling/derivedFvPatchFields/filmPyrolysisRadiativeCoupledMixed/filmPyrolysisRadiativeCoupledMixedFvPatchScalarField.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2013-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2013-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -35,7 +35,7 @@ Description
type filmPyrolysisRadiativeCoupledMixed;
Tnbr T;
kappa fluidThermo;
Qr Qr;
qr qr;
kappaName none;
filmDeltaDry 0.0;
filmDeltaWet 3e-4;
@ -58,7 +58,7 @@ Description
\verbatim
qConv = ratio*htcwfilm*(Tfilm - *this);
qRad = (1.0 - ratio)*Qr;
qRad = (1.0 - ratio)*qr;
\endverbatim
Then the solid can gain or loose energy through radiation or conduction
@ -67,7 +67,7 @@ Description
Notes:
- kappa and \c kappaName are inherited from temperatureCoupledBase.
- Qr is the radiative flux defined in the radiation model.
- qr is the radiative flux defined in the radiation model.
See also
@ -124,7 +124,7 @@ private:
const word TnbrName_;
//- Name of the radiative heat flux
const word QrName_;
const word qrName_;
//- Convective Scaling Factor (as determined by Prateep's tests)
const scalar convectiveScaling_;

8
src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/constantRadiation/constantRadiation.C
View File

@ -56,11 +56,11 @@ constantRadiation::constantRadiation
)
:
filmRadiationModel(typeName, film, dict),
QrConst_
qrConst_
(
IOobject
(
typeName + ":QrConst",
typeName + ":qrConst",
film.time().timeName(),
film.regionMesh(),
IOobject::MUST_READ,
@ -125,10 +125,10 @@ tmp<volScalarField> constantRadiation::Shs()
if ((time >= timeStart_) && (time <= timeStart_ + duration_))
{
scalarField& Shs = tShs.ref();
const scalarField& Qr = QrConst_;
const scalarField& qr = qrConst_;
const scalarField& alpha = filmModel_.alpha();
Shs = mask_*Qr*alpha*absorptivity_;
Shs = mask_*qr*alpha*absorptivity_;
}
return tShs;

4
src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/constantRadiation/constantRadiation.H
View File

@ -63,7 +63,7 @@ private:
// Private data
//- Constant radiative flux [kg/s3]
volScalarField QrConst_;
volScalarField qrConst_;
//- Radiation mask
volScalarField mask_;
@ -115,7 +115,7 @@ public:
virtual void correct();
//- Return the radiation sensible enthalpy source
// Also updates QrNet
// Also updates qrNet
virtual tmp<volScalarField> Shs();
};

2
src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/primaryRadiation/primaryRadiation.H
View File

@ -101,7 +101,7 @@ public:
virtual void correct();
//- Return the radiation sensible enthalpy source
// Also updates QrNet
// Also updates qrNet
virtual tmp<volScalarField> Shs();
};

12
src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/standardRadiation/standardRadiation.C
View File

@ -71,11 +71,11 @@ standardRadiation::standardRadiation
dimensionedScalar("zero", dimMass/pow3(dimTime), 0.0),
film.mappedPushedFieldPatchTypes<scalar>()
),
QrNet_
qrNet_
(
IOobject
(
"QrNet",
"qrNet",
film.time().timeName(),
film.regionMesh(),
IOobject::NO_READ,
@ -100,7 +100,7 @@ standardRadiation::~standardRadiation()
void standardRadiation::correct()
{
// Transfer Qr from primary region
// Transfer qr from primary region
QinPrimary_.correctBoundaryConditions();
}
@ -131,9 +131,9 @@ tmp<volScalarField> standardRadiation::Shs()
Shs = beta_*QinP*alpha*(1.0 - exp(-kappaBar_*delta));
// Update net Qr on local region
QrNet_.primitiveFieldRef() = QinP - Shs;
QrNet_.correctBoundaryConditions();
// Update net qr on local region
qrNet_.primitiveFieldRef() = QinP - Shs;
qrNet_.correctBoundaryConditions();
return tShs;
}

4
src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/standardRadiation/standardRadiation.H
View File

@ -63,7 +63,7 @@ private:
volScalarField QinPrimary_;
//- Remaining radiative flux after removing local contribution
volScalarField QrNet_;
volScalarField qrNet_;
// Model coefficients
@ -112,7 +112,7 @@ public:
virtual void correct();
//- Return the radiation sensible enthalpy source
// Also updates QrNet
// Also updates qrNet
virtual tmp<volScalarField> Shs();
};

4
src/regionModels/thermalBaffleModels/derivedFvPatchFields/thermalBaffle/thermalBaffleFvPatchScalarField.H
View File

@ -54,8 +54,8 @@ Usage
Tnbr T;
kappa fluidThermo; // or solidThermo
KappaName none;
QrNbr Qr;//or none.Name of Qr field on neighbourregion
Qr none;// or none.Name of Qr field on localregion
qrNbr qr;//or none.Name of qr field on neighbourregion
qr none;// or none.Name of qr field on localregion
value uniform 300;
// Baffle region name

4
src/thermophysicalModels/radiation/derivedFvPatchFields/greyDiffusiveRadiation/greyDiffusiveRadiationMixedFvPatchScalarField.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -185,7 +185,7 @@ updateCoeffs()
const scalarField nAve(n & ray.dAve());
ray.Qr().boundaryFieldRef()[patchi] += Iw*nAve;
ray.qr().boundaryFieldRef()[patchi] += Iw*nAve;
const scalarField temissivity = emissivity();

14
src/thermophysicalModels/radiation/derivedFvPatchFields/greyDiffusiveViewFactor/greyDiffusiveViewFactorFixedValueFvPatchScalarField.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -40,7 +40,7 @@ greyDiffusiveViewFactorFixedValueFvPatchScalarField
:
fixedValueFvPatchScalarField(p, iF),
radiationCoupledBase(patch(), "undefined", scalarField::null()),
Qro_(p.size(), 0.0)
qro_(p.size(), 0.0)
{}
@ -60,7 +60,7 @@ greyDiffusiveViewFactorFixedValueFvPatchScalarField
ptf.emissivityMethod(),
ptf.emissivity_
),
Qro_(ptf.Qro_)
qro_(ptf.qro_)
{}
@ -74,7 +74,7 @@ greyDiffusiveViewFactorFixedValueFvPatchScalarField
:
fixedValueFvPatchScalarField(p, iF, dict, false),
radiationCoupledBase(p, dict),
Qro_("Qro", dict, p.size())
qro_("qro", dict, p.size())
{
if (dict.found("value"))
{
@ -104,7 +104,7 @@ greyDiffusiveViewFactorFixedValueFvPatchScalarField
ptf.emissivityMethod(),
ptf.emissivity_
),
Qro_(ptf.Qro_)
qro_(ptf.qro_)
{}
@ -122,7 +122,7 @@ greyDiffusiveViewFactorFixedValueFvPatchScalarField
ptf.emissivityMethod(),
ptf.emissivity_
),
Qro_(ptf.Qro_)
qro_(ptf.qro_)
{}
@ -163,7 +163,7 @@ write
{
fixedValueFvPatchScalarField::write(os);
radiationCoupledBase::write(os);
Qro_.writeEntry("Qro", os);
qro_.writeEntry("qro", os);
}

14
src/thermophysicalModels/radiation/derivedFvPatchFields/greyDiffusiveViewFactor/greyDiffusiveViewFactorFixedValueFvPatchScalarField.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -29,12 +29,12 @@ Group
Description
This boundary condition provides a grey-diffuse condition for radiative
heat flux, \c Qr, for use with the view factor model
heat flux, \c qr, for use with the view factor model
Usage
\table
Property | Description | Required | Default value
Qro | external radiative heat flux | yes |
qro | external radiative heat flux | yes |
emissivityMode | emissivity mode: solidRadiation or lookup | yes |
\endtable
@ -43,7 +43,7 @@ Usage
<patchName>
{
type greyDiffusiveRadiationViewFactor;
Qro uniform 0;
qro uniform 0;
emissivityMode solidRadiation;
value uniform 0;
}
@ -85,7 +85,7 @@ class greyDiffusiveViewFactorFixedValueFvPatchScalarField
// Private data
//- External radiative heat flux
scalarField Qro_;
scalarField qro_;
public:
@ -165,9 +165,9 @@ public:
// Access
//- Return external radiative heat flux
const scalarList& Qro()
const scalarList& qro()
{
return Qro_;
return qro_;
}

4
src/thermophysicalModels/radiation/derivedFvPatchFields/wideBandDiffusiveRadiation/wideBandDiffusiveRadiationMixedFvPatchScalarField.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -186,7 +186,7 @@ updateCoeffs()
const scalarField nAve(n & ray.dAve());
ray.Qr().boundaryFieldRef()[patchi] += Iw*nAve;
ray.qr().boundaryFieldRef()[patchi] += Iw*nAve;
const scalarField Eb
(

18
src/thermophysicalModels/radiation/radiationModels/P1/P1.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -62,18 +62,18 @@ Foam::radiation::P1::P1(const volScalarField& T)
),
mesh_
),
Qr_
qr_
(
IOobject
(
"Qr",
"qr",
mesh_.time().timeName(),
mesh_,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh_,
dimensionedScalar("Qr", dimMass/pow3(dimTime), 0.0)
dimensionedScalar("qr", dimMass/pow3(dimTime), 0.0)
),
a_
(
@ -132,18 +132,18 @@ Foam::radiation::P1::P1(const dictionary& dict, const volScalarField& T)
),
mesh_
),
Qr_
qr_
(
IOobject
(
"Qr",
"qr",
mesh_.time().timeName(),
mesh_,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh_,
dimensionedScalar("Qr", dimMass/pow3(dimTime), 0.0)
dimensionedScalar("qr", dimMass/pow3(dimTime), 0.0)
),
a_
(
@ -242,14 +242,14 @@ void Foam::radiation::P1::calculate()
- 4.0*(e_*physicoChemical::sigma*pow4(T_) ) - E_
);
volScalarField::Boundary& QrBf = Qr_.boundaryFieldRef();
volScalarField::Boundary& qrBf = qr_.boundaryFieldRef();
// Calculate radiative heat flux on boundaries.
forAll(mesh_.boundaryMesh(), patchi)
{
if (!G_.boundaryField()[patchi].coupled())
{
QrBf[patchi] =
qrBf[patchi] =
-gamma.boundaryField()[patchi]
*G_.boundaryField()[patchi].snGrad();
}

4
src/thermophysicalModels/radiation/radiationModels/P1/P1.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -65,7 +65,7 @@ class P1
volScalarField G_;
//- Total radiative heat flux [W/m2]
volScalarField Qr_;
volScalarField qr_;
//- Absorption coefficient
volScalarField a_;

18
src/thermophysicalModels/radiation/radiationModels/fvDOM/fvDOM/fvDOM.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -222,18 +222,18 @@ Foam::radiation::fvDOM::fvDOM(const volScalarField& T)
mesh_,
dimensionedScalar("G", dimMass/pow3(dimTime), 0.0)
),
Qr_
qr_
(
IOobject
(
"Qr",
"qr",
mesh_.time().timeName(),
mesh_,
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
),
mesh_,
dimensionedScalar("Qr", dimMass/pow3(dimTime), 0.0)
dimensionedScalar("qr", dimMass/pow3(dimTime), 0.0)
),
Qem_
(
@ -309,18 +309,18 @@ Foam::radiation::fvDOM::fvDOM
mesh_,
dimensionedScalar("G", dimMass/pow3(dimTime), 0.0)
),
Qr_
qr_
(
IOobject
(
"Qr",
"qr",
mesh_.time().timeName(),
mesh_,
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
),
mesh_,
dimensionedScalar("Qr", dimMass/pow3(dimTime), 0.0)
dimensionedScalar("qr", dimMass/pow3(dimTime), 0.0)
),
Qem_
(
@ -490,7 +490,7 @@ void Foam::radiation::fvDOM::updateBlackBodyEmission()
void Foam::radiation::fvDOM::updateG()
{
G_ = dimensionedScalar("zero",dimMass/pow3(dimTime), 0.0);
Qr_ = dimensionedScalar("zero",dimMass/pow3(dimTime), 0.0);
qr_ = dimensionedScalar("zero",dimMass/pow3(dimTime), 0.0);
Qem_ = dimensionedScalar("zero", dimMass/pow3(dimTime), 0.0);
Qin_ = dimensionedScalar("zero", dimMass/pow3(dimTime), 0.0);
@ -498,7 +498,7 @@ void Foam::radiation::fvDOM::updateG()
{
IRay_[rayI].addIntensity();
G_ += IRay_[rayI].I()*IRay_[rayI].omega();
Qr_.boundaryFieldRef() += IRay_[rayI].Qr().boundaryField();
qr_.boundaryFieldRef() += IRay_[rayI].qr().boundaryField();
Qem_.boundaryFieldRef() += IRay_[rayI].Qem().boundaryField();
Qin_.boundaryFieldRef() += IRay_[rayI].Qin().boundaryField();
}

6
src/thermophysicalModels/radiation/radiationModels/fvDOM/fvDOM/fvDOM.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -89,7 +89,7 @@ class fvDOM
volScalarField G_;
//- Total radiative heat flux [W/m2]
volScalarField Qr_;
volScalarField qr_;
//- Emmited radiative heat flux [W/m2]
volScalarField Qem_;
@ -228,7 +228,7 @@ public:
inline const volScalarField& G() const;
//- Const access to total radiative heat flux field
inline const volScalarField& Qr() const;
inline const volScalarField& qr() const;
//- Const access to incident radiative heat flux field
inline const volScalarField& Qin() const;

6
src/thermophysicalModels/radiation/radiationModels/fvDOM/fvDOM/fvDOMI.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -86,9 +86,9 @@ inline const Foam::volScalarField& Foam::radiation::fvDOM::G() const
}
inline const Foam::volScalarField& Foam::radiation::fvDOM::Qr() const
inline const Foam::volScalarField& Foam::radiation::fvDOM::qr() const
{
return Qr_;
return qr_;
}
inline const Foam::volScalarField& Foam::radiation::fvDOM::Qin() const

10
src/thermophysicalModels/radiation/radiationModels/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -69,18 +69,18 @@ Foam::radiation::radiativeIntensityRay::radiativeIntensityRay
mesh_,
dimensionedScalar("I", dimMass/pow3(dimTime), 0.0)
),
Qr_
qr_
(
IOobject
(
"Qr" + name(rayId),
"qr" + name(rayId),
mesh_.time().timeName(),
mesh_,
IOobject::NO_READ,
IOobject::NO_WRITE
),
mesh_,
dimensionedScalar("Qr", dimMass/pow3(dimTime), 0.0)
dimensionedScalar("qr", dimMass/pow3(dimTime), 0.0)
),
Qin_
(
@ -206,7 +206,7 @@ Foam::radiation::radiativeIntensityRay::~radiativeIntensityRay()
Foam::scalar Foam::radiation::radiativeIntensityRay::correct()
{
// Reset boundary heat flux to zero
Qr_.boundaryFieldRef() = 0.0;
qr_.boundaryFieldRef() = 0.0;
scalar maxResidual = -GREAT;

8
src/thermophysicalModels/radiation/radiationModels/fvDOM/radiativeIntensityRay/radiativeIntensityRay.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -80,7 +80,7 @@ private:
volScalarField I_;
//- Total radiative heat flux on boundary
volScalarField Qr_;
volScalarField qr_;
//- Incident radiative heat flux on boundary
volScalarField Qin_;
@ -176,10 +176,10 @@ public:
inline const volScalarField& I() const;
//- Return const access to the boundary heat flux
inline const volScalarField& Qr() const;
inline const volScalarField& qr() const;
//- Return non-const access to the boundary heat flux
inline volScalarField& Qr();
inline volScalarField& qr();
//- Return non-const access to the boundary incident heat flux
inline volScalarField& Qin();

10
src/thermophysicalModels/radiation/radiationModels/fvDOM/radiativeIntensityRay/radiativeIntensityRayI.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -31,15 +31,15 @@ Foam::radiation::radiativeIntensityRay::I() const
inline const Foam::volScalarField&
Foam::radiation::radiativeIntensityRay::Qr() const
Foam::radiation::radiativeIntensityRay::qr() const
{
return Qr_;
return qr_;
}
inline Foam::volScalarField& Foam::radiation::radiativeIntensityRay::Qr()
inline Foam::volScalarField& Foam::radiation::radiativeIntensityRay::qr()
{
return Qr_;
return qr_;
}
inline const Foam::volScalarField& Foam::radiation::

64
src/thermophysicalModels/radiation/radiationModels/viewFactor/viewFactor.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -49,17 +49,17 @@ namespace Foam
void Foam::radiation::viewFactor::initialise()
{
const polyBoundaryMesh& coarsePatches = coarseMesh_.boundaryMesh();
const volScalarField::Boundary& Qrp = Qr_.boundaryField();
const volScalarField::Boundary& qrp = qr_.boundaryField();
label count = 0;
forAll(Qrp, patchi)
forAll(qrp, patchi)
{
//const polyPatch& pp = mesh_.boundaryMesh()[patchi];
const fvPatchScalarField& QrPatchi = Qrp[patchi];
const fvPatchScalarField& qrPatchi = qrp[patchi];
if ((isA<fixedValueFvPatchScalarField>(QrPatchi)))
if ((isA<fixedValueFvPatchScalarField>(qrPatchi)))
{
selectedPatches_[count] = QrPatchi.patch().index();
selectedPatches_[count] = qrPatchi.patch().index();
nLocalCoarseFaces_ += coarsePatches[patchi].size();
count++;
}
@ -264,11 +264,11 @@ Foam::radiation::viewFactor::viewFactor(const volScalarField& T)
mesh_,
finalAgglom_
),
Qr_
qr_
(
IOobject
(
"Qr",
"qr",
mesh_.time().timeName(),
mesh_,
IOobject::MUST_READ,
@ -322,11 +322,11 @@ Foam::radiation::viewFactor::viewFactor
mesh_,
finalAgglom_
),
Qr_
qr_
(
IOobject
(
"Qr",
"qr",
mesh_.time().timeName(),
mesh_,
IOobject::MUST_READ,
@ -394,7 +394,7 @@ void Foam::radiation::viewFactor::insertMatrixElements
void Foam::radiation::viewFactor::calculate()
{
// Store previous iteration
Qr_.storePrevIter();
qr_.storePrevIter();
scalarField compactCoarseT(map_->constructSize(), 0.0);
scalarField compactCoarseE(map_->constructSize(), 0.0);
@ -407,7 +407,7 @@ void Foam::radiation::viewFactor::calculate()
DynamicList<scalar> localCoarseEave(nLocalCoarseFaces_);
DynamicList<scalar> localCoarseHoave(nLocalCoarseFaces_);
volScalarField::Boundary& QrBf = Qr_.boundaryFieldRef();
volScalarField::Boundary& qrBf = qr_.boundaryFieldRef();
forAll(selectedPatches_, i)
{
@ -416,17 +416,17 @@ void Foam::radiation::viewFactor::calculate()
const scalarField& Tp = T_.boundaryField()[patchID];
const scalarField& sf = mesh_.magSf().boundaryField()[patchID];
fvPatchScalarField& QrPatch = QrBf[patchID];
fvPatchScalarField& qrPatch = qrBf[patchID];
greyDiffusiveViewFactorFixedValueFvPatchScalarField& Qrp =
greyDiffusiveViewFactorFixedValueFvPatchScalarField& qrp =
refCast
<
greyDiffusiveViewFactorFixedValueFvPatchScalarField
>(QrPatch);
>(qrPatch);
const scalarList eb = Qrp.emissivity();
const scalarList eb = qrp.emissivity();
const scalarList& Hoi = Qrp.Qro();
const scalarList& Hoi = qrp.qro();
const polyPatch& pp = coarseMesh_.boundaryMesh()[patchID];
const labelList& coarsePatchFace = coarseMesh_.patchFaceMap()[patchID];
@ -499,23 +499,23 @@ void Foam::radiation::viewFactor::calculate()
// Create global size vectors
scalarField T(totalNCoarseFaces_, 0.0);
scalarField E(totalNCoarseFaces_, 0.0);
scalarField QrExt(totalNCoarseFaces_, 0.0);
scalarField qrExt(totalNCoarseFaces_, 0.0);
// Fill lists from compact to global indexes.
forAll(compactCoarseT, i)
{
T[compactGlobalIds[i]] = compactCoarseT[i];
E[compactGlobalIds[i]] = compactCoarseE[i];
QrExt[compactGlobalIds[i]] = compactCoarseHo[i];
qrExt[compactGlobalIds[i]] = compactCoarseHo[i];
}
Pstream::listCombineGather(T, maxEqOp<scalar>());
Pstream::listCombineGather(E, maxEqOp<scalar>());
Pstream::listCombineGather(QrExt, maxEqOp<scalar>());
Pstream::listCombineGather(qrExt, maxEqOp<scalar>());
Pstream::listCombineScatter(T);
Pstream::listCombineScatter(E);
Pstream::listCombineScatter(QrExt);
Pstream::listCombineScatter(qrExt);
// Net radiation
scalarField q(totalNCoarseFaces_, 0.0);
@ -538,7 +538,7 @@ void Foam::radiation::viewFactor::calculate()
if (i==j)
{
C(i, j) = invEj - (invEj - 1.0)*Fmatrix_()(i, j);
q[i] += (Fmatrix_()(i, j) - 1.0)*sigmaT4 - QrExt[j];
q[i] += (Fmatrix_()(i, j) - 1.0)*sigmaT4 - qrExt[j];
}
else
{
@ -594,7 +594,7 @@ void Foam::radiation::viewFactor::calculate()
if (i==j)
{
q[i] += (Fmatrix_()(i, j) - 1.0)*sigmaT4 - QrExt[j];
q[i] += (Fmatrix_()(i, j) - 1.0)*sigmaT4 - qrExt[j];
}
else
{
@ -614,7 +614,7 @@ void Foam::radiation::viewFactor::calculate()
}
}
// Scatter q and fill Qr
// Scatter q and fill qr
Pstream::listCombineScatter(q);
Pstream::listCombineGather(q, maxEqOp<scalar>());
@ -625,7 +625,7 @@ void Foam::radiation::viewFactor::calculate()
const polyPatch& pp = mesh_.boundaryMesh()[patchID];
if (pp.size() > 0)
{
scalarField& Qrp = QrBf[patchID];
scalarField& qrp = qrBf[patchID];
const scalarField& sf = mesh_.magSf().boundaryField()[patchID];
const labelList& agglom = finalAgglom_[patchID];
label nAgglom = max(agglom)+1;
@ -646,8 +646,8 @@ void Foam::radiation::viewFactor::calculate()
{
label facei = fineFaces[k];
Qrp[facei] = q[globalCoarse];
heatFlux += Qrp[facei]*sf[facei];
qrp[facei] = q[globalCoarse];
heatFlux += qrp[facei]*sf[facei];
}
globCoarseId ++;
}
@ -656,11 +656,11 @@ void Foam::radiation::viewFactor::calculate()
if (debug)
{
forAll(QrBf, patchID)
forAll(qrBf, patchID)
{
const scalarField& Qrp = QrBf[patchID];
const scalarField& qrp = qrBf[patchID];
const scalarField& magSf = mesh_.magSf().boundaryField()[patchID];
scalar heatFlux = gSum(Qrp*magSf);
scalar heatFlux = gSum(qrp*magSf);
InfoInFunction
<< "Total heat transfer rate at patch: "
@ -669,8 +669,8 @@ void Foam::radiation::viewFactor::calculate()
}
}
// Relax Qr if necessary
Qr_.relax();
// Relax qr if necessary
qr_.relax();
}

6
src/thermophysicalModels/radiation/radiationModels/viewFactor/viewFactor.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -80,7 +80,7 @@ class viewFactor
singleCellFvMesh coarseMesh_;
//- Net radiative heat flux [W/m2]
volScalarField Qr_;
volScalarField qr_;
//- View factor matrix
autoPtr<scalarSquareMatrix> Fmatrix_;
@ -168,7 +168,7 @@ public:
// Access
//- Const access to total radiative heat flux field
inline const volScalarField& Qr() const;
inline const volScalarField& qr() const;
};

7
src/thermophysicalModels/radiation/radiationModels/viewFactor/viewFactorI.H
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -23,11 +23,10 @@ License
\*---------------------------------------------------------------------------*/
inline const Foam::volScalarField& Foam::radiation::viewFactor::Qr() const
inline const Foam::volScalarField& Foam::radiation::viewFactor::qr() const
{
return Qr_;
return qr_;
}
// ************************************************************************* //

4
tutorials/combustion/fireFoam/LES/flameSpreadWaterSuppressionPanel/0/T
View File

@ -53,8 +53,8 @@ boundaryField
filmRegion filmRegion;
Tnbr T;
kappaMethod fluidThermo;
QrNbr none;
Qr Qr;
qrNbr none;
qr qr;
filmDeltaDry 0.0;
filmDeltaWet 2e-4;
value $internalField;

4
tutorials/combustion/fireFoam/LES/flameSpreadWaterSuppressionPanel/0/pyrolysisRegion/T
View File

@ -38,8 +38,8 @@ boundaryField
filmRegion filmRegion;
Tnbr T;
kappaMethod solidThermo;
QrNbr Qr;
Qr none;
qrNbr qr;
qr none;
filmDeltaDry 0.0;
filmDeltaWet 2e-4;
value uniform 298;

4
tutorials/combustion/fireFoam/LES/flameSpreadWaterSuppressionPanel/0/pyrolysisRegion/Qr → tutorials/combustion/fireFoam/LES/flameSpreadWaterSuppressionPanel/0/pyrolysisRegion/qr
View File

@ -11,7 +11,7 @@ FoamFile
format ascii;
class volScalarField;
location "0";
object Qr;
object qr;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -38,7 +38,7 @@ boundaryField
sampleMode nearestPatchFace;
samplePatch region0_to_pyrolysisRegion_coupledWall;
offset (0 0 0);
field Qr; // this is the name of Qr field in region0
field qr; // this is the name of qr field in region0
setAverage no;
average 0;
value uniform 0;

4
tutorials/combustion/fireFoam/LES/flameSpreadWaterSuppressionPanel/constant/pyrolysisZones
View File

@ -28,9 +28,9 @@ pyrolysis
filmCoupled true;
gasHSource true;
QrHSource false;
qrHSource false;
radFluxName Qr;
radFluxName qr;
moveMesh false; // true;
minimumDelta 1e-6;

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

@ -51,8 +51,8 @@ boundaryField
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
kappaMethod fluidThermo;
QrNbr none;
Qr Qr;
qrNbr none;
qr qr;
value $internalField;
}
}

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

@ -42,8 +42,8 @@ boundaryField
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
kappaMethod solidThermo;
QrNbr Qr;
Qr none;
qrNbr qr;
qr none;
value uniform 298.15;
}

6
tutorials/combustion/fireFoam/LES/oppositeBurningPanels/0/panelRegion/Qr → tutorials/combustion/fireFoam/LES/oppositeBurningPanels/0/panelRegion/qr
View File

@ -11,7 +11,7 @@ FoamFile
format ascii;
class volScalarField;
location "0";
object Qr;
object qr;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -43,7 +43,7 @@ boundaryField
sampleMode nearestPatchFace;
samplePatch region0_to_panelRegion_fLeft_zone;
offset (0 0 0);
field Qr;
field qr;
setAverage no;
average 0;
value uniform 0;
@ -56,7 +56,7 @@ boundaryField
sampleMode nearestPatchFace;
samplePatch region0_to_panelRegion_fRight_zone;
offset (0 0 0);
field Qr;
field qr;
setAverage no;
average 0;
value uniform 0;

2
tutorials/combustion/fireFoam/LES/oppositeBurningPanels/constant/pyrolysisZones
View File

@ -26,7 +26,7 @@ pyrolysis
reactingOneDimCoeffs
{
gasHSource false; //Energy source term due to pyrolysis gas
QrHSource false; //Energy source term due in depht radiation
qrHSource false; //Energy source term due in depht radiation
minimumDelta 1e-12;
reactionDeltaMin 1e-12;

2
tutorials/combustion/fireFoam/LES/smallPoolFire2D/constant/pyrolysisZones
View File

@ -25,7 +25,7 @@ pyrolysis
reactingOneDimCoeffs
{
radFluxName Qr;
radFluxName qr;
minimumDelta 1e-8;

2
tutorials/combustion/fireFoam/LES/smallPoolFire3D/constant/pyrolysisZones
View File

@ -25,7 +25,7 @@ pyrolysis
reactingOneDimCoeffs
{
radFluxName Qr;
radFluxName qr;
minimumDelta 1e-8;

2
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/0/Qr → tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/0/qr
View File

@ -10,7 +10,7 @@ FoamFile
version 2.0;
format ascii;
class volScalarField;
object Qr;
object qr;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

2
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/Allrun.pre
View File

@ -12,7 +12,7 @@ runApplication splitMeshRegions -cellZones -overwrite
# remove fluid fields from solid regions (important for post-processing)
for i in heater leftSolid rightSolid
do
rm -f 0*/$i/{rho,nut,alphat,epsilon,k,U,p_rgh,Qr,G,IDefault}
rm -f 0*/$i/{rho,nut,alphat,epsilon,k,U,p_rgh,qr,G,IDefault}
done
for i in bottomAir topAir heater leftSolid rightSolid

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

@ -54,8 +54,8 @@ T
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
kappaMethod fluidThermo;
QrNbr none;
Qr Qr;
qrNbr none;
qr qr;
value uniform 300;
}
}
@ -113,7 +113,7 @@ p
}
}
Qr
qr
{
internalField uniform 0;
boundaryField
@ -122,7 +122,7 @@ Qr
{
type greyDiffusiveRadiationViewFactor;
emissivityMode lookup;
Qro uniform 0;
qro uniform 0;
emissivity uniform 1.0;
value uniform 0;
}
@ -131,7 +131,7 @@ Qr
{
type greyDiffusiveRadiationViewFactor;
emissivityMode solidRadiation;
Qro uniform 0;
qro uniform 0;
value uniform 0;
}
}

2
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/system/bottomAir/fvSolution
View File

@ -74,7 +74,7 @@ relaxationFactors
"(k|epsilon|omega)" 0.7;
G 0.7;
"ILambda.*" 0.7;
Qr 0.7;
qr 0.7;
}
}

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

@ -47,8 +47,8 @@ T
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
kappaMethod solidThermo;
QrNbr Qr;
Qr none;
qrNbr qr;
qr none;
value uniform 300;
}
@ -57,8 +57,8 @@ T
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
kappaMethod solidThermo;
QrNbr Qr;
Qr none;
qrNbr qr;
qr none;
value uniform 300;
}

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

@ -43,8 +43,8 @@ T
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
kappaMethod solidThermo;
QrNbr Qr;
Qr none;
qrNbr qr;
qr none;
value uniform 300;
}
@ -53,8 +53,8 @@ T
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
kappaMethod solidThermo;
QrNbr Qr;
Qr none;
qrNbr qr;
qr none;
value uniform 300;
}

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

@ -43,8 +43,8 @@ T
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
kappaMethod solidThermo;
QrNbr Qr;
Qr none;
qrNbr qr;
qr none;
value uniform 300;
}
@ -53,8 +53,8 @@ T
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
kappaMethod solidThermo;
QrNbr Qr;
Qr none;
qrNbr qr;
qr none;
value uniform 300;
}

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

@ -65,8 +65,8 @@ T
type compressible::turbulentTemperatureRadCoupledMixed;
Tnbr T;
kappaMethod fluidThermo;
QrNbr none;
Qr Qr;
qrNbr none;
qr qr;
value uniform 300;
}
}
@ -160,7 +160,7 @@ p
}
}
Qr
qr
{
internalField uniform 0;
boundaryField
@ -169,7 +169,7 @@ Qr
{
type greyDiffusiveRadiationViewFactor;
emissivityMode lookup;
Qro uniform 0;
qro uniform 0;
emissivity uniform 1.0;
value uniform 0;
}
@ -178,7 +178,7 @@ Qr
{
type greyDiffusiveRadiationViewFactor;
emissivityMode solidRadiation;
Qro uniform 0;
qro uniform 0;
value uniform 0;
}
}

2
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/system/topAir/fvSolution
View File

@ -74,7 +74,7 @@ relaxationFactors
"(k|epsilon|omega)" 0.7;
G 0.7;
"ILambda.*" 0.7;
Qr 0.7;
qr 0.7;
}
}