Prerpints.org logo
Preprint
Article

Measuring Absorbed Energy in the Human Auditory System Using Finite Element Models: A Comparison with Experimental Results

Altmetrics

Downloads

134

Views

79

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

18 October 2023

Posted:

19 October 2023

Read the latest preprint version here

Alerts
Abstract
here are different ways to analyse energy absorbance (EA) in the human auditory system. In previous research, we developed a complete finite element model (FEM) of the human auditory system. In this mentioned work, the external auditory canal (EAC), middle ear, and inner ear (spiral cochlea, vestibule, and semi-circular canals) were modelled based on human temporal bone histological sections. Multiple acoustic, structure and fluid-coupled analyses were conducted using the FEM to perform harmonic analyses in the 0.1–10 kHz range. Once the FEM had been validated with published experimental data, the FEM numerical results were used to calculate the EA or energy reflected (ER) by the tympanic membrane. This EA was also measured in clinical audiology tests which were used as a diagnostic parameter. A mathematical approach was developed to calculate the EA and ER, with numerical and experimental results showing adequate correlation up to 1 kHz. Another published FEM had adapted its boundary conditions to replicate experimental results. Here, we recalculated those numerical results by applying the natural boundary conditions of human hearing and found that the results almost totally agreed with our FEM. This boundary problem is frequent and problematic in experimental hearing test protocols: the more invasive they are, the more the results are affected. One of the main objectives of using FEMs is to explore how the experimental test conditions influence the results. Further work will still be required to uncover the relationship between the middle ear structure and EA to clarify how to best use FEMs. Moreover, the FEM boundary conditions must be more representative in future work to ensure their adequate interpretation.
Keywords: 
Subject: Engineering  -   Bioengineering

1. Introduction

The human auditory system (HAS) is an essential tool in performing numerous daily tasks. Therefore, it is important to be able to assess the proper functioning of the HAS in the most efficient way possible. The HAS is more complex than one might think at first glance, and there are still mechanisms whose exact functioning is unknown, or for which there are only hypotheses, such as the process of feedback from the inner and outer hair cells in the Basilar Membrane. In addition to this, there is the difficulty of accessing the HAS, which is embedded in the skull. Any intervention beyond the eardrum requires perforation or opening of the skull, with the corresponding consequences. These special circumstances have led medicine and engineering to develop non-invasive methods to determine the correct functioning of the HAS and, in case of malfunction, detect any associated pathology. The first and simplest test is audiometry, which involves subjecting the patient to the identification of beeps at different frequencies, increasing the intensity until surpassing the hearing threshold, at which point the patient informs the technician of the beep. Audiometry determines our hearing capacity in terms of frequency and decibels. Despite its simplicity, this test requires the patient’s cooperation and is not applicable to babies or individuals who are not conscious at the time. The measurement of impedance and/or Absorbed or Reflected Energy has been presented as an optimal solution.
Moller [1] first introduced the use of calibrated sound sources as a means for measuring impedances. The technique involves introducing a transmitter and receiver of waves into the auditory canal, so that the difference in energy between the emission and reception is evaluated. This difference can be related to the impedances of the different subsystems of the middle and inner ear, allowing for the inference of possible pathologies. This method, with some modifications or performances, is now extensively used. The following research [2,3,4,5] establishes an identical foundation in terms of the pressure source and receiver in the auditory canal. The difference between them lies in the technique used to approximate the calculation of eardrum impedance. Some of the studies from this time [6,7] begin to relate variations in impedance and, consequently, absorbed and reflected energies, to potential hearing loss.
In recent years, numerous studies have been conducted to improve the measurement of acoustic impedances [8,9,10,11,12], introducing Wideband Tympanometry (WBT) to replace the original Tympanometry. Nowadays, WBT is widely accepted and used for the calculation of impedances and reflected and absorbed energies. To do this, a passive test is performed that does not require any action from the patient. Some of these studies are more descriptive of the process and its proper setup [8,9], while others focus more on the relationship between the results and potential pathologies and hearing loss [10,11,12]. Norton or Thevenin equivalent sources are used for measuring device calibration, with only a single pressure measurement then being needed to calculate the external auditory canal (EAC) entrance impedance [13]. Other research has used the transfer function method, which relies on measuring the influence of the termination impedance on duct pressure [14], thereby requiring robust coupling of the device and the EAC [4,6]. Lanoye et al. [15] proposed a third method using impedance probes containing separated pressure and volume velocity sensors, although high levels of disturbance of the sound field at the measuring position are a problem when using this technique.
Most studies published to date have employed the two first methods where a sound (pressure) source and measuring device are usually placed inside the EAC. This is to avoid discontinuity between the device and the EAC to avoid exciting higher-order modes [16,17,18,19]. The frequencies of the minimum and maximum input impedances are mainly affected by the EAC length and its cross-sectional area. This means that the sound source and measurement device directly affect the impedance calculations because they change the canal length. One possible solution to this problem could be the application of inverse procedures to derive the EAC shape from its input impedance [12]. Thus, knowing the calculated cross-sectional area of the EAC allows the required eardrum impedance transformation for the energy calculations to be estimated.
To obtain middle ear diagnostics, the input impedances of the EAC measurements are transformed to represent the position of the eardrum and thus, record reflected energy (RE) and the energy absorbed, abbreviated as EA [5,7,8,13,14,15,16,17,18]. However, these impedance and energy results are often inaccurate. Thus, the main objective of this current work was to clarify the origin of possible errors in these measurements and their transformation as well as the consequences of these errors in terms of the impedances and energy calculated from them. In this sense, the use of finite element models (FEMs) is now widely accepted as an adequate and complementary research tool. Thus, this paper was based on several numerical simulations conducted using a FEM previously described and validated in the academic literature. This original work was first conducted for the outer and middle ear [20,21], was then applied to an inner ear FEM by employing a semiautomatic algorithm [22] and was finally coupled to the previous outer and middle ear FEMs [23]. Therefore, this FEM presents a complete fluid–structure interaction between the EAC, tympanic membrane (TM), and the oval window.
TM modelling was a crucial innovation with respect to previous models because the elements it used were better formulated and thereby eliminated problems associated with ‘shear locking’ elements used in thin membranes. Thus, together with adequate mesh convergence analysis, TM modelling provides sufficient guarantees of correct results. This means that, apart from geometric uncertainties caused by natural variability, most inaccuracy in this type of modelling comes from the difficulty of discerning the mechanical properties of some components such as the TM, tensors, and joints, among others. In this current article, we used assumed values for these components without attempting to discuss their accuracy. We simulated different combinations to discern the impact of each subsystem in the human auditory system (AS), building on a previous paper with a similar methodology [24] that determined how the AS influences pressure distribution in the EAC.
In this context, over the past century, Rong Gan has led a research group that has published several FEMs that have become an important source of inspiration for other researchers developing FEMs [25,26,27,28,29]. Gan’s work has focused on calculating EAC and eardrum impedances as well as the EA and ER from the TM [29]. Our work differs from that of Gan et al. in two clear ways. First, their main objective was to calculate the relationship between three middle ear disorders (otitis media, otosclerosis, and ossicular chain disarticulation), as simulated in their FEM, and any changes in the EA. Second, in contrast, our objective was to determine how the AS subsystems affect EA. We employed a range of FEMs, with the most basic one comprising an EAC and eardrum, to calculate the impedances and EAs. This methodology has already been applied successfully in previous work which gave us a better understanding of the mutual influence of each part of the AS on these factors.
It is important to understand the differences in our strategic goals and those of Gan et al. because these differing objectives affect the simulation conditions employed. Gan and colleagues aimed to offer a useful tool for the diagnosis of middle ear disorders by testing EAs through pneumatic otoscopic and wideband absorbance audiometry tests. Therefore, they literally reproduced the experimental setup boundary conditions and placed the sound source 20 mm from the eardrum to simulate the experimental test conditions in the FEM (see Figure 1). They determined that this approach gave them the best correlation between their numerical and experimental results. However, this does not mean that these numerical results most closely resemble reality. Indeed, controversy remains regarding the accuracy of the experimental methods used to measure impedance and in turn, the subsequent energy calculation results [30,31,32,33,34]. Therefore, more theoretical background work will still be required to increase our knowledge in this area.

2. Theoretical Background

The EA is calculated based on the ER. It describes the fraction of incident acoustic power reflected by the TM, where a reflective power of 1 corresponds to complete reflection of all the acoustic power and a reflectance of 0 corresponds to the condition in which all the power is absorbed by the TM [35]. First, the characteristic impedance of the EAC was calculated as:
Z c = ρ c S
where ρ is the density of the air contained in the EAC, c is the speed of sound in air, and S is the cross-sectional area of the EAC. In this study, air density and sound speed values of 1.21 kg/m3 and 343 m/s were used, respectively. To calculate the TM impedance, the impedance of each of the elements of the FEM for the TM was determined as a function of the acoustic pressure p T M i ( f ) , with v T M i f   being the velocity in element i of the TM, and where A T M i was the element area. Both the velocity and pressure depend on the frequency f, so we used values obtained for a frequency range of 0.1–10 kHz according to the FEM:
Z T M i f = p T M i ( f ) A T M i v T M i ( f )
The total TM impedance, Z T M f , was obtained by adding the impedances of each of the elements of the TM surface in parallel:
Z T M f = 1 ∑ i 1 Z T M i f
The acoustic impedance in the EAC was calculated based on Z T M f and Z c :
Z E C f = Z c Z T M f + j Z c tan ⁡ ( k L ) Z c + j Z T M f tan ⁡ ( k L )
where k is the wave number and L is the distance between the TM and the location of the measurement points in the EAC, which was at 30 mm in our study. The reflected acoustic pressure is obtained using the expression:
R P f = Z E C − Z C Z E C + Z C
Thus, the ER is calculated based on the reflected acoustic pressure R P f :
E R f = R P f 2 = Z E C − Z C Z E C + Z C 2
Finally, the EA was obtained as a function of the frequency calculated as:
E A f = 1 − E R f

3. Materials Models

Finite Element Model

The construction and validation of the FEM applied in this section, as well as its material properties, have been previously published [20,21,22,23,24]. However, in this current work, the EAC was re-meshed with hexahedral elements [36,37] to help improve the accuracy of the FEM post-process calculations (Figure 2).
As summarised in Table 1, three different FEM combinations were simulated to discern the impact of the different AS subsystems.

4. Results

The TM was modelled using 7880 hexahedral elements, with each one itself comprising 8 nodes, although only the 4 that formed the surface in contact with EAC air were valid in our calculations. Thus, a total of 8193 nodes formed the TM surface. With harmonic analysis, we obtained 16,386 real and imaginary values each for pressure and speed. Thus, when we performed the harmonic calculation for the 0.1–10 kHz range, we produced matrices comprising a total of 1,638,600 data points each for pressure and speed. In turn, the area was composed of only the 7880 element data points. We then used ANSYS 2022R2 engineering simulation software to obtain a total of 3,285,080 values for use in impedance and EA calculations using MATLAB R2023A software.

4.1. Tympanic membrane velocity

The average of all the elements in the TM velocity module is shown in Figure 3. There were three clearly differentiated resonance zones. One was at around 800–1000 Hz because of the eardrum itself, while the other two, located at 4 kHz and 9 kHz, were the result of the characteristic resonances of the EAC [24]. These results highlight the importance of modelling the ossicular chain, as there is a significant difference in the outcomes when the ossicular chain and cochlea are not modelled. On the other hand, there are no significant differences observed between modelling the simplified cochlea and the realistic spiral model, resulting in significant computational savings.

4.2. External auditory canal characteristic impedance

The characteristic impedance of the EAC, Zc depends on the air density ( ρ ), speed of sound in air (c), and cross-sectional area of the canal (S). The cross-sectional area was measured at the entrance of the canal where the pressure source was located and was 89,652 mm2. The density and speed of sound values were 1.21 kg/m3 and 343 m/s, respectively, and so the characteristic impedance of the EAC was 4.6294 × 103 Pa s/mm3. This value was used for all the FEMs and was frequency independent.

4.3. Tympanic membrane impedance (direct calculation)

Figure 4 represents the module and phase TM impedance in a frequency range of 0.1 to 10 kHz. Figure 4A shows a comparison of the module impedance for the three studied FEMs and showed a decrease in the impedance with frequency up to 800–1000 Hz, after which it increased. This minimum coincided with the first eardrum resonance frequency in that range. Figure 4B shows how the system formed by the EAC and TM opposes the system of least resistance to wave propagation and so was where the highest speeds occurred. The phase impedance started at a value of −80° for 100 Hz and increased with the frequency until it reached 0° for a range of 700–900 Hz. For frequency values between 8000–9000 Hz, the phase impedance reached a maximum between 60–80°. Therefore, the phase advanced until it reached the resonance frequency at around 900 Hz, with the phase also changing from lagging to leading at this point. In addition, there was a tendency towards an asymptotic phase shift at 180°, meaning that this was a second order system.

4.4. External auditory canal impedance (backward calculation) based on tympanic membrane impedance.

Figure 5A shows the module EAC impedance for each of the studied FEMs, which strongly coincided, especially at high frequencies. The EAC curves presented two minima, one at 3000 Hz and the other at 9000 Hz, with a maximum at around 6000 Hz. In turn, Figure 5B shows the phase EAC impedance curves for each FEM, which also strongly coincided at high frequencies. The TM stopped influencing the results in both the module and phase EAC impedance at around 1000 Hz [40,41], with curve disturbances at a lower range than those present at 3000, 6000, or 9000 Hz. The first and second resonance frequencies of the EAC were at 3000 and 9000 Hz, respectively, while the anti-resonance of the canal was at 6000 Hz.

4.5. Energy absorbance at the tympanic membrane

The results of EA for our models are presented In Figure 6, the EA curve presented a maximum for a frequency value around 700–900 Hz in all three tested FEMs, reaching a minimum for low frequency values between 100–200 Hz and a maximum at 9000–10,000 Hz. Figure 6 shows how the maximum EA values coincided with the minimum of the TM impedance, as shown in Figure 4B.
In Figure 7 are shown the EA results published. Zhang [29] presented results from numerical simulations by FEM. Other results [5,35,43,44] belongs to experimental tests. These experimental tests were carried out placing the pressure source inside the EAC, changing the natural EAC boundary conditions. There are differences between the experimental results, but It could be separate in two groups: on the one hand, Feeney [45] has the maximum before 1 kHz and fast decay; on the other hand, rest of experimental results have the maximum at 4 kHz. Our FEM results are more coincident with those presented by Feeney [45]

5. Discussion

To calculate the impedances, it appears that a well-defined cochlea model is unnecessary, especially at high frequencies. To validate this finding, we compared our data with those provided in previous academic publications. As shown in Figure 8 for module (A) and phase (B) EAC impedance, the values obtained with our FEMs corresponded to the data obtained by Zhang and Gan [36]. The blue curve represents the impedance of the EAC calculated according to Zhang and Gan to which a correction method [42] was applied to simulate our test condition (L = 0.03 m) and our EAC section (S), which thereby modified the value of the characteristic impedance (Zc). These results strongly concurred with the complete FEM model (CATIOSCO) for this correction.

6. Conclusions

Here we established a process for the numerical calculation of impedances and acoustic EA by the TM. This process starts with post-processing in ANSYS software followed by data exportation and processing in MATLAB. Regarding the impedance of the EAC, it is worth noting that the results obtained in this study agreed with those published by Gan and colleagues once the length of the canal was modified to 0.03 m, with values ranging from 10−4 to 10−1 Pa·s/mm³, valleys at 3 and 9 kHz, and a peak at 7 kHz. Modelling of the ossicular chain and cochlea were shown to be crucial to determine the impedance of the TM. Without these elements, the measured impedance dropped to 0.6·10−2 Pa·s/mm³ with a minimum around 800 Hz, while modelling the complete system reduced the minimum to 2·10−2 Pa·s/mm³, shifting the valley frequency to approximately 1 kHz. Of note, the experimental values previously reported [2,5,36] closely matched those obtained in this current study.
In this study were significantly different when the ossicular chain and cochlea were introduce ed into the models. Without these elements, EA reached a value of almost 1 at around 900 Hz. This is quite logical because this frequency coincides with the first natural frequency of the eardrum at which almost all incoming energy is absorbed. When the ossicular chain and cochlea were modelled, this maximum value reduced to 0.6 and the frequency slightly increased to 1 kHz. Nonetheless, experimental results showed discrepancies in these figures. Previous work [46] exhibited a similar maximum EA frequency but with a value close to 1 while other results [5,23,25,26,36] showed that the first part was nearly identical to our findings in this current study, with a maximum around 1 kHz and a value of 0.6, except that the EA continued to increase to 3–4 kHz, with values between 0.7 and 1.
Considering the equation used to calculate EA (Equation 7) and referring to its root equations (Equations 4– 6) involving ZTM, ZC, and ZEC (the TM and EAC impedances), we first showed that the ZTM we obtained coincided with the findings from experimental work presented elsewhere. Thus, we deduced that ZEC caused the differences between the modelled and experimental results mentioned above. This was because placement of the measuring device in the EAC considerably reduced its length, thereby affecting these impedance values, as demonstrated in Figure 8. This may be because it is difficult to establish the conditions of a ‘normal ear’ as a FEM [23].

Author Contributions

A.G.-G. and C.C.-E. developed all the FEMs in ANSYS software. J.A.A.-G. and P.M.C. carried out all the post-processing in MATLAB software. A.G.-H. was the project coordinator and supervised all the processes. All authors have read and agreed to the published version of the manuscript.

Funding

Not applicable.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Raw data of this article are available upon request from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Møller, A.R. Improved Technique for Detailed Measurements of the Middle Ear Impedance. J. Acoust. Soc. Am. 1960, 32, 250–257. [Google Scholar] [CrossRef]
  2. Rabinowitz, W.M. Measurement of the acoustic input immittance of the human ear. J. Acoust. Soc. Am. 1981, 70, 1025–1035. [Google Scholar] [CrossRef] [PubMed]
  3. Larson, V.D.; Nelson, J.A.; Cooper, W.A.; Egolf, D.P. Measurements of acoustic impedance at the input to the occluded ear canal. J. Rehabil. Res. Dev. 1993, 30, 129. [Google Scholar]
  4. Keefe, D.H.; Ling, R.; Bulen, J.C. Method to measure acoustic impedance and reflection coefficient. J. Acoust. Soc. Am. 1992, 91, 470–485. [Google Scholar] [CrossRef]
  5. Voss, S.E.; Allen, J.B. Measurement of acoustic impedance and reflectance in the human ear canal. J. Acoust. Soc. Am. 1994, 95, 372–384. [Google Scholar] [CrossRef] [PubMed]
  6. Sanborn, P.E. Predicting hearing aid response in real ears. J. Acoust. Soc. Am. 1998, 103, 3407–3417. [Google Scholar] [CrossRef] [PubMed]
  7. Shahnaz, N.; Bork, K.; Polka, L.; Longridge, N.; Bell, D.; Westerberg, B.D. Energy Reflectance and Tympanometry in Normal and Otosclerotic Ears. Ear Hear. 2009, 30, 219–233. [Google Scholar] [CrossRef]
  8. Jønsson, S.; Schuhmacher, A.; Ingerslev, H. Wideband impedance measurement in the human ear canal; In vivo study on 32 subjects. Brüel & Kjær Sound & Vibration A/S, Skodsborgvej 307, DK-2850, Denmark 2018, arXiv:1811.03389.
  9. Niemczyk, E.; Lachowska, M.; Tataj, E.; Kurczak, K.; Niemczyk, K. Wideband acoustic immitance—Absorbance measurements in ears after stapes surgery. Auris Nasus Larynx 2020, 47, 909–923. [Google Scholar] [CrossRef]
  10. Chris, S.; Lisa, H.; Patrick, F.; Hideko, H. Wideband Acoustic Immittance: Tympanometric Measures. Ear Hear. 2013, 34, 65s–71s. [Google Scholar]
  11. Ellison, J.C.; Gorga, M.; Cohn, E.; Fitzpatrick, D.; Sanford, C.A.; Keefe, D.H. Wideband acoustic transfer functions predict middle-ear effusion. Laryngoscope 2012, 122, 887–894. [Google Scholar] [CrossRef]
  12. Keefe, D.H.; Sanford, C.A.; Ellison, J.C.; Fitzpatrick, D.F.; Gorga, M.P. Wideband aural acoustic absorbance predicts conductive hearing loss in children. Int. J. Audiol. 2012, 51, 880–891. [Google Scholar] [CrossRef] [PubMed]
  13. Lawton, B.W.; Shaw, E.A.G. Estimation of acoustical energy reflectance at the eardrum from measurements of pressure distribution in the human ear canal. J. Acoust. Soc. Am. 1982, 72, 766–773. [Google Scholar]
  14. Utsuno, H.; Tanaka, T.; Fujikawa, T.; Seybert, A.F. Transfer function method for measuring characteristic impedance and propagation constant of porous materials. J. Acoust. Soc. Am. 1989, 86, 637–643. [Google Scholar] [CrossRef]
  15. Lanoye, R.; Vermeir, G.; Lauriks, W.; Kruse, R.; Mellert, V. Measuring the free field acoustic impedance and absorption coefficient of sound absorbing materials with a combined particle velocity-pressure sensor. J. Acoust. Soc. Am. 2006, 119, 2826–2831. [Google Scholar] [CrossRef]
  16. Hudde, H.; Letens, U. Scattering matrix of a discontinuity with a nonrigid wall in a lossless circular duct. J. Acoust. Soc. Am. 1985, 78, 1826–1837. [Google Scholar] [CrossRef]
  17. Fletcher, N.H.; Smith, J.; Tarnopolsky, A.Z.; Wolfe, J. Acoustic impedance measurements—Correction for probe geometry mismatch. J. Acoust. Soc. Am. 2005, 117, 2889–2895. [Google Scholar] [CrossRef]
  18. Stinson, M.R.; Daigle, G.A. Transverse pressure distributions in a simple model ear canal occluded by a hearing aid test fixture. J. Acoust. Soc. Am. 2007, 121, 3689–3702. [Google Scholar] [CrossRef]
  19. Brass, D.; Locke, A. The effect of the evanescent wave upon acoustic measurements in the human ear canal. J. Acoust. Soc. Am. 1997, 101, 2164–2175. [Google Scholar] [CrossRef]
  20. Garcia-Gonzalez, A.; Castro-Egler, C.; Gonzalez-Herrera, A. Analysis of the mechano-acoustic influence of the tympanic cavity in the auditory system. Biomed. Eng. Online 2016, 15, 33. [Google Scholar] [CrossRef]
  21. Garcia-Gonzalez, A.; Gonzalez-Herrera, A. Effect of the middle ear cavity on the response of the human auditory system. J. Acoust. Soc. Am. 2013, 133. [Google Scholar] [CrossRef]
  22. Castro-Egler, C.; García-González, A. Semiautomatic algorithm for 3D modelling of finite elements of the cochlea. J. Med. Imaging Health Inform. 2017, 7, 1049–1056. [Google Scholar] [CrossRef]
  23. García-González, A.; Durán-Escalante, A.; Castro-Egler, C. 3D Modelling and Numerical Analysis of Human Inner Ear by Means of Finite Elements Method. Int. J. Comput. Inf. Eng. 2016, 10, 596–606. [Google Scholar]
  24. Garcia-Gonzalez, A.; Castro-Egler, C.; Gonzalez-Herrera, A. Influence of the auditory system on pressure distribution in the ear canal. J. Mech. Med. Biol. 2018, 18, 1850021. [Google Scholar] [CrossRef]
  25. Gan, R.Z.; Reeves, B.P.; Wang, X. Modeling of Sound Transmission from Ear Canal to Cochlea. Ann. Biomed. Eng. 2007, 14, 35. [Google Scholar] [CrossRef] [PubMed]
  26. Gan, R.Z.; Sun, Q.; Feng, B.; Wood, M.W. Acoustic–structural coupled finite element analysis for sound transmission in human ear—Pressure distributions. Med. Eng. Phys. 2006, 28, 395–404. [Google Scholar] [CrossRef] [PubMed]
  27. Gan, R.Z.; Feng, B.; Sun, Q. Three-Dimensional Finite Element Modeling of Human Ear for Sound Transmission. Ann. Biomed. Eng. 2004, 32, 847–859. [Google Scholar] [CrossRef] [PubMed]
  28. Brown, M.A.; Bradshaw, J.J.; Gan, R.Z. Three-Dimensional Finite Element Modeling of Blast Wave Transmission From the External Ear to a Spiral Cochlea. J. Biomech. Eng. 2022, 144, 014503. [Google Scholar] [CrossRef] [PubMed]
  29. Zhang, X.; Gan, R.Z. Finite element modeling of energy absorbance in normal and disordered human ears. Hear Res. 2013, 301, 146–155. [Google Scholar] [CrossRef] [PubMed]
  30. Nørgaard, K.R.; Fernandez-Grande, E.; Laugesen, S. Compensating for oblique ear-probe insertions in ear-canal reflectance measurements. J. Acoust. Soc. Am. 2019, 145, 3499–3509. [Google Scholar] [CrossRef]
  31. Hudde, H.; Engel, A. Measuring and Modeling Basic Properties of the Human Middle Ear and Ear Canal. Part III: Eardrum Impedances, Transfer Functions and Model Calculations. In Acta Acustica United with Acustica; S. Hirzel Verlag: Stuttgart, Germany, 1998; Volume 84. [Google Scholar]
  32. Hudde, H.; Engel, A. Measuring and Modeling Basic Properties of the Human Middle Ear and Ear Canal. Part II: Ear Canal, Middle Ear Cavities, Eardrum, and Ossicles. In Acta Acustica United with Acustica; S. Hirzel Verlag: Stuttgart, Germany, 1998; Volume 84. [Google Scholar]
  33. Hudde, H.; Engel, A. Measuring and Modeling Basic Properties of the Human Middle Ear and Ear Canal. Part I: Model Structure and Measuring Techniques. In Acta Acustica United with Acustica; S. Hirzel Verlag: Stuttgart, Germany, 1998; Volume 84. [Google Scholar]
  34. Hudde, H. Measurement of the eardrum impedance of human ears. J. Acoust. Soc. Am. 1983, 73, 242–247. [Google Scholar] [CrossRef] [PubMed]
  35. Rosowski, J.J.; Nakajima, H.H.; Hamade, M.A.; Mahfoud, L.; Merchant, G.R.; Halpin, C.F. ; Ear-Canal Reflectance, Umbo Velocity, and Tympanometry in Normal-Hearing Adults. Ear Hear. 2012, 33, 19–34. [Google Scholar] [CrossRef] [PubMed]
  36. Ramos Valero, J.P. Modelado y Análisis Numérico Con el Método de los Elementos Finitos del Canal Auditivo Externo y Membrana Timpánica; Riuma, Málaga, Spain, 2013.
  37. Castro Egler, C. Modelado 3D & Análisis Numérico del Oído Interno Mediante Elementos Finitos. 2020. Available online: https://hdl.handle.net/10630/19420 (accessed on 28 July 2022).
  38. Stepp, C.E.; Voss, S.E. Acoustics of the human middle-ear air space. J. Acoust. Soc. Am 2005, 118, 861–871. [Google Scholar] [CrossRef]
  39. Voss, S.E.; Rosowski, J.J.; Merchant, S; N., *!!! REPLACE !!!*; Peake, W.T. Acoustic responses of the human middle ear. Hear. Res. 2000, 150, 43–69. [Google Scholar] [CrossRef] [PubMed]
  40. Rasetshwane, D.M.; Neely, S.T. Inverse solution of ear-canal area function from reflectance. J. Acoust. Soc. Am. 2011, 130, 3873–3881. [Google Scholar] [CrossRef]
  41. Allen, J.B.; Jeng, P.S.; Levitt, H. Evaluation of human middle ear function via an acoustic power assessment. J. Rehabil. Res. Dev. 2005, 42, 63–78. [Google Scholar] [CrossRef]
  42. Sebothoma, B.; Khoza-Shangase, K.; Mol, D.; Masege, D. The sensitivity and specificity of wideband absorbance measure in identifying pathologic middle ears in adults living with HIV. South Afr. J. Commun. Disord. 2021, 30, 68. [Google Scholar] [CrossRef]
  43. Piskorski, P.; Keefe, D.H.; Simmons, J.L.; Gorga, M.P. Prediction of conductive hearing loss based on acoustic ear-canal response using a multivariate clinical decision theory. J. Acoust. Soc. Am. 1999, 105, 1749–1764. [Google Scholar] [CrossRef] [PubMed]
  44. Keefe, D.H.; Simmons, J.L. Energy transmittance predicts conductive hearing loss in older children and adults. J. Acoust. Soc. Am. 2003, 114, 3217–3238. [Google Scholar] [CrossRef]
  45. Feeney, M.P.; Grant, I.L.; Marryott, L.P. Wideband energy reflectance measurements in adults with middle-ear disorders. Journal of Speech, Language, and Hearing Research 2003, 46, 901–911. [Google Scholar] [CrossRef]
  46. Yilmaz, N.; Soylemez, E.; Sanuc, M.B.; Bayrak, M.H.; Sener, V. Sound energy absorbance changes in the elderly with presbycusis with normal outer and middle ear. Eur. Arch. Oto-Rhino-Laryngol. 2023, 280, 2265–2271. [Google Scholar] [CrossRef]
Preprints 88122 g001
Figure 1. A cross-section of the ear illustrating the probe used to measure impedance when inserted into the ear canal. The probe emits a sound pressure wave that is incident to the tympanic membrane. Some of the incident sound pressure is reflected and this is then measured by the probe. The remaining incident sound pressure is absorbed by the tympanic membrane and the structures behind it.
Figure 1. A cross-section of the ear illustrating the probe used to measure impedance when inserted into the ear canal. The probe emits a sound pressure wave that is incident to the tympanic membrane. Some of the incident sound pressure is reflected and this is then measured by the probe. The remaining incident sound pressure is absorbed by the tympanic membrane and the structures behind it.
Preprints 88122 g001
Preprints 88122 g002
Figure 2. Hexahedral finite element model (FEM) for the external auditory canal based on the tympanic membrane FEM.
Figure 2. Hexahedral finite element model (FEM) for the external auditory canal based on the tympanic membrane FEM.
Preprints 88122 g002
Preprints 88122 g003
Figure 3. The means for the tympanic membrane speed module for the three finite element models studied, CATIOSCO, CATI, and CATIOS. The CATI comprised the external auditory canal (EAC) and tympanic membrane (TM); CATIOS comprised the EAC, TM, ossicular chain (OC), and simplified cochlea; and CATIOSCO comprised the EAC, TM, OC, and the entire cochlea.
Figure 3. The means for the tympanic membrane speed module for the three finite element models studied, CATIOSCO, CATI, and CATIOS. The CATI comprised the external auditory canal (EAC) and tympanic membrane (TM); CATIOS comprised the EAC, TM, ossicular chain (OC), and simplified cochlea; and CATIOSCO comprised the EAC, TM, OC, and the entire cochlea.
Preprints 88122 g003
Preprints 88122 g004
Figure 4. Comparison of the three studied finite element models with experimental test [2,38,39] and FEM [29] for the module (A) and phase (B) tympanic membrane impedance (ZTM).
Figure 4. Comparison of the three studied finite element models with experimental test [2,38,39] and FEM [29] for the module (A) and phase (B) tympanic membrane impedance (ZTM).
Preprints 88122 g004
Preprints 88122 g005
Figure 5. Comparison of the three studied finite element models for the module (A) and phase (B) external auditory canal impedance (ZEC).
Figure 5. Comparison of the three studied finite element models for the module (A) and phase (B) external auditory canal impedance (ZEC).
Preprints 88122 g005
Preprints 88122 g006
Figure 6. Energy absorbed by the tympanic membrane for the different finite element models studied.
Figure 6. Energy absorbed by the tympanic membrane for the different finite element models studied.
Preprints 88122 g006
Preprints 88122 g007
Figure 7. Energy absorbed by the tympanic membrane for the finite element models and experimental results described in other publications [4,5,29,45,46].
Figure 7. Energy absorbed by the tympanic membrane for the finite element models and experimental results described in other publications [4,5,29,45,46].
Preprints 88122 g007
Preprints 88122 g008
Figure 8. Phase comparison of the module (A) and phase (B) external auditory canal impedance with the findings published by Zhang [29].
Figure 8. Phase comparison of the module (A) and phase (B) external auditory canal impedance with the findings published by Zhang [29].
Preprints 88122 g008
Table 1. Combinations of different finite elements simulated in the CATI, CATIOS, and CATIOSCO finite element models.
Table 1. Combinations of different finite elements simulated in the CATI, CATIOS, and CATIOSCO finite element models.
Model Subsystem Modelled Name
Preprints 88122 i001 External auditory canal
Tympanic membrane		 CATI
Preprints 88122 i002 External auditory canal
Tympanic membrane
Ossicular chain
Cochlea simplified		 CATIOS
Preprints 88122 i003 External auditory canal
Tympanic membrane
Ossicular chain
Cochlea
Vestibuli
Semicircular canals		 CATIOSCO
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated