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.

The presented filtering technique is proposed to detect errors and correct outliers inside the acoustic sources, respectively, the first time derivative of the incompressible pressure obtained from large eddy simulations with prescribed vocal fold motion using overlay mesh methods. Regarding the perturbed convective wave equation, the time derivative of the incompressible pressure is the primary sound source in the human phonation process. However, the incompressible pressure can be erroneous and have outliers when fulfilling the divergence-free constraint of the velocity field. This error is primarily occurring for non-conserving prescribed vocal fold motions. Therefore, the method based on a continuous stationary random process was designed to detect rare events in the time derivative of the pressure. The detected events are then localized and treated by a defined window function to increase their probability. As a consequence, the data quality of the non-linearly filtered data is enhanced significantly. Furthermore, the proposed method can also be used to assess convergence of the aeroacoustic source terms, and detect regions and time intervals, which show a non-converging behavior by an impulse-like structure.

Porous acoustic absorbers have excellent properties in the low-frequency range when positioned in room edges, therefore they are a common method for reducing low-frequency reverberation. However, standard room acoustic simulation methods such as ray tracing and mirror sources are invalid for low frequencies in general which is a consequence of using geometrical methods, yielding a lack of simulation tools for these so-called edge absorbers. In this article, a validated finite element simulation model is presented, which is able to predict the effect of an edge absorber on the acoustic field. With this model, the interaction mechanisms between room and absorber can be studied by high-resolved acoustic field visualizations in both room and absorber. The finite element model is validated against transfer function data computed from impulse response measurements in a reverberation chamber in style of ISO 354. The absorber made of Basotect^® is modeled using the Johnson-Champoux-Allard-Lafarge model, which is fitted to impedance tube measurements using the four-microphone transfer matrix method. It is shown that the finite element simulation model is able to predict the influence of different edge absorber configurations on the measured transfer functions to a high degree of accuracy. The evaluated third-octave band error exhibits deviations of 3.3–4.1 dB computed from third-octave band averaged spectra.

Sound generation in human phonation and the underlying fluid–structure–acoustic interaction that describes the sound production mechanism are not fully understood. A previous experimental study, with a silicone made vocal fold model connected to a straight vocal tract pipe of fixed length, showed that vibroacoustic coupling can cause a deviation in the vocal fold vibration frequency. This occurred when the fundamental frequency of the vocal fold motion was close to the lowest acoustic resonance frequency of the pipe. What is not fully understood is how the vibroacoustic coupling is influenced by a varying vocal tract length. Presuming that this effect is a pure coupling of the acoustical effects, a numerical simulation model is established based on the computation of the mechanical-acoustic eigenvalue. With varying pipe lengths, the lowest acoustic resonance frequency was adjusted in the experiments and so in the simulation setup. In doing so, the evolution of the vocal folds’ coupled eigenvalues and eigenmodes is investigated, which confirms the experimental findings. Finally, it was shown that for normal phonation conditions, the mechanical mode is the most efficient vibration pattern whenever the acoustic resonance of the pipe (lowest formant) is far away from the vocal folds’ vibration frequency. Whenever the lowest formant is slightly lower than the mechanical vocal fold eigenfrequency, the coupled vocal fold motion pattern at the formant frequency dominates.

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.

Impedance tube measurements are a widespread method for determining the sound absorption coefficient (SAC) of porous materials for normal sound incidence. The measurement method is standardised in ISO 10534-2. However, the standards offer limited guidance on sample preparation and mounting. Many impedance tubes have a circular cross-section and can be mounted at various angles. Therefore, it is crucial to investigate the influence of the mounting angle on the SAC and its interaction with diameter imperfections and different thicknesses using different materials. Following the dataset FOAM 01 this work documents the creation of a new dataset (called FOAM 02) containing SAC measurements from the ISO 10534-2 two-microphone method. The following factors are considered: two different materials, four thicknesses, three diameters, and four rotation angles. For each combination of material, thickness, and diameter, three specimens are produced, resulting in 72 specimens. Each specimen is measured three times at each rotation angle, yielding 864 SAC measurements. The dataset contains the SAC measurements and one-hot encoded label vectors, and is publicly available. A cutting device is proposed to saw cylindrical specimens accurately on a band saw with variable thickness and diameter. The workshop drawings of the cutting device are available as supplementary material.