[florian_krause]

Dear all,
I use a modified simpleFoam in OpenFOAM-1.7.x. I simulate a laminar pipe flow with uniform fixed value inlet boundary conditons and zeroGradient outlet boundary conditions. I've set up the case as a wedge using wedge boundary conditions. so far so good...

Now I want to stop the simulations if a pre-defined convergence criterion is satisfied for pressure p and the axial (x) and wall-normal (y) velocity component.

If I use

tmp<fvVectorMatrix> UEqn
(
fvm::div(phi, U)
+ turbulence->divDevReff(U)
);

UEqn().relax();

eqnResidualU = solve
(
UEqn() == -fvc::grad(p)
).initialResidual();

maxResidualU = max(eqnResidualU, maxResidualU);

then I get the initial residuals for all three velocity components. And using

if (maxResidualU < convergenceCriterionU && maxResidualP < convergenceCriterionP)
{
Info<< "reached convergence criterion for U and p"<< endl;
runTime.writeAndEnd();
Info<< "latestTime = " << runTime.timeName() << endl;
}

I perform the actual convergence check on all three velocity components.

As mentioned above, for this wedge case, I want to stop the calculation if only the x and y component of U satisfy the pre-defined convergence criterion.

Does anyone know how to do that? As always, any help is appreciated!

Best,
Florian