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

Writing out of a pointField along with the mesh - allows restarts. Added FCC and BCC initial conditions (commented), their duals are rhombic dodecahedra and bitruncated octahedra respectively.

Browse Source
This commit is contained in:
graham
2009-02-12 19:06:22 +00:00
parent 20349d2eaf
commit 10dfba5795
1 changed files with 174 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

176
applications/utilities/mesh/generation/CV3DMesher/CV3D.C
View File

@ -400,7 +400,7 @@ void Foam::CV3D::insertPoints(const fileName& pointFileName)
) )
); );
insertPoints(points, 0.5*controls_.minCellSize2); insertPoints(points, 0.25*controls_.minCellSize2);
} }
@ -427,6 +427,8 @@ void Foam::CV3D::insertGrid()
vector delta(xR/ni, yR/nj, zR/nk); vector delta(xR/ni, yR/nj, zR/nk);
// SC lattice, also a simple cubic lattice
delta *= pow((1.0),-(1.0/3.0)); delta *= pow((1.0),-(1.0/3.0));
Random rndGen(1321); Random rndGen(1321);
@ -457,12 +459,167 @@ void Foam::CV3D::insertGrid()
if (qSurf_.wellInside(p, 0.5*controls_.minCellSize2)) if (qSurf_.wellInside(p, 0.5*controls_.minCellSize2))
{ {
initialPoints.push_back(Point(p.x(), p.y(), p.z())); initialPoints.push_back(Point(p.x(), p.y(), p.z()));
// insert(Point(p.x(), p.y(), p.z()))->index() = nVert++;
} }
} }
} }
} }
// BCC lattice, the Voronoi diagram of which is a lattice of bitruncated
// octahedron
// delta *= pow((1.0/2.0),-(1.0/3.0));
// Random rndGen(1321);
// scalar pert = controls_.randomPerturbation*cmptMin(delta);
// std::vector<Point> initialPoints;
// for (int i=0; i<ni; i++)
// {
// for (int j=0; j<nj; j++)
// {
// for (int k=0; k<nk; k++)
// {
// point p;
// p = point
// (
// x0 + i*delta.x(),
// y0 + j*delta.y(),
// z0 + k*delta.z()
// );
// if (controls_.randomiseInitialGrid)
// {
// p.x() += pert*(rndGen.scalar01() - 0.5);
// p.y() += pert*(rndGen.scalar01() - 0.5);
// p.z() += pert*(rndGen.scalar01() - 0.5);
// }
// if (qSurf_.wellInside(p, 0.5*controls_.minCellSize2))
// {
// initialPoints.push_back(Point(p.x(), p.y(), p.z()));
// }
// p = point
// (
// x0 + (i + 0.5)*delta.x(),
// y0 + (j + 0.5)*delta.y(),
// z0 + (k + 0.5)*delta.z()
// );
// if (controls_.randomiseInitialGrid)
// {
// p.x() += pert*(rndGen.scalar01() - 0.5);
// p.y() += pert*(rndGen.scalar01() - 0.5);
// p.z() += pert*(rndGen.scalar01() - 0.5);
// }
// if (qSurf_.wellInside(p, 0.5*controls_.minCellSize2))
// {
// initialPoints.push_back(Point(p.x(), p.y(), p.z()));
// }
// }
// }
// }
// FCC lattice, the Voronoi diagram of which is a lattice of rhombic
// dodecahedra
// delta *= pow((1.0/4.0),-(1.0/3.0));
// Random rndGen(1321);
// scalar pert = controls_.randomPerturbation*cmptMin(delta);
// std::vector<Point> initialPoints;
// for (int i=0; i<ni; i++)
// {
// for (int j=0; j<nj; j++)
// {
// for (int k=0; k<nk; k++)
// {
// point p;
// p = point
// (
// x0 + i*delta.x(),
// y0 + j*delta.y(),
// z0 + k*delta.z()
// );
// if (controls_.randomiseInitialGrid)
// {
// p.x() += pert*(rndGen.scalar01() - 0.5);
// p.y() += pert*(rndGen.scalar01() - 0.5);
// p.z() += pert*(rndGen.scalar01() - 0.5);
// }
// if (qSurf_.wellInside(p, 0.5*controls_.minCellSize2))
// {
// initialPoints.push_back(Point(p.x(), p.y(), p.z()));
// }
// p = point
// (
// x0 + i*delta.x(),
// y0 + (j + 0.5)*delta.y(),
// z0 + (k + 0.5)*delta.z()
// );
// if (controls_.randomiseInitialGrid)
// {
// p.x() += pert*(rndGen.scalar01() - 0.5);
// p.y() += pert*(rndGen.scalar01() - 0.5);
// p.z() += pert*(rndGen.scalar01() - 0.5);
// }
// if (qSurf_.wellInside(p, 0.5*controls_.minCellSize2))
// {
// initialPoints.push_back(Point(p.x(), p.y(), p.z()));
// }
// p = point
// (
// x0 + (i + 0.5)*delta.x(),
// y0 + j*delta.y(),
// z0 + (k + 0.5)*delta.z()
// );
// if (controls_.randomiseInitialGrid)
// {
// p.x() += pert*(rndGen.scalar01() - 0.5);
// p.y() += pert*(rndGen.scalar01() - 0.5);
// p.z() += pert*(rndGen.scalar01() - 0.5);
// }
// if (qSurf_.wellInside(p, 0.5*controls_.minCellSize2))
// {
// initialPoints.push_back(Point(p.x(), p.y(), p.z()));
// }
// p = point
// (
// x0 + (i + 0.5)*delta.x(),
// y0 + (j + 0.5)*delta.y(),
// z0 + k*delta.z()
// );
// if (controls_.randomiseInitialGrid)
// {
// p.x() += pert*(rndGen.scalar01() - 0.5);
// p.y() += pert*(rndGen.scalar01() - 0.5);
// p.z() += pert*(rndGen.scalar01() - 0.5);
// }
// if (qSurf_.wellInside(p, 0.5*controls_.minCellSize2))
// {
// initialPoints.push_back(Point(p.x(), p.y(), p.z()));
// }
// }
// }
// }
Info<< nl << initialPoints.size() << " vertices to insert." << endl; Info<< nl << initialPoints.size() << " vertices to insert." << endl;
// using the range insert (it is faster than inserting points one by one) // using the range insert (it is faster than inserting points one by one)
@ -1178,6 +1335,21 @@ void Foam::CV3D::relaxPoints(const scalar relaxation)
if (runTime_.outputTime()) if (runTime_.outputTime())
{ {
writeMesh(true); writeMesh(true);
pointIOField internalDVs
(
IOobject
(
"internalDelaunayVertices",
runTime_.path()/runTime_.timeName(),
runTime_,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
internalDelaunayVertices
);
internalDVs.write();
} }
// Remove the entire triangulation // Remove the entire triangulation