Too much iterations for k, epsilon with Pointwise mesh
2 Attachment(s)
Hello, Foamers
I'm trying to solve flow past a floating square box. The solver is interFoam(Multiphase water, air). The domain, suface mesh, volume mesh are generated by Pointwise and I set boundary conditions as below : Inlet, oulet, atmosphere : patch bottom, midplane, side : symmetry box : wall Volume mesh is composed of tetrahedron mesh. Mesh quality checked by checkMesh utility was fine with some skew cells. The problem is that Iterations for alpha.water and k, epsilon is 1000 and they blow up later... I already faced this situation when the mesh was structured mesh.. I thought unstructured mesh will be a solution for this situation but... nothing changed Here is terminal output Code:
Time = 0.000119976 Have anyone experienced this too much iteration problem? If so, how can I solve this problem? :confused: Best regards, tigger |
This is quite unusual... the linear solver doesn't converge at all.
Did you try to run a different solver for k and epsilon? Code:
"(U|k|epsilon|omega|UFinal|kFinal|epsilonFinal|omegaFinal)" |
Just additionally, this is probably not the reason for your problems but you should always start a case with the uncorrected laplacian scheme, which is numerically more stable.
Code:
laplacianSchemes |
There is no reason why an unstructured tet mesh should work better than a structured mesh. I would never prever a tet mesh for free surface flows. Can you plot some screenshots of your mesh and from your initialized solution (especially alpha field)?
|
Quote:
|
Quote:
|
Thanks! Philipp
1 Attachment(s)
Quote:
I tried to use PBiCG solver you suggested Now the K, epsilon iterations become stable Code:
DILUPBiCG: Solving for alpha.water, Initial residual = 9.41257114753943e-07, Final residual = 6.69250718355886e-11, No Iterations 1 You can see the residuals in attached picture. I have to find solutions for this situation.. Anyway, Thanks again! |
Here is the picture
Quote:
I attached alpha field screenshot. I'm not sure this picture is that you wanted to see... Well, by the way, I had same trouble(too many iterations for k, epsilon..) when I used Hex mesh(Structured grid)... Thanks! |
You need to post the log output of the blowup (and some iterations before)...
|
Jason, I can't see the screenshots. Can you upload again?
|
1 Attachment(s)
Quote:
Actually, I stopped the calculation because the residual is getting weird.. Thanks, tigger |
1 Attachment(s)
Quote:
Here it is Thanks, tigger |
ok, can you do one more with mesh grid on?
|
Can you post the fvSchemes and fvSolution that you used for to get that log output?
|
2 Attachment(s)
Quote:
|
2 Attachment(s)
Quote:
|
Several things:
1) You did not change the laplacian scheme to uncorrected. 2) You should use numerically safe settings for the other schemes as well, such as "Gauss upwind" instead of linearUpwind and vanLeer. If they work, you can start to introduce better schemes, one by one. 3) Using p_rghFinal with relTol 0 is a waste of time, I guess. I would commend out that "p_rghFinal" block. 4) This is a PIMPLE based solver, right? You use "nOuterCorrectors 1" which basically means, that you don't run PIMPLE but PISO. But PISO needs a Courant number of less than 1, which is already violated in the 3rd or 4th time step. Thus, you need to reduce the time step or increase the numer of outer (PIMPLE-)iterations per time step. If you do the second, you should set turbOnFinalIterOnly to "no" and also use some safer relaxation factors for the pressure (such as p 0.3) for the beginning. Try to set nOuterCorrectors 15 or so and see if that runs stable. |
Quote:
3) relTol 0 on p_rghFinal will be removed! 4) I'm using 'interFoam'solver and this solver may be based on PIMPLE algorithm.. I missed that the 'nOuterCorrector 1' means calculation is running PISO algorithm.. Thanks for your kind advice :D Those 4 advice will be adopted to my case and I'll update results! |
3 Attachment(s)
Well,:rolleyes:
every advices applied to new case! 1) uncorrected laplacian scheme, 2) Gauss upwind scheme for all fvSchemes 3) relTol 0 was removed 4) nOuterCorrectors 15, turbOnFinalIterOnly no Now, I can see that bounding k is getting bigger and bigger Finally, the calculation blow up and openfoam do Iteration 1000 suddenly in final pimple iteration (Iteration 15)... Here is my log file and fvSolution & fvSchemes... * When I run simulation with turbOnFinalIterOnly yes the volume fraction went above 1.. Can someone explain how turbOnFinalIterOnly affect the pimple iteration or recommend some paper to study? |
Jason, I don't know the solver you use and also don't know what it solves for ;)
But: Is there any differential equation for "alpha" that is solved? Like a transport equation? You need to set an under-relaxation factor for PIMPLE for all the values (such as k, epsilon, p, alpha), otherwise this will probably be very unstable. Please try that. For debugging: You can have a look at the residuals "before" and "after" setting the new under-relaxation factor in the log file. If they change, you set the right value. |
All times are GMT -4. The time now is 12:45. |