[chenye1717]

Dear all,

I am new to OpenFOAM and I have met a problem of accessing the scalarField element values in my code. Please help to take a look.

I used a 3D fluent mesh file generated by ICEM and concerted to OpenFOAM by 'Fluent3DMeshToFoam'. The simulation works well but I want to also output the surface element area at the inlet. So I tried to modify ./0/U using 'codedFixedValue'.


According to some online resources, I have done these:
const fvPatch& boundaryPatch = patch();
const fvBoundaryMesh& boundaryMesh = boundaryPatch.boundaryMesh();
const fvMesh& mesh = boundaryMesh.mesh();
label patchID = mesh.boundaryMesh().findPatchID("INLET_R");
const surfaceScalarField& magSf = mesh.magSf();
const scalarField& magSf0 = magSf.boundaryField()[patchID];

As you may see, I have found the patchID of the particular inlet I am looking at and the value of patchID is correct. My question is how can I output the values in either magSf or magSf0 one by one?

I have tried to output the sum of both magSf and magSf0 by 'gSum(magSf.boundaryField()[patchID])' and 'gSum(magSf0)'. The values are both correct.

I tried to use:

forAll(magSf.boundaryField()[patchID],ii)
{
Info << "ii "<<ii<<" Area "<<magSf.boundaryField()[patchID][ii]<<endl;
}

But it never outputs. Could you guys help to find out the reason?

Thanks, and best regards.

[chenye1717]

Hi all,

This issue has been fixed. The problem is caused by the inconsistent mesh information between fluent mesh and the initial simulation file of ./constant/polyMesh/boundary.

Thanks!