# An Analytical Model for Bandwidth Enhancement of Air-Coupled Unsealed Helmholtz Structural CMUTs

Authors: Xin Liu; Yuanyu Yu; Jiujiang Wang; Sio Hang Pun; Mang I Vai; Peng Un Mak
Journal: Journal of Sensors
 (2019)

---

## Abstract

Capacitive micromachined ultrasonic transducers (CMUTs) were reported to own high potential in air-coupled ultrasonic applications such as noncontact nondestructive examination and gas flow measurement. The unsealed CMUTs which utilized the squeeze film effect were reported to overcome the narrow output pressure bandwidth of the conventional sealed CMUTs in air operation. This kind of unsealed CMUTs can also be regarded as Helmholtz resonators. In this work, we present the air-coupled unsealed Helmholtz structural CMUTs which utilize both the squeeze film effect and the Helmholtz resonant effect to enhance the output pressure bandwidth. Based on the mechanism of vibration coupling between membrane and air pistons in membrane holes, we propose an analytical model to aid the design process of this kind of CMUTs. We also use finite element method (FEM) to investigate this kind of CMUTs for our analytical model validation. The FEM results show that the significant bandwidth enhancement can be achieved when the Helmholtz resonant frequency is designed close to the fundamental resonant frequency of the CMUT membrane. Compared with the conventional sealed CMUT cell, the 4-hole unsealed Helmholtz structural CMUT cell improves both the 3-dB fractional bandwidth and SPL-bandwidth product around 35 times. Furthermore, it is found that, with more holes under the same hole area ratio or with a smaller ratio of the cavity height to the viscous boundary layer thickness, the Helmholtz resonant effect becomes weaker and thus the output pressure bandwidth decreases.


---

## Body

## 1. Introduction

Capacitive micromachined ultrasonic transducers (CMUTs) were invented as an alternative to conventional piezoelectric ultrasonic transducers in early 1990s. Compared with the piezoelectric ultrasonic transducers, CMUTs have many advantages such as they are suitable for batch production and can be easily integrated with modern electronics. Conventional piezoelectric ultrasonic transducers typically require matching layers to compensate the acoustic impedance mismatch between piezoelectric materials and surrounding gaseous or liquid media. CMUTs feature a better impedance matching with gaseous and liquid media, and therefore they can transmit/receive ultrasound more efficiently into gas and liquid without the matching layer [1].In many ultrasonic applications, the transducers were used in immersion [2, 3]. However, in some applications, some objects cannot contact with liquid couplants, such as paper products and electronic products [4]. Also there are many other air-borne ultrasonic applications such as air-coupled ultrasonic imaging [4], ultrasonic distance measurement [5], anthropometrical measurement of scoliosis [6], and gas flow measurement [7]. CMUTs ideally are suited for air-coupled applications because of their good acoustic impedance matching with air.For air-coupled ultrasonic applications, transducers with wide bandwidth have sharp ultrasonic pulses and short decay time of the vibration amplitude, thus providing high resolution and a small “dead zone” where objects cannot be detected [8]. Wide bandwidth air-coupled transducers are needed for applications such as echo profile evaluation and distance measurement with high resolution [8]. However, the conventional vacuum sealed air-coupled CMUTs still own a narrow bandwidth. For instance, for the CMUT array with only one cell radius, the 6-dB fractional bandwidth is about 1% [9]. To overcome this, unsealed CMUTs were proposed to achieve wider bandwidth by utilizing the squeeze film effects [10–13]. For unsealed CMUTs, the vented cavity not only introduces the squeeze film effects, but also constructs a Helmholtz resonator with the vented holes. The Helmholtz resonator comprised a cavity with an opening denoted as “neck”. When an excitation with a specific frequency range is applied, a volume of air in the neck would vibrate due to the elasticity of the air inside the cavity and thus radiates sound waves into surroundings. This is called Helmholtz resonant effect [14]. The Helmholtz resonant effect can also exert an influence on the device performance; this phenomenon was reported, but there was a lack of analytical explanation in [11]. In [15], the Helmholtz resonant effect was used to enlarge the bandwidth of the CMUT’s output, but, in its FEM simulation, the squeeze film damping was neglected. Actually, with proper adjustment of volume ratio between vented hole and cavity, the performance of the device would be benefited by the vibration coupling between the air pistons vibration system of Helmholtz resonator and the membrane vibration system of CMUT cell. This vibration coupling was widely used to enhance the low frequency response of musical instrument such as acoustic guitar [16] and vented-box loudspeaker systems [17].In this paper, we present an air-coupled unsealed CMUT cell with Helmholtz resonant cavity, which utilizes both the squeeze film effect and Helmholtz resonant effect to enhance the output pressure bandwidth of the device. For convenience, we hereafter call it as Helmholtz structural CMUT. We propose an analytical model to explain the vibration coupling working mechanism of this kind of CMUT and aid the design process of this kind of CMUT. In order to validate this analytical model and provide deeper investigation on this kind of Helmholtz structural CMUT, the FEM model is built with COMSOL Multiphysics version 5.2a (COMSOL Inc., Stockholm, Sweden) to simulate three cases of 4-hole CMUTs with different hole radii.


## 2. Methodologies

Figure1(a) shows an air-coupled unsealed Helmholtz structural CMUT cell with 4 holes in its membrane, the cross-sectional illustration of this CMUT cell along the red dashed line in (a) is shown in Figure 1(b). Major components of this CMUT cell are a circular membrane with radius R and thicknesst, a cavity with heighth, holes with radiusa, the thin metal layer as full top electrode, and one insulator layer with thicknesstins to avoid short circuit between electrodes. The bottom electrode is located between the insulator layer and the substrate. The rim of the membrane is clamped by support. Using Helmholtz resonant effect and squeeze film effect to enhance bandwidth, several circular holes need to be punched through the membrane. These holes and the cavity form a Helmholtz resonator. If a DC voltage is applied between top and bottom electrodes, the membrane will be activated by the electrostatic force induced by the voltage and the restoring spring force due to the membrane’s stiffness. When these two forces are equal, the membrane will reach an equilibrium position. Then when a small amplitude AC voltage is superimposed, the unbalance of these two forces will be caused by the small perturbation of the electrostatic force, which leads to the membrane’s vibration. The membrane’s vibration will drive the air pistons in the holes to vibrate due to the elastic property of the air inside the cavity. The vibration of the air pistons will also induce a pressure change in the cavity, which leads to a force acted on the membrane. Therefore, the membrane’s vibration and the air pistons’ vibration interact with each other. During the operation of this air-coupled unsealed Helmholtz structural CMUT cell, the ultrasonic wave will be generated by both vibrations of the membrane and the air pistons.Figure 1
(a) A 4-hole unsealed circular CMUT cell with Helmholtz resonator structure. (b) Cross-sectional diagram of this cell along dashed line of (a).
(a) (b)


### 2.1. Lumped-Parameter Analytical Model

As Figure2 shows, inspired by the simple model of acoustic guitar [16], the Helmholtz structural CMUT cell can be simplified as a Helmholtz resonator with a movable top plate. The Helmholtz resonator is formed from the cavity of CMUT cell with volume V and an equivalent air piston with mass ma and top surface areaS, where the top surface area S is the sum of top surface areas of all small holes. And we can simplify the membrane as a movable top plate with mass mp, top surface area A, and stiffness kp.Figure 2
Simplified model of Helmholtz structural CMUT cell.Then referencing from [16], we can develop a lumped-parameter analytical model of the air-coupled unsealed Helmholtz structural CMUT cell as follows.When the DC and AC bias voltages are superimposed, the top plate vibrates as a harmonic oscillator and drives the equivalent air piston to vibrate due to the elastic property of air cushion in the cavity. We assume that the deflected shape of the top plate can be neglected and the uniform displacement of the top plate is denoted asxp. The uniform displacement of the equivalent air piston is taken as xa. Then we take outward movement as positive; the volume change of the cavity is(1)ΔV=Axp+SxaThus the pressure change as adiabatic compression in cavity is(2)Δp=-μΔVwhereμ=c2ρ/V,c is sound velocity in air, and ρ is the density of air. The forces acted on plate and equivalent air piston due to the pressure change can be represented byAΔpandSΔp, and the equations of the vibrations of the top plate and the equivalent air piston can be written as(3a)mpd2xpdt2=F-Rpdxpdt-kp+μA2xp-μSAxa(3b)mad2xadt2=-Radxadt-μS2xa-μSAxpwhere F is the driving force of top plate, Ra and Rp are resistances to movement of the top plate and the equivalent air piston, respectively. From (3a) and (3b), we can find that the top plate’s vibration and the equivalent air piston’s vibration are coupled, because they both depend on xa and xp. Namely, the vibration of one of them can have effect on the other.Uncoupled Helmholtz resonator and uncoupled top plate: firstly, if we setS=0 in (3a) and xp=0 in (3b), we can get the equations which describe the vibrations of the uncoupled top plate and uncoupled Helmholtz resonator, respectively. (4a)mpd2xpdt2=F-Rpdxpdt-kp+μA2xp(4b)mad2xadt2=-Radxadt-μS2xaFrom (4a), we can get the resonant frequency of the uncoupled top plate as(5)ωp=kp+μA2mpReferencing from [18], this frequency can be calculated more accurately as the fundamental resonant frequency of a clamped circular membrane as(6)ωp=2πfp=2π·0.47tR2Yρm1-σ2where ρm,σ, and Y are membrane material’s density, Poisson’s ratio and Young’s modulus, respectively. From (4b), the resonant frequency of an uncoupled Helmholtz resonator can be solved as (7)ωh=2πfh=μS2maWe substituteμ=c2ρ/V, ma=Slρ into (7), and we can get (8)ωh=cSVlAccording to [19], the length of the Helmholtz resonator’s neck l should take an effective valueleff=l+δ for better estimation; thus,(9)ωh=cSVl+δwhereδ is the end correction. Referencing from [19], the end correction δ should take the sum of the interior correction δi and the exterior correction δe as follows:(10)δ=δi+δe=0.48S1-1.25nholeaR+0.48Swherenhole is the number of holes in the membrane of the CMUT cell.Coupled Helmholtz resonator and top plate: now, we focus on the coupled vibration system, assuming in harmonic oscillation state, and we solve (3a) and (3b) in frequency domain, by writing in phasor format xp=Xpeiωt, xa=Xaeiωt, F=F0eiωt, the amplitude of top plate’s velocity Up=iωXp, and the amplitude of equivalent air piston’s velocity Ua=iωXa. Similar to [16], we can get (11a)Up=iωF0ωh2-ω2+iωγampD(11b)Ua=-iωF0Aωh2mpSDwhere(12)D=ωp2-ω2+iωγpωh2-ω2+iωγa-ωph4(13)γa=Rama,γp=Rpmp,ωph4=ωh2ωa2,ωa2=μA2mpωp and ωh are the resonant frequencies of uncoupled top plate and uncoupled Helmholtz resonator, respectively.From (11a) and (11b), the coupled vibration system shown in Figure 2 resonates while the moduli of Ua and Up arrive maximum (D=0). If we neglect the damping terms ofiωγa andiωγp in (12), then this coupled vibration system’s resonant frequencies can be solved as(14)ω±2=ωp2+ωh22±ωp2-ωh22+4ωph42As shown in (6) and (9), for the uncoupled top plate and the uncoupled Helmholtz resonator, each of them only has one resonant frequency. From (14) we can see that, unlike the uncoupled situation, both the top plate and the equivalent air piston have two resonant frequencies in the coupled vibration system.For an air-coupled unsealed Helmholtz structural CMUT cell, when calculatingωh by (9), the length of the neck l is equal to the membrane’s thickness t. And the lumped parameters are calculated as follows:(15)S=nholeπa2(16)A=πR2-S(17)V=πR2h(18)mp=1.8Atρm+mrp(19)ma=Stρ+mrawhere the factor 1.8 in (18) is from [20], mrp is the acoustic radiation mass of the plate, and mra in (19) is the acoustic radiation mass of the equivalent air piston. Referenced from [21, 22], the CMUT cell should be considered as a clamped plate source located on an infinite baffle, and the acoustic radiation mass of the plate is (20)mrp=-ρcAω20F22ω/cA/πω/cA/π9,F2x=-x4-91x2+504H1x-14xx2-18H0x+14x415π-168x2πwhere Hn(x) is the nth order Struve function. The equivalent air piston is considered as a baffled piston source; the radiation mass of the equivalent air piston is [14](21)mra=nholeρcπa2ω2H12ω/ca2ω/caIn the simplification in Figure 2, the driving force F is a harmonic force; hence, we take F as the harmonic electrostatic force, according to [1],(22)F0=Fac=ε0AVdcVach2where Vdc, Vac are DC voltage and AC voltage applied to the cell, and ε0 is the dielectric permittivity of air.The squeeze film effect includes squeeze film stiffening effect and squeeze film damping effect. The proportion of the stiffening effect and damping effect is determined by the squeeze number [20](23)σsqeeze=12μairωR2Pah2whereμair is the dynamic viscosity of air, Pa is the ambient pressure, and ω is the plate’s vibrating angular frequency. For the air-coupled unsealed CMUT with perforated membrane, the squeeze damping effect dominates the squeeze film effect because of the small squeeze number, and the squeeze stiffening effect can be neglected [20, 23]. From [20], the squeeze damping can be represented by the lump resistance as (24)Rsqeeze=3πμairR4κβ2h3nholewhere(25)κβ=4β2-β4-4ln⁡β-3(26)β=a1.05r1/2where r1 represents the distance between center of hole and center of membrane andnhole is the total number of holes.For the plate, although there are other damping mechanisms, the squeeze film damping is the most noticeable once it is presented in MEMS devices [23]. In addition, CMUT is an acoustic device; the loss of acoustic radiation cannot be ignored [10]. According to [21, 22], the resistance of the acoustic radiation of the clamped plate is(27)Rrad_plate=ρcA1-20F12ω/cA/πω/cA/π9,F1x=x4-91x2+504J1x+14xx2-18J0x-x516-x7768where Jn(x) is the Bessel function of the first kind of order n. Thus, the damping resistance of the top plate Rp is (28)Rp=Rsqueeze+Rrad_platewhereRsqueeze is calculated as in (24) to (26).Rrad_plate is calculated as in (27).For the equivalent air piston, not only the acoustic radiation damping but also the viscous damping needs to be considered. When the air pistons vibrate, due to the viscosity of the air, they are damped by the viscous force at the side walls of the holes. From [24], the resistance of viscous boundary layer per unit area can be expressed as (29)Rvisc=ρμairω2When calculating the viscous damping of the equivalent air piston, the area of viscous boundary layer is taken as the area of the side walls of the holes and a viscous end correction [19]. The resistance of acoustic radiationRrad_air piston is calculated as the baffled piston source [14].(30)Rrad_airpiston=nholeρcπa21-2J12ωa/c2ωa/cThus, the damping resistance of the equivalent air piston is (31)Ra=nhole2πat+4πa2Rvisc+Rrad_airpistonAfter gettingUa and Up, the complex amplitude of total volume velocity U can be written as(32)U=AUp+SUaFrom [21, 25], if we consider a CMUT cell as a clamped edge radiator source, the complex amplitude of sound pressure of the point at distance d in the far field above the center of the membrane can be calculated as (33)p=iωρU2πdthen the sound pressure level (SPL) at this point can be calculated.


### 2.2. Finite Element Method Model

In order to investigate the air-coupled unsealed Helmholtz structural CMUTs and validate the analytical model, we have developed the finite element method (FEM) model to simulate the device. Unsealed CMUTs would induce squeeze film effect; usually the squeeze film effect in MEMS devices is described by nonlinear compressible Reynolds equation, which is simplified from the continuity equation and the Navier-Stokes equation when the modified Reynolds number satisfies [13, 20](34)Remodified=ωρh2μair≪1This condition is not effective for CMUTs which have slightly higher gaps or higher operating frequencies. Thus, we need to start with the continuity equation and the Navier-Stokes equation to model the squeeze film effect in our unsealed Helmholtz structural CMUTs [13].The commercial software COMSOL Multiphysics version 5.2a is used to build our FEM model. The Electromechanical (emi) module and Thermoviscous Acoustics, Frequency Domain (ta) module are used in our FEM model of a CMUT cell. The governing equations ofta include the Navier-Stokes equation, the continuity equation, and the energy equation formulated using the Fourier heat law [26]. Thus, the squeeze film effect is included and the thermal loss is also taken into account.We use four symmetric distributed holes here to avoid the damage of the center part of membrane which makes main contribution to membrane’s vibration. Table1 lists the material properties in simulation. The 3D model of a quarter of a CMUT cell in octant spherical air domain is built, as shown in Figures 3(a) and 3(b). A “perfectly matched layer (PML)” is used to eliminate the acoustic reflection and to simulate an infinite spread medium. Inemi module, for the “domain selection”, it includes all domains. The membrane, the insulator layer, and the support are selected as “linear elastic dielectric”. The rim of the membrane, the support, and the insulator layer are applied with “fixed constraint” boundary condition. The substrate can be ignored because of the fixed insulator layer. For simplification in FEM, the top and bottom electrodes can be regarded as infinitely thin [27]; thus, we directly use the membrane’s top surface as the top electrode, and the bottom electrode is located underneath the insulator layer. A constant DC bias voltage and a harmonic AC voltage (amplitude as 1% of DC voltage) are applied on the top surface of membrane by the “terminal” with “harmonic perturbation” conditions, and the bottom surface of insulator layer is set as “ground”. The “symmetry” is applied on the vertical sections of the membrane, insulator layer, and support. Other boundary conditions keep the default setting in COMSOL. Inta module, for the “domain selection”, only domains filled with air (cavity, holes, spherical air domain, PML) are selected. To couple two modules, “velocity” and “isothermal” boundary conditions are applied inta on the top surface of membrane and inner wall of cavity and of hole, which are interfaces betweenta andemi [26]. The “velocity” inta is prescribed as inemi. This coupling setting is a two-way coupling, which means when the membrane vibrates to generate sound waves in medium, a reaction force would be added to the membrane [28]. Thus we do not need to add other coupling condition inemi module. The “symmetry” boundary condition is applied on the vertical sections of cavity and the vertical sections of octant spherical air and PML. Other boundary conditions keep the default setting in COMSOL. In the mesh setup, the mesh size should be sufficiently small to resolve the wavelength by at least 5-6 elements per wavelength [29], and the maximum element size of the acoustic medium is set as one-sixth of wavelength of operating frequency. Also the “boundary layers mesh” is needed at the inner walls of both cavity and holes; the thickness of “boundary layers mesh” is equal to the viscous boundary layer thickness under fundamental resonant frequency of CMUTs membrane [26]. The meshing details of vertical section of the cell are shown in Figure 3(c), where the “boundary layers meshes” are highlighted. The “prestressed analysis” is applied in calculating the average pressure over different frequencies. Firstly, a stationary study is used to calculate the membrane deflection under a certain DC bias voltage. Then an AC voltage is biased on the membrane using harmonic analysis over a frequency range. The pressure is calculated by averaging the pressure between the membrane and the air medium.Table 1
Material properties used in COMSOL (∗: the factor 0.6 is from [30]; others are the preset values/functions in COMSOL).
Parameter Value/Preset function in COMSOL Air: (used in cavity, holes, spherical air domain, PML) Relative permittivity 1 Dynamic viscosityµair [Pa.s] eta(T[1/K]) Heat capacity at constant pressureCp [J/(kg.K)] Cp(T[1/K]) Densityρair [kg/m3] rho(pA[1/Pa],T[1/K]) Thermal conductivityK [W/(m.K)] k(T[1/K]) ∗Bulk viscosity µb [Pa.s] 0.6∗eta(T[1/K]) Typical Wave Speed for PML [m/s] 343  Si: (used in membrane, support beam) Young’s modulusY [GPa] 170 Poisson ratioσ 0.28 DensityρSi [kg/m3] 2329 Relative permittivity 11.7  SiO 2 : (used in insulator layer) Young’s modulus [GPa] 70 Poisson ratio 0.17 Density [kg/m3] 2200 Relative permittivity 4.2Figure 3
A 3D model in COMSOL: (a) a quarter of a CMUT cell and spherical air domain; (b) details of the unsealed Helmholtz structural CMUT cell; (c) meshing details of the vertical section of (b).In order to make a comparison between the unsealed Helmholtz structural CMUTs and the conventional vacuum sealed CMUTs, we also build an FEM model for the latter. The model building process of the conventional cell is slightly different from the unsealed Helmholtz structural CMUT cell. The solid-air interfaces for the conventional cell only contain the top surface of the membrane; thus, the “velocity” and “isothermal” are only added on the membrane’s top surface. The cavity is vacuum; thus, for a conventional cell, we do not take the acoustic field of the cavity intota module’s calculation. Also, the “boundary layers mesh” is not needed for the vacuum cavity.


## 2.1. Lumped-Parameter Analytical Model

As Figure2 shows, inspired by the simple model of acoustic guitar [16], the Helmholtz structural CMUT cell can be simplified as a Helmholtz resonator with a movable top plate. The Helmholtz resonator is formed from the cavity of CMUT cell with volume V and an equivalent air piston with mass ma and top surface areaS, where the top surface area S is the sum of top surface areas of all small holes. And we can simplify the membrane as a movable top plate with mass mp, top surface area A, and stiffness kp.Figure 2
Simplified model of Helmholtz structural CMUT cell.Then referencing from [16], we can develop a lumped-parameter analytical model of the air-coupled unsealed Helmholtz structural CMUT cell as follows.When the DC and AC bias voltages are superimposed, the top plate vibrates as a harmonic oscillator and drives the equivalent air piston to vibrate due to the elastic property of air cushion in the cavity. We assume that the deflected shape of the top plate can be neglected and the uniform displacement of the top plate is denoted asxp. The uniform displacement of the equivalent air piston is taken as xa. Then we take outward movement as positive; the volume change of the cavity is(1)ΔV=Axp+SxaThus the pressure change as adiabatic compression in cavity is(2)Δp=-μΔVwhereμ=c2ρ/V,c is sound velocity in air, and ρ is the density of air. The forces acted on plate and equivalent air piston due to the pressure change can be represented byAΔpandSΔp, and the equations of the vibrations of the top plate and the equivalent air piston can be written as(3a)mpd2xpdt2=F-Rpdxpdt-kp+μA2xp-μSAxa(3b)mad2xadt2=-Radxadt-μS2xa-μSAxpwhere F is the driving force of top plate, Ra and Rp are resistances to movement of the top plate and the equivalent air piston, respectively. From (3a) and (3b), we can find that the top plate’s vibration and the equivalent air piston’s vibration are coupled, because they both depend on xa and xp. Namely, the vibration of one of them can have effect on the other.Uncoupled Helmholtz resonator and uncoupled top plate: firstly, if we setS=0 in (3a) and xp=0 in (3b), we can get the equations which describe the vibrations of the uncoupled top plate and uncoupled Helmholtz resonator, respectively. (4a)mpd2xpdt2=F-Rpdxpdt-kp+μA2xp(4b)mad2xadt2=-Radxadt-μS2xaFrom (4a), we can get the resonant frequency of the uncoupled top plate as(5)ωp=kp+μA2mpReferencing from [18], this frequency can be calculated more accurately as the fundamental resonant frequency of a clamped circular membrane as(6)ωp=2πfp=2π·0.47tR2Yρm1-σ2where ρm,σ, and Y are membrane material’s density, Poisson’s ratio and Young’s modulus, respectively. From (4b), the resonant frequency of an uncoupled Helmholtz resonator can be solved as (7)ωh=2πfh=μS2maWe substituteμ=c2ρ/V, ma=Slρ into (7), and we can get (8)ωh=cSVlAccording to [19], the length of the Helmholtz resonator’s neck l should take an effective valueleff=l+δ for better estimation; thus,(9)ωh=cSVl+δwhereδ is the end correction. Referencing from [19], the end correction δ should take the sum of the interior correction δi and the exterior correction δe as follows:(10)δ=δi+δe=0.48S1-1.25nholeaR+0.48Swherenhole is the number of holes in the membrane of the CMUT cell.Coupled Helmholtz resonator and top plate: now, we focus on the coupled vibration system, assuming in harmonic oscillation state, and we solve (3a) and (3b) in frequency domain, by writing in phasor format xp=Xpeiωt, xa=Xaeiωt, F=F0eiωt, the amplitude of top plate’s velocity Up=iωXp, and the amplitude of equivalent air piston’s velocity Ua=iωXa. Similar to [16], we can get (11a)Up=iωF0ωh2-ω2+iωγampD(11b)Ua=-iωF0Aωh2mpSDwhere(12)D=ωp2-ω2+iωγpωh2-ω2+iωγa-ωph4(13)γa=Rama,γp=Rpmp,ωph4=ωh2ωa2,ωa2=μA2mpωp and ωh are the resonant frequencies of uncoupled top plate and uncoupled Helmholtz resonator, respectively.From (11a) and (11b), the coupled vibration system shown in Figure 2 resonates while the moduli of Ua and Up arrive maximum (D=0). If we neglect the damping terms ofiωγa andiωγp in (12), then this coupled vibration system’s resonant frequencies can be solved as(14)ω±2=ωp2+ωh22±ωp2-ωh22+4ωph42As shown in (6) and (9), for the uncoupled top plate and the uncoupled Helmholtz resonator, each of them only has one resonant frequency. From (14) we can see that, unlike the uncoupled situation, both the top plate and the equivalent air piston have two resonant frequencies in the coupled vibration system.For an air-coupled unsealed Helmholtz structural CMUT cell, when calculatingωh by (9), the length of the neck l is equal to the membrane’s thickness t. And the lumped parameters are calculated as follows:(15)S=nholeπa2(16)A=πR2-S(17)V=πR2h(18)mp=1.8Atρm+mrp(19)ma=Stρ+mrawhere the factor 1.8 in (18) is from [20], mrp is the acoustic radiation mass of the plate, and mra in (19) is the acoustic radiation mass of the equivalent air piston. Referenced from [21, 22], the CMUT cell should be considered as a clamped plate source located on an infinite baffle, and the acoustic radiation mass of the plate is (20)mrp=-ρcAω20F22ω/cA/πω/cA/π9,F2x=-x4-91x2+504H1x-14xx2-18H0x+14x415π-168x2πwhere Hn(x) is the nth order Struve function. The equivalent air piston is considered as a baffled piston source; the radiation mass of the equivalent air piston is [14](21)mra=nholeρcπa2ω2H12ω/ca2ω/caIn the simplification in Figure 2, the driving force F is a harmonic force; hence, we take F as the harmonic electrostatic force, according to [1],(22)F0=Fac=ε0AVdcVach2where Vdc, Vac are DC voltage and AC voltage applied to the cell, and ε0 is the dielectric permittivity of air.The squeeze film effect includes squeeze film stiffening effect and squeeze film damping effect. The proportion of the stiffening effect and damping effect is determined by the squeeze number [20](23)σsqeeze=12μairωR2Pah2whereμair is the dynamic viscosity of air, Pa is the ambient pressure, and ω is the plate’s vibrating angular frequency. For the air-coupled unsealed CMUT with perforated membrane, the squeeze damping effect dominates the squeeze film effect because of the small squeeze number, and the squeeze stiffening effect can be neglected [20, 23]. From [20], the squeeze damping can be represented by the lump resistance as (24)Rsqeeze=3πμairR4κβ2h3nholewhere(25)κβ=4β2-β4-4ln⁡β-3(26)β=a1.05r1/2where r1 represents the distance between center of hole and center of membrane andnhole is the total number of holes.For the plate, although there are other damping mechanisms, the squeeze film damping is the most noticeable once it is presented in MEMS devices [23]. In addition, CMUT is an acoustic device; the loss of acoustic radiation cannot be ignored [10]. According to [21, 22], the resistance of the acoustic radiation of the clamped plate is(27)Rrad_plate=ρcA1-20F12ω/cA/πω/cA/π9,F1x=x4-91x2+504J1x+14xx2-18J0x-x516-x7768where Jn(x) is the Bessel function of the first kind of order n. Thus, the damping resistance of the top plate Rp is (28)Rp=Rsqueeze+Rrad_platewhereRsqueeze is calculated as in (24) to (26).Rrad_plate is calculated as in (27).For the equivalent air piston, not only the acoustic radiation damping but also the viscous damping needs to be considered. When the air pistons vibrate, due to the viscosity of the air, they are damped by the viscous force at the side walls of the holes. From [24], the resistance of viscous boundary layer per unit area can be expressed as (29)Rvisc=ρμairω2When calculating the viscous damping of the equivalent air piston, the area of viscous boundary layer is taken as the area of the side walls of the holes and a viscous end correction [19]. The resistance of acoustic radiationRrad_air piston is calculated as the baffled piston source [14].(30)Rrad_airpiston=nholeρcπa21-2J12ωa/c2ωa/cThus, the damping resistance of the equivalent air piston is (31)Ra=nhole2πat+4πa2Rvisc+Rrad_airpistonAfter gettingUa and Up, the complex amplitude of total volume velocity U can be written as(32)U=AUp+SUaFrom [21, 25], if we consider a CMUT cell as a clamped edge radiator source, the complex amplitude of sound pressure of the point at distance d in the far field above the center of the membrane can be calculated as (33)p=iωρU2πdthen the sound pressure level (SPL) at this point can be calculated.


## 2.2. Finite Element Method Model

In order to investigate the air-coupled unsealed Helmholtz structural CMUTs and validate the analytical model, we have developed the finite element method (FEM) model to simulate the device. Unsealed CMUTs would induce squeeze film effect; usually the squeeze film effect in MEMS devices is described by nonlinear compressible Reynolds equation, which is simplified from the continuity equation and the Navier-Stokes equation when the modified Reynolds number satisfies [13, 20](34)Remodified=ωρh2μair≪1This condition is not effective for CMUTs which have slightly higher gaps or higher operating frequencies. Thus, we need to start with the continuity equation and the Navier-Stokes equation to model the squeeze film effect in our unsealed Helmholtz structural CMUTs [13].The commercial software COMSOL Multiphysics version 5.2a is used to build our FEM model. The Electromechanical (emi) module and Thermoviscous Acoustics, Frequency Domain (ta) module are used in our FEM model of a CMUT cell. The governing equations ofta include the Navier-Stokes equation, the continuity equation, and the energy equation formulated using the Fourier heat law [26]. Thus, the squeeze film effect is included and the thermal loss is also taken into account.We use four symmetric distributed holes here to avoid the damage of the center part of membrane which makes main contribution to membrane’s vibration. Table1 lists the material properties in simulation. The 3D model of a quarter of a CMUT cell in octant spherical air domain is built, as shown in Figures 3(a) and 3(b). A “perfectly matched layer (PML)” is used to eliminate the acoustic reflection and to simulate an infinite spread medium. Inemi module, for the “domain selection”, it includes all domains. The membrane, the insulator layer, and the support are selected as “linear elastic dielectric”. The rim of the membrane, the support, and the insulator layer are applied with “fixed constraint” boundary condition. The substrate can be ignored because of the fixed insulator layer. For simplification in FEM, the top and bottom electrodes can be regarded as infinitely thin [27]; thus, we directly use the membrane’s top surface as the top electrode, and the bottom electrode is located underneath the insulator layer. A constant DC bias voltage and a harmonic AC voltage (amplitude as 1% of DC voltage) are applied on the top surface of membrane by the “terminal” with “harmonic perturbation” conditions, and the bottom surface of insulator layer is set as “ground”. The “symmetry” is applied on the vertical sections of the membrane, insulator layer, and support. Other boundary conditions keep the default setting in COMSOL. Inta module, for the “domain selection”, only domains filled with air (cavity, holes, spherical air domain, PML) are selected. To couple two modules, “velocity” and “isothermal” boundary conditions are applied inta on the top surface of membrane and inner wall of cavity and of hole, which are interfaces betweenta andemi [26]. The “velocity” inta is prescribed as inemi. This coupling setting is a two-way coupling, which means when the membrane vibrates to generate sound waves in medium, a reaction force would be added to the membrane [28]. Thus we do not need to add other coupling condition inemi module. The “symmetry” boundary condition is applied on the vertical sections of cavity and the vertical sections of octant spherical air and PML. Other boundary conditions keep the default setting in COMSOL. In the mesh setup, the mesh size should be sufficiently small to resolve the wavelength by at least 5-6 elements per wavelength [29], and the maximum element size of the acoustic medium is set as one-sixth of wavelength of operating frequency. Also the “boundary layers mesh” is needed at the inner walls of both cavity and holes; the thickness of “boundary layers mesh” is equal to the viscous boundary layer thickness under fundamental resonant frequency of CMUTs membrane [26]. The meshing details of vertical section of the cell are shown in Figure 3(c), where the “boundary layers meshes” are highlighted. The “prestressed analysis” is applied in calculating the average pressure over different frequencies. Firstly, a stationary study is used to calculate the membrane deflection under a certain DC bias voltage. Then an AC voltage is biased on the membrane using harmonic analysis over a frequency range. The pressure is calculated by averaging the pressure between the membrane and the air medium.Table 1
Material properties used in COMSOL (∗: the factor 0.6 is from [30]; others are the preset values/functions in COMSOL).
Parameter Value/Preset function in COMSOL Air: (used in cavity, holes, spherical air domain, PML) Relative permittivity 1 Dynamic viscosityµair [Pa.s] eta(T[1/K]) Heat capacity at constant pressureCp [J/(kg.K)] Cp(T[1/K]) Densityρair [kg/m3] rho(pA[1/Pa],T[1/K]) Thermal conductivityK [W/(m.K)] k(T[1/K]) ∗Bulk viscosity µb [Pa.s] 0.6∗eta(T[1/K]) Typical Wave Speed for PML [m/s] 343  Si: (used in membrane, support beam) Young’s modulusY [GPa] 170 Poisson ratioσ 0.28 DensityρSi [kg/m3] 2329 Relative permittivity 11.7  SiO 2 : (used in insulator layer) Young’s modulus [GPa] 70 Poisson ratio 0.17 Density [kg/m3] 2200 Relative permittivity 4.2Figure 3
A 3D model in COMSOL: (a) a quarter of a CMUT cell and spherical air domain; (b) details of the unsealed Helmholtz structural CMUT cell; (c) meshing details of the vertical section of (b).In order to make a comparison between the unsealed Helmholtz structural CMUTs and the conventional vacuum sealed CMUTs, we also build an FEM model for the latter. The model building process of the conventional cell is slightly different from the unsealed Helmholtz structural CMUT cell. The solid-air interfaces for the conventional cell only contain the top surface of the membrane; thus, the “velocity” and “isothermal” are only added on the membrane’s top surface. The cavity is vacuum; thus, for a conventional cell, we do not take the acoustic field of the cavity intota module’s calculation. Also, the “boundary layers mesh” is not needed for the vacuum cavity.


## 3. Results

### 3.1. FEM Results

We take three (Standard, Double, and Half) cases in FEM simulation; Table2 shows the parameters of physical structural dimension of the CMUT cell used. The main target of this paper is to demonstrate synergy properties of both the squeeze film effect and the Helmholtz resonant effect in CMUTs. In order to fulfill this, we set the corresponding parameters and they turn out to be a bit large for Standard case. In the Standard case, the radius R and thickness t of membrane are referenced from [12] as a representative one. The cavity height h is designed to achieve a properh/2σvisc ratio to generate the obvious Helmholtz effect, whereσvisc is the thickness of the viscous boundary layer, as we will explain in the discussion section. For air-coupled CMUTs, designs and fabrications with such large cell dimensions can be found in [31, 32]. The fabrication method of the vented holes can be found in [12, 13]. In the Double case, the physical dimension parameters are 2 times of the Standard case. In the frequency response, the vibration coupling characteristic of two resonant peaks is more obvious with such large dimension cells, because the proportion of thermal viscous loss energy to the total energy gets lower when the cell dimension increases from Half case to Double case, as we shall show in discussion section. The results of the Double case with obvious characteristic are helpful to clearly investigate the working mechanism of Helmholtz structural CMUT. In the Half case, the physical dimension parameters are taken as half of the Standard case. As we shall show in discussion section, this case is to investigate influence factor of the Helmholtz effect.Table 2
Parameters of CMUT cells in our design.
Parameters Double Standard Half R membrane radius [µm] 1000 500 250 h cavity height [µm] 60 30 15 t membrane thickness [µm] 20 10 5 a ini initial hole radius [µm] 110 55 27.5 a hole radius [µm] 40, 80, 120, 160, 200 20, 40, 60, 80, 100 10, 20, 30, 40, 50 t ins insulator layer thickness [µm] 6 3 1.5 r 1 hole position [µm] 750 375 187.5 n hole holes number 4 4 4 V dc DC voltage [V] 6720 3360 1680For hole radius, we use (6) to calculate fp; let fh=fp, then use (9), (10), and (15) to solve an initial hole radiusaini of each hole as(35)aini=0.96nholeπqa+0.962nholeπqa2-41+qaqb-qat21+qaqbwhere(36)qa=2πfpc2Vnholeπ(37)qb=0.6πnholeRThen choose hole radius aroundaini. In each case in Table 2, compared with the no-hole cell, the collapse voltage of the CMUTs with 4 holes decreases a little. For simplicity, the DC voltage is taken as about 70% of the collapse voltage of the no-hole cell.Firstly, the Double case is considered to show the vibration coupling working mechanism of the air-coupled unsealed Helmholtz structural CMUTs.The membrane and air piston’s motions are mainly on vertical direction; thus, in COMSOL, we take the amplitude of the membrane’s velocityUp and the amplitude of the air piston’s velocity Ua as the z directional average velocity complex amplitude of the membrane’s top surface and of the hole’s top surface, respectively.Figure4 shows magnitudes of Up and Ua. As shown in Figure 4, for unsealed Helmholtz structural CMUTs with different hole sizes, as mentioned in methodologies section, the magnitudes of Up and Ua have two resonant frequencies corresponding to f-=ω-/2π, f+=ω+/2π, respectively. It is noticeable that between these two resonant peaks there is a trough, whose corresponding frequency is called the antiresonant frequency. Referring to [16], this antiresonant frequency is equal to the resonant frequency of uncoupled Helmholtz resonator fh.Figure 4
FEM results of the Double case: (a) magnitude of complex amplitude of air piston velocity; (b) magnitude of complex amplitude of membrane velocity.
(a) (b)The FEM complex amplitude of total volume velocityU is calculated as (32). Figure 5 shows its magnitude U. Also, U of the unsealed Helmholtz structural CMUTs has two resonant frequencies, locating at two sides of the peak of the vacuum sealed cell at fundamental resonant frequency fp = 83 kHz.Figure 5
FEM results of the Double case: magnitude of complex amplitude of total volume velocity.Figure6 illustrates the air pistons’ and the membrane’s contributions to the magnitudes of total volume velocity U of the unsealed Helmholtz structural CMUT cells. Contributions of air pistons and the membrane are represented by the magnitudes SUa and AUp, respectively. In order to mark the fundamental resonant frequency fp, we also add the magnitudes of total volume velocity of the vacuum sealed cell. From Figures 6(a) and 6(b), for the unsealed Helmholtz structural CMUT cells, when fh < fp, the air pistons make major contribution to the left peak of U at f-, and to the right peak atf+ while fh > fp. In simpler terms, for the unsealed Helmholtz structural CMUT cells, the air pistons mainly contribute to the resonant peak whose corresponding frequency is farther from fp than the other resonant peak.Figure 6
FEM results of the Double case: (a)SUa, AUp, U ofa=40μm dimension (fh<fp); (b) SUa, AUp, U ofa=160μm dimension (fh>fp).
(a) (b)The CMUTs from the Double case to the Half case have different working frequencies; thus, at the point with same distance from the CMUT cells, the acoustic attenuation differs among cases. In order to fairly compare the output pressure from the Double case to the Half case, the SPL is taken at the top surfaces of membrane and hole. We first take the area average sound pressure of the top surfaces of membrane and hole then calculate the SPL. Figure7 shows the SPLs over the membrane and hole’s top surface. From (33), the sound pressure is proportional toU; on this account, each curve in Figure 7 has two resonant peaks at resonant frequencies of its corresponding curve in Figure 5. As shown in Figure 7, compared with the sealed one, the SPLs of unsealed Helmholtz structural CMUTs realize the dual-band and extend range of frequency response towards wider frequencies.Figure 7
FEM results of the Double case: sound pressure level (SPL) over the membrane and hole’s top surface (dB).Secondly, for the Standard case in Table2, the FEM results of SPLs are shown in Figure 8. Two resonant frequencies due to the vibration coupling also can be observed for the unsealed cells. As we can see from Figure 8, for SPLs, the bandwidths of unsealed Helmholtz structural CMUT cells are obviously larger than the conventional sealed CMUT cell. However, compared with the sealed cell, the energy loss in the form of squeeze damping loss, viscous loss, and thermal loss is introduced in the cavity and holes of the unsealed Helmholtz structural CMUT cell. Under the same Vdc and Vac, for the sealed cell, the electrical energy is converted into the acoustic energy, while for the unsealed cell, the electrical energy is converted into the acoustic energy and the energy loss in the cavity and holes. Thus the pressure values of the Helmholtz structural CMUT cells are lower than the sealed one. Referring to [33], we define the SPL-bandwidth product as bandwidth [MHz]×peak value of SPL [dB]: both bandwidth and pressure values performance of the CMUTs can be reflected by this product. We calculate the 3-dB fractional bandwidth and SPL-bandwidth product of the CMUT cells in Figure 8; the results are listed in Table 3.Table 3
3-dB fractional bandwidth and SPL-bandwidth product of the Standard case as shown in Figure8.
Hole radius [µm] a=20 a=40 a=60 a=80 a=100 no hole f h [kHz] 92.3 138.5 176.5 211.3 244.7 f p=167.5 Separation betweenfp and fh [kHz] 75.2 29 9 43.8 77.2 3-dB fractional bandwidth 4.8% 11.6% 35.7% 7.3% 3.8% 0.9% SPL-bandwidth product 1.22 2.88 8.01 1.51 0.77 0.23Figure 8
FEM results of the Standard case: sound pressure level (SPL) over the membrane and holes (dB).As shown in Table3, the fp is calculated from (6) and the fh is calculated from (9) and (10); when fh is close to fp, the unsealed Helmholtz structural CMUT cell witha= 60 μm achieves a significant improvement around 35 times of the conventional sealed cell in both 3-dB fractional bandwidth and SPL-bandwidth product. Other unsealed dimensions also have performance enhancement compared with the conventional sealed one. Therefore, compared with the conventional sealed cell, unsealed Helmholtz structural CMUT cells not only enhance the bandwidth of SPL but also keep the pressure value at a reasonable level.Last, for the Half case in Table2, the FEM results of the SPLs over the top surfaces of the membrane and hole are shown in Figure 9, and the corresponding 3-dB fractional bandwidth and the SPL-bandwidth product are listed in Table 4. The characteristic of two resonant peaks is not so obvious as the Double case and the Standard case. However, the significant enhancement around 17 times of the conventional sealed cell still is achieved while fh is close to fp (a= 30 μm cell). Unsealed Helmholtz structural CMUT cells still have a performance enhancement over the conventional sealed CMUT cell.Table 4
3-dB fractional bandwidth and SPL-bandwidth product of the Half case as shown in Figure9.
Hole radius [µm] a=10 a=20 a=30 a=40 a=50 no hole f h [kHz] 184.6 277.1 353.1 422.6 489.6 f p=334.9 Separation betweenfp and fh [kHz] 150.3 57.8 18.2 87.7 154.7 3-dB fractional bandwidth 5.4% 10.9% 24.0% 14.7% 6.5% 1.3% SPL-bandwidth product 2.67 5.29 11.06 6.08 2.62 0.67Figure 9
FEM results of the Half case: sound pressure level (SPL) over the membrane and holes (dB).As shown above in Tables3 and 4, compared with the conventional sealed CMUT cell, the air-coupled unsealed 4-hole Helmholtz structural CMUT cells achieved the enhancement of both 3-dB fractional bandwidth and SPL-bandwidth product, and the significant enhancement was achieved while fh is close to fp.


### 3.2. Analytical Results Compared with FEM Results

In order to validate our analytical model, we make a comparison between analytical results and FEM results from all cases. For the analytical results, we first use (15) to (22), (28) and (31) to calculate values of basic parameters; then, we use (6) to get ωp, and we use (9) and (10) to get ωh. After that we use (11a) to (13) to get Up and Ua, then we substitute back (32) to get U and use (33) to get p at a point above (d= 3600 μm in the Double case,d=1800 μm in the Standard case,d= 900 μm in the Half case) the center of the membrane. Then we calculate SPL at this point. In this section, the FEM SPL results also are calculated at the same corresponding point above the center of the membrane. The comparison results are shown in Figure 10.Figure 10
Comparison of analytical and FEM results of SPL at a point above the center of the membrane.One point we would like to mention, in (20) and (27), we calculate the acoustic radiation impedance of the CMUT membrane by the formulas of the baffled clamped plate source. From [22], the piston plate radiation and the clamped plate radiation are almost same forkRsource<1, where k is the wave number andRsource is the source radius. In our design, the value ofkRsource is small; thus, we can also use baffled piston source to calculate the acoustic radiation impedance of the CMUT membrane; the analytical results are almost same by using these two ways to calculate acoustic radiation impedance of the CMUT membrane.In order to evaluate the similarity between the analytical SPL results and the FEM SPL results, we calculate the normalized mean square error (NMSE) between our analytical results and FEM results [34].(38)NMSEvalue=1-x→ref-x→2x→ref-meanx→ref2The NMSE values vary between negative infinity (not fit) and 1 (perfect fit). In each subgraph in Figure10, we take FEM results and analytical results as x→ref and x→, respectively. The calculated NMSE values in Table 5 indicate our analytical results have a good matching with FEM results.Table 5
Normalized mean square error (NMSE) values between analytical results and FEM results in Figure10.
Double case a=80µm a=120µm NMSE 0.9173 0.9347  Standard case a=40µm a=60µm NMSE 0.9355 0.9513  Half case a =20µm a=30µm NMSE 0.9523 0.9689The deviations between analytical results and FEM results could be caused by the simplifications of the analytical model such as negligence of the membrane’s deformation and assumption as a plate with uniform displacement. In (22), the analytical driving force is only taken as the harmonic part of the electrical force. In (28) and (31), the loss taken into consideration is only the squeezing, viscous and radiation loss, not including the thermal loss which is included in the FEM model. And in the analytical model we use the equivalent air piston instead of four small air pistons in the holes; thus, the interaction among holes is neglected. Compared with the no-hole cell, the membrane stiffness and the uncoupled natural resonant frequency of the 4-hole CMUT cell decrease. In [15], the natural resonant frequency of the membrane with holes was approximated by the formula of the no-hole membrane. For our design, using (6) of the no-hole membrane to approximate the natural resonant frequency of the 4-hole membrane would not have a significant influence of the analytical results. In the future work, the accurate formula of the natural resonant frequency of the membrane with holes will be investigated.However, this analytical model gives the explanation of the most significant vibration coupling mechanism of the air-coupled unsealed Helmholtz structural CMUTs and can be used to predict the characteristics such asU, SPL andf+ and f- of this kind of CMUTs during the design process. Figure 11 shows the flow chart on how the analytical model can be used to aid the design of Helmholtz structural CMUTs. In the design phase of Helmholtz structural CMUTs, for the membrane with inputted thickness t and radiusR, we firstly select an h value within the interval [4σvisc, hmax], where the lower limit 4σvisc is for the obvious Helmholtz resonance, and the consideration of the limitation of DC voltage and maximum allowed gap height is included in hmax. Then sweep hole radius a around the calculated initial hole radiusaini to get the point SPL and perform the FEM simulation of the dimension neighbored around the hole radius with significant performance enhancement to search for more accurate fine results. If there is no significant performance enhancement for all of the swept hole radii, then we need to reselect the h value, and sweep hole radius a around the new calculatedaini to find performance enhanced cell dimension. If there is no properh, we need to go back to reinput t and R. Using this design process, the analytical model can help us to initially pick out the cell dimensions with performance enhancement. Then we only need to perform the FEM simulations of these dimensions to check more accurate results. This simplifies the design process and shortens the design period, saving a lot of computation time with our reasonable guesstimates of analytical design procedure.Figure 11
The flow chart of using analytical model to aid the design process of Helmholtz structural CMUTs.


## 3.1. FEM Results

We take three (Standard, Double, and Half) cases in FEM simulation; Table2 shows the parameters of physical structural dimension of the CMUT cell used. The main target of this paper is to demonstrate synergy properties of both the squeeze film effect and the Helmholtz resonant effect in CMUTs. In order to fulfill this, we set the corresponding parameters and they turn out to be a bit large for Standard case. In the Standard case, the radius R and thickness t of membrane are referenced from [12] as a representative one. The cavity height h is designed to achieve a properh/2σvisc ratio to generate the obvious Helmholtz effect, whereσvisc is the thickness of the viscous boundary layer, as we will explain in the discussion section. For air-coupled CMUTs, designs and fabrications with such large cell dimensions can be found in [31, 32]. The fabrication method of the vented holes can be found in [12, 13]. In the Double case, the physical dimension parameters are 2 times of the Standard case. In the frequency response, the vibration coupling characteristic of two resonant peaks is more obvious with such large dimension cells, because the proportion of thermal viscous loss energy to the total energy gets lower when the cell dimension increases from Half case to Double case, as we shall show in discussion section. The results of the Double case with obvious characteristic are helpful to clearly investigate the working mechanism of Helmholtz structural CMUT. In the Half case, the physical dimension parameters are taken as half of the Standard case. As we shall show in discussion section, this case is to investigate influence factor of the Helmholtz effect.Table 2
Parameters of CMUT cells in our design.
Parameters Double Standard Half R membrane radius [µm] 1000 500 250 h cavity height [µm] 60 30 15 t membrane thickness [µm] 20 10 5 a ini initial hole radius [µm] 110 55 27.5 a hole radius [µm] 40, 80, 120, 160, 200 20, 40, 60, 80, 100 10, 20, 30, 40, 50 t ins insulator layer thickness [µm] 6 3 1.5 r 1 hole position [µm] 750 375 187.5 n hole holes number 4 4 4 V dc DC voltage [V] 6720 3360 1680For hole radius, we use (6) to calculate fp; let fh=fp, then use (9), (10), and (15) to solve an initial hole radiusaini of each hole as(35)aini=0.96nholeπqa+0.962nholeπqa2-41+qaqb-qat21+qaqbwhere(36)qa=2πfpc2Vnholeπ(37)qb=0.6πnholeRThen choose hole radius aroundaini. In each case in Table 2, compared with the no-hole cell, the collapse voltage of the CMUTs with 4 holes decreases a little. For simplicity, the DC voltage is taken as about 70% of the collapse voltage of the no-hole cell.Firstly, the Double case is considered to show the vibration coupling working mechanism of the air-coupled unsealed Helmholtz structural CMUTs.The membrane and air piston’s motions are mainly on vertical direction; thus, in COMSOL, we take the amplitude of the membrane’s velocityUp and the amplitude of the air piston’s velocity Ua as the z directional average velocity complex amplitude of the membrane’s top surface and of the hole’s top surface, respectively.Figure4 shows magnitudes of Up and Ua. As shown in Figure 4, for unsealed Helmholtz structural CMUTs with different hole sizes, as mentioned in methodologies section, the magnitudes of Up and Ua have two resonant frequencies corresponding to f-=ω-/2π, f+=ω+/2π, respectively. It is noticeable that between these two resonant peaks there is a trough, whose corresponding frequency is called the antiresonant frequency. Referring to [16], this antiresonant frequency is equal to the resonant frequency of uncoupled Helmholtz resonator fh.Figure 4
FEM results of the Double case: (a) magnitude of complex amplitude of air piston velocity; (b) magnitude of complex amplitude of membrane velocity.
(a) (b)The FEM complex amplitude of total volume velocityU is calculated as (32). Figure 5 shows its magnitude U. Also, U of the unsealed Helmholtz structural CMUTs has two resonant frequencies, locating at two sides of the peak of the vacuum sealed cell at fundamental resonant frequency fp = 83 kHz.Figure 5
FEM results of the Double case: magnitude of complex amplitude of total volume velocity.Figure6 illustrates the air pistons’ and the membrane’s contributions to the magnitudes of total volume velocity U of the unsealed Helmholtz structural CMUT cells. Contributions of air pistons and the membrane are represented by the magnitudes SUa and AUp, respectively. In order to mark the fundamental resonant frequency fp, we also add the magnitudes of total volume velocity of the vacuum sealed cell. From Figures 6(a) and 6(b), for the unsealed Helmholtz structural CMUT cells, when fh < fp, the air pistons make major contribution to the left peak of U at f-, and to the right peak atf+ while fh > fp. In simpler terms, for the unsealed Helmholtz structural CMUT cells, the air pistons mainly contribute to the resonant peak whose corresponding frequency is farther from fp than the other resonant peak.Figure 6
FEM results of the Double case: (a)SUa, AUp, U ofa=40μm dimension (fh<fp); (b) SUa, AUp, U ofa=160μm dimension (fh>fp).
(a) (b)The CMUTs from the Double case to the Half case have different working frequencies; thus, at the point with same distance from the CMUT cells, the acoustic attenuation differs among cases. In order to fairly compare the output pressure from the Double case to the Half case, the SPL is taken at the top surfaces of membrane and hole. We first take the area average sound pressure of the top surfaces of membrane and hole then calculate the SPL. Figure7 shows the SPLs over the membrane and hole’s top surface. From (33), the sound pressure is proportional toU; on this account, each curve in Figure 7 has two resonant peaks at resonant frequencies of its corresponding curve in Figure 5. As shown in Figure 7, compared with the sealed one, the SPLs of unsealed Helmholtz structural CMUTs realize the dual-band and extend range of frequency response towards wider frequencies.Figure 7
FEM results of the Double case: sound pressure level (SPL) over the membrane and hole’s top surface (dB).Secondly, for the Standard case in Table2, the FEM results of SPLs are shown in Figure 8. Two resonant frequencies due to the vibration coupling also can be observed for the unsealed cells. As we can see from Figure 8, for SPLs, the bandwidths of unsealed Helmholtz structural CMUT cells are obviously larger than the conventional sealed CMUT cell. However, compared with the sealed cell, the energy loss in the form of squeeze damping loss, viscous loss, and thermal loss is introduced in the cavity and holes of the unsealed Helmholtz structural CMUT cell. Under the same Vdc and Vac, for the sealed cell, the electrical energy is converted into the acoustic energy, while for the unsealed cell, the electrical energy is converted into the acoustic energy and the energy loss in the cavity and holes. Thus the pressure values of the Helmholtz structural CMUT cells are lower than the sealed one. Referring to [33], we define the SPL-bandwidth product as bandwidth [MHz]×peak value of SPL [dB]: both bandwidth and pressure values performance of the CMUTs can be reflected by this product. We calculate the 3-dB fractional bandwidth and SPL-bandwidth product of the CMUT cells in Figure 8; the results are listed in Table 3.Table 3
3-dB fractional bandwidth and SPL-bandwidth product of the Standard case as shown in Figure8.
Hole radius [µm] a=20 a=40 a=60 a=80 a=100 no hole f h [kHz] 92.3 138.5 176.5 211.3 244.7 f p=167.5 Separation betweenfp and fh [kHz] 75.2 29 9 43.8 77.2 3-dB fractional bandwidth 4.8% 11.6% 35.7% 7.3% 3.8% 0.9% SPL-bandwidth product 1.22 2.88 8.01 1.51 0.77 0.23Figure 8
FEM results of the Standard case: sound pressure level (SPL) over the membrane and holes (dB).As shown in Table3, the fp is calculated from (6) and the fh is calculated from (9) and (10); when fh is close to fp, the unsealed Helmholtz structural CMUT cell witha= 60 μm achieves a significant improvement around 35 times of the conventional sealed cell in both 3-dB fractional bandwidth and SPL-bandwidth product. Other unsealed dimensions also have performance enhancement compared with the conventional sealed one. Therefore, compared with the conventional sealed cell, unsealed Helmholtz structural CMUT cells not only enhance the bandwidth of SPL but also keep the pressure value at a reasonable level.Last, for the Half case in Table2, the FEM results of the SPLs over the top surfaces of the membrane and hole are shown in Figure 9, and the corresponding 3-dB fractional bandwidth and the SPL-bandwidth product are listed in Table 4. The characteristic of two resonant peaks is not so obvious as the Double case and the Standard case. However, the significant enhancement around 17 times of the conventional sealed cell still is achieved while fh is close to fp (a= 30 μm cell). Unsealed Helmholtz structural CMUT cells still have a performance enhancement over the conventional sealed CMUT cell.Table 4
3-dB fractional bandwidth and SPL-bandwidth product of the Half case as shown in Figure9.
Hole radius [µm] a=10 a=20 a=30 a=40 a=50 no hole f h [kHz] 184.6 277.1 353.1 422.6 489.6 f p=334.9 Separation betweenfp and fh [kHz] 150.3 57.8 18.2 87.7 154.7 3-dB fractional bandwidth 5.4% 10.9% 24.0% 14.7% 6.5% 1.3% SPL-bandwidth product 2.67 5.29 11.06 6.08 2.62 0.67Figure 9
FEM results of the Half case: sound pressure level (SPL) over the membrane and holes (dB).As shown above in Tables3 and 4, compared with the conventional sealed CMUT cell, the air-coupled unsealed 4-hole Helmholtz structural CMUT cells achieved the enhancement of both 3-dB fractional bandwidth and SPL-bandwidth product, and the significant enhancement was achieved while fh is close to fp.


## 3.2. Analytical Results Compared with FEM Results

In order to validate our analytical model, we make a comparison between analytical results and FEM results from all cases. For the analytical results, we first use (15) to (22), (28) and (31) to calculate values of basic parameters; then, we use (6) to get ωp, and we use (9) and (10) to get ωh. After that we use (11a) to (13) to get Up and Ua, then we substitute back (32) to get U and use (33) to get p at a point above (d= 3600 μm in the Double case,d=1800 μm in the Standard case,d= 900 μm in the Half case) the center of the membrane. Then we calculate SPL at this point. In this section, the FEM SPL results also are calculated at the same corresponding point above the center of the membrane. The comparison results are shown in Figure 10.Figure 10
Comparison of analytical and FEM results of SPL at a point above the center of the membrane.One point we would like to mention, in (20) and (27), we calculate the acoustic radiation impedance of the CMUT membrane by the formulas of the baffled clamped plate source. From [22], the piston plate radiation and the clamped plate radiation are almost same forkRsource<1, where k is the wave number andRsource is the source radius. In our design, the value ofkRsource is small; thus, we can also use baffled piston source to calculate the acoustic radiation impedance of the CMUT membrane; the analytical results are almost same by using these two ways to calculate acoustic radiation impedance of the CMUT membrane.In order to evaluate the similarity between the analytical SPL results and the FEM SPL results, we calculate the normalized mean square error (NMSE) between our analytical results and FEM results [34].(38)NMSEvalue=1-x→ref-x→2x→ref-meanx→ref2The NMSE values vary between negative infinity (not fit) and 1 (perfect fit). In each subgraph in Figure10, we take FEM results and analytical results as x→ref and x→, respectively. The calculated NMSE values in Table 5 indicate our analytical results have a good matching with FEM results.Table 5
Normalized mean square error (NMSE) values between analytical results and FEM results in Figure10.
Double case a=80µm a=120µm NMSE 0.9173 0.9347  Standard case a=40µm a=60µm NMSE 0.9355 0.9513  Half case a =20µm a=30µm NMSE 0.9523 0.9689The deviations between analytical results and FEM results could be caused by the simplifications of the analytical model such as negligence of the membrane’s deformation and assumption as a plate with uniform displacement. In (22), the analytical driving force is only taken as the harmonic part of the electrical force. In (28) and (31), the loss taken into consideration is only the squeezing, viscous and radiation loss, not including the thermal loss which is included in the FEM model. And in the analytical model we use the equivalent air piston instead of four small air pistons in the holes; thus, the interaction among holes is neglected. Compared with the no-hole cell, the membrane stiffness and the uncoupled natural resonant frequency of the 4-hole CMUT cell decrease. In [15], the natural resonant frequency of the membrane with holes was approximated by the formula of the no-hole membrane. For our design, using (6) of the no-hole membrane to approximate the natural resonant frequency of the 4-hole membrane would not have a significant influence of the analytical results. In the future work, the accurate formula of the natural resonant frequency of the membrane with holes will be investigated.However, this analytical model gives the explanation of the most significant vibration coupling mechanism of the air-coupled unsealed Helmholtz structural CMUTs and can be used to predict the characteristics such asU, SPL andf+ and f- of this kind of CMUTs during the design process. Figure 11 shows the flow chart on how the analytical model can be used to aid the design of Helmholtz structural CMUTs. In the design phase of Helmholtz structural CMUTs, for the membrane with inputted thickness t and radiusR, we firstly select an h value within the interval [4σvisc, hmax], where the lower limit 4σvisc is for the obvious Helmholtz resonance, and the consideration of the limitation of DC voltage and maximum allowed gap height is included in hmax. Then sweep hole radius a around the calculated initial hole radiusaini to get the point SPL and perform the FEM simulation of the dimension neighbored around the hole radius with significant performance enhancement to search for more accurate fine results. If there is no significant performance enhancement for all of the swept hole radii, then we need to reselect the h value, and sweep hole radius a around the new calculatedaini to find performance enhanced cell dimension. If there is no properh, we need to go back to reinput t and R. Using this design process, the analytical model can help us to initially pick out the cell dimensions with performance enhancement. Then we only need to perform the FEM simulations of these dimensions to check more accurate results. This simplifies the design process and shortens the design period, saving a lot of computation time with our reasonable guesstimates of analytical design procedure.Figure 11
The flow chart of using analytical model to aid the design process of Helmholtz structural CMUTs.


## 4. Discussion

### 4.1. Explanation of the Decreased Phenomenon of the Two Resonant Peaks

In this section, we make a comparison of the FEM results of SPLs from the Double case to the Half case and give an explanation of the decreased phenomenon of the two resonant peaks. As we can see from Figure8, for the SPL of the unsealed Helmholtz structural CMUT cells in the Standard case, the left peaks of cells with fh < fp and the right peaks of cells with fh > fp are not as evident as the Double case in Figure 7. For the Half case in Figure 9, for the unsealed Helmholtz structural CMUT cells, the peaks majorly contributed by air pistons appear further weakly and even almost disappear. This is due to the thermal viscous loss, which is caused by viscosity and heat conduction.When sound propagates in structures with small geometric dimensions, the interaction between air and solid boundaries gives rise to the viscous boundary layer and the thermal boundary layer. In these two boundary layers, the loss is much greater than in the free field which is far away from the solid boundaries. The thicknesses of these two boundary layers are calculated as [24](39)σvisc=2μairρairω≈2.2×10-3f(40)σtherm=2KρairCpω≈2.5×10-3fFor viscous boundary layer, as Figure12 shows, we can regard the flow in the cavity and the holes as shear flow; it can be divided into two parts: the part withinσvisc is heavily damped and with great mechanical loss and the other part far fromσvisc is free and lossless. As listed in Table 6, the ratio h/2σvisc reduces from the Double case to the Half case. Therefore, for the whole flow in the cavity and the holes, the proportion of the part with heavy loss increases and the proportion of the lossless part decreases. That means more portion of the energy is dissipated by the boundary layers, and portion of the energy provided for the air pistons’ vibration gets weaker; as a consequence, the peak mostly contributed by air pistons gets weaker. For thermal boundary layer, also the proportion of the total energy to the thermal loss energy can be reflected by the ratioh/2σtherm. For our designed CMUTs parameters, the thermal loss can be neglected compared with the squeeze film loss and the radiation loss. Thus in our analytical model, we did not take the thermal loss into account.Table 6
Ratio of cavity height to 2σvisc, taking f as fp in (6).
Double case Standard case Half case h/2σvisc 3.94 2.79 1.97Figure 12
Illustration of viscous boundary layer and flow in cavity and holes.We have also simulated other two cases whose dimensions are 10 times and 0.1 times of the Standard case in Table2. The results are in accordance with our explanation above. For the 10 times case, the two resonant peaks are obvious. But for the 0.1 times case, the Helmholtz resonant effect tends to disappear; only one resonant peak is observed in the SPL frequency response. We also did FEM simulations and analytical analyses for aR = 650 μm,t = 10 μm, andh = 30 μm 4-hole cell with different hole radii. Compared with the cell in Standard case, the Helmholtz effect gets weaker because of the decreasedh/2σvisc ratio. Thus when we sweep hole radii from 25 μm to 45 μm aroundaini = 36 μm, the two resonant peaks cannot be within 3 dB level and cannot achieve the 3 dB bandwidth improvement likea = 60 μm cell in Figure 8.


### 4.2. Contribution of Helmholtz Resonant Effect to the Output Enhancement of the 4-Hole Cell

Here we check the bandwidth and SPL-bandwidth product enhancement of the unsealed Helmholtz structural CMUTs is indeed affected by the Helmholtz resonator, although the vibration coupling feature of two resonant peaks is not obvious in Figures8 and 9. The hole area ratio is defined as the total hole area divided by the total membrane area, i.e.,nholeπa2/πR2. We simulate an unsealed cell with 68 holes uniformly distributed on the membrane and having the same hole area ratio as the cell (a=60μm, 4 holes) with significant enhancement in the Standard case. Other dimension parameters of this 68-hole cell remain the same as the Standard case in Table 2. Thus from (9), (10), and (15), for the 4-hole cell and the 68-hole cell, they have the same value of the uncoupled Helmholtz resonant frequency fh. The simulation results are shown in Figure 13. We can see from Figure 13, as we explained in Figure 12, in spite of the same fh, there is more air-solid contact area contributed by holes’ side wall in the 68-hole cell. Hence, there are more lossy flows which weaken the vibration of air pistons of Helmholtz resonator. Therefore, the performance of 68-hole CMUT with decreased Helmholtz resonator effect is much less than the 4-hole one’s.Figure 13
FEM results of the Standard case: SPL comparison of two cells with same hole area ratio,a=60μm 4 holes &a=14.5μm 68 holes (FBW: 3-dB fractional bandwidth, SPL-BWP: SPL-bandwidth product).We also simulate another 68-hole cell which has the same hole area ratio as the cell (a=30μm, 4 holes) with a significant enhancement in the Half case; other dimension parameters of this 68-hole cell are the same as the Half case in Table 2. As shown in Figure 14, the 4-hole cell still performs better than the 68-hole one, which means Helmholtz resonant effect still has a great influence for the performance improvement of the 4-hole unsealed Helmholtz structural CMUT cells in the Half case, although the characteristic of two resonant peaks is not obvious.Figure 14
FEM results of the Half case: SPL comparison of two cells with same hole area ratio,a=30μm 4 holes &a=7.25μm 68 holes (FBW: 3-dB fractional bandwidth, SPL-BWP: SPL-bandwidth product).


### 4.3. Range of Effectiveness for Analytical Model

From discussion above, when theh/2σvisc ratio decreases, or the hole number increases under same hole area ratio, the Helmholtz resonant effect gets weaker and even disappears. For our analytical model, the existence of Helmholtz resonant effect is a basic assumption. In Figure 11, when choosing the cavity heighth, the conditionh≥4σvisc is an FEM empirical value that the Helmholtz resonant effect is strong enough to benefit the bandwidth of the CMUT device. Our analytical model is effective under these situations where the Helmholtz resonant effect is strong. For smallerh/2σvisc ratio, the Helmholtz resonant effect tends to disappear, and our analytical model loses efficacy. Also, for 4-hole cells in Figures 13 and 14, the Helmholtz resonant effect is strong and our analytical model is effective. For 68-holes cells, the analytical model loses efficacy. In the future work, the research attempt is to investigate the effective range of our analytical model.


### 4.4. The “Best” Hole Radius for Bandwidth Enhancement

From Tables3 and 4, we can see, when fh is close to fp, the unsealed Helmholtz structural CMUT achieves a significant output pressure bandwidth performance enhancement. We can see a tendency that when fh=fp, the improvement should achieve the “best”. We use 5 μm step size to simulate Standard case in Figure 8 with hole radius changes from 40 μm to 80 μm; the results of SPL over the membrane and holes are shown below. As Figure 15 shows, thea=55 μm,a=60 μm, anda=65 μm cells have larger 3-dB bandwidth than other dimensions. The 3-dB fractional bandwidth are 29.4%, 35.7%, and 29.7% for a=55 μm, 60 μm, and 65 μm, respectively. The calculatedaini which let fh=fp is 55 μm; however, its fractional bandwidth enhancement is not the “best” one. For the “best” fractional bandwidth enhancement witha=60 um, as we described above, the fh is close to fp. Up to now, we cannot give a quantitative analysis for how close are the fh and fp for the “best” bandwidth performance enhancement. This still needs to be further investigated.Figure 15
FEM results of the Standard case: sound pressure level (SPL) over the membrane and holes (dB), hole radius changes from 40μm to 80 μm with 5 μm step size.


## 4.1. Explanation of the Decreased Phenomenon of the Two Resonant Peaks

In this section, we make a comparison of the FEM results of SPLs from the Double case to the Half case and give an explanation of the decreased phenomenon of the two resonant peaks. As we can see from Figure8, for the SPL of the unsealed Helmholtz structural CMUT cells in the Standard case, the left peaks of cells with fh < fp and the right peaks of cells with fh > fp are not as evident as the Double case in Figure 7. For the Half case in Figure 9, for the unsealed Helmholtz structural CMUT cells, the peaks majorly contributed by air pistons appear further weakly and even almost disappear. This is due to the thermal viscous loss, which is caused by viscosity and heat conduction.When sound propagates in structures with small geometric dimensions, the interaction between air and solid boundaries gives rise to the viscous boundary layer and the thermal boundary layer. In these two boundary layers, the loss is much greater than in the free field which is far away from the solid boundaries. The thicknesses of these two boundary layers are calculated as [24](39)σvisc=2μairρairω≈2.2×10-3f(40)σtherm=2KρairCpω≈2.5×10-3fFor viscous boundary layer, as Figure12 shows, we can regard the flow in the cavity and the holes as shear flow; it can be divided into two parts: the part withinσvisc is heavily damped and with great mechanical loss and the other part far fromσvisc is free and lossless. As listed in Table 6, the ratio h/2σvisc reduces from the Double case to the Half case. Therefore, for the whole flow in the cavity and the holes, the proportion of the part with heavy loss increases and the proportion of the lossless part decreases. That means more portion of the energy is dissipated by the boundary layers, and portion of the energy provided for the air pistons’ vibration gets weaker; as a consequence, the peak mostly contributed by air pistons gets weaker. For thermal boundary layer, also the proportion of the total energy to the thermal loss energy can be reflected by the ratioh/2σtherm. For our designed CMUTs parameters, the thermal loss can be neglected compared with the squeeze film loss and the radiation loss. Thus in our analytical model, we did not take the thermal loss into account.Table 6
Ratio of cavity height to 2σvisc, taking f as fp in (6).
Double case Standard case Half case h/2σvisc 3.94 2.79 1.97Figure 12
Illustration of viscous boundary layer and flow in cavity and holes.We have also simulated other two cases whose dimensions are 10 times and 0.1 times of the Standard case in Table2. The results are in accordance with our explanation above. For the 10 times case, the two resonant peaks are obvious. But for the 0.1 times case, the Helmholtz resonant effect tends to disappear; only one resonant peak is observed in the SPL frequency response. We also did FEM simulations and analytical analyses for aR = 650 μm,t = 10 μm, andh = 30 μm 4-hole cell with different hole radii. Compared with the cell in Standard case, the Helmholtz effect gets weaker because of the decreasedh/2σvisc ratio. Thus when we sweep hole radii from 25 μm to 45 μm aroundaini = 36 μm, the two resonant peaks cannot be within 3 dB level and cannot achieve the 3 dB bandwidth improvement likea = 60 μm cell in Figure 8.


## 4.2. Contribution of Helmholtz Resonant Effect to the Output Enhancement of the 4-Hole Cell

Here we check the bandwidth and SPL-bandwidth product enhancement of the unsealed Helmholtz structural CMUTs is indeed affected by the Helmholtz resonator, although the vibration coupling feature of two resonant peaks is not obvious in Figures8 and 9. The hole area ratio is defined as the total hole area divided by the total membrane area, i.e.,nholeπa2/πR2. We simulate an unsealed cell with 68 holes uniformly distributed on the membrane and having the same hole area ratio as the cell (a=60μm, 4 holes) with significant enhancement in the Standard case. Other dimension parameters of this 68-hole cell remain the same as the Standard case in Table 2. Thus from (9), (10), and (15), for the 4-hole cell and the 68-hole cell, they have the same value of the uncoupled Helmholtz resonant frequency fh. The simulation results are shown in Figure 13. We can see from Figure 13, as we explained in Figure 12, in spite of the same fh, there is more air-solid contact area contributed by holes’ side wall in the 68-hole cell. Hence, there are more lossy flows which weaken the vibration of air pistons of Helmholtz resonator. Therefore, the performance of 68-hole CMUT with decreased Helmholtz resonator effect is much less than the 4-hole one’s.Figure 13
FEM results of the Standard case: SPL comparison of two cells with same hole area ratio,a=60μm 4 holes &a=14.5μm 68 holes (FBW: 3-dB fractional bandwidth, SPL-BWP: SPL-bandwidth product).We also simulate another 68-hole cell which has the same hole area ratio as the cell (a=30μm, 4 holes) with a significant enhancement in the Half case; other dimension parameters of this 68-hole cell are the same as the Half case in Table 2. As shown in Figure 14, the 4-hole cell still performs better than the 68-hole one, which means Helmholtz resonant effect still has a great influence for the performance improvement of the 4-hole unsealed Helmholtz structural CMUT cells in the Half case, although the characteristic of two resonant peaks is not obvious.Figure 14
FEM results of the Half case: SPL comparison of two cells with same hole area ratio,a=30μm 4 holes &a=7.25μm 68 holes (FBW: 3-dB fractional bandwidth, SPL-BWP: SPL-bandwidth product).


## 4.3. Range of Effectiveness for Analytical Model

From discussion above, when theh/2σvisc ratio decreases, or the hole number increases under same hole area ratio, the Helmholtz resonant effect gets weaker and even disappears. For our analytical model, the existence of Helmholtz resonant effect is a basic assumption. In Figure 11, when choosing the cavity heighth, the conditionh≥4σvisc is an FEM empirical value that the Helmholtz resonant effect is strong enough to benefit the bandwidth of the CMUT device. Our analytical model is effective under these situations where the Helmholtz resonant effect is strong. For smallerh/2σvisc ratio, the Helmholtz resonant effect tends to disappear, and our analytical model loses efficacy. Also, for 4-hole cells in Figures 13 and 14, the Helmholtz resonant effect is strong and our analytical model is effective. For 68-holes cells, the analytical model loses efficacy. In the future work, the research attempt is to investigate the effective range of our analytical model.


## 4.4. The “Best” Hole Radius for Bandwidth Enhancement

From Tables3 and 4, we can see, when fh is close to fp, the unsealed Helmholtz structural CMUT achieves a significant output pressure bandwidth performance enhancement. We can see a tendency that when fh=fp, the improvement should achieve the “best”. We use 5 μm step size to simulate Standard case in Figure 8 with hole radius changes from 40 μm to 80 μm; the results of SPL over the membrane and holes are shown below. As Figure 15 shows, thea=55 μm,a=60 μm, anda=65 μm cells have larger 3-dB bandwidth than other dimensions. The 3-dB fractional bandwidth are 29.4%, 35.7%, and 29.7% for a=55 μm, 60 μm, and 65 μm, respectively. The calculatedaini which let fh=fp is 55 μm; however, its fractional bandwidth enhancement is not the “best” one. For the “best” fractional bandwidth enhancement witha=60 um, as we described above, the fh is close to fp. Up to now, we cannot give a quantitative analysis for how close are the fh and fp for the “best” bandwidth performance enhancement. This still needs to be further investigated.Figure 15
FEM results of the Standard case: sound pressure level (SPL) over the membrane and holes (dB), hole radius changes from 40μm to 80 μm with 5 μm step size.


## 5. Conclusions

In this paper, we present the air-coupled unsealed Helmholtz structural CMUTs which utilize both Helmholtz resonant effect and squeeze film effect to enhance the device’s output pressure bandwidth performance in transmit mode. We have developed an analytical model to explain the vibration coupling between the vibration of membrane and the vibration of air pistons in the holes, which is the basic working mechanism of the Helmholtz structural CMUTs. In order to validate the analytical model, we use FEM model to simulate three cases of the 4-hole unsealed Helmholtz structural CMUT cells. From the FEM results, for the Double case, the two resonant peaks in the frequency response clearly show the vibration coupling working mechanism of the Helmholtz structural CMUT cell. For an unsealed Helmholtz structural CMUT cell with various hole radii, compared with the conventional sealed one, the improvements of 3-dB fractional SPL bandwidth and SPL-bandwidth product are around 35 times in the Standard case and 17 times in the Half case, which are achieved with the hole radius which letfh being close to fp. Compared with FEM results, within reasonable deviations, the analytical model can predict the characteristics such as U, SPL andf+ and f- of the unsealed Helmholtz structural CMUTs. In the design phase of Helmholtz structural CMUTs, this analytical model can help to initially select the cell dimensions with output pressure bandwidth performance enhancement, thus simplifying the design process and shortening the design period. The smallerh/2σvisc ratio or the more holes under the same hole area ratio, the weaker the peak majorly contributed by air pistons, which means the weaker Helmholtz resonant effect. Our analytical model loses efficacy when the Helmholtz effect tends to disappear, and the exact effective range of our analytical model still needs to be further investigated. In summary, the air-coupled 4-hole Helmholtz structural CMUT cell, which is designed based on the vibration coupling mechanism, can greatly enhance the output pressure bandwidth and keep the output pressure at a reasonable level. The analytical model can simplify the design process of this kind of CMUTs.



---

*Source: 3896965-2019-04-17.xml*