diff --git a/src/turbulenceModels/RAS/compressible/RASModel/RASModel.C b/src/turbulenceModels/RAS/compressible/RASModel/RASModel.C
index 6a4a6c055b62170f983195e5810eb3cc96efbb04..61dbc1fc277ea9b2c9d65d01d0b9cbb0e34e7b15 100644
--- a/src/turbulenceModels/RAS/compressible/RASModel/RASModel.C
+++ b/src/turbulenceModels/RAS/compressible/RASModel/RASModel.C
@@ -89,14 +89,24 @@ RASModel::RASModel
kappa_
(
- subDict("wallFunctionCoeffs").lookupOrAddDefault<scalar>
+ dimensioned<scalar>::lookupOrAddToDict
(
"kappa",
+ subDict("wallFunctionCoeffs"),
0.4187
)
),
- E_(subDict("wallFunctionCoeffs").lookupOrAddDefault<scalar>("E", 9.0)),
- yPlusLam_(yPlusLam(kappa_, E_)),
+ E_
+ (
+ dimensioned<scalar>::lookupOrAddToDict
+ (
+ "E",
+ subDict("wallFunctionCoeffs"),
+ 9.0
+ )
+ ),
+
+ yPlusLam_(yPlusLam(kappa_.value(), E_.value())),
k0_("k0", dimVelocity*dimVelocity, SMALL),
epsilon0_("epsilon", k0_.dimensions()/dimTime, SMALL),
@@ -171,10 +181,10 @@ bool RASModel::read()
lookup("turbulence") >> turbulence_;
coeffDict_ = subDict(type() + "Coeffs");
- subDict("wallFunctionCoeffs").readIfPresent<scalar>("kappa", kappa_);
- subDict("wallFunctionCoeffs").readIfPresent<scalar>("E", E_);
+ kappa_.readIfPresent(subDict("wallFunctionCoeffs"));
+ E_.readIfPresent(subDict("wallFunctionCoeffs"));
- yPlusLam_ = yPlusLam(kappa_, E_);
+ yPlusLam_ = yPlusLam(kappa_.value(), E_.value());
readIfPresent("k0", k0_);
readIfPresent("epsilon0", epsilon0_);
diff --git a/src/turbulenceModels/RAS/compressible/RASModel/RASModel.H b/src/turbulenceModels/RAS/compressible/RASModel/RASModel.H
index 06e3b6b8fb788f1f7859fd1cf7e18f9f2e9ba579..39960659a4a5bb2900450387328022c97d3ee2af 100644
--- a/src/turbulenceModels/RAS/compressible/RASModel/RASModel.H
+++ b/src/turbulenceModels/RAS/compressible/RASModel/RASModel.H
@@ -91,8 +91,8 @@ protected:
Switch printCoeffs_;
dictionary coeffDict_;
- scalar kappa_;
- scalar E_;
+ dimensionedScalar kappa_;
+ dimensionedScalar E_;
scalar yPlusLam(const scalar kappa, const scalar E);
scalar yPlusLam_;
@@ -226,13 +226,13 @@ public:
//- Return kappa for use in wall-functions
- scalar kappa() const
+ dimensionedScalar kappa() const
{
return kappa_;
}
//- Return E for use in wall-functions
- scalar E() const
+ dimensionedScalar E() const
{
return E_;
}
diff --git a/src/turbulenceModels/RAS/compressible/kOmegaSST/kOmegaWallFunctionsI.H b/src/turbulenceModels/RAS/compressible/kOmegaSST/kOmegaWallFunctionsI.H
index 1f34c5be9f6bfaa25e2e05257fce4db7615d923d..e7b9db40cd1b18a21e17e37de7123e2601258e5d 100644
--- a/src/turbulenceModels/RAS/compressible/kOmegaSST/kOmegaWallFunctionsI.H
+++ b/src/turbulenceModels/RAS/compressible/kOmegaSST/kOmegaWallFunctionsI.H
@@ -89,7 +89,7 @@ Description
omega_[faceCelli] +=
sqrt(k_[faceCelli])
- /(Cmu25*kappa_*y_[faceCelli]);
+ /(Cmu25*kappa_.value()*y_[faceCelli]);
if (yPlus > yPlusLam_)
{
@@ -97,7 +97,7 @@ Description
(mutw[facei] + muw[facei])
*magFaceGradU[facei]
*Cmu25*sqrt(k_[faceCelli])
- /(kappa_*y_[faceCelli]);
+ /(kappa_.value()*y_[faceCelli]);
}
}
}
diff --git a/src/turbulenceModels/RAS/compressible/kOmegaSST/kOmegaWallViscosityI.H b/src/turbulenceModels/RAS/compressible/kOmegaSST/kOmegaWallViscosityI.H
index 4aae893e977d889351880a16be4c28e26ede42a4..68ce16201cf3489b464b8e52bcaf681a7c5aa8f4 100644
--- a/src/turbulenceModels/RAS/compressible/kOmegaSST/kOmegaWallViscosityI.H
+++ b/src/turbulenceModels/RAS/compressible/kOmegaSST/kOmegaWallViscosityI.H
@@ -58,7 +58,7 @@ Description
{
mutw[facei] =
muw[facei]
- *(yPlus*kappa_/log(E_*yPlus) - 1);
+ *(yPlus*kappa_.value()/log(E_.value()*yPlus) - 1);
}
else
{
diff --git a/src/turbulenceModels/RAS/compressible/wallFunctions/mutWallFunctions/mutStandardRoughWallFunction/mutStandardRoughWallFunctionFvPatchScalarField.C b/src/turbulenceModels/RAS/compressible/wallFunctions/mutWallFunctions/mutStandardRoughWallFunction/mutStandardRoughWallFunctionFvPatchScalarField.C
index 5238bcda209369193ea7d04a48c006026dee49d6..b019bfdd516b613df46f722ad7951e2cb985e9bd 100644
--- a/src/turbulenceModels/RAS/compressible/wallFunctions/mutWallFunctions/mutStandardRoughWallFunction/mutStandardRoughWallFunctionFvPatchScalarField.C
+++ b/src/turbulenceModels/RAS/compressible/wallFunctions/mutWallFunctions/mutStandardRoughWallFunction/mutStandardRoughWallFunctionFvPatchScalarField.C
@@ -123,14 +123,14 @@ void mutStandardRoughWallFunctionFvPatchScalarField::evaluate
const RASModel& rasModel
= db().lookupObject<RASModel>("RASProperties");
- const scalar kappa = rasModel.kappa();
- const scalar E = rasModel.E();
+ const scalar kappa = rasModel.kappa().value();
+ const scalar E = rasModel.E().value();
const scalar yPlusLam = 11.225;
// The reciprical of the distance to the adjacent cell centre.
const scalarField& ry = patch().deltaCoeffs();
- const fvPatchVectorField& U =
+ const fvPatchVectorField& U =
patch().lookupPatchField<volVectorField, vector>("U");
const fvPatchScalarField& rho =
@@ -139,7 +139,7 @@ void mutStandardRoughWallFunctionFvPatchScalarField::evaluate
// The flow velocity at the adjacent cell centre.
scalarField magUp = mag(U.patchInternalField() - U);
- const scalarField& muw =
+ const scalarField& muw =
patch().lookupPatchField<volScalarField, scalar>("mu");
scalarField& mutw = *this;
@@ -155,7 +155,7 @@ void mutStandardRoughWallFunctionFvPatchScalarField::evaluate
//if (KsPlusBasedOnYPlus_)
{
- // If KsPlus is based on YPlus the extra term added to the law
+ // If KsPlus is based on YPlus the extra term added to the law
// of the wall will depend on yPlus.
forAll(mutw, facei)
{
@@ -206,7 +206,7 @@ void mutStandardRoughWallFunctionFvPatchScalarField::evaluate
const scalar sint_2 = sin(t_2);
const scalar logt_1 = log(t_1);
G = logt_1*sint_2;
- yPlusGPrime =
+ yPlusGPrime =
(c_1*sint_2*KsPlus/t_1) + (c_3*logt_1*cos(t_2));
}
@@ -268,7 +268,7 @@ void mutStandardRoughWallFunctionFvPatchScalarField::evaluate
mag(ryPlusLam*(yPlus - yPlusLast)) > 0.0001
&& ++iter < 10
);
-
+
if (yPlus > yPlusLam)
{
mutw[facei] = muw[facei]*(yPlus*yPlus/Re - 1);
@@ -291,7 +291,7 @@ void mutStandardRoughWallFunctionFvPatchScalarField::write(Ostream& os) const
<< roughnessConstant_ << token::END_STATEMENT << nl;
os.writeKeyword("roughnessFudgeFactor")
<< roughnessFudgeFactor_ << token::END_STATEMENT << nl;
-}
+}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/src/turbulenceModels/RAS/compressible/wallFunctions/wallFunctionsI.H b/src/turbulenceModels/RAS/compressible/wallFunctions/wallFunctionsI.H
index 19ae89d4ca43042671e0cf80eaf707849294c89f..35b1a401cfb4353f6272cf60d0a54937d6838961 100644
--- a/src/turbulenceModels/RAS/compressible/wallFunctions/wallFunctionsI.H
+++ b/src/turbulenceModels/RAS/compressible/wallFunctions/wallFunctionsI.H
@@ -90,7 +90,7 @@ Description
epsilon_[faceCelli] +=
Cmu75*pow(k_[faceCelli], 1.5)
- /(kappa_*RASModel::y_[patchi][facei]);
+ /(kappa_.value()*RASModel::y_[patchi][facei]);
if (yPlus > yPlusLam_)
{
@@ -98,7 +98,7 @@ Description
(mutw[facei] + muw[facei])
*magFaceGradU[facei]
*Cmu25*sqrt(k_[faceCelli])
- /(kappa_*RASModel::y_[patchi][facei]);
+ /(kappa_.value()*RASModel::y_[patchi][facei]);
}
}
}
diff --git a/src/turbulenceModels/RAS/compressible/wallFunctions/wallViscosityI.H b/src/turbulenceModels/RAS/compressible/wallFunctions/wallViscosityI.H
index 938555e860675d50fc4c45758cdf1b49475737c9..659e3286c2c3f1f99a8349d56bf96a945e3aa69f 100644
--- a/src/turbulenceModels/RAS/compressible/wallFunctions/wallViscosityI.H
+++ b/src/turbulenceModels/RAS/compressible/wallFunctions/wallViscosityI.H
@@ -58,7 +58,7 @@ Description
{
mutw[facei] =
muw[facei]
- *(yPlus*kappa_/log(E_*yPlus) - 1);
+ *(yPlus*kappa_.value()/log(E_.value()*yPlus) - 1);
}
else
{
diff --git a/src/turbulenceModels/RAS/incompressible/LienCubicKELowRe/LienCubicKELowReSetWallDissipation.H b/src/turbulenceModels/RAS/incompressible/LienCubicKELowRe/LienCubicKELowReSetWallDissipation.H
index 76d0a2021171f4d1f4589c57640bcf4d7f616141..ef216dcf68765654ef1d86c72b077c71ff542f57 100644
--- a/src/turbulenceModels/RAS/incompressible/LienCubicKELowRe/LienCubicKELowReSetWallDissipation.H
+++ b/src/turbulenceModels/RAS/incompressible/LienCubicKELowRe/LienCubicKELowReSetWallDissipation.H
@@ -43,7 +43,7 @@
epsilon_[faceCelli] +=
Cmu75*pow(k_[faceCelli], 1.5)
/(
- kappa_*y_[faceCelli]
+ kappa_.value()*y_[faceCelli]
*(1.0 - exp(-Aepsilon_.value()*yStar_[faceCelli]))
)
*exp(-Amu_.value()*sqr(yStar_[faceCelli]));
diff --git a/src/turbulenceModels/RAS/incompressible/LienLeschzinerLowRe/LienLeschzinerLowReSetWallDissipation.H b/src/turbulenceModels/RAS/incompressible/LienLeschzinerLowRe/LienLeschzinerLowReSetWallDissipation.H
index cfda062b594bf64966001211a22cfc0b01ea98d3..a998d3af8e092ecb639f6fac1fae615ad0bd828a 100644
--- a/src/turbulenceModels/RAS/incompressible/LienLeschzinerLowRe/LienLeschzinerLowReSetWallDissipation.H
+++ b/src/turbulenceModels/RAS/incompressible/LienLeschzinerLowRe/LienLeschzinerLowReSetWallDissipation.H
@@ -43,7 +43,7 @@
epsilon_[faceCelli] +=
Cmu75*pow(k_[faceCelli], 1.5)
/(
- kappa_*y_[faceCelli]
+ kappa_.value()*y_[faceCelli]
*(1.0 - exp(-Aepsilon*yStar_[faceCelli]))
)
*exp(-Amu*sqr(yStar_[faceCelli]));
diff --git a/src/turbulenceModels/RAS/incompressible/RASModel/RASModel.C b/src/turbulenceModels/RAS/incompressible/RASModel/RASModel.C
index 36567e7b8c5778d66646f4644b637ff5cabc16c4..34be1e6e39165d0235d699c5470d539b793f18d2 100644
--- a/src/turbulenceModels/RAS/incompressible/RASModel/RASModel.C
+++ b/src/turbulenceModels/RAS/incompressible/RASModel/RASModel.C
@@ -85,14 +85,24 @@ RASModel::RASModel
kappa_
(
- subDict("wallFunctionCoeffs").lookupOrAddDefault<scalar>
+ dimensioned<scalar>::lookupOrAddToDict
(
"kappa",
+ subDict("wallFunctionCoeffs"),
0.4187
)
),
- E_(subDict("wallFunctionCoeffs").lookupOrAddDefault<scalar>("E", 9.0)),
- yPlusLam_(yPlusLam(kappa_, E_)),
+ E_
+ (
+ dimensioned<scalar>::lookupOrAddToDict
+ (
+ "E",
+ subDict("wallFunctionCoeffs"),
+ 9.0
+ )
+ ),
+
+ yPlusLam_(yPlusLam(kappa_.value(), E_.value())),
k0_("k0", dimVelocity*dimVelocity, SMALL),
epsilon0_("epsilon", k0_.dimensions()/dimTime, SMALL),
@@ -168,14 +178,14 @@ bool RASModel::read()
lookup("turbulence") >> turbulence_;
coeffDict_ = subDict(type() + "Coeffs");
- subDict("wallFunctionCoeffs").readIfPresent<scalar>("kappa", kappa_);
- subDict("wallFunctionCoeffs").readIfPresent<scalar>("E", E_);
+ kappa_.readIfPresent(subDict("wallFunctionCoeffs"));
+ E_.readIfPresent(subDict("wallFunctionCoeffs"));
- yPlusLam_ = yPlusLam(kappa_, E_);
+ yPlusLam_ = yPlusLam(kappa_.value(), E_.value());
- readIfPresent("k0", k0_);
- readIfPresent("epsilon0", epsilon0_);
- readIfPresent("epsilonSmall", epsilonSmall_);
+ k0_.readIfPresent(*this);
+ epsilon0_.readIfPresent(*this);
+ epsilonSmall_.readIfPresent(*this);
return true;
}
diff --git a/src/turbulenceModels/RAS/incompressible/RASModel/RASModel.H b/src/turbulenceModels/RAS/incompressible/RASModel/RASModel.H
index 77ed3a35a4d4be5bf84fef2e3b7026f3bcf1a77f..40e087612e3640320e627151f9d5bf1df9cf2c9a 100644
--- a/src/turbulenceModels/RAS/incompressible/RASModel/RASModel.H
+++ b/src/turbulenceModels/RAS/incompressible/RASModel/RASModel.H
@@ -88,8 +88,8 @@ protected:
Switch printCoeffs_;
dictionary coeffDict_;
- scalar kappa_;
- scalar E_;
+ dimensionedScalar kappa_;
+ dimensionedScalar E_;
scalar yPlusLam(const scalar kappa, const scalar E);
scalar yPlusLam_;
@@ -213,13 +213,13 @@ public:
//- Return kappa for use in wall-functions
- scalar kappa() const
+ dimensionedScalar kappa() const
{
return kappa_;
}
//- Return E for use in wall-functions
- scalar E() const
+ dimensionedScalar E() const
{
return E_;
}
diff --git a/src/turbulenceModels/RAS/incompressible/kOmegaSST/kOmegaWallFunctionsI.H b/src/turbulenceModels/RAS/incompressible/kOmegaSST/kOmegaWallFunctionsI.H
index 3b3a8a59bff1a711d07f6c23c9cb4bbc501d88ff..ab2bcc148e6f34dd211cb428594cbdbfdb92c0ab 100644
--- a/src/turbulenceModels/RAS/incompressible/kOmegaSST/kOmegaWallFunctionsI.H
+++ b/src/turbulenceModels/RAS/incompressible/kOmegaSST/kOmegaWallFunctionsI.H
@@ -87,7 +87,7 @@ Description
omega_[faceCelli] +=
sqrt(k_[faceCelli])
- /(Cmu25*kappa_*y_[faceCelli]);
+ /(Cmu25*kappa_.value()*y_[faceCelli]);
if (yPlus > yPlusLam_)
{
@@ -95,7 +95,7 @@ Description
(nutw[facei] + nuw[facei])
*magFaceGradU[facei]
*Cmu25*sqrt(k_[faceCelli])
- /(kappa_*y_[faceCelli]);
+ /(kappa_.value()*y_[faceCelli]);
}
}
}
diff --git a/src/turbulenceModels/RAS/incompressible/kOmegaSST/kOmegaWallViscosityI.H b/src/turbulenceModels/RAS/incompressible/kOmegaSST/kOmegaWallViscosityI.H
index 3231565d7a38a3cca60601ccad1a11ce12d1fdba..f2667bc675705215b39ab967ec5aee0297d4190c 100644
--- a/src/turbulenceModels/RAS/incompressible/kOmegaSST/kOmegaWallViscosityI.H
+++ b/src/turbulenceModels/RAS/incompressible/kOmegaSST/kOmegaWallViscosityI.H
@@ -55,7 +55,7 @@ Description
{
nutw[facei] =
nuw[facei]
- *(yPlus*kappa_/log(E_*yPlus) - 1);
+ *(yPlus*kappa_.value()/log(E_.value()*yPlus) - 1);
}
else
{
diff --git a/src/turbulenceModels/RAS/incompressible/wallFunctions/nonLinearWallFunctionsI.H b/src/turbulenceModels/RAS/incompressible/wallFunctions/nonLinearWallFunctionsI.H
index 07967c87d95bce1fa283e825d7cfd3776b1111c3..f286ad0b67031f29438e8ca3b2b5223100756bc4 100644
--- a/src/turbulenceModels/RAS/incompressible/wallFunctions/nonLinearWallFunctionsI.H
+++ b/src/turbulenceModels/RAS/incompressible/wallFunctions/nonLinearWallFunctionsI.H
@@ -89,7 +89,7 @@ Description
epsilon_[faceCelli] +=
Cmu75*pow(k_[faceCelli], 1.5)
- /(kappa_*y_[patchi][facei]);
+ /(kappa_.value()*y_[patchi][facei]);
if (yPlus > yPlusLam_)
{
@@ -97,7 +97,7 @@ Description
(nutw[facei] + nuw[facei])
*magFaceGradU[facei]
*Cmu25*sqrt(k_[faceCelli])
- /(kappa_*y_[patchi][facei])
+ /(kappa_.value()*y_[patchi][facei])
- (nonlinearStress_[faceCelli] && gradU_[faceCelli]);
}
}
diff --git a/src/turbulenceModels/RAS/incompressible/wallFunctions/nutWallFunctions/nutStandardRoughWallFunction/nutStandardRoughWallFunctionFvPatchScalarField.C b/src/turbulenceModels/RAS/incompressible/wallFunctions/nutWallFunctions/nutStandardRoughWallFunction/nutStandardRoughWallFunctionFvPatchScalarField.C
index 2f242c1a4feccc66b2f6308ab3828726c5300efb..dcc849510eaf88320a7c7ff1ba8fecde33942028 100644
--- a/src/turbulenceModels/RAS/incompressible/wallFunctions/nutWallFunctions/nutStandardRoughWallFunction/nutStandardRoughWallFunctionFvPatchScalarField.C
+++ b/src/turbulenceModels/RAS/incompressible/wallFunctions/nutWallFunctions/nutStandardRoughWallFunction/nutStandardRoughWallFunctionFvPatchScalarField.C
@@ -123,8 +123,8 @@ void nutStandardRoughWallFunctionFvPatchScalarField::evaluate
const RASModel& rasModel
= db().lookupObject<RASModel>("RASProperties");
- const scalar kappa = rasModel.kappa();
- const scalar E = rasModel.E();
+ const scalar kappa = rasModel.kappa().value();
+ const scalar E = rasModel.E().value();
const scalar yPlusLam = 11.225;
// The reciprical of the distance to the adjacent cell centre.
diff --git a/src/turbulenceModels/RAS/incompressible/wallFunctions/nutWallFunctions/nutStandardWallFunction/nutStandardWallFunctionFvPatchScalarField.C b/src/turbulenceModels/RAS/incompressible/wallFunctions/nutWallFunctions/nutStandardWallFunction/nutStandardWallFunctionFvPatchScalarField.C
index 3eee19f990f8a3a24f9902158d99882d97fe1d0b..f0d6eb714b7a274d50e8668b783ebfc4df762ab7 100644
--- a/src/turbulenceModels/RAS/incompressible/wallFunctions/nutWallFunctions/nutStandardWallFunction/nutStandardWallFunctionFvPatchScalarField.C
+++ b/src/turbulenceModels/RAS/incompressible/wallFunctions/nutWallFunctions/nutStandardWallFunction/nutStandardWallFunctionFvPatchScalarField.C
@@ -108,9 +108,9 @@ void nutStandardWallFunctionFvPatchScalarField::evaluate
const RASModel& rasModel
= db().lookupObject<RASModel>("RASProperties");
- scalar kappa = rasModel.kappa();
- scalar E = rasModel.E();
- scalar yPlusLam = rasModel.E();
+ scalar kappa = rasModel.kappa().value();
+ scalar E = rasModel.E().value();
+ scalar yPlusLam = rasModel.E().value();
const scalarField& ry = patch().deltaCoeffs();
diff --git a/src/turbulenceModels/RAS/incompressible/wallFunctions/nutWallFunctions/nutWallFunction/nutWallFunctionFvPatchScalarField.C b/src/turbulenceModels/RAS/incompressible/wallFunctions/nutWallFunctions/nutWallFunction/nutWallFunctionFvPatchScalarField.C
index 12dc46e96d0e8827dbfd0c67a4e3abc3029cfb40..7ec3523fd7dcb6cb994dc9dfb19378d466a1d234 100644
--- a/src/turbulenceModels/RAS/incompressible/wallFunctions/nutWallFunctions/nutWallFunction/nutWallFunctionFvPatchScalarField.C
+++ b/src/turbulenceModels/RAS/incompressible/wallFunctions/nutWallFunctions/nutWallFunction/nutWallFunctionFvPatchScalarField.C
@@ -100,8 +100,8 @@ void nutWallFunctionFvPatchScalarField::evaluate(const Pstream::commsTypes)
const RASModel& rasModel
= db().lookupObject<RASModel>("RASProperties");
- scalar kappa = rasModel.kappa();
- scalar E = rasModel.E();
+ scalar kappa = rasModel.kappa().value();
+ scalar E = rasModel.E().value();
const scalarField& ry = patch().deltaCoeffs();
diff --git a/src/turbulenceModels/RAS/incompressible/wallFunctions/wallFunctionsI.H b/src/turbulenceModels/RAS/incompressible/wallFunctions/wallFunctionsI.H
index 52c4b0c0e14491ecadd0d779a0aaf6a0e602dc6d..c5f97fc043e9b1f300505d9828df0f3dda75b5e0 100644
--- a/src/turbulenceModels/RAS/incompressible/wallFunctions/wallFunctionsI.H
+++ b/src/turbulenceModels/RAS/incompressible/wallFunctions/wallFunctionsI.H
@@ -88,7 +88,7 @@ Description
epsilon_[faceCelli] +=
Cmu75*pow(k_[faceCelli], 1.5)
- /(kappa_*RASModel::y_[patchi][facei]);
+ /(kappa_.value()*RASModel::y_[patchi][facei]);
if (yPlus > yPlusLam_)
{
@@ -96,7 +96,7 @@ Description
(nutw[facei] + nuw[facei])
*magFaceGradU[facei]
*Cmu25*sqrt(k_[faceCelli])
- /(kappa_*RASModel::y_[patchi][facei]);
+ /(kappa_.value()*RASModel::y_[patchi][facei]);
}
}
}
diff --git a/src/turbulenceModels/RAS/incompressible/wallFunctions/wallNonlinearViscosityI.H b/src/turbulenceModels/RAS/incompressible/wallFunctions/wallNonlinearViscosityI.H
index d1d3617e9257eaede933e20816429350a5c9d818..05bc8411abbb2417102d939a7188f0e00c70d07a 100644
--- a/src/turbulenceModels/RAS/incompressible/wallFunctions/wallNonlinearViscosityI.H
+++ b/src/turbulenceModels/RAS/incompressible/wallFunctions/wallNonlinearViscosityI.H
@@ -55,7 +55,7 @@ Description
if (yPlus > yPlusLam_)
{
nutw[facei] = nuw[facei]
- *(yPlus*kappa_/log(E_*yPlus) - 1);
+ *(yPlus*kappa_.value()/log(E_.value()*yPlus) - 1);
}
else
{
diff --git a/src/turbulenceModels/RAS/incompressible/wallFunctions/wallViscosityI.H b/src/turbulenceModels/RAS/incompressible/wallFunctions/wallViscosityI.H
index 82a47f39b9efb9003dad55847f40d75fb93db1b7..9a7fe62d3e58c88a50ea1100828a48b71307d30d 100644
--- a/src/turbulenceModels/RAS/incompressible/wallFunctions/wallViscosityI.H
+++ b/src/turbulenceModels/RAS/incompressible/wallFunctions/wallViscosityI.H
@@ -55,7 +55,7 @@ Description
{
nutw[facei] =
nuw[facei]
- *(yPlus*kappa_/log(E_*yPlus) - 1);
+ *(yPlus*kappa_.value()/log(E_.value()*yPlus) - 1);
}
else
{