[pvtpilot]

I am currently working on an independent study in CFD based around using OpenFOAM and my task right now is to model flow around a SD7003 airfoil at Re=60000. So far, I have been successful with my laminar models leading up to this using icoFoam but this is my first experience with the turbulent solvers. I have created the airfoil and free stream and meshed it using Salome imported it. For BC I have inlet, outlet, symmetry on the free stream upper, lower, inner and outer boundaries and wall functions (since wall gives the error of you should use wall functions). Some of the problems/questions that have come up are:

I would like to solve this as a 2D case but I always run into a problem when trying to use empty for the faces bounding the z axis. I get an error saying essentially due to the faces empty is the wrong condition so I set it to symmetry (which is probably wrong).

An overall question would be: what is the best solver to use (k-epsilon or LES)? Right now I'm just trying to get sensible results in simpleFoam but I end up getting results that just diverge once I turn turbulence on.

For an airfoil what is the best Div scheme to use?

I'm sure there are some issues I'm missing right now and as you can tell these are some newbie problems I'm having but I appreciate any help you can provide.