|
[Sponsors] |
June 26, 2013, 13:54 |
are residulas real or normalized?
|
#1 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
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?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. Last edited by immortality; June 26, 2013 at 15:18. |
|
June 26, 2013, 14:35 |
|
#2 |
Senior Member
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 16 |
Dear Ehsan
Can you tell me how to see the residuals in ParaFoam? (btw I don't know the answer of your question) Thanks
__________________
“If you have an apple and I have an apple and we exchange these apples then you and I will still each have one apple. But if you have an idea and I have an idea and we exchange these ideas, then each of us will have two ideas.” |
|
June 26, 2013, 15:21 |
|
#3 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
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?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
|
June 27, 2013, 03:20 |
|
#4 |
Senior Member
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 16 |
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.
__________________
“If you have an apple and I have an apple and we exchange these apples then you and I will still each have one apple. But if you have an idea and I have an idea and we exchange these ideas, then each of us will have two ideas.” |
|
June 27, 2013, 06:41 |
|
#5 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
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.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
|
January 1, 2014, 12:21 |
|
#6 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
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? Code:
U { solver smoothSolver; smoother DICGaussSeidel; //nSweeps 2; tolerance 1e-14; relTol 0; maxIter 25000; }
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
|
January 1, 2014, 13:09 |
|
#7 | |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Hi Ehsan,
Quote:
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
__________________
|
||
January 2, 2014, 02:29 |
|
#8 |
Senior Member
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 21 |
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.
|
|
January 2, 2014, 06:01 |
|
#9 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
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.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
error: uninitialized local variable 't' used | MASOUD | Fluent UDF and Scheme Programming | 5 | October 17, 2016 04:24 |
Problems in air flow udf - divergence | PJT | Fluent UDF and Scheme Programming | 0 | May 28, 2013 10:01 |
view factor | yoollooz | Fluent UDF and Scheme Programming | 0 | March 1, 2013 00:44 |
defining a term for a domain using DEFINE_ADJUST | MASOUD | Fluent UDF and Scheme Programming | 1 | September 24, 2010 05:08 |
enum | MASOUD | Fluent UDF and Scheme Programming | 0 | June 5, 2010 00:49 |