Restraining vocal fold vertical motion reduces source-filter interaction in a two-mass model.

Previous experimental studies suggested that restraining the vocal fold vertical motion may reduce the coupling strength between the voice source and vocal tract. In this study, the effects of vocal fold vertical motion on source-filter interaction were systematically examined in a two-dimensional two-mass model coupled to a compressible flow simulation. The results showed that when allowed to move vertically, the vocal folds exhibited subharmonic vibration due to entrainment to the first vocal tract acoustic resonance. Restraining the vertical motion suppressed this entrainment. This indicates that the vertical mobility of the vocal folds may play a role in regulating source-filter interaction.

Characterizing infection risk in a restaurant environment due to airborne diseases using discrete droplet dispersion simulations.

The use of masks as a measure to control the spread of respiratory viruses has been widely acknowledged. However, there are instances where wearing a mask is not possible, making these environments potential vectors for virus transmission. Such environments can contain multiple sources of infection and are challenging to characterize in terms of infection risk. To address this issue, we have developed a methodology to investigate the role of ventilation in reducing the infection risk in such environments. We use a restaurant setting as a representative scenario to demonstrate the methodology. Using implicit large eddy simulations along with discrete droplet dispersion modeling we investigate the impact of ventilation and physical distance on the spread of respiratory viruses and the risk of infection. Our findings show that operating ventilation systems, such as mechanical mixing and increasing physical distance between subjects, can significantly reduce the average room infection risk and number of newly infected subjects. However, this observation is subject to the transmissibility of the airborne viruses. In the case of a highly transmissible virus, the use of mechanical mixing may be inconsequential when compared to only fresh air ventilation. These findings provide valuable insights into the mitigation of infection risk in situations where the use of masks is not possible.

Comparison of one-dimensional and three-dimensional glottal flow models in left-right asymmetric vocal fold conditions.

While the glottal flow is often simplified as one-dimensional (1D) in computational models of phonation to reduce computational costs, the 1D flow model has not been validated in left-right asymmetric vocal fold conditions, as often occur in both normal and pathological voice production. In this study, we performed three-dimensional (3D) and 1D flow simulations coupled to a two-mass model of adult male vocal folds and compared voice production at different degrees of left-right stiffness asymmetry. The flow and acoustic fields in 3D were obtained by solving the compressible Navier-Stokes equations using the volume penalization method with the moving vocal fold wall as an immersed boundary. Despite differences in the predicted flow pressure on vocal fold surface between the 1D and 3D flow models, the results showed reasonable agreement in vocal fold vibration patterns and selected voice outcome measures between the 1D and 3D models for the range of left-right asymmetric conditions investigated. This indicates that vocal fold properties play a larger role than the glottal flow in determining the overall pattern of vocal fold vibration and the produced voice, and the 1D flow simplification is sufficient in modeling phonation, at least for the simplified glottal geometry of this study.

Numerical investigation of effects of tongue articulation and velopharyngeal closure on the production of sibilant [s].

A numerical simulation of sibilant /s/ production with the realistically moving vocal tract was conducted to investigate the flow and acoustic characteristics during the articulation process of velopharyngeal closure and tongue movement. The articulation process was simulated from the end of /u/ to the middle of /s/ in the Japanese word /usui/, including the tongue elevation and the velopharyngeal valve closure. The time-dependent vocal tract geometry was reconstructed from the computed tomography scan. The moving immersed boundary method with the hierarchical structure grid was adopted to approach the complex geometry of the human speech organs. The acoustic characteristics during the co-articulation process were observed and consistent with the acoustic measurement for the subject of the scan. The further simulations with the different closing speeds of the velopharyngeal closure showed that the far-field sound during the co-articulation process was amplified with the slower closing case, and the velum closure speed was inverse proportional to the sound amplitude with the slope value of - 35.3 dB s/m. This indicates possible phonation of indistinguishable aeroacoustics sound between /u/ and /s/ with slower velopharyngeal closure.

Estimation of sibilant groove formation and sound generation from early hominin jawbones.

The speech production capability of sibilant fricatives of early hominin was assessed by interpolating the modern human vocal tract to an Australopithecine specimen based on the jawbone landmarks, and then simulating the airflow and sound generation. The landmark interpolation demonstrates the possibility to form the sibilant groove in the anterior part of the oral tract, and results of the aeroacoustic simulation indicate that the early hominins had the potential to produce the fricative broadband noise with a constant supply of airflow to the oral cavity, although the ancestor's tongue deformation ability is still uncertain, and the results are highly speculative.

Quantifying the COVID19 infection risk due to droplet/aerosol inhalation.

The dose-response model has been widely used for quantifying the risk of infection of airborne diseases like COVID-19. The model has been used in the room-average analysis of infection risk and analysis using passive scalars as a proxy for aerosol transport. However, it has not been employed for risk estimation in numerical simulations of droplet dispersion. In this work, we develop a framework for the evaluation of the probability of infection in droplet dispersion simulations using the dose-response model. We introduce a version of the model that can incorporate the higher transmissibility of variant strains of SARS-CoV2 and the effect of vaccination in evaluating the probability of infection. Numerical simulations of droplet dispersion during speech are carried out to investigate the infection risk over space and time using the model. The advantage of droplet dispersion simulations for risk evaluation is demonstrated through the analysis of the effect of ambient wind, humidity on infection risk, and through a comparison with risk evaluation based on passive scalars as a proxy for aerosol transport.

Numerical investigation of effects of incisor angle on production of sibilant /s/.

The effects of the inclination angle of the incisor on the speech production of the fricative consonant /s/ was investigated using an implicit compressible flow solver. The hierarchical structure grid was applied to reduce the grid generation time for the vocal tract geometry. The airflow and sound during the pronunciation of /s/ were simulated using the adaptively switched time stepping scheme, and the angle of the incisor in the vocal tract was changed from normal position up to 30°. The results showed that increasing the incisor angle affected the flow configuration and moved the location of the high turbulence intensity region thereby decreased the amplitudes of the sound in the frequency range from 8 to 12 kHz. Performing the Fourier transform on the velocity fluctuation, we found that the position of large magnitudes of the velocity at 10 kHz shifted toward the lip outlet when the incisor angle was increased. In addition, separate acoustic simulations showed that the shift in the potential sound source position decreased the far-field sound amplitudes above 8 kHz. These results provide the underlying insights necessary to design dental prostheses for the production of sibilant fricatives.

Aeroacoustic differences between the Japanese fricatives [ɕ] and [ç].

To elucidate the linguistic similarity between the alveolo-palatal sibilant [ɕ] and palatal non-sibilant [ç] in Japanese, the aeroacoustic differences between the two consonants were explored via experimentation with participants and analysis using simplified vocal tract models. The real-time magnetic resonance imaging (rtMRI) observations of articulatory movements demonstrated that some speakers use a nearly identical place of articulation for /si/ [ɕi] and /hi/ [çi]. Simplified vocal tract models were then constructed based on the data captured by static MRI, and the model-generated synthetic sounds were compared with speaker data producing [ɕ] and [ç]. Speaker data demonstrated that the amplitude of the broadband noise of [ç] was weaker than that of [ɕ]; the characteristic peak amplitude at approximately 4 kHz was greater in [ç] than in [ɕ], although the mid-sagittal vocal tract profiles were nearly identical for three of ten subjects in the rtMRI observation. These acoustic differences were reproduced by the proposed models, with differences in the width of the coronal plane constriction and the flow rate. The results suggest the need to include constriction width and flow rate as parameters for articulatory phonetic descriptions of speech sounds.

Global numerical simulation of fluid-structure-acoustic interaction in a single-reed instrument.

A numerical simulation of a single-reed instrument with a pressure chamber is conducted to examine the interaction among the flow, reed oscillation, and acoustic propagation. The flow and acoustic fields are predicted using the three-dimensional compressible Navier-Stokes equations, whereas the one-dimensional dynamic beam equation is solved for reed oscillation. The deforming geometry in the aeroacoustic field is expressed by the volume penalization method as an immersed boundary technique. The results showed that the waveforms of the tip opening and far-field acoustic spectra agreed well with those measured experimentally. The three-dimensional flow configuration near the tip opening was visualized, and the measurement of the instantaneous volume flow rate at the tip opening revealed that 30%-40% of the total flow rate passed through the side opening. The spectral tendencies of the time derivatives of the flow rate for different tip openings were consistent with that of the far-field sound, indicating that the slope of the flow rate waveform significantly affects the generated sound's harmonics.

Analysis of jet oscillations with acoustic radiation in the recorder by direct aeroacoustic simulations.

To elucidate the selection mechanism of a predominant mode in acoustic radiation from the recorder, the jet oscillations predicted by direct aeroacoustic simulations are analyzed based on a proposed formula for hydrodynamic and acoustic jet displacements. The displacements are well represented by the formula, taking into account the non-zero initial amplitude around the windway exit and variations of oscillation center with streamwise position for jet displacement. The analysis is applied to the jet oscillations in two different recorders with a straight- and an arch-shaped windway, where the shift of the predominant mode from the first to the second mode occurs for the straight-shaped recorder at a lower jet velocity compared to the arch-shaped recorder. The analytical results present the influence of the recorder shape on the amplification rate of the hydrodynamic jet displacement, the acoustic feedback effects and the phase relation between the hydrodynamic jet displacement and acoustic pressure in the resonator, along with the jet offset to the edge. Compared to the arch-shaped recorder, the convex curve of the amplification rate with the non-dimensional frequency based on the windway height for the straight-shaped recorder is located in the higher-frequency region, which contributes to the predominance of the higher mode.

Direct numerical simulation of fluid-acoustic interactions in a recorder with tone holes.

To clarify fluid-acoustic interactions in an actual recorder with opened and closed tone holes, flow and acoustic fields were directly numerically simulated on the basis of the compressible Navier-Stokes equations. To validate the simulation accuracy, the flow field around the windway and sound pressure above the window were measured. The predicted acoustic fields clarify changes of the positions of pressure nodes and anti-nodes in accordance with the state of the tone holes and the Mach number of the jet velocity. The fundamental mechanism of the self-sustained oscillations in a three-dimensional actual recorder is visualized by the predicted acoustic and flow fields. This result is also consistent with the relationship between the jet behaviors and pressure fluctuations based on the jet-drive model. Moreover, the effects of the fine vortices in the jet, which appear at the high Mach number of jet velocity of 0.099, on the sound are discussed. The time difference between the generation of the disturbances and the most intense deflection of the jet is identified and compared with the time delay of acoustic reflection around the window.