[mkhm]

Hi foamers,

My question is simple: U for each cell has 3 components for a 3D case and U is calculated for the whole domain, so many cell. Imagine that I want to find out the maximum value of the difference of U(iteration2) and U(iteration1) over the domain.


- for openfoam, which one is correct ?
1)

max(mag(U_iteration1)-mag(U_iteration2))

2)

max(mag(U_iteration1-U_iteration2))

Does the following notations are used by OpenFoam ?
U_cell(i)=[Ux Uy Uz]_cell(i):
U_iteration1=[ [Ux Uy Uz]_cell(1)
[Ux Uy Uz]_cell(2)
.

.
.
[Ux Uy Uz]_cell(n)];


mag(U_iteration1)=[ mag([Ux Uy Uz]_cell(1))
mag([Ux Uy Uz]_cell(2))
.

.
.
mag([Ux Uy Uz]_cell(n))];

Thanks for your help,
Mary

[mkhm]

In other terms, could we replace maxValU[/CODE] simply by:

max(mag(U_f)-mag(U_0))

?

 forAll(U_f,i)
        {
            DiffRhoU =mag(mag(U_f[i])-mag(U_0[i]));
            if(DiffU>maxValU)
            {
                maxValU=DiffU;
                posMaxU=i;
            }
            sumU_f  += mag(U_f[i]);
            sumDiffU += DiffU;
         }