CFD Online Discussion Forums - Access volume field in the boundary condition of a point field

Dear Foamers

I have developed a new boundary condition for a point field. In this boundary condition, I need to calculate forces that act on a specific surface; therefore, I have loaded pressure from the object registry. The code was compiled successfully but whenever I run a case with this boundary and the code wants to access to pressure field (please see below), I receive the floating-point error.

here is the code:

Code:

        const polyMesh&  pointMeshIF = this->dimensionedInternalField().mesh()();

        

        const volScalarField& p_ = pointMeshIF.thisDb().lookupObject<volScalarField>("p");



        const fvMesh& mesh = p_.mesh();



        label IDts = mesh.boundaryMesh().findPatchID(targetSurface_);

  

         scalarField targetPatchPressureField  = p_.boundaryField()[IDts];

I could not find the actual problem, but I am sure that floating-point error is because of loading volume fields (pressure). I think it should be possible to use a volume Field in the calculation of a point field.

I will appreciate it if anyone has any idea or suggestion about this problem.

Kind regards,
Peyman