Schoder Stefan, Weitz Michael, Maurerlehner Paul, Hauser Alexander, Falk Sebastian, Kniesburges Stefan, Döllinger Michael, Kaltenbacher Manfred
Institute of Mechanics and Mechatronics, TU Wien, Austria.
Division of Phoniatrics and Pediatric Audiology at the Department of Otorhinolaryngology Head & Neck Surgery, University Hospital Erlangen, Friedrich-Alexander-University Erlangen-Nürnberg, Germany.
J Acoust Soc Am. 2020 Feb;147(2):1179. doi: 10.1121/10.0000785.
A hybrid aeroacoustic approach was developed for the efficient numerical computation of human phonation. In the first step, an incompressible flow simulation on a three-dimensional (3 D) computational grid, which is capable of resolving all relevant turbulent scales, is performed using STARCCM+ and finite volume method. In the second step, the acoustic source terms on the flow grid are computed and a conservative interpolation to the acoustic grid is performed. Finally, the perturbed convective wave equation is solved to obtain the acoustic field in 3 D with the finite element solver CFS++. Thereby, the conservative transformation of the acoustic sources from the flow grid to the acoustic grid is a key step to allow coarse acoustic grids without reducing accuracy. For this transformation, two different interpolation strategies are compared and grid convergence is assessed. Overall, 16 simulation setups are compared. The initial (267 000 degrees of freedom) and the optimized (21 265 degrees of freedom) simulation setup were validated by measurements of a synthetic larynx model. To conclude, the total computational time of the acoustic simulation is reduced by 95% compared to the initial simulation setup without a significant reduction of accuracy, being 7%, in the frequency range of interest.
为了高效地对人类发声进行数值计算,开发了一种混合气动声学方法。第一步,使用STARCCM+和有限体积法在能够解析所有相关湍流尺度的三维(3D)计算网格上进行不可压缩流模拟。第二步,计算流场上的声源项,并对声学网格进行守恒插值。最后,使用有限元求解器CFS++求解扰动对流波动方程,以获得三维声场。因此,将声源从流场网格保守转换到声学网格是允许使用粗声学网格而不降低精度的关键步骤。对于这种转换,比较了两种不同的插值策略并评估了网格收敛性。总体而言,比较了16种模拟设置。通过对合成喉模型的测量验证了初始(267000自由度)和优化(21265自由度)模拟设置。总之,与初始模拟设置相比,声学模拟的总计算时间减少了95%,在感兴趣的频率范围内精度仅降低7%,精度没有显著降低。