anisotropicFilter: Corrected formulation of the anisotropic coefficient

This commit is contained in:
Henry
2013-03-21 17:04:42 +00:00
parent af0378317e
commit 3f09e6e3b3
3 changed files with 18 additions and 10 deletions

View File

@ -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
@ -66,8 +66,12 @@ Foam::anisotropicFilter::anisotropicFilter
coeff_.internalField().replace
(
d,
(2.0/widthCoeff_)*mesh.V()
(1/widthCoeff_)*
sqr
(
2.0*mesh.V()
/fvc::surfaceSum(mag(mesh.Sf().component(d)))().internalField()
)
);
}
}
@ -99,8 +103,12 @@ Foam::anisotropicFilter::anisotropicFilter
coeff_.internalField().replace
(
d,
(2.0/widthCoeff_)*mesh.V()
(1/widthCoeff_)*
sqr
(
2.0*mesh.V()
/fvc::surfaceSum(mag(mesh.Sf().component(d)))().internalField()
)
);
}
}

View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -57,7 +57,7 @@ Foam::laplaceFilter::laplaceFilter(const fvMesh& mesh, scalar widthCoeff)
calculatedFvPatchScalarField::typeName
)
{
coeff_.internalField() = pow(mesh.V(), 2.0/3.0)/widthCoeff_;
coeff_.dimensionedInternalField() = pow(mesh.V(), 2.0/3.0)/widthCoeff_;
}
@ -78,7 +78,7 @@ Foam::laplaceFilter::laplaceFilter(const fvMesh& mesh, const dictionary& bd)
calculatedFvPatchScalarField::typeName
)
{
coeff_.internalField() = pow(mesh.V(), 2.0/3.0)/widthCoeff_;
coeff_.dimensionedInternalField() = pow(mesh.V(), 2.0/3.0)/widthCoeff_;
}

View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -66,7 +66,7 @@ volScalarField dynOneEqEddy::ck
const volSymmTensorField MM
(
simpleFilter_(-2.0*delta()*pow(KK, 0.5)*filter_(D))
simpleFilter_(-2.0*delta()*sqrt(KK)*filter_(D))
);
const volScalarField ck