Fluid damping arising from squeeze film flow of air or some inert gas trapped between an oscillating micro mechanical structure, such as a beam or a plate, and a fixed substrate often dominates the other energy dissipation mechanisms in silicon based MEM devices. As a consequence, it has maximum effect on the resonant response or dynamic response of the device. Unfortunately, modelling of the squeeze film flow in most MEMS devices is quite complex because of several factors unique to MEMS structures. In this thesis, we set out to study the effect of these factors on squeeze film flow. First we list these factors and study each of them in the context of a particular application, using experimental measurements, extensive numerical simulations, and analytical modelling for all chosen factors. We consider five important factors. The most important factor perhaps is the effect of rarefaction that is dominant when a device is vacuum packed with low to moderate vacuum, typical for MEMS packaging. The second problem is to investigate and model the effect of perforations which are usually provided for efficient etching of the sacrificial layer during fabrication of the suspended structures. The third problem is to consider the effect of non-uniform deflection of the structure such as those in MEMS cantilever beams and model its effect on the squeeze film. The fourth effect studied is the influence of different boundary conditions such as simple, fully open and partially closed boundaries around the vibrating structure on the characteristics of the squeeze film flow. The fifth problem undertaken is to analyze the effect of high operating frequencies on the squeeze film damping. In the first problem, the rarefaction effect is studied by performing experiments under varying pressures. Depending on the ambient pressure or the size of the gap between the vibrating and the fixed structure, the fluid flow may fall in any of the flow regimes, ranging from continuum flow to molecular flow, and giving a wide range of dissipation. The relevant fluid flow characteristics are determined by the Knudsen number, which is the ratio of the mean free path of the gas molecule to the characteristic flow length of the device. This number is very small for continuum flow and reasonably big for molecular flow. Here, we study the effect of fluid pressure on the squeeze film damping by carrying out experiments on a MEMS device that consists of a double gimbaled torsional resonator. Such devices are commonly used in optical cross-connects and switches. We vary fluid pressure to make the Knudsen number go through the entire range of continuum flow, slip flow, transition flow, and molecular flow. We experimentally determine the quality factor of the torsional resonator at different air pressures ranging from 760 torr to 0.001 torr. The variation of this pressure over six orders of magnitude ensures the required rarefaction to range over all flow conditions. Finally, we get the variation of the quality factor with pressure. The result indicates that the quality factor, Q, follows a power law, Q P-r, with different values of the exponent r in different flow regimes. To numerically model the effect of rarefaction, we propose the use of effective viscosity in Navier-Stokes equation. This concept is validated with analytical results for a simple case. It is then compared with the experimental results presented in this thesis. The study shows that the effective viscosity concept can be used effectively even for the molecular regime if the air-gap to length ratio is sufficiently small (h0/L < 0.01). However, as this ratio increases, the range of validity decreases. Next, a semianalytical approach is presented to model the rarefaction effect in double-gimballed MEMS torsion mirror. In this device, the air gap thickness is 80 µm which is comparable to the lateral dimension 400 µm of the oscillating plate and thus giving the air-gap to length ratio of 0.2. As the ratio 0.2 is much greater than 0.01, the conventional Reynolds equation cannot be used to compute the squeeze film damping. Consequently, we find the effective length of an equivalent simple mirror corresponding to the motion about the two axes of the mirror such that the Reynolds equation still holds. After finding the effective length, we model the rarefaction effect by incorporating effective viscosity which is based on different models including the one proposed in this paper. Then we compare the analytical solution with the experimental result and find that the proposed model not only captures the rarefaction effect in the slip, transition and molecular regimes but also couples well with the non-fluid damping in the intrinsic regime. For the second problem, several analytical models exist for evaluating squeeze film damping in rigid rectangular perforated MEMS structures. These models vary in their treatment of losses through perforations and squeezed film, in their assumptions of compressibility, rarefaction and inertia, and their treatment of various second order corrections. We present a model that improves upon previously reported models by incorporating more accurate losses through holes proposed by Veijola and treating boundary cells and interior cells differently as proposed by Mohite et al. The proposed model is governed by a modified Reynolds equation that includes compressibility and rarefaction effect. This equation is linearized and transformed to the standard two-dimensional diffusion equation using a simple mapping function. The analytical solution is then obtained using Green’s function. The solution thus obtained adds an additional term Γ to the damping and spring force expressions derived by Blech for compressible squeeze flow through non-perforated plates. This additional term contains several parameters related to perforations and rarefaction. Setting Γ = 0, one recovers Blech’s formulas. We benchmark all the models against experimental results obtained for a typical perforated MEMS structure with geometric parameters (e.g., perforation geometry, air gap, plate thickness) that fall well within the acceptable range of parameters for these models (with the sole exception of Blech’s model that does not include perforations but is included for historical reasons). We compare the results and discuss the sources of errors. We show that the proposed model gives the best result by predicting the damping constant within 10% of the experimental value. The approximate limit of maximum frequencies under which the formulas give reasonable results is also discussed. In the third problem, we study the effect of elastic modeshape during vibration on the squeeze film flow. We present an analytical model that gives the values of squeeze film damping and spring coefficients for MEMS cantilever resonators taking into account the effect of flexural modes of the resonator. We use the exact modeshapes of a 2D cantilever plate to solve for pressure in the squeeze film and then derive the equivalent damping and spring coefficient relations from the back force calculations. The relations thus obtained can be used for any flexural mode of vibration of the resonators. We validate the analytical formulas by comparing the results with numerical simulations carried out using coupled finite element analysis in ANSYS, as well as experimentally measured values from MEMS cantilever resonators of various sizes and vibrating in different modes. The analytically predicted values of damping are, in the worst case, within less than 10% of the values obtained experimentally or numerically. We also compare the results with previously reported analytical formulas based on approximate flexural modeshapes and show that the proposed model gives much better estimates of the squeeze film damping. From the analytical model presented here, we find that the squeeze film damping drops by 84% from the first mode to the second mode in a cantilever resonator, thus improving the quality factor by a factor of six to seven. This result has significant implications in using cantilever resonators for mass detection where a significant increase in quality factor is obtained only by using vacuum. In the fourth and fifth problem, the effects of partially blocked boundary condition and high operating frequencies on squeeze films are studied in a MEMS torsion mirror with different boundary conditions. For the structures with narrow air-gap, Reynolds equation is used for calculating squeeze film damping, generally with zero pressure boundary conditions on the side walls. This procedure, however, fails to give satisfactory results for structures under two important conditions: (a) for an air-gap thickness comparable to the lateral dimensions of the micro structure, and (b) for non-trivial pressure boundary conditions such as fully open boundaries on an extended substrate or partially blocked boundaries that provide side clearance to the fluid flow. Several formulas exist to account for simple boundary conditions. In practice, however, there are many micromechanical structures, such as torsional MEMS structures, that have non-trivial boundary conditions arising from partially blocked boundaries. The most common example is the double-gimballed MEMS torsion mirror of rectangular, circular, or hexagonal shape. Such boundaries usually have clearance parameters that can vary due to fabrication. These parameters, however, can also be used as design parameters if we understand their role on the dynamics of the structure. We take a MEMS torsion mirror as an example device that has large air-gap and partially blocked boundaries due to static frames. Next we model the same structure in ANSYS and carry out CFD (computational fluid dynamics) analysis to evaluate the stiffness constant K, the damping constant C, as well as the quality factor Q due to the squeeze film. We compare the computational results with experimental results and show that without taking care of the partially blocked boundaries properly in the computational model, we get unacceptably large errors. Subsequently, we use the CFD calculations to study the effect of two important boundary parameters, the side clearance c, and the flow length s, that specify the partial blocking. We show the sensitivity of K and C on these boundary design parameters. The results clearly show that the effect of these parameters on K and C is substantial, especially when the frequency of excitation becomes close to resonant frequency of the oscillating fluid and high enough for inertial and compressibility effects to be significant. Later, we present a compact model to capture the effect of side boundaries on the squeeze film damping in a simple rectangular torsional structure with two sides open and two sides closed. The analytical model matches well with the numerical results. However, the proposed analytical model is limited to low operating frequencies such that the inertial effect is negligible. The emphasis of this work has been towards developing a comprehensive understanding of different significant factors on the squeeze film damping in MEMS devices. We have proposed various ways of modelling these effects, both numerically as well as analytically, and shown the efficacy of these models by comparing their predictive results with experimental results. In particular, we think that the proposed analytical models can help MEMS device designers by providing quick estimates of damping while incorporating complex effects in the squeeze film flow. The contents of the thesis may also be of interest to researchers working in the area of microfluidics and nanofluidics.

The fast-paced growth in microelectromechanical systems (MEMS), microfluidic fabrication, porous media applications, biomedical assemblies, space propulsion, and vacuum technology demands accurate and practical transport equations for rarefied gas flows. It is well-known that in rarefied situations, due to strong deviations from the continuum regime, traditional fluid models such as Navier-Stokes-Fourier (NSF) fail. The shortcoming of continuum models is rooted in nonequilibrium behavior of gas particles in miniaturized and/or low-pressure devices, where the Knudsen number (Kn) is sufficiently large. Since kinetic solutions are computationally very expensive, there has been a great desire to develop macroscopic transport equations for dilute gas flows, and as a result, several sets of extended equations are proposed for gas flow in nonequilibrium states. However, applications of many of these extended equations are limited due to their instabilities and/or the absence of suitable boundary conditions. In this work, we concentrate on regularized 13-moment (R13) equations, which are a set of macroscopic transport equations for flows in the transition regime, i.e., Kn≤1. The R13 system provides a stable set of equations in Super-Burnett order, with a great potential to be a powerful CFD tool for rarefied flow simulations at moderate Knudsen numbers. The goal of this research is to implement the R13 equations for problems of practical interest in arbitrary geometries. This is done by transformation of the R13 equations and boundary conditions into general curvilinear coordinate systems. Next steps include adaptation of the transformed equations in order to solve some of the popular test cases, i.e., shear-driven, force-driven, and temperature-driven flows in both planar and curved flow passages. It is shown that inexpensive analytical solutions of the R13 equations for the considered problems are comparable to expensive numerical solutions of the Boltzmann equation. The new results present a wide range of linear and nonlinear rarefaction effects which alter the classical flow patterns both in the bulk and near boundary regions. Among these, multiple Knudsen boundary layers (mechanocaloric heat flows) and their influence on mass and energy transfer must be highlighted. Furthermore, the phenomenon of temperature dip and Knudsen paradox in Poiseuille flow; Onsager's reciprocity relation, two-way flow pattern, and thermomolecular pressure difference in simultaneous Poiseuille and transpiration flows are described theoretically. Through comparisons it is shown that for Knudsen numbers up to 0.5 the compact R13 solutions exhibit a good agreement with expensive solutions of the Boltzmann equation.