From 3f09e6e3b3df8cd4fb4e172675a827075d9ab121 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Thu, 21 Mar 2013 17:04:42 +0000
Subject: [PATCH] anisotropicFilter: Corrected formulation of the anisotropic
coefficient
---
.../anisotropicFilter/anisotropicFilter.C | 18 +++++++++++++-----
.../LESfilters/laplaceFilter/laplaceFilter.C | 6 +++---
.../LES/dynOneEqEddy/dynOneEqEddy.C | 4 ++--
3 files changed, 18 insertions(+), 10 deletions(-)
diff --git a/src/turbulenceModels/LES/LESfilters/anisotropicFilter/anisotropicFilter.C b/src/turbulenceModels/LES/LESfilters/anisotropicFilter/anisotropicFilter.C
index 060ed5bda3c..08ce6197c65 100644
--- a/src/turbulenceModels/LES/LESfilters/anisotropicFilter/anisotropicFilter.C
+++ b/src/turbulenceModels/LES/LESfilters/anisotropicFilter/anisotropicFilter.C
@@ -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()
- /fvc::surfaceSum(mag(mesh.Sf().component(d)))().internalField()
+ (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()
- /fvc::surfaceSum(mag(mesh.Sf().component(d)))().internalField()
+ (1/widthCoeff_)*
+ sqr
+ (
+ 2.0*mesh.V()
+ /fvc::surfaceSum(mag(mesh.Sf().component(d)))().internalField()
+ )
);
}
}
diff --git a/src/turbulenceModels/LES/LESfilters/laplaceFilter/laplaceFilter.C b/src/turbulenceModels/LES/LESfilters/laplaceFilter/laplaceFilter.C
index 231ae2737fa..24db120ffa9 100644
--- a/src/turbulenceModels/LES/LESfilters/laplaceFilter/laplaceFilter.C
+++ b/src/turbulenceModels/LES/LESfilters/laplaceFilter/laplaceFilter.C
@@ -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_;
}
diff --git a/src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.C b/src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.C
index 40c4deaf33a..5e125d53108 100644
--- a/src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.C
+++ b/src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.C
@@ -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
--
GitLab