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
|