diff --git a/src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.C b/src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.C
index 0c96555aa582c5f4f6ddf6feb6bd53838f337275..36e562882c5226193434f5fb5aefbd0f84814532 100644
--- a/src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.C
+++ b/src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.C
@@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
- Copyright (C) 2015-2020 OpenCFD Ltd.
+ Copyright (C) 2015-2021 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@@ -120,6 +120,7 @@ Foam::humidityTemperatureCoupledMixedFvPatchScalarField::thicknessField
}
+
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
Foam::humidityTemperatureCoupledMixedFvPatchScalarField::
@@ -159,7 +160,9 @@ humidityTemperatureCoupledMixedFvPatchScalarField
fluid_(false),
cp_(patch().size(), Zero),
thickness_(patch().size(), Zero),
- rho_(patch().size(), Zero)
+ rho_(patch().size(), Zero),
+ thicknessLayers_(0),
+ kappaLayers_(0)
{
this->refValue() = 0.0;
this->refGrad() = 0.0;
@@ -199,7 +202,9 @@ humidityTemperatureCoupledMixedFvPatchScalarField
fluid_(psf.fluid_),
cp_(psf.cp_, mapper),
thickness_(psf.thickness_, mapper),
- rho_(psf.rho_, mapper)
+ rho_(psf.rho_, mapper),
+ thicknessLayers_(psf.thicknessLayers_),
+ kappaLayers_(psf.kappaLayers_)
{}
@@ -234,7 +239,9 @@ humidityTemperatureCoupledMixedFvPatchScalarField
fluid_(false),
cp_(patch().size(), Zero),
thickness_(patch().size(), Zero),
- rho_(patch().size(), Zero)
+ rho_(patch().size(), Zero),
+ thicknessLayers_(0),
+ kappaLayers_(0)
{
if (!isA<mappedPatchBase>(this->patch().patch()))
{
@@ -254,6 +261,11 @@ humidityTemperatureCoupledMixedFvPatchScalarField
fluid_ = true;
}
+ if (dict.readIfPresent("thicknessLayers", thicknessLayers_))
+ {
+ dict.readEntry("kappaLayers", kappaLayers_);
+ }
+
if (fluid_)
{
switch(mode_)
@@ -357,7 +369,9 @@ humidityTemperatureCoupledMixedFvPatchScalarField
fluid_(psf.fluid_),
cp_(psf.cp_),
thickness_(psf.thickness_),
- rho_(psf.rho_)
+ rho_(psf.rho_),
+ thicknessLayers_(psf.thicknessLayers_),
+ kappaLayers_(psf.kappaLayers_)
{}
@@ -468,6 +482,16 @@ void Foam::humidityTemperatureCoupledMixedFvPatchScalarField::updateCoeffs()
myKDelta_ = K*patch().deltaCoeffs();
+ if (thicknessLayers_.size() > 0)
+ {
+ myKDelta_ = 1.0/myKDelta_;
+ forAll(thicknessLayers_, iLayer)
+ {
+ myKDelta_ += thicknessLayers_[iLayer]/kappaLayers_[iLayer];
+ }
+ myKDelta_ = 1.0/myKDelta_;
+ }
+
scalarField dm(patch().size(), Zero);
// Fluid Side
@@ -484,6 +508,8 @@ void Foam::humidityTemperatureCoupledMixedFvPatchScalarField::updateCoeffs()
scalarField hfg(patch().size(), Zero);
scalarField htc(patch().size(), GREAT);
scalarField liquidRho(patch().size(), Zero);
+ scalarField Tdew(patch().size(), Zero);
+ scalarField RH(patch().size(), Zero);
fixedGradientFvPatchField<scalar>& Yp =
const_cast<fixedGradientFvPatchField<scalar>&>
@@ -537,25 +563,26 @@ void Foam::humidityTemperatureCoupledMixedFvPatchScalarField::updateCoeffs()
const scalar invMwmean =
Yi[faceI]/Mv + (1.0 - Yi[faceI])/Mcomp_;
const scalar Xv = Yi[faceI]/invMwmean/Mv;
- const scalar RH = min(Xv*pf/pSat, 1.0);
+ RH[faceI] = min(Xv*pf/pSat, 1.0);
scalar RHmin = 0.01;
- scalar Tdew = -GREAT;
+ Tdew[faceI] = -GREAT;
- if (RH > RHmin)
+ if (RH[faceI] > RHmin)
{
scalar b = 243.5;
scalar c = 17.65;
scalar TintDeg = Tint - 273;
- Tdew =
- b*(log(RH) + (c*TintDeg)/(b + TintDeg))
- /(c - log(RH) - ((c*TintDeg)/(b + TintDeg))) + 273;
+ Tdew[faceI] =
+ b*(log(RH[faceI]) + (c*TintDeg)/(b + TintDeg))
+ /(c - log(RH[faceI]) - ((c*TintDeg)/(b + TintDeg)))
+ + 273;
}
if
(
- Tf < Tdew
- && RH > RHmin
+ Tf < Tdew[faceI]
+ && RH[faceI] > RHmin
&& (
mode_ == mtCondensation
|| mode_ == mtCondensationAndEvaporation
@@ -613,7 +640,7 @@ void Foam::humidityTemperatureCoupledMixedFvPatchScalarField::updateCoeffs()
htc[faceI] = htcCondensation(TSat, Re)*nbrK[faceI]/L_;
}
- else if (Tf > Tdew && Tf < Tvap_ && mass_[faceI] > 0.0)
+ else if (Tf > Tdew[faceI] && Tf < Tvap_ && mass_[faceI] > 0.0)
{
htc[faceI] = htcCondensation(TSat, Re)*nbrK[faceI]/L_;
}
@@ -649,6 +676,15 @@ void Foam::humidityTemperatureCoupledMixedFvPatchScalarField::updateCoeffs()
// Heat flux due to change of phase [W/m2]
dmHfg_ = dm*hfg;
+
+ // Output RH and Tdew
+ scalarField& bRH =
+ thicknessField("RH", refCast<const fvMesh>(mesh));
+ bRH = RH;
+
+ scalarField& bTdew =
+ thicknessField("Tdew", refCast<const fvMesh>(mesh));
+ bTdew = Tdew;
}
else
{
@@ -763,6 +799,12 @@ void Foam::humidityTemperatureCoupledMixedFvPatchScalarField::write
liquidDict_.write(os);
}
+ if (thicknessLayers_.size())
+ {
+ thicknessLayers_.writeEntry("thicknessLayers", os);
+ kappaLayers_.writeEntry("kappaLayers", os);
+ }
+
temperatureCoupledBase::write(os);
}
diff --git a/src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.H b/src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.H
index 2f6d5801d9cf3549eda0ba168aa1174080139008..bc5d081cdb6e943be8e64ca117dd54ad32524a60 100644
--- a/src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.H
+++ b/src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.H
@@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
- Copyright (C) 2015-2016 OpenCFD Ltd.
+ Copyright (C) 2015-2021 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@@ -91,6 +91,8 @@ Description
- The boundary condition of this specie on the coupled wall must be
fixedGradient in order to allow condensation or evaporation of the
vapour in or out of this wall
+ - Addition of extra layers in possible using thicknessLayers and
+ kappaLayers
Example usage:
@@ -127,6 +129,9 @@ Description
}
}
+ thicknessLayers (0.1 0.2 0.3 0.4);
+ kappaLayers (1 2 3 4);
+
// thickness, density and cp required for inert and condensation
// modes
@@ -267,6 +272,12 @@ private:
//- Density field for inert mode
scalarField rho_;
+ //- Thickness of layers
+ scalarList thicknessLayers_;
+
+ //- Conductivity of layers
+ scalarList kappaLayers_;
+
// Private members