Tổng quan luận án

Luận án tiến sĩ "On the quality factors of MEMS Resonators" (Nghiên cứu hệ số phẩm chất của các bộ cộng hưởng MEMS) được thực hiện bởi nghiên cứu sinh Chi-Cuong Nguyen (Nguyễn Chí Cường) dưới sự hướng dẫn của Giáo sư Wang-Long Li (Lý Vượng Long) tại Khoa Khoa học và Kỹ thuật Vật liệu, Trường Đại học Quốc gia Thành Công (National Cheng Kung University - NCKU), Đài Loan, hoàn thành vào tháng 5 năm 2017.

Tính cấp thiết và khoảng trống nghiên cứu

Các bộ cộng hưởng vi cơ điện tử (MEMS resonators) đóng vai trò trung tâm trong cấu tạo của nhiều thiết bị đo lường và cảm biến như gia tốc kế, con quay vi cơ (micro gyroscopes), gương vi cơ xoắn (torsional micromirrors), micro thu âm, cũng như các đầu dò kính hiển vi lực nguyên tử (AFM tips) và thiết bị thu hoạch năng lượng rung động. Hai đặc tính động lực học cốt lõi của bộ cộng hưởng MEMS là tần số cộng hưởng ($\omega_n$) và hệ số phẩm chất ($Q$-factor). Hệ số phẩm chất cao tương ứng với mức độ tiêu tán năng lượng thấp, là yêu cầu then chốt để nâng cao độ nhạy, độ phân giải và độ ổn định tổng thể của các hệ thống cảm biến MEMS.

Khi các cấu trúc vi cơ bị thu nhỏ kích thước xuống thang micromet, tỷ số giữa diện tích bề mặt và thể tích tăng lên đáng kể. Trong môi trường chất khí xung quanh, hiện tượng cản màng nén (squeeze film damping - SFD) xuất hiện do lớp khí bị ép trong khe hở hẹp giữa cấu trúc dao động và đế cố định, trở thành nguồn tiêu tán năng lượng chủ đạo làm suy giảm nghiêm trọng hệ số phẩm chất ở mode dao động cơ bản. Để giảm thiểu tổn hao do SFD, các cấu trúc vi dầm thường được vận hành ở áp suất thấp hoặc trong khe hở siêu mỏng. Khi đó, hiệu ứng loãng khí (gas rarefaction) và độ nhám bề mặt (surface roughness) phát sinh và tương tác phức tạp với nhau. Tuy nhiên, các nghiên cứu trước đây chưa giải quyết đồng thời và toàn diện tác động kết hợp của độ loãng khí (bao gồm số Knudsen nghịch đảo $D$ và hệ số thích ứng bề mặt ACs), độ nhám bề mặt (tỷ số bề dày màng $H_s$ và số Peklenik $\gamma$) cùng với các mode dao động uốn bậc cao lên tổng hệ số phẩm chất $Q_T$, vốn bao gồm cả tổn hao màng nén (SFD), tổn hao nhiệt đàn hồi nội tại (thermoelastic damping - TED) và tổn hao qua mấu ngàm (anchor loss).

Mục tiêu nghiên cứu

  1. Thiết lập mô hình giải tích và số trị toàn diện giải đồng thời phương trình bôi trơn khí phân tử hiệu chỉnh (MMGL) và phương trình dao động ngang của vi dầm thông qua bài toán giá trị riêng phức (complex eigenvalue problem) bằng phương pháp phần tử hữu hạn (FEM) để tính toán tần số cộng hưởng và hệ số phẩm chất $Q_{SFD}$.
  2. Đánh giá định lượng cơ chế tiêu tán năng lượng nội tại do cản nhiệt đàn hồi ($Q_{TED}$) thông qua việc giải hệ phương trình nhiệt kết hợp và dao động ngang bằng FEM, đồng thời xác định tổn hao màng ngàm ($Q_{anchor}$) theo mô hình sóng đàn hồi giải tích.
  3. Phân tích sự phụ thuộc của tổng hệ số phẩm chất ($Q_T$) và tỷ trọng đóng góp của SFD ($W_{SFD}$) dưới tác động đơn lẻ và kết hợp của độ loãng khí ($0{,}01 \le D \le 100$; $0{,}1 \le \alpha_1, \alpha_2 \le 1{,}0$), độ nhám bề mặt ($H_s = 3, 4, 6, 8, 10$; $\gamma = 9, 1, 1/9$) trên các mode dao động uốn khác nhau (từ mode 1 đến các mode bậc cao hơn 3).

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Các cấu trúc vi dầm cộng hưởng MEMS dạng tiết diện chữ nhật (micro-cantilever beams, clamped-clamped micro-beams) dao động uốn ngang trong môi trường chất khí vi mô.
  • Phạm vi không gian và thông số hình học: Khe hở màng khí $h$ ở mức micromet ($h_0 = 4, 6, 8, 10\ \mu\text{m}$), chiều dài dầm $L$, bề rộng $W$, bề dày $H$ ở thang vi mô; dải độ nhám bề mặt với số Peklenik $\gamma \in [1/9, 9]$ (đại diện cho bề mặt vân ngang, đẳng hướng và vân dọc).
  • Phạm vi vật lý và môi trường: Dải áp suất từ chân không, áp suất thấp đến áp suất khí quyển ($p \in [0{,}01\ \text{atm}, 1\ \text{atm}]$ hoặc biểu diễn qua số Knudsen nghịch đảo $D \in [0{,}01, 100]$); hệ số thích ứng động lượng tiếp tuyến đối xứng ($\alpha_1 = \alpha_2$) và phi đối xứng ($\alpha_1 \neq \alpha_2$) trong khoảng $[0{,}1; 1{,}0]$; nhiệt độ môi trường chuẩn $T_0 = 300\ \text{K}$.

Tổng quan tài liệu và vị trí của luận án

Luận án tiến hành điểm luận và phân loại có hệ thống các công trình nghiên cứu về các cơ chế tiêu tán năng lượng trong cấu trúc vi cơ:

                            CƠ CHẾ TIÊU TÁN NĂNG LƯỢNG MEMS
Tổn hao cản màng nén (SFD)      Tổn hao nhiệt đàn hồi (TED)        Tổn hao qua mấu ngàm (Anchor)
- Newell (1968)                 - Zener (1937, 1938)               - Jimbo & Itao (1968)
- Christian (1966)              - Lifshitz & Roukes (2000)         - Cross & Lifshitz (2001)
- Bao et al. (2002)             - Srikar & Senturia (2002)         - Hao et al. (2003)
- Blech (1983)                  - Nayfeh & Younis (2004)           - Park & Park (2004)
- Nayfeh & Younis (2004)        - De & Aluru (2006)                - Photiadis et al. (2004)
- Pandey & Pratap (2007)        - Rezazadeh et al. (2009)          - Chouvion (2010)
- Li (1999, 2003, 2008)         - Li et al. (2012)

Các hướng nghiên cứu trước đây

  1. Nghiên cứu về cản màng nén (SFD) và hiệu ứng trượt/loãng khí:

    • Newell (1968) phân chia ảnh hưởng của áp suất khí lên dao động cơ học thành 3 vùng: chân không tự do (áp suất rất thấp, cản khí không đáng kể), vùng chuyển tiếp phân tử và vùng nhớt liên tục (gần áp suất khí quyển).
    • Christian (1966) và Bao et al. (2002) xây dựng các mô hình giải tích dựa trên sự truyền động lượng phân tử trong chế độ dòng phân tử tự do, tuy nhiên giả định dịch chuyển màng khí đồng nhất khiến sai số lớn khi áp dụng cho kết cấu dầm uốn linh hoạt.
    • Burgdorfer (1958), Hsia và Domoto (1983), Mitsuya (1993) lần lượt đề xuất các điều kiện biên trượt bậc 1, bậc 2 và bậc 1.5 để hiệu chỉnh phương trình Reynolds cho dòng khí loãng.
    • Fukui và Kaneko (1988, 1990) giải phương trình Boltzmann tuyến tính hóa trên mô hình Bhatnagar-Gross-Krook (BGK) để dẫn xuất phương trình bôi trơn khí phân tử (MGL) và lập cơ sở dữ liệu hệ số lưu lượng dòng Poiseuille $Q_p(D)$.
    • Veijola et al. (1995, 1997) và Li & Hughes (2000) đưa ra các công thức xấp xỉ cho độ nhớt hiệu dụng và hệ số lưu lượng Poiseuille có tính đến hệ số thích ứng bề mặt (ACs).
    • Li (1999, 2002, 2003, 2004, 2008) phát triển cơ sở dữ liệu hoàn chỉnh cho hệ số hiệu chỉnh dòng Poiseuille $Q_p(D, \alpha_1, \alpha_2)$ và dòng Couette $Q_c(D, \alpha_1, \alpha_2)$ cho dải $0{,}01 \le D \le 100$ và $0{,}1 \le \alpha_1, \alpha_2 \le 1{,}0$, ứng dụng cho đầu đọc ổ đĩa cứng và vi gương MEMS, nhưng chưa tính toán hệ số phẩm chất cho vi dầm uốn ở nhiều mode.
    • Pandey và Pratap (2007) đề xuất mô hình giải tích tính toán hệ số cản SFD và độ cứng lò xo có xét đến ảnh hưởng của các mode uốn của dầm công xôn vi cơ, sử dụng mô hình hiệu chỉnh của Veijola.

  2. Nghiên cứu về cản nhiệt đàn hồi (TED):

    • Zener (1937, 1938) đặt nền tảng lý thuyết cho TED trong vật liệu đẳng hướng dựa trên mô hình phi đàn hồi chuẩn và phương pháp tiếp cận nhiệt - năng lượng.
    • Lifshitz và Roukes (2000) dẫn xuất công thức giải tích chính xác hơn cho các vi dầm hình chữ nhật chịu dao động uốn điều hòa.
    • Nayfeh và Younis (2004), De và Aluru (2006), Sun et al. (2009), Rezazadeh et al. (2009), Li et al. (2012) mở rộng phân tích TED cho vi tấm và dầm chịu tải tĩnh điện bằng mô hình truyền nhiệt 1D và 2D hyperbolic kết hợp mô phỏng phần tử hữu hạn.

  3. Nghiên cứu về tổn hao qua mấu ngàm (Anchor loss):

    • Jimbo và Itao (1968), Cross và Lifshitz (2001) xây dựng biểu thức truyền sóng đàn hồi 2D từ dầm sang đế vô hạn ở mode cơ bản.
    • Hao et al. (2003) thiết lập mô hình giải tích dựa trên nền tảng của Miller và Pursey (1954), cung cấp công thức tính $Q_{anchor}$ phụ thuộc tường minh vào tần số cộng hưởng và mode dao động cho cả cấu trúc dầm ngàm-tự do và ngàm-ngàm.
    • Park và Park (2004a, b), Chouvion (2010), Photiadis et al. (2004), Judge et al. (2007) tiếp tục phát triển các mô hình truyền sóng và lý thuyết đế có bề dày hữu hạn.

  4. Nghiên cứu về độ nhám bề mặt:

    • Patir và Cheng (1978), Bhushan và Tonder (1989a, b) đề xuất phương trình Reynolds trung bình thông qua các hệ số dòng áp suất và dòng trượt.
    • Li et al. (1995), Li và Weng (1997), Li (1999) mở rộng phương trình bôi trơn khí phân tử hiệu chỉnh (MMGL) kết hợp các hệ số dòng áp suất $\phi_x, \phi_y$ cho bề mặt nhám dưới điều kiện dòng khí trượt và khí loãng.

Khoảng trống nghiên cứu được luận án giải quyết

Mặc dù các cơ chế SFD, TED và Anchor loss đã được nghiên cứu riêng lẻ, các mô hình trước đây thường bỏ qua hệ số thích ứng bề mặt phi đối xứng ($\alpha_1 \neq \alpha_2$), chưa tích hợp đồng thời cơ sở dữ liệu lưu lượng khí phân tử chính xác với độ nhám bề mặt ngẫu nhiên có hướng (số Peklenik $\gamma$), và chưa đánh giá sự chuyển dịch vai trò chi phối (weighting of SFD) trên dải mode dao động rộng (từ mode 1 đến các mode bậc cao) trong bài toán giá trị riêng phức tổng thể. Luận án này giải quyết trọn vẹn khoảng trống đó.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án sử dụng phương pháp mô hình hóa toán lý định lượng kết hợp phân tích số trị phần tử hữu hạn (FEM).

1. Hệ phương trình chi phối dòng màng khí nén (SFD)

Phương trình chuyển động của chất lưu nhớt Navier-Stokes và phương trình liên tục được viết trong hệ tọa độ Descartes $(x, y, z)$:

$$\rho \left( \frac{\partial u}{\partial t} + u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} + w\frac{\partial u}{\partial z} \right) = -\frac{\partial p}{\partial x} + \frac{\partial \sigma_x}{\partial x} + \frac{\partial \tau_{yx}}{\partial y} + \frac{\partial \tau_{zx}}{\partial z} + \rho X_b$$

$$\frac{\partial \rho}{\partial t} + \frac{\partial (\rho u)}{\partial x} + \frac{\partial (\rho v)}{\partial y} + \frac{\partial (\rho w)}{\partial z} = 0$$

Trong đó, ứng suất pháp $\sigma_i$ và ứng suất cắt $\tau_{ij}$ được xác định qua độ nhớt động lực học $\mu$ và hệ số độ nhớt thứ hai $\lambda_B = -2\mu/3$.

Khi áp dụng lý thuyết bôi trơn cho màng khí nén mỏng giữa hai bề mặt nhám có tính đến hiệu ứng loãng khí, phương trình bôi trơn khí phân tử hiệu chỉnh (MMGL) được dẫn xuất dưới dạng:

$$\frac{\partial}{\partial x} \left( \phi_x Q_p(D, \alpha_1, \alpha_2) p h^3 \frac{\partial p}{\partial x} \right) + \frac{\partial}{\partial y} \left( \phi_y Q_p(D, \alpha_1, \alpha_2) p h^3 \frac{\partial p}{\partial y} \right) = 12 \mu \frac{\partial (p h)}{\partial t}$$

Trong đó:

  • $D = \frac{\sqrt{\pi}}{2 K_n} = \frac{\sqrt{\pi} p_0 h}{2 \lambda_0 p}$ là số Knudsen nghịch đảo, đặc trưng cho mức độ loãng khí.
  • $Q_p(D, \alpha_1, \alpha_2)$ là hệ số hiệu chỉnh lưu lượng Poiseuille được nội suy từ cơ sở dữ liệu số giải phương trình Boltzmann tuyến tính hóa (Li, 2003).
  • $\phi_x, \phi_y$ là các hệ số dòng áp suất mô tả hướng và độ cao nhám bề mặt, phụ thuộc vào tỷ số bề dày màng $H_s = h/\sigma$ và số Peklenik $\gamma = \gamma_x / \gamma_y$.
Thông số Ký hiệu Ý nghĩa vật lý / Giá trị khảo sát
Số Knudsen nghịch đảo $D$ Mức độ loãng khí ($0{,}01 \le D \le 100$)
Hệ số thích ứng bề mặt $\alpha_1, \alpha_2$ Tỷ lệ va chạm khuếch tán trên bề mặt động và tĩnh ($0{,}1 \le \alpha \le 1{,}0$)
Tỷ số bề dày màng $H_s = h/\sigma$ Mức độ gồ ghề của bề mặt so với khe hở ($H_s = 3, 4, 6, 8, 10$)
Số Peklenik $\gamma = \gamma_x/\gamma_y$ Định hướng vân nhám: $\gamma = 9$ (vân dọc), $\gamma = 1$ (đẳng hướng), $\gamma = 1/9$ (vân ngang)

2. Phương trình dao động ngang của vi dầm

Dao động uốn ngang $w(x,t)$ của vi dầm đàn hồi dưới tác dụng của lực khí động phân bố $\Delta p(x, y, t)$ được mô tả bởi phương trình:

$$E I \frac{\partial^4 w}{\partial x^4} + \rho_b A \frac{\partial^2 w}{\partial t^2} = \int_{-W/2}^{W/2} \Delta p(x, y, t) , dy$$

Với $E$ là mô đun Young, $I$ là mô men quán tính mặt cắt ngang, $\rho_b$ là khối lượng riêng vật liệu dầm, $A = W \cdot H$ là diện tích mặt cắt.

3. Phương trình cản nhiệt đàn hồi (TED) và mô hình tổn hao ngàm (Anchor loss)

  • TED: Biến dạng uốn tuần hoàn sinh ra gradien nhiệt độ qua chiều dày dầm, dẫn đến dòng nhiệt truyền không thuận nghịch. Phương trình truyền nhiệt ghép với chuyển vị ngang được giải số theo phương pháp Zener và Lifshitz-Roukes.
  • Anchor loss: Được tính theo công thức giải tích của Hao et al. (2003) cho từng mode dao động $n$:

$$Q_{anchor}^{(n)} = C_n \left( \frac{L}{H} \right)^3 f(\nu, n)$$

Trong đó $C_n$ là hệ số phụ thuộc điều kiện biên ngàm và mode dao động, $\nu$ là hệ số Poisson.

4. Quy trình giải bài toán giá trị riêng phức bằng FEM

Tuyến tính hóa các biến động điều hòa: $w(x, t) = \bar{w}(x) e^{\lambda t}$, $p(x, y, t) = p_0 + \bar{p}(x, y) e^{\lambda t}$, với $\lambda = -\delta + i \omega_n$ là giá trị riêng phức. Rời rạc hóa hệ phương trình vi phân liên kết bằng FEM tạo thành hệ phương trình đại số ma trận giá trị riêng:

$$(\mathbf{K} + \lambda \mathbf{C} + \lambda^2 \mathbf{M}) \mathbf{U} = \mathbf{0}$$

Từ nghiệm $\lambda$, các thông số động lực học được trích xuất trực tiếp:

  • Tần số cộng hưởng: $\omega_n = \text{Im}(\lambda)$ (rad/s) hay $f_n = \omega_n / (2\pi)$ (Hz).
  • Hệ số cản: $\delta = -\text{Re}(\lambda)$.
  • Hệ số phẩm chất của từng cơ chế: $Q = \frac{\omega_n}{2 \delta} = \frac{\text{Im}(\lambda)}{2 |\text{Re}(\lambda)|}$.
  • Tổng hệ số phẩm chất:

$$\frac{1}{Q_T} = \frac{1}{Q_{SFD}} + \frac{1}{Q_{TED}} + \frac{1}{Q_{anchor}} = \frac{1}{Q_{SFD}} + \frac{1}{Q_{TA}}$$

  • Tỷ trọng đóng góp của cản màng nén: $W_{SFD} = \frac{Q_T}{Q_{SFD}} \times 100%$.

Nội dung chính theo từng chương

                                  CẤU TRÚC LUẬN ÁN
Chương 1    Chương 2    Chương 3                        Chương 4    Chương 5    Chương 6
Tổng quan   Mô hình     Phương pháp phần tử hữu hạn     Hiệu ứng    Ảnh hưởng   Kết luận &
& Bối cảnh  lý thuyết   & Bài toán giá trị riêng phức   loãng khí   kết hợp     Định hướng

Chương 1: Giới thiệu (Introduction)

Chương này tổng thuật chi tiết cơ sở lý thuyết dao động của hệ một bậc tự do (SDOF) tương đương (mô hình khối lượng - lò xo - giảm chấn với các tham số tập trung $m, k, c$), định nghĩa tỷ số cản $\zeta = c / (2\sqrt{km}) = 1/(2Q)$ và phương pháp đo $Q$-factor thực nghiệm bằng băng thông nửa công suất ($3\text{dB}$ bandwidth: $Q = \omega_r / \Delta\omega$).

Chương 1 phân tích cặn kẽ 4 cơ chế tiêu tán năng lượng chính trong cấu trúc vi dầm MEMS:

  1. Cản màng nén (SFD) trong môi trường khí quyển và áp suất thấp.
  2. Cản nhiệt đàn hồi nội tại (TED) do gradient nhiệt sinh ra trong quá trình co giãn tuần hoàn.
  3. Tổn hao mấu ngàm (Anchor loss) do sóng đàn hồi truyền qua các điểm liên kết vào đế cố định.
  4. Tổn hao bề mặt (Surface loss) do sai hỏng mạng, tạp chất và lớp hấp phụ bề mặt (có thể khử bằng xử lý nhiệt luyện/annealing).

Chương 1 xác định rõ mục tiêu và phạm vi nghiên cứu của luận án, tập trung vào việc mô hình hóa SFD bằng phương trình MMGL kết hợp cơ sở dữ liệu dòng khí phân tử và các hệ số dòng áp suất nhám bề mặt.

Chương 2: Các mô hình lý thuyết (Theoretical Models)

Chương 2 trình bày chi tiết các bước thiết lập toán học:

  • Xuất phát từ hệ phương trình Navier-Stokes 3 chiều cho chất lưu nhớt nén được và phương trình liên tục bảo toàn khối lượng, áp dụng các giả thiết bôi trơn để suy dẫn phương trình Reynolds nén được.
  • Tích hợp hệ số hiệu chỉnh lưu lượng Poiseuille $Q_p(D, \alpha_1, \alpha_2)$ và hệ số dòng áp suất $\phi_x, \phi_y$ để thu được phương trình bôi trơn khí phân tử hiệu chỉnh (MMGL).
  • Thiết lập phương trình dao động ngang của vi dầm chữ nhật chịu kích thích điều hòa dưới tác dụng của lực màng khí nén.
  • Thiết lập phương trình truyền nhiệt vi phân ghép với trường biến dạng đàn hồi để xác định tổn hao nhiệt đàn hồi TED.
  • Giới thiệu mô hình giải tích dạng tường minh tính toán tổn hao qua mấu ngàm $Q_{anchor}$ phụ thuộc vào từng mode dao động uốn của vi dầm theo công trình của Hao et al. (2003).

Chương 3: Quy trình bài toán giá trị riêng và tính toán hệ số phẩm chất (Eigenvalue Problems Procedures and Quality Factors Calculations)

Chương 3 tập trung vào phương pháp luận số trị:

  • Dẫn xuất dạng biến phân yếu (weak form) của hệ phương trình vi phân liên kết giữa chuyển vị uốn của dầm và biến thiên áp suất màng khí.
  • Áp dụng phương pháp phần tử hữu hạn (FEM) để rời rạc hóa không gian cho vi dầm và màng khí, thiết lập bài toán giá trị riêng phức tổng quát $(\mathbf{K} + \lambda \mathbf{C} + \lambda^2 \mathbf{M}) \mathbf{U} = \mathbf{0}$.
  • Xây dựng thuật toán tính toán giá trị riêng phức $\lambda_n$, từ đó trích xuất tần số cộng hưởng $\omega_n = \text{Im}(\lambda_n)$ và hệ số suy giảm $\delta_n = -\text{Re}(\lambda_n)$.
  • Xây dựng quy trình tổng hợp hệ số phẩm chất tổng $Q_T$ và xác định tỷ trọng cản màng nén $W_{SFD}$.

Chương 4: Ảnh hưởng của hiệu ứng loãng khí đến hệ số phẩm chất của bộ cộng hưởng vi dầm (Effects of Gas Rarefaction on the Quality Factors of Micro-beam Resonators)

Chương này công bố các kết quả tính toán chi tiết khi xét bề mặt nhẵn trên dải số Knudsen nghịch đảo $D_0 \in [0{,}01; 100]$ và hệ số thích ứng $\alpha_1, \alpha_2 \in [0{,}1; 1{,}0]$:

  • Kiểm tra sự hội tụ của nghiệm phần tử hữu hạn theo số phần tử lưới $N' = (N_h \times N_w \times N_l)$, khẳng định tính ổn định của thuật toán.
  • So sánh kiểm chứng tỷ số cản $\zeta = 1/(2Q)$ với kết quả mô phỏng và thực nghiệm của Pandey và Pratap (2007) cho thấy sự tương thích chính xác với sai số phần trăm rất nhỏ.
  • Phân tích chi tiết hệ số phẩm chất $Q_{SFD}$, $Q_{TED}$, $Q_{anchor}$ và tổng hệ số phẩm chất $Q_T$ trên 3 mode dao động đầu tiên (mode 1, mode 2, mode 3) và các mode bậc cao hơn 3.
  • Đánh giá sự phụ thuộc của tỷ trọng $W_{SFD}$ vào áp suất môi trường ($p$), hệ số thích ứng đối xứng ($\alpha_1 = \alpha_2$) và phi đối xứng ($\alpha_1 \neq \alpha_2$).
Mode dao động ($n$) $Q_{SFD}$ (ở áp suất khí quyển $D_0=10$) $Q_{TED}$ $Q_{anchor}$ Cơ chế tiêu tán chủ đạo
Mode 1 Thấp (bị chi phối mạnh bởi SFD) Rất cao Rất cao SFD chiếm ưu thế ($W_{SFD} \approx 100%$)
Mode 2 Tăng đáng kể so với Mode 1 Cao Cao SFD giảm dần tỷ trọng
Mode 3 Tăng cao Giảm dần Giảm dần Bắt đầu có sự cạnh tranh giữa SFD và TED/Anchor
Mode > 3 Rất cao (SFD suy giảm mạnh) Thấp (chi phối) Thấp (chi phối) TED và Anchor loss chiếm ưu thế tuyệt đối

Chương 5: Ảnh hưởng kết hợp của độ nhám bề mặt và hiệu ứng loãng khí đến hệ số phẩm chất của bộ cộng hưởng vi dầm (Coupled Effects of Surface Roughness and Gas Rarefaction on the Quality Factors of Micro-beam Resonators)

Chương 5 nghiên cứu toàn diện sự tương tác đa trường giữa độ nhám bề mặt và hiệu ứng loãng khí:

  • Khảo sát các hệ số dòng áp suất $\phi_x, \phi_y$ theo số Peklenik $\gamma \in {1/9, 1, 9}$ và tỷ số bề dày màng $H_s \in {3, 4, 6}$ ở các mức loãng khí $D_0$ và hệ số ACs khác nhau.
  • Phân tích ảnh hưởng của bề dày khe hở ban đầu $h_0 = 4, 6, 8, 10\ \mu\text{m}$ đến phần thực $\delta$ và phần ảo $\omega_n$ của giá trị riêng phức.
  • Đánh giá tỷ số $Q_{sro} / Q_{ssm}$ (hệ số phẩm chất của bề mặt nhám so với bề mặt nhẵn) theo tỷ số bề dày màng $H_s$ và số Peklenik $\gamma$ tại các điều kiện $D_0 = 0{,}1; 1; 10$.
  • Chứng minh hiện tượng hiệu ứng loãng khí làm suy giảm (pha loãng - diluted) ảnh hưởng của độ nhám bề mặt. Tuy nhiên, khi vận hành ở các mode bậc cao và độ loãng khí cao, độ nhám bề mặt quay trở lại ảnh hưởng đáng kể đến tổng hệ số phẩm chất.
  • Khảo sát sự phân bố tỷ trọng $W_{SFD}(%)$ theo áp suất môi trường đối với các dạng vân nhám bề mặt khác nhau ($\gamma = 9$ vân dọc, $\gamma = 1$ đẳng hướng, $\gamma = 1/9$ vân ngang).

Chương 6: Kết luận và hướng nghiên cứu tương lai (Conclusions and Future Work)

Tóm tắt các kết luận cốt lõi thu được từ toàn bộ luận án và đề xuất các phương hướng mở rộng nghiên cứu tiếp theo.

Kết quả và những đóng góp mới

1. Đóng góp mới về mặt lý luận và phương pháp học thuật

  • Thiết lập thành công khung tính toán số trị giải đồng thời phương trình bôi trơn khí phân tử hiệu chỉnh (MMGL) và phương trình dao động cơ học của vi dầm dưới dạng bài toán giá trị riêng phức (complex eigenvalue problem) bằng phương pháp phần tử hữu hạn (FEM), khắc phục hoàn toàn nhược điểm giả định dịch chuyển màng khí đồng nhất của các mô hình giải tích cổ điển.
  • Tích hợp trọn vẹn cơ sở dữ liệu chính xác về hệ số lưu lượng Poiseuille $Q_p(D, \alpha_1, \alpha_2)$ cho dải số Knudsen nghịch đảo rộng ($0{,}01 \le D \le 100$) và hệ số thích ứng bề mặt tùy ý ($0{,}1 \le \alpha_1, \alpha_2 \le 1{,}0$), bao gồm cả trường hợp đối xứng và phi đối xứng.
  • Lần đầu tiên đánh giá đồng thời tác động kết hợp của độ nhám bề mặt định hướng (thông qua hệ số dòng áp suất $\phi_x, \phi_y$, tỷ số $H_s$ và số Peklenik $\gamma$) với hiệu ứng loãng khí lên hệ số phẩm chất của vi dầm MEMS ở nhiều mode dao động uốn khác nhau.
                                 TỔNG HỆ SỐ PHẨM CHẤT Q_T
         Áp suất cao / Mode thấp                       Áp suất thấp / Mode cao
         - SFD chiếm ưu thế tuyệt đối                  - SFD suy giảm mạnh
         - W_SFD xấp xỉ 100%                           - TED và Anchor loss chi phối chính
         - Ảnh hưởng nhám bị loãng hóa                 - Nhám bề mặt (H_s, gamma) tác động rõ rệt

2. Phát hiện và đóng góp thực tiễn

  • Quy luật chi phối của các cơ chế cản:
    • Ở các mode dao động thấp (đặc biệt là mode 1), áp suất môi trường cao và hệ số ACs lớn, cản màng nén (SFD) là nguồn tiêu tán năng lượng chủ đạo quyết định hệ số phẩm chất ($W_{SFD} \approx 100%$).
    • Khi tăng bậc mode dao động uốn (từ mode 2 trở lên) hoặc giảm áp suất môi trường (tăng độ loãng khí, giảm $D$ và giảm $\alpha$), hệ số cản SFD giảm mạnh, làm tỷ trọng $W_{SFD}$ suy giảm nhanh chóng.
    • Ở các mode bậc cao (>3) và trong môi trường áp suất thấp/chân không, cản nhiệt đàn hồi nội tại (TED) và tổn hao qua mấu ngàm (Anchor loss) trở thành các nguồn cản chủ đạo khống chế giới hạn trên của tổng hệ số phẩm chất $Q_T$.
  • Ảnh hưởng tương hỗ giữa độ nhám và độ loãng khí: Hiệu ứng loãng khí có xu hướng làm suy giảm (pha loãng) tác động của độ nhám bề mặt lên tổn hao SFD ở áp suất thấp. Tuy nhiên, đối với các bộ cộng hưởng thiết kế ở tần số cao (mode bậc cao) với khe hở siêu mỏng, thông số độ nhám bề mặt ($H_s$ và $\gamma$) vẫn gây ra sự thay đổi đáng kể đến $Q$-factor và cần được kiểm soát chặt chẽ trong chế tạo.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án

  • Luận án tập trung chủ yếu vào kết cấu vi dầm chữ nhật đồng nhất (micro-beam resonators dạng công xôn hoặc ngàm hai đầu); chưa mở rộng mô hình hóa cho các cấu trúc vi cơ có hình học phức tạp hơn như màng lục giác (cMUT), vi dầm gấp khúc dạng lược (comb-drive folded beams) hay con quay vi cơ có khối quán tính treo phức tạp.
  • Chưa trực tiếp thực hiện chuỗi thực nghiệm đo đạc độc lập để khảo sát hệ số thích ứng bề mặt $\alpha$ trong các điều kiện xử lý bề mặt cụ thể, mà chủ yếu kiểm chứng mô hình dựa trên tập dữ liệu thực nghiệm đã công bố của Pandey & Pratap (2007) và Legtenberg & Tilmans (1994).
  • Phân tích cản nhiệt đàn hồi TED chủ yếu áp dụng cho vật liệu dầm đồng chất, đẳng hướng ở nhiệt độ phòng.

Hướng nghiên cứu tương lai

  • Mở rộng mô hình bài toán giá trị riêng phức MMGL-FEM cho các cấu trúc vi cơ 2D và 3D phức tạp (màng vi cơ, vi gương xoắn đa trục, đĩa cộng hưởng MEMS).
  • Nghiên cứu thực nghiệm xác định chính xác hệ số thích ứng động lượng tiếp tuyến (TMAC) và các thông số nhám thực tế sau các quy trình chế tạo vi cơ (khắc axit, lắng đọng, ủ nhiệt).
  • Tích hợp mô hình dẫn nhiệt 2D/3D phi Fourier (non-Fourier heat conduction) để đánh giá TED trong các bộ cộng hưởng kích thước nano (NEMS) vận hành ở tần số siêu cao (VHF/UHF).

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị học thuật và ứng dụng kỹ thuật cao trong các lĩnh vực:

  1. Đối với nghiên cứu sinh và học viên sau đại học chuyên ngành Cơ điện tử vi mô (MEMS/NEMS), Khoa học Vật liệu và Cơ học chất lưu:

    • Cung cấp phương pháp luận chi tiết và chặt chẽ về việc dẫn xuất phương trình Reynolds nén được từ hệ phương trình Navier-Stokes và phương trình liên tục (Chương 2).
    • Quy trình biến đổi và giải bài toán giá trị riêng phức bằng phương pháp phần tử hữu hạn (FEM) để xác định đồng thời tần số cộng hưởng và hệ số cản (Chương 3).

  2. Đối với kỹ sư thiết kế và chuyên gia R&D vi cảm biến:

    • Cung cấp cơ sở dữ liệu và công cụ tính toán để tối ưu hóa hệ số phẩm chất $Q$ cho các thiết bị như đầu dò AFM, vi gia tốc kế, con quay vi cơ thông qua việc lựa chọn chế độ áp suất vận hành, chiều dày khe hở màng khí và mode dao động phù hợp.
    • Giúp dự báo chính xác ngưỡng chuyển tiếp mà tại đó việc hút chân không sâu không còn làm tăng thêm $Q$-factor do sự khống chế của TED và Anchor loss.

  3. Đối với các nhà nghiên cứu công nghệ bề mặt và chế tạo vi cơ:

    • Cung cấp các biểu đồ định lượng về ảnh hưởng của định hướng vân nhám (số Peklenik $\gamma$) và độ cao nhám ($H_s$) tương tác với hệ số thích ứng bề mặt (ACs), định hướng cho các yêu cầu hoàn thiện bề mặt cấu trúc sau chế tạo.

Câu hỏi thường gặp

1. Tại sao hệ số phẩm chất (Q-factor) lại là thông số quan trọng hàng đầu trong thiết kế bộ cộng hưởng MEMS?

Hệ số phẩm chất $Q$ là đại lượng không thứ nguyên biểu thị tỷ số giữa năng lượng dự trữ cực đại trên năng lượng tiêu tán trong một chu kỳ dao động ($Q = 2\pi W_{stored} / \Delta W_{cycle} = \omega_n / (2\delta) = 1/(2\zeta)$). Giá trị $Q$ càng cao đồng nghĩa với mức độ tiêu tán năng lượng càng thấp, giúp bộ cộng hưởng đạt được biên độ dao động lớn tại điểm cộng hưởng, độ chọn lọc tần số cao (băng thông $\Delta\omega$ hẹp), độ nhạy vượt trội khi đo lường các đại lượng cơ học/sinh hóa và giảm thiểu nhiễu nhiệt pha trong hệ thống vi cảm biến.

2. Phương trình MMGL khác gì so với phương trình Reynolds cổ điển trong chất lưu?

Phương trình Reynolds cổ điển dựa trên giả thuyết dòng chất lưu liên tục và không trượt tại biên solid-gas, chỉ đúng khi số Knudsen $K_n \to 0$. Phương trình bôi trơn khí phân tử hiệu chỉnh (MMGL) được bổ sung hai yếu tố quan trọng:

  • Hệ số hiệu chỉnh lưu lượng dòng Poiseuille $Q_p(D, \alpha_1, \alpha_2)$ được trích xuất từ nghiệm của phương trình động học chất khí Boltzmann tuyến tính hóa trên mô hình BGK, mô tả chính xác hiện tượng trượt biên và loãng khí từ chế độ dòng trượt, chuyển tiếp đến phân tử tự do ($0{,}01 \le D \le 100$).
  • Các hệ số dòng áp suất $\phi_x, \phi_y$ mô tả ảnh hưởng của độ nhám ngẫu nhiên có định hướng (số Peklenik $\gamma$) và tỷ số bề dày màng khí $H_s$.

3. Quy luật dịch chuyển vai trò giữa cản màng nén (SFD) và các cơ chế cản nội tại (TED, Anchor loss) diễn ra như thế nào khi thay đổi mode dao động?

Ở mode dao động cơ bản (mode 1), chuyển động uốn có dạng uốn đơn thuần tạo ra sự thay đổi thể tích màng khí nén lớn trên toàn bộ diện tích, khiến tổn hao cản màng nén SFD chiếm ưu thế tuyệt đối ($W_{SFD} \approx 100%$). Khi chuyển dịch lên các mode dao động bậc cao (mode 2, mode 3 và $>3$), biến dạng uốn cục bộ xen kẽ các nút dao động làm giảm hiệu ứng nén ép thể tích chất khí, khiến tổn hao SFD suy giảm rõ rệt. Đồng thời, tần số dao động tăng cao và biến dạng uốn uốn cong cục bộ lớn làm tăng mạnh gradien nhiệt độ qua chiều dày dầm (tăng TED) và tăng cường phát xạ sóng đàn hồi qua điểm ngàm (tăng Anchor loss), biến TED và Anchor loss thành hai cơ chế khống chế chính hệ số phẩm chất tổng $Q_T$.

4. Hệ số thích ứng bề mặt (Accommodation Coefficients - ACs) ảnh hưởng như thế nào đến lưu lượng Poiseuille và hệ số phẩm chất?

Hệ số thích ứng bề mặt $\alpha$ phản ánh mức độ trao đổi động lượng tiếp tuyến trung bình khi các phân tử khí va chạm với bề mặt rắn. Giá trị $\alpha = 1{,}0$ biểu thị phản xạ khuếch tán hoàn toàn (diffusive reflection), trong khi $\alpha = 0$ biểu thị phản xạ gương hoàn toàn (specular reflection). Khi $\alpha$ giảm (bề mặt trơn nhẵn về mặt phân tử, ít trao đổi động lượng), hiện tượng trượt biên tăng lên, hệ số lưu lượng Poiseuille $Q_p$ tăng mạnh (ví dụ: tại $D_0 = 1$, $Q_p$ biến thiên từ $105{,}23$ xuống $0{,}97$ khi thay đổi $\alpha$), làm giảm lực cản nhớt SFD và giúp nâng cao hệ số phẩm chất $Q_{SFD}$ của vi dầm.

Kết luận

Luận án tiến sĩ của Chi-Cuong Nguyen đã giải quyết trọn vẹn bài toán xác định hệ số phẩm chất của vi dầm cộng hưởng MEMS thông qua việc phát triển mô hình phần tử hữu hạn giải bài toán giá trị riêng phức kết hợp phương trình MMGL với cơ sở dữ liệu khí phân tử và các hệ số nhám bề mặt định hướng. Công trình đã làm sáng tỏ quy luật tương tác đa trường giữa độ loãng khí, thông số độ nhám và bậc mode dao động, đồng thời lượng hóa chính xác sự chuyển dịch vai trò chi phối giữa cản màng nén khí động ngoại vi (SFD) và các cơ chế tiêu tán nội tại (TED, Anchor loss). Các kết quả và phương pháp số trị được thiết lập trong luận án cung cấp cơ sở khoa học quan trọng phục vụ phân tích, thiết kế và tối ưu hóa hiệu năng động lực học cho các thế hệ linh kiện vi cơ điện tử cộng hưởng độ nhạy cao.