|
[Sponsors] |
October 10, 2012, 15:00 |
Unstable Spallart-Allmaras
|
#1 |
New Member
Peter Doings
Join Date: Sep 2012
Location: Portugal
Posts: 4
Rep Power: 13 |
Hello,
I'm trying to implement the Spalart-Allmaras model in my code, but I ran into some difficulties. It seems the diffusion coefficient, which depends on the eddy viscosity, makes the solution process very unstable. Atm I'm using a structured grid with Finite Volume discretization and linear interpolation at the boundaries to calculate the fluxes, with GMRES to solve the system of eqs. Will using a higher order scheme help? Is the use of flux limiters recommended? If so, which one is the best? I plan on applying the method to unstructured grids later on, by the way. Thank you for your comments. |
|
October 10, 2012, 15:11 |
|
#2 |
Senior Member
Chris DeGroot
Join Date: Nov 2011
Location: Canada
Posts: 414
Rep Power: 17 |
Higher order schemes are almost always more difficult to converge. Don't move on to higher order schemes until you can get a first order upwind scheme to converge. That being said, a flux limiter will be helpful when you are working with the higher order schemes.
Also it sounds like your problem is with diffusion so changing the advection discretization probably will not help. Maybe you need to apply under-relaxation to the eddy viscosity to smooth things out and keep things stable. That is to say, you should blend in some of the previous estimate to prevent the eddy viscosity from changing too abruptly between iterations. |
|
October 10, 2012, 20:06 |
|
#3 |
New Member
Peter Doings
Join Date: Sep 2012
Location: Portugal
Posts: 4
Rep Power: 13 |
Well, I'm already using quite low values for the under-relaxation factor.
Another question: limiting the viscosity to positive values seems to help. Is this frowned upon or too unrealistic to be used in practice? |
|
October 11, 2012, 03:12 |
|
#4 |
Member
Serge A. Suchkov
Join Date: Oct 2011
Location: Moscow, Russia
Posts: 74
Blog Entries: 5
Rep Power: 14 |
Take a look http://turbmodels.larc.nasa.gov/spalart.html especially notes.
For test your implementation take a look http://turbmodels.larc.nasa.gov/flatplate_sa.html http://turbmodels.larc.nasa.gov/shear_sa.html http://turbmodels.larc.nasa.gov/bump_sa.html http://turbmodels.larc.nasa.gov/bump3d_sa.html
__________________
OpenHyperFLOW2D Project |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Spalart and Allmaras Paper | Mohamed | Main CFD Forum | 7 | July 10, 2014 16:39 |
CFX Spalart Allmaras turbulence models | aweizazuji | CFX | 9 | September 24, 2013 10:53 |
Convergence on unstable steady states | ekor | FLUENT | 0 | August 24, 2010 05:50 |
spalart allmaras paper | hafiz | Main CFD Forum | 1 | April 10, 2007 03:47 |
extracting K and Epsilon from Spalart Allmaras | Sundar | Main CFD Forum | 4 | July 31, 2006 02:02 |