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

First order simulation better than second order

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 1, 2012, 13:42
Default First order simulation better than second order
  #1
Member
 
Join Date: Jan 2012
Location: Edmonton, CA
Posts: 87
Rep Power: 14
Torque_Converter is on a distinguished road
Send a message via AIM to Torque_Converter
I am finding that the first order solution is much closer to test than second order. I am using the same convergence criteria (max res 5e-4), same boundary conditions, same fluid, etc. Just changing the convergence order. I start with first order/upwind, then from those results go high resolution on both turbulence and pressure, momentum. First order is within 8% of test torque, and second order is 115% over pedicted torque. Mesh changes definitely affect the first order result, but not second order.

This is all a little concerning, as I don't want to just use first order with a mesh that happens to match test as this serves little purpose in development when without test data.
Torque_Converter is offline   Reply With Quote

Old   May 1, 2012, 20:11
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,871
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
The apparent accuracy of the first order solution is probably just luck. You need to refine the mesh such that you get a mesh independant solution - you will probably find this occur for the second order system with a coarser mesh than the first order one.

Your comment that the first order system is more sensitive to mesh changes supports this.
ghorrocks is offline   Reply With Quote

Old   May 1, 2012, 23:16
Default
  #3
Member
 
Join Date: Jan 2012
Location: Edmonton, CA
Posts: 87
Rep Power: 14
Torque_Converter is on a distinguished road
Send a message via AIM to Torque_Converter
Thanks, getting lucky seems right since it obviously not consistent.
I have a question thou about the coarseness of the mesh thou. The mesh is pretty fine right now, with inflation layer first layer of .003 mm, 15 layers, and body mesh of .4 mm with the domain being approximately 80x80 mm at its largest dimensions. This produces a mesh of about 2.5 million. This probably overdoes it. However, how can I know when the solution is mesh independent. Let's say 500,000 cells gives a second order solution that is correct. But then when i jump to 1.5 million cells it goes back to being over 100% overpredicted. That seems mesh dependent.
Torque_Converter is offline   Reply With Quote

Old   May 2, 2012, 02:05
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,871
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
When a mesh becomes too fine the following problems occur:
1) numerical round off means spurious errors start creeping in. This can be helped by using double precision numeric.
2) Sometimes the reduced dissipation of a finer mesh is not fully offset by the turbulence model. This means that psuedo-turbulent fluctuations start occuring and affects the result. You can normally see this in the the results file as the residuals do not achieve tolerance (for steady state runs) or you can see transient structures in the flow which (for transient runs).
3) For transient simulations a finer mesh will probably need a smaller timestep to keep the Courant number about the same. You might need to do a new time step sensitivity check for the finer mesh.
4) Likewise, the CFX convergence tolerance is pretty good at being independant of mesh size, but sometimes it is not perfect. You might need to do a new convergence tolerance sensitivity check for the finer mesh.

This is why sensitivity checks are slow. Once you have checked convergence tolerance, time step size and mesh you have probably found a new time step size, tolerance and mesh is required and you have to do the exercise again. It is an iterative process as you converge on the best mesh, convergence tolerance and time step size.
ghorrocks is offline   Reply With Quote

Old   May 2, 2012, 13:15
Default
  #5
Member
 
Join Date: Jan 2012
Location: Edmonton, CA
Posts: 87
Rep Power: 14
Torque_Converter is on a distinguished road
Send a message via AIM to Torque_Converter
Thank you for all your help. Starting with a coarser mesh and looking at first and second order I see that second order is consistently about 100% high, while first order is all over the place. This shows me that as you said first order is just luck if you happen to stop at some mesh count, while second order is mesh independent past a certain point.

I think perhaps the 100% overprediction is not a facet of first being somehow better than second order (which I did not feel good about before!) nor is it a facet of mesh, but somehow the physics isn't right. Perhaps there is cavitation present, or some other factor that reduces the torque production.
Torque_Converter is offline   Reply With Quote

Old   May 2, 2012, 19:12
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,871
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Quote:
somehow the physics isn't right
I agree. Something important has been missed.

Good work for actually looking into this and figuring this out - many people don't bother and they fool themselves into thinking their model is accurate.
ghorrocks is offline   Reply With Quote

Old   May 2, 2012, 22:00
Default
  #7
Member
 
Join Date: Jan 2012
Location: Edmonton, CA
Posts: 87
Rep Power: 14
Torque_Converter is on a distinguished road
Send a message via AIM to Torque_Converter
Since you have been so helpful, for myself and on the site in general I wanted to let you know any interesting find on this in case you were curious.

The parasolid I have has sharp perpendicular wall intersections and sharp blade geometry. I then got a look at the actual device and noticed everything was rounded off to 3mm fillets and has drafts. I asked around and found this was necessary since it is sand cast and needs these features to be manufacterable. Tomorrow my new direction will be to study the effects of these relatively large variations from "ideal" sharp geometry. I imagine it will be significant considering it is in high pressure and key momentum change regions.

Thanks again
Torque_Converter is offline   Reply With Quote

Old   January 8, 2013, 05:07
Default
  #8
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27
immortality is on a distinguished road
hi glenn
I have a 3mm duct with 6cm in length.i have set courant number to .1 for convergency but results are not as expected.how should i do?my grids are 20*500.should i make them finer?how much?and CO lower?
Thanks.
immortality is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Why is my fluent simulation giving Cp values above 1??? wils Main CFD Forum 3 March 26, 2010 10:12
Replacing mesh while running a simulation akultane CFX 1 November 15, 2009 14:46
What's the best order to run this simulation in? siw CFX 1 November 4, 2009 20:42
FSI TWO-WAY SIMULATION Smagmon CFX 1 March 6, 2009 14:24
2nd order conservative schemes taw Main CFD Forum 1 September 16, 2008 08:05


All times are GMT -4. The time now is 03:46.