[Martin_K_lalelu]

Hi all,

I added a volume force rho*AF to UEqn in an interDyMFoam solver to handle a moving reference frame.
After running successfully in non-parallel mode, the parallel calculations of this volume force fails.

AF and alpha1 are volVectorFields, x is a vector.

x=gSum(alpha1.internalField()*
            mesh.C().internalField()*mesh.V())/
            gSum(alpha1.internalField()*mesh.V());

AF=alpha1*f(x)

Is there a specific global formulation for parallelized calculations?

The attachement shows the AF value over iterations. Up to iteration 170 in parallel mode (wrong solution), after that starting from the same initial condition in non-parallelized mode (correct solution) .


Any advice will be appreciated!

best regards

Martin

Edit: Solved, AF was calculated correctly.