why pimple foam used under relaxation in all iteration of pimple Loop ?
P { margin-bottom: 0.08in; } Hi all
the under relaxation factor is used in steady simulation for better convergencey and I think it can not be used in unsteady flow in pimpleFoam1.6ext the under relaxation is not applied at the last iteration of pimple loop therefor its OK but in pimpleFoam2.x in all iteration pimple loop the flow is under relaxed. I will be so grateful if anybody can tell me that pimplefoam2.x is true or not ? Regards |
Hello,
Neither SIMPLE (where under-relaxation is applied) nor PISO (where its not necessary) is designed particularly for steady or unsteady flows. In unsteady RANS simlations you just solve to steady state for each timestep. PIMPLE works the same as PISO, but you can start the iteration process from the beginning with the final results of the PISO-alorithm (so you can also use under-relaxation, but its not necessary). This is a very short explanation, but i hope i could help you. Best regards, Andre |
Quote:
if we used URF we average the given P and U from current ddt and previous time step and it is not true |
Mechy you are right that in a transient simulation, you should not use underrelaxation on the final iteration (using PISO/PIMPLE). And thus, OpenFOAM doesn't :) In older OpenFOAM versions, you could see an "if" clause preventing underrelaxation on the last iteration (if I remember correctly). In newer versions, GeometricFields are smart:
Code:
00889 template<class Type, template<class> class PatchField, class GeoMesh> |
Quote:
you are right and it is true about p.relax() what about UEqn.relax ? in icoFoam the UEqn does not relaxed but in pisoFoam and all iteration of pimpe loop it relaxed ? do you know its reason ? Best Regards |
Relaxing a fvMatrix is different. You boost the diagonal term for the linear solver and iterate until convergence. This should always be ok in my opinion, no matter what you do with the outer loops. If you really want to be careful, you can relax all outer iterations except the final one by using the appropriate settings in the fvSolution file (e.g. relax U, but not UFinal).
- Anton |
OK thanks
in pimpleFoam OF-1.6ext at last iteration of pimple loop the RF for UEqn is set to 1. however, I agree with you Best Regards |
Hi All,
In the OF211 tutorials, the under relaxation factors are still there in controls dictionary. So does it affect the relaxation of the outer loop (PIMPLE loops)? When nOuterCorrectors=1, the under-relaxation factor should be 1? Thank you very much for any comments. Quote:
|
Dear akidess,
So based on what you mentioned in this post, we can still use the under-relaxation for the unsteady simulations? Besides, do you know the relation between under-relaxation factors and iteration times? Actually in openfoam, there seems no such criterions to judge if the solutions are convegent or not within one time step. For example, if I set nOuterCorrectors=5, the solutions are not convengent while if I set nOuterCorrectors=50, actually the solution will not be convengent either. So how to set the iteration times and reasonable under-relaxation factors? Thank you very much. Quote:
|
Yes, you should be able to use underrelaxation even in unsteady simulations. That's pretty much the idea of the PIMPLE algorithm as far as I know.
I'm not quite sure what you mean with the relation between URF and iteration times. Assuming you mean the number of iterations you'll need to get convergence and the URF minimizing that number, that's something you can't easily predict (if possible at all). PISO was developed with the goal to get rid of the underrelaxation necessary using SIMPLE, and thus attempting to achieve optimal convergence. - Anton Quote:
|
Hi,
When I use pimpleFoam with "nOuterCorrectors 1" setting, this should be the same as pisoFoam. But I figured out, that pimple doesn't seem to recognize the relaxationFactors I set. For pimple I get the same results as for piso with all relaxation factors set to 1. Do I have to tell pimple explicitly to use them? Philipp. Edit: Ok, I solved it. To everyone who encounters this: It seems, that pimpleFoam needs a declaration that looks like: Code:
relaxationFactors Code:
relaxationFactors |
Quote:
use "(p)" instead or "(U|k|omega)" |
Allright, I better change it in my original post! Thanks a lot.
... for some reason there is no "EDIT" button in that post, so I can't. Be carefull. |
As I understood: Piso=no underrelaxation,
Pimple= underrelaxation except for final iteration. Sometimes one does not need underrelaxation at all. Does underrelaxation speed up pimple or just make it more stable? Which values are usually used fur the underrelaxation and what kind of field is "nuSgs"? Do I need this one for incompressible flow with k-omega-modell? Best regards Tobi |
Hi Tobias,
The first two things are correct. As far as I know, underrelaxation is always used to make an algorithm more stable. Sometimes it can happen that you become faster using it. But that is just by pure chance, because you "accidentally" hit the right value. The last point: I wondered about that, too. I guess, only the actual variables (such as u, p, omega, k) will use the underrelaxation, but maybe someone can correct me. (Perhaps just run a case with and without it and see, if the residuals change). |
Thanks for your answer :-)
Firstly I´m going to try without nuSg. Would you say, that the relaxation-factors from above are good, or would you use the factors used for simple: p=0.3 (U|k|omega)=0.7 Best regards Tobi |
I normally run pimple and simple with the latter, but as I said: Just try. If it is not stable, you need to reduce. If it is stable it should usually converge faster the higher the underrelaxation factors are.
|
Quote:
Therefore you dont need it when using RANSE k-omega-modell. |
Yes Hans, thats true, but what about setting a "nut" underrelaxation parameter? Does that have any influence ? And if so, why would you do that, if already "k" and "epsilon" or "k" and "omega" are underrelaxated.
|
In k-omega model, nut is calculated from k and omega in an algebraic equation. So I think defining a relaxation factor for nut makes no sense.
Just my thoughts. |
All times are GMT -4. The time now is 19:11. |