[AbhishekTAMU]

I want to know if it is possible to get the value of cell skewness (please see the code attached from cellQuality.C) when the simulation itself is running. skewness is shown as an output to 'checkMesh' command and generally it is done after the simulation has been run.

However, I would like to know if I can get the value of skewness when the simulation itself is running and I would also like to know if this value of skewness can be stored as a 'volScalar Field'.


#include "cellQuality.H"
#include "unitConversion.H"
#include "SubField.H"

//This is "cellQuality.C"


Foam::tmp<Foam::scalarField> Foam::cellQuality::skewness() const
{
tmp<scalarField> tresult
(
new scalarField
(
mesh_.nCells(), 0.0
)
);
scalarField& result = tresult.ref();

scalarField sumArea(mesh_.nCells(), 0.0);

const vectorField& cellCtrs = mesh_.cellCentres();
const vectorField& faceCtrs = mesh_.faceCentres();
const vectorField& areas = mesh_.faceAreas();

const labelList& own = mesh_.faceOwner();
const labelList& nei = mesh_.faceNeighbour();

forAll(nei, facei)
{
scalar dOwn = mag
(
(faceCtrs[facei] - cellCtrs[own[facei]]) & areas[facei]
)/mag(areas[facei]);

scalar dNei = mag
(
(cellCtrs[nei[facei]] - faceCtrs[facei]) & areas[facei]
)/mag(areas[facei]);

point faceIntersection =
cellCtrs[own[facei]]
+ (dOwn/(dOwn+dNei))*(cellCtrs[nei[facei]] - cellCtrs[own[facei]]);

scalar skewness =
mag(faceCtrs[facei] - faceIntersection)
/(mag(cellCtrs[nei[facei]] - cellCtrs[own[facei]]) + vSmall);

result[own[facei]] = max(skewness, result[own[facei]]);

result[nei[facei]] = max(skewness, result[nei[facei]]);
}

forAll(mesh_.boundaryMesh(), patchi)
{
const labelUList& faceCells =
mesh_.boundaryMesh()[patchi].faceCells();

const vectorField::subField faceCentres =
mesh_.boundaryMesh()[patchi].faceCentres();

const vectorField::subField faceAreas =
mesh_.boundaryMesh()[patchi].faceAreas();

forAll(faceCentres, facei)
{
vector n = faceAreas[facei]/mag(faceAreas[facei]);

point faceIntersection =
cellCtrs[faceCells[facei]]
+ ((faceCentres[facei] - cellCtrs[faceCells[facei]])&n)*n;

scalar skewness =
mag(faceCentres[facei] - faceIntersection)
/(
mag(faceCentres[facei] - cellCtrs[faceCells[facei]])
+ vSmall
);

result[faceCells[facei]] = max(skewness, result[faceCells[facei]]);
}
}

return tresult;
}