[shipman]

Hi to All,

I am trying to plot the residual by using gnuplot. One of the iteration saved in my log file is as follows:

DILUPBiCG: Solving for alpha1, Initial residual = 4.38942e-05, Final residual = 1.79432e-11, No Iterations 1
Liquid phase volume fraction = 0.999624 Min(alpha1) = 0.0663787 Max(alpha1) = 1
PIMPLE: iteration 1
GAMG: Solving for p_rgh, Initial residual = 0.0390795, Final residual = 0.00328955, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 0.00474834, Final residual = 0.000251044, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 0.00695737, Final residual = 0.000493905, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 0.00151405, Final residual = 7.09018e-05, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 0.00310708, Final residual = 0.000183107, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 0.000543072, Final residual = 1.95473e-05, No Iterations 2
PIMPLE: iteration 2
GAMG: Solving for p_rgh, Initial residual = 0.00142025, Final residual = 6.74854e-05, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 0.000205943, Final residual = 7.01413e-06, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 0.000606703, Final residual = 2.57123e-05, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 8.15515e-05, Final residual = 3.07118e-06, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 0.000268793, Final residual = 1.01818e-05, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 3.37112e-05, Final residual = 1.29354e-06, No Iterations 2
PIMPLE: iteration 3
GAMG: Solving for p_rgh, Initial residual = 0.000123156, Final residual = 4.13567e-06, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1.46604e-05, Final residual = 5.41923e-07, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 5.54985e-05, Final residual = 1.6922e-06, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 6.40674e-06, Final residual = 2.2905e-07, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 2.58142e-05, Final residual = 7.05381e-07, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 2.96623e-06, Final residual = 1.05465e-07, No Iterations 2
PIMPLE: iteration 4
GAMG: Solving for p_rgh, Initial residual = 1.23282e-05, Final residual = 3.02816e-07, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1.46768e-06, Final residual = 5.97165e-08, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 5.82797e-06, Final residual = 1.35998e-07, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 7.50777e-07, Final residual = 3.27265e-08, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 2.83949e-06, Final residual = 6.4069e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 4.13769e-07, Final residual = 1.63602e-08, No Iterations 2
PIMPLE: iteration 5
GAMG: Solving for p_rgh, Initial residual = 1.43059e-06, Final residual = 3.14307e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 2.44785e-07, Final residual = 8.03719e-09, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 7.19056e-07, Final residual = 1.59301e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1.48655e-07, Final residual = 1.32698e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 3.65799e-07, Final residual = 1.3061e-08, No Iterations 1
GAMGPCG: Solving for p_rgh, Initial residual = 9.6084e-08, Final residual = 9.6084e-08, No Iterations 0
DILUPBiCG: Solving for omega, Initial residual = 1.05021e-07, Final residual = 2.1834e-12, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 4.41147e-06, Final residual = 6.37142e-12, No Iterations 1
ExecutionTime = 59.58 s ClockTime = 60 s

And i set the my Residual file as follows:

plot "< cat log | grep 'Solving for alpha1' | cut -d' ' -f9 | tr -d ','" title 'alpha1' with lines,\
     "< cat log | grep 'Solving for omega' | cut -d' ' -f9 | tr -d ','" title 'omega' with lines,\
     "< cat log | grep 'Solving for k' | cut -d' ' -f9 | tr -d ','" title 'k' with lines,\
 "< cat log | grep 'Solving for Uk' | cut -d' ' -f9 | tr -d ','" title 'U' with lines,\
     "< cat log | grep 'Solving for p_rgh' | cut -d' ' -f9 | tr -d ','" title 'p_rgh' with lines
pause 1
reread

But the problem is in the log file i cant see solving of the U. I am using interphasechangeFoam. Could someone tell me that how can i see the solving of the U in my .log file.

Thanks in advance.

[ahmmedshakil]

Hi,
In my cases, I prefer to use pyFoamPlotWatcher.py (which comes with pyFoam: http://openfoamwiki.net/index.php/Contrib/PyFoam). It's easy to use, and no need of any extra programming.

cheers,
#shakil

[shipman]

Hi Shakil,

Thank you for your email.

Actually, I already used pyFoamPlotWatcher.py. But as you can see in my one iteration,

GAMG: Solving for p_rgh, Initial residual = 1.43059e-06, Final residual = 3.14307e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 2.44785e-07, Final residual = 8.03719e-09, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 7.19056e-07, Final residual = 1.59301e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1.48655e-07, Final residual = 1.32698e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 3.65799e-07, Final residual = 1.3061e-08, No Iterations 1
GAMGPCG: Solving for p_rgh, Initial residual = 9.6084e-08, Final residual = 9.6084e-08, No Iterations 0
DILUPBiCG: Solving for omega, Initial residual = 1.05021e-07, Final residual = 2.1834e-12, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 4.41147e-06, Final residual = 6.37142e-12, No Iterations 1
ExecutionTime = 59.58 s ClockTime = 60 s

I cant see the solving of the U like this:

smoothSolver: [B]Solving for Ux[/B], Initial residual = 0.001088847, Final residual = 8.176141e-05, No Iterations 6 
smoothSolver: [B]Solving for Uy[/B], Initial residual = 0.003432152, Final residual = 0.0002506418, No Iterations 6 
smoothSolver: [B]Solving for Uz[/B], Initial residual = 0.004479452, Final residual = 0.0003201872, No Iterations 6 
GAMG: Solving for p, Initial residual = 0.02337899, Final residual = 0.0008043688, No Iterations 2 
GAMG: Solving for p, Initial residual = 0.007815368, Final residual = 0.0005850772, No Iterations 1 
GAMG: Solving for p, Initial residual = 0.001613976, Final residual = 0.0001103631, No Iterations 2 
GAMG: Solving for p, Initial residual = 0.001157861, Final residual = 8.520689e-05, No Iterations 1 
GAMG: Solving for p, Initial residual = 0.0002718626, Final residual = 2.288669e-05, No Iterations 2 
time step continuity errors : sum local = 1.079882e-05, global = 3.552218e-08, cumulative = 0.002264786 
smoothSolver: Solving for epsilon, Initial residual = 0.0006579817, Final residual = 5.450659e-05, No

So, my equation is that how can i see the velocity solving inside my iteration?

Thank you

[wyldckat]

Greetings to all!

@Baris: The first problem here is that your initial description is seriously misleading The problem isn't a plotting issue, the problem is that the residuals for the velocity simply don't exist in the log file.

The first approach on this is to search (instead of asking ) on the forum if someone else has seen this issue in the past. Let's see... I searched for:

interPhaseChangeFoam velocity residuals

and the very first post that appeared on Google for me was this: http://www.cfd-online.com/Forums/ope...hangefoam.html
And the answer is right there on the second post:

Originally Posted by abe

The main reason could be that the velocity equation UEqn is not solved.
You may change the momentumPredictor condition in System/fvSolution file--PIMPLE entity from 'no' to 'yes'. So, the velocity equations will be solved and you can see their results.

Best regards,
Bruno

[shipman]

Dear Bruno,

Thank you so much for your answer. Yeah you are right. sorry for my misleading.I should be more careful before asking.

B.