[yhy20081016]

Hello everybody,
I am recently working on modifying the compressible solver "rhoCentralFoam".
Here I want to retrieve the boundary field of "e" and "p", then substitute them into a user-defined function named "ep_d", then assigned the result to the boundary field "rho".

...
p.correctBoundaryConditions();
forAll(rho.boundaryField(),i)
{
rho.boundaryField()[i] = pp.ep_d(e.boundaryField()[i], p.boundaryField()[i]);

}
runTime.write();
...

However I got a compilation error:


nistCentralFoam.C: In function ‘int main(int, char**)’:
nistCentralFoam.C:248:85: error: no matching function for call to ‘phys_property::ep_d(Foam::fvPatchField<double>&, Foam::fvPatchField<double>&)’
rho.boundaryField()[i] = pp.ep_d(e.boundaryField()[i], p.boundaryField()[i]);
^
In file included from nistCentralFoam.C:46:0:
property.H:248:6: note: candidate: double phys_property::ep_d(double, double)
REAL phys_property::ep_d(REAL e, REAL p)
^
property.H:248:6: note: no known conversion for argument 1 from ‘Foam::fvPatchField<double>’ to ‘double’

So what is the solution? Thank you very much.

[adhiraj]

Will this work:

forAll(rho.boundaryField(), patchI)
{
    fvPatchScalarField& rhoPatch = rho.boundaryField()[patchI];
    fvPatchScalarField& ePatch = e.boundaryField()[patchI];
    fvPatchScalarField& pPatch = p.boundaryField()[patchI];

    forAll(pPatch, faceI)
    {
        const scalar ei = ePatch[faceI];
        const scalar pi = pPatch[faceI];
        rhoPatch[faceI] = pp.ep_d(ei, pi);
    }
}

[yhy20081016]

Originally Posted by adhiraj

Will this work:

forAll(rho.boundaryField(), patchI)
{
    fvPatchScalarField& rhoPatch = rho.boundaryField()[patchI];
    fvPatchScalarField& ePatch = e.boundaryField()[patchI];
    fvPatchScalarField& pPatch = p.boundaryField()[patchI];

    forAll(pPatch, faceI)
    {
        const scalar ei = ePatch[faceI];
        const scalar pi = pPatch[faceI];
        rhoPatch[faceI] = pp.ep_d(ei, pi);
    }
}

Yes, it works. Thank you.