[immortality]

the residual of equations that we can see also by pyFoam(thanks to dear Bernhard Gschaider) are real or normalized?if the latter by what?

[cfdonline2mohsen]

Dear Ehsan
Can you tell me how to see the residuals in ParaFoam?
(btw I don't know the answer of your question)
Thanks

[immortality]

Hi Mohsen
I corrected it.I meant pyFoam
you also can see residuals after run by gnuplot if you had not it be saved by pyFoam.there is a script code for that but pyFoam way is better(and easier) in my opinion.
do you know how to work with pyFoam?

[cfdonline2mohsen]

Thanks for the reply.
Yeah, I've used gnuplot but it works (as far as I know and used) only for a single processor run but I couldn't utilize it for my parallel run with multiple processors. do you how to use it in parallel?
And no, I am not familiar with pyFoam but I will search the forums for it.
thanks for the hint.

[immortality]

see openfoamwiki and search for pyFoam contrib (as I remember)
after installing I'll tell you which commands you use.it shows residuals during run and also save it to case folder and also save residuals for next uses.its very useful.

[immortality]

Hi,
I have noticed and encountered another question about residuals,does anyone know?
for example I have set the tolerance for U to be 1e-14, but in residuals that are shown,its far more of 1e-14.its around 1e-5.why?

U
    {
        solver          smoothSolver;
        smoother        DICGaussSeidel;
        //nSweeps         2;
        tolerance       1e-14;
        relTol          0;
        maxIter        25000;
    }

[wyldckat]

Hi Ehsan,

Originally Posted by immortality

for example I have set the tolerance for U to be 1e-14, but in residuals that are shown,its far more of 1e-14.its around 1e-5.why?

Did you look at the output from the solver?

The idea is that the residual that matters to us for considering if it the solution is as good as "converged", is the Initial residual, which should be lower than 1e-2 to 1e-4 or smaller, depending on the kind of simulation you're performing. That is the residual value being controled inside the "SIMPLE", "PISO" and "PIMPLE" settings in "system/fvSolution" (not certain if all of them have got it). And that is the one being plotted.

The residual being controlled in the matrix/equation solver is the Final residual, which is also present in the output of the solver (Initial ... Final ... No iterations...).

Best regards,
Bruno

[Bernhard]

Also, I suppose you are hitting the maxIter limit? You can easily check that with your log-file. In any case, especially a transient one, you do not want 25000 iterations, because this is mostly related to a bad case. Also, you are showing the solver settings for U. Which means that you are using the momentum predictor, I guess. Why would you want to solve it down to this accuracy? I don't think it makes much sense.

[immortality]

Oh Bruno,I thought that I was knowing that before but had forgotten about initial and final residuals.I remembered all things now. thanks for help.
@Bernhard:thanks to your consideration dear Bernhard,I used that maxIter when I was encountering convergence problems that was because of boundary conditions,after the problems be solved,I didn't managed to change it,thought it may be better to be there.I wanted reaching to convergence be sure.