The flow-induced sound sources of human voice production are investigated based on a validated voice model. This analysis is performed using a hybrid aeroacoustic workflow based on the perturbed convective wave equation. In the first step, the validated 3D incompressible turbulent flow simulation is computed by the finite volume method using STARCCM+. In a second step, the aeroacoustic sources are evaluated and studied in detail. The formulation of the sound sources is compared to the simplification (neglecting the convective sources) systematically using time-domain and Fourier-space analysis. Additionally, the wave equation is solved with the finite element solver openCFS to obtain the 3D sound field in the acoustic far-field. During the detailed effect analysis, the far-field sound spectra are compared quantitatively, and the flow-induced sound sources are visualized within the larynx. In this contribution, it is shown that the convective part of the sources dominates locally near the vocal folds (VFs) while the local time derivative of the incompressible pressure is distributed in the whole supra-glottal area. Although the maximum amplitude of the time derivative is lower, the integral contribution dominates the sound spectrum. As a by-product of the detailed perturbed convective wave equation source study, we show that the convective source term can be neglected since it only reduces the validation error by 0.6%. Neglecting the convective part reduces the algorithmic complexity of the aeroacoustic source computation of the perturbed convective wave equation and the stored flow data. From the source visualization, we learned how the VF motion transforms into specific characteristics of the aeroacoustic sources. We found that if the VFs are fully closing, the aeroacoustic source terms yield the highest dynamical range. If the VFs are not fully closing, VFs motion does not provide as much source energy to the flow-induced sound sources as in the case of a healthy voice. As a consequence of not fully closing VFs, the cyclic pulsating velocity jet is not cut off entirely and therefore turbulent structures are permanently present inside the supraglottal region. These turbulent structures increase the broadband component of the voice signal, which supports research results of previous studies regarding glottis closure and insufficient voice production.

The hybrid aeroacoustic approach is an efficient way to address the issue of the disparity of scales in Computational AeroAcoustics (CAA) at low Mach numbers. In the present paper, three wave equations governing propagation of flow-induced sound of low Mach number flows, namely the Perturbed Convective Wave Equation (PCWE), Ribner’s Dilatation (RIB) equation, and Lighthill’s wave equation, are applied using the Finite Element Method (FEM). An airflow through a circular pipe with a half-moon-shaped orifice at three operating flow speeds is considered, where validation data from measurements on a dedicated test rig is available. An extensive analysis of the flow field is provided based on the results of the incompressible flow simulation. The resulting acoustic source terms are investigated, and the relevant source term contributions are determined. The results of the acoustic propagation simulations revealed that the PCWE and RIB are best suited for the present task. The overall deviation of the predicted pressure spectra from the measured mean values amounted to 2.26 and 2.13 times the standard deviation of the measurement compared to 3.55 for Lighthill’s wave equation. Besides reliably predicting the flow-induced sound, the numerical procedure of source term computation is straightforward for PCWE and RIB, where the source term contributions, shown to be relevant, solely consist of time derivatives of the incompressible pressure. In contrast, the Lighthill source term involves spatial derivatives and, thus, is strongly dependent on the spatial resolution and the numerical method actually used for approximating these terms.

The process of voice production is a complex process and depends on the correct interaction of the vocal folds and the glottal airstream inducing the primary voice source, which is subsequently modulated by the vocal tract. Due to the restricted access to the glottis, not all aspects of the three-dimensional process can be captured by measurements without influencing the measurement object. Hence, the application of a numerical tool capturing the physical process of phonation can provide an extended database for voice treatment and, therefore, can contribute to an increased effectiveness of voice treatment. However, such numerical models involve complex and demanding procedures to model the material behavior and the mechanical contact of the vocal folds and to realize moving boundaries of the involved physical domains. The present paper proposes a numerical model calledsimVoice, which circumvents these computational expenses by prescribing the experimentally obtained vocal fold motion within the simulation. Additionally, a hybrid approach for sound computation further enhances the computational efficiency and yields good agreement with acoustic measurements. An analysis of the computational workloads suggests that the key factor for a further increase in efficiency is an optimized flow simulation and source term computation.

This article deals with large-eddy simulations of three-dimensional incompressible laryngeal flow followed by acoustic simulations of human phonation of five cardinal English vowels, /ɑ, æ, i, o, u/. The flow and aeroacoustic simulations were performed in OpenFOAM and in-house code openCFS, respectively. Given the large variety of scales in the flow and acoustics, the simulation is separated into two steps: (1) computing the flow in the larynx using the finite volume method on a fine moving grid with 2.2 million elements, followed by (2) computing the sound sources separately and wave propagation to the radiation zone around the mouth using the finite element method on a coarse static grid with 33 000 elements. The numerical results showed that the anisotropic minimum dissipation model, which is not well known since it is not available in common CFD software, predicted stronger sound pressure levels at higher harmonics, and especially at first two formants, than the wall-adapting local eddy-viscosity model. The model on turbulent flow in the larynx was employed and a positive impact on the quality of simulated vowels was found.

Acoustic emissions play a major role in the usability of many product categories. Therefore, mitigating the emitted sound directly at the source is paramount to improve usability and customer satisfaction. To reliably predict acoustic emissions, numerical methods such as the boundary element method (BEM) are employed, which allow for predicting, e.g., the acoustic emission into the free field. BEM algorithms need appropriate boundary conditions to couple the sound field with the structural motion of the vibrating body. In this contribution, firstly, an interpolation scheme is presented, which allows for appropriate interpolation of arbitrary velocity data to the computational grid of the BEM. Secondly, the free-field Helmholtz problem is solved with the open-source BEM software framework NiHu. The forward coupling between the device of interest and BEM is based on the surface normal velocities (i.e., a Neumann boundary condition). The BEM simulation results are validated using a previously established aeroacoustic benchmark problem. Furthermore, an application to a medical device (knee prosthesis frame) is presented. Furthermore, the radiated sound power is evaluated and contextualized with other low-cost approximations. Regarding the validation example, very good agreements are achieved between the measurements and BEM results, with a mean effective pressure level error of $0.63$ dB averaged across three microphone positions. Applying the workflow to a knee prosthesis frame, the simulation is capable of predicting the acoustic radiation to four microphone positions with a mean effective pressure level error of $1.52$ dB.