[bernarde]

Hi

Thanks for this excellent script.

Is it possible to zoom into a section of the residual plot?

[ata]

Hi every one
When I try to plot residuals with gnuplot I got this error:
plot "< cat log | grep 'Solving for Ux' | cut -d' ' -f9 | tr -d ','" title 'Ux' with lines,"< cat log | grep 'Solving for Uy' | cut -d' ' -f9 | tr -d ','" title 'Uy' with lines,"< cat log | grep 'Solving for Uz' | cut -d' ' -f9 | tr -d ','" title 'Uz' 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 p' | cut -d' ' -f9 | tr -d ','" title 'p' with lines
^
"Residuals", line 10: warning: Skipping data file with no valid points

plot "< cat log | grep 'Solving for Ux' | cut -d' ' -f9 | tr -d ','" title 'Ux' with lines,"< cat log | grep 'Solving for Uy' | cut -d' ' -f9 | tr -d ','" title 'Uy' with lines,"< cat log | grep 'Solving for Uz' | cut -d' ' -f9 | tr -d ','" title 'Uz' 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 p' | cut -d' ' -f9 | tr -d ','" title 'p' with lines
^
"Residuals", line 10: warning: Skipping data file with no valid points

plot "< cat log | grep 'Solving for Ux' | cut -d' ' -f9 | tr -d ','" title 'Ux' with lines,"< cat log | grep 'Solving for Uy' | cut -d' ' -f9 | tr -d ','" title 'Uy' with lines,"< cat log | grep 'Solving for Uz' | cut -d' ' -f9 | tr -d ','" title 'Uz' 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 p' | cut -d' ' -f9 | tr -d ','" title 'p' with lines
^
"Residuals", line 10: warning: Skipping data file with no valid points

plot "< cat log | grep 'Solving for Ux' | cut -d' ' -f9 | tr -d ','" title 'Ux' with lines,"< cat log | grep 'Solving for Uy' | cut -d' ' -f9 | tr -d ','" title 'Uy' with lines,"< cat log | grep 'Solving for Uz' | cut -d' ' -f9 | tr -d ','" title 'Uz' 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 p' | cut -d' ' -f9 | tr -d ','" title 'p' with lines
Can any body help me?
Regards

Ata

[ata]

Hi every one
When I try to plot residuals with gnuplot I got this error:
plot "< cat log | grep 'Solving for Ux' | cut -d' ' -f9 | tr -d ','" title 'Ux' with lines,"< cat log | grep 'Solving for Uy' | cut -d' ' -f9 | tr -d ','" title 'Uy' with lines,"< cat log | grep 'Solving for Uz' | cut -d' ' -f9 | tr -d ','" title 'Uz' 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 p' | cut -d' ' -f9 | tr -d ','" title 'p' with lines
^
"Residuals", line 10: warning: Skipping data file with no valid points

plot "< cat log | grep 'Solving for Ux' | cut -d' ' -f9 | tr -d ','" title 'Ux' with lines,"< cat log | grep 'Solving for Uy' | cut -d' ' -f9 | tr -d ','" title 'Uy' with lines,"< cat log | grep 'Solving for Uz' | cut -d' ' -f9 | tr -d ','" title 'Uz' 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 p' | cut -d' ' -f9 | tr -d ','" title 'p' with lines
^
"Residuals", line 10: warning: Skipping data file with no valid points

plot "< cat log | grep 'Solving for Ux' | cut -d' ' -f9 | tr -d ','" title 'Ux' with lines,"< cat log | grep 'Solving for Uy' | cut -d' ' -f9 | tr -d ','" title 'Uy' with lines,"< cat log | grep 'Solving for Uz' | cut -d' ' -f9 | tr -d ','" title 'Uz' 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 p' | cut -d' ' -f9 | tr -d ','" title 'p' with lines
^
"Residuals", line 10: warning: Skipping data file with no valid points

plot "< cat log | grep 'Solving for Ux' | cut -d' ' -f9 | tr -d ','" title 'Ux' with lines,"< cat log | grep 'Solving for Uy' | cut -d' ' -f9 | tr -d ','" title 'Uy' with lines,"< cat log | grep 'Solving for Uz' | cut -d' ' -f9 | tr -d ','" title 'Uz' 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 p' | cut -d' ' -f9 | tr -d ','" title 'p' with lines
Can any body help me?
Regards

Ata

[RBJ]

It just means you are trying to extract fields from the 'log' file, which don't exist. May be you don't have a turbulence model and try to extract k and eps in example. Change your residual.txt file accordingly.

[ata]

Hello RBj
Thanks. It seems so.
Best Regards

Ata

[vinz]

Dear Openfoamers,

Wolfgang tool is very useful.
However, I have a question. I usually start my computation with 3 nonOrthogonalCorector loops and then decreased to 0.
if I use the script of wolfgang with the following line:
"< cat log | grep 'Solving for p' | cut -d' ' -f9 | sed -n 'p;N;N' | tr -d ','" title 'p' with lines

the line of p residual just stop on the graph when I decrease the number of p loops.
how would it be possible to handle that?

regards,

Vincent

[lth]

Dear Wolfgang,

want to add to the list, thank you for this script!

[sarahk]

Thanks for the code, its great ...

can i use it as well to plot the time step continuity errors?
tried to exchange it, but i always get an error ...

thanks in advance ...

[alberto]

You might want to take a look at:

http://openfoamwiki.net/index.php/Contrib_PyFoam

[anger]

Originally Posted by sarahk

Thanks for the code, its great ...

can i use it as well to plot the time step continuity errors?
tried to exchange it, but i always get an error ...

thanks in advance ...

sarahk,

What's your script and what error do you get?

Regards,
-Thomas

[mvoss]

hi,

i´m using the gnuplot approach quite a while and it is really helpful.

But how about this: i want to know where the residuals are high and not convergent and want to look at the residual-fields (u,v,w, final p e.g.) in paraview.
Basically it´s easy, just create field for each residualvalue and store them BUT i do not know how to manage smth. like residuals (u,v,w,p).write.


Any ideas?

neewbie

[afo3]

hello,
I've followed thew instruction of the tutorial to plot the forces.dat file with the forces of my simulation, generated by simpleFoam.

here is my problem: gnuplot, only plots some of the columns, for example 1:3 (time versus Fy) but not 1:2 (time versus Fx) it gives me the following:

Warning : skipping data file with no valid points

attached there's the file.


any ideas??

thanks!

Michele

ps thanks to Wolfgang for the thread!

[afo3]

I think the problem is due to the parentesis in my file..there's a way to convert the .dat file in order to be read properly by gnuplot?

EDIT:

solved here:

http://www.cfd-online.com/Forums/ope...-brackets.html

[mvoss]

hi,

good to know about the previous post. I solved it with
plot "..." using ... '%lf ( %lf %lf %lf) ( %lf %lf %lf)' but this only works up to 10 x %lf ...so max. 3 Vectors + first row.

neewbie

[cfdonline2mohsen]

hi
thank you

[Amir]

Hi guys,
thank you Wolfgang !

[megacrout]

Awesome.
Thx Wolfgang.
Tibo

[calim_cfd]

it's an old topic but still i feel like i need to thank, it surely speeds things up

[suh]

Thanks for the script,

One additional information is, the sequence of residuals reading in command line for gnuplot also to be same as solved by that particular solver. Hope it helps someone.

ThankYou.
Regards
Suhas

[robingilbert]

Hi everyone,

I tried this:

set logscale y
set title "Residuals"
set ylabel 'Residual'
set xlabel 'Iteration'
plot "< cat log | grep 'Solving for Ux' | cut -d' ' -f9 | sed -n 'Ux;N' | tr -d ','" title 'Ux' with lines,\
     "< cat log | grep 'Solving for Uy' | cut -d' ' -f9 | sed -n 'Uy;N' | tr -d ','" title 'Uy' with lines,\
     "< cat log | grep 'Solving for Uz' | cut -d' ' -f9 | sed -n 'Uz;N' | tr -d ','" title 'Uz' with lines,\
     "< cat log | grep 'Solving for k' | cut -d' ' -f9 | sed -n 'k;N'| tr -d ','" title 'k' with lines,\
     "< cat log | grep 'Solving for p' | cut -d' ' -f9 | sed -n 'p;N;N;N' | tr -d ','" title 'p' with lines 
pause 1
reread

But this plots only p but gives errors for U and k. Can anybody help me?