Tổng quan về luận án

Nghiên cứu công nghệ tách khí và làm giàu nitơ đóng vai trò chiến lược trong nền công nghiệp hiện đại và kỹ thuật quân sự. Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Hóa học (Mã số: 9520301) với đề tài "Nghiên cứu mô phỏng và tối ưu hóa máy tạo khí nitơ sử dụng công nghệ hấp phụ dao động áp suất" do nghiên cứu sinh Phạm Văn Chính thực hiện tại Học viện Khoa học và Công nghệ Quân sự (2021), dưới sự hướng dẫn của PGS.TS. Vũ Đình Tiến và TS. Lê Quang Tuấn, đã giải quyết các nút thắt cốt lõi về mặt lý thuyết và thực nghiệm trong việc tự chủ công nghệ sản xuất khí nitơ độ tinh khiết cao.

Trước đây, công nghệ sản xuất nitơ chủ yếu dựa vào phương pháp chưng cất phân đoạn không khí hóa lỏng ở nhiệt độ siêu thấp (Cryogenic Air Separation) do Carl Von Linde phát minh từ năm 1895. Dù tạo ra sản lượng lớn, phương pháp hóa lỏng đòi hỏi chi phí đầu tư hạ tầng khổng lồ, áp suất vận hành rất cao và phụ thuộc hoàn toàn vào hệ thống vận chuyển bình chứa lỏng. Khoảng trống nghiên cứu (research gap) đặt ra là: Trong bối cảnh công nghệ hấp phụ dao động áp suất (Pressure Swing Adsorption - PSA) trên thế giới đã phát triển mạnh, việc xây dựng một mô hình toán học tích hợp đầy đủ quy luật động học, truyền khối, biến thiên áp suất và khuếch tán dọc trục để mô phỏng chính xác tháp hấp phụ đơn và hệ hai tháp sử dụng vật liệu rây phân tử carbon (Carbon Molecular Sieve - CMS) tại Việt Nam vẫn chưa được hệ thống hóa toàn diện. Luận án tập trung giải quyết 3 câu hỏi nghiên cứu và 3 giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế khuếch tán và biến thiên áp suất riêng phần của $O_2$ trong lớp vật liệu CMS-240 biến đổi như thế nào theo thời gian và chiều cao tháp hấp phụ?
  • Câu hỏi nghiên cứu 2 (RQ2): Mô hình toán học cân bằng động học nào phản ánh chính xác nhất tương tác pha khí - pha rắn trong chu trình PSA 4 bước (Skarstrom) và 6 bước (Berlin)?
  • Câu hỏi nghiên cứu 3 (RQ3): Bộ thông số công nghệ tối ưu (áp suất nạp, lưu lượng, thời gian chu kỳ) nào cho phép hệ thống PSA đạt năng suất cực đại $F_{p,\max}$ với độ tinh khiết $N_2 \ge 99.5%$ phục vụ sản xuất công nghiệp và quốc phòng?

Các giả thuyết nghiên cứu tương ứng:

  • Giả thuyết 1 (H1): Quá trình hấp phụ chọn lọc động học của $O_2$ trên rây phân tử CMS-240 bị chi phối bởi khuếch tán Knudsen ($D_k$) và khuếch tán phân tử ($D_{AB}$), tạo ra vùng truyền khối (Mass Transfer Zone - MTZ) di chuyển tịnh tiến dọc trục tháp.
  • Giả thuyết 2 (H2): Mô hình động học lực truyền tuyến tính (Linear Driving Force - LDF) kết hợp phương trình đẳng nhiệt Bi-Langmuir mô tả chính xác sự sụt giảm áp suất riêng phần $p_i(z,t)$ của $O_2$ và $N_2$.
  • Giả thuyết 3 (H3): Tối ưu hóa chu trình PSA 6 bước trên phần mềm Aspen Adsorption, Matlab và Presto cho phép thiết lập chế độ vận hành cân bằng áp suất, nâng hiệu suất thu hồi nitơ lên mức $21 - 48%$ ở dải áp suất $5 - 8\text{ bar}$.

Phạm vi nghiên cứu bao quát từ đặc tính cấu trúc nano của vật liệu CMS-240 (kích thước mao quản $\approx 4\text{ \AA}$), chế tạo hệ thiết bị thực nghiệm quy mô phòng thí nghiệm đến nâng cấp quy mô công nghiệp (scale-up) cho máy tạo khí $N_2$ công suất $50\text{ L/phút}$ ($3\text{ m}^3\text{/h}$) tại điều kiện tiêu chuẩn, ứng dụng trực tiếp vào dây chuyền sản xuất thuốc nổ RDX (Cyclotrimethylenetrinitramine).

Literature Review và Positioning

Nghiên cứu công nghệ phân tách không khí trên thế giới được chia thành ba nhánh công nghệ chính: chưng cất phân đoạn hóa lỏng sâu (Linde, 1895; Barron, 1985), công nghệ màng polymer thẩm thấu chọn lọc (Spillman, 1989; Baker, 2004) và công nghệ hấp phụ chu kỳ (Ruthven, 1984; Yang, 1987; Sircar, 2002). Trong lĩnh vực hấp phụ, các chu trình dao động nhiệt độ (TSA), dao động chân không (VSA) và dao động áp suất (PSA) đại diện cho các phương thức giải hấp khác nhau dựa trên động lực học áp suất và nhiệt độ.

                     ┌─────────────────────────────────────────────────────────┐
                     │          CÁC CÔNG NGHỆ TÁCH KHÍ TỪ KHÔNG KHÍ            │
                     └────────────────────────────┬────────────────────────────┘
                                                  │
         ┌────────────────────────────────────────┼────────────────────────────────────────┐
         │                                        │                                        │
┌────────┴────────┐                      ┌────────┴────────┐                      ┌────────┴────────┐
│ Chưng cất hóa   │                      │ Thẩm thấu qua   │                      │ Kỹ thuật hấp    │
│ lỏng sâu        │                      │ màng polymer    │                      │ phụ phân tử     │
│ (Cryogenic)     │                      │ (Membrane)      │                      │ (Adsorption)    │
└─────────────────┘                      └─────────────────┘                      └────────┬────────┘
                                                                                           │
                                                  ┌────────────────────────────────────────┼────────────────────────────────────────┐
                                                  │                                        │                                        │
                                         ┌────────┴────────┐                      ┌────────┴────────┐                      ┌────────┴────────┐
                                         │  Dao động nhiệt │                      │ Dao động chân   │                      │ Dao động áp     │
                                         │   độ (TSA)      │                      │ không (VSA)     │                      │ suất (PSA)      │
                                         └─────────────────┘                      └─────────────────┘                      └────────┬────────┘
                                                                                                                                    │
                                                                                           ┌────────────────────────────────────────┴────────┐
                                                                                           │                                                 │
                                                                                  ┌────────┴────────┐                               ┌────────┴────────┐
                                                                                  │ Chu trình 4 bước│                               │ Chu trình 6 bước│
                                                                                  │ (Skarstrom)     │                               │ (Berlin)        │
                                                                                  └─────────────────┘                               └─────────────────┘

Trong dòng y văn quốc tế, tồn tại hai luồng quan điểm đối lập về cơ chế phân tách khí trên vật liệu xốp:

  1. Trường phái chọn lọc cân bằng pha (Equilibrium selectivity): Tiêu biểu với vật liệu Zeolite (3A, 4A, 5A, 13X hình thành từ khung tứ diện $AlO_4$ và $SiO_4$ tạo lồng sodalite $\beta$-cage). Zeolite ưu tiên hấp phụ $N_2$ mạnh hơn $O_2$ do tương tác tĩnh điện giữa cation kim loại ($Na^+$, $K^+$, $Ca^{2+}$) với mômen tứ cực (quadrupole moment) của phân tử $N_2$.
  2. Trường phái chọn lọc động học khuếch tán (Kinetic selectivity): Tiêu biểu với vật liệu rây phân tử carbon (CMS). Mặc dù dung lượng hấp phụ cân bằng của $N_2$ và $O_2$ trên CMS gần như tương đương ở trạng thái bão hòa vô hạn, đường kính động học của $O_2$ ($3.46 \times 3.98\text{ \AA}$, kích thước hiệu dụng $\approx 3.8\text{ \AA}$) nhỏ hơn $N_2$ ($3.64 \times 4.1\text{ \AA}$, kích thước hiệu dụng $\approx 4.3\text{ \AA}$). Do đó, tốc độ khuếch tán của $O_2$ vào các vi mao quản $4\text{ \AA}$ nhanh hơn $N_2$ từ 20 đến 30 lần.

So sánh với các công trình quốc tế tiêu biểu:

  • Nghiên cứu của Ruthven, Farooq và Knaebel (1987-1993) đã đặt nền móng cho mô hình LDF trên tháp PSA nhưng chủ yếu giả định tháp đẳng nhiệt, bỏ qua biến thiên hệ số khuếch tán dọc trục ($D_L$) phụ thuộc vận tốc dòng chảy thực tế.
  • Công trình của Sircar và Golden (2000) cùng các nghiên cứu tại Trung Quốc (Jiangxi Xintao, Dalian University of Technology) tập trung vào chế tạo vật liệu CMS nhưng chưa tích hợp mô hình mô phỏng chuyển tiếp tự động giữa các bước nạp áp, hấp phụ, cân bằng áp suất và thổi rửa ngược dòng (countercurrent purge).

Luận án của Phạm Văn Chính định vị chính xác tại điểm giao thoa: kết hợp thực nghiệm phân tích cấu trúc vi mô (SEM, FE-SEM 200k, EDS, DSC, TG, BET) của CMS-240 với mô hình toán học giải tích - số trị mô tả áp suất riêng phần $p_i(z,t)$, từ đó liên kết phần mềm Presto và Aspen Adsorption để điều khiển hệ SCADA công nghiệp.

Đó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 thuyết va chạm chất khí và mô hình hấp phụ đơn lớp của Irving Langmuir (1918) sang hệ hỗn hợp khí đa cấu tử dưới áp suất động học. Theo lý thuyết động học chất khí, tốc độ va chạm bề mặt $R_s$ được xác định bởi: $$R_s = \frac{P}{\sqrt{2\pi M R_g T}}$$

Luận án đã tích hợp hệ số dính bám $\alpha$ và độ che phủ bề mặt $\theta = q/q_0$ để xây dựng tốc độ hấp phụ thực tế $R_a$: $$R_a = \alpha R_s (1 - \theta) = \frac{\alpha P}{\sqrt{2\pi M R_g T}}(1 - \theta)$$

Kết hợp với phương trình tốc độ giải hấp Arrhenius $R_d = k_d \theta = (k_{d\infty} e^{-\frac{E_d}{R_g T}})\theta$, luận án chứng minh hằng số cân bằng hấp phụ $K_L$ tỷ lệ nghịch với căn bậc hai của khối lượng phân tử $M$ và nhiệt độ $T$, khẳng định bản chất hấp phụ vật lý thuần túy không có rào cản năng lượng hoạt hóa ($E_d \approx \Delta H$).

             ┌─────────────────────────────────────────────────────────┐
             │       HỆ THỐNG PHƯƠNG TRÌNH ĐỘNG HỌC & TRUYỀN KHỐI      │
             └────────────────────────────┬────────────────────────────┘
                                          │
         ┌────────────────────────────────┼────────────────────────────────┐
         │                                │                                │
┌────────┴────────┐              ┌────────┴────────┐              ┌────────┴────────┐
│ Cân bằng đẳng   │              │ Động học lực    │              │ Truyền khối     │
│ nhiệt Langmuir  │              │ truyền LDF      │              │ dòng khí        │
│    mở rộng      │              │                 │              │   dọc trục      │
└────────┬────────┘              └────────┬────────┘              └────────┬────────┘
         │                                │                                │
         │   qm,i · b_i · P_i             │   ∂qi                          │   ∂Ci      ∂^2Ci     ∂(u·Ci)
  q*_i = ─────────────────                │   ─── = ω_i·(q*_i - q_i)       │   ─── = DL·───── - ────────
         1 + Σ (b_j · P_j)                │   ∂t                           │   ∂t        ∂z^2       ∂z
         │                                │                                │
         └────────────────────────────────┼────────────────────────────────┘
                                          │
                                          ▼
             ┌─────────────────────────────────────────────────────────┐
             │     MÔ HÌNH TOÁN HỌC ÁP SUẤT RIÊNG PHẦN p_i(z, t)       │
             │           TRONG THÁP HẤP PHỤ CỐ ĐỊNH                    │
             └─────────────────────────────────────────────────────────┘

Mô hình lý thuyết được thiết lập dựa trên hệ 3 phương trình vi phân đạo hàm riêng kết nối pha khí và pha rắn:

  1. Phương trình đẳng nhiệt hấp phụ Bi-Langmuir mở rộng: $$q_i^* = \frac{a_{i,1} b_{i,1} P_i}{1 + \sum b_{j,1} P_j} + \frac{a_{i,2} b_{i,2} P_i}{1 + \sum b_{j,2} P_j}$$
  2. Phương trình động học LDF (Glueckauf, 1955): $$\frac{\partial q_i}{\partial t} = \omega_i (q_i^* - q_i) = \frac{15 D_{e,i}}{R_p^2}(q_i^* - q_i)$$
  3. Phương trình bảo toàn vật chất pha khí dọc trục tháp: $$\frac{\partial C_i}{\partial t} + \frac{\partial (u_c C_i)}{\partial z} + \left(\frac{1-\varepsilon_t}{\varepsilon_t}\right)\rho_s \frac{\partial q_i}{\partial t} = D_L \frac{\partial^2 C_i}{\partial z^2}$$

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

Khung phân tích của luận án tích hợp đồng thời 4 trụ cột lý thuyết: Thuyết hấp phụ Langmuir, Thuyết khuếch tán Knudsen/phân tử, Phương trình động học LDF và Phương trình động lực học dòng chảy Ergun. Cách tiếp cận mới lạ nằm ở việc chuyển đổi biến nồng độ $C_i(z,t)$ sang trường áp suất riêng phần $p_i(z,t)$ thông qua phương trình trạng thái khí lý tưởng $P_i = C_i R_g T$, giúp việc kiểm soát quá trình nạp và xả áp trên thiết bị thực tế trở nên trực quan.

Điều kiện biên và điều kiện ban đầu được chuẩn hóa nghiêm ngặt:

  • Tại $t = 0$: $P(z,0) = P_{atm}$, $y_{O_2}(z,0) = 0.2095$, $y_{N_2}(z,0) = 0.7808$.
  • Tại đầu vào ($z = 0$): Dòng khí cấp $u_c C_{i,in} = u_c C_i|{z=0} - D_L \left.\frac{\partial C_i}{\partial z}\right|{z=0}$.
  • Tại đầu ra ($z = L$): Điều kiện Danckwerts $\left.\frac{\partial C_i}{\partial z}\right|_{z=L} = 0$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ triết lý thực chứng (positivism) với phương pháp tiếp cận đa tầng (multi-level design), kết hợp giữa phân tích vi cấu trúc vật liệu xốp, mô hình hóa toán học giải tích, mô phỏng số trên máy tính và kiểm chứng thực nghiệm trên hệ pilot.

Vật liệu nghiên cứu là Carbon Molecular Sieve CMS-240 (xuất xứ JXT - Trung Quốc), dạng viên nén hình trụ màu đen. Thiết kế tháp hấp phụ đơn và hệ hai tháp hấp phụ song song được chế tạo bằng thép không gỉ chịu áp suất làm việc từ $0$ đến $8\text{ bar}$ ($0.8\text{ MPa}$).

┌────────────────────────────────────────────────────────────────────────┐
│             QUY TRÌNH THỰC NGHIỆM VÀ PHÂN TÍCH ĐA PHƯƠNG PHÁP          │
└───────────────────────────────────┬────────────────────────────────────┘
                                    │
    ┌───────────────────────────────┼───────────────────────────────┐
    │                               │                               │
┌───┴──────────────────────┐  ┌─────┴────────────────────┐  ┌───────┴──────────────────────┐
│  Đặc trưng hóa CMS-240   │  │   Hệ thiết bị PSA Pilot  │  │   Phần mềm mô phỏng & SCADA  │
│  - SEM, FE-SEM (200k)    │  │   - Tháp đôi B1 & B2     │  │   - Matlab (ODE/PDE solver)  │
│  - EDS-JED 2300          │  │   - Van V1 - V8, VN1-VN2 │  │   - Presto & Aspen Adsorption│
│  - DSC, TG & BET         │  │   - Cảm biến P, F, %N2   │  │   - SCADA thu thập dữ liệu   │
└───┬──────────────────────┘  └─────┬────────────────────┘  └───────┬──────────────────────┘
    │                               │                               │
    └───────────────────────────────┼───────────────────────────────┘
                                    │
                                    ▼
┌────────────────────────────────────────────────────────────────────────┐
│  TỐI ƯU HÓA THÔNG SỐ: P_feed = 5-8 bar, t_cycle = 120s, N2 >= 99.5%   │
└────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

  1. Khảo sát đặc tính cấu trúc hạt CMS-240:

    • Kính hiển vi điện tử quét phát xạ trường FE-SEM ở độ phóng đại $1\times 200.000$ lần nhằm xác định cấu trúc mạng lưới vi mao quản siêu nhỏ $0.4 - 0.9\text{ nm}$.
    • Phổ tán xạ năng lượng tia X (EDS-JED 2300) kiểm tra thành phần nguyên tố tại lớp vỏ bề mặt và lớp lõi bên trong hạt, xác nhận thành phần carbon nguyên chất chiếm tỷ trọng tuyệt đối.
    • Phân tích nhiệt quét vi sai (DSC) và nhiệt trọng lượng (TG) đánh giá độ bền nhiệt và độ ổn định của vật liệu trong môi trường khí nén.
    • Đo bề mặt riêng đẳng nhiệt hấp phụ theo phương pháp BET.

  2. Quy trình đo đạc và thu thập dữ liệu hệ thống:

    • Hệ thống van điện từ điều khiển tự động gồm $V_1, V_2, V_3, V_4, V_5, V_6, V_7, V_8$ cùng hai van vi chỉnh lưu lượng xả ngược $VN_1, VN_2$.
    • Giám sát thời gian thực thông qua phần mềm SCADA chuyên dụng với giao diện cài đặt tham số chu kỳ, ghi nhận liên tục áp suất đỉnh tháp, đáy tháp, lưu lượng dòng nạp $F_1$, dòng sản phẩm $P_1$ và nồng độ cảm biến $O_2/N_2$ đầu ra.

Data và phân tích

Các thông số động học và đặc tính vật lý của hạt chất hấp phụ CMS-240 và tháp đệm được xác định chính xác qua tính toán và thực nghiệm:

Thông số ký hiệu Tên thông số và đơn vị Giá trị thực nghiệm / tính toán
$d_p$ ($R_p$) Đường kính hạt (bán kính hạt) CMS-240 $1.2\text{ mm}$ ($0.6\text{ mm}$)
$\rho_b$ Khối lượng thể tích xốp (Bulk density) $650 - 680\text{ g/L}$ ($0.68\text{ g/cm}^3$)
$\rho_p$ Khối lượng thể tích hạt (Particle density) $1.05 - 1.15\text{ g/cm}^3$
Độ bền cơ học Khả năng chịu nén của hạt $\ge 60 - 100\text{ N/hạt}$
$\varepsilon_i$ Độ rỗng ngoài hạt (khoảng trống giữa các hạt) $0.36 - 0.42$
$\varepsilon_p$ Độ rỗng trong hạt (mao quản hạt) $0.32 - 0.38$
$\varepsilon_t$ Độ rỗng tổng cộng của lớp đệm $0.60 - 0.68$
$u_c$ Vận tốc dòng khí qua lớp đệm $0.05 - 0.25\text{ m/s}$
$D_{AB}$ Hệ số khuếch tán phân tử $O_2 - N_2$ $0.18 - 0.22\text{ cm}^2\text{/s}$
$D_k$ Hệ số khuếch tán Knudsen trong mao quản $1.2 \times 10^{-3}\text{ cm}^2\text{/s}$
$D_L$ Hệ số khuếch tán dọc trục tháp $1.5 - 4.2\text{ cm}^2\text{/min}$
$\Delta P_{\max}$ Độ giảm áp lớn nhất qua lớp đệm $0.15 - 0.35\text{ bar}$

Phương trình Ergun được áp dụng để tính toán độ sụt áp qua lớp đệm hạt: $$\frac{\Delta P}{L} = 150 \frac{(1-\varepsilon_i)^2}{\varepsilon_i^3}\frac{\mu u_c}{d_p^2} + 1.75 \frac{(1-\varepsilon_i)}{\varepsilon_i^3}\frac{\rho_g u_c^2}{d_p}$$

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

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

  1. Quy luật phân bố áp suất riêng phần $p_{O_2}(z,t)$ theo không gian và thời gian: Kết quả mô phỏng trên phần mềm Presto và Matlab khớp với thực nghiệm với sai số $< 4.5%$. Tại thời điểm $t = 40\text{ s}$ ở áp suất nạp $5\text{ bar}$, nồng độ oxy giảm mạnh từ $20.95%$ tại đáy tháp ($z = 0$) xuống mức $< 0.5%$ tại đỉnh tháp ($z = 0.65\text{ m}$), chứng minh vùng truyền khối MTZ được giữ hoàn toàn bên trong lớp vật liệu CMS-240.
  2. Xác lập chu trình Berlin 6 bước tối ưu hơn chu trình Skarstrom 4 bước: Việc bổ sung bước cân bằng áp suất (Pressure Equalization) giữa hai tháp $B_1$ và $B_2$ giúp tận dụng dòng khí giàu nitơ có áp suất cao từ tháp vừa kết thúc hấp phụ để nạp sơ bộ cho tháp vừa giải hấp, giảm $28%$ năng lượng tiêu thụ của máy nén khí và tăng tỷ lệ thu hồi sản phẩm từ $21%$ lên $42%$.
  3. Hiện tượng giải hấp ngược dòng tăng cường (Countercurrent Purging Effect): Khi áp suất tháp giảm từ $5\text{ bar}$ về áp suất khí quyển kết hợp thổi rửa bằng $5 - 8%$ dòng sản phẩm $N_2$ tinh khiết từ tháp đối diện, lượng $O_2$ bị giữ trong các vi mao quản được giải phóng hoàn toàn chỉ trong $35 - 40\text{ s}$, phục hồi $99.2%$ dung lượng hấp phụ ban đầu của CMS-240.
  4. Giới hạn tới hạn của vận tốc dòng nạp: Khi vận tốc dòng khí $u_c > 0.3\text{ m/s}$, hiện tượng xuyên thủng vùng truyền khối (breakthrough) xảy ra sớm trước $50\text{ s}$, khiến độ tinh khiết $N_2$ tụt dốc từ $99.5%$ xuống dưới $95.0%$.
Áp suất riêng phần P_O2 (bar)
1.0 ┼───────┐
    │       │ (z = 0 m - Khí nạp vào)
0.8 │       │
    │        \
0.6 │         \
    │          \  (Vùng truyền khối MTZ di chuyển)
0.4 │           \
    │            \
0.2 │             \
    │              └──────────────────────── (z = 0.65 m - Khí thành phẩm N2)
0.0 ┼──────────────────────────────────────────► Chiều cao tháp z (m)
    0       0.15     0.30     0.45     0.65

Implications đa chiều

  • Về mặt học thuật: Cung cấp bộ dữ liệu thực nghiệm chuẩn xác về hệ số khuếch tán dọc trục $D_L$, dung lượng hấp phụ cân bằng $q^*$ và thông số LDF của vật liệu CMS-240, phục vụ việc giảng dạy và phát triển các mô hình truyền khối nâng cao trong công nghệ hóa học.
  • Về mặt công nghệ: Xây dựng thành công quy trình thiết kế, mô phỏng và phóng to quy mô (scale-up) máy tạo khí nitơ PSA từ $10\text{ L/phút}$ lên $50\text{ L/phút}$ và $200\text{ L/phút}$ mà không cần trải qua quá trình thử - sai tốn kém.
  • Về mặt ứng dụng thực tiễn: Ứng dụng thành công hệ thống PSA công suất $50\text{ L/phút}$ ($N_2 \ge 99.5%$) trực tiếp vào trục vít cấp liệu Hexamine trong dây chuyền sản xuất thuốc nổ quân dụng RDX, loại bỏ nguy cơ cháy nổ, bảo quản khí tài quân sự và ứng dụng đóng gói bảo quản nông sản thực phẩm (Modified Atmosphere Packaging - MAP).

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  • Giả định đẳng nhiệt trong mô hình toán: Quá trình hấp phụ thực tế là quá trình tỏa nhiệt ($\Delta H < 0$). Nhiệt độ cục bộ trong tháp có thể tăng từ $2 - 5^\circ\text{C}$ khi nạp áp và giảm khi xả áp. Mặc dù ảnh hưởng này không làm biến dạng lớn đường cong nồng độ ở quy mô nhỏ, mô hình phi đẳng nhiệt (non-isothermal) cần được tích hợp khi thiết kế tháp công nghiệp siêu lớn ($> 500\text{ m}^3\text{/h}$).
  • Vật liệu rây phân tử: Luận án tập trung khảo sát trên mẫu hạt CMS-240 thương mại của JXT; chưa mở rộng so sánh song song với các dòng vật liệu rây phân tử carbon biến tính bằng lắng đọng hơi hóa học (CVD) hoặc vật liệu khung cơ kim (MOFs - Metal-Organic Frameworks).
  • Điều kiện khí đầu vào: Nghiên cứu giả định không khí đã qua xử lý lọc sạch bụi, dầu và tách ẩm hoàn toàn bằng tháp đệm Silicagel. Trong thực tế, độ ẩm dư có thể làm suy giảm tuổi thọ của CMS-240 sau 3 - 5 năm vận hành.

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

  1. Phát triển mô hình 3D động lực học chất lưu tính toán (CFD) tích hợp truyền nhiệt - truyền khối phi đẳng nhiệt và phi đoạn nhiệt.
  2. Thử nghiệm chu trình PSA kết hợp chân không (VPSA) nhằm giảm áp suất nạp xuống $2 - 3\text{ bar}$, tiết kiệm điện năng tiêu thụ.
  3. Nghiên cứu biến tính vật liệu CMS trong nước từ nguồn than antraxit Quảng Ninh bằng phương pháp hoạt hóa nhiệt than kết hợp cracking hydrocarbon.
  4. Tích hợp thuật toán trí tuệ nhân tạo (AI/Machine Learning) vào hệ SCADA để tự động điều chỉnh chu kỳ đóng mở van theo độ suy giảm của vật liệu qua thời gian.

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

Luận án tạo ra tác động sâu rộng trên nhiều phương diện:

  • Ảnh hưởng học thuật: Công bố 10 công trình khoa học chuyên ngành, bao gồm 01 bài báo quốc tế chuẩn ISI, 01 bài báo tại Hội nghị Quốc tế CASEAN-6 và 08 bài báo trên các tạp chí khoa học chuyên ngành uy tín trong nước (Tạp chí Hóa học, Tạp chí Khoa học và Công nghệ Quân sự).
  • Tự chủ công nghệ quốc phòng: Chấm dứt sự phụ thuộc vào thiết bị nhập khẩu đắt đỏ từ nước ngoài, chủ động cung cấp nguồn khí trơ tinh khiết cao phục vụ sản xuất thuốc phóng, thuốc nổ RDX, bảo dưỡng tên lửa, vũ khí công nghệ cao trong quân đội.
  • Hiệu quả kinh tế xã hội: Giảm giá thành sản xuất khí nitơ tại chỗ từ $40 - 60%$ so với việc mua bình khí nén hóa lỏng truyền thống, giảm thiểu rủi ro vận chuyển áp lực cao và thúc đẩy hiện đại hóa công nghệ bảo quản nông sản sau thu hoạch tại Việt Nam.

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

  • Nghiên cứu sinh và học viên cao học chuyên ngành Hóa kỹ thuật: Tiếp cận phương pháp luận mô hình hóa quá trình hấp phụ động học và kỹ thuật lập trình giải hệ phương trình đạo hàm riêng trên Matlab/Presto.
  • Kỹ sư vận hành và thiết kế thiết bị Hóa chất - Dầu khí: Nhận chuyển giao quy trình tính toán kích thước tháp đệm, lựa chọn van điều khiển, thiết lập chu kỳ Berlin 6 bước và cấu hình SCADA.
  • Doanh nghiệp quốc phòng và công nghiệp chế biến: Trực tiếp ứng dụng thiết bị máy tạo nitơ PSA tự động, an toàn, ổn định với chi phí đầu tư và bảo trì thấp.

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 va chạm động học chất khí Langmuir (1918) và mô hình Linear Driving Force (LDF) sang mô hình trường áp suất riêng phần $p_i(z,t)$ hai chiều không gian - thời gian, phản ánh chính xác cơ chế khuếch tán chọn lọc động học của $O_2$ trong vi mao quản $4\text{ \AA}$ của rây phân tử carbon CMS-240.

2. Điểm mới về phương pháp nghiên cứu so với các công trình trước đây?

Khác với các nghiên cứu của Ruthven (1984) hay Sircar (2000) vốn chỉ giải bài toán giải tích đơn giản hoặc khảo sát thực nghiệm thuần túy, luận án đã thiết lập quy trình tích hợp 4 bước khép kín: Đặc trưng hóa vật lý nano hạt CMS-240 (FE-SEM 200k, EDS, DSC/TG) $\rightarrow$ Mô hình hóa toán học $\rightarrow$ Mô phỏng đa phần mềm (Presto, Aspen Adsorption, Matlab) $\rightarrow$ Thực nghiệm kiểm chứng trên hệ SCADA tự động hóa hoàn toàn.

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

Tại áp suất làm việc $5\text{ bar}$, quá trình cân bằng áp suất giữa hai tháp trong $3 - 5\text{ s}$ không chỉ giúp tiết kiệm năng lượng nén mà còn làm sắc nét ranh giới vùng truyền khối MTZ, ngăn chặn hiện tượng xáo trộn nồng độ ở đầu tháp và nâng độ tinh khiết $N_2$ thành phẩm từ $98.2%$ lên $\ge 99.5%$.

4. Giao thức nhân bản thực nghiệm (Replication Protocol) có được cung cấp rõ ràng không?

Luận án cung cấp chi tiết toàn bộ thông số hình học tháp ($D \times H$), đặc tính lớp đệm ($\rho_b = 680\text{ g/L}$, $\varepsilon_t = 0.64$), lưu lượng dòng nạp $F_1$, bảng trạng thái đóng mở van theo thời gian của chu trình 6 bước (Berlin cycle) và thuật toán điều khiển SCADA, cho phép tái lập thực nghiệm hoàn toàn.

5. Định hướng nghiên cứu 10 năm tới được vạch ra như thế nào?

Lộ trình 10 năm tập trung vào: (1) Tự chủ tổng hợp vật liệu CMS chất lượng cao từ than antraxit Việt Nam; (2) Phát triển hệ thiết bị PSA di động siêu nhỏ gọn phục vụ các đơn vị quân sự dã chiến; (3) Tối ưu hóa chu trình VPSA tích hợp màng lọc hybrid nhằm nâng công suất lên quy mô hàng ngàn $m^3\text{/h}$.

Kết luận

Luận án tiến sĩ của NCS. Phạm Văn Chính đã đạt được 5 đóng góp học thuật và thực tiễn cốt lõi:

  1. Xác định toàn diện bộ thông số vi cấu trúc, độ bền nhiệt, diện tích bề mặt và thành phần nguyên tố của rây phân tử carbon CMS-240.
  2. Thiết lập thành công hệ phương trình toán học mô tả động học hấp phụ và biến thiên áp suất riêng phần $p_i(z,t)$ của $O_2$ và $N_2$ trong tháp đệm cố định.
  3. Ứng dụng hiệu quả các phần mềm mô phỏng tiên tiến (Aspen Adsorption, Presto, Matlab) để tối ưu hóa chu trình PSA 6 bước, nâng cao tỷ lệ thu hồi khí $N_2$.
  4. Thiết kế, chế tạo và thử nghiệm thành công hệ máy tạo khí nitơ PSA công suất $50\text{ L/phút}$, độ tinh khiết $N_2 \ge 99.5%$ tại điều kiện áp suất $5 - 8\text{ bar}$.
  5. Ứng dụng thực tế hệ thống vào dây chuyền sản xuất thuốc nổ quân dụng RDX, tạo tiền đề vững chắc cho việc tự chủ công nghệ phân tách khí và phát triển ngành kỹ thuật hóa chất công nghiệp tại Việt Nam.