[titio]

Hi foamers,

I am starting to feel desesperate guys. I am trying to calculate the mass flow through a plane that it is not a patch, as the massCalcFlow does. Is there any way of this the easy way, say using foamCalc or some other utility that I am not aware of. If not possible to do that easily, I will do it in the code, and I will contribute it to the community, if I manage to do it.

Can anyone help me out. Thanks in advance,

António Martins

PS: OpenFOAM really needs a good documentation, but it seems capitalism and common good do not mix very well... and the problems with the documentation site only helped the commercial software even more....

[olesen]

Originally Posted by titio

Hi foamers,

I am starting to feel desesperate guys. I am trying to calculate the mass flow through a plane that it is not a patch, as the massCalcFlow does. Is there any way of this the easy way, say using foamCalc or some other utility that I am not aware of. If not possible to do that easily, I will do it in the code, and I will contribute it to the community, if I manage to do it.

Can anyone help me out. Thanks in advance,

I'd guess that there just isn't a generic place for such a massflow sampling and that's why it hasn't made its way into the code.
Note that this is all much easier when sampling on a patch, since you can simply use the phi directly there.

Nonetheless, this code fragment should give you a few ideas (I hope).

// defined somewhere in the surrounding class:
PtrList<sampledPlane> gammaPlanes_;


void Foam::massFlowInfo::calculateAndWrite
(
    const volVectorField& U,
    const volScalarField& rho
) const
{
    forAll(gammaPlanes_, planeI)
    {
        const sampledPlane& pln = gammaPlanes_[planeI];

        label nFaces = pln.faces().size();
        reduce(nFaces, sumOp<label>());

        if (nFaces)
        {
            // the velocity normal to the plane
            scalarField Uproj = pln.project
            (
                pln.sample(U)
            );

            scalar U_mean = pln.average(Uproj);

            scalarField rhoU = pln.sample(rho) * Uproj;
            scalar massFlow = pln.integrate(rhoU);

            Info<< pln.name() << nl
                << "  area       = " << pln.area() << nl
                << "  massFlow   = " << massFlow << "  ("
                << (massFlow * 3600) << " kg/h)" << nl
                << "  U_mean     = " << U_mean
                << endl;
        }
        else
        {
            Info<< pln.name() << nl
                << "  area       = " << 0
                << endl;
        }
    }
}

[titio]

Thanks olesen. I will try the code and I will let you know if its worked.

Thanks,

Titio

[chegdan]

Was this resolved? if so, please share how you made this possible.

[robingilbert]

try swak4Foam:

http://openfoamwiki.net/index.php/Contrib/swak4Foam

to calculate flow rate through internal faceZone this works. i have tried it. to calculate the flow through a patch u can use calcMassFlow:

http://openfoamwiki.net/index.php/Contrib_calcMassFlow

[maysmech]

hi,


How can i calculate mass flow rate of a patch by using paraView software?

[misakagan]

Originally Posted by olesen

I'd guess that there just isn't a generic place for such a massflow sampling and that's why it hasn't made its way into the code.
Note that this is all much easier when sampling on a patch, since you can simply use the phi directly there.

Nonetheless, this code fragment should give you a few ideas (I hope).

// defined somewhere in the surrounding class:
PtrList<sampledPlane> gammaPlanes_;


void Foam::massFlowInfo::calculateAndWrite
(
    const volVectorField& U,
    const volScalarField& rho
) const
{
    forAll(gammaPlanes_, planeI)
    {
        const sampledPlane& pln = gammaPlanes_[planeI];

        label nFaces = pln.faces().size();
        reduce(nFaces, sumOp<label>());

        if (nFaces)
        {
            // the velocity normal to the plane
            scalarField Uproj = pln.project
            (
                pln.sample(U)
            );

            scalar U_mean = pln.average(Uproj);

            scalarField rhoU = pln.sample(rho) * Uproj;
            scalar massFlow = pln.integrate(rhoU);

            Info<< pln.name() << nl
                << "  area       = " << pln.area() << nl
                << "  massFlow   = " << massFlow << "  ("
                << (massFlow * 3600) << " kg/h)" << nl
                << "  U_mean     = " << U_mean
                << endl;
        }
        else
        {
            Info<< pln.name() << nl
                << "  area       = " << 0
                << endl;
        }
    }
}

Dear Mr. Olesen,

The code is very helpful, but I wonder also how the parts related to construction of sampledPlane's is implemented. Normally sampledPlane's are constructed on runtime using the sampleDict with runtimeselectiontable. I don't know how to construct them in a regular way with the direction vectors and basepoints.

Any suggestion would be a lot of help.

Best,

Asim

[misakagan]

Ok I found the answer. I have to create cuttingPlane's and construct sampledPlane's using those cuttingPlane's...

[krndv]

Originally Posted by maysmech

hi,


How can i calculate mass flow rate of a patch by using paraView software?

hello,

have you got answer for this? Please share if you know/

Thank you

[Flowkersma]

Originally Posted by krndv

hello,

have you got answer for this? Please share if you know/

Thank you

1. Calculate rho * U_N with "Calculator".
2. Use "IntegrateVariables" filter to calculate the massflow.