|
[Sponsors] |
Divergence when implementing Spalart-Allmaras turbulence model |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
December 4, 2023, 06:04 |
Divergence when implementing Spalart-Allmaras turbulence model
|
#1 |
New Member
Join Date: Nov 2021
Posts: 8
Rep Power: 4 |
Hi,
I was trying to implement the Spalart-Allmaras turbulence model into my FVM rotor solver, which can work well with the Baldwin-Lomax model. The version of the S-A model I used is the "Standard" model on https://turbmodels.larc.nasa.gov/spalart.html#sa. I used the LUSGS scheme for time marching, and convective and viscous fluxes in the S-A model are 1st upwind and 2nd central scheme, respectively. The freestream nut is set to 3*nu as recommended. The case I am running is a hovering state of a 2-bladed rotor. Now the computation suddenly diverges after 5 revolutions, when the force coefficients seem to have achieved a periodic solution. It is weird to me. Lowering the CFL number helps to delay the divergence but it still crashes finally. Has anyone encountered similar problems? What are the possible causes of such a sudden divergence? Any ideas are welcomed! Thanks in advance!!! |
|
December 4, 2023, 07:27 |
|
#2 | |
Senior Member
|
Quote:
|
||
December 4, 2023, 07:55 |
|
#3 |
New Member
Join Date: Nov 2021
Posts: 8
Rep Power: 4 |
Thanks, prior to the rotor case I have already tested my implementation with some 2D airfoils and 3D wing test cases provided in the NASA validation archive. In these cases my computation results show good agreement with experiments and converge well. When it is applied to this rotor case, the sudden divergence occurs
|
|
December 4, 2023, 08:07 |
|
#4 | |
Senior Member
|
Quote:
Do you have any other equation based turbulence model? Is the model implementation based, say, on previously used components or is it just a whole new implementation? If you were developing this for a commercial product you should have tens of cases in between so, no, two or three cases is not typically enough to be confident on a model implementation. One thing that usually requires attention in the SA implementation is the non linear diffusive term, which requires rearrangement to be better implemented in a conservative FV code, but I have doubts it actually is your problem. My suggestion is to carefully assess what are the elements of your failing case and test the code with intermediate mixes of ingredients. In the end, SA is among the most stable models around and should not fail if other options don't as well. |
||
December 4, 2023, 08:29 |
|
#5 | |
New Member
Join Date: Nov 2021
Posts: 8
Rep Power: 4 |
Quote:
Actually I am trying a whole new implementation without any other equation-based model, for academic use. And I did rearrange the diffusive term as advised. I truly need to test the intermediate elements of my code first, before entering the mixed case. |
||
December 4, 2023, 09:51 |
|
#7 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,782
Rep Power: 71 |
Just to test your moving mesh, does it works without any turbulence model?
|
|
December 4, 2023, 20:58 |
|
#8 | |
New Member
Join Date: Nov 2021
Posts: 8
Rep Power: 4 |
Quote:
I output the rotor flow fields of several time steps before it diverges, as attached. The density and pressure contour exhibit some singularities in the wake region(cut boundary), which is the region of large turbulent viscosities. And this is the region where divergence occurs first. |
||
December 5, 2023, 08:54 |
|
#9 | |
Senior Member
|
Quote:
Also, does the same code and case work when mesh rotation is set to 0? |
||
December 5, 2023, 20:39 |
|
#10 | |
New Member
Join Date: Nov 2021
Posts: 8
Rep Power: 4 |
Quote:
And thanks for that hint, I will examine my code setting the rotation to 0, to check if the moving mesh part works well. |
||
Tags |
spalart-allmaras model, sudden diverge |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Error: WorkBench Error: Could not handle event: SolutionStatusUpdate | Kieyo | Fluent Multiphase | 0 | November 9, 2022 23:58 |
Implementing a new turbulence model | Sakun | OpenFOAM | 15 | November 11, 2021 13:56 |
Fail to converge when solving with a fabricated solution | zizhou | FLUENT | 0 | March 22, 2021 06:33 |
Spalart Allmaras turbulence model | Venkatesh Devaraj | Main CFD Forum | 0 | January 20, 2015 23:22 |
Spalart Allmaras Turbulence model - Need article! | CFD Student | Main CFD Forum | 0 | December 27, 2004 04:28 |