CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CFX

Multigrid cycles and numerical scheme

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   January 9, 2014, 10:01
Default Multigrid cycles and numerical scheme
  #1
Nkl
Member
 
Nikola
Join Date: Sep 2013
Location: Madrid, Spain
Posts: 60
Rep Power: 4
Nkl is on a distinguished road
Hi everyone!

My supervisor is very cryptic...a solution of one of my models is not reaching convergence after more than 10k iterations. After telling me to keep on with patience for many days, today he told me to "play with computational parameters" such as type and number of numerical scheme and multigrid cycles. Where should I look? How should I change these parameters?
Thanks in advance
Nkl is offline   Reply With Quote

Old   January 9, 2014, 10:20
Default
  #2
Senior Member
 
Lance
Join Date: Mar 2009
Posts: 589
Rep Power: 12
Lance is on a distinguished road
The default multigrid settings should be fine, but the choice of numerical scheme can definitely affect convergence. You'll find them under solver control...
Besides numerical schemes there are many factors affecting convergence, including mesh quality, boundary conditions, time step size, etc...
Lance is offline   Reply With Quote

Old   January 9, 2014, 10:29
Default
  #3
Nkl
Member
 
Nikola
Join Date: Sep 2013
Location: Madrid, Spain
Posts: 60
Rep Power: 4
Nkl is on a distinguished road
I have created two models and the only difference is the inlet geometry. In the second one, the one where convergence is slow, the inlet is divided in 16 circular entrances (so mass flow rate has been divided by 16). For the rest the models are identical. Where lies the problem then?
I will try to change numerical scheme, thank you!
Nkl is offline   Reply With Quote

Old   January 12, 2014, 06:56
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,238
Rep Power: 96
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
You are probably getting a complex flow feature occurring in the difficult to converge case. Maybe shock waves getting somewhere tricky, maybe a difficult separation.

What to do is discussed in the FAQ: http://www.cfd-online.com/Wiki/Ansys...gence_criteria

And Lance is correct in that you should not fiddle with the multigrid parameters. The problem will be more fundamental - eg mesh quality, time step size, physical setup, boundary conditions.
ghorrocks is offline   Reply With Quote

Old   January 12, 2014, 08:01
Default
  #5
Nkl
Member
 
Nikola
Join Date: Sep 2013
Location: Madrid, Spain
Posts: 60
Rep Power: 4
Nkl is on a distinguished road
Is it maybe a mistake dividing the mass flow rate by 16 because I have 16 inlets instead of only 1? I was thinking that maybe the mass flow rate is referred to the whole sum of inlets. Thanks for the advices!
Nkl is offline   Reply With Quote

Old   January 12, 2014, 17:43
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,238
Rep Power: 96
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
The mass flow is evenly distributed over the inlet, whatever its size, and no matter how many surfaces it is spread over.
ghorrocks is offline   Reply With Quote

Old   January 13, 2014, 08:39
Default
  #7
Nkl
Member
 
Nikola
Join Date: Sep 2013
Location: Madrid, Spain
Posts: 60
Rep Power: 4
Nkl is on a distinguished road
Thank you very much, you are far more helpful than my supervisor!
Now I have run the simulation with a physical timescale of 5s (it was 0.25s) and the correct mass flow rate. I have reached convergence after only 160 iterations and the plots seem to be ok. The only problem now lays in the yplus values, which are in the order of magnitude of e+5. This sounds strange because the mesh settings are the same I used for the model with one inlet, for which I have good yplus values (less than 100). Thank you in advance again!
Nkl is offline   Reply With Quote

Old   January 13, 2014, 17:50
Default
  #8
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,238
Rep Power: 96
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
As a general note: The coupled solver in CFX generally converges very quickly. For most single phase flows, if you are not converged (or at least well on the way there) by 100 iterations then something is wrong. For multiphase simulations and other simulations with tricky coupling this could be ore like 1000 iterations.

But forget about running it for 10000 iterations like you used to have to do with SIMPLE based solvers. With SIMPLE solvers if you had convergence difficulties you fiddled with the solver parameters (under relaxation, linear solver etc) to get convergence.

With CFX you almost never adjust the under relaxation or solver parameters. If you have convergence difficulties you look at your mesh quality, time step size, initial condition and double precision numerics.
ghorrocks is offline   Reply With Quote

Old   January 14, 2014, 17:48
Default
  #9
Nkl
Member
 
Nikola
Join Date: Sep 2013
Location: Madrid, Spain
Posts: 60
Rep Power: 4
Nkl is on a distinguished road
Now that convergence is good I have one last question. My two models are identical from all points of view (fluid domain, boundary conditions, solver controls), they only differ for the geometry of inlet. When I plot the results, they are very similar. How is it possible that for my second model (with 16inlets) the yplus values are around 3000, while in the other model they are below 100??
Thanks for all the help!
Nkl is offline   Reply With Quote

Old   January 14, 2014, 17:59
Default
  #10
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,238
Rep Power: 96
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
Where is the high value of y+? Post an image of what you are seeing.
ghorrocks is offline   Reply With Quote

Old   January 14, 2014, 18:07
Default
  #11
Nkl
Member
 
Nikola
Join Date: Sep 2013
Location: Madrid, Spain
Posts: 60
Rep Power: 4
Nkl is on a distinguished road
They are even around 200 000. Here is the image!
Attached Images
File Type: jpg Immagine.jpg (43.5 KB, 8 views)
Nkl is offline   Reply With Quote

Old   January 14, 2014, 18:16
Default
  #12
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,238
Rep Power: 96
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
y+ is a function of mesh size and flow velocity. So either the mesh size or the flow velocity has changed. You need to use the post processor to find out which.

Also your y+ plot looks very blotchy. This is a sure sign of inadequate mesh resolution. You have not done a mesh sensitivity study, have you? If that is the case your results are little better than random numbers at the moment.
ghorrocks is offline   Reply With Quote

Old   January 14, 2014, 18:24
Default
  #13
Nkl
Member
 
Nikola
Join Date: Sep 2013
Location: Madrid, Spain
Posts: 60
Rep Power: 4
Nkl is on a distinguished road
I have assigned face sizings, each with element size of 0.01 m. The inflation option is total thickness with 15 layers, 1.7 growth rate and a maximum thickness of 0.004 m. With orthogonal quality I have an average of 0.85. I used these settings because with the other model I achieved values of yplus lower than 100.
I have tried to increase both element size and growth rate, but yplus did not change. Should I try mesh adaptation?
Nkl is offline   Reply With Quote

Old   January 14, 2014, 18:27
Default
  #14
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,238
Rep Power: 96
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
But did you look at the mesh you generated?

My comment on mesh sensitivity applies to the whole mesh, not just the inflation layers. For more details see the FAQ http://www.cfd-online.com/Wiki/Ansys..._inaccurate.3F and this is an excellent description of what to do http://journaltool.asme.org/Template...umAccuracy.pdf
ghorrocks is offline   Reply With Quote

Old   January 14, 2014, 18:46
Default
  #15
Nkl
Member
 
Nikola
Join Date: Sep 2013
Location: Madrid, Spain
Posts: 60
Rep Power: 4
Nkl is on a distinguished road
What does it mean exactly "have take a look at your mesh?". From the beginning of my simulations I have changed advection scheme from high resolution to upwind, turbulence model from k-epsilon to SST, boundary condition type from outlet to opening in order to achieve convergence. I don't understand why with the other model I have good results with the same parameters. I will try to change other parameters. Thanks for all the help!
Nkl is offline   Reply With Quote

Old   January 14, 2014, 19:05
Default
  #16
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,238
Rep Power: 96
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
I mean use the post processor to look at the mesh and see its quality and size. Especially compare the two simulations which you report are giving different results despite similar setups.

Also be careful of using upwind differencing. You are going to struggle to get accurate answers with that.
ghorrocks is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 17:21.