[saba_saeb]

Hi,

I need to take the average of a field on a boundary during the simulation. I am doing that using the normal approach:

scalar tmp= average(frictionvelocity);

where frictionvelocity is a scalar field. Everything is running like a charm in serial, but the simulation would stop if I run it in parallel. I tried to tackle the problem using this approach:

scalar tmp= average(frictionVelocity);
scalar tauw = returnReduce(tmp, sumOp<scalar>());

But the output is not correct since the calculated averages on every processors are being added together which make tauw that seems incorrect. Any idea would be appreciated.

Tnx,
Saba

[akidess]

Something along these lines might work:

int numprocs;
MPI_Comm_size(MPI_COMM_WORLD, &numprocs);

tauw /= numprocs;

Info<< tauw << endl;

Ninja edit - this is better:

label numprocs=Pstream::nProcs();
tauw /= numprocs;

[saba_saeb]

Originally Posted by akidess

Something along these lines might work:

int numprocs;
MPI_Comm_size(MPI_COMM_WORLD, &numprocs);

tauw /= numprocs;

Info<< tauw << endl;

Ninja edit - this is better:

label numprocs=Pstream::nProcs();
tauw /= numprocs;

Hi Anton,

Tnx for your reply.
That was the first idea that came to my mind, not however division by the total number of processors, but, division by no. of processors in x direction time no. of processors in z direction since friction velocity is just being calculated on the boundary located at the bottom of a channel, in my case. However, this is not still giving me the expected results since friction velocity is not being spread uniformly along the wall!
If anything seems unclear to me, please let me know.

Cheers,
Saba

[alexeym]

There is gAverage function.

Explanation and examples can be found for example at - http://www.cfd-online.com/Forums/ope...tml#post202817.

[saba_saeb]

Originally Posted by alexeym

There is gAverage function.

Explanation and examples can be found for example at - http://www.cfd-online.com/Forums/ope...tml#post202817.

Hi Alexey,

gAverage is apparently solving the problem. Tnx for your hint.

Cheers,
Saba