From 6663702c457a17cb59430c8f5754d969f5b63f32 Mon Sep 17 00:00:00 2001 From: andy Date: Fri, 25 Jul 2008 18:01:06 +0100 Subject: [PATCH 1/2] cleaned up output message (removed dims etc) --- .../utilities/postProcessing/wall/yPlusLES/yPlusLES.C | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/applications/utilities/postProcessing/wall/yPlusLES/yPlusLES.C b/applications/utilities/postProcessing/wall/yPlusLES/yPlusLES.C index c5524896c1..d7c145b953 100644 --- a/applications/utilities/postProcessing/wall/yPlusLES/yPlusLES.C +++ b/applications/utilities/postProcessing/wall/yPlusLES/yPlusLES.C @@ -107,11 +107,12 @@ int main(int argc, char *argv[]) *mag(U.boundaryField()[patchi].snGrad()) ) /sgsModel->nu().boundaryField()[patchi]; + const scalarField& Yp = yPlus.boundaryField()[patchi]; Info<< "Patch " << patchi << " named " << currPatch.name() - << " y+ : min: " << min(yPlus) << " max: " << max(yPlus) - << " average: " << average(yPlus) << nl << endl; + << " y+ : min: " << min(Yp) << " max: " << max(Yp) + << " average: " << average(Yp) << nl << endl; } } From e97b721db785ecd4ceb7fee8c5342d86e09f0129 Mon Sep 17 00:00:00 2001 From: andy Date: Wed, 30 Jul 2008 15:32:51 +0100 Subject: [PATCH 2/2] bugfix for lambda_ --- .../kineticTheoryModel/kineticTheoryModel.C | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C b/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C index 54590f55a3..724602009a 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C @@ -256,8 +256,7 @@ void Foam::kineticTheoryModel::solve() /(alpha*rhoa_*sqrtPi*(ThetaSqrt + TsmallSqrt)); // bulk viscosity p. 45 (Lun et al. 1984). - volScalarField lambda_ = - (4.0/3.0)*sqr(alpha_)*rhoa_*da_*gs0_*(1.0+e_)*ThetaSqrt/sqrtPi; + lambda_ = (4.0/3.0)*sqr(alpha_)*rhoa_*da_*gs0_*(1.0+e_)*ThetaSqrt/sqrtPi; // stress tensor, Definitions, Table 3.1, p. 43