zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

mirrorMesh: Added support for cellLevel and pointLevel to support dynamic mesh refinement

Browse Source 
Patch contributed by Mattijs Janssens
Resolves bug-report https://bugs.openfoam.org/view.php?id=2712
This commit is contained in:
Henry Weller
2017-10-27 14:58:26 +01:00
parent 3ad3db0c86
commit 19da21ed93
6 changed files with 166 additions and 38 deletions
Download Patch File Download Diff File Expand all files Collapse all files

72
applications/utilities/mesh/manipulation/mirrorMesh/mirrorMesh.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -32,6 +32,8 @@ Description
#include "argList.H"
#include "Time.H"
#include "mirrorFvMesh.H"
#include "mapPolyMesh.H"
#include "hexRef8Data.H"
using namespace Foam;
@ -39,18 +41,12 @@ using namespace Foam;
int main(int argc, char *argv[])
{
#include "addOverwriteOption.H"
#include "setRootCase.H"
#include "createTime.H"
const bool overwrite = args.optionFound("overwrite");
if (!overwrite)
{
runTime++;
}
mirrorFvMesh mesh
(
IOobject
@ -61,11 +57,71 @@ int main(int argc, char *argv[])
)
);
hexRef8Data refData
(
IOobject
(
"dummy",
mesh.facesInstance(),
polyMesh::meshSubDir,
mesh,
IOobject::READ_IF_PRESENT,
IOobject::NO_WRITE,
false
)
);
if (!overwrite)
{
runTime++;
mesh.setInstance(runTime.timeName());
}
// Set the precision of the points data to 10
IOstream::defaultPrecision(max(10u, IOstream::defaultPrecision()));
// Generate the mirrorred mesh
const fvMesh& mMesh = mesh.mirrorMesh();
const_cast<fvMesh&>(mMesh).setInstance(mesh.facesInstance());
Info<< "Writing mirrored mesh" << endl;
mesh.mirrorMesh().write();
mMesh.write();
// Map the hexRef8 data
mapPolyMesh map
(
mesh,
mesh.nPoints(), //nOldPoints,
mesh.nFaces(), //nOldFaces,
mesh.nCells(), //nOldCells,
mesh.pointMap(), //pointMap,
List<objectMap>(0), // pointsFromPoints,
labelList(0), //faceMap,
List<objectMap>(0), //facesFromPoints,
List<objectMap>(0), //facesFromEdges,
List<objectMap>(0), //facesFromFaces,
mesh.cellMap(), //cellMap,
List<objectMap>(0), //cellsFromPoints,
List<objectMap>(0), //cellsFromEdges,
List<objectMap>(0), //cellsFromFaces,
List<objectMap>(0), //cellsFromCells,
labelList(0), //reversePointMap,
labelList(0), //reverseFaceMap,
labelList(0), //reverseCellMap,
labelHashSet(0), //flipFaceFlux,
labelListList(0), //patchPointMap,
labelListList(0), //pointZoneMap,
labelListList(0), //faceZonePointMap,
labelListList(0), //faceZoneFaceMap,
labelListList(0), //cellZoneMap,
pointField(0), //preMotionPoints,
labelList(0), //oldPatchStarts,
labelList(0), //oldPatchNMeshPoints,
autoPtr<scalarField>() //oldCellVolumesPtr
);
refData.updateMesh(map);
refData.write();
Info<< "End" << endl;