From 2aaf0fb67051bfb1e00540a345045828df2473d5 Mon Sep 17 00:00:00 2001 From: henry Date: Thu, 30 Jul 2009 15:56:41 +0100 Subject: [PATCH] Force construction of mesh.deltaCoeffs() which may be needed for the wall-function BCs. --- src/turbulenceModels/compressible/LES/LESModel/LESModel.C | 4 ++++ src/turbulenceModels/compressible/RAS/RASModel/RASModel.C | 6 +++++- src/turbulenceModels/incompressible/LES/LESModel/LESModel.C | 4 ++++ .../nuSgsWallFunction/nuSgsWallFunctionFvPatchScalarField.C | 4 +++- src/turbulenceModels/incompressible/RAS/RASModel/RASModel.C | 6 +++++- 5 files changed, 21 insertions(+), 3 deletions(-) diff --git a/src/turbulenceModels/compressible/LES/LESModel/LESModel.C b/src/turbulenceModels/compressible/LES/LESModel/LESModel.C index 5602e44d..dd166f26 100644 --- a/src/turbulenceModels/compressible/LES/LESModel/LESModel.C +++ b/src/turbulenceModels/compressible/LES/LESModel/LESModel.C @@ -84,6 +84,10 @@ LESModel::LESModel delta_(LESdelta::New("delta", U.mesh(), *this)) { readIfPresent("k0", k0_); + + // Force the construction of the mesh deltaCoeffs which may be needed + // for the construction of the derived models and BCs + mesh_.deltaCoeffs(); } diff --git a/src/turbulenceModels/compressible/RAS/RASModel/RASModel.C b/src/turbulenceModels/compressible/RAS/RASModel/RASModel.C index 8d18cec8..b390bca5 100644 --- a/src/turbulenceModels/compressible/RAS/RASModel/RASModel.C +++ b/src/turbulenceModels/compressible/RAS/RASModel/RASModel.C @@ -88,7 +88,11 @@ RASModel::RASModel omegaSmall_("omegaSmall", omega0_.dimensions(), SMALL), y_(mesh_) -{} +{ + // Force the construction of the mesh deltaCoeffs which may be needed + // for the construction of the derived models and BCs + mesh_.deltaCoeffs(); +} // * * * * * * * * * * * * * * * * Selectors * * * * * * * * * * * * * * * * // diff --git a/src/turbulenceModels/incompressible/LES/LESModel/LESModel.C b/src/turbulenceModels/incompressible/LES/LESModel/LESModel.C index 3e3c952f..3458e03c 100644 --- a/src/turbulenceModels/incompressible/LES/LESModel/LESModel.C +++ b/src/turbulenceModels/incompressible/LES/LESModel/LESModel.C @@ -82,6 +82,10 @@ LESModel::LESModel delta_(LESdelta::New("delta", U.mesh(), *this)) { readIfPresent("k0", k0_); + + // Force the construction of the mesh deltaCoeffs which may be needed + // for the construction of the derived models and BCs + mesh_.deltaCoeffs(); } diff --git a/src/turbulenceModels/incompressible/LES/derivedFvPatchFields/wallFunctions/nuSgsWallFunctions/nuSgsWallFunction/nuSgsWallFunctionFvPatchScalarField.C b/src/turbulenceModels/incompressible/LES/derivedFvPatchFields/wallFunctions/nuSgsWallFunctions/nuSgsWallFunction/nuSgsWallFunctionFvPatchScalarField.C index 5555fd30..1d66afbc 100644 --- a/src/turbulenceModels/incompressible/LES/derivedFvPatchFields/wallFunctions/nuSgsWallFunctions/nuSgsWallFunction/nuSgsWallFunctionFvPatchScalarField.C +++ b/src/turbulenceModels/incompressible/LES/derivedFvPatchFields/wallFunctions/nuSgsWallFunctions/nuSgsWallFunction/nuSgsWallFunctionFvPatchScalarField.C @@ -133,8 +133,8 @@ void nuSgsWallFunctionFvPatchScalarField::evaluate const scalarField& nuw = patch().lookupPatchField(nuName_); - scalarField& nuSgsw = *this; + scalarField& nuSgsw = *this; scalarField magFaceGradU = mag(U.snGrad()); @@ -178,6 +178,8 @@ void nuSgsWallFunctionFvPatchScalarField::evaluate nuSgsw[facei] = 0; } } + + fixedValueFvPatchScalarField::evaluate(); } diff --git a/src/turbulenceModels/incompressible/RAS/RASModel/RASModel.C b/src/turbulenceModels/incompressible/RAS/RASModel/RASModel.C index 6b300390..d3366d82 100644 --- a/src/turbulenceModels/incompressible/RAS/RASModel/RASModel.C +++ b/src/turbulenceModels/incompressible/RAS/RASModel/RASModel.C @@ -87,7 +87,11 @@ RASModel::RASModel omegaSmall_("omegaSmall", omega0_.dimensions(), SMALL), y_(mesh_) -{} +{ + // Force the construction of the mesh deltaCoeffs which may be needed + // for the construction of the derived models and BCs + mesh_.deltaCoeffs(); +} // * * * * * * * * * * * * * * * * * Selectors * * * * * * * * * * * * * * * // -- 2.11.4.GIT