Modeling of the flow over aerofoil profiles at low Reynolds numbers is difficult due to the complex physics associated with the laminar flow separation mechanism. Two major problems arise in the estimation of profile drag: (1) the drag force at low Reynolds numbers is extremely small to be measured in a wind tunnel by force balance techniques, (2) the profile drag is usually calculated by pressure integration, hence the skin friction component of drag is excluded. In the present work, three different 4-digit NACA aerofoils are investigated. Measurements are conducted in an open-ended subsonic wind tunnel, while numerical work is performed by time Reynolds-averaged Navier Stokes (RANS) coupled with the laminar-kinetic-energy ( K-kl-w) turbulence model. The influence of the flow separation bubbles and transition locations on the profile drag is discussed and addressed. This paper gives important insights into importance of measurements at low Reynolds numbers for better aerodynamic loads predictions.