diff --git a/src/turbulenceModels/LES/LESdeltas/maxhxhyhzDelta/maxhxhyhzDelta.C b/src/turbulenceModels/LES/LESdeltas/maxhxhyhzDelta/maxhxhyhzDelta.C deleted file mode 100644 index 0ba1d63b1aa42642dcadf10a0e81f86d734096c5..0000000000000000000000000000000000000000 --- a/src/turbulenceModels/LES/LESdeltas/maxhxhyhzDelta/maxhxhyhzDelta.C +++ /dev/null @@ -1,142 +0,0 @@ -/*---------------------------------------------------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | - \\ / A nd | Copyright (C) 2010-2010 OpenCFD Ltd. - \\/ M anipulation | -------------------------------------------------------------------------------- -License - This file is part of OpenFOAM. - - OpenFOAM is free software: you can redistribute it and/or modify it - under the terms of the GNU General Public License as published by - the Free Software Foundation, either version 3 of the License, or - (at your option) any later version. - - OpenFOAM is distributed in the hope that it will be useful, but WITHOUT - ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or - FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License - for more details. - - You should have received a copy of the GNU General Public License - along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>. - -\*---------------------------------------------------------------------------*/ - -#include "maxhxhyhzDelta.H" -#include "addToRunTimeSelectionTable.H" - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -namespace Foam -{ - -// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * // - -defineTypeNameAndDebug(maxhxhyhzDelta, 0); -addToRunTimeSelectionTable(LESdelta, maxhxhyhzDelta, dictionary); - -// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * // - -void maxhxhyhzDelta::calcDelta() -{ - label nD = mesh().nGeometricD(); - - tmp<volScalarField> hmax - ( - new volScalarField - ( - IOobject - ( - "hmax", - mesh().time().timeName(), - mesh(), - IOobject::NO_READ, - IOobject::NO_WRITE - ), - mesh(), - dimensionedScalar("zrero", dimLength, 0.0) - ) - ); - - const cellList& cells = mesh().cells(); - - forAll(cells,cellI) - { - scalar deltaMaxTmp = 0.0; - const labelList& cFaces = mesh().cells()[cellI]; - const point& centrevector = mesh().cellCentres()[cellI]; - - forAll(cFaces, cFaceI) - { - label faceI = cFaces[cFaceI]; - const point& facevector = mesh().faceCentres()[faceI]; - scalar tmp = mag(facevector - centrevector); - if (tmp > deltaMaxTmp) - { - deltaMaxTmp = tmp; - } - } - hmax()[cellI] = deltaCoeff_*deltaMaxTmp; - } - - if (nD == 3) - { - delta_.internalField() = hmax(); - } - else if (nD == 2) - { - WarningIn("maxhxhyhzDelta::calcDelta()") - << "Case is 2D, LES is not strictly applicable\n" - << endl; - - delta_.internalField() = hmax(); - } - else - { - FatalErrorIn("maxhxhyhzDelta::calcDelta()") - << "Case is not 3D or 2D, LES is not applicable" - << exit(FatalError); - } -} - - -// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // - -maxhxhyhzDelta::maxhxhyhzDelta -( - const word& name, - const fvMesh& mesh, - const dictionary& dd -) -: - LESdelta(name, mesh), - deltaCoeff_(readScalar(dd.subDict(type() + "Coeffs").lookup("deltaCoeff"))) -{ - calcDelta(); -} - - -// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // - -void maxhxhyhzDelta::read(const dictionary& dd) -{ - dd.subDict(type() + "Coeffs").lookup("deltaCoeff") >> deltaCoeff_; - calcDelta(); -} - - -void maxhxhyhzDelta::correct() -{ - if (mesh_.changing()) - { - calcDelta(); - } -} - - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -} // End namespace Foam - -// ************************************************************************* // diff --git a/src/turbulenceModels/LES/LESdeltas/maxhxhyhzDelta/maxhxhyhzDelta.H b/src/turbulenceModels/LES/LESdeltas/maxhxhyhzDelta/maxhxhyhzDelta.H deleted file mode 100644 index 30f26e6bd09ef902fe2f9dbc2f12ee6b3e307d78..0000000000000000000000000000000000000000 --- a/src/turbulenceModels/LES/LESdeltas/maxhxhyhzDelta/maxhxhyhzDelta.H +++ /dev/null @@ -1,111 +0,0 @@ -/*---------------------------------------------------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | - \\ / A nd | Copyright (C) 2010-2010 OpenCFD Ltd. - \\/ M anipulation | -------------------------------------------------------------------------------- -License - This file is part of OpenFOAM. - - OpenFOAM is free software: you can redistribute it and/or modify it - under the terms of the GNU General Public License as published by - the Free Software Foundation, either version 3 of the License, or - (at your option) any later version. - - OpenFOAM is distributed in the hope that it will be useful, but WITHOUT - ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or - FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License - for more details. - - You should have received a copy of the GNU General Public License - along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>. - -Class - Foam::maxhxhyhzDelta - -Description - maxhxhyhzDelta takes the maximum of the three dimensions per cell: - max(hx, hy, hz). Valid for structures hexahedral cells only. - - -SourceFiles - maxhxhyhzDelta.C - -\*---------------------------------------------------------------------------*/ - -#ifndef maxhxhyhzDeltaDelta_H -#define maxhxhyhzDeltaDelta_H - -#include "LESdelta.H" - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -namespace Foam -{ - -/*---------------------------------------------------------------------------*\ - Class maxhxhyhzDelta Declaration -\*---------------------------------------------------------------------------*/ - -class maxhxhyhzDelta -: - public LESdelta -{ - // Private data - - scalar deltaCoeff_; // - - - // Private Member Functions - - //- Disallow default bitwise copy construct and assignment - maxhxhyhzDelta(const maxhxhyhzDelta&); - void operator=(const maxhxhyhzDelta&); - - // Calculate the delta values - void calcDelta(); - - -public: - - //- Runtime type information - TypeName("maxhxhyhzDelta"); - - - // Constructors - - //- Construct from name, mesh and IOdictionary - maxhxhyhzDelta - ( - const word& name, - const fvMesh& mesh, - const dictionary& - ); - - - //- Destructor - virtual ~maxhxhyhzDelta() - {} - - - // Member Functions - - //- Read the LESdelta dictionary - virtual void read(const dictionary&); - - // Correct values - virtual void correct(); -}; - - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -} // End namespace Foam - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - - -#endif - -// ************************************************************************* //