Research Article  Open Access
PullIn Analysis of the Flat Circular CMUT Cell Featuring Sealed Cavity
Abstract
Capacitive micromachined ultrasonic transducers (CMUTs) are one of the appealing MEMS devices. Most studies treat CMUTs as rigid plates vibrating in open air, ignoring the mechanical boundary conditions for simplification and resulting in cumulative errors in coupled fields. This paper presents a new analytical model for the pullin characteristics of the flat circular CMUT cell featuring sealed cavity. Utilizing the plate theory coupled with Boyle’s law, the paper establishes a strong relation between the pressures inside the sealed cavity and the pullin characteristics for the first time. Not only did we point out that the existence of the pressure inside the sealed cavity cannot be omitted, but we also quantified the direct effect of the pressure ratios on the pullin phenomenon. The pullin voltages increase while the pullin ratios decrease with the pressure ratios of the pressure inside the sealed cavity to the ambient pressure. The proposed calculation process delivers a good approximation of the pullin voltages and displacements, which are consistent with COMSOL simulation results. Particularly, the percentage error of our calculation process is 6.986% for the worst case. Therefore, our proposed analytical model accurately and efficiently predicts the pullin characteristics and this paper offers new perspectives and reference value in designing and modeling the CMUTs.
1. Introduction
Research on MEMS (microelectromechanical systems) has seen an amazing growth during the last three decades, stimulated both by their interesting physical properties and their attractive application potential. One of those appealing MEMS devices is capacitive micromachined ultrasonic transducers (CMUTs) [1, 2]. Advantages of CMUTs over conventional piezoelectric transducers include better impedance matching with fluids, broader bandwidth, stable device properties, the ease of fabrication, and the possibility of onchip integration with electronics [3]. Thus CMUTs are turning out to be a promising alternative in the fields of medical ultrasonic imaging, flow measurement, and aircoupled nondestructive applications [4].
Basically, the CMUT converts electrical energy into mechanical energy and vice versa [5]. Throughout the electromechanical coupling process, the pullin characteristics are of great significance [6]. Different analyses and modeling approaches can be found in the literatures [7, 8]. Most studies [5, 9, 10] treat CMUTs as rigid plates in open air in their analysis, with a uniform velocity profile across the CMUT surface. Combining the electrostatic force and the elastic force, the approximate pullin voltage could be obtained. This assumption, first issued by [5], was adopted because of its ease in modeling and the similarity to the theory used for piezoelectric transducers. However, considering the relative small dimensions of CMUTs, the cumulative effect of errors in coupled fields may lead to a huge deviation, and the reasons are as follows.
Firstly, the model does not take the top electrode into consideration for the static calculation, and the flexural rigidity depends on the dimensions and properties of the membrane only. Unlike the singlelayer plate, the flexural rigidity of the moving parts depends on both the top electrode and the membrane, resulting in a homogenized flexural rigidity. Secondly, the model does not take into account the mechanical boundary conditions of the fixed circumference and simplifies the moving parts to be a uniform plate with exactly the same deflection and velocity throughout. Thirdly, the model neglects the significant effects of the pressure inside the sealed cavity. Once the cavity is hermetically sealed during the fabrication, the pressure inside is set. Under the external pressure and the biased voltage, the thin film inside the cavity will be squeezed and have an impact on the pullin characteristics.
Our aim was to build up a calculation process of pullin characteristics considering the effects of the top electrode, the fixed boundary conditions, and the pressure inside the sealed cavity, which would be accurate enough for designing while avoiding a massive electromechanical Finite Element Modeling (FEM). Combining the plate theory featuring fixed circumference and Boyle’s law inside the sealed cavity, this paper, for the first time, obtained a relation of the pullin characteristics with respect to the pressure ratios of the pressures inside the cavity to the ambient pressure.
Section 1 introduces the device and issues the main problems with former research. Details of the development of pullin analysis are provided. Section 2 takes the top electrode into consideration and treats the CMUT as a flexible thin plate. Based on the fixed boundary conditions and Boyle’s law, the mathematical descriptions of pullin characteristics are obtained. In Section 3, the theoretical results provided by this model are compared with FEM results using COMSOL Multiphysics. The discussion involves an investigation of the key features of pullin characteristics, including both pullin ratios and pullin voltages. Section 4 concludes the paper and assesses the regime of validity and limitations of the present model.
2. Problem Formulation and Solution
The typical CMUT cell is built with a circular, square, or hexagonal membrane separated from a fixed substrate by a small air gap [11]. The circular flat CMUT cell is chosen for its prevalence. Several thousand such cells electrically connected in parallel to form a CMUT array. The geometry of the CMUT cell is shown as Figure 1. The moving part of the CMUT cell is made from a silicon nitride membrane covered with a gold electrode [12]. A layer of silicon dioxide is deposited in order to prevent the electrical shortcut between the two electrodes. The highly doped silicon substrate is also utilized as the bottom electrode. Using waferbonding technology [13, 14], the cavity and the membrane are defined on separate wafers and then are bonded under vacuum conditions.
CMUTs are efficient transmitters and receivers for aircoupled nondestructive evaluation applications, generating or detecting ultrasonic waves by vibrating membrane featuring fixed circumference, respectively [15]. When a DC voltage is applied between the two electrodes, the membrane is attracted to the substrate by electrostatic forces. When an AC voltage is superimposed over the DC voltage, the membrane will move in response to the applied signal and an ultrasonic wave is generated and launched into the ambient fluid. Similarly, when an ultrasonic wave approaches the membrane kept under a DC bias, it will vibrate and an output AC voltage can be measured from the CMUT electrodes. The applied DC voltage is necessary because it supplies the charge that is modulated by the membrane movement to create an AC signal.
2.1. Mechanical Modeling of a CMUT Cell
Taking the electrode into consideration, the flexural rigidity depends on the properties and dimensions of both the membrane and the electrode [16, 17]. Thus, the flexural rigidity is given bywhere and are Young’s Modulus of the membrane and the top electrode, respectively. and are their thicknesses. and are their Poisson’s ratios. is the location of the median plane and is the location of the upper surface of the th layer, as shown in Figure 2.
Together with the top electrode and the membrane, the moving part of the CMUT cell could be taken as a flexible thin plate with radius . The equation governing deflection of a thin flexible circular plate under a uniform loading is given asIntegrating several times [18] leads to where are constants that can be evaluated from the boundary conditions.
The CMUT cell is assumed to be clamped on the circumference; thus the boundary conditions are
Taking the boundary conditions (4) into (3) leads to
The vertical deflection is obtained as and .
Over the entire plate area, the average deflection could be defined as
2.2. Electrical Modeling of a CMUT Cell
As shown in Figure 1, the membrane and the insulation layer make up a capacitor in series with the gap capacitance. Therefore, the effective gap distance is expressed aswhere is the thickness of the insulation layer and the relative permittivity of the silicon dioxide.
When a DC voltage is applied between the electrodes, as shown in Figure 3, the electrical capacitance inside the sealed cavity could be calculated as where is the permittivity in the vacuum gap and = 8.85 × 10^{−12} F/m. For simplicity, we denote . Then (9) could be rewritten as
The first derivative of with respect to is given by
The electrical force on the membrane caused by an applied voltage is given as [19]
2.3. Calculation of PullIn Characteristics
Usually, the cavity between the membrane and the substrate is hermetically sealed during the fabrication, and the pressure inside the cavity has great effect on the deviation of the pullin characteristics [13].
Suppose the ambient pressure is . The initial pressure inside the cavity is and . Before the DC voltage is applied between two electrodes, the moving part has an initial deflection as shown in Figure 4(a). The equilibrium equation is
(a)
(b)
After the DC voltage is applied, the equilibrium equation changes towhere is the pressure inside the cavity after the DC voltage, as shown in Figure 4(b).
Based on the analysis above, the process before and after the biased voltage should be isothermal yet not adiabatic. According to Boyle’s law, the product of the pressure inside the cavity and the volume of the cavity are a constant; that is,
Combining (7) and (13), the spring coefficient for the circular plate with fixed circumference could be obtained as
When a biased voltage is applied between the CMUT’s electrodes, regardless of its polarity, it deflects the thin plate towards the bottom electrode. If the applied voltage is approaching the pullin voltage, the plate will reach a stable state eventually. After several steps of calculation and iteration, the relation among the biased voltage , the ratio of the central deflection to the effective gap distance, and the ratio of the pressure inside the cavity to the ambient pressure could be obtained as
This expression points out the relations among , , and . Under every pressure ratio , the equation has two digital roots. As the voltage increases, two solutions approach each other and coalesce together, as in Figure 5. This is the pullin point where the pullin phenomenon occurs and the pullin characteristics could be obtained. In Figure 5, since the pressure inside the cavity is not zero (), there will be initial vertical deflections; thus the ratio does not start from 0.
This paper focuses on the pullin analysis of the flat circular CMUT cell featuring sealed cavity. In Section 3, in order to verify the new model for CMUTs without biased voltage, (6) and (7) will be used to calculate the maximum and average deflections of the CMUT’s moving part as a flexible thin plate. Equation (17) will be used to calculate the pullin characteristics under different pressure ratios for the validation of the proposed analytical model.
3. Validation and Discussion
For the validation of the present analytical model for the pullin characteristics, this section is organized as follows.
Firstly, in order to verify the accuracy of the new model in mechanical field only, the maximum deflections and the average deflections without biased voltage have been analyzed and modeled by the analytical model and FEM, respectively. For every pressure ratio , and obtained from the new model are compared with the FEM results.
Secondly, for the validation and accuracy of the new model in electromechanical coupled fields, the pullin characteristics for different pressures inside the sealed cavity are calculated and modeled by the analytical model and FEM, respectively. For every pressure ratio , the pullin voltage and the pullin ratio predicted by the new model are compared with the FEM simulations. Besides, the pullin characteristics obtained from the analytical model in [5] are also used for comparison.
The dimensions of the CMUT cell used in this section are listed in Table 1. The material properties needed for calculations and simulations are listed in Table 2.


Featuring small dimensions, the MEMS devices often involve complex interactions in coupled domains, and FEM is commonly used to deal with such nonlinear models. In this section, FEM simulations are taken as reference, with whom the analytical results are compared for the model’s validation and accuracy.
We constructed a Finite Element Model of a CMUT cell in COMSOL Multiphysics (COMSOL Inc., Burlington, MA, USA) software, coupling the structural mechanics subdomain and the electrostatics subdomain. The 2D axisymmetric electromechanical coupling model is used, as shown in Figure 6. The free triangular mesh is used and the minimum mesh size is 0.00625 μm as illustrated in Figure 7. The average mesh quality factor is 0.9274, which is high enough for accurate simulations. As the bottom electrode, the fixed substrate of CMUT does not create any displacements, so it could be omitted when modeling. Instead, the interface between the bottom electrode and the insulation layer could be set for zero voltage for the electrostatics subdomain. The CPU of the workstation for FEM simulations in this paper is Inter(R) Xeon(R) CPU E52680 V3 @2.5 GHz, and the memory is 192 GB.
3.1. Displacement Comparison without Biased Voltage
Featuring the aforementioned dimensions, the CMUT cell is simulated under different pressure inside the sealed cavity. Without biased voltage, still, the membrane together with the top electrode would deflect due to the difference between the ambient pressure and the pressure inside the sealed cavity, as shown in Figure 8.
The maximum deflections and the average deflections for different pressures inside the sealed cavity are calculated by (6), (7) and simulated by FEM. The deflections and the relative errors between the analytical and FEM simulations are demonstrated in Figures 9 and 10. The black line stands for the data given by the new model, the red line stands for the FEM results, and the blue line represents the relative errors between the new model and the FEM. Relevant data are listed in Table 3.

From Figures 9 and 10, it can be obtained that the maximum deflections and the average deflections decrease as the pressure ratio increases. Since is the ratio of the pressure inside the sealed cavity to the ambient pressure (), it also represents the pressure difference across the moving part. When the pressure inside the sealed cavity is smaller ( is smaller), the deflections are bigger. Therefore, the moving part of the CMUT cell has bigger initial deflections with a smaller .
Together with Table 3, it can be observed that, for the mechanical field only, the new model could accurately predict the vertical deflections of the CMUT’s moving part without biased voltage. The analytical predictions are in great agreement with the FEM results, yet the analytical deflections are slightly bigger. Also, the FEM results are taken as reference. The relative errors between the new model and FEM results for the maximum deflections are 0.850%~2.964%, and the relative errors for the average deflections are 0.006%~2.073%.
3.2. PullIn Comparison under Different Pressure Ratios
3.2.1. PullIn Ratios Comparison
To perform the parametric analysis in coupled fields, a range of pressure ratios are applied and the model is run through a sequence of static analyses for different parameter values. Every sequence follows through a certain number of iterations and the convergence is inspected. The process is repeated and the voltage is incremented until the membrane is about to contact the insulation layer. At this point, the solution ceases to converge and the simulation terminates. The last voltage value before divergence of the solution is the pullin voltage. When the pullin voltage is applied to the electrodes, we call the then to be the pullin gap and the then to be the pullin ratio. Under different pressure ratios, the pullin ratios and the relative errors are shown in Figure 11. Relevant data are listed in Table 4.

In Figure 11, the black line stands for the data given by our new model, the red line stands for FEM results, and the blue line represents the relative errors between the new model and FEM.
From Figure 11, it can be obtained that pullin ratios decrease as the pressure ratios increase for both analytical and simulation results. When the pressure inside the sealed cavity is bigger ( is bigger), the pullin ratio is smaller, indicating that the moving part of the CMUT cell is harder to approach the pullin point with a bigger pressure inside the sealed cavity. Also, the analytical results are in great agreement with FEM results, producing the relative errors ranging from −6.986% to 1.612% (Table 4). As a matter of fact, FEMs like COMSOL are precise yet require more computing power and need more time to get the results, while the proposed new model is just several steps of calculations and more efficient to get the accurate results with acceptable relative errors. However, when compared to the relative errors produced in mechanical field only, the relative errors in electromechanical coupled fields are slightly bigger, indicating that the nonlinearity of the electrical force has brought bigger errors.
Notice that the tendency for analytical results is basically linear from to , while the simulations show a smaller slope after . From to , the pullin ratios given by the new model decrease from 0.502 to 0.466, yet the FEM results stay around 0.501. This nonlinearity is maybe due to the electromechanical coupling effect aroused by the applied voltage. When the pressure inside the cavity is bigger than a certain value, the electromechanical coupling effect has bigger influence on the pullin ratios, resulting in a series of similar pullin ratios after that. This coupling effect also leads to bigger relative errors between the new model and FEM, which needs to be taken into consideration afterwards.
3.2.2. PullIn Voltages Comparison
Commonly, when treated as rigid plates in open air, the CMUT cell has a uniform velocity profile across the surface. Combining the electrostatic force and the elastic force, the approximate pullin voltage could be obtained. A case in point is the frequently cited literature [5]. Referring to [5], the pullin voltage under such assumption is given by , without consideration of the pressure ratios. For the dimensions and material properties listed in Tables 1 and 2, the pullin voltage calculated by [5] is about 266.773 V. The pullin voltages obtained from our new model, the FEM, and [5] under different pressure ratios are all shown in Figure 12.
Relevant data are listed in Table 4.
In Figure 12, the black line stands for the data given by our new model, the red line stands for FEM results, the blue line represents the relative errors between the new model and FEM, and the pink line represents the pullin voltage calculated by [5].
From Figure 12, it can be observed that pullin voltages obtained from the new model and the FEM simulations both increase along with the pressure ratios, while the pullin voltage calculated by [5] does not change with since the pressures inside the sealed cavity are not taken into consideration. The analytical results predicted by the new model are in great agreement with FEM results, producing the relative errors ranging from −2.410% to 5.028% (Table 4). However, the pullin voltage calculated by [5] is much larger than them. The analytical model used in [5] takes the CMUT as rigid plates in open air, resulting in the invariable and overestimated pullin voltage that could lead to improper characterization and experimental preparations. Again, notice that the tendency of the analytical results predicted by our new model is basically linear from to , while the FEM simulations show fluctuations around the analytical results. Thus, the aforementioned electromechanical coupling effect does need to be taken into consideration for the accuracy of the present model.
4. Conclusion
This paper proposed a new analytical model for the pullin characteristics of the flat circular CMUT cell featuring sealed cavity based on the plate theory and Boyle’s law. Not only did we point out that the existence of the pressure inside the sealed cavity cannot be omitted, but we also quantified the direct effect of the pressure ratios on the pullin phenomenon. The computational results of our analytical model were compared with the COMSOL simulations. The good agreement between them shows that our analytical model can approximate the CMUT’s behavior fairly well. The work in this paper is summarized as follows.
(1) An improved analytical model for the pullin characteristics of the flat circular CMUT cell featuring the sealed cavity has been presented. The strong relation between the pressures inside the sealed cavity and the pullin characteristics is established for the first time. The derivation of this new model takes the top electrode into consideration and treats the CMUT as a flexible plate with fixed boundary conditions. The relations of the biased voltage , the ratio of the central deflection to the effective gap distance, and the ratio of the pressure inside the cavity to the ambient pressure are obtained.
(2) The validation and accuracy of the new model in mechanical field only have been obtained. The maximum and average deflections without biased voltage have been calculated by the new model and simulated by FEM. The maximum relative error is 2.964% for the maximum deflection and 2.073% for the average deflection, respectively. Thus, the mechanical analysis of the new model could accurately and efficiently predict the CMUT’s vertical deflections.
(3) The validation and accuracy of the new model in coupled fields have been obtained. The pullin characteristics under different pressure ratios have been verified by both the new model and FEM simulations. The maximum relative error is 6.986% for the pullin ratios and 5.028% for the pullin voltages, respectively. In particular, the comparison of pullin voltages among the present model, the FEM simulations, and the literature [5] that treats the CMUT as rigid plates leads to the conclusion that the fixed boundary conditions of the flexible thin plate are crucial to the model’s validation and accuracy. Also the comparison does point out that the consideration of the pressures inside the sealed cavity is significant and necessary for the designing and modeling of the CMUT’s pullin analysis.
(4) Although our proposed process is quite precise during the analysis and validation, it has some limitations to be improved. The complex coupling nonlinearity causes bigger relative errors between the analytical pullin characteristics and FEM results. In order to make the analytical model deliver smaller relative errors, the electromechanical coupling effect should be taken into consideration afterwards.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work has been supported by the Young Scientists Fund of the National Natural Science Foundation of China (Grant no. 61201039).
References
 Ö. Oralkan, A. S. Ergun, J. A. Johnson et al., “Capacitive micromachined ultrasonic transducers: nextgeneration arrays for acoustic imaging?” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 49, no. 11, pp. 1596–1610, 2002. View at: Publisher Site  Google Scholar
 D. Gross, C. Coutier, M. Legros, A. Bouakaz, and D. Certon, “A cMUT probe for ultrasoundguided focused ultrasound targeted therapy,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 62, no. 6, pp. 1145–1160, 2015. View at: Publisher Site  Google Scholar
 T. A. Emadi and D. A. Buchanan, “Design and characterization of a capacitive micromachined tranducer with a defectable bottom electrode,” IEEE Electron Device Letters, vol. 36, no. 6, pp. 612–614, 2015. View at: Google Scholar
 B. T. KhuriYakub and Ö. Oralkan, “Capacitive micromachined ultrasonic transducers for medical imaging and therapy,” Journal of Micromechanics and Microengineering, vol. 21, no. 5, Article ID 054004, 2011. View at: Publisher Site  Google Scholar
 G. G. Yaralioglu, A. S. Ergun, B. Bayram, E. Hæggström, and B. T. KhuriYakub, “Calculation and measurement of electromechanical coupling coefficient of capacitive micromachined ultrasonic transducers,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 50, no. 4, pp. 449–456, 2003. View at: Publisher Site  Google Scholar
 B. Ahmad and R. Pratap, “Elastoelectrostatic analysis of circular microplates used in capacitive micromachined ultrasonic transducers,” IEEE Sensors Journal, vol. 10, no. 11, pp. 1767–1773, 2010. View at: Publisher Site  Google Scholar
 E. Aydoǧdu, A. Ozgurluk, A. Atalar, and H. Köymen, “Parametric nonlinear lumped element model for circular CMUTs in collapsed mode,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 61, no. 1, pp. 173–181, 2014. View at: Publisher Site  Google Scholar
 N. Apte, K. K. Park, and B. T. KhuriYakub, “Finite element analysis of CMUTs with pressurized cavities,” in Proceedings of the IEEE International Ultrasonics Symposium (IUS '12), pp. 979–982, IEEE, Dresden, Germany, October 2012. View at: Publisher Site  Google Scholar
 B. Ahmad and R. Pratap, “Analytical evaluation of squeeze film forces in a CMUT with sealed airfilled cavity,” IEEE Sensors Journal, vol. 11, no. 10, pp. 2426–2431, 2011. View at: Publisher Site  Google Scholar
 I. O. Wygant, M. Kupnik, and B. T. KhuriYakub, “Analytically calculating membrane displacement and the equivalent circuit model of a circular CMUT cell,” in Proceedings of the IEEE International Ultrasonics Symposium (IUS '08), pp. 2111–2114, Beijing, China, November 2008. View at: Publisher Site  Google Scholar
 C. B. Doody, X. Y. Cheng, C. A. Rich, D. F. Lemmerhirt, and R. D. White, “Modeling and characterization of CMOSfabricated capacitive micromachined ultrasound transducers,” Journal of Microelectromechanical Systems, vol. 20, no. 1, pp. 104–118, 2011. View at: Publisher Site  Google Scholar
 M. L. Kuntzman, D. Kim, and N. A. Hall, “Microfabricated and experimental evaluation of a rotational capacitive micromachined ultrasonic transducers,” Journal of Microelectromechanical Systems, vol. 24, no. 2, pp. 404–413, 2015. View at: Publisher Site  Google Scholar
 Y. Huang, A. Sanli Ergun, E. Hæggström, M. H. Badi, and B. T. KhuriYakub, “Fabricating capacitive micromachined ultrasonic transducers with waferbonding technology,” Journal of Microelectromechanical Systems, vol. 12, no. 2, pp. 128–137, 2003. View at: Publisher Site  Google Scholar
 K. K. Park, H. Lee, M. Kupnik, and B. T. KhuriYakub, “Fabrication of capacitive micromachined ultrasonic transducers via local oxidation and direct wafer bonding,” Journal of Microelectromechanical Systems, vol. 20, no. 1, pp. 95–103, 2011. View at: Publisher Site  Google Scholar
 W. Zhang, H. Zhang, Y. Wang, F. Du, S. Jin, and Z. Zeng, “Simulation characterization of CMUT with vented square membrane,” in International Conference on Optical Instruments and Technology: Micro/Nano Photonics and Fabrication, vol. 9624 of Proceedings of SPIE, Beijing, China, May 2015. View at: Publisher Site  Google Scholar
 S. P. Timoshenko and S. WoinowskyKrieger, Theory of Plates and Shells, McGrawHill, New York, NY, USA, 1959.
 D. Certon, F. Teston, and F. Patat, “A finite difference model for CMUT devices,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 52, no. 12, pp. 2199–2210, 2005. View at: Publisher Site  Google Scholar
 E. Ventsel and T. Krauthammer, Thin Plates and Shells, Marcel Deker, New York, NY, USA, 1st edition, 2001.
 Z. Wang, Microsystem Design and Fabrication, Tsinghua University Press, Beijing, China, 2008.
Copyright
Copyright © 2015 Wen Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.