[BSengupta]

Hi!


I wish to create a solver that will give me coordinate of cell center and boundary face center as a solution file. I am presenting the file which I wish to get:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     4.0                                |
|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       vectorField;
    location    "0.5";
    object      coordinate;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 0 0 0 0 0];

internalField   nonuniform List<vector> 
400
(
    (interCell0_x interCell0_y interCell0_z)
    (interCell1_x interCell1_y interCell1_z)
    (interCell2_x interCell2_y interCell2_z)
    (interCell3_x interCell3_y interCell3_z)
    ...........
    ...........
    ...........
    (interCell399_x interCell399_y interCell399_z)
)
;

boundaryField
{
    movingWall
    {
        type            fixedValue;
        value           nonuniform List<vector>
        19
        (
        (face0Center0_x face0Center0_y face0Center0_z)
        (face0Center1_x face0Center1_y face0Center1_z)
        (face0Center2_x face0Center2_y face0Center2_z)
        ...........
        ...........
        ...........
        (face0Center18_x face0Center18_y face0Center18_z)
        )
    }
    fixedWalls
    {
        type            fixedValue;
        value           nonuniform List<vector>
        57
        (
        (face1Center0_x face1Center0_y face1Center0_z)
        (face1Center1_x face1Center1_y face1Center1_z)
        (face1Center2_x face1Center2_y face1Center2_z)
        ...........
        ...........
        ...........
        (face1Center56_x face1Center56_y face1Center56_z)
        )
    }
}


// ************************************************************************* //

I tried something by declaring

vectorField cellCenter
(
    IOobject
    (
        "cellCenter",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    mesh
);

for the cell center and have no idea how to get face center coordinate in the generic manner. I am aware we can get these coordinates by :

label patchI = mesh.boundaryMesh().findPatchID("movingWall");

    forAll(mesh.boundaryMesh()[patchI].faceCentres(), faceI)
    {
        scalar x = mesh.boundaryMesh()[patchI].faceCentres()[faceI].x();
        scalar y = mesh.boundaryMesh()[patchI].faceCentres()[faceI].y();
        scalar z = mesh.boundaryMesh()[patchI].faceCentres()[faceI].z();
        Info<<faceI<<" "<<x<<" "<<y<<" "<<z<<" "<<endl;
    }

something like this. But I don't want this way. I wish to have the coordinates as a result file as indicated in the first snippet.


I will be grateful to receive suggestion to achieve what I wish.


Thank you,
Bidesh

[jurich]

Hello Bidesh,

You can use

postProcess -func writeCellCentres

to get that. You can see the code to do it in src/functionObjects/field/writeCellCentres

[BSengupta]

Originally Posted by jurich

Hello Bidesh,

You can use

postProcess -func writeCellCentres

to get that. You can see the code to do it in src/functionObjects/field/writeCellCentres

Thank you. It was very useful.