For the clinical analysis of underlying mechanisms of voice disorders, we developed a numerical aeroacoustic larynx model, calledsimVoice, that mimics commonly observed functional laryngeal disorders as glottal insufficiency and vibrational left-right asymmetries. The model is a combination of the Finite Volume (FV) CFD solver Star-CCM+ and the Finite Element (FE) aeroacoustic solver CFS++.simVoicemodels turbulence using Large Eddy Simulations (LES) and the acoustic wave propagation with the perturbed convective wave equation (PCWE). Its geometry corresponds to a simplified larynx and a vocal tract model representing the vowel /a/. The oscillations of the vocal folds are externally driven. In total, 10 configurations with different degrees of functional-based disorders were simulated and analyzed. The energy transfer between the glottal airflow and the vocal folds decreases with an increasing glottal insufficiency and potentially reflects the higher effort during speech for patients being concerned. This loss of energy transfer may also have an essential influence on the quality of the sound signal as expressed by decreasing sound pressure level (SPL), Cepstral Peak Prominence (CPP), and Vocal Efficiency (VE). Asymmetry in the vocal fold oscillations also reduces the quality of the sound signal. However,simVoiceconfirmed previous clinical and experimental observations that a high level of glottal insufficiency worsens the acoustic signal quality more than oscillatory left-right asymmetry. Both symptoms in combination will further reduce the quality of the sound signal. In summary,simVoiceallows for detailed analysis of the origins of disordered voice production and hence fosters the further understanding of laryngeal physiology, including occurring dependencies. A current walltime of 10 h/cycle is, with a prospective increase in computing power, auspicious for a future clinical use ofsimVoice.

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 enormous computational power and time required for simulating the complex phonation process preclude the effective clinical use of computational larynx models. The aim of this study was to evaluate the potential of a numerical larynx model, considering the computational time and resources required. Using Large Eddy Simulations (LES) in a 3D numerical larynx model with prescribed motion of vocal folds, the complicated fluid-structure interaction problem in phonation was reduced to a pure flow simulation with moving boundaries. The simulated laryngeal flow field is in good agreement with the experimental results obtained from authors’ synthetic larynx model. By systematically decreasing the spatial and temporal resolutions of the numerical model and optimizing the computational resources of the simulations, the elapsed simulation time was reduced by 90% to less than 70 h for 10 oscillation cycles of the vocal folds. The proposed computational larynx model with reduced mesh resolution is still able to capture the essential laryngeal flow characteristics and produce results with sufficiently good accuracy in a significant shorter time-to-solution. The reduction in computational time achieved is a promising step towards the clinical application of these computational larynx models in the near future.

Modeling fluid-structure-acoustic interaction is a highly challenging task and currently there is a lack of available benchmark datasets. For instance, to understand the processes of human phonation, it is essential to fully understand the fluid-structure-acoustic interaction process. In this article, a synthetic human phonation model is presented for benchmarking numerical methods in the field of aerodynamics, aeroacoustics and the highly-complex fluid-structure-acoustic interaction process of the voice production. The objective is to present the wide range of experimental data available, including the model geometry, the material characteristics, quantitiesreferred to the vocal fold dynamics, fluid mechanical quantities and the acoustic field. Firstly, the experimental setup is specified. Secondly, the experimental data is described in detail, accompanied by illustrations. Having the dataset in hand, finally, it is shown how to use this data to validate a computational phonation model called simVoice successfully. In conclusion, this benchmarking dataset offers the opportunity to validate structural dynamics, aerodynamics and aeroacoustics of a highly-complex fluid-structure-acoustic interaction simulation, obtained from different mathematical formulations and numerical procedures.

Most sports like surfing are highly developed. It is necessary to tease the last percentages out of the competitors and equipment—in the case of surfing the surfboard-fin-system—to win competitions or championships. In this computational investigation, a parameter study of the positioning of the two rear fins within a 4-fin configuration with fixed front fins on a surfboard is executed to find appropriate fin positions for specific surf situations. Four different inflow velocities are investigated. The RANS and URANS models combined with the SST $k-\omega$ turbulence model, which is available within the computational fluid dynamics (CFD) package STAR-CCM+, are used to simulate the flow field around the fins for angles of attack (AoA) between 0° and 45°. The simulation results show that shifting the rear fins toward the longitudinal axis of the surfboard lowers the maximum lift. Surfboards with 4-fin configurations are slower in nearly the whole range of AoA due to a higher drag force but produce a higher lift force compared to the 3-fin configuration. The lift and drag forces increase significantly with increasing inflow velocity. This study shows a significant influence of the rear fin positioning and the inflow velocity on lift and drag performance characteristics.

Understanding the risk of infection by routine medical examination is important for the protection of the medical personnel. In this study we investigated respiratory particles emitted by patients during routine otolaryngologic procedures and assessed the risks for the performing physician. We developed two experimental setups to measure aerosol and droplet emission during rigid/flexible laryngoscopy, rhinoscopy, pharyngoscopy, otoscopy, sonography and patient interview for subjects with and without masks. A high-speed-camera setup was used to detect ballistic droplets (approx. > 100 µm) and an aerosol-particle-sizer was used to detect aerosol particles in the range of 0.3 µm to 10 µm. Aerosol particle counts were highly increased for coughing and slightly increased for heavy breathing in subjects without masks. The highest aerosol particle counts occurred during rigid laryngoscopy. During laryngoscopy and rhinoscopy, the examiner was exposed to increased particle emission due to close proximity to the patient’s face and provoked events such as coughing. However, even during sonography or otoscopy without a mask, aerosol particles were expelled close to the examiner. The physician’s exposure to respiratory particles can be reduced by deliberate choice of examination technique depending on medical indication and the use of appropriate equipment for the examiners and the patients (e.g., FFP2 masks for both).