First version of dynamicPolyRefinementFvMesh
Refactored version of dynamicRefinePolyFvMesh in terms of topoChangerFvMesh
This commit is contained in:
Vuko Vukcevic
parent
15fc87e5dd
commit
ae02d50179
2 changed files with 296 additions and 0 deletions
|
|
@ -0,0 +1,178 @@
|
||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | foam-extend: Open Source CFD
|
||||||
|
\\ / O peration | Version: 4.0
|
||||||
|
\\ / A nd | Web: http://www.foam-extend.org
|
||||||
|
\\/ M anipulation | For copyright notice see file Copyright
|
||||||
|
-------------------------------------------------------------------------------
|
||||||
|
License
|
||||||
|
This file is part of foam-extend.
|
||||||
|
|
||||||
|
foam-extend 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.
|
||||||
|
|
||||||
|
foam-extend 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 foam-extend. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
#include "dynamicPolyRefinementFvMesh.H"
|
||||||
|
#include "addToRunTimeSelectionTable.H"
|
||||||
|
#include "refinementSelection.H"
|
||||||
|
#include "polyhedralRefinement.H"
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
namespace Foam
|
||||||
|
{
|
||||||
|
|
||||||
|
defineTypeNameAndDebug(dynamicPolyRefinementFvMesh, 0);
|
||||||
|
|
||||||
|
addToRunTimeSelectionTable
|
||||||
|
(
|
||||||
|
topoChangerFvMesh,
|
||||||
|
dynamicPolyRefinementFvMesh,
|
||||||
|
IOobject
|
||||||
|
);
|
||||||
|
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||||
|
|
||||||
|
void Foam::dynamicPolyRefinementFvMesh::readDict()
|
||||||
|
{
|
||||||
|
// Read and check refinement interval
|
||||||
|
refineInterval_ = readLabel(refinementDict_.lookup("refineInterval"));
|
||||||
|
if (refineInterval_ < 1)
|
||||||
|
{
|
||||||
|
FatalErrorIn("dynamicPolyRefinementFvMesh::readDict()")
|
||||||
|
<< "Illegal refineInterval found: " << refineInterval << nl
|
||||||
|
<< "The refineInterval controls the refinement/unrefinement"
|
||||||
|
<< " trigerring within a certain time step and should be > 0".
|
||||||
|
<< exit(FatalError);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
dynamicPolyRefinementFvMesh::dynamicPolyRefinementFvMesh(const IOobject& io)
|
||||||
|
:
|
||||||
|
dynamicFvMesh(io),
|
||||||
|
refinementDict_
|
||||||
|
(
|
||||||
|
IOdictionary
|
||||||
|
(
|
||||||
|
IOobject
|
||||||
|
(
|
||||||
|
"dynamicMeshDict",
|
||||||
|
time().constant(),
|
||||||
|
*this,
|
||||||
|
IOobject::MUST_READ,
|
||||||
|
IOobject::NO_WRITE,
|
||||||
|
false
|
||||||
|
)
|
||||||
|
).subDict(typeName + "Coeffs")
|
||||||
|
),
|
||||||
|
refineInterval_(readLabel(refinementDict_.lookup("refineInterval"))),
|
||||||
|
refinementSelectionPtr_(refinementSelection::New(*this, refinementDict_))
|
||||||
|
{
|
||||||
|
// Add the topology modifier engine
|
||||||
|
Info<< "Adding polyhedralRefinement topology modifier" << endl;
|
||||||
|
|
||||||
|
topoChanger_.setSize(1);
|
||||||
|
topoChanger_.set
|
||||||
|
(
|
||||||
|
0,
|
||||||
|
new polyhedralRefinement
|
||||||
|
(
|
||||||
|
"polyhedralRefinement",
|
||||||
|
refinementDict_,
|
||||||
|
0,
|
||||||
|
topoChanger_
|
||||||
|
)
|
||||||
|
);
|
||||||
|
|
||||||
|
// Write mesh and modifiers
|
||||||
|
topoChanger_.writeOpt() = IOobject::AUTO_WRITE;
|
||||||
|
topoChanger_.write();
|
||||||
|
write();
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
dynamicPolyRefinementFvMesh::~dynamicPolyRefinementFvMesh()
|
||||||
|
{}
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
bool dynamicPolyRefinementFvMesh::update()
|
||||||
|
{
|
||||||
|
// Re-read the data from dictionary for on-the-fly changes
|
||||||
|
readDict();
|
||||||
|
|
||||||
|
// Performing refinement/unrefinement when:
|
||||||
|
// 1. We are at the first time step
|
||||||
|
// 2. Time step is a multiplier of specified refineInterval
|
||||||
|
if
|
||||||
|
(
|
||||||
|
time().timeIndex() > 0
|
||||||
|
&& time().timeIndex() % refineInterval == 0
|
||||||
|
)
|
||||||
|
{
|
||||||
|
// Get reference to polyhedralRefinement polyMeshModifier
|
||||||
|
polyhedralRefinement& polyRefModifier =
|
||||||
|
refCast<polyhedralRefinement>(topoChanger_[0]);
|
||||||
|
|
||||||
|
// Get refinement candidates from refinement selection algorithm
|
||||||
|
labelList cellsToRefine
|
||||||
|
(
|
||||||
|
refinementSelectionPtr_->cellsToRefine()
|
||||||
|
);
|
||||||
|
|
||||||
|
// Set cells to refine.
|
||||||
|
// Notes:
|
||||||
|
// 1. Invalidates the argument list
|
||||||
|
// 2. polyhedralRefinement ensures that face and point consistent
|
||||||
|
// refinement is performed
|
||||||
|
polyRefModifier.setCellsToRefine(cellsToRefine);
|
||||||
|
|
||||||
|
// Get unrefinement split points from refinement selection algorithm
|
||||||
|
labelList splitPointsToUnrefine
|
||||||
|
(
|
||||||
|
refinementSelectionPtr_->splitPointsToUnrefine()
|
||||||
|
);
|
||||||
|
|
||||||
|
// Set split points to unrefine around.
|
||||||
|
// Notes:
|
||||||
|
// 1. Invalidates the argument list
|
||||||
|
// 2. polyhedralRefinement ensures that only a consistent set of
|
||||||
|
// split points is used for unrefinement
|
||||||
|
polyRefModifier.setSplitPointsToUnrefine(splitPointsToUnrefine);
|
||||||
|
|
||||||
|
// Perform refinement and unrefinement in one go
|
||||||
|
autoPtr<mapPolyMesh> topoChangeMap = topoChanger_.changeMesh();
|
||||||
|
|
||||||
|
return topoChangeMap->morphing();
|
||||||
|
}
|
||||||
|
|
||||||
|
// No refinement/unrefinement at this time step. Return false
|
||||||
|
return false;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
} // End namespace Foam
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
|
|
@ -0,0 +1,118 @@
|
||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | foam-extend: Open Source CFD
|
||||||
|
\\ / O peration | Version: 4.0
|
||||||
|
\\ / A nd | Web: http://www.foam-extend.org
|
||||||
|
\\/ M anipulation | For copyright notice see file Copyright
|
||||||
|
-------------------------------------------------------------------------------
|
||||||
|
License
|
||||||
|
This file is part of foam-extend.
|
||||||
|
|
||||||
|
foam-extend 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.
|
||||||
|
|
||||||
|
foam-extend 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 foam-extend. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
|
||||||
|
Class
|
||||||
|
Foam::dynamicPolyRefinementFvMesh
|
||||||
|
|
||||||
|
Description
|
||||||
|
Adaptive mesh refinement for arbitrary polyhedral cells.
|
||||||
|
|
||||||
|
SourceFiles
|
||||||
|
dynamicPolyRefinementFvMesh.C
|
||||||
|
|
||||||
|
Author
|
||||||
|
Vuko Vukcevic, Wikki Ltd. All rights reserved
|
||||||
|
|
||||||
|
Notes
|
||||||
|
Generalisation and refactorisation of dynamicRefineMesh for polyhedral cells
|
||||||
|
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
#ifndef dynamicPolyRefinementFvMesh_H
|
||||||
|
#define dynamicPolyRefinementFvMesh_H
|
||||||
|
|
||||||
|
#include "topoChangerFvMesh.H"
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
namespace Foam
|
||||||
|
{
|
||||||
|
|
||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
Class dynamicPolyRefinementFvMesh Declaration
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
class dynamicPolyRefinementFvMesh
|
||||||
|
:
|
||||||
|
public topoChangerFvMesh
|
||||||
|
{
|
||||||
|
// Private data
|
||||||
|
|
||||||
|
//- Dictionary containing dynamic mesh coefficient controls
|
||||||
|
dictionary refinementDict_;
|
||||||
|
|
||||||
|
//- Refinement interval
|
||||||
|
label refineInterval_;
|
||||||
|
|
||||||
|
//- Refinement selection algorithm that has two jobs:
|
||||||
|
// 1. Selects cells to refine from all cells
|
||||||
|
// 2. Selects split points that are ok to refine based on given set of
|
||||||
|
// all split points
|
||||||
|
autoPtr<refinementSelection> refinementSelectionPtr_;
|
||||||
|
|
||||||
|
|
||||||
|
// Private Member Functions
|
||||||
|
|
||||||
|
//- Helper function for reading the dictionary and updating the data
|
||||||
|
void readDict();
|
||||||
|
|
||||||
|
// Copy control
|
||||||
|
|
||||||
|
//- Disallow default bitwise copy construct
|
||||||
|
dynamicPolyRefinementFvMesh(const dynamicPolyRefinementFvMesh&);
|
||||||
|
|
||||||
|
//- Disallow default bitwise assignment
|
||||||
|
void operator=(const dynamicPolyRefinementFvMesh&);
|
||||||
|
|
||||||
|
public:
|
||||||
|
|
||||||
|
//- Runtime type information
|
||||||
|
TypeName("dynamicPolyRefinementFvMesh");
|
||||||
|
|
||||||
|
|
||||||
|
// Constructors
|
||||||
|
|
||||||
|
//- Construct from IOobject
|
||||||
|
explicit dynamicPolyRefinementFvMesh(const IOobject& io);
|
||||||
|
|
||||||
|
|
||||||
|
//- Destructor
|
||||||
|
virtual ~dynamicPolyRefinementFvMesh();
|
||||||
|
|
||||||
|
|
||||||
|
// Member Functions
|
||||||
|
|
||||||
|
//- Update the mesh for topology change
|
||||||
|
virtual bool update();
|
||||||
|
};
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
} // End namespace Foam
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
#endif
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
Reference in a new issue