This paper presents the first numerical investigation via direct numerical simulation of the tone noise phenomenon occurring in the flow past laminar airfoils. This phenomenon corresponds to the radiation of discrete acoustic tones in some specific flow conditions, and has received much attention since the 1970s, and several experimental studies have been carried out to identify and understand the underlying physical mechanisms. However, several points remain to be clarified in order to provide a complete explanation of its origin. The flow around a two-dimensional NACA0012 airfoil is considered in order to have a deeper understanding of the tone noise phenomenon. Consistently with previous experimental studies, it is shown that depending on the Reynolds number and angle of attack, two different types of acoustic spectrum are observed: one which exhibits a broadband contribution with a dominant frequency together with a sequence of regularly spaced discrete frequencies, while the other one is only characterized by a simple broadband contribution. The first configuration is typical of the tone noise phenomenon. The present work shows that in this case, the mean flow on the pressure side of the airfoil exhibits a separation bubble near the trailing edge and the main tone frequency is close to the most amplified frequency of the boundary layer. The mechanism proposed in previous works for the main tone generation – which implies the existence of a separation bubble at the pressure side – is therefore validated by numerical simulation. On the other hand, the analysis of the suction side boundary layer reveals that there is no separation and that the most amplified frequency is different from the main tonal one. However, the suction side boundary layer is highly receptive to the tone frequency. Finally, an original explanation for the existence of the secondary discrete frequencies observed in the radiated pressure spectrum is given. They are associated to a bifurcation of the airfoil wake from a symmetric to a non-symmetric vortex pattern. A possible explanation for the existence of this bifurcation is the interaction between the disturbances which are the most amplified by the suction side boundary layer and those originating in the forcing of the suction side flow by the main tone noise mechanism.

A lattice Boltzmann method (LBM) with enhanced stability and accuracy is presented for various Hermite tensor-based lattice structures. The collision operator relies on a regularization step, which is here improved through a recursive computation of nonequilibrium Hermite polynomial coefficients. In addition to the reduced computational cost of this procedure with respect to the standard one, the recursive step allows to considerably enhance the stability and accuracy of the numerical scheme by properly filtering out second- (and higher-) order nonhydrodynamic contributions in under-resolved conditions. This is first shown in the isothermal case where the simulation of the doubly periodic shear layer is performed with a Reynolds number ranging from ${10}^4$ to ${10}^6$, and where a thorough analysis of the case at $\text{Re}=3\times {10}^4$ is conducted. In the latter, results obtained using both regularization steps are compared against the Bhatnagar-Gross-Krook LBM for standard (D2Q9) and high-order (D2V17 and D2V37) lattice structures, confirming the tremendous increase of stability range of the proposed approach. Further comparisons on thermal and fully compressible flows, using the general extension of this procedure, are then conducted through the numerical simulation of Sod shock tubes with the D2V37 lattice. They confirm the stability increase induced by the recursive approach as compared with the standard one.