CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   FLUENT (http://www.cfd-online.com/Forums/fluent/)
-   -   3D Separation Model using K-Omega SST Divergence (http://www.cfd-online.com/Forums/fluent/113080-3d-separation-model-using-k-omega-sst-divergence.html)

TempestCFD February 12, 2013 01:27

3D Separation Model using K-Omega SST Divergence
 
1 Attachment(s)
Hey guys,

I'm trying to model a 3-D wing in Fluent to study the effects of boundary restabilisation via Vortex Generators.

I started by using the k-e model, before realising that this was inappropriate for the high angles of attack I would be studying. I have now tried to switch to k-Omega SST to better capture boundary layer separation, and refined my mesh (50 layers of inflation, first off wall thickness of 1e-06 due to a fairly high Reynolds number of ~ 10^5). My domain is about 20 chord lengths radius (I have a cylinder) of my wing, and I'm using a pressure-far-field condition.

My solution doesnt last too long before it starts to diverge, and I have no clue why. I have been searching all over the forums for a good two days trying to solve this issue, but without any luck.
Any ideas?



-------------------------------------------------------------------------
Also, here are some of the parameters that I have set (just to reference):

1) Pressure Based, Steady
2) Energy - On, SST k-omega
3) Air - ideal gas (for pressure-far-field condition)
4) Velocity - 13.44 m/s
5) Courant # of 5


Attached is a picture of part of my mesh.


Thank you for your help

cfd seeker February 12, 2013 08:47

A common problem, not a big issue. Jut reduce your courant number to 1 at the start of solution and slowly increase it. By doing so, I hope your problem will be solved.

TempestCFD February 12, 2013 13:22

Progress but a new fault
 
1 Attachment(s)
cfd seeker, thank you for your reply.


The initialisation now works, however the solution begins to diverge soon after, giving new errors (see screenshot attached).

I'm curious about the behavior of the continuity residuals at the beginning as well (staying at exactly 1 initially)?


I will try reducing my courant number further, but am open to other ideas as well.


Thanks in advance.


Frustrated but still determined,


TempestCFD

cfd seeker February 13, 2013 12:57

what is the mach no? in which sotware you made your mesh? what is the minimum qulity of your mesh? also trying with 1st order upwind for all equztions at the start of the solution.

TempestCFD February 13, 2013 13:24

The mach number is low (M = 0.0386).

I meshed the setup in Ansys Mesher.

My elements are fairly good quality and I have a fine mesh of about 6.5 million elements. Anything specific you want to know? Skewness ect?

I re-ran the simulation with a courant number of 0.1 and relaxed the energy from 1 -> 0.9 and the solution runs okay. The continuity residuals seem to flat line at 10e-02 though which is not satisfactory in my eyes. I'm going to look at improving the elements near the trailing edge of the airfoil as its sharp and not a great transition to the tet elements behind that section.

Any additional concepts or ideas would be greatly appreciated.


Thank you

cfd seeker February 14, 2013 08:35

Continuity convergence is very slow and secondly courant no and convergence are directly propotional, so let the solution run as it is and check the convergence of global parameters like CL,CD,CM etc.

TempestCFD February 14, 2013 08:39

I tried running the solution with a velocity inlet and pressure outlet instead, and the continuity went down to 10e-04. Much better.

My Cl and Cd seem exceptionally high however. I can't seem to find any issues with my reference values that would cause this.


Thank you for your guidence by the way, it is greatly appreciated.

cfd seeker February 15, 2013 14:21

Quote:

My Cl and Cd seem exceptionally high however
Have you resolved the boundary layer properly? what's the range of wall y+ ?
Have you scaled the mesh after importing it in Fluent?

TempestCFD February 15, 2013 14:29

My boundary layer (in my opinion) has been resolved well, even too well.

My y+ is ~ 0.15 - 0.89


What do you mean by "scaling the mesh after importing it into fluent" ?

cfd seeker February 15, 2013 15:10

Wall y+ 0.15 - 0.89 Wow it's great...what is the total mesh size? How much RAM does you machine have....
What is the units in which you created geometry and then meshed it? If it's other than "meters" then you have to scale it in Fluent, follow Mesh>Scale

TempestCFD February 15, 2013 15:25

Total mesh size is huge (about 5.5 million elements), 12 GB of ram, lol. It lags a bit :D

Everything was done in metres

Far February 16, 2013 08:44

Why energy is on? Are you using pressure based -coupled solver?

TempestCFD February 16, 2013 18:33

Sorry, I forgot to mention, I switched energy off, switching from pressure-far-field condition to a velocity-inlet/pressure-outlet setup.

Yes, pressure based, coupled.


All times are GMT -4. The time now is 15:33.