[jqtranus]

Hello,

I'm a beginner at C++, and coding in general (I'm a first year), so I'm struggling to figure this out. I've been attempting to add the contact angle model from multiphaseInterFoam into icoReactingMultiphaseInterFoam.

I'm trying to get the wall velocity by subtracting the surface velocity from the internal field velocity at the boundary.

// Calculated the component of the velocity parallel to the wall
                 vectorField Uwall
                 (
                     U().boundaryField()[patchi].patchInternalField()
                   - U().boundaryField()[patchi]
                 );
                 Uwall -= (AfHatPatch & Uwall)*AfHatPatch;

Where U() is this function:

Foam::tmp<Foam::volVectorField> Foam::phaseSystem::U() const
{
    auto tstf = tmp<volVectorField>::New
    (
        IOobject
        (
            "U",
            mesh_.time().timeName(),
            mesh_
        ),
        mesh_,
        dimensionedVector(dimVelocity, Zero)
    );

    auto& stf = tstf.ref();

    forAllConstIters(phaseModels_, iter)
    {
        stf += iter()() * iter()->U();
    }

    return tstf;
}

My problem is that I get this error:

phaseSystem/phaseSystem.C:1189:26: error: ‘class Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> >’ has no member named ‘boundaryField’
                      U().boundaryField()[patchi].patchInternalField()

From what I can understand, I believe U() is returning the internal velocity field, so why can't I get the boundary field from it?

Thanks in advance.