[GRAUPS]

Hey guys, I'm looking to calculate a mass weighted average over a faceZone for pressure using swak4foam. Specifically, this is the equation I'm trying to implement...

Capture.PNG

My question comes when trying to implement. Is it possible in swak4foam to divide a sum by a sum within one expression? Can I do two sums separately and divide them? So far I've only seen expressions that can do one accumulation, such as the one below for mass flow rate...

mass-in
    {
        type swakExpression;
        valueType faceZone;
        zoneName interior-inlet-1;
        expression "phi*flip()";
        accumulations (
            sum
        );
        verbose true;
	autoInterpolate	true;
    }

Also, (if possible) what is the syntax for getting the absolute value and dot product within a swakExpression? Thanks in advance!

[GRAUPS]

Ok, I think I found an equivalent expression. The following function I believe does what I need it to do...

p_stat2
    {
        type swakExpression;
        valueType faceZone;
        zoneName interior-inlet-1;
        expression "(p*phi*flip())/sum(phi*flip())";
        accumulations (
            sum
        );
        verbose true;
	autoInterpolate	true;
    }

However, I would still be interested in knowing if you can perform a dot product or an absolute value within a swakExpression?

Also, does anyone know if when you have backflow on a face, does it register as negative mass flow using the phi*flip() formulation?

[gschaider]

Originally Posted by GRAUPS

Ok, I think I found an equivalent expression. The following function I believe does what I need it to do...

p_stat2
    {
        type swakExpression;
        valueType faceZone;
        zoneName interior-inlet-1;
        expression "(p*phi*flip())/sum(phi*flip())";
        accumulations (
            sum
        );
        verbose true;
	autoInterpolate	true;
    }

However, I would still be interested in knowing if you can perform a dot product or an absolute value within a swakExpression?

Also, does anyone know if when you have backflow on a face, does it register as negative mass flow using the phi*flip() formulation?

Dot product is according to the conventions of OpenFOAM the &-operator.

Yes. You should get the correct backflow if you use flip() (the sign of phi depends on the orientation of the face. Multiplying with flip() "removes" that dependence)

[GRAUPS]

Originally Posted by gschaider

Dot product is according to the conventions of OpenFOAM the &-operator.

Awesome, thank you. And using the mag() function on a scalar I found is basically the equivalent of taking the absolute value.

Originally Posted by gschaider

Yes. You should get the correct backflow if you use flip() (the sign of phi depends on the orientation of the face. Multiplying with flip() "removes" that dependence)

Fantastic, just wanted to make sure. Thanks for your replies.

[Alex2]

Dear Foamers,

I want to use groovybc, in order to create a T bc for a reciculation device. To do this I have do calculate the mass-flux-weighted temperature of a face. I use the following code:

PHP Code:

    Auslasslinks_1_5    {        type            groovyBC;        valueExpression "Tein";    variables    "massenstrom=1;Taus{Ansaug_3}=sum(T*phi)/sum(phi);cp=1007;Qdot=100;Tein=Taus-Qdot/(massenstrom*cp);";    fractionExpression "1";        value           uniform 300;        refValue        $value;    } 


Unfurtunately this bc doesn't work and my simulation crashes because of the Taus{Ansaug_3} calculation. If I calculate a faceweighted temperature instead of the mass flux weighted temperature for Taus{Ansaug_3} or if I use Taus instead of Taus{Ansaug_3}, the simulation runs. So my question: Is this a bug in groovybc or is there a bug in my code?

Thank's in advance

[TRIDEV G]

Hi guys
I am trying to calculate mass average weighted kinetic energy for Tank sloshing case.
The equation is similar to the one discussed in this thread.

(summation(mass*0.5*U*U)/ summation (mass))

Can someone help me with the syntax to write this expression using swak4Foam ?

Iam new to openFoam. any help will be highly apreciated.

Thank you.

[gschaider]

Originally Posted by TRIDEV G

Hi guys
I am trying to calculate mass average weighted kinetic energy for Tank sloshing case.
The equation is similar to the one discussed in this thread.

(summation(mass*0.5*U*U)/ summation (mass))

Can someone help me with the syntax to write this expression using swak4Foam ?

Iam new to openFoam. any help will be highly apreciated.

Thank you.

The expression for the total mass would be (assuming you have the denisty in a field rho) "sum(vol()*rho)". The sum of the kinetic energy would be "sum(vol()*rho*0.5*(U&U))" (using & for the scalar product). Now with swakExpression you could write the complete expression as "(vol()*rho*0.5*(U&U))/sum(vol()*rho)" (without the sum) and in the accumulations-list use "sum". Or you write the equation with the sum and in the accumulations use one of "min", "max", "average". But the first option needs a little less computation

[TRIDEV G]

Originally Posted by gschaider

The expression for the total mass would be (assuming you have the denisty in a field rho) "sum(vol()*rho)". The sum of the kinetic energy would be "sum(vol()*rho*0.5*(U&U))" (using & for the scalar product). Now with swakExpression you could write the complete expression as "(vol()*rho*0.5*(U&U))/sum(vol()*rho)" (without the sum) and in the accumulations-list use "sum". Or you write the equation with the sum and in the accumulations use one of "min", "max", "average". But the first option needs a little less computation

Hello gschaider, Thank you very much for the reply

The equation is writing results perfectly.

But I dont understand the logic behind using summation in the numerator for calculating max, min, average and not using summation for calculating sum (accumulation).

the equation calculates mass average weighted kinetic energy for each cell in the domain for each timestep and calculates sum. where is the difference coming by introducing the summation in the numerator ?

There is something that I am missing out in understanding here.
Can u please explain this??

Thanks in advance

[anthony761]

Hello Foamers,

I'm new in OF, I'm studying compressible flows and for post-processing, I just want to have confirmation. I'm not sure about syntax, but I would like to calculate the mass-weighted average of one variable, defines as at my "r1-outblock-outflow" (patch name)' so is it correct to write in my

controlDict (OpenFOAM v2006)

    weightAvgOut
    {
        #includeEtc "caseDicts/postProcessing/surfaceFieldValue/surfaceRegion.cfg"
        regionType      patch;
        name            "r1-outblock-outflow";
        operation       weightedAverage;
        weightField     phi;
        fields          (p T Ma total(p));
        writeInterval   5;
        log             true;
    }

Do you think this is right?

Also, does anyone know if OF can provide the total temperature?