[linch]

Hi,

I have a question to following lines in src/finiteVolume/cfdTools/incompressible/continuityErrs.H

volScalarField contErr = fvc::div(phi);
scalar globalContErr = runTime.deltaT().value()*
        contErr.weightedAverage(mesh.V()).value();

globalContErr should be an integral of the error over the complete domain volume and time step. fvc::div(phi) is a volume integral for each cell. Thus, time integration and summation over all cells is needed to obtain the global error. We already have the time integration *runTime.deltaT(), but why do we have averaging instead of summation over all cells?

Best regards,
Ilya

[linch]

No ideas?

[boger]

I don't have any answers for you, except to add my own observation that the localContErr would, I think, be equivalent to the L1-norm of the pressure equation residual were it not for the volume-averaging.