BUG: changing from turbulent to lamninar viscocity in the calculation of the

Darcy coefficient for the porousBafflePressure
This commit is contained in:
sergio
2014-03-14 14:56:11 +00:00
committed by Andrew Heather
parent 5844801027
commit 354197bf33
10 changed files with 2293 additions and 24 deletions

View File

@ -34,7 +34,7 @@ Description
The porous baffle introduces a pressure jump defined by: The porous baffle introduces a pressure jump defined by:
\f[ \f[
\Delta p = -(I \mu U + 0.5 D \rho |U|^2 )L \Delta p = -(D \mu U + 0.5 I \rho |U|^2 )L
\f] \f]
where where
@ -42,7 +42,7 @@ Description
\vartable \vartable
p | pressure [Pa] p | pressure [Pa]
\rho | density [kg/m3] \rho | density [kg/m3]
\mu | viscosity [Pa s] \mu | laminar viscosity [Pa s]
I | inertial coefficient I | inertial coefficient
D | Darcy coefficient D | Darcy coefficient
L | porous media length in the flow direction L | porous media length in the flow direction
@ -66,8 +66,8 @@ Description
type porousBafflePressure; type porousBafflePressure;
patchType cyclic; patchType cyclic;
jump uniform 0; jump uniform 0;
D 1000000; I 1000000;
I 0.001; D 0.001;
L 0.1; L 0.1;
value uniform 0; value uniform 0;
} }

View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -76,9 +76,9 @@ void Foam::porousBafflePressureFvPatchField<Foam::scalar>::updateCoeffs()
"turbulenceModel" "turbulenceModel"
); );
const scalarField nuEffw = turbModel.nuEff()().boundaryField()[patchI]; const scalarField nu = turbModel.nu()().boundaryField()[patchI];
jump_ = -sign(Un)*(I_*nuEffw + D_*0.5*magUn)*magUn*length_; jump_ = -sign(Un)*(D_*nu + I_*0.5*magUn)*magUn*length_;
} }
else else
{ {
@ -88,14 +88,14 @@ void Foam::porousBafflePressureFvPatchField<Foam::scalar>::updateCoeffs()
"turbulenceModel" "turbulenceModel"
); );
const scalarField muEffw = turbModel.muEff()().boundaryField()[patchI]; const scalarField mu = turbModel.mu().boundaryField()[patchI];
const scalarField rhow = const scalarField rhow =
patch().lookupPatchField<volScalarField, scalar>("rho"); patch().lookupPatchField<volScalarField, scalar>("rho");
Un /= rhow; Un /= rhow;
jump_ = -sign(Un)*(I_*muEffw + D_*0.5*rhow*magUn)*magUn*length_; jump_ = -sign(Un)*(D_*mu + I_*0.5*rhow*magUn)*magUn*length_;
} }
if (debug) if (debug)

View File

@ -1,7 +1,7 @@
/*--------------------------------*- C++ -*----------------------------------*\ /*--------------------------------*- C++ -*----------------------------------*\
| ========= | | | ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 | | \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org | | \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | | | \\/ M anipulation | |
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/

View File

@ -1,7 +1,7 @@
/*--------------------------------*- C++ -*----------------------------------*\ /*--------------------------------*- C++ -*----------------------------------*\
| ========= | | | ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 | | \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org | | \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | | | \\/ M anipulation | |
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/

View File

@ -1,7 +1,7 @@
/*--------------------------------*- C++ -*----------------------------------*\ /*--------------------------------*- C++ -*----------------------------------*\
| ========= | | | ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 | | \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org | | \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | | | \\/ M anipulation | |
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
@ -28,7 +28,6 @@ boundaryField
Cmu 0.09; Cmu 0.09;
kappa 0.41; kappa 0.41;
E 9.8; E 9.8;
value uniform 0.1;
} }
rightWall rightWall
{ {
@ -37,7 +36,6 @@ boundaryField
Cmu 0.09; Cmu 0.09;
kappa 0.41; kappa 0.41;
E 9.8; E 9.8;
value uniform 0.1;
} }
lowerWall lowerWall
{ {
@ -46,7 +44,6 @@ boundaryField
Cmu 0.09; Cmu 0.09;
kappa 0.41; kappa 0.41;
E 9.8; E 9.8;
value uniform 0.1;
} }
atmosphere atmosphere
{ {

View File

@ -1,7 +1,7 @@
/*--------------------------------*- C++ -*----------------------------------*\ /*--------------------------------*- C++ -*----------------------------------*\
| ========= | | | ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 | | \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org | | \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | | | \\/ M anipulation | |
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/

View File

@ -1,7 +1,7 @@
/*--------------------------------*- C++ -*----------------------------------*\ /*--------------------------------*- C++ -*----------------------------------*\
| ========= | | | ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 | | \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org | | \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | | | \\/ M anipulation | |
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/

View File

@ -1,7 +1,7 @@
/*--------------------------------*- C++ -*----------------------------------*\ /*--------------------------------*- C++ -*----------------------------------*\
| ========= | | | ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 | | \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org | | \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | | | \\/ M anipulation | |
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/

View File

@ -1,7 +1,7 @@
/*--------------------------------*- C++ -*----------------------------------*\ /*--------------------------------*- C++ -*----------------------------------*\
| ========= | | | ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 | | \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org | | \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | | | \\/ M anipulation | |
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
@ -59,7 +59,7 @@ boundaryField
jump uniform 0; jump uniform 0;
value uniform 0; value uniform 0;
D 1000; D 1000;
I 500; I 1000;
length 0.15; length 0.15;
} }
porous_half1 porous_half1
@ -68,7 +68,7 @@ boundaryField
patchType cyclic; patchType cyclic;
value uniform 0; value uniform 0;
D 1000; D 1000;
I 500; I 1000;
length 0.15; length 0.15;
} }
} }