|
[Sponsors] |
February 20, 2013, 07:18 |
Divergence problem
|
#1 |
Member
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 14 |
"Already posted this yesterday but no response that's why"
Hello, I am getting this problem with the 3D VOF (air jet impinging over water). The simulation is good till 1.4e-2s but after that it gives the following error. The time step size i am using is 1e-6 with 40 time steps. The setting are as follows: 1. Double with SST 2. Surface tension 0.073 n/m 3. inlet velocity 96.672 m/s 4. PISO Scheme with VF(Mod HRIC) Momentum (2nd ord) 5. URF Pressure 0.3 Density 0.5 Body Forces 1 VF 0.5 In 2D for the same thing, i had mass loss problem but i was able to find the solution with the help of few people here. Would be thankful for your kind reply. Regards, Smaras ===>>>>turbulent viscosity limited to viscosity ratio of 1.000000e+05 in 240257 cells # Divergence detected in AMG solver: x-momentum -> Increasing relaxation sweeps! Divergence detected in AMG solver: x-momentum Divergence detected in AMG solver: y-momentum Divergence detected in AMG solver: z-momentum Divergence detected in AMG solver: vof-1 Divergence detected in AMG solver: x-momentum Divergence detected in AMG solver: y-momentum Divergence detected in AMG solver: z-momentum Divergence detected in AMG solver: vof-1 Divergence detected in AMG solver: x-momentum Divergence detected in AMG solver: y-momentum Divergence detected in AMG solver: z-momentum Divergence detected in AMG solver: vof-1 Divergence detected in AMG solver: x-momentum Divergence detected in AMG solver: y-momentum Divergence detected in AMG solver: z-momentum Divergence detected in AMG solver: vof-1 Divergence detected in AMG solver: x-momentum Divergence detected in AMG solver: y-momentum Divergence detected in AMG solver: z-momentum Divergence detected in AMG solver: vof-1 Divergence detected in AMG solver: x-momentum Divergence detected in AMG solver: y-momentum Divergence detected in AMG solver: z-momentum Divergence detected in AMG solver: vof-1 turbulent viscosity limited to viscosity ratio of 1.000000e+05 in 251038 cells Primitive Error at Node 0: floating point exception Primitive Error at Node 1: floating point exception Primitive Error at Node 2: floating point exception Primitive Error at Node 3: floating point exception Primitive Error at Node 4: floating point exception Primitive Error at Node 5: floating point exception<<<<==== |
|
February 20, 2013, 07:33 |
|
#2 |
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
1) Does each time step converge until these errors appear? Can you show a residual plot?
2) Did you make a solution animation? Can you see something curious in the results until the errors?
__________________
The skeleton ran out of shampoo in the shower. |
|
February 20, 2013, 08:43 |
|
#3 | |
Member
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 14 |
Quote:
1. Yes, from very start it was converging, but just before the diverging occurs three or four time steps has no convergence rather turbulent viscosity. 2. Yes, i did only did for the penetration depth. Please see the attachments. I thinks it's more clear in that. Thanks once again and looking forward for your reply. Regards, Azam |
||
February 20, 2013, 08:47 |
hi
|
#4 |
Senior Member
shoeb khan
Join Date: Nov 2011
Posts: 179
Rep Power: 14 |
Sometimes divergence problem is due to bad meshing that you have used .
Try to refine your mesh and run the simulation and let me know. try to decrease the under relaxation factors also tell what is the courunt number in your simulation. Regards shk |
|
February 20, 2013, 09:19 |
|
#5 |
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
I think shk is right here. It looks like the fluid firstly moves through a domain of "good mesh" where the solver is able to converge sufficiently, but then enters some bad cells, where the mess starts.
You should a) post a picture of your mesh b) show us the courant number (i think Fluent can show cell courant number) c) in the meantime let the simulation run with a smaller time step. Edit: before you do c) please try, if you can reduce the residuals right from the start each time step better than now. It looks like you get "continuity" to about 1.0e-3. Can you get better convergence? How? By increasing the number of iterations per time step or only by reducing the time step?
__________________
The skeleton ran out of shampoo in the shower. |
|
February 20, 2013, 09:58 |
|
#6 | |
Member
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 14 |
Quote:
This is my mesh a) Now just as shk as said i rechecked in fluent i am getting mesh Rodriguez i have used smaller steps but this is stall point where it crashes. No matter the step size is small or large. I even tried number of iteration but no use. Now what should i do???? regards, Smaras |
||
February 20, 2013, 10:10 |
|
#7 |
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
1) Post a picture of the Courant number shortly before it crashes.
2) Why do you use a prism layer at the wall? I don't think you expect any flow along that wall, right? As I understand it, a prism layer isn't helpful in such a case. 3) Can you use SIMPLE algorithm for your model? In my experience this is the most robust one. 4) Do you have ICEM for meshing?
__________________
The skeleton ran out of shampoo in the shower. |
|
February 20, 2013, 10:32 |
|
#8 | |
Member
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 14 |
Quote:
2,3) Ok i will try that without prism layers. (after your reply) 4) Yes i am using ICEM for meshing. |
||
February 21, 2013, 02:57 |
|
#9 |
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
If you have ICEM, you could make an excellent hexa-mesh. Your geometry is really simple and blocking would be straight forward. Don't use tet meshs unless you are not able to do the blocking. They are numerically low-grade.
__________________
The skeleton ran out of shampoo in the shower. |
|
February 21, 2013, 04:33 |
|
#10 | |
Member
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 14 |
Quote:
I know Rodriguez but this isn't the real geometry i will be working on more complex geometry, and the requirement for my thesis is tetra meshing. This is just a mock-up or example to make practice which are based on previous research papers. I am today re-meshing it without the prism layer and plus more refined mesh. Further i have 2 questions: 1. I need to know one more thing i am using densities to get precise result for the jet flow and its impingement. as shown previously. Is it creating the problem???? 2. Further after opening opening the geometry in fluent there is seperation line within the air region. and when i saw the velocity profile @ 1.4e-2sec the velocity is reaching the line. Is this meshing problem?????? Regards, Smaras |
||
February 21, 2013, 04:40 |
|
#11 |
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
1) I don't understand the question. Do you mean, if a variable density can cause problems? Then, the answer is yes.
2) What do you mean by "separation line"?
__________________
The skeleton ran out of shampoo in the shower. |
|
February 21, 2013, 04:47 |
|
#12 | |
Member
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 14 |
Quote:
2.) The line that starts from top till the water bed. |
||
February 21, 2013, 04:54 |
|
#13 |
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27 |
But what kind of line is this? Did you paint it?
Also I am curious why the courant number at the top, next to the needle is so high. I guess the whole top area is a pressure outlet, right? Why would you have such a high courant number there?
__________________
The skeleton ran out of shampoo in the shower. |
|
February 21, 2013, 05:03 |
|
#14 | |
Member
Smaras
Join Date: Sep 2012
Location: Germany
Posts: 78
Rep Power: 14 |
Quote:
2. The top area is outlet. And that i don't know. It's appeared just after divergence this is the image before divergence |
||
Tags |
vof model |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Quarter Burner mesh with periosic condition | SamCanuck | FLUENT | 2 | August 31, 2011 11:34 |
Meshing and divergence problem | ghost | FLUENT | 9 | February 5, 2010 12:24 |
divergence problem | vincent | FLUENT | 3 | August 3, 2006 15:44 |
strange divergence when solving multiphase problem | tanghao | FLUENT | 2 | July 27, 2006 19:47 |
divergence problem | Ayyappan.T | FLUENT | 2 | May 16, 2005 12:10 |