Simulation-Based Design and Optimization of Rectangular Micro-Cantilever-Based Aerosols Mass Sensor

Micro-Cantilever (MCL) is a thin film structure that is applied for aerosol particle mass sensing. Several modifications to the rectangular MCL (length-to-width ratio, slots at the anchor, serrations at its side edges) are made to deduce the role and influence of the shape of rectangular MCL-based aerosol mass sensors and reduce gas damping. A finite element fluid-structure interaction model was used to investigate the performance of MCL. It is found that (I) the mass sensitivity and quality factor decline with the increasing of length-to-width ratio which alters the resonant frequency of the MCL. The optimum conditions, including the length-to-width ratio (σlw = 5) and resonant frequency (f0 = 540.7 kHz) of the MCL, are obtained with the constant surface area (S = 45,000 μm2) in the frequency domain ranging from 0 to 600 kHz. (II) The slots can enhance the read-out signal and bring a small Q factor drop. (III) The edge serrations on MCL significantly reduce the gas damping. The results provide a reference for the design of aerosol mass sensor, which makes it possible to develop aerosol mass sensor with high frequency, sensitivity, and quality.


Introduction
Particulate matter (PM) in air, considered to be one of the air pollutants, often results in the reduction of visibility and agricultural crop yield due to light dispersion and dissipation effect. Due to their small sizes and high surface-to-volume ratio, fine and ultrafine PM can also penetrate deeply into human bronchi and alveoli, causing/aggravating respiratory diseases Researchers have paid much attention to develop instruments and sensors for measuring aerosol mass concentration (due to the present air quality standards). The aerosol mass sensing method has evolved through two generation according to the volume of instruments [1].
Instruments in the early generation are designed for scientific studies. They are typically bulky in the package, such as standard gravimetric measurement, Beta Gauge [2], Tapered Element Oscillating Microbalance (TEOM) [3][4][5], Berner Low Pressure Impactor (ELPI) [6], and Nephelometric [7]. While promising to have great potential for reducing the damping of the MCL. Instead of time-consuming and costly MEMS fabrication process, the simulation-based method will minimize the cost and time in the design and optimization. At a small scale, resonators [18,19], gyroscopes [20], accelerometers [21,22], and actuators [23] must consider the physical phenomena involved in their operation in order to design and optimize the structure of MCL. A general theoretical model for the frequency response of an MCL executing torsional vibration in a fluid of arbitrary viscosity and density was presented [24]. However, the principal assumptions must be satisfied that the MCL length greatly exceeds its width, the amplitude of vibration is small, and the fluid is incompressible in nature. Another model of the resonant frequencies of rectangular cantilever beams immersed in fluid was presented [25], which rigorously valid for arbitrary mode number and for cases where the fluid maybe considered to be in viscid in nature. This model is directly applicable to cantilever beams of macroscopic size. The hydrodynamic loading of elastic MCLs vibrating in viscous fluids is analyzed computationally using a finite element fluid-structure interaction model [26]. The results indicate that MCL damping arises from localized fluid shear near the edges of the MCL. The MCL-based aerosols mass sensor is operated in air, which is easier to compress than liquid-phase fluids. In measurement of atmospheric or compressible fluid, as the mode number of MCL increases and passes a "coincidence point", the Q factor will finally start to decrease [27,28]. However, the result is different in the incompressible fluid [29][30][31][32]. Since the compression of the gas leads to additional acoustic damping, it was previously assumed to be negligible in the incompressible fluid [33].
In the present study, the solid MCLs and surrounding air are modeled respectively using shell elements and Thermoviscous Acoustics and Pressure Acoustic. The effects of length-to-width ratios (σ lw ) on the aerosol mass sensor with a constant thickness is critically investigated. Furthermore, it is the first time to add edge serration to MCL. Combined with a slotted structure, the read-out signal strength can be enhanced. Meanwhile, the edge serration on the MCL-based aerosols sensors can reduce the damping in air fluid. Thus, it will greatly reduce the damping limitation of the MCL in fluid detection under high frequency mode.

Principle of the MCL-Based Mass Sensors
The actuate methods of MCL are piezoelectric [34][35][36], electric [37,38], flexoelectric [39], atomic force [40], magnetic [41], ultrasound radiation force [42], and laser irradiation [43]. The methods of reading out the output signal of MCL include optics [44,45], piezoresistance [46,47], piezoelectricity [48], capacitive [49], and tunneling [50]. As an example, there is a magnetic-actuating MCL, be used as a sensing element as shown in Figure 1, which can be made of Silicon-On-Insulator (SOI). The MCL placed in the magnetic field can be actuated, when an alternating current is applied to the A and H pins. Making piezoresistance on the surface of MCL is a typical MEMS technology. For example, we can make it with Silicon-On-Insulator (SOI), which's thickness of the top layer is 5 µm, and the resistivity is 1-10 Ω·cm. Piezoresistance on the MCL obtained by boron (B)-doping of the n-type silicon. The injection energy is 60 keV, the dose is 8 × 10 14 /cm 3 , and the photoresist is used as mask. After the photoresist is removed, it is annealed at 950 • C in nitrogen atmosphere for 30 min. A piezoresistance of about 380 Ω is obtained. The output signal of the MCL is read through the output pins of the Diagonal Wheatstone Bridge (B/F and C/G). An input DC voltage, V in , is across the E and D terminals and the output in the form of voltage, V out , is measured across C/G and B/F terminals. In dynamic mode, the deflects of the vibration resulting in a change of resistance of the piezoresistor placed on MCL (R 2 &R 4 ). The output voltage can be expressed as: Sensors 2020, 20, 626 where the η is the adsorption efficiency of airborne particles, T l and Q l are the sampling/deposition time of the airborne particles and flow rate during sampling, respectively. Therefore, the mass concentration of particles can be obtained by detecting the resonant frequency shift of MCL.

Theory
The resonant frequency of MCL in vacuum can be calculated by Equation (3). The resonant frequency of MCL in air will be affected by air fluid. Each resonant frequency of MCL in vacuum can be calculated by the formula [51]: Accordingly, the frequency of the MCL can be converted into the frequency of the V out. The mass of particles loaded on the surface of the MCL can be calculated by measuring the resonant frequency shift of the MCL. The resonant behavior of MCL can be described by Hooke's law. For rectangular MCL, its effective stiffness/spring constant k eff can be obtained from the effective mass (m eff ) is: and the resonant frequency of the MCL (f 0 ) is: where E, b, W, and L are the elastic modulus, thickness, width and length of MCL respectively [51]. The π is the ratio of a circle's circumference. The ρ is the density of the MCL. When the MCL is loaded with particles, the resonant frequency of the MCL becomes f 1 The mass changes (∆m) can calculated from: The airborne particulate mass concentration (X) is: where the η is the adsorption efficiency of airborne particles, T 1 and Q 1 are the sampling/deposition time of the airborne particles and flow rate during sampling, respectively. Therefore, the mass concentration of particles can be obtained by detecting the resonant frequency shift of MCL.

Theory
The resonant frequency of MCL in vacuum can be calculated by Equation (3). The resonant frequency of MCL in air will be affected by air fluid. Each resonant frequency of MCL in vacuum can be calculated by the formula [51]: where the f n is resonant frequency corresponding to different modes, β n is coefficient related to the eigenvalue of each harmonic, β n = π(n − 0.5), n = 1, 2, . . . , n. When the MCL vibrates in the air, the resonant frequency ( f n_Air ) is [24,52]: The Q-factor determine the performance of the MCL-based aerosols mass sensor which influenced by fluid damping and the shapes. The Q-factor refers to the resolution for the MCL-based aerosols mass sensor. As we known, it can be defined as the ratio of the total energy stored in the vibration system (E m ) to the energy lost (E c ) in the vibration period. The Q-factor can be characterized as follows: Thermal and viscous damping of the MCL also needs to be considered. It is necessary to include thermal conduction effects and viscous losses explicitly in the governing equation. The governing equations for the fluctuation (wave equation or Helmholtz equation) are derived by perturbing or linearizing the Navier-Stokes equations, momentum equation, continuity equation and energy equation [53]. This is achieved by solving the thermoviscous acoustics equations in the Thermoviscous Acoustics interfaces in the COMSOL Acoustics Module.
The continuity equation for Thermoviscous Acoustics: The conservation equations for momentum and energy for acoustic perturbation: And in the fluid domains, the pressure distribution is governed by the wave equation: The dependent variables take the form: where p, u, and T are the pressure, velocity and temperature of air around MCL, respectively, and w is the angular frequency. Where α 0 is the coefficient of thermal expansion (isobaric), Q h and Q m are a possible heat source and the monopole source, respectively, and the q d is the dipole source. Variables with sign (') are acoustic variables, and variables with subscript (0) represent background mean flow quantities.

Design
To investigate the effects of σ lw , slots, and edge serrations on the detection performance, the MCLs with different σ lw and special structure are set up. In the simulation, modeled using shell elements, the MCL can be simplified to a spring structure with anchor fixed. It should be noted that the particles were considered to load on the upper surface of the MCL approximately uniformly. The mechanical properties of silicon material embedded in COMSOL were compiled in Table 1, they come from the default properties of silicon dynamic analysis in the MEMS module of COMSOL.

The Shapes of MCLs
Different shapes affect the detection performance of micro-cantilever. For rectangular shell structure, the σ lw is mainly considered. Different length to width ratios (σ lw = 1, 2, 3, 4, 5, 6,7,8,9) were set to investigate the effects of different shapes as shown in Figure 2. To eliminate the influence of the effective mass of the MCL, the thickness of the MCLs are set to 5 µm and the effective area (tunable surface area, S) is 4.5 × 10 4 µm 2 .
Sensors 2020, 20, x FOR PEER REVIEW 6 of 19 were set to investigate the effects of different shapes as shown in Figure 2. To eliminate the influence of the effective mass of the MCL, the thickness of the MCLs are set to 5 μm and the effective area (tunable surface area, S) is 4.5 × 10 4 μm 2 .

The Slots of MCLs
The rate of resistance change of the piezoresistive resistance can be described as: where the E is the elastic modulus, the µ is Poisson's ratio, and the ε is strain. High rate of resistance change is beneficial to improve the signal intensity of the detection system, especially in MEMS sensors. In dynamic mode, the strain ε is positively correlated with the deflection of MCL in piezoresistive attachment area (0 μm < L < 50 μm). In addition, the deflection is

The Slots of MCLs
The rate of resistance change of the piezoresistive resistance can be described as: where the E is the elastic modulus, the µ is Poisson's ratio, and the ε is strain. High rate of resistance change is beneficial to improve the signal intensity of the detection system, especially in MEMS sensors. In dynamic mode, the strain ε is positively correlated with the deflection of MCL in piezoresistive attachment area (0 µm < L < 50 µm). In addition, the deflection is also positively correlated with other signal readout methods. Hence, the read-out signal strength increases with the deflection. As shown in Figure 3a, three slots are arranged at the anchor end of the MCL. The length of the slots is 50 µm, the width of the slots is W/7 for each shape. The deflection, mass sensitivity, power dissipation, and Q-factor of the MCL will be investigated.

The Edge Serrations of MCLs
The edge serrations structure can change the flow state of the fluid to bring about some improvements. For example, the leading-edge serrations on the owl's wings allow him to achieve remarkably low noise gliding and flapping flights [54,55]. It is assumed that the edge serrations also reduce the disturbance to the fluid in the micro length scale. The effect of edge serrations on MCLs were investigated by simulation. As shown in Figure 4, edge serrations of 15 μm by 2 μm are added to the both free edge connected to the anchor end of MCL.

Material Property of Air
According to Newton's law of viscosity, the dynamic viscosity ( ) of air is significantly correlated with temperature, but almost independent of pressure. When the temperature T is less than 2000 K, the dynamic viscosity can be calculated by Sutherland formula [56]: where is the reference dynamic viscosity corresponding to the reference temperature ( ). Where B is a constant related to the type of gas (B of air is 110.4 K). We need to simulate thermal and viscous damping, so we need to consider the influence of temperature on the heat capacity and thermal conductivity of air. The heat capacity (C P ) of air was described by empirical formulas:

The Edge Serrations of MCLs
The edge serrations structure can change the flow state of the fluid to bring about some improvements. For example, the leading-edge serrations on the owl's wings allow him to achieve remarkably low noise gliding and flapping flights [54,55]. It is assumed that the edge serrations also reduce the disturbance to the fluid in the micro length scale. The effect of edge serrations on MCLs were investigated by simulation. As shown in Figure

The Edge Serrations of MCLs
The edge serrations structure can change the flow state of the fluid to bring about some improvements. For example, the leading-edge serrations on the owl's wings allow him to achieve remarkably low noise gliding and flapping flights [54,55]. It is assumed that the edge serrations also reduce the disturbance to the fluid in the micro length scale. The effect of edge serrations on MCLs were investigated by simulation. As shown in Figure 4, edge serrations of 15 μm by 2 μm are added to the both free edge connected to the anchor end of MCL.

Material Property of Air
According to Newton's law of viscosity, the dynamic viscosity ( ) of air is significantly correlated with temperature, but almost independent of pressure. When the temperature T is less than 2000 K, the dynamic viscosity can be calculated by Sutherland formula [56]: where is the reference dynamic viscosity corresponding to the reference temperature ( ). Where B is a constant related to the type of gas (B of air is 110.4 K). We need to simulate thermal and viscous damping, so we need to consider the influence of temperature on the heat capacity and thermal conductivity of air. The heat capacity (C P ) of air was described by empirical formulas:

Material Property of Air
According to Newton's law of viscosity, the dynamic viscosity (µ) of air is significantly correlated with temperature, but almost independent of pressure. When the temperature T is less than 2000 K, the dynamic viscosity can be calculated by Sutherland formula [56]: where µ B is the reference dynamic viscosity corresponding to the reference temperature (T 0 ). Where B is a constant related to the type of gas (B of air is 110.4 K). We need to simulate thermal and viscous damping, so we need to consider the influence of temperature on the heat capacity and thermal conductivity of air. The heat capacity (C P ) of air was described by empirical formulas: Sensors 2020, 20, 626 where a, b, c, d, and e are empirical constants: a = 1047.63657 J/(kg·K); b = −0.372589265 J/(kg·K 3 ); c = 9.45304214 × 10 −4 J/(kg·K 3 ); d = −6.02409443 × 10 −7 J/(kg·K 4 ); e = 1.2858961 × 10 −10 J/(kg·K 5 ). Similarly, the thermal conductivity (λ) of air can be obtained by empirical formula: .The relationship between air density, temperature and pressure can be described by the Clapeyron equation [57]: where ρ t is the density of air corresponding to the temperature (T), P is pressure. Where k and R constants related to the type of gas (k of air is 0.02897, R c of air is about 8314 Pa·L/mol·K). In the flow field where the temperature T is not constant, the sound velocity at each point is different. The expressions of sound velocity (u) and temperature are as follows: where λ 0 and R are the specific heat ratio and the universal gas constant (λ 0 of air is 1.4, R of air is 287.14 J ·kg − 1 ·K − 1 ).

Modeling and Validation
The model is established and calculated by using the COMSOL Multiphysics 5.3a. A sketch of MCL modeled are shown in Figures 2-4. The solid MCLs are modeled using shell elements and the surrounding air is modeled using Thermoviscous Acoustics and Pressure Acoustic [58]. When vibration waves propagate in a fluid bounded by the surface of MCL, so-called viscous and thermal boundary layers are created at the MCL surfaces. The temperature of the atmosphere is 20 • C, and the pressure is standard atmospheric pressure (1 atm). There is a typical spherical geometry is taken as an approximation of the actual geometry [26,53]. The MCL is submerged in a 3D enclosure of Navier-Stokes elements of an incompressible viscous fluid which is set as the "Thermoviscous Acoustics, frequency domain" boundary (TA). These Navier-Stokes fluid elements are then surrounded by Pressure Acoustic 3D fluid elements which is set as the "Pressure Acoustic, frequency domain" boundary (ACPR). For the case of MCL submerged in air of infinite extent and ensure the convergence of the model [26], as shown in Figure 5, the air geometry is plotted as a sphere with a diameter of 800 µm. The central sphere with an inner diameter of 500 µm. field where the temperature T is not constant, the sound velocity at each point is different. The expressions of sound velocity (u) and temperature are as follows: where λ 0 and R are the specific heat ratio and the universal gas constant (λ 0 of air is 1.4, R of air is 287.14 J· kg ·K ).

Modeling and Validation
The model is established and calculated by using the COMSOL Multiphysics 5.3a. A sketch of MCL modeled are shown in Figures 2-4. The solid MCLs are modeled using shell elements and the surrounding air is modeled using Thermoviscous Acoustics and Pressure Acoustic [58]. When vibration waves propagate in a fluid bounded by the surface of MCL, so-called viscous and thermal boundary layers are created at the MCL surfaces. The temperature of the atmosphere is 20 °C, and the pressure is standard atmospheric pressure (1 atm). There is a typical spherical geometry is taken as an approximation of the actual geometry [26,53]. The MCL is submerged in a 3D enclosure of Navier-Stokes elements of an incompressible viscous fluid which is set as the "Thermoviscous Acoustics, frequency domain" boundary (TA). These Navier-Stokes fluid elements are then surrounded by Pressure Acoustic 3D fluid elements which is set as the "Pressure Acoustic, frequency domain" boundary (ACPR). For the case of MCL submerged in air of infinite extent and ensure the convergence of the model [26], as shown in Figure 5, the air geometry is plotted as a sphere with a diameter of 800 μm. The central sphere with an inner diameter of 500 μm. The validation data are used to evaluate the simulation method based on the practical example of these fabricated rectangular MCL for detecting airborne particles. The performance of these fabricated MCLs with different dimensions was reported. We compared the experimental data with the theoretical model simulation data from this work. As the goal was to verify the accuracy of the simulation model, the focus of the validation was on the basic performance of MCL in the air at the micro scale rather than other performance aspects. The validation data are used to evaluate the simulation method based on the practical example of these fabricated rectangular MCL for detecting airborne particles. The performance of these fabricated MCLs with different dimensions was reported. We compared the experimental data with the theoretical model simulation data from this work. As the goal was to verify the accuracy of the simulation model, the focus of the validation was on the basic performance of MCL in the air at the micro scale rather than other performance aspects.

Simulation-Model Validation
To verify the accuracy of the simulation model, the model of this work is used to simulate the basic performance of MCLs as listed in Table 2. The data from experimental and simulation are compared. To verify the validity of the method and the obtained results of the theoretical model. As expected, most of the experimental and simulation data have minor differences (<10%). It is noted that there is only relative error more than 20% (24%) in the Q factor of MCL with dimension 2750 µm × 100 µm × 50 µm in first vibration mode. The difference is more obvious in the larger amplitude vibration mode (first vibration mode) of the smaller MCL. These errors may come from the difference between the actual parameters of the measurement environment and the simulated settings. As mentioned above, the Q-factor will finally start to decrease in measurement of atmospheric as the mode number of MCL increases and passes a "coincidence point". This is also confirmed on the third vibration mode MCL with size 2750 µm × 100 µm × 50 µm. Nonetheless, the relative error between experimental data and simulation data is very small. It is verified that the multi physical field coupling model can precisely simulate the performance of micro-cantilever in detecting airborne particles.

The Effects of MCLs' Shapes
As we know, the dimension of resonator will directly affect its detection performance. To explore the influence of the length width ratio of the rectangular MCL on the detection performance, it is impossible to compare the MCL with different volume or even with different area. We therefore eliminate the influence factors of volume and effective area in constant volume, and analyze the influence of the shapes of MCLs on the detection of airborne particles. Frequency search range is set to 0-600 kHz.

The Vibration Mode
MCL have different actuating, sampling, and signal read-out methods. It is necessary to refer to the vibration mode to evaluate whether it can perform normally when determining the working frequency and sampling method of the MCL-based aerosols sensors. Too high frequency and too large deflection will lead to the particle separation from the resonator. To investigate the vibration mode, a set of simulations are used to investigate the performance of the MCLs with different shapes.
As shown in Figure 6, there are three types of vibration modes of MCLs with different shapes. The first vibration mode is the standard bending mode (z-Bending mode), different waveform numbers are the characteristics that distinguish them. The second mode is torsional mode. The former two types of vibration modes can adopt the conventional actuating and sampling methods. The vibration direction of the x-Bending is along the direction of the shell surface. This is bad for sampling, since the transverse vibration will make the particles separate from the effective area easily. The actuating and signal read-out methods applicable to different vibration modes of MCLs are shown in Table 3. The actuating and signal read-out methods in torsional mode are limited. As the example of MCL proposed above, the Lorentz force, which is periodically changing in the magnetic field, is used to actuate the torsional vibration. Therefore, the bending mode of MCL is more general. Besides, the order and frequency position of each vibration mode are different between MCL with different σ lw , The σ lw affects the order and resonant frequency of various modes of MCL.

The Detection Performance
To evaluate the detection performance, the mass sensitivity of MCL structure, Q-factor and power dissipation of MCLs structure are simulated. The mass sensitivity is calculated by simulating the frequency shift of the MCL loaded with 1 ng and 2 ng particles, the Q-factor is simulated through frequency-domain sweep, and the power dissipation of one vibration period are calculated by integral power density. Figure 7a shows the linear fit of mass sensitivity and working frequency of MCLs. Compared to the continuous line, it is evident that the linear relationship between resonant frequency and mass sensitivity. In other words, the influence of of MCLs on the linear relationship between

The Detection Performance
To evaluate the detection performance, the mass sensitivity of MCL structure, Q-factor and power dissipation of MCLs structure are simulated. The mass sensitivity is calculated by simulating the frequency shift of the MCL loaded with 1 ng and 2 ng particles, the Q-factor is simulated through frequency-domain sweep, and the power dissipation of one vibration period are calculated by integral power density. Figure 7a shows the linear fit of mass sensitivity and working frequency of MCLs. Compared to the continuous line, it is evident that the linear relationship between resonant frequency and mass sensitivity. In other words, the influence of σ lw of MCLs on the linear relationship between frequency and mass sensitivity is almost negligible. However, it is safe to say that different shapes have their resonant frequencies, which determine the sensitivity. The mass sensitivity is inversely proportional to the σ lw of MCLs in each vibration mode as shown in Figure 7b. The Q-factors of MCLs in different vibration modes are depicted in Figure 8. The relationship between mass sensitivity, , and Q-factor is demonstrated. In the z-Bending vibration mode, there are three separate curves in the figure, which represent the former three vibration modes respectively. There is a noticeable downward trend with the increase of . As we can see, the YZ plane, Q-factor is linearly proportional to mass sensitivity. In other words, the Q-factor is proportional to the resonant frequency. Therefore, it is easy to determine the optimal parameters of MCL in the frequency domain. For example, the best of MCL with the same surface area is 5, and the working frequency is 540.7 kHz in the frequency domain of 0-600 kHz. In the torsional vibration mode, the range of Q-factor and sensitivity is narrow. It was clear that the mass sensitivity, and Q-factor of the MCL have poor correlation. It can be inferred that the thickness of the MCL is the parameter that determines the detection performance of the torsional vibration mode. To further evaluate the performance of this vibration mode, we can compare the power dissipation in the air. In the x-Bending vibration mode, the noticeable trend is similar to the z-Bending vibration mode. In addition, there are higher Q-factor than other modes. The line graph of the power dissipation corresponding to varying and resonant frequency is shown in Figure 9. The power dissipation is divided into thermal and viscous dissipation, which be denoted by solid line and dotted line. The viscous dissipation increased with increment in of MCLs, while decreased with increment in resonant frequency. In contrast, there is an obviously The Q-factors of MCLs in different vibration modes are depicted in Figure 8. The relationship between mass sensitivity, σ lw , and Q-factor is demonstrated. In the z-Bending vibration mode, there are three separate curves in the figure, which represent the former three vibration modes respectively. There is a noticeable downward trend with the increase of σ lw . As we can see, the YZ plane, Q-factor is linearly proportional to mass sensitivity. In other words, the Q-factor is proportional to the resonant frequency. Therefore, it is easy to determine the optimal parameters of MCL in the frequency domain. For example, the best σ lw of MCL with the same surface area is 5, and the working frequency is 540.7 kHz in the frequency domain of 0-600 kHz. In the torsional vibration mode, the range of Q-factor and sensitivity is narrow. It was clear that the mass sensitivity, σ lw and Q-factor of the MCL have poor correlation. It can be inferred that the thickness of the MCL is the parameter that determines the detection performance of the torsional vibration mode. To further evaluate the performance of this vibration mode, we can compare the power dissipation in the air. In the x-Bending vibration mode, the noticeable trend is similar to the z-Bending vibration mode. In addition, there are higher Q-factor than other modes. The Q-factors of MCLs in different vibration modes are depicted in Figure 8. The relationship between mass sensitivity, , and Q-factor is demonstrated. In the z-Bending vibration mode, there are three separate curves in the figure, which represent the former three vibration modes respectively. There is a noticeable downward trend with the increase of . As we can see, the YZ plane, Q-factor is linearly proportional to mass sensitivity. In other words, the Q-factor is proportional to the resonant frequency. Therefore, it is easy to determine the optimal parameters of MCL in the frequency domain. For example, the best of MCL with the same surface area is 5, and the working frequency is 540.7 kHz in the frequency domain of 0-600 kHz. In the torsional vibration mode, the range of Q-factor and sensitivity is narrow. It was clear that the mass sensitivity, and Q-factor of the MCL have poor correlation. It can be inferred that the thickness of the MCL is the parameter that determines the detection performance of the torsional vibration mode. To further evaluate the performance of this vibration mode, we can compare the power dissipation in the air. In the x-Bending vibration mode, the noticeable trend is similar to the z-Bending vibration mode. In addition, there are higher Q-factor than other modes. The line graph of the power dissipation corresponding to varying and resonant frequency is shown in Figure 9. The power dissipation is divided into thermal and viscous dissipation, which be denoted by solid line and dotted line. The viscous dissipation increased with increment in of MCLs, while decreased with increment in resonant frequency. In contrast, there is an obviously The line graph of the power dissipation corresponding to varying σ lw and resonant frequency is shown in Figure 9. The power dissipation is divided into thermal and viscous dissipation, which be denoted by solid line and dotted line. The viscous dissipation increased with increment in σ lw of MCLs, while decreased with increment in resonant frequency. In contrast, there is an obviously reversed trend of thermal dissipation. At first, it is crystal clear that the power dissipation primarily comes from viscous dissipation in the search frequency domain. It can be inferred that thermal dissipation will be dominant with the increase of resonant frequency. Graphs confirm the linear relation between the σ lw and logarithm of both thermal and viscous dissipation. There is a logarithmic relation between resonant frequency and both thermal and viscous dissipation. In addition, it can be concluded that the thermal dissipation has little correlation with the vibration mode, and its change mainly depends on the value of resonant frequency. However, viscous dissipation is not only related to resonant frequency but also to the vibration mode.
Sensors 2020, 20, x FOR PEER REVIEW 12 of 19 reversed trend of thermal dissipation. At first, it is crystal clear that the power dissipation primarily comes from viscous dissipation in the search frequency domain. It can be inferred that thermal dissipation will be dominant with the increase of resonant frequency. Graphs confirm the linear relation between the and logarithm of both thermal and viscous dissipation. There is a logarithmic relation between resonant frequency and both thermal and viscous dissipation. In addition, it can be concluded that the thermal dissipation has little correlation with the vibration mode, and its change mainly depends on the value of resonant frequency. However, viscous dissipation is not only related to resonant frequency but also to the vibration mode. In this section, the shape of conventional MCLs effects on the vibration mode and detection performance is simulated. The of conventional MCLs affect the frequency of various vibration modes and sensitivity. The relationship between Q-factor and is different in each vibration mode. Overall, the Q-factor is inversely proportional to the in the bending modes. The viscous dissipation is increasing and the thermal is decreasing with increment in .

Effect of MCL Slots
In this section, a set of simulations are used to investigate the slotted effects on deflection of MCLs.
The 1:1 line graph of variation in slotted and conventional MCLs is shown in Figure 10. Compared with the 1:1 line, which is represented as one continuous line. The green dotted line is a linear fit curve. We can see that only the fitting curve of sensitivity is close to the line of 1:1, and others have different degrees of offset. Among them, deflection has the largest deviation. This shows that the deflection of micro cantilevers increases significantly. Besides, there is a significant increase in the power dissipation of the slotted MCLs relative to the conventional MCLs is shown in Figure 10c. This may be due to the fact that the slotted structure increases the deflection of the MCLs. As expected, In this section, the shape of conventional MCLs effects on the vibration mode and detection performance is simulated. The σ lw of conventional MCLs affect the frequency of various vibration modes and sensitivity. The relationship between Q-factor and σ lw is different in each vibration mode. Overall, the Q-factor is inversely proportional to the σ lw in the bending modes. The viscous dissipation is increasing and the thermal is decreasing with increment in σ lw .

Effect of MCL Slots
In this section, a set of simulations are used to investigate the slotted effects on deflection of MCLs. The 1:1 line graph of variation in slotted and conventional MCLs is shown in Figure 10. Compared with the 1:1 line, which is represented as one continuous line. The green dotted line is a linear fit curve. We can see that only the fitting curve of sensitivity is close to the line of 1:1, and others have different degrees of offset. Among them, deflection has the largest deviation. This shows that the deflection of micro cantilevers increases significantly. Besides, there is a significant increase in the power dissipation of the slotted MCLs relative to the conventional MCLs is shown in Figure 10c. This may be due to the fact that the slotted structure increases the deflection of the MCLs. As expected, the maximum deflection at the piezoresistive resistance area with length =50 µm of the slotted MCLs are larger, which is clearly shown in Figure 10d. It is worth noting that the growth rate corresponding to the green arrow in the figure is the largest one. This point corresponds to a torsional mode with σ lw = 1. The line graph of variation in deflection of slotted and conventional MCLs with distance to anchor end is shown in Figure 11a. Figure 11b shows 3D visualization of slotted and conventional MCLs deflection in torsional mode. In the torsional mode with σ lw = 1, the wider width and the longer anchor end result in an enormous deflection change.
Sensors 2020, 20, x FOR PEER REVIEW 13 of 19 the maximum deflection at the piezoresistive resistance area with length =50 μm of the slotted MCLs are larger, which is clearly shown in Figure 10d. It is worth noting that the growth rate corresponding to the green arrow in the figure is the largest one. This point corresponds to a torsional mode with = 1. The line graph of variation in deflection of slotted and conventional MCLs with distance to anchor end is shown in Figure 11a. Figure 11b shows 3D visualization of slotted and conventional MCLs deflection in torsional mode. In the torsional mode with = 1, the wider width and the longer anchor end result in an enormous deflection change. As a conclusion, there is an evident effect of slots on increasing the deflection of MCLs. This will increase the sensitivity of piezoresistive materials. However, it brings the disadvantages of increasing the power dissipation. In general, it is acceptable to enhance the readout weak signals of MCL-based mass sensors at the expense of increased power dissipation in the low frequency range. As a conclusion, there is an evident effect of slots on increasing the deflection of MCLs. This will increase the sensitivity of piezoresistive materials. However, it brings the disadvantages of increasing the power dissipation. In general, it is acceptable to enhance the readout weak signals of MCL-based mass sensors at the expense of increased power dissipation in the low frequency range.

Effect of Edge Serration
As we know, viscous and thermal damping are important factors that cannot be ignored when design the MCL-based aerosols mass sensors. We can measure damping amount by the power dissipation caused by damping which is divided into thermal and viscous dissipation. In this section, a set of simulations are used to verify the effect of the micro edge serrations on damping based on the slotted MCL described above ( = 6, in 1st Bending vibration mode and torsional vibration mode). The temperature variation of slotted MCL in 1st Bending mode and Torsional mode are shown in Figure 12a,c, respectively. The position of maximum temperature variation is opposite to that of minimum temperature variation. As we know, temperature variation is closely related to pressure variation through equation of state. The temperature gradient near the surface of the MCL under isothermal conditions leads to thermal dissipation. As an obvious contrast, it can be seen that the MCL with edge serrations has a smaller temperature variation as shown in Figure 12b,d. The max temperature variation in slotted MCL and slotted MCL with edge serrations are summarized in Table 4. Lower thermal dissipation is observed in slotted MCL with edge serrations.

Effect of Edge Serration
As we know, viscous and thermal damping are important factors that cannot be ignored when design the MCL-based aerosols mass sensors. We can measure damping amount by the power dissipation caused by damping which is divided into thermal and viscous dissipation. In this section, a set of simulations are used to verify the effect of the micro edge serrations on damping based on the slotted MCL described above (σ lw = 6, in 1st Bending vibration mode and torsional vibration mode).
The temperature variation of slotted MCL in 1st Bending mode and Torsional mode are shown in Figure 12a,c, respectively. The position of maximum temperature variation is opposite to that of minimum temperature variation. As we know, temperature variation is closely related to pressure variation through equation of state. The temperature gradient near the surface of the MCL under isothermal conditions leads to thermal dissipation. As an obvious contrast, it can be seen that the MCL with edge serrations has a smaller temperature variation as shown in Figure 12b,d. The max temperature variation in slotted MCL and slotted MCL with edge serrations are summarized in Table 4. Lower thermal dissipation is observed in slotted MCL with edge serrations. There is an obvious difference between slotted MCL and slotted MCL with edge serration in terms of viscous dissipation density, as shown in Figure 13. The maximum viscous dissipation density is concentrated at the edge of the slotted MCL surface. However, it is evenly distributed around the edge serrations of the serrated MCL. The edge serration-based mechanism is in terms of passive control of the airflow, the serrations play a crucial role in reducing flow fluctuations and making changes in the lamina-turbulent transition as well as the turbulence distribution around the edge of MCL. This may because the edge serration on the MCL are equivalent to the air vortex generator that slow down the transition from laminar to turbulent flow. mode).
The temperature variation of slotted MCL in 1st Bending mode and Torsional mode are shown in Figure 12a,c, respectively. The position of maximum temperature variation is opposite to that of minimum temperature variation. As we know, temperature variation is closely related to pressure variation through equation of state. The temperature gradient near the surface of the MCL under isothermal conditions leads to thermal dissipation. As an obvious contrast, it can be seen that the MCL with edge serrations has a smaller temperature variation as shown in Figure 12b,d. The max temperature variation in slotted MCL and slotted MCL with edge serrations are summarized in Table 4. Lower thermal dissipation is observed in slotted MCL with edge serrations.  There is an obvious difference between slotted MCL and slotted MCL with edge serration in terms of viscous dissipation density, as shown in Figure 13. The maximum viscous dissipation density is concentrated at the edge of the slotted MCL surface. However, it is evenly distributed around the edge serrations of the serrated MCL. The edge serration-based mechanism is in terms of passive control of the airflow, the serrations play a crucial role in reducing flow fluctuations and making changes in the lamina-turbulent transition as well as the turbulence distribution around the edge of MCL. This may because the edge serration on the MCL are equivalent to the air vortex generator that slow down the transition from laminar to turbulent flow.  Table 5 presents the variation in dissipation between slotted MCL and slotted MCL with edge serrations. As concluded above that the edge serration reduces viscous and thermal dissipation, it is found that the viscous damping plays the leading role for the performance of MCLs in constant dimensions working in the frequency domain of 0-600 kHz. In the manufacture of the MCL-based aerosol mass sensors, we should evaluate the influence of different dimensions of edge serration.
The key criteria in this part of analysis is not to substantially raise the manufacturing cost. density is concentrated at the edge of the slotted MCL surface. However, it is evenly distributed around the edge serrations of the serrated MCL. The edge serration-based mechanism is in terms of passive control of the airflow, the serrations play a crucial role in reducing flow fluctuations and making changes in the lamina-turbulent transition as well as the turbulence distribution around the edge of MCL. This may because the edge serration on the MCL are equivalent to the air vortex generator that slow down the transition from laminar to turbulent flow.  Table 5 presents the variation in dissipation between slotted MCL and slotted MCL with edge serrations. As concluded above that the edge serration reduces viscous and thermal dissipation, it is found that the viscous damping plays the leading role for the performance of MCLs in constant dimensions working in the frequency domain of 0-600 kHz. In the manufacture of the MCL-based aerosol mass sensors, we should evaluate the influence of different dimensions of edge serration. The key criteria in this part of analysis is not to substantially raise the manufacturing cost.

Conclusions
A finite element fluid-structure interaction model was used to investigate the designed and optimized the structure of MCL-based aerosol mass sensors. The solid MCLs are modeled using shell elements and the surrounding air is modeled using Thermoviscous Acoustics and Pressure Acoustic.

Conclusions
A finite element fluid-structure interaction model was used to investigate the designed and optimized the structure of MCL-based aerosol mass sensors. The solid MCLs are modeled using shell elements and the surrounding air is modeled using Thermoviscous Acoustics and Pressure Acoustic. The effects of σ lw on the rectangular MCL-based aerosols mass sensor with a constant thickness is investigated in detail. The vibration modes and performance of MCLs in different shapes (σ lw = 1, 2, 3, 4, 5, 6,7,8,9) were modeled. It is found that the σ lw of conventional MCLs affect the frequency of various vibration modes and sensitivity; the Q-factor is inversely proportional to σ lw in the bending modes; the viscous dissipation is increased and the thermal is decreased with increment in σ lw . Different frequency domain and volume of MCL were varied to optimize the MCL performance. Therefore, the optimal σ lw (5) and resonant frequency (540.7 kHz) of the MCL with the constant surface area (45,000 µm 2 ) in the frequency domain from 0 kHz to 600 kHz were obtained, the Q factor is 2860.43, the mass sensitivity is 511.78 Hz/ng. Combined with a slotted structure, the read-out signal can be enhanced. We further found that the MCL slots have obvious effect on increasing the deflection of MCLs, but it introduces the disadvantage of increasing the power dissipation. Therefore, as far as we know, it is the first time researchers added edge serration to the edge of MCL and carried out a qualitative simulation. As expected, the edge serration on the MCL-based aerosols sensors can reduce the damping in fluid. It will greatly reduce the damping limitation of the MCL in fluid detection under high frequency mode.