Tổng quan về luận án

Quá trình đô thị hóa và tăng trưởng công nghiệp toàn cầu kéo theo sự gia tăng áp lực chất thải rắn. Theo thống kê của DEFRA (Anh), Eurostat (Châu Âu) và EPA (Hoa Kỳ), khối lượng chất thải rắn phát sinh hàng năm dao động từ hàng trăm triệu đến hàng tỷ tấn, trong đó thành phần giấy và cactông chiếm tỷ trọng lớn từ 20% đến 35,2%. Sau các công đoạn thu gom, phân loại và tái chế sơ cấp, một lượng lớn bã thải hữu cơ không thể tái chế tiếp tục hình thành—được định danh học thuật là "chất thải của chất thải" (waste of waste). Phương thức xử lý tối ưu cho dòng vật liệu này là công nghệ thiêu đốt thu hồi năng lượng (Waste-to-Energy) trên thiết bị lò đốt có vỉ lò di chuyển (travelling bed) hoặc lò đốt tĩnh hai cấp.

"Tiết kiệm sử dụng nhiên liệu hoá thạch, giảm thiểu lượng rác thải, tái chế, tái sử dụng và tái sinh năng lượng từ rác thải hiện nay đang là vấn đề được quan tâm ở mức toàn cầu hoá."

Mặc dù lò đốt tầng ghi di chuyển được ứng dụng rộng rãi trong công nghiệp nhiệt, các mô hình tính toán truyền thống thường xem lớp nhiên liệu như một khối đồng nhất hoặc môi trường xốp liên tục tĩnh, bỏ qua sự bất đồng nhất cấu trúc, hiện tượng co rút hình học, sự thay đổi độ rỗng cục bộ và động thái sụp lún cơ học của lớp hạt trong quá trình nhiệt phân. Khoảng trống nghiên cứu (research gap) cốt lõi nằm ở việc thiếu vắng một khung lý thuyết phân tích hệ thống hoàn chỉnh kết hợp mô tả toán học đa cấp độ (từ hạt đơn lẻ đến tập đoàn hạt đa phân tán) để giải mã chính xác các hiện tượng truyền nhiệt, truyền khối và động học phản ứng dị thể.

Luận án "Phân tích hệ thống quá trình đốt nhiệt phân chất thải rắn công nghiệp đặc trưng" của NCS. Lê Anh Kiên (hướng dẫn bởi PGS.TSKH. Lê Xuân Hải, GS. Vida Sharifi và GS. Jim Swithenbank) giải quyết trọn vẹn khoảng trống trên thông qua các câu hỏi và giả thuyết nghiên cứu tường minh:

  • Research Question 1 (RQ1): Các thông số cấu trúc vật liệu, độ rỗng ($\varepsilon$) và suất cấp khí sơ cấp tác động như thế nào đến động học cháy và cơ chế lan truyền màng bắt cháy trong thiết bị dạng cột nhồi?
  • Research Question 2 (RQ2): Làm thế nào để xây dựng mô hình toán học giải tích - số trị biểu diễn đồng thời quá trình biến đổi pha, truyền nhiệt bức xạ/đối lưu và sụp lún của hệ dị thể đa phân tán?
  • Hypothesis 1 (H1): Quá trình cháy của hỗn hợp chất thải rắn trên vỉ di chuyển có thể quy đổi tương đương về động thái biến thiên theo thời gian của một cột nhồi 1D chuyển động theo phương ngang mà không làm mất đi tính đại diện của trường nhiệt động.
  • Hypothesis 2 (H2): Việc phân hoạch hệ thống 5 lớp kết hợp phương trình cân bằng tính chất tập đoàn hạt cho phép dự báo chính xác diễn biến nhiệt độ, tốc độ giảm khối lượng và nồng độ khí thải với sai số kiểm định thống kê $\chi^2$ nằm trong giới hạn cho phép ($p > 0{,}05$).

Quy mô thực nghiệm của nghiên cứu được triển khai đồng thời trên hệ thống phân tích nhiệt vi sai TGA (NETZSCH TG 209, mẫu 15 mg, tốc độ gia nhiệt $10^\circ\text{C/phút}$) và tháp thực nghiệm cột nhồi đường kính $D = 200\text{ mm}$, chiều cao $H = 1{,}5\text{ m}$, cách nhiệt chịu nhiệt $1300^\circ\text{C}$, với lưu lượng không khí sơ cấp biến thiên từ $468\text{ kg/m}^2\text{h}$ đến $1404\text{ kg/m}^2\text{h}$ trên 5 loại vật liệu đặc trưng: khăn giấy, cactông, gỗ thông, cỏ tranh và rơm.

Literature Review và Positioning

Nghiên cứu về quá trình cháy và khí hóa chất thải rắn sinh khối trong thiết bị lớp đệm cố định và vỉ chuyển động đã ghi nhận nhiều công trình kinh điển:

  • Động học phân hủy nhiệt và khí hóa than: Henrich và cộng sự, Bassilakis và cộng sự tập trung vào động học phản ứng oxy hóa của than bán phân hủy từ sinh khối. Di Blasi phát triển các mô hình chuyển pha ngược chiều mô tả các mặt sóng phản ứng cục bộ.
  • Truyền nhiệt và khí động học cột nhồi: Fjellerup và cộng sự phân tích truyền nhiệt dẫn xạ qua lớp than rơm; Ryu và cộng sự, Horttanainen và cộng sự khảo sát tương quan giữa kích thước hạt, giới hạn cấp khí và tốc độ bắt cháy; Grammelis và cộng sự nghiên cứu quá trình đồng cháy sinh khối - than đá.
  • Mô phỏng số trị trên vỉ lò và cột nhồi: Yang, Swithenbank và Sharifi (Đại học Sheffield) thiết lập nền tảng tính toán trường nhiệt và phát thải pha khí cho lò đốt sinh khối công nghiệp. Shin và cộng sự (Viện KIST Hàn Quốc) phân chia vỉ lò di động thành chuỗi các cột nhồi độc lập để tính toán cân bằng năng lượng từng thể tích kiểm soát. Zhou và cộng sự mô phỏng cháy rơm trong buồng đốt tầng cố định, trong khi Van der Lans và cộng sự triển khai mô phỏng 2D cho lò ghi chuyển động thông thường.
┌────────────────────────────────────────────────────────────────────────┐
│                        TIẾN TRÌNH LITERATURE                           │
├────────────────────────────────────────────────────────────────────────┤
│  [Mô hình hạt đơn/tĩnh]          [Mô hình CFD liên tục]                │
│  Henrich, Bassilakis, Di Blasi   Shin et al. (KIST), Yang et al.       │
│  - Động học Arrhenius cơ bản     - Cân bằng pha khí/rắn                │
│  - Bỏ qua tương tác tập đoàn     - Giả định môi trường xốp đồng nhất   │
│                 │                                  │                   │
│                 └─────────────────┬────────────────┘                   │
│                                   ▼                                    │
│       [KHUNG PHÂN TÍCH HỆ THỐNG ĐA CẤP ĐỘ - LUẬN ÁN LÊ ANH KIÊN]       │
│       - Tiếp cận Hệ thống (Lớp 1 -> Lớp 5)                             │
│       - Động học biến thiên độ rỗng ε(t) & Hiện tượng sụp lún khối hạt │
│       - Tích hợp FVM + SIMPLE + TDMA trong phần mềm CIS chuyên dụng    │
└────────────────────────────────────────────────────────────────────────┘

Trong bức tranh tổng quan đó, tồn tại một tranh luận học thuật sâu sắc giữa hai trường phái:

  1. Trường phái Continuum Homogeneous: Coi lớp hạt là môi trường đồng thể xốp cố định với độ rỗng không đổi, tập trung giải quyết phương trình Navier-Stokes trung bình hóa (Yang et al., Van der Lans et al.).
  2. Trường phái Discrete/Dynamic Bed: Đòi hỏi phải tính đến sự biến dạng thể tích, sụt lún cơ học do tiêu hao pha rắn và sự hình thành các kênh khí ưu tiên (Ryu et al., Shin et al.).

Luận án định vị nghiên cứu vào nhánh cấu trúc động lực học, song tạo bước đột phá bằng cách không chỉ xem xét cơ chế thực nghiệm thuần túy mà thiết lập một khung lý thuyết phân tích hệ thống phân cấp hoàn chỉnh. So với mô hình của Shin et al. (chỉ cân bằng nhiệt - khối đơn giản trên 1 cell) và mô hình của Yang et al. (thiên về CFD pha khí buồng trên), công trình của Lê Anh Kiên giải quyết triệt để tương tác vi mô - vĩ mô của hệ dị thể đa phân tán, đặc biệt là cơ chế lan truyền sóng nhiệt bức xạ kết hợp sụp lún lớp hạt khi cháy hỗn hợp phức tạp.

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

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

Luận án tạo ra những đóng góp bản lề cho lý thuyết kỹ thuật hóa học và lý thuyết các quá trình nhiệt động học:

  • Mở rộng lý thuyết cân bằng tính chất tập đoàn hạt (Population Balance Equations): Ứng dụng phương trình vi phân chuyển vận tổng quát cho hệ dị thể đa phân tán để biểu diễn mật độ phân bố kích thước hạt, diện tích bề mặt tiếp xúc riêng ($A_s$) và độ rỗng thay đổi theo thời gian thực: $$\frac{\partial (\rho \phi)}{\partial t} + \text{div}(\rho \phi \mathbf{v}) = S_\phi$$ trong đó $\phi$ là đại lượng chuyên chở đơn vị (khối lượng, động lượng, năng lượng, cấu tử), $\mathbf{v}$ là trường vận tốc và $S_\phi$ là số hạng nguồn tương tác liên pha.
  • Hình thành hệ chỉ số động học cháy lớp hạt toàn diện: Định nghĩa và lượng hóa tường minh các thông số kiểm soát:
    • Tốc độ bắt cháy tức thời ($IR$ - Ignition Rate, $\text{kg/m}^2\text{h}$)
    • Tốc độ di chuyển màng lửa ($IFS$ - Ignition Front Speed, $\text{mm/s}$)
    • Tốc độ cháy tức thời ($BR$ - Burning Rate, $\text{kg/m}^2\text{h}$)
    • Tốc độ cháy trung bình ($ABR$ - Average Burning Rate, $\text{kg/m}^2\text{h}$)
    • Tốc độ cháy nghiêm ngặt ($SBR$ - Strict Burning Rate, $\text{kg/m}^2\text{h}$)
  • Chuyển dịch hệ hình (Paradigm Shift) trong mô tả cơ chế sụp lún: Chứng minh rằng hiện tượng mất cân bằng động lực học lớp đệm không phải là một nhiễu loạn ngẫu nhiên vô tổ chức, mà là hệ quả tất yếu của sự tiêu hao pha rắn dẫn đến vượt ngưỡng bền cấu trúc lỗ xốp dưới tác động của trọng lực và áp lực khí xuyên lớp.

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

Khung phân tích của luận án dựa trên phương pháp luận Tiếp cận Hệ thống, phân chia đối tượng lò đốt vỉ di chuyển / cột nhồi thành cấu trúc phân tầng 5 lớp (5-Level Hierarchical System Framework):

┌────────────────────────────────────────────────────────────────────────┐
│               HỆ THỐNG PHÂN CẤP 5 LỚP TRONG LUẬN ÁN                   │
├────────────────────────────────────────────────────────────────────────┤
│ Lớp 5: NGUYÊN HỆ (Toàn bộ thiết bị lò đốt vỉ di chuyển / tĩnh 2 cấp)   │
│   ▲                                                                    │
│ Lớp 4: HỆ DỊ THỂ ĐA PHÂN TÁN (Cột nhồi lớp hạt chất thải rắn đa cấu tử)│
│   ▲                                                                    │
│ Lớp 3: HỆ DỊ THỂ ĐƠN HẠT (Tương tác truyền nhiệt, truyền khối biên pha)│
│   ▲                                                                    │
│ Lớp 2: HỆ ĐỒNG THỂ (Pha khí liên tục: N2, O2, CO, CO2, H2O, CnHm)      │
│   ▲                                                                    │
│ Lớp 1: TIỂU PHÂN TỬ / MÔI TRƯỜNG BIÊN NGOÀI (Cân bằng nhiệt - áp suất) │
└────────────────────────────────────────────────────────────────────────┘
  • Lớp 1: Các điều kiện môi trường vi mô và trạng thái biên nhiệt động.
  • Lớp 2 (Hệ đồng thể): Pha khí chuyển động trong khe rỗng, tuân theo các định luật bảo toàn khối lượng, động lượng (Navier-Stokes) và phản ứng pha khí đồng thể.
  • Lớp 3 (Hệ dị thể đơn hạt): Quá trình sấy, thoát ẩm, nhiệt phân giải phóng chất bốc, phản ứng oxy hóa bề mặt hạt than (char oxidation) và truyền nhiệt bức xạ/đối lưu bề mặt.
  • Lớp 4 (Hệ dị thể đa phân tán): Tập hợp tập đoàn hạt trong cột nhồi với sự phân bố kích thước, độ xốp khối hạt $\rho_{bulk}$ và thông lượng bức xạ khối hạt.
  • Lớp 5 (Nguyên hệ): Toàn bộ không gian thiết bị buồng đốt sơ cấp và thứ cấp, liên kết dòng nhiệt và vật chất với vỉ lò di chuyển.

Điều kiện biên và ràng buộc (Boundary Conditions): Luận án xác lập biên bức xạ đáy vỉ lò với hệ số phản xạ/hấp thụ đặc trưng, dòng khí sơ cấp đồng nhất qua 700 lỗ ghi lò ($D=3\text{ mm}$), bỏ qua thất thoát nhiệt dẫn truyền ngang giữa các thanh ghi lân cận dựa trên phân tích tỷ số kích thước chiều cao/đường kính.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường thực chứng và hiện thực phê phán (Positivism & Scientific Realism), kết hợp chặt chẽ giữa quy trình thực nghiệm đối chứng nghiêm ngặt với mô hình hóa toán học số trị tiên tiến.

┌────────────────────────────────────────────────────────────────────────┐
│                     SƠ ĐỒ PHƯƠNG PHÁP NGHIÊN CỨU                       │
├────────────────────────────────────────────────────────────────────────┤
│  [THỰC NGHIỆM VI MÔ]               [THỰC NGHIỆM VĨ MÔ]                 │
│  TGA NETZSCH TG 209 (15 mg)        Cột nhồi H=1.5m, D=200mm            │
│  - Xác định Ea, k0 (Arrhenius)     - 11 Cặp nhiệt Ni-Cr/Ni-Al          │
│  - Phân tích Proximate & Ultimate  - Cân điện tử ghi khối lượng liên tục│
│                 │                                  │                   │
│                 └─────────────────┬────────────────┘                   │
│                                   ▼                                    │
│                       [MÔ PHỎNG SỐ TRỊ (CIS)]                          │
│                       - Phương pháp Thể tích Hữu hạn (FVM)             │
│                       - Thuật toán SIMPLE & TDMA                       │
│                       - Kiểm định thống kê Chi-Square (χ2)             │
└────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

  1. Phân tích cơ bản và nhiệt vi sai (Micro-scale):

    • Phân tích thành phần tổng quát (độ ẩm, tro, chất bốc, cacbon cố định) và phân tích nguyên tố (C, H, O, N, S) theo sổ tay CPE 01 (Đại học Sheffield).
    • Đo nhiệt trị bằng bom nhiệt lượng tiêu chuẩn.
    • Phân tích nhiệt trọng lượng vi sai (TGA) trên máy NETZSCH TG 209 trong môi trường khí nén oxy hóa với tốc độ tăng nhiệt chuẩn xác $10^\circ\text{C/phút}$, ghi nhận đường cong TG và DTG để trích xuất năng lượng hoạt hóa $E_a$ và hằng số tốc độ $k_0$.

  2. Thực nghiệm trên mô hình vật thể cột nhồi (Macro-scale):

    • Buồng đốt hình trụ đứng cao $1{,}5\text{ m}$, đường kính trong $200\text{ mm}$, vỏ ngoài thép không gỉ dày $30\text{ mm}$, lớp lót cách nhiệt dày $100\text{ mm}$ chịu nhiệt độ cực đại $1300^\circ\text{C}$.
    • Hệ thống ghi lò: đĩa vỉ khoan 700 lỗ thông khí ($D = 3\text{ mm}$, cao $8\text{ mm}$, dày $10\text{ mm}$).
    • Đo nhiệt độ: 11 đầu dò cặp nhiệt điện loại K (Ni-Cr/Ni-Al) bố trí dọc thân lò (từ đáy $y=0$ đến $y=560\text{ mm}$ cách nhau $80\text{ mm}$; các điểm kiểm soát trên tại $y=900, 1000, 1400\text{ mm}$).
    • Theo dõi khối lượng: toàn bộ cụm phản ứng được treo trên hệ thống cảm biến cân tải trọng điện tử kết nối máy tính ghi nhận liên tục theo chu kỳ 5 phút.
    • Đo thành phần khí: trích mẫu tại vị trí $y = 40\text{ cm}$ và $y = 56\text{ cm}$, phân tích trực tiếp nồng độ các cấu tử $\text{O}_2, \text{CO}, \text{CO}_2, \text{C}_n\text{H}_m, \text{NO}_x$.

Data và phân tích

Mô hình toán học được rời rạc hóa bằng Phương pháp Thể tích Hữu hạn (Finite Volume Method - FVM) trên lưới so le (staggered grid). Hệ phương trình phi tuyến liên kết bảo toàn khối lượng, động lượng, năng lượng và nồng độ cấu tử được giải bằng thuật toán SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) kết hợp bộ giải ma trận ba đường chéo TDMA (Tri-Diagonal Matrix Algorithm).

Toàn bộ thuật toán được đóng gói trong bộ chương trình mô phỏng chuyên dụng CIS (Combustion in Incinerator Simulator).

Độ tương thích giữa kết quả mô phỏng số trị và chuỗi dữ liệu thực nghiệm được đánh giá định lượng bằng kiểm định Khi bình phương ($\chi^2$ test): $$\chi^2 = \sum_{i=1}^n \frac{(O_i - E_i)^2}{E_i}$$ với $O_i$ là giá trị thực nghiệm đo đạc và $E_i$ là giá trị tính toán từ mô hình. Mọi kết quả kiểm định trên các chuỗi dữ liệu biến thiên nhiệt độ, độ suy giảm khối lượng và phát thải CO ở các chế độ khí cấp $Q = 468, 936, 1404\text{ kg/m}^2\text{h}$ đều thỏa mãn giả thuyết không ($H_0$), khẳng định tính vững chắc (robustness) và độ tin cậy của mô hình toán học.

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

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

  1. Quy luật phi tuyến của tốc độ màng bắt cháy ($IFS$) theo suất cấp khí: Khi lưu lượng khí sơ cấp tăng từ $Q = 468\text{ kg/m}^2\text{h}$ lên $936\text{ kg/m}^2\text{h}$ và $1404\text{ kg/m}^2\text{h}$, tốc độ bắt cháy ($IR$) và tốc độ lan truyền màng lửa ($IFS$) tăng mạnh do sự gia tăng cung cấp oxy cho phản ứng dị thể. Tuy nhiên, khi vượt qua ngưỡng tới hạn, hiệu ứng làm mát đối lưu (convective cooling) của dòng khí lạnh từ đáy bắt đầu chiếm ưu thế, làm suy giảm tốc độ gia nhiệt bức xạ ngược chiều, kìm hãm tốc độ lan truyền màng nhiệt.

  2. Cơ chế cân bằng nhiệt bức xạ trong pha rắn: Xác lập bằng chứng thực nghiệm cho thấy phương thức truyền nhiệt chủ đạo chi phối quá trình sấy và nhiệt phân lớp hạt bên dưới màng lửa là bức xạ nhiệt định hướng dọc trục từ vùng than cháy cục bộ bên trên xuống vùng nhiên liệu tươi:

    "Nhiệt thu được trong pha rắn trong một đơn vị thể tích kiểm soát bằng nhiệt sinh ra do phản ứng hoá học cộng với nhiệt mất mát do bay hơi ẩm, nhiệt trao đổi với thể tích kiểm soát bên trên và bên dưới và nhiệt trao đổi với pha khí."

  3. Hiện tượng sụp lún đột ngột của khối hạt (Bed Collapse Dynamics): Trong giai đoạn cacbon hóa mạnh, độ rỗng lớp hạt tăng cục bộ dẫn đến hiện tượng sụt lún cơ học. Dữ liệu thực nghiệm từ cặp nhiệt điện cho thấy sự dịch chuyển gián đoạn của profile nhiệt độ dọc trục tại các vị trí $y = 32\text{ cm}$ và $y = 56\text{ cm}$, tương ứng với sự giảm đột ngột chiều cao cột liệu.

  4. Động thái phát thải khí CO và tỷ lệ cấp khí tương đương ($ER$): Nồng độ khối lượng của CO đạt đỉnh trong pha nhiệt phân phân rã chất bốc cực đại trước khi giảm dần khi phản ứng oxy hóa than chiếm ưu thế. Tỷ lệ tương đương $ER$ (Equivalence Ratio) biến thiên mạnh theo thời gian, chứng minh tính cấp thiết của việc cấp khí thứ cấp tại vùng trên vỉ ($y = 75\text{ cm}$) để thiêu đốt hoàn toàn các hợp chất hydrocarbon chưa cháy.

┌────────────────────────────────────────────────────────────────────────┐
│          BẢNG TỔNG HỢP CÁC THÔNG SỐ VẬN HÀNH & KẾT QUẢ ĐẠT ĐƯỢC        │
├──────────────────────────┬─────────────────────────────────────────────┤
│ Thông số                 │ Giá trị / Quy chuẩn thực nghiệm             │
├──────────────────────────┼─────────────────────────────────────────────┤
│ Suất cấp khí sơ cấp (Q)  │ 468 - 936 - 1404 kg/m²h                     │
│ Khối lượng nạp mẫu (m)   │ 1,0 kg (hỗn hợp giấy, cactông, sinh khối)   │
│ Kích thước buồng phản ứng│ Cao H = 1,5 m; Đường kính D = 200 mm        │
│ Cấu tạo ghi đáy lò       │ 700 lỗ thông khí (đường kính mỗi lỗ d=3 mm) │
│ Kiểm định thống kê       │ Chuẩn số Chi-square (χ²) thỏa mãn p > 0,05  │
│ Độ chính xác mô phỏng    │ Khớp biên nhiệt độ, khối lượng và khí CO    │
└──────────────────────────┴─────────────────────────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp cơ sở khoa học hoàn chỉnh để tích hợp động học phản ứng phân hủy nhiệt Arrhenius vào hệ phương trình bảo toàn chuyển vận 2 pha khí - rắn, đặt nền móng cho các nghiên cứu chuyên sâu về công nghệ nhiệt phân - khí hóa tại Đông Nam Á.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình nghiên cứu tích hợp từ phân tích nhiệt vi sai TGA quy mô micro-gram đến mô hình cột nhồi bench-scale và phần mềm mô phỏng FVM số trị, có thể chuyển giao trực tiếp cho nghiên cứu các loại chất thải nguy hại, bùn thải hoặc chất thải y tế.
  • Về mặt thực tiễn công nghệ:

    "Thiết bị lò đốt kiểu cột nhồi được nghiên cứu trong luận án có thể được triển khai ứng dụng cho thiết bị lò đốt có vỉ lò di chuyển và thiết bị lò đốt tĩnh hai cấp. Đây là hai loại lò đốt chất thải phổ biến trên thế giới và ở Việt Nam." Cho phép tối ưu hóa chiều dài ghi lò, vận tốc dịch chuyển băng tải ghi và định tuyến phân phối tỷ lệ gió cấp 1 / cấp 2 nhằm tối đa hóa hiệu suất sinh hơi phát điện và triệt tiêu dioxin/furan.

  • Về mặt chính sách & môi trường: Đóng góp luận cứ kỹ thuật cho Chiến lược quốc gia về Quản lý tổng hợp chất thải rắn, định hướng chuyển dịch từ chôn lấp truyền thống sang công nghệ Đốt rác phát điện (Waste-to-Energy), giảm áp lực quỹ đất cho các bãi chôn lấp đô thị tại TP.HCM và Hà Nội.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án cũng ghi nhận các giới hạn biên cần được tiếp tục hoàn thiện:

  1. Giả định mô hình hình học 1 chiều (1D Axial): Mô hình CIS xem xét trường truyền nhiệt và nồng độ cấu tử biến thiên chủ yếu theo phương trục đứng $y$, bỏ qua bức xạ nhiệt theo phương ngang và gradient nhiệt độ cục bộ tại sát vách lò cách nhiệt.
  2. Xấp xỉ dạng hạt đẳng hướng: Các cấu tử rác thải (khăn giấy, cactông, rơm) trong thực tế có dạng dẹt, phi hình cầu phức tạp; việc quy đổi về kích thước đặc trưng hình cầu $l_s$ có thể dẫn đến sai số cục bộ trong tính toán diện tích bề mặt riêng $A_s$.
  3. Giản lược chuỗi phản ứng cracking hắc ín (Tar cracking): Các phản ứng nhiệt phân cracking bậc hai của hắc ín trong pha khí buồng trên chưa được mô tả chi tiết theo từng cơ chế gốc tự do.

Chương trình nghiên cứu 10 năm tiếp theo:

  • Nâng cấp mô hình CIS lên không gian 3 chiều (3D CFD - Discrete Phase Model / DEM) mô phỏng chính xác tương tác va chạm cơ học giữa các mảnh rác dị hình khi sụp lún.
  • Khảo sát quá trình đồng thiêu đốt (Co-incineration) giữa bã thải rắn công nghiệp và bùn thải công nghiệp có độ ẩm cao.
  • Nghiên cứu cơ chế chuyển hóa và kết tụ tro xỉ (Ash slagging/fouling) đối với sinh khối giàu kiềm (KCl, $K_2O$) trên dàn ống sinh hơi.
  • Đánh giá vòng đời phát thải (LCA) tích hợp hệ thống thu hồi nhiệt và xử lý khí thải khô/bán khô.

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

  • Tác động học thuật quốc tế: Luận án là minh chứng cho sự hợp tác nghiên cứu hiệu quả giữa Đại học Quốc gia TP.HCM và Đại học Sheffield (Vương quốc Anh). Các phát hiện khoa học đã được công bố và khẳng định vị thế học thuật tại các diễn đàn uy tín: Hội nghị Kỹ thuật Hóa học Khu vực (RSCE 2006, Singapore), Hội nghị Kỹ thuật Hóa học Malaysia lần thứ 21 (SOMCHE 2007), Tạp chí Kỹ thuật Hóa học ASEAN (ASEAN Journal of Chemical Engineering, 2008), và Hội nghị Kỹ thuật Hóa học ASEAN lần 2 (2009).
  • Chuyển giao và ứng dụng thực tiễn:

    "Phương pháp nghiên cứu này đã được triển khai ứng dụng xây dựng mô tả toán học trong đề tài 'Nghiên cứu ứng dụng công nghệ đốt nhiệt phân để xử lý chất thải rắn nguy hại tại TP.HCM' và đề tài 'Nghiên cứu ứng dụng quá trình nhiệt phân để xử lý thành phần hữu cơ trơ trong chất thải rắn sinh hoạt tại TP.HCM theo hướng sản xuất nhiên liệu'." Kết quả tính toán đã hỗ trợ thiết kế chế tạo lò đốt chất thải công nghiệp $4\text{ tấn/ngày}$ tại TP.HCM và cung cấp thông số vận hành cho các dự án xử lý rác công nghiệp trọng điểm tại Hà Nội (Nam Sơn - Sóc Sơn) và Bà Rịa - Vũng Tàu (Công ty Sông Xanh).

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận một cấu trúc phương pháp luận hoàn chỉnh về Tiếp cận Hệ thống (System Approach), kỹ thuật mô hình hóa FVM và phương pháp kiểm định thống kê $\chi^2$ trong cơ học chất lưu và nhiệt động lực học phản ứng.
  • Giảng viên & Chuyên gia Kỹ thuật Hóa học / Nhiệt năng: Sở hữu bộ số liệu thực nghiệm chuẩn hóa về động học phân hủy nhiệt (TGA), nhiệt trị và thông số cháy của 5 nhóm sinh khối và rác thải điển hình.
  • Kỹ sư R&D và Nhà máy Chế tạo Lò đốt: Ứng dụng thuật toán CIS để tính toán thiết kế kích thước buồng đốt, phân vùng cấp gió sơ cấp/thứ cấp, tối ưu hóa quá trình thu hồi nhiệt tạo hơi phát điện.
  • Cơ quan Quản lý Nhà nước & Quy hoạch Đô thị: Có cơ sở khoa học để ban hành quy chuẩn kỹ thuật vận hành lò đốt chất thải sinh hoạt/công nghiệp, định giá công nghệ và hoạch định chính sách phân loại rác tại nguồn.

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ì và đã mở rộng lý thuyết nào? Trả lời: Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Cân bằng Tính chất Tập đoàn Hạt (Population Balance Equations) kết hợp với Phương pháp luận Tiếp cận Hệ thống 5 lớp cho hệ dị thể đa phân tán. Luận án đã liên kết động học phân hủy nhiệt Arrhenius vi mô với phương trình bảo toàn chuyển vận vĩ mô trong điều kiện độ rỗng $\varepsilon(t)$ và thể tích khối hạt biến thiên phi tuyến.

2. Đột phá phương pháp luận của nghiên cứu so với các công trình quốc tế trước đây? Trả lời: Khác với Shin et al. (Viện KIST, Hàn Quốc) chỉ lập mô hình cân bằng nhiệt - khối đơn giản trên thể tích kiểm soát tĩnh, và Yang et al. (Đại học Sheffield) chủ yếu tập trung vào CFD pha khí buồng trên, nghiên cứu của Lê Anh Kiên đã lượng hóa đồng thời trường bức xạ nhiệt dọc trục trong lớp hạt, cơ chế sụp lún cơ học của lớp đệm và kiểm định toàn diện bằng chuẩn số Khi bình phương ($\chi^2$).

3. Phát hiện thực nghiệm nào gây bất ngờ nhất về mặt nhiệt động học? Trả lời: Sự tồn tại của ngưỡng cấp khí tới hạn chi phối tốc độ lan truyền màng bắt cháy ($IFS$). Tăng lưu lượng khí sơ cấp $Q$ ban đầu làm tăng mạnh $IFS$, nhưng khi $Q > 936\text{ kg/m}^2\text{h}$, hiệu ứng làm mát đối lưu của dòng khí lạnh bắt đầu triệt tiêu nhiệt bức xạ truyền ngược, làm suy giảm tốc độ bắt cháy cục bộ.

4. Nghiên cứu có cung cấp giao thức tái lập (Replication Protocol) hoàn chỉnh không? Trả lời: Có. Toàn bộ thông số hình học lò ($H=1{,}5\text{ m}, D=200\text{ mm}$), cấu tạo ghi lò (700 lỗ, $\phi 3\text{ mm}$), quy trình vận hành TGA NETZSCH TG 209 ($10^\circ\text{C/phút}$), tọa độ 11 cặp nhiệt điện Ni-Cr/Ni-Al và mã nguồn/lưu đồ thuật toán FVM-SIMPLE-TDMA đều được công bố chi tiết trong thuyết minh và phụ lục luận án.

5. Chương trình nghị sự 10 năm tới giải quyết vấn đề cốt lõi nào? Trả lời: Mở rộng mô hình lên không gian 3D Discrete Element Method (DEM-CFD) để giải quyết bài toán biến dạng phi đối xứng của các mảnh rác thải công nghiệp phức tạp và kiểm soát hiện tượng kết tụ tro kiềm (slagging) khi đồng thiêu đốt sinh khối với chất thải nguy hại.

Kết luận

"Đây là nghiên cứu đầu tiên ở trong và ngoài nước cho đến nay đã vận dụng nhất quán phương pháp Phân tích Hệ thống với tư cách là một tác vụ chiến lược của phương pháp luận Tiếp cận Hệ thống để nghiên cứu quá trình đốt vật liệu rắn 'chất thải của chất thải' trên đối tượng thiết bị lò đốt kiểu cột nhồi."

Tựu trung lại, công trình của NCS. Lê Anh Kiên đã xác lập những dấu ấn khoa học then chốt:

  1. Xác lập thành công phương pháp luận Tiếp cận Hệ thống 5 lớp, giải quyết thấu đáo tính phức tạp của lò đốt vỉ di chuyển và lò đốt tĩnh hai cấp.
  2. Xây dựng hệ phương trình vi phân đạo hàm riêng mô tả chính xác quá trình truyền nhiệt, truyền khối và biến đổi pha trong hệ dị thể đa phân tán.
  3. Khám phá và lượng hóa các quy luật động học cháy thông qua hệ chỉ số mới: $IFS, IR, BR, ABR, SBR$ trên nhiều hệ vật liệu rác thải đặc trưng.
  4. Phát triển thành công phần mềm mô phỏng số trị CIS dựa trên phương pháp Thể tích Hữu hạn (FVM) và thuật toán SIMPLE/TDMA, được chứng minh tương thích hoàn hảo qua kiểm định thống kê $\chi^2$.
  5. Mở ra hướng nghiên cứu mới về công nghệ thiêu đốt thu hồi năng lượng (Waste-to-Energy) cho các dòng bã thải sau phân loại tại Việt Nam và khu vực ASEAN.
  6. Đóng góp trực tiếp vào công tác thiết kế chế tạo thiết bị lò đốt công nghiệp thực tế, tạo nền tảng vững chắc cho các chính sách bảo vệ môi trường và phát triển kinh tế tuần hoàn.