From e0e369d8889b1a78d877d25ca9754811da705a9e Mon Sep 17 00:00:00 2001 From: Sergio Ferraris Date: Tue, 14 May 2013 09:57:32 +0100 Subject: [PATCH] BUG:Correcting expression for the inertial contribution to the pressure drop in the bafflePressureDrop --- .../porousBafflePressure/porousBafflePressureFvPatchField.H | 4 ++-- .../porousBafflePressure/porousBafflePressureFvPatchFields.C | 4 ++-- 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/src/turbulenceModels/derivedFvPatchFields/porousBafflePressure/porousBafflePressureFvPatchField.H b/src/turbulenceModels/derivedFvPatchFields/porousBafflePressure/porousBafflePressureFvPatchField.H index a52e8e5b34..ed22575785 100644 --- a/src/turbulenceModels/derivedFvPatchFields/porousBafflePressure/porousBafflePressureFvPatchField.H +++ b/src/turbulenceModels/derivedFvPatchFields/porousBafflePressure/porousBafflePressureFvPatchField.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -34,7 +34,7 @@ Description The porous baffle introduces a pressure jump defined by: \f[ - \Delta p = -(I \mu U + 0.5 D \rho |U|^2 L) + \Delta p = -(I \mu U + 0.5 D \rho |U|^2 )L \f] where diff --git a/src/turbulenceModels/derivedFvPatchFields/porousBafflePressure/porousBafflePressureFvPatchFields.C b/src/turbulenceModels/derivedFvPatchFields/porousBafflePressure/porousBafflePressureFvPatchFields.C index 08855fa7fb..32317c0ad1 100644 --- a/src/turbulenceModels/derivedFvPatchFields/porousBafflePressure/porousBafflePressureFvPatchFields.C +++ b/src/turbulenceModels/derivedFvPatchFields/porousBafflePressure/porousBafflePressureFvPatchFields.C @@ -78,7 +78,7 @@ void Foam::porousBafflePressureFvPatchField::updateCoeffs() const scalarField nuEffw = turbModel.nuEff()().boundaryField()[patchI]; - jump_ = -sign(Un)*(I_*nuEffw + D_*0.5*magUn*length_)*magUn; + jump_ = -sign(Un)*(I_*nuEffw + D_*0.5*magUn)*magUn*length_; } else { @@ -95,7 +95,7 @@ void Foam::porousBafflePressureFvPatchField::updateCoeffs() Un /= rhow; - jump_ = -sign(Un)*(I_*muEffw + D_*0.5*rhow*magUn*length_)*magUn; + jump_ = -sign(Un)*(I_*muEffw + D_*0.5*rhow*magUn)*magUn*length_; } if (debug)