CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   failed validating Cd for 2D airfoil (https://www.cfd-online.com/Forums/openfoam-solving/214240-failed-validating-cd-2d-airfoil.html)

Sereff January 24, 2019 09:52

failed validating Cd for 2D airfoil
 
Dear Foamers:

The case I was simulating is air flow around a 2D airfoil with Re = 1e6, with following:
  • Solver: pimpleFoam
  • Model: k-Omega SST
  • Meshing: snappHexMesh

I've been getting pretty covincing results for Cl, but not for Cd. I've read from some source that Cd is far more sensitive to mesh quality than Cl does. I assume that 'high mesh quality' means more refined mesh as in more cells and finer boundary layers?

And talking about boundary layers, I took a closer look at the yPlus value for my simulations on different mesh, please correct me if I was using wrong BCs,
case 1,30<y+<70, with boundary conditions on airfoil:
  • U: type noSlip;
  • turbulence properties: type xxWallFunction; value $internalField
case 2,0.09<y+<1.2, with boundary conditions on airfoil:
  • U: type noSlip;
  • turbulence properties: type fixedValue; value uniform 1e-15
But all my simulations are having trouble with Cd. Any idea?

Kind regards,

Sereff January 24, 2019 10:50

Attached are my Cd and Cl results for different aoa.

yambanshee January 25, 2019 01:41

Drag is a difficult property to get right in airfoil simulations. It is very dependent on the wall boundary layer, and so using a wall function will likely not work.


What is meant by a high mesh quality is details such as density, minimal grading near the walls, low non-orthogonality and skewness, and cells aligned to the direction of flow as much as possible.


Also, as the drag is so dependent on the boundary layer, accurate turbulence modelling at the boundary is also very important. Unfortunately, all RANS formulations are simply an approximation of turbulence effects, so this introduces enhanced difficulty in getting a correct boundary layer

Sereff January 25, 2019 02:44

Hi, thank you for your anwer.

So I assume with k-Omega-SST, it is expected to have a rather bad Cd prediction. But do you know what boundary condition should I use for velocity on the surface of airfoil? I've been told if wall function were to apply then I should use 'zeroGradient' for airfoil surface instead of 'noSlip'.

Kind regards,

yambanshee January 25, 2019 04:10

AFAIK if you are interested in the effect of the wall, a no slip condition is always needed


All times are GMT -4. The time now is 07:06.