This contribution benchmarks the aeroacoustic workflow of the perturbed convective wave equation and the Ffowcs Williams and Hawkings analogy in Farassat’s 1A version for a low-pressure axial fan. ...Thereby, we focus on the turbulence modeling of the flow simulation and mesh convergence concerning the complete aeroacoustic workflow. During the validation, good agreement has been found with the efficiency, the wall pressure sensor signals, and the mean velocity profiles in the duct. The analysis of the source term structures shows a strong correlation to the sound pressure spectrum. Finally, both acoustic sound propagation models are compared to the measured sound field data.
For the clinical analysis of underlying mechanisms of voice disorders, we developed a numerical aeroacoustic larynx model, called
, 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++.
models 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,
confirmed 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,
allows 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 of
.
In the human phonation process, acoustic standing waves in the vocal tract can influence the fluid flow through the glottis as well as vocal fold oscillation. To investigate the amount of acoustic ...back-coupling, the supraglottal flow field has been recorded via high-speed particle image velocimetry (PIV) in a synthetic larynx model for several configurations with different vocal tract lengths. Based on the obtained velocity fields, acoustic source terms were computed. Additionally, the sound radiation into the far field was recorded via microphone measurements and the vocal fold oscillation via high-speed camera recordings. The PIV measurements revealed that near a vocal tract resonance frequency
, the vocal fold oscillation frequency
(and therefore also the flow field's fundamental frequency) jumps onto
. This is accompanied by a substantial relative increase in aeroacoustic sound generation efficiency. Furthermore, the measurements show that
-
-coupling increases vocal efficiency, signal-to-noise ratio, harmonics-to-noise ratio and cepstral peak prominence. At the same time, the glottal volume flow needed for stable vocal fold oscillation decreases strongly. All of this results in an improved voice quality and phonation efficiency so that a person phonating with
-
-coupling can phonate longer and with better voice quality.
Numerical modeling of the human phonatory process has become more and more in focus during the last two decades. The increase in computational power and the use of high-performance computation (HPC) ...yielded more complex models being closer to the actual fluid-structure-acoustic interaction (FSAI) within the human phonatory process. However, several different simulation approaches with varying mathematical complexity and focus on certain parts of the phonatory process exist. Currently, models are suggested based on ordinary differential equations (reduced order models) but also on partial differential equations based on continuum mechanics as e.g. the Navier–Stokes equations for the flow discretized by Finite-Volume or Finite-Element-Methods. This review will illuminate current trends and recent progress within the area. In summary, the ultimate simulation model satisfying all physiological needs and scientific opinions still has to be developed.
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.
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.
In an aeroacoustic simulation of human voice production, the effect of the sub-grid scale (SGS) model on the acoustic spectrum was investigated. In the first step, incompressible airflow in a 3D ...model of larynx with vocal folds undergoing prescribed two-degree-of-freedom oscillation was simulated by laminar and Large-Eddy Simulations (LES), using the One-Equation and Wall-Adaptive Local-Eddy (WALE) SGS models. Second, the aeroacoustic sources and the sound propagation in a domain composed of the larynx and vocal tract were computed by the Perturbed Convective Wave Equation (PCWE) for vowels u: and i:. The results show that the SGS model has a significant impact not only on the flow field, but also on the spectrum of the sound sampled 1 cm downstream of the lips. With the WALE model, which is known to handle the near-wall and high-shear regions more precisely, the simulations predict significantly higher peak volumetric flow rates of air than those of the One-Equation model, only slightly lower than the laminar simulation. The usage of the WALE SGS model also results in higher sound pressure levels of the higher harmonic frequencies.
In this paper, a transient and compressible solver based on the Lattice-Boltzmann Method/Very Large Eddy Simulation (LBM/VLES) approach is employed to predict unsteady flow physics and flow-induced ...noise generation of a low-pressure axial fan. Aerodynamic and aeroacoustic measurements provided by the European Acoustics Association (EAA) benchmark platform are used for validation purposes. Boundary and design operating conditions are applied to the numerical model to replicate the experimental setup. Simulation and experimental data are compared, showing an excellent agreement in terms of fan efficiency with less than 1% deviation, as well as broadband and tonal noise within 0.7 dBA in terms of overall sound pressure level. An advanced post-processing analysis is performed to shed light on the noise generation mechanism in tip clearance. It is observed that both fine random turbulent structures and large coherent vortices are generated in the tip gap. The continuous impingement of the fine turbulence with the following blades and the blade itself is responsible for the radiation of broadband noise, while the interaction between the large coherent tip vortices, spinning at a lower angular velocity with respect to the fan shaft, and the following blades leads to the generation of narraowband peaks at sub-harmonics of the blade-passing frequencies. Finally, a beamforming analysis further confirms that the main noise sources are located in the blade tip clearance and tip regions.
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.
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.