[CFDanielGER]

Hello everybody,


for a new hybrid turbulence model I want to implement I need the dimensions of each cell to calculate a new delta. I found a way of calculating the dimensions in each direction (x, y, z) in meter. See code attached:

  const faceList & ff = this->mesh().faces();
  const pointField & pp = this->mesh().points();
  volScalarField deltax = this->y_;
  volScalarField deltay = this->y_;
  volScalarField deltaz = this->y_;

  forAll(this->mesh().C(), celli)
  {
    const cell & cc = this->mesh().cells()[celli];
    labelList pLabels(cc.labels(ff));
    pointField pLocal(pLabels.size(), vector::zero);

    forAll(pLabels, pointi)
    {
      pLocal[pointi] = pp[pLabels[pointi]];
    }

    dimensionedScalar Convert ("Convert", dimensionSet(0, 1, 0, 0, 0, 0, 0), scalar(1));

    deltax = (Foam::max(pLocal & vector(1,0,0)) - Foam::min(pLocal & vector(1,0,0))) * Convert;
    deltay = (Foam::max(pLocal & vector(0,1,0)) - Foam::min(pLocal & vector(0,1,0))) * Convert;
    deltaz = (Foam::max(pLocal & vector(0,0,1)) - Foam::min(pLocal & vector(0,0,1))) * Convert;
  }

I think the code is working fine, at least it is compiling and a few iterations are running but I have to check further.


But my question is:


Is there a more effective way of computing the cell dimensions in x, y and z than using a forAll loop? I have the impression that this is rather time-consuming when using a large mesh.



Many thanks for your help!


Daniel