diff --git a/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C b/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
index 4c16ae0e541b6fdccd4adaa916ceca13583beef8..18fa6a61ef6cf6294c2eb905c28a165af914e3ca 100644
--- a/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
+++ b/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -104,13 +104,15 @@ int main(int argc, char *argv[])
Info<< "Setting boundary layer velocity" << nl << endl;
scalar yblv = ybl.value();
- forAll(U, celli)
+ forAll(U, cellI)
{
- if (y[celli] <= yblv)
+ if (y[cellI] <= yblv)
{
- U[celli] *= ::pow(y[celli]/yblv, (1.0/7.0));
+ mask[cellI] = 1;
+ U[cellI] *= ::pow(y[cellI]/yblv, (1.0/7.0));
}
}
+ mask.correctBoundaryConditions();
Info<< "Writing U\n" << endl;
U.write();
@@ -128,11 +130,15 @@ int main(int argc, char *argv[])
if (isA<incompressible::RASModel>(turbulence()))
{
- // Calculate nut
+ // Calculate nut - reference nut is calculated by the turbulence model
+ // on its construction
tmp<volScalarField> tnut = turbulence->nut();
volScalarField& nut = tnut();
volScalarField S(mag(dev(symm(fvc::grad(U)))));
- nut = sqr(kappa*min(y, ybl))*::sqrt(2)*S;
+ nut = (1 - mask)*nut + mask*sqr(kappa*min(y, ybl))*::sqrt(2)*S;
+
+ // do not correct BC - wall functions will 'undo' manipulation above
+ // by using nut from turbulence model
if (args.optionFound("writenut"))
{
@@ -140,9 +146,6 @@ int main(int argc, char *argv[])
nut.write();
}
- // Create G field - used by RAS wall functions
- volScalarField G(turbulence().GName(), nut*2*sqr(S));
-
//--- Read and modify turbulence fields
@@ -150,8 +153,11 @@ int main(int argc, char *argv[])
tmp<volScalarField> tk = turbulence->k();
volScalarField& k = tk();
scalar ck0 = pow025(Cmu)*kappa;
- k = sqr(nut/(ck0*min(y, ybl)));
- k.correctBoundaryConditions();
+ k = (1 - mask)*k + mask*sqr(nut/(ck0*min(y, ybl)));
+
+ // do not correct BC - operation may use inconsistent fields wrt these
+ // local manipulations
+ // k.correctBoundaryConditions();
Info<< "Writing k\n" << endl;
k.write();
@@ -161,8 +167,11 @@ int main(int argc, char *argv[])
tmp<volScalarField> tepsilon = turbulence->epsilon();
volScalarField& epsilon = tepsilon();
scalar ce0 = ::pow(Cmu, 0.75)/kappa;
- epsilon = ce0*k*sqrt(k)/min(y, ybl);
- epsilon.correctBoundaryConditions();
+ epsilon = (1 - mask)*epsilon + mask*ce0*k*sqrt(k)/min(y, ybl);
+
+ // do not correct BC - wall functions will use non-updated k from
+ // turbulence model
+ // epsilon.correctBoundaryConditions();
Info<< "Writing epsilon\n" << endl;
epsilon.write();
@@ -181,12 +190,12 @@ int main(int argc, char *argv[])
if (omegaHeader.headerOk())
{
volScalarField omega(omegaHeader, mesh);
- omega =
- epsilon
- /(
- Cmu*k+dimensionedScalar("VSMALL", k.dimensions(), VSMALL)
- );
- omega.correctBoundaryConditions();
+ dimensionedScalar k0("VSMALL", k.dimensions(), VSMALL);
+ omega = (1 - mask)*omega + mask*epsilon/(Cmu*k + k0);
+
+ // do not correct BC - wall functions will use non-updated k from
+ // turbulence model
+ // omega.correctBoundaryConditions();
Info<< "Writing omega\n" << endl;
omega.write();
@@ -207,7 +216,9 @@ int main(int argc, char *argv[])
{
volScalarField nuTilda(nuTildaHeader, mesh);
nuTilda = nut;
- nuTilda.correctBoundaryConditions();
+
+ // do not correct BC
+ // nuTilda.correctBoundaryConditions();
Info<< "Writing nuTilda\n" << endl;
nuTilda.write();
diff --git a/applications/utilities/preProcessing/applyBoundaryLayer/createFields.H b/applications/utilities/preProcessing/applyBoundaryLayer/createFields.H
index e7d72692fc1f640e0b25c6023d0eaccf729f61b4..22515ac83270fdd31dd712aa50e3dd4b214208ab 100644
--- a/applications/utilities/preProcessing/applyBoundaryLayer/createFields.H
+++ b/applications/utilities/preProcessing/applyBoundaryLayer/createFields.H
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -23,39 +23,55 @@ License
\*---------------------------------------------------------------------------*/
- Info<< "Reading field U\n" << endl;
- volVectorField U
+Info<< "Reading field U\n" << endl;
+volVectorField U
+(
+ IOobject
(
- IOobject
- (
- "U",
- runTime.timeName(),
- mesh,
- IOobject::MUST_READ,
- IOobject::NO_WRITE
- ),
- mesh
- );
-
- Info<< "Calculating wall distance field" << endl;
- volScalarField y(wallDist(mesh).y());
-
- // Set the mean boundary-layer thickness
- dimensionedScalar ybl("ybl", dimLength, 0);
-
- if (args.optionFound("ybl"))
- {
- // If the boundary-layer thickness is provided use it
- ybl.value() = args.optionRead<scalar>("ybl");
- }
- else if (args.optionFound("Cbl"))
- {
- // Calculate boundary layer thickness as Cbl * mean distance to wall
- ybl.value() = gAverage(y) * args.optionRead<scalar>("Cbl");
- }
-
- Info<< "\nCreating boundary-layer for U of thickness "
- << ybl.value() << " m" << nl << endl;
+ "U",
+ runTime.timeName(),
+ mesh,
+ IOobject::MUST_READ,
+ IOobject::NO_WRITE
+ ),
+ mesh
+);
+
+Info<< "Calculating wall distance field" << endl;
+volScalarField y(wallDist(mesh).y());
+
+// Set the mean boundary-layer thickness
+dimensionedScalar ybl("ybl", dimLength, 0);
+
+if (args.optionFound("ybl"))
+{
+ // If the boundary-layer thickness is provided use it
+ ybl.value() = args.optionRead<scalar>("ybl");
+}
+else if (args.optionFound("Cbl"))
+{
+ // Calculate boundary layer thickness as Cbl*mean distance to wall
+ ybl.value() = gAverage(y)*args.optionRead<scalar>("Cbl");
+}
+
+Info<< "\nCreating boundary-layer for U of thickness "
+ << ybl.value() << " m" << nl << endl;
+
+Info<< "Creating mask field" << endl;
+volScalarField mask
+(
+ IOobject
+ (
+ "mask",
+ runTime.timeName(),
+ mesh,
+ IOobject::NO_READ,
+ IOobject::NO_WRITE
+ ),
+ mesh,
+ dimensionedScalar("zero", dimless, 0.0),
+ zeroGradientFvPatchScalarField::typeName
+);
// ************************************************************************* //