[Shibi]

Hello to all,


I am using the solver pimpleFoam in openFoam v7 to do some calculations.


After I achieve convergence of the main fields I would like to solve another equation.


So, I place the following:

 if (pimple.finalPimpleIter())
             {
                bool switch_solver(true);
                label counter(0);
                while (switch_solver)
                {
                    #include "solveMyEquation.H" //has the solve statement inside

                    counter+=1;
                    if (counter==3)
                    {
                        switch_solver=false;
                    }

                }


            }

This works fine but I would like to end the loop by a predefined value of the residuals and not the number of cycles.
Something like:

   while (residue > myPredefinedResidue)
                {
                    #include "solveMyEquation.H"
                }

Does anyone know how to do this?


Best regards

[Shibi]

Update:


in solve statement within #include "solveMyEquation.H" I made

SymmTensor<scalar> Residue_ = solve(MyEquation).initialResidual();

Outside of the include:

counter+=1;



SymmTensor<scalar> tolerance_ (1e-05,1e-05,1e-05,1e-05,1e-05,1e-05);  
 
scalar checkConvergence_(1);

forAll(Residue_,i)
{
    checkConvergence_*=pos( tolerance_[i]-Residue_[i]);
}


if (checkConvergence_>0)
{
      switch_solver=false; // if checkConvergence_ is positive, all multiplications will result in the value 1 and the residue is smaller than the tolerance
                
}
else
{
      switch_solver=true;
}



// If solver takes more than 3 iteration to converge, leave the loop
 if ( counter==3)
 {
      end_orientation_loop=false;            

 }

It does the first pimple cycle and everything is Ok. During the 2nd cycle, I get the following error:

--> FOAM FATAL ERROR: 
index 3 out of range 0 ... 2

    From function void Foam::UList<T>::checkIndex(Foam::label) const [with T = Foam::SolverPerformance<Foam::SymmTensor<double> >; Foam::label = int]
      in file /home/pc/OpenFOAM/OpenFOAM-7/src/OpenFOAM/lnInclude/UListI.H at line 106.

Any suggestions on what to do?