mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-11-28 03:28:01 +00:00
ENH: createZeroFolder renamed to createZeroDirectory
This commit is contained in:
andy
@ -0,0 +1,452 @@
|
||||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2014 OpenFOAM Foundation
|
||||
\\/ 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 "solverTemplate.H"
|
||||
#include "Time.H"
|
||||
#include "IOPtrList.H"
|
||||
#include "polyMesh.H"
|
||||
#include "regionProperties.H"
|
||||
|
||||
using namespace Foam;
|
||||
|
||||
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
template<>
|
||||
const char* Foam::NamedEnum<Foam::solverTemplate::solverType, 4>::names[] =
|
||||
{
|
||||
"compressible",
|
||||
"incompressible",
|
||||
"buoyant",
|
||||
"unknown"
|
||||
};
|
||||
}
|
||||
|
||||
const Foam::NamedEnum<Foam::solverTemplate::solverType, 4>
|
||||
Foam::solverTemplate::solverTypeNames_;
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||
|
||||
Foam::word Foam::solverTemplate::readFromDict
|
||||
(
|
||||
IOobject& dictHeader,
|
||||
const word& entryName
|
||||
) const
|
||||
{
|
||||
if (!dictHeader.headerOk())
|
||||
{
|
||||
FatalErrorIn
|
||||
(
|
||||
"Foam::word Foam::solverTemplate::readFromDict"
|
||||
"("
|
||||
"IOobject&, "
|
||||
"const word&"
|
||||
") const"
|
||||
)
|
||||
<< "Unable to open file "
|
||||
<< dictHeader.objectPath()
|
||||
<< exit(FatalError);
|
||||
}
|
||||
|
||||
IOdictionary dict(dictHeader);
|
||||
return dict.lookup(entryName);
|
||||
}
|
||||
|
||||
|
||||
Foam::dictionary Foam::solverTemplate::readFluidFieldTemplates
|
||||
(
|
||||
const word& regionName,
|
||||
const fileName& baseDir,
|
||||
const dictionary& solverDict,
|
||||
const Time& runTime
|
||||
) const
|
||||
{
|
||||
Info<< " Reading fluid field templates";
|
||||
|
||||
if (regionName == word::null)
|
||||
{
|
||||
Info<< endl;
|
||||
}
|
||||
else
|
||||
{
|
||||
Info<< " for region " << regionName << endl;
|
||||
}
|
||||
|
||||
dictionary fieldTemplates = solverDict.subDict("fluidFields");
|
||||
|
||||
const fileName turbModelDir(baseDir/"models"/"turbulence");
|
||||
word turbulenceModel("laminar"); // default to laminar
|
||||
|
||||
const dictionary fieldModels(solverDict.subDict("fluidModels"));
|
||||
|
||||
word turbulenceType = "none";
|
||||
if (fieldModels.readIfPresent("turbulenceModel", turbulenceType))
|
||||
{
|
||||
word simulationType(word::null);
|
||||
|
||||
IOobject RASHeader
|
||||
(
|
||||
"RASProperties",
|
||||
runTime.constant(),
|
||||
regionName,
|
||||
runTime,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::NO_WRITE
|
||||
);
|
||||
|
||||
IOobject LESHeader
|
||||
(
|
||||
"LESProperties",
|
||||
runTime.constant(),
|
||||
regionName,
|
||||
runTime,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::NO_WRITE
|
||||
);
|
||||
|
||||
|
||||
if (turbulenceType == "turbulenceModel")
|
||||
{
|
||||
IOobject turbulencePropertiesHeader
|
||||
(
|
||||
"turbulenceProperties",
|
||||
runTime.constant(),
|
||||
regionName,
|
||||
runTime,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::NO_WRITE
|
||||
);
|
||||
|
||||
if (turbulencePropertiesHeader.headerOk())
|
||||
{
|
||||
const IOdictionary turbulenceProperties
|
||||
(
|
||||
turbulencePropertiesHeader
|
||||
);
|
||||
|
||||
turbulenceProperties.lookup("simulationType") >> simulationType;
|
||||
|
||||
if (simulationType == "RASModel")
|
||||
{
|
||||
turbulenceModel = readFromDict(RASHeader, "RASModel");
|
||||
}
|
||||
else if (simulationType == "LESModel")
|
||||
{
|
||||
turbulenceModel = readFromDict(LESHeader, "LESModel");
|
||||
}
|
||||
}
|
||||
else
|
||||
{
|
||||
FatalError<< "Unable to open file "
|
||||
<< turbulencePropertiesHeader.objectPath()
|
||||
<< exit(FatalError);
|
||||
}
|
||||
}
|
||||
else if (turbulenceType == "RASModel")
|
||||
{
|
||||
turbulenceModel = readFromDict(RASHeader, "RASModel");
|
||||
}
|
||||
else if (turbulenceType == "LESModel")
|
||||
{
|
||||
turbulenceModel = readFromDict(LESHeader, "LESModel");
|
||||
}
|
||||
else
|
||||
{
|
||||
FatalErrorIn
|
||||
(
|
||||
"Foam::dictionary Foam::solverTemplate::readFluidFieldTemplates"
|
||||
"("
|
||||
"const word&, "
|
||||
"const fileName&, "
|
||||
"const dictionary&, "
|
||||
"const Time&"
|
||||
") const"
|
||||
)
|
||||
<< "Unhandled turbulence model option. Valid options are "
|
||||
<< "turbulenceModel, RASModel, LESModel"
|
||||
<< abort(FatalError);
|
||||
}
|
||||
}
|
||||
|
||||
Info<< " Selecting " << turbulenceType << ": " << turbulenceModel
|
||||
<< endl;
|
||||
|
||||
IOdictionary turbModelDict
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
fileName(turbModelDir/turbulenceModel),
|
||||
runTime,
|
||||
IOobject::MUST_READ
|
||||
)
|
||||
);
|
||||
|
||||
switch (solverType_)
|
||||
{
|
||||
case stIncompressible:
|
||||
{
|
||||
fieldTemplates.merge(turbModelDict.subDict("incompressibleFields"));
|
||||
break;
|
||||
}
|
||||
case stCompressible:
|
||||
case stBuoyant:
|
||||
{
|
||||
fieldTemplates.merge(turbModelDict.subDict("compressibleFields"));
|
||||
break;
|
||||
}
|
||||
default:
|
||||
{
|
||||
// do nothing
|
||||
}
|
||||
}
|
||||
|
||||
return fieldTemplates;
|
||||
}
|
||||
|
||||
|
||||
Foam::dictionary Foam::solverTemplate::readSolidFieldTemplates
|
||||
(
|
||||
const word& regionName,
|
||||
const dictionary& solverDict
|
||||
) const
|
||||
{
|
||||
Info<< " Reading solid field templates for region " << regionName
|
||||
<< endl;
|
||||
|
||||
// nothing more to do for solids
|
||||
return solverDict.subDict("solidFields");
|
||||
}
|
||||
|
||||
|
||||
void Foam::solverTemplate::setRegionProperties
|
||||
(
|
||||
const dictionary& fieldDict,
|
||||
const word& regionType,
|
||||
const word& regionName,
|
||||
const label regionI
|
||||
)
|
||||
{
|
||||
regionTypes_[regionI] = regionType,
|
||||
regionNames_[regionI] = regionName,
|
||||
fieldNames_[regionI] = fieldDict.toc();
|
||||
fieldTypes_[regionI].setSize(fieldNames_[regionI].size());
|
||||
fieldDimensions_[regionI].setSize(fieldNames_[regionI].size());
|
||||
|
||||
forAll(fieldNames_[regionI], i)
|
||||
{
|
||||
const word& fieldName = fieldNames_[regionI][i];
|
||||
const dictionary& dict = fieldDict.subDict(fieldName);
|
||||
|
||||
dict.lookup("type") >> fieldTypes_[regionI][i];
|
||||
fieldDimensions_[regionI].set
|
||||
(
|
||||
i,
|
||||
new dimensionSet(dict.lookup("dimensions"))
|
||||
);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::solverTemplate::solverTemplate
|
||||
(
|
||||
const fileName& baseDir,
|
||||
const Time& runTime,
|
||||
const word& solverName
|
||||
)
|
||||
:
|
||||
solverType_(stUnknown),
|
||||
multiRegion_(false),
|
||||
regionTypes_(),
|
||||
fieldNames_(),
|
||||
fieldTypes_(),
|
||||
fieldDimensions_()
|
||||
{
|
||||
IOdictionary solverDict
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
fileName(baseDir/"solvers"/solverName),
|
||||
runTime,
|
||||
IOobject::MUST_READ
|
||||
)
|
||||
);
|
||||
|
||||
Info<< "Selecting " << solverName << ": ";
|
||||
|
||||
solverType_ = solverTypeNames_.read(solverDict.lookup("solverType"));
|
||||
Info<< solverTypeNames_[solverType_];
|
||||
|
||||
multiRegion_ = readBool(solverDict.lookup("multiRegion"));
|
||||
if (multiRegion_)
|
||||
{
|
||||
Info<< ", multi-region";
|
||||
}
|
||||
else
|
||||
{
|
||||
Info<< ", single-region";
|
||||
}
|
||||
|
||||
Info<< " case" << endl;
|
||||
|
||||
if (multiRegion_)
|
||||
{
|
||||
// read regionProperties
|
||||
regionProperties rp(runTime);
|
||||
|
||||
const wordList fluidNames(rp["fluid"]);
|
||||
const wordList solidNames(rp["solid"]);
|
||||
|
||||
const label nRegion = fluidNames.size() + solidNames.size();
|
||||
|
||||
regionTypes_.setSize(nRegion);
|
||||
regionNames_.setSize(nRegion);
|
||||
fieldNames_.setSize(nRegion);
|
||||
fieldTypes_.setSize(nRegion);
|
||||
fieldDimensions_.setSize(nRegion);
|
||||
|
||||
// read templates for solver lists of avaliable
|
||||
// - fields and dimensions required for the solver
|
||||
// - models
|
||||
label regionI = 0;
|
||||
forAll(fluidNames, i)
|
||||
{
|
||||
const dictionary fieldDict
|
||||
(
|
||||
readFluidFieldTemplates
|
||||
(
|
||||
fluidNames[i],
|
||||
baseDir,
|
||||
solverDict,
|
||||
runTime
|
||||
)
|
||||
);
|
||||
|
||||
setRegionProperties(fieldDict, "fluid", fluidNames[i], regionI++);
|
||||
}
|
||||
|
||||
forAll(solidNames, i)
|
||||
{
|
||||
const dictionary fieldDict
|
||||
(
|
||||
readSolidFieldTemplates
|
||||
(
|
||||
solidNames[i],
|
||||
solverDict
|
||||
)
|
||||
);
|
||||
setRegionProperties(fieldDict, "solid", solidNames[i], regionI++);
|
||||
}
|
||||
}
|
||||
else
|
||||
{
|
||||
regionTypes_.setSize(1);
|
||||
regionNames_.setSize(1);
|
||||
fieldNames_.setSize(1);
|
||||
fieldTypes_.setSize(1);
|
||||
fieldDimensions_.setSize(1);
|
||||
|
||||
// read templates for solver lists of avaliable
|
||||
// - fields and dimensions required for the solver
|
||||
// - models
|
||||
const dictionary fieldDict
|
||||
(
|
||||
readFluidFieldTemplates
|
||||
(
|
||||
word::null, // use region = null for single-region cases
|
||||
baseDir,
|
||||
solverDict,
|
||||
runTime
|
||||
)
|
||||
);
|
||||
|
||||
setRegionProperties(fieldDict, "fluid", word::null, 0);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
|
||||
|
||||
Foam::word Foam::solverTemplate::type() const
|
||||
{
|
||||
return solverTypeNames_[solverType_];
|
||||
}
|
||||
|
||||
|
||||
bool Foam::solverTemplate::multiRegion() const
|
||||
{
|
||||
return multiRegion_;
|
||||
}
|
||||
|
||||
|
||||
label Foam::solverTemplate::nRegion() const
|
||||
{
|
||||
return regionTypes_.size();
|
||||
}
|
||||
|
||||
|
||||
const Foam::word& Foam::solverTemplate::regionType(const label regionI) const
|
||||
{
|
||||
return regionTypes_[regionI];
|
||||
}
|
||||
|
||||
|
||||
const Foam::word& Foam::solverTemplate::regionName(const label regionI) const
|
||||
{
|
||||
return regionNames_[regionI];
|
||||
}
|
||||
|
||||
|
||||
const Foam::wordList& Foam::solverTemplate::fieldNames
|
||||
(
|
||||
const label regionI
|
||||
) const
|
||||
{
|
||||
return fieldNames_[regionI];
|
||||
}
|
||||
|
||||
|
||||
const Foam::wordList& Foam::solverTemplate::fieldTypes
|
||||
(
|
||||
const label regionI
|
||||
) const
|
||||
{
|
||||
return fieldTypes_[regionI];
|
||||
}
|
||||
|
||||
|
||||
const Foam::PtrList<Foam::dimensionSet>& Foam::solverTemplate::fieldDimensions
|
||||
(
|
||||
const label regionI
|
||||
) const
|
||||
{
|
||||
return fieldDimensions_[regionI];
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
Reference in New Issue
Block a user