|
[Sponsors] |
February 6, 2012, 23:01 |
Numerical Schemes for Vortex Capturing
|
#1 |
Member
|
Hi all,
I am trying to model the wing tip vortices from a round tip wing ( NACA 0012) using OpenFOAM. I tried both the Spalart Allmaras and KW-SST models and realized that the results are no where near the experimental results. When I tweaked the fvSchemes file, I got results which are at least comparable to the experiment but they are still quite a bit away. I understand that the diffusion of the vortex depends majorly on the diverging schemes. I am using Gauss QUICK and when I use MUSCL, the solution blows up. Can anyone suggest what schemes I can employ to get lower diffusion? Kalyan |
|
February 7, 2012, 01:48 |
|
#2 |
Senior Member
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 283
Rep Power: 18 |
You might also want to look at your mesh in this region since it will play an important role.
|
|
February 7, 2012, 09:38 |
|
#4 |
Senior Member
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 283
Rep Power: 18 |
Out of curiosity, do you have some pictures of the mesh in this region?
|
|
February 7, 2012, 10:45 |
|
#5 |
Member
|
Vinz,
Here are some pics of the mesh. The first mesh is around the wing with a cylinder over the round tip. The second mesh is in the wake behind the wing, with the fine mesh in the region of the vortex. Please excuse the weirdness of the cells, since paraview doesn't render well. The mesh was made in snappyHex wing.png wake.png |
|
February 7, 2012, 11:01 |
|
#6 |
Senior Member
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 283
Rep Power: 18 |
Ok, thanks for images. A small tip for paraview, you can use "extract cell by region" instead of "Slice". It will keep the cells instead of cutting them.
When you say "the results are no where near the experimental results." What are the parameters that your are looking at? forces? coefficients? flow topology? vortex shedding? |
|
February 7, 2012, 11:08 |
|
#7 |
Member
|
Vinz,
Thanks for the tip. The main parameters which I was comparing were the Normalized axial velocity through the vortex core and the cross flow velocity through the vortex core. Immediately behind the wing, the norm. axial velocity is around 1.77*U, I get around 1.74*U which is very reasonable but at around x/c=0.6 downstream of the wing, the experimental results have 1.69*U and I have around 0.2-0.3 *U which is a clear indicator of over diffusion either due to the turbulence modeling or the numerical method. I used both the SA and the KW-SST model with Gauss linear and Gauss QUICK conditions. Gauss QUICK gives slightly better results compared to the Gauss linear but not drastically different. One thing to note, the residuals have flat lined in just around 600 iterations are about 10e-2 which is not really a good convergence criterion considering the basic geometry involved,and they dont even change. I am running a steady state simulation by the way. Do you think I need to run a transient ? |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Numerical schemes for free surface flows (VOF) | botp | OpenFOAM | 2 | March 11, 2011 15:27 |
numerical schemes | Valerio | FLUENT | 2 | July 28, 2008 04:02 |
Kinetic schemes and numerical dissipation | Praveen | Main CFD Forum | 0 | September 6, 2002 07:09 |
Behaviour of Numerical Schemes | Ravi. B. R | Main CFD Forum | 3 | September 22, 2001 10:18 |
Standard for checking and testing numerical schemes? | X. Ye | Main CFD Forum | 7 | August 31, 1999 17:05 |