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

ENH: When overridding cell centres, change face centre to midpoint.

Browse Source
This commit is contained in:
graham
2010-12-02 18:04:03 +00:00
parent 7335b82d52
commit aa26d6703f
1 changed files with 25 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

25
src/OpenFOAM/meshes/primitiveMesh/primitiveMeshCellCentresAndVols.C
View File

@ -203,6 +203,31 @@ void Foam::primitiveMesh::overrideCellCentres
calcCellCentresAndVols();
*cellCentresPtr_ = newCellCtrs;
// Set internal face centres to the midpoint of the cell-centre delta vector
if (debug)
{
Pout<< "void Foam::primitiveMesh::overrideCellCentres"
<< "(const vectorField& newCellCtrs) const : "
<< "overriding internal face centres." << endl;
}
deleteDemandDrivenData(faceCentresPtr_);
deleteDemandDrivenData(faceAreasPtr_);
calcFaceCentresAndAreas();
vectorField& fCtrs = *faceCentresPtr_;
const vectorField& C = cellCentres();
const labelUList& owner = faceOwner();
const labelUList& neighbour = faceNeighbour();
forAll(neighbour, faceI)
{
fCtrs[faceI] = 0.5*(C[neighbour[faceI]] + C[owner[faceI]]);
}
}