[Cyp]

hello!

I am facing a new problem : I do not know how to evaluate the surface of a patchField labeled patchID.

I was thinking about

p.boundaryField[patchID].magSf();

but it doesn't work.


Thank you very much for your help,
Cyp

[MartinB]

Hi Cyp,

try:

p.boundaryField()[patchID].magSf();

Martin

[Cyp]

Hi Martin!

It doesn't work neither...

(in fact it was what I tried... my first post was a typo)

[MartinB]

May be

mesh.magSf().boundaryField()[patchID];

?

Can you post some more code?

Martin

[Cyp]

It still doesn't work...

This is the piece of code I consider:

    surfaceScalarField phiKsa = fvc::interpolate(kSa);
    surfaceScalarField phiKsb = fvc::interpolate(kSb);

    surfaceScalarField phiP
    (
        IOobject
        (
            "phiP",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::NO_WRITE
        ),
        mesh,
        dimensionedScalar("zero", dimensionSet(1, -2, -2, 0, 0), 0.0)
    );

    phiP = ((-fvc::interpolate(Ub+Ua)&mesh.Sf())+((rhoa*phiKsa+rhob*phiKsb)*g&mesh.Sf()))
            / (phiKsa + phiKsb)/mesh.magSf().boundaryField()[patchID];

[MartinB]

phiP.boundaryField() = ((-fvc::interpolate(Ub+Ua)&mesh.Sf())+((rhoa*phiKsa+rhob*phiKsb)*g&mesh.Sf()))/(phiKsa + phiKsb)/mesh.magSf().boundaryField();

?

What's the error message? Is it one of the standard solvers of OpenFOAM, so I can check it?

May be you can loop over the patches like here:

volVectorField tau
        (
            IOobject
            (
                "tau",
                runTime.timeName(),
                mesh,
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            mesh,
            dimensionedVector
            (
                "tau",
                Reff.dimensions(),
                vector::zero
            )
        );

        forAll(tau.boundaryField(), patchi)
        {
            tau.boundaryField()[patchi] =
            (
                -mesh.Sf().boundaryField()[patchi]
                /mesh.magSf().boundaryField()[patchi]
            ) & Reff.boundaryField()[patchi];
        }

[Cyp]

thank you very much for your help. It works now!

    scalarField fluxP = -(phi.boundaryField()[patchID])
                /(kSa.boundaryField()[patchID]+kSb.boundaryField()[patchID])
                /mesh.magSf().boundaryField()[patchID];