zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Adding a pointField write to the writeTargetCellSize function to allow

Browse Source 
comparision of the isoSurfaces created from the interpolated cell
centre values on the polyhedral dual mesh and the directly calculated
point values on the tetMesh (i.e. the raw Delaunay).
This commit is contained in:
graham
2009-10-21 11:30:16 +01:00
parent 481852be9a
commit 8b4b531357
1 changed files with 116 additions and 70 deletions
Download Patch File Download Diff File Expand all files Collapse all files

186
src/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshIO.C
View File

@ -28,6 +28,8 @@ License
#include "IOstreams.H"
#include "OFstream.H"
#include "zeroGradientPointPatchField.H"
#include "pointMesh.H"
#include "pointFields.H"
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
@ -134,39 +136,39 @@ void Foam::conformalVoronoiMesh::writeMesh(bool writeToConstant)
{
writeInternalDelaunayVertices(writeToConstant);
// {
// pointField points;
// faceList faces;
// labelList owner;
// labelList neighbour;
// wordList patchNames;
// labelList patchSizes;
// labelList patchStarts;
{
pointField points;
faceList faces;
labelList owner;
labelList neighbour;
wordList patchNames;
labelList patchSizes;
labelList patchStarts;
// calcTetMesh
// (
// points,
// faces,
// owner,
// neighbour,
// patchNames,
// patchSizes,
// patchStarts
// );
calcTetMesh
(
points,
faces,
owner,
neighbour,
patchNames,
patchSizes,
patchStarts
);
// writeMesh
// (
// "tetDualMesh",
// writeToConstant,
// points,
// faces,
// owner,
// neighbour,
// patchNames,
// patchSizes,
// patchStarts
// );
// }
writeMesh
(
"tetDualMesh",
writeToConstant,
points,
faces,
owner,
neighbour,
patchNames,
patchSizes,
patchStarts
);
}
{
pointField points;
@ -317,48 +319,92 @@ void Foam::conformalVoronoiMesh::writeObjMesh
void Foam::conformalVoronoiMesh::writeTargetCellSize() const
{
Info<< nl << "Create fvMesh" << endl;
fvMesh fMesh
(
IOobject
(
Foam::polyMesh::defaultRegion,
runTime_.constant(),
runTime_,
IOobject::MUST_READ
)
);
timeCheck();
Info<< nl << "Create targetCellSize volScalarField" << endl;
volScalarField targetCellSize
(
IOobject
(
"targetCellSize",
runTime_.timeName(),
runTime_,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
fMesh,
dimensionedScalar("cellSize", dimLength, 0),
zeroGradientPointPatchField<scalar>::typeName
);
scalarField& cellSize = targetCellSize.internalField();
const vectorField& C = fMesh.cellCentres();
forAll(cellSize, i)
{
cellSize[i] = cellSizeControl().cellSize(C[i]);
Info<< nl << "Create fvMesh" << endl;
fvMesh fMesh
(
IOobject
(
Foam::polyMesh::defaultRegion,
runTime_.constant(),
runTime_,
IOobject::MUST_READ
)
);
timeCheck();
Info<< nl << "Create targetCellSize volScalarField" << endl;
volScalarField targetCellSize
(
IOobject
(
"targetCellSize",
runTime_.timeName(),
runTime_,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
fMesh,
dimensionedScalar("cellSize", dimLength, 0),
zeroGradientPointPatchField<scalar>::typeName
);
scalarField& cellSize = targetCellSize.internalField();
const vectorField& C = fMesh.cellCentres();
forAll(cellSize, i)
{
cellSize[i] = cellSizeControl().cellSize(C[i]);
}
targetCellSize.write();
}
{
polyMesh tetMesh
(
IOobject
(
"tetDualMesh",
runTime_.constant(),
runTime_,
IOobject::MUST_READ
)
);
pointMesh ptMesh(tetMesh);
pointScalarField ptTargetCellSize
(
IOobject
(
"ptTargetCellSize",
runTime_.timeName(),
tetMesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
ptMesh,
dimensionedScalar("ptTargetCellSize", dimLength, 0),
pointPatchVectorField::calculatedType()
);
scalarField& cellSize = ptTargetCellSize.internalField();
const vectorField& P = tetMesh.points();
forAll(cellSize, i)
{
cellSize[i] = cellSizeControl().cellSize(P[i]);
}
ptTargetCellSize.write();
}
targetCellSize.write();
}