From d02c8d1a52b6005e766353bc7340b35f520153a9 Mon Sep 17 00:00:00 2001 From: andy Date: Thu, 13 Aug 2009 09:44:09 +0100 Subject: [PATCH] updated Up calc to use mag() --- .../solvers/incompressible/boundaryFoam/evaluateNearWall.H | 7 +++---- 1 file changed, 3 insertions(+), 4 deletions(-) diff --git a/applications/solvers/incompressible/boundaryFoam/evaluateNearWall.H b/applications/solvers/incompressible/boundaryFoam/evaluateNearWall.H index b9683b9e68..8844f21e6e 100644 --- a/applications/solvers/incompressible/boundaryFoam/evaluateNearWall.H +++ b/applications/solvers/incompressible/boundaryFoam/evaluateNearWall.H @@ -7,8 +7,7 @@ scalar epsilon = turbulence->epsilon()()[cellId]; // scalar omega = turbulence->omega()()[cellId]; scalar k = turbulence->k()()[cellId]; - scalar Up = - flowDirection & (U[cellId] - U.boundaryField()[patchId][faceId]); + scalar magUp = mag(U[cellId] - U.boundaryField()[patchId][faceId]); scalar tauw = flowDirection & R & wallNormal; @@ -16,7 +15,7 @@ scalar yPlus = uTau*y[cellId]/(nu + ROOTVSMALL); - scalar uPlus = Up/(uTau + ROOTVSMALL); + scalar uPlus = magUp/(uTau + ROOTVSMALL); scalar nutPlus = nut/nu; @@ -26,7 +25,7 @@ // scalar omegaPlus = omega*nu/(sqr(uTau) + ROOTVSMALL); - scalar Rey = Up*y[cellId]/nu; + scalar Rey = magUp*y[cellId]/nu; Info<< "Rey = " << Rey << ", uTau = " << uTau << ", nut+ = " << nutPlus << ", y+ = " << yPlus << ", u+ = " << uPlus