CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   capturing cell centers (https://www.cfd-online.com/Forums/openfoam-programming-development/243058-capturing-cell-centers.html)

saicharan662000@gmail.com May 27, 2022 00:31

capturing cell center to center distances
 
Hi guys,
I added a code for getting cell center to center distance. And I want only the center to center distances on x direction. But i am not getting the expected result. I am attaching the code

Code:

const unallocLabelList& neighbour = mesh.neighbour();
const unallocLabelList& owner = mesh.owner();
 const volVectorField& C = mesh.C();
forAll(owner, facei)
 {
  vector delta = C[neighbour[facei]] - C[owner[facei]];

  const dimensionedScalar len("len",dimLength,delta.component(vector::X));
 
  Info << "len: " << len.value() << " m" <<endl;
 
}

as my cell size is (0.01,0.0025,0.05) and uniform mesh I should get len = 0.01 all over the timesteps,
But i am getting worng result. I am attaching the result
Code:

len: 0.01 m
len: -8.67362e-19 m
len: 0.01 m
len: 1.73472e-18 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 6.93889e-18 m
len: 0.01 m
len: 2.08167e-17 m
len: 0.01 m
len: 6.93889e-18 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: -1.38778e-17 m
len: 0.01 m
len: 1.38778e-17 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 2.77556e-17 m
len: 0.01 m
len: -1.38778e-17 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: -2.77556e-17 m
len: 0.01 m
len: -2.77556e-17 m
len: 0.01 m
len: -2.77556e-17 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: -2.77556e-17 m
len: 0.01 m
len: 2.77556e-17 m
len: 0.01 m
len: -2.77556e-17 m
len: 0.01 m
len: -5.55112e-17 m
len: 0.01 m
len: -8.32667e-17 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 5.55112e-17 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: -5.55112e-17 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: -5.55112e-17 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 5.55112e-17 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 0 m
len: 0.01 m
len: 5.55112e-17 m
len: 0.01 m
len: -5.55112e-17 m

Can anyone tell me whats wrong ? open to accept any advice.

Thanks in advance


All times are GMT -4. The time now is 22:09.