Some effects of domain size and boundary conditions on the accuracy of airfoil simulations

This paper investigates a specific case of one of the most popular fluid dynamic simulations, the incompressible flow around an airfoil (NACA 0012 here) at a high Reynolds number ( $$6 \times 10^6$$ ). OpenFOAM software was used to study the effect of domain size and four common choices of boundary conditions on airfoil lift, drag, surface friction, and pressure. We also examine the relation between boundary conditions and the velocity, pressure, and vorticity distributions throughout the domain. In addition to the common boundary conditions, we implement the “point vortex” boundary condition that was introduced many years ago but is now rarely used. We also applied the point vortex condition for the outlet pressure instead of using the traditional Neumann condition. With the airfoil generating significant lift at incidence angles of $$5^\circ , 10^\circ$$ , and $$14^\circ$$ , we confirm a previous finding that the boundary conditions combine with domain size to produce an induced (pressure) drag. The change in the pressure drag with domain size is significant for the commonly-used boundary conditions but is much smaller for the point vortex alternative. The point vortex boundary condition increases the execution time, but this is more than offset by the reduction in domain size needed to achieve a specified accuracy in the lift and drag. This study also estimates the error in total drag and lift due to domain size and shows it can be almost eliminated using the point vortex boundary condition. We also used the impulse form of the momentum equations to study the relation between drag and lift and spurious vorticity, which is generated as a result of using non-exact boundary conditions. These equations reveal that the spurious vorticity throughout the domain is associated with cancelling circulation around the domain boundaries.