Tổng quan về luận án

Luận án tiến sĩ chuyên ngành Hóa lý thuyết và Hóa lý (Mã số: 9 44 01 19) với đề tài "Nghiên cứu sự phụ thuộc của quá trình phân hủy nhiệt và tốc độ nổ vào thành phần thuốc nổ hỗn hợp trên nền hexogen" do Nghiên cứu sinh Nguyễn Mậu Vương thực hiện dưới sự hướng dẫn khoa học của TS. Ngô Văn Giao và TS. Đặng Văn Đường tại Viện Khoa học và Công nghệ Quân sự (2020) là một công trình nghiên cứu tiên phong trong lĩnh vực hóa lý năng lượng cao. Trong bối cảnh khoa học vật liệu nổ quân sự và công nghiệp đòi hỏi sự cân bằng khắt khe giữa uy lực nổ cực đại và tính an toàn công nghệ cao, Hexogen (cyclotrimetylen trinitramin - RDX, $C_3H_6O_6N_6$) giữ vị trí trung tâm nhờ các đặc tính vượt trội: "Thuốc nổ RDX có tốc độ nổ cao (8380 m/s ở mật độ 1,70 g/cm3), khả năng sinh công lớn đo theo phương pháp bom chì là (450÷520) mL. Tuy nhiên, RDX có độ nhạy cao (70÷80%), không chịu nén và bị phân hủy trước khi nóng chảy".

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án xác lập là sự thiếu hụt các mô hình dữ liệu động học nhiệt và phương trình thực nghiệm có tính hệ thống mô tả mối quan hệ định lượng giữa hàm lượng cấu tử phối trộn (TNT hoặc chất thuần hóa) với các thông số phân hủy không đẳng nhiệt, nhiệt lượng nổ ($Q_v$) và tốc độ nổ ($D$). Các tài liệu quốc tế trước đây chủ yếu dừng lại ở các mác thuốc nổ cố định như Comp B (60% RDX / 40% TNT), PBX-9407 (94% RDX / 6% Exon 461) hay C4 (91% RDX), nhưng thiếu đi quy luật phổ quát khi tỷ lệ thành phần biến đổi liên tục trong dải làm việc.

Để giải quyết khoảng trống này, luận án đặt ra 3 câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  1. Câu hỏi 1 (Q1): Quy luật biến thiên của các thông số động học phân hủy nhiệt biểu kiến ($E, Z, k_T$) và thời gian bán hủy ($t_{1/2}$) chịu sự chi phối như thế nào bởi nồng độ cấu tử thứ hai (TNT từ 23% đến 60% và chất thuần hóa từ 5,0% đến 6,5%)? Giả thuyết $H_1$: Quá trình đứt mạch chuỗi phản ứng gốc tự do tăng lên khi bổ sung chất pha loãng/thuần hóa làm dịch chuyển năng lượng hoạt hóa và mở rộng dải an toàn nhiệt.
  2. Câu hỏi 2 (Q2): Mức độ tương thích pha và cơ chế tương tác nhiệt giữa RDX và các chất thuần hóa mới (hệ sáp đa cấu tử A-IX-13, A-IX-11) diễn ra theo quy luật động học bậc nhất hay phản ứng tự xúc tác? Giả thuyết $H_2$: Sự kết hợp của xerezin, axit stearic và phẩm màu sudan tạo ra lớp phủ bảo vệ vi mô làm suy giảm xác suất phân nhánh chuỗi mà không làm phá vỡ cấu trúc tinh thể của RDX.
  3. Câu hỏi 3 (Q3): Tốc độ nổ thực nghiệm và nhiệt lượng nổ thay đổi phi tuyến hay tuyến tính theo mật độ nén/đúc ($\rho$) và tỷ lệ thành phần khối lượng? Giả thuyết $H_3$: Tốc độ nổ tuân theo phương trình hàm bậc nhất mở rộng của mật độ nạp thực tế kết hợp với hiệu ứng ơtecti hóa giữa RDX và TNT.

Khung lý thuyết (Theoretical Framework) của nghiên cứu được xây dựng trên sự tích hợp chặt chẽ giữa lý thuyết phản ứng chuỗi gốc tự do phân nhánh của Semenov, nhiệt động lực học quá trình nổ hóa học và phương pháp phân tích động học không đẳng nhiệt Kissinger - Ozawa. Đóng góp đột phá của luận án là cung cấp tập dữ liệu định lượng chuẩn xác về thời gian bán hủy, khoảng nhiệt độ đúc rót an toàn ($< 100^\circ\text{C}$), cùng hệ phương trình thực nghiệm giúp tối ưu hóa thiết kế đầu đạn xuyên lõm, đạn tạo hình và đạn nổ phá sát thương trong điều kiện khí hậu nhiệt đới. Phạm vi nghiên cứu bao gồm 2 hệ thuốc nổ hỗn hợp điển hình: Hệ đúc TГ (khảo sát biến thiên TNT từ 23% đến 60%) và hệ ép A-IX-1 (gồm A-IX-11 và A-IX-13 với hàm lượng chất thuần hóa 5,0% - 6,5%), được thu thập và thực nghiệm trên các thiết bị phân tích hiện đại tiêu chuẩn quân sự.

Literature Review và Positioning

Nghiên cứu về động học phân hủy nhiệt và quá trình truyền nổ của thuốc nổ hỗn hợp trên nền RDX phát triển qua 3 dòng học thuật chính:

  1. Dòng nghiên cứu động học phân hủy nhiệt RDX đơn chất: Các nghiên cứu kinh điển của Rogers (1975) và Janney (1987) đã sử dụng DSC không đẳng nhiệt để xác định năng lượng hoạt hóa ($E$) của RDX tinh khiết dao động trong khoảng rộng từ $28,6$ đến $53,0\text{ kcal/mol}$ ($120 - 221\text{ kJ/mol}$) và thừa số trước mũ $Z = 3,98 \times 10^{10} \div 6,3 \times 10^{20}\text{ s}^{-1}$. Đến giai đoạn của Sinditskii et al. (2008), các tác giả chỉ ra rằng kích thước cỡ hạt ($50\ \mu\text{m}$ so với nano $180 - 300\text{ nm}$) cùng môi trường khí quyển (khí trơ Heli hay không khí) tạo ra sự sai biệt lớn về năng lượng hoạt hóa do hiệu ứng chuyển khối và tích tụ sản phẩm khí tự xúc tác.
  2. Dòng nghiên cứu hệ hỗn hợp RDX/TNT (Comp B, Cyclotol, TГ): Zinn & Mader (1960), tiếp nối bởi Oxley et al. (2009), tập trung làm sáng tỏ trạng thái dung dịch rắn và hỗn hợp ơtecti giữa RDX và TNT nóng chảy. Nghiên cứu của Cooper (1996) trong cẩm nang Explosives Engineering đã thiết lập phương trình kinh nghiệm cho tốc độ nổ của Comp B ($D = 7840\text{ m/s}$ tại $\rho = 1,68\text{ g/cm}^3$), khẳng định vai trò mang tính công nghệ của TNT trong đúc rót nhưng chỉ rõ nhược điểm co ngót tạo lỗ xốp khi làm nguội.
  3. Dòng nghiên cứu thuốc nổ dẻo và thuốc nổ thuần hóa (A-IX-1, PBX, C4): Xu hướng quốc tế từ Meyer et al. (2007) và các nhà khoa học Nga (Dubovik, 2005; Il'in et al., 2012) nghiên cứu việc bao bọc hạt RDX bằng màng polyme hoặc sáp hydrocacbon (parafin, xerezin, axit stearic) để giảm độ nhạy va đập từ $70 - 80%$ xuống còn $25 - 35%$.
                       ┌─────────────────────────────────────────────────────────────┐
                       │               CÁC DÒNG TIẾP CẬN KHOA HỌC                   │
                       └──────────────────────────────┬──────────────────────────────┘
                                                      │
         ┌────────────────────────────────────────────┼────────────────────────────────────────────┐
         │                                            │                                            │
         ▼                                            ▼                                            ▼
┌─────────────────────────────────┐   ┌─────────────────────────────────┐   ┌─────────────────────────────────┐
│     Động học RDX Đơn chất       │   │    Hệ Hỗn hợp Đúc RDX/TNT       │   │     Hệ Thuần hóa A-IX-1 / PBX   │
│   (Rogers 1975, Sinditskii 2008)│   │   (Cooper 1996, Oxley 2009)     │   │   (Meyer 2007, Dubovik 2005)    │
│  - E = 120-221 kJ/mol           │   │  - Trạng thái ơtecti nóng chảy  │   │  - Bao phủ vi mô hạt RDX        │
│  - Ảnh hưởng kích thước hạt     │   │  - Điểm mác cố định (Comp B)    │   │  - Giảm độ nhạy va đập cơ học   │
└────────────────┬────────────────┘   └────────────────┬────────────────┘   └────────────────┬────────────────┘
                 │                                     │                                     │
                 └─────────────────────────────────────┼─────────────────────────────────────┘
                                                       │
                                                       ▼
                       ┌─────────────────────────────────────────────────────────────┐
                       │                 KHOẢNG TRỐNG KHOA HỌC (GAP)                 │
                       │ - Thiếu mô hình hàm liên tục E(c), Z(c), D(c) khi biến đổi  │
                       │   nồng độ cấu tử đúc (TГ: 23-60%) và nén (A-IX-1: 5.0-6.5%).│
                       │ - Chưa định lượng tương quan giữa thời gian bán hủy t(1/2)  │
                       │   và nhiệt độ lưu kho nhiệt đới với dữ liệu phổ thực nghiệm. │
                       └──────────────────────────────┬──────────────────────────────┘
                                                      │
                                                      ▼
                       ┌─────────────────────────────────────────────────────────────┐
                       │               ĐỊNH VỊ CỦA LUẬN ÁN (2020)                    │
                       │ - Tích hợp DSC/DTA Kissinger + SEM/EDX + Bom nổ thực nghiệm. │
                       │ - Xây dựng sổ tay dữ liệu số hóa hoàn chỉnh cho hệ TГ và     │
                       │   A-IX-13 phục vụ tối ưu hóa mô phỏng nổ và an toàn lưu kho.│
                       └─────────────────────────────────────────────────────────────┘

Trong bức tranh tổng quan đó, những tranh luận khoa học (debates) diễn ra gay gắt giữa hai luồng quan điểm:

  • Quan điểm 1: Cho rằng quá trình phân hủy nhiệt của RDX trong môi trường TNT nóng chảy bị chi phối hoàn toàn bởi phản ứng tự xúc tác pha lỏng của sản phẩm phân hủy trung gian $NO_2$ và formaldehyde.
  • Quan điểm 2: Khẳng định hiệu ứng pha tạp chất làm giảm rào cản năng lượng hoạt hóa thông qua sự hòa tan một phần RDX vào mạng tinh thể TNT (độ tan từ $4,5\text{ g/100g TNT}$ ở $81^\circ\text{C}$ tăng lên $8,2\text{ g/100g TNT}$ ở $110^\circ\text{C}$), biến đổi cơ chế phản ứng thành cơ chế chuỗi gốc tự do đồng thể.

Luận án của Nguyễn Mậu Vương đã định vị xuất sắc ở điểm giao thoa này, giải quyết triệt để sự thiếu hụt dữ liệu biến thiên liên tục. So sánh với nghiên cứu quốc tế của Badgujar et al. (2008) tại Ấn Độ về hệ RDX/sáp polyme và công trình của Song et al. (2014) tại Trung Quốc về động học phân hủy Comp B, luận án đã mở rộng khảo sát trên ma trận mẫu thực tế với 15 tỷ lệ phối trộn hệ TГ (từ TГ-23,0 đến TГ-60,0) và 8 mẫu hệ A-IX-1, kết hợp đồng thời việc phân tích định tính - định lượng sản phẩm pha khí ($CO, CO_2, N_2, H_2, H_2O$) và cặn rắn sau nổ bằng sắc ký khí NARL8514, hồng ngoại JASCO 4600 và kính hiển vi điện tử quét phân giải cao kèm đầu dò EDX.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và kiểm chứng thực nghiệm 3 lý thuyết hóa lý cốt lõi:

  1. Lý thuyết phản ứng chuỗi gốc tự do phân nhánh của N.N. Semenov: Luận án áp dụng hoàn chỉnh hệ phương trình vi phân mô tả biến thiên nồng độ gốc tự do $n$: $$\frac{dn}{dt} = n_0 + (f - g)n = n_0 + \varphi n$$ Trong đó $f$ là thừa số phân nhánh, $g$ là thừa số đứt mạch. Luận án chứng minh rằng khi đưa chất thuần hóa (axit stearic, xerezin) hoặc TNT vào cấu trúc vi mô, xác suất đứt mạch tăng vọt ($\beta > \delta$), đưa hệ thống từ miền nổ ($\varphi > 0$) về miền ổn định nhiệt ($\varphi < 0$) trong điều kiện gia nhiệt dưới điểm bùng cháy ($T_{bc} = 230^\circ\text{C}$).
  2. Lý thuyết động học phản ứng phân hủy trạng thái rắn không đẳng nhiệt: Luận án chuẩn hóa việc áp dụng phương trình Kissinger: $$\ln\left(\frac{\beta}{T_p^2}\right) = \ln\left(\frac{Z \cdot R}{E}\right) - \frac{E}{R \cdot T_p}$$ và phương trình Ozawa để xác định năng lượng hoạt hóa biểu kiến $E$ và thừa số trước mũ $Z$, chứng minh phản ứng phân hủy nhiệt của hỗn hợp RDX tuân theo động học phản ứng bậc 1: $$\frac{dx}{dt} = k_T (1 - x) = Z \exp\left(-\frac{E}{RT}\right)(1 - x)$$
  3. Mô hình hóa thời gian bán hủy ($t_{1/2}$): Từ hằng số tốc độ riêng $k_T$, luận án thiết lập công thức dự báo tuổi thọ bền nhiệt: $$t_{1/2} = \frac{\ln 2}{k_T} = \frac{0,693}{Z \exp\left(-\frac{E}{RT}\right)}$$ cung cấp bằng chứng định lượng về độ suy giảm năng lượng theo thời gian bảo quản tại nhiệt độ môi trường.
                  ┌─────────────────────────────────────────────────────────────┐
                  │                 KHUNG PHÂN TÍCH TÍCH HỢP                    │
                  └──────────────────────────────┬──────────────────────────────┘
                                                 │
        ┌────────────────────────────────────────┼────────────────────────────────────────┐
        │                                        │                                        │
        ▼                                        ▼                                        ▼
┌───────────────────────────────┐┌───────────────────────────────┐┌───────────────────────────────┐
│  Lý thuyết Phản ứng Chuỗi Gốc ││  Lý thuyết Động học Không Đẳng││  Nhiệt Động lực học & Cơ chế  │
│  (Semenov / Hinshelwood)      ││  nhiệt (Kissinger - Ozawa)    ││  Truyền nổ Thủy động học      │
│  - Phân nhánh mạch (f)        ││  - Năng lượng hoạt hóa E      ││  - Cân bằng oxi Kb, hệ số A   │
│  - Đứt mạch bảo vệ (g)        ││  - Thừa số tần số Z           ││  - Mật độ nén/đúc ρ           │
│  - Trạng thái tới hạn phi > 0 ││  - Tốc độ gia nhiệt beta      ││  - Vận tốc sóng nổ Chapman-J. │
└───────────────┬───────────────┘└───────────────┬───────────────┘└───────────────┬───────────────┘
                │                                │                                │
                └────────────────────────────────┼────────────────────────────────┘
                                                 │
                                                 ▼
┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                 MÔ HÌNH DỰ BÁO VÀ TỐI ƯU HÓA                                    │
│   1. Động học bền nhiệt: t(1/2) = ln(2) / [Z * exp(-E/RT)] -> Dự báo tuổi thọ bảo quản kho      │
│   2. Năng lượng nổ: Qv(thực tế) = f(Q_RDX, Q_TNT, Delta_H_trộn) -> Tối ưu hóa áp suất nổ       │
│   3. Tốc độ nổ: D = a + b*rho -> Chuẩn hóa thông số đầu vào cho thủy mã mô phỏng số            │
└─────────────────────────────────────────────────────────────────────────────────────────────────┘

Khung phân tích độc đáo

Khung phân tích của luận án kết hợp đa chiều giữa:

  • Công thức phân tử giả định (CTPTGĐ): Tính toán lại tỷ lệ mol tương đương $C_aH_bO_cN_d$ cho từng tỷ lệ khối lượng hỗn hợp phức tạp.
  • Cân bằng oxi ($K_b$) và hệ số oxi ($A$): Đánh giá mức độ thiếu hụt oxy để dự báo thành phần sản phẩm nổ sơ cấp và thứ cấp theo mô hình Kistiakowsky - Wilson.
  • Biên giới hạn (Boundary conditions): Thiết lập rõ ràng khoảng mật độ đúc $\rho = (1,55 \div 1,70)\text{ g/cm}^3$ cho hệ TГ và mật độ ép $\rho = (1,60 \div 1,68)\text{ g/cm}^3$ cho hệ A-IX-1, dưới ngưỡng giới hạn đường kính tới hạn $d_{th} = 5\text{ mm}$.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án tuân thủ triệt để thế giới quan thực chứng (positivism) với phương pháp định lượng thực nghiệm nghiêm ngặt. Thiết kế nghiên cứu đa mức (multi-level design) được cấu trúc hóa từ vi mô (cấu trúc tinh thể, cỡ hạt, năng lượng liên kết) đến vĩ mô (nhiệt lượng nổ, tốc độ nổ thỏi thuốc). Cỡ mẫu thực nghiệm được thiết kế bao phủ toàn diện:

  • Đối với hệ TГ: 15 tỷ lệ thành phần biến thiên với bước nhảy $2,5%$ hàm lượng TNT (từ TГ-23,0; TГ-25,0; TГ-27,5... đến TГ-60,0).
  • Đối với hệ A-IX-1: 8 mẫu thí nghiệm phân nhóm rõ ràng giữa chất thuần hóa truyền thống A-IX-11 (sử dụng 100% xerezin) và chất thuần hóa cải tiến A-IX-13 (gồm $60%$ xerezin, $38,8%$ axit stearic, $1,2%$ sudan) với hàm lượng chất thuần hóa từ $5,0%$ đến $6,5%$.

Quy trình nghiên cứu rigorous

Quy trình chế tạo mẫu và thu thập dữ liệu được chuẩn hóa nghiêm ngặt:

  • Xử lý cỡ hạt: RDX tinh thể được phân tích phân bố cỡ hạt bằng thiết bị tán xạ laze Partica LA-950V2 và kiểm tra hình thái bề mặt qua kính hiển vi điện tử quét phát xạ trường FEI Nova NanoSEM 450.
  • Phối trộn công nghệ: Hệ TГ được nấu chảy trong nồi cách thủy ở $85 - 95^\circ\text{C}$, khuấy cơ học đồng tốc để RDX phân tán đồng đều trong TNT nóng chảy trước khi đúc rót. Hệ A-IX-1 được bọc hạt trong môi trường nước có gia nhiệt, sau đó sấy khô chân không và ép định hình ở áp suất thủy lực $2000\text{ kG/cm}^2$.
  • Độ tin cậy và chuẩn hóa: Các phép đo DTA và DSC được lặp lại ở 4 tốc độ gia nhiệt khác nhau ($\beta = 5; 10; 15; 20^\circ\text{C/min}$) trên thiết bị phân tích nhiệt quét vi sai DSC 131. Độ lặp lại của các đỉnh pic $T_{onset}$ và $T_p$ có sai số tuyệt đối $< 0,5^\circ\text{C}$.
┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                QUY TRÌNH THỰC NGHIỆM ĐỒNG BỘ                                    │
└────────────────────────────────────────────────┬────────────────────────────────────────────────┘
                                                 │
         ┌───────────────────────────────────────┴───────────────────────────────────────┐
         │                                                                               │
         ▼                                                                               ▼
┌─────────────────────────────────┐                             ┌─────────────────────────────────┐
│      CHẾ TẠO MẪU ĐỘC QUYỀN      │                             │    ĐO ĐẠC HÓA LÝ ĐỘNG HỌC       │
│ - Hệ TГ: Nấu chảy 85-95°C       │                             │ - DSC 131: Đo nhiệt lượng phân  │
│   Khuấy cơ học đồng tốc         │                             │   hủy, xác định Tonset, Tp      │
│ - Hệ A-IX-1: Ép thủy lực        │                             │ - DTA: Khảo sát tại 4 tốc độ    │
│   Áp suất 2000 kG/cm2           │                             │   beta = 5, 10, 15, 20 °C/min   │
└────────────────┬────────────────┘                             └────────────────┬────────────────┘
                 │                                                               │
                 └───────────────────────────────┬───────────────────────────────┘
                                                 │
                                                 ▼
┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│                              ĐO ĐẠC UY LỰC VÀ SẢN PHẨM NỔ THỰC TẾ                              │
│ - VOD-8: Đo tốc độ nổ thỏi thuốc ghép ống định chuẩn bằng cảm biến quang điện / xung điện       │
│ - DCA-5: Đo nhiệt lượng nổ trong bom nhiệt lượng kế dung tích kín                               │
│ - NARL8514 + JASCO 4600: Sắc ký khí và phổ FTIR phân tích thành phần sản phẩm nổ              │
│ - JSM-6510LV / EDX: Phân tích vi cấu trúc và phân bố nguyên tố cặn rắn                         │
└─────────────────────────────────────────────────────────────────────────────────────────────────┘

Data và phân tích

Thiết bị và phần mềm chuyên dụng được tích hợp đồng bộ:

  • Thiết bị đo tốc độ nổ VOD-8: Đo thời gian truyền sóng nổ qua các thỏi thuốc ghép trong ống nhựa định chuẩn, kích nổ bằng kíp nổ số 8 tiêu chuẩn với khoảng cách cảm biến quang học chính xác đến micro-giây.
  • Bom nhiệt lượng kế DCA-5: Xác định nhiệt lượng nổ thực tế $Q_v$ trong môi trường khí quyển chuẩn.
  • Hệ thống phân tích khí nổ NARL8514 & Phổ hồng ngoại JASCO 4600: Định lượng phần trăm thể tích các cấu tử khí nhẹ $H_2, N_2, CO, CO_2, CH_4$ và phổ hấp thụ của chất lỏng ngưng tụ ($H_2O$).
  • Kính hiển vi JSM-6510LV tích hợp đầu dò tán xạ năng lượng tia X (EDX): Xác định thành phần carbon tự do và tạp chất vô cơ trong sản phẩm rắn sau nổ.
  • Xử lý số liệu: Hồi quy tuyến tính phương trình Kissinger bằng phương pháp bình phương bé nhất, hệ số tương quan $R^2 > 0,995$ trên toàn bộ các chuỗi dữ liệu thực nghiệm.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Quy luật tương tác nhiệt và tính tương hợp hoàn hảo: Kết quả đường cong DSC chứng minh tính tương hợp cao giữa RDX và các cấu tử phối trộn. Trên giản đồ DSC của hệ TГ-60/40, pic nóng chảy của TNT xuất hiện rõ rệt ở $80,8^\circ\text{C}$ mà không làm thay đổi đặc tính đỉnh phân hủy chính của RDX ở vùng $205 - 235^\circ\text{C}$, khẳng định RDX không bị suy thoái hóa học trong quá trình gia nhiệt đúc rót.
  2. Sự biến thiên phi tuyến của năng lượng hoạt hóa ($E$) theo thành phần TNT: Đối với hệ TГ, khi hàm lượng TNT tăng từ $23%$ lên $60%$, năng lượng hoạt hóa biểu kiến $E$ và thừa số trước mũ $Z$ biến đổi có quy luật: Năng lượng hoạt hóa của TГ-60 đạt giá trị cao hơn đáng kể so với RDX nguyên chất, chứng minh TNT đóng vai trò chất cản trở không gian, kìm hãm tốc độ phát triển chuỗi gốc tự do ở giai đoạn đầu phân hủy.
  3. Ưu thế vượt trội của hệ chất thuần hóa đa cấu tử A-IX-13 so với A-IX-11: Thành phần chất thuần hóa gồm: $$60%\text{ xerezin} + 38,8%\text{ axit stearic} + 1,2%\text{ sudan}$$ cho khả năng bao phủ vi mô hạt RDX đồng đều hơn hẳn sáp xerezin đơn chất (A-IX-11). Ảnh SEM cho thấy màng film mỏng phân bố đều trên bề mặt tinh thể RDX, giúp mẫu A-IX-13 giữ vững độ bền nhiệt, tăng thời gian bán hủy $t_{1/2}$ ở $50^\circ\text{C}$ lên hàng chục năm mà không làm suy giảm đột ngột tốc độ nổ ($D$ duy trì ở mức cao $7900 - 8100\text{ m/s}$).
  4. Hệ phương trình thực nghiệm xác định tốc độ nổ ($D$) và nhiệt lượng nổ ($Q_v$): Luận án đã xây dựng tập phương trình hồi quy thực nghiệm có độ chính xác cao:
    • Đối với hệ TГ tại mật độ $\rho = 1,60\text{ g/cm}^3$: Tốc độ nổ $D$ giảm tuyến tính theo tỷ lệ phần trăm khối lượng TNT: $$D = D_{\text{RDX}} - k \cdot (%\text{TNT})$$ biến thiên trong dải $7500 - 7900\text{ m/s}$.
    • Đối với hệ A-IX-1: Tốc độ nổ phụ thuộc chặt chẽ vào mật độ nén ép $\rho$: $$D = 2,56 + 3,47\rho\quad (\text{km/s})$$ đạt cực đại $8380\text{ m/s}$ tại $\rho = 1,70\text{ g/cm}^3$.
  5. Định lượng sai biệt giữa nhiệt lượng nổ tính toán lý thuyết và đo đạc thực tế: Luận án phát hiện nhiệt lượng nổ thực tế $Q_v$ đo bằng bom DCA-5 luôn cao hơn tính toán lý thuyết Kistiakowsky từ $3 - 5%$ đối với các mẫu giàu RDX, nguyên nhân do phản ứng tỏa nhiệt thứ cấp của các gốc tự do dạng cacbon rắn vô định hình với hơi nước ở nhiệt độ cao ($T_{\text{nổ}} = 3380^\circ\text{C}$).
┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                 MA TRẬN PHÁT HIỆN THỰC NGHIỆM                                    │
├───────────────────┬───────────────────────────────┬─────────────────────────────────────────────┤
│   Đối tượng mẫu   │   Thông số động học nhiệt     │          Đặc tính uy lực nổ                 │
├───────────────────┼───────────────────────────────┼─────────────────────────────────────────────┤
│ RDX nguyên chất   │ E = 160-180 kJ/mol            │ D = 8380 m/s (tại rho = 1,70 g/cm3)         │
│ (d_tb = 50 um)    │ Z = 10^16 - 10^18 s^-1        │ Qv = 1250 kcal/kg; Độ nhạy: 72%             │
├───────────────────┼───────────────────────────────┼─────────────────────────────────────────────┤
│ Hệ đúc TГ-50      │ E tăng so với RDX đơn chất    │ D = 7570 m/s (tại rho = 1,70 g/cm3)         │
│ (50% RDX/50% TNT) │ T_onset nóng chảy = 80,8°C    │ Qv = 1180 kcal/kg; Độ co ngót thấp          │
├───────────────────┼───────────────────────────────┼─────────────────────────────────────────────┤
│ Hệ ép A-IX-13     │ T_p dịch chuyển an toàn       │ D = 7950-8100 m/s (tại rho = 1,62-1,65 g/cm3│
│ (RDX + 5% sáp hỗn)│ t(1/2) ở 50°C kéo dài vượt bậc│ Độ nhạy va đập giảm mạnh còn 25-35%         │
└───────────────────┴───────────────────────────────┴─────────────────────────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình định lượng liên tục về sự chuyển pha và động học phân hủy trong hệ thuốc nổ hỗn hợp hai hay nhiều cấu tử, làm cầu nối giữa lý thuyết hóa lý bề mặt và động lực học nổ khí quyển cao áp.
  • Về phương pháp luận: Hoàn thiện quy trình phân tích nhiệt kết hợp (DSC-DTA-SEM-EDX-GC-FTIR) có thể chuyển giao để đánh giá độ ổn định của bất kỳ hệ vật liệu năng lượng cao thế hệ mới nào (CL-20, FOX-7, HMX).
  • Về ứng dụng thực tiễn và công nghệ quốc phòng:
    • Xác lập ngưỡng nhiệt độ rót an toàn cho hệ TГ ở mức $85 - 95^\circ\text{C}$, ngăn ngừa triệt để nguy cơ quá nhiệt gây phân hủy tự xúc tác.
    • Cung cấp dữ liệu chuẩn cho các nhà máy sản xuất vũ khí trong việc chế tạo đầu đạn chống tăng B-41M, đạn pháo nổ phá và mìn lõm.
  • Về quản lý và chính sách kho tàng: Cung cấp công thức tính toán thời gian bán hủy $t_{1/2}$ theo nhiệt độ kho bảo quản, đặt nền tảng khoa học cho việc kéo dài niên hạn sử dụng và quy định điều kiện thông gió, làm mát cho các kho đạn dược quốc gia tại vùng khí hậu nhiệt đới gió mùa nóng ẩm.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án vẫn ghi nhận 3 giới hạn nghiên cứu:

  1. Giới hạn môi trường áp suất cao: Các phép đo DTA/DSC chủ yếu thực hiện ở áp suất khí quyển hoặc dòng khí trơ đẳng áp, chưa mô phỏng hết động học phân hủy ở điều kiện giam hãm áp suất siêu cao ($> 100\text{ MPa}$) như trong buồng đốt kín khi bắt đầu kích nổ.
  2. Kích thước cỡ hạt RDX: Nghiên cứu tập trung vào cỡ hạt RDX kỹ thuật công nghiệp tiêu chuẩn ($d_{tb} \approx 50\ \mu\text{m}$), chưa khảo sát so sánh sâu với RDX siêu mịn nano ($< 1\ \mu\text{m}$) vốn có cơ chế phân hủy bề mặt khác biệt.
  3. Ảnh hưởng của chu kỳ lão hóa ẩm: Dữ liệu độ bền nhiệt $t_{1/2}$ được ngoại suy từ mô hình Arrhenius nhiệt khô, chưa tích hợp đồng thời biến số độ ẩm tương đối cao ($RH > 85%$) trong điều kiện kho thực tế lâu năm.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) mở ra 4 hướng đi chiến lược:

  • Khảo sát động học phân hủy nhiệt dưới áp suất cao (High-Pressure DSC).
  • Nghiên cứu cơ chế suy giảm năng lượng của thuốc nổ hỗn hợp chứa bột nhôm nano (hệ A-IX-2, ТГАФ-5М).
  • Ứng dụng trí tuệ nhân tạo (Machine Learning) để dự báo vận tốc nổ $D$ và áp suất nổ $P_{CJ}$ từ cấu trúc mạng tinh thể 3D của các chất điều biến sáp mới.
  • Xây dựng mô hình suy thoái vật liệu nổ dưới tác động đồng thời của nhiệt độ - độ ẩm - ứng suất rung xóc cơ học.

Tác động và ảnh hưởng

Luận án của TS. Nguyễn Mậu Vương đã tạo nên bước chuyển biến quan trọng trên nhiều phương diện:

  • Tác động học thuật: Trở thành tài liệu tham khảo chuẩn mực trong các chương trình đào tạo sau đại học ngành Hóa học, Kỹ thuật Hóa liệu nổ và Công nghệ Vũ khí tại Viện Khoa học và Công nghệ Quân sự, Đại học Bách khoa Hà Nội, Học viện Kỹ thuật Quân sự. Tiềm năng trích dẫn cao trong các nghiên cứu chuyên ngành vật liệu năng lượng cao trong và ngoài nước.
  • Hiện đại hóa công nghiệp quốc phòng: Đóng góp trực tiếp vào việc làm chủ dây chuyền sản xuất RDX quy mô công nghiệp (công nghệ chuyển giao từ Italia đạt tiêu chuẩn MIL-248D), giải quyết triệt để vấn đề điều chỉnh độ nhạy va đập phù hợp với Tiêu chuẩn Quốc gia TCVN 1274:2010 và TCVN sửa đổi năm 2017.
  • Đóng góp cho an ninh quốc gia và kinh tế: Tối ưu hóa đơn công nghệ nhồi nạp hệ A-IX-13 và TГ giúp giảm thiểu tỷ lệ phế phẩm trong đúc rót đầu đạn, nâng cao độ tin cậy và uy lực xuyên phá của vũ khí nội địa, tiết kiệm ngân sách nhập khẩu vật liệu nổ quân dụng hàng triệu USD.

Đối tượng hưởng lợi

  • Nghiên cứu sinh & Nhà khoa học chuyên ngành Vật liệu Năng lượng cao: Tiếp cận bộ dữ liệu thực nghiệm toàn diện về động học nhiệt Kissinger và phương pháp luận tích hợp phân tích cấu trúc - uy lực nổ.
  • Kỹ sư R&D tại các nhà máy quốc phòng & khai khoáng: Ứng dụng trực tiếp thông số nhiệt độ đúc an toàn và tỷ lệ sáp A-IX-13 tối ưu để nâng cao chất lượng đạn dược và khối mồi nổ công nghiệp chất lượng cao (như dòng mồi nổ VE-05A).
  • Cơ quan quản lý đạn dược & Quy hoạch kho vũ khí: Sở hữu công cụ tính toán khoa học về thời hạn bảo quản an toàn dựa trên năng lượng hoạt hóa $E$ và thời gian bán hủy $t_{1/2}$, từ đó ban hành tiêu chuẩn kiểm định định kỳ chính xác.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?

Luận án đã mở rộng thành công lý thuyết phân hủy trạng thái rắn không đẳng nhiệt Kissinger vào hệ thuốc nổ đa cấu tử chứa chất thuần hóa dạng sáp phức hợp (A-IX-13). Nghiên cứu đã chứng minh bằng thực nghiệm rằng năng lượng hoạt hóa biểu kiến $E$ không phải là một hằng số cố định mà biến thiên phụ thuộc vào nồng độ chất hoạt động bề mặt, làm sáng tỏ cơ chế ức chế phân nhánh chuỗi gốc tự do của màng sáp trên bề mặt tinh thể RDX.

2. Điểm đổi mới phương pháp luận so với các nghiên cứu trước đây?

So với các nghiên cứu đơn lẻ của Cooper (1996) hay Badgujar (2008), luận án đã thiết lập quy trình kiểm chứng chéo (triangulation) đa tầng: Đo thông số nhiệt học (DSC 131/DTA) $\rightarrow$ Quan sát vi cấu trúc pha (FEI NanoSEM / EDX) $\rightarrow$ Phân tích sản phẩm cháy nổ thực tế (NARL8514 / JASCO 4600) $\rightarrow$ Đo uy lực nổ thực tế (VOD-8 / DCA-5). Sự kết hợp đồng bộ này loại bỏ triệt để các sai số giả định trong mô hình hóa.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm?

Đó là sự gia tăng vượt bậc của độ ổn định nhiệt trên mẫu A-IX-13 so với A-IX-11 dù cùng hàm lượng chất thuần hóa ($5,0 - 6,5%$). Sự có mặt của $38,8%$ axit stearic trong hỗn hợp sáp đóng vai trò chất tạo màng phân cực liên kết chặt chẽ với các tâm nitramin của RDX, làm tăng thời gian bán hủy dự báo ở $50^\circ\text{C}$ cao hơn $1,8$ lần so với việc chỉ dùng xerezin đơn thuần.

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) không?

Có. Luận án mô tả chi tiết từng bước: Từ thông số công nghệ nhiệt độ phối trộn ($85 - 95^\circ\text{C}$), áp suất ép nén định hình ($2000\text{ kG/cm}^2$), kích thước mẫu đo DTA ($1 - 2\text{ mg}$), tốc độ gia nhiệt ($\beta = 5, 10, 15, 20^\circ\text{C/min}$) đến quy cách ghép thỏi thuốc định chuẩn trên thiết bị VOD-8, cho phép các phòng thí nghiệm tiêu chuẩn quân sự tái lập kết quả với độ chính xác cao.

5. Khung chương trình nghiên cứu 10 năm được phác thảo như thế nào?

Luận án định hướng chương trình 10 năm tập trung vào: Hoàn thiện thư viện thông số nổ số hóa tích hợp vào các phần mềm mô phỏng thủy động lực học chuyên dụng (ANSYS Autodyn, LS-DYNA); nghiên cứu chế tạo các hệ thuốc nổ PBX thế hệ mới chịu nhiệt siêu cao trên nền Hexogen biến tính bề mặt; và làm chủ công nghệ sản xuất thuốc nổ composite phục vụ vũ khí công nghệ cao dẫn đường chính xác.

Kết luận

  1. Hệ thống hóa toàn diện: Luận án đã giải quyết trọn vẹn mối quan hệ phụ thuộc giữa thành phần khối lượng với động học phân hủy nhiệt và tốc độ nổ của 2 hệ thuốc nổ hỗn hợp nòng cốt trên nền RDX: Hệ đúc TГ (RDX/TNT) và hệ ép A-IX-1 (RDX/chất thuần hóa).
  2. Đột phá dữ liệu thực nghiệm: Xác lập hệ thống thông số động học nhiệt chính xác ($E, Z, k_T$) trên thiết bị DSC 131 và DTA theo phương pháp Kissinger - Ozawa cho toàn bộ dải nồng độ phối trộn.
  3. Cung cấp công cụ dự báo tuổi thọ: Thiết lập phương trình tính toán thời gian bán hủy $t_{1/2}$, mở ra phương pháp định lượng khoa học để xác định niên hạn bảo quản và điều kiện an toàn cháy nổ trong hệ thống kho tàng quân sự nhiệt đới.
  4. Chuẩn hóa công thức thực nghiệm tốc độ nổ: Xây dựng hệ phương trình thực nghiệm liên hệ giữa mật độ $\rho$, hàm lượng cấu tử và tốc độ nổ $D$, cung cấp đầu vào số liệu chuẩn xác cho việc mô phỏng hiệu ứng nổ lõm trên các phần mềm Autodyn và LS-DYNA.
  5. Khẳng định ưu thế của đơn pha chế A-IX-13: Chứng minh thành công hiệu quả vượt trội của hệ chất thuần hóa ba cấu tử ($60%$ xerezin, $38,8%$ axit stearic, $1,2%$ sudan) trong việc giảm độ nhạy va đập xuống mức an toàn ($25 - 35%$) trong khi vẫn bảo tồn tốc độ nổ cao ($> 8000\text{ m/s}$).
  6. Mở ra các hướng nghiên cứu hiện đại: Đặt nền móng vững chắc cho các công trình nghiên cứu tiếp theo về vật liệu nổ dẻo, thuốc nổ nhôm hóa năng lượng cao và công nghệ kiểm soát an toàn vũ khí tự động trong kỷ nguyên số.