Tổng quan luận án

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

Khí sinh học (biogas) thu được từ quá trình lên men kỵ khí các chất thải hữu cơ (rác thải sinh hoạt, phân gia súc, phụ phẩm nông nghiệp) là nguồn năng lượng tái tạo chứa thành phần chính gồm 55–70 vol% $\text{CH}_4$ và 30–45 vol% $\text{CO}_2$, cùng một lượng nhỏ tạp chất ($\text{H}_2\text{S}$, $\text{NH}_3$, $\text{H}_2\text{O}$, $\text{N}_2$, $\text{O}_2$). Việc chuyển hóa biogas thành điện năng qua các động cơ nhiệt truyền thống có hiệu suất bị giới hạn ở mức 25–30% do phải trải qua chuỗi biến đổi năng lượng đa bước (hóa năng $\rightarrow$ nhiệt năng $\rightarrow$ cơ năng $\rightarrow$ điện năng). Ngược lại, pin nhiên liệu oxit rắn (Solid Oxide Fuel Cell - SOFC) vận hành ở nhiệt độ cao (600–1000 °C, phổ biến là 700–800 °C) cho phép chuyển hóa trực tiếp hóa năng thành điện năng với hiệu suất cao, đồng thời tận dụng nhiệt lượng chất lượng cao sinh ra từ phản ứng điện hóa.

Khi cấp trực tiếp biogas vào anode gốc Niken (Ni) của SOFC, quá trình hoàn nguyên khô (dry reforming với $\text{CO}_2$) và hoàn nguyên ướt (steam reforming với $\text{H}_2\text{O}$) xảy ra đồng thời trên bề mặt hạt xúc tác kim loại Ni để tạo thành khí tổng hợp ($\text{H}_2$ và $\text{CO}$). Quá trình này được gọi là vận hành hoàn nguyên nội tại trực tiếp (Direct Internal Reforming - DIR), và chuỗi phản ứng đồng thời được định danh là quá trình đa hoàn nguyên methane (Methane Multiple-Reforming - MMR). Tuy nhiên, việc ứng dụng DIR-SOFC chạy biogas gặp các thách thức kỹ thuật:

  1. Phản ứng hoàn nguyên methane thu nhiệt mạnh gây sụt giảm nhiệt độ cục bộ tại vùng cửa nạp nhiên liệu, tạo ra gradient nhiệt và ứng suất nhiệt lớn dẫn tới nguy cơ nứt vỡ màng điện phân (electrolyte fracture).
  2. Sự suy giảm hoạt tính anode do đóng cặn carbon (coking) và ngộ độc tạp chất lưu huỳnh ($\text{H}_2\text{S}$).
  3. Thiếu hụt một mô hình động học chính xác mô tả quá trình MMR; các nghiên cứu trước đây thường coi MMR là tổng tuyến tính độc lập của phản ứng hoàn nguyên khô và hoàn nguyên ướt (bỏ qua ảnh hưởng đồng thời của $\text{CO}_2$ và $\text{H}_2\text{O}$ lên độ chuyển hóa xúc tác của $\text{CH}_4$), dẫn đến sai lệch khi mô phỏng nhiệt độ, nồng độ và ứng suất nhiệt trong tế bào pin.

Mục tiêu và các nhiệm vụ nghiên cứu

Mục tiêu chính của luận án là phát triển mô hình toán học và công cụ mô phỏng số động lực học chất lưu tính toán (Computational Fluid Dynamics - CFD) hoàn chỉnh cho tế bào đơn DIR-SOFC vận hành bằng biogas, từ đó phân tích cơ chế vận hành và đề xuất các giải pháp kỹ thuật nâng cao độ bền cơ-nhiệt và mật độ công suất. Các nhiệm vụ cụ thể được tác giả đánh số và triển khai như sau:

  1. Khảo sát thực nghiệm hành vi điện hóa và đặc tính phản ứng reforming nội tại của anode cermet gốc Ni khi cấp hỗn hợp khí biogas mô phỏng ($\text{CH}_4\text{-}\text{CO}_2\text{-}\text{N}_2$ và $\text{CH}_4\text{-}\text{CO}_2\text{-}\text{H}_2\text{O}$).
  2. Xây dựng mô hình động học mới cho quá trình đa hoàn nguyên methane (MMR) trong cấu trúc xốp của anode gốc Ni bằng phương pháp tiếp cận quy nạp kết hợp mạng nơ-ron nhân tạo (Artificial Neural Network - ANN) và hệ suy luận mờ (Fuzzy Inference System - FIS).
  3. Thiết lập mô hình CFD ba chiều (3D-CFD) toàn diện kết hợp mô hình MMR và xây dựng chiến lược thẩm định mô hình ba bước (three-step model validation) đối chiếu với dữ liệu thực nghiệm ở cả chế độ không DIR (non-DIR với $\text{H}_2$) và chế độ DIR (với biogas).
  4. Ứng dụng mô hình 3D-CFD đã thẩm định để đánh giá hiệu quả của hai khái niệm DIR cải tiến (Advanced DIR concepts) trên các cấu trúc tế bào có độ dày lớp nền anode khác nhau.

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

  • Đối tượng nghiên cứu: Tế bào pin nhiên liệu oxit rắn đơn lớp (planar single cell) dạng hỗ trợ bởi anode (Anode-Supported Cell - ASC) và dạng hỗ trợ bởi chất điện phân (Electrolyte-Supported Cell - ESC); các phản ứng hóa học (hoàn nguyên khô, hoàn nguyên ướt, chuyển dịch khí-nước WGS, phân hủy methane hình thành carbon) và phản ứng điện hóa oxy hóa $\text{H}_2$, $\text{CO}$ tại ranh giới ba pha (Triple-Phase Boundary - TPB).
  • Phạm vi không gian:
    • Tế bào dạng nút ESC diện tích anode $8 \times 8\ \text{mm}^2$, độ dày anode $\sim 38\ \mu\text{m}$ ($\text{Ni-10ScSZ}$).
    • Tế bào nửa bán thành phẩm (half-cell) hỗ trợ bởi anode đường kính $20\ \text{mm}$ ($\text{Ni-8YSZ}$).
    • Tế bào phẳng ASC kích thước $14 \times 14\ \text{mm}^2$ và $20 \times 50\ \text{mm}^2$ với hai biến thể độ dày chất nền anode: ASC-A ($d_a = 950\ \mu\text{m}$) và ASC-B ($d_a = 200\ \mu\text{m}$).
  • Phạm vi thời gian và điều kiện vận hành: Vận hành ở dải nhiệt độ 700 °C, 750 °C và 800 °C; nhiên liệu là hỗn hợp biogas mô phỏng với tỷ lệ $\text{CH}_4/\text{CO}_2 = 1$ hoặc thay đổi lưu lượng từ $40$ đến $200\ \text{mL}\cdot\text{min}^{-1}$; áp suất khí quyển ($1\ \text{atm}$).

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

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

Luận án hệ thống hóa các mốc phát triển công nghệ và các công trình nghiên cứu tiền đề:

  • Lịch sử và phát triển hệ thống SOFC: Từ nguyên lý của Sir W. Nernst (1897) và E. Baur & H. Preis (1937) với chất điện phân $\text{ZrO}_2\text{-}\text{Y}_2\text{O}_3$, công nghệ SOFC đã được thương mại hóa qua các hệ thống lớn của Siemens Westinghouse (1998, 140 kW), Mitsubishi Heavy Industries (hệ lai SOFC-MGT 200 kW đạt hiệu suất 52,1% LHV), Rolls-Royce (1 MW Integrated-Planar SOFC), dự án trình diễn SOFC-MGT 250 kW của liên doanh Mitsubishi Hitachi tại Đại học Kyushu (2015), các hệ dân dụng CHP 1–2 kW (Kyocera, Tokyo Gas, TOTO, Ceramic Fuel Cells - BlueGen đạt hiệu suất điện ban đầu 60% LHV) và xe nguyên mẫu SOFC 5 kW chạy ethanol của Nissan (2016).
  • Nghiên cứu vận hành DIR với hydrocarbon và biogas: Các công trình thực nghiệm của Yentekakis (2006), Leone (2010), Sasaki (2011), Kameshima (2011), Shiratori và cộng sự (2012–2015) chứng minh tính khả thi của việc nạp trực tiếp biogas, cồn và nhiên liệu sinh học biodiesel vào SOFC.
  • Nghiên cứu kiểm soát đóng cặn carbon và ngộ độc tạp chất: Sử dụng màng cản khuếch tán khí trơ (Barnett và cộng sự, 2005), bổ sung lớp xúc tác chống đóng cặn (Georges, 2006; Fonseca, 2007), phát triển vật liệu anode kháng carbon như Ceria-doped YSZ (Campbell, 2010; Vohs, 2010), anode gốc Cu (Ahn, 2007; Conghua, 2006), kim loại quý Ru/Pd (Shao, 2011; Watanabe, 2013; Li, 2016), và chất xúc tác cấu trúc giấy (Paper-Structured Catalyst - PSC) tích hợp khoáng hydrotalcite phân tán (Shiratori, Kitaoka, Sakamoto và cộng sự, 2013–2016).
  • Mô hình hóa động học phản ứng reforming và CFD trong SOFC:
    • Động học hoàn nguyên ướt đơn lẻ: Lehnert và cộng sự (2000), Achenbach & Riensche (1994), Haberman & Young (2004), Xu & Froment (1989).
    • Động học chi tiết 42 phản ứng sơ cấp vi mô: Deutschmann và cộng sự (2000, 2006), Janardhanan & Deutschmann (2006, 2007). Cách tiếp cận này có chi phí tính toán rất lớn khi đưa vào mô phỏng đa chiều 3D.
    • Cách tiếp cận hoàn nguyên song song (Parallel-reforming approach): Meng Ni (2012, 2013) và Girona và cộng sự (2012) tính tốc độ hoàn nguyên của hỗn hợp bằng cách cộng tách rời tốc độ hoàn nguyên ướt ($\text{CH}_4\text{-}\text{H}_2\text{O}$) và hoàn nguyên khô ($\text{CH}_4\text{-}\text{CO}_2$).
Hướng tiếp cận mô hình hóa Tác giả tiêu biểu Đặc điểm chính Hạn chế ghi nhận trong văn bản
Động học hoàn nguyên ướt đơn phân tử Lehnert et al. (2000), Haberman & Young (2004) Chỉ tính phản ứng giữa $\text{CH}_4$ và $\text{H}_2\text{O}$ Bỏ qua hoàn nguyên khô khi $\text{CO}_2$ chiếm tỷ lệ lớn
Cơ chế phản ứng sơ cấp đa bước (42 phản ứng) Deutschmann et al. (2000, 2006) Mô tả chi tiết hóa học bề mặt trên Ni Phức tạp cao, tốn tài nguyên tính toán trong CFD 3D
Hoàn nguyên song song (Parallel-reforming) Meng Ni (2012, 2013), Girona et al. (2012) Lấy tổng đại số tốc độ hoàn nguyên khô và ướt riêng rẽ Bỏ qua tương tác đồng thời giữa $\text{CO}_2$ và $\text{H}_2\text{O}$ lên tâm xúc tác
Mô hình quy nạp ANN/FIS (Luận án này) Tran Dang Long (2017) Huấn luyện từ dữ liệu thực nghiệm đa thành phần Mô tả chính xác tốc độ tiêu thụ/sinh khí ở mọi dải tỷ lệ

Khoảng trống luận án giải quyết

Luận án chỉ ra rằng cách tiếp cận hoàn nguyên song song truyền thống không phản ánh đúng thực tế vì khi cả $\text{CO}_2$ và $\text{H}_2\text{O}$ cùng hiện diện trong anode, chúng cạnh tranh và tương tác phức tạp trên bề mặt xúc tác Ni. Do đó, luận án lấp đầy khoảng trống này bằng cách:

  1. Thiết lập một mô hình toán học thực nghiệm dựa trên mạng nơ-ron ANN và hệ mờ FIS, tính toán chính xác tốc độ thực của các cấu tử trong quá trình MMR mà không cần giải hệ 42 phản ứng vi mô phức tạp.
  2. Thiết lập quy trình thẩm định 3 bước kết hợp dữ liệu điện hóa thực nghiệm, khắc phục tình trạng các công trình trước đây (Wang et al., Colpan et al., Razbani et al., Santarelli et al.) chỉ làm khớp đường cong dòng-áp ($i\text{-}V$) một cách cục bộ mà không kiểm chứng ở chế độ nạp $\text{H}_2$ cũng như phân bố nhiệt độ và thành phần khí thực tế.

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

Khung lý thuyết và các phản ứng hóa học

Quá trình vận hành DIR-SOFC sử dụng nhiên liệu gốc methane bao gồm các phản ứng hóa học và điện hóa diễn ra đồng thời trong cấu trúc xốp của anode và bề mặt phân cách:

$$\text{O}_2 + 4\text{e}^- \rightarrow 2\text{O}^{2-} \quad (\text{Phản ứng khử điện hóa tại Cathode TPB})$$

$$\text{H}_2 + \text{O}^{2-} \rightarrow \text{H}_2\text{O} + 2\text{e}^- \quad (\text{Oxy hóa điện hóa } \text{H}_2 \text{ tại Anode TPB})$$

$$\text{CO} + \text{O}^{2-} \rightarrow \text{CO}_2 + 2\text{e}^- \quad (\text{Oxy hóa điện hóa CO tại Anode TPB})$$

$$\text{CH}_4 + \text{CO}_2 \rightleftharpoons 2\text{H}_2 + 2\text{CO} \quad (\text{Hoàn nguyên khô trong Anode, } \Delta H > 0)$$

$$\text{CH}_4 + \text{H}_2\text{O} \rightleftharpoons 3\text{H}_2 + \text{CO} \quad (\text{Hoàn nguyên ướt trong Anode, } \Delta H > 0)$$

$$\text{H}_2\text{O} + \text{CO} \rightleftharpoons \text{H}_2 + \text{CO}_2 \quad (\text{Phản ứng chuyển dịch khí-nước WGS, } \Delta H < 0)$$

$$\text{CH}_4 \rightarrow \text{C} + 2\text{H}_2 \quad (\text{Nhiệt phân methane tạo carbon, } \text{coking})$$

$$\text{C} + \text{CO}_2 \rightleftharpoons 2\text{CO} \quad (\text{Phản ứng Boudouard nghịch})$$

$$\text{C} + \text{H}_2\text{O} \rightleftharpoons \text{H}_2 + \text{CO} \quad (\text{Khí hóa carbon})$$

Sức điện động lý thuyết (Nernst voltage) được xác định dựa trên chênh lệch thế hóa học của oxy giữa cathode và anode:

$$E = -\frac{\Delta G^\circ}{2F} + \frac{RT}{2F}\ln\left(\frac{p_{\text{H}2,a}\cdot p{\text{O}2,c}^{1/2}}{p{\text{H}_2\text{O},a}}\right)$$

Điện áp thực tế của tế bào pin $V_{\text{cell}}$ khi có dòng điện $i$ chạy qua bị suy giảm bởi các tổn hao nội:

$$V_{\text{cell}} = E - i(R_a + R_e + R_c) - (\eta_{\text{act},a} + \eta_{\text{act},c}) - (\eta_{\text{conc},a} + \eta_{\text{conc},c})$$

Phương pháp nghiên cứu thực nghiệm

  • Đo phổ tổng trở điện hóa (Electrochemical Impedance Spectroscopy - EIS): Thực hiện trên tế bào ESC ($\text{Ni-10ScSZ}$) ở 800 °C dưới điều kiện hở mạch (Open-Circuit Condition - OCV) với lưu lượng khí $80\ \text{mL}\cdot\text{min}^{-1}$, nhằm tách biệt điện trở truyền khối ($R_{\text{mt},a}$) ở dải tần số thấp và điện trở truyền điện tích ($R_{\text{ct},a}$) ở dải tần số trung/cao.
  • Thử nghiệm phản ứng Reforming trên Half-cell: Sử dụng đĩa cermet $\text{Ni-8YSZ}$ (đường kính $20\ \text{mm}$) đặt trong lò phản ứng thạch anh ở các mức nhiệt độ 700 °C, 750 °C và 800 °C với lưu lượng $100\ \text{mL}\cdot\text{min}^{-1}$ để thu thập ma trận tốc độ phản ứng tiêu thụ $\text{CH}4$ ($r{\text{CH}_4}$) và sinh $\text{H}2$ ($r{\text{H}_2}$) dưới các tỷ lệ khí $\text{CH}_4\text{-}\text{CO}_2\text{-}\text{H}_2\text{O}\text{-}\text{N}_2$ khác nhau.
  • Đo đạc đường đặc trưng dòng-áp ($i\text{-}V$): Đo đạc trên trạm thử nghiệm tế bào đơn ASC ($14 \times 14\ \text{mm}^2$ và $20 \times 50\ \text{mm}^2$) gắn trên vỏ đỡ alumina và hàn kín bằng keo gốm ceramic bond.

Phương pháp mô hình hóa và mô phỏng số

  • Mô hình MMR dựa trên ANN/FIS: Mạng nơ-ron ANN cấu trúc 5 đầu vào (tỷ lệ phần mol các khí) sinh ra các tỷ số tốc độ $q_{\text{CH}4}^T$ và $q{\text{H}2}^T$ tại 3 mức nhiệt độ chuẩn ($T_1=700\ ^\circ\text{C}, T_2=750\ ^\circ\text{C}, T_3=800\ ^\circ\text{C}$). Sau đó, hệ mờ FIS thực hiện nội suy để xác định chính xác tốc độ thực $r{\text{CH}4}$ và $r{\text{H}_2}$ tại nhiệt độ bất kỳ $T$.
  • Mô hình 3D-CFD đa vật lý (Multiphysics CFD): Triển khai trên nền tảng phần mềm COMSOL Multiphysics sử dụng phương pháp phần tử hữu hạn (FEM).
    • Truyền khối: Phương trình dòng chảy Darcy-Brinkman trong môi trường xốp và phương trình Navier-Stokes trong kênh dẫn khí; khuếch tán đa cấu tử giải bằng mô hình Stefan-Maxwell (SMM) kết hợp số hạng khuếch tán Knudsen hoặc mô hình khí bụi (Dusty-Gas Model - DGM).
    • Truyền nhiệt: Truyền nhiệt dẫn trong các pha rắn, đối lưu nhiệt trong dòng khí, nguồn nhiệt hóa học từ phản ứng reforming (thu nhiệt $Q_{\text{chem}}$) và nguồn nhiệt điện hóa/Joule (tỏa nhiệt $Q_{\text{elec}}$).
    • Ứng suất cơ-nhiệt: Tính toán phân bố ứng suất chính thứ nhất (first principal stress $\sigma_T$) trên mặt phẳng màng điện phân dựa trên gradient nhiệt độ.

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

Chương 1: Giới thiệu chung (General introduction)

Chương này tổng quan về bối cảnh năng lượng tái tạo, nguồn gốc và thành phần của khí biogas, cũng như nguyên lý vận hành của pin nhiên liệu oxit rắn (SOFC). Tác giả so sánh chi tiết ưu nhược điểm giữa cấu trúc tế bào dạng ống (tubular) và dạng phẳng (planar), cùng 5 cấu hình tế bào đơn: hỗ trợ bởi anode (anode-supported), hỗ trợ bởi chất điện phân (electrolyte-supported), hỗ trợ bởi cathode (cathode-supported), hỗ trợ bởi interconnector và hỗ trợ bởi đế xốp (porous substrate).

Cấu hình tế bào Ưu điểm Nhược điểm
Hỗ trợ bởi Anode (ASC) Anode dẫn điện tốt; hạ thấp nhiệt độ vận hành nhờ dùng màng điện phân mỏng. Nguy cơ tái oxy hóa Ni; giới hạn truyền khối qua lớp anode dày.
Hỗ trợ bởi Điện phân (ESC) Độ bền cơ học cao từ lớp điện phân đặc; ít rủi ro do tái oxy hóa anode. Điện trở thuần lớn do độ dẫn ion của điện phân thấp; yêu cầu nhiệt độ vận hành cao.
Hỗ trợ bởi Cathode Không gặp vấn đề oxy hóa kim loại; điện phân mỏng giúp giảm nhiệt độ. Độ dẫn điện cathode thấp hơn anode; giới hạn truyền khối qua cathode dày.
Hỗ trợ bởi Interconnect Các lớp điện cực/điện phân rất mỏng; cấu trúc cơ học vững chắc từ kim loại. Nguy cơ oxy hóa lớp interconnect; hạn chế trong thiết kế kênh dẫn dòng.
Đế xốp (Porous substrate) Các thành phần tế bào mỏng; có thể dùng vật liệu rẻ tiền làm giá đỡ. Tăng độ phức tạp vật liệu; bề mặt không bằng phẳng dễ gây đoản mạch.

Chương 1 cũng tính toán vùng nhiệt động học hình thành carbon bằng phần mềm HSC Chemistry 9.0 (Outotec, Finland) cho hỗn hợp biogas điển hình ($60%\ \text{CH}_4 - 40%\ \text{CO}_2$). Kết quả chỉ ra rằng tỷ lệ hơi nước trên carbon $\text{S/C} = 1,1$ là điều kiện đủ để ngăn chặn hoàn toàn hiện tượng đóng cặn carbon trong dải nhiệt độ 600–1000 °C. Tuy nhiên, việc tăng S/C làm giảm sức điện động hở mạch Nernst và tiêu tốn năng lượng hóa hơi nước, do đó cần duy trì S/C ở mức thấp nhất có thể.

Chương 2: Hành vi điện hóa của DIR-SOFC vận hành bằng khí sinh học (Electrochemical behavior of DIR-SOFCs operating with biogas)

Chương 2 trình bày các kết quả thực nghiệm trên tế bào nút ESC ($\text{Ni-10ScSZ}$, diện tích $8 \times 8\ \text{mm}^2$, chiều dày anode $\sim 38\ \mu\text{m}$) ở 800 °C:

  • Khi nạp $80\ \text{mL}\cdot\text{min}^{-1}$ hỗn hợp biogas mô phỏng ($\text{CH}_4:\text{CO}_2:\text{N}2 = 20:f{\text{CO}2}:(60-f{\text{CO}_2})$), độ chuyển hóa $\text{CH}_4$ và tốc độ sinh $\text{H}2, \text{CO}$ tăng tuyến tính theo lưu lượng $f{\text{CO}_2}$ nạp vào. Tỷ lệ mol $\text{H}_2/\text{CO}$ trong dòng khí sau phản ứng biến thiên từ 1,0 đến 1,6.
  • Kết quả đo phổ tổng trở EIS dưới điều kiện hở mạch OCV chỉ ra rằng: khi tăng $f_{\text{CO}2}$ từ $5$ lên $20\ \text{mL}\cdot\text{min}^{-1}$, điện trở truyền điện tích của anode ($R{\text{ct},a}$) giảm mạnh từ $0,38\ \Omega\cdot\text{cm}^2$ xuống $0,15\ \Omega\cdot\text{cm}^2$, trong khi điện trở truyền khối ($R_{\text{mt},a}$) duy trì ổn định ở mức thấp ($\sim 0,05\ \Omega\cdot\text{cm}^2$). Tổng trở phân cực của anode ($R_{p,a} = R_{\text{mt},a} + R_{\text{ct},a}$) giảm từ $0,43\ \Omega\cdot\text{cm}^2$ xuống $0,20\ \Omega\cdot\text{cm}^2$, tiến sát giá trị đo được khi chạy bằng $\text{H}_2$ ẩm ($0,18\ \Omega\cdot\text{cm}^2$).
  • Điều này chứng minh rằng sự gia tăng nồng độ $\text{CO}_2$ thúc đẩy phản ứng reforming sinh ra lượng lớn $\text{H}_2$ cục bộ, làm giảm đáng kể trở kháng hoạt hóa của anode.

Chương 3: Mô hình hóa quá trình đa hoàn nguyên methane trong anode gốc Ni của SOFC (Modeling of methane multiple-reforming within the Ni-based anode of an SOFC)

Chương 3 trình bày chi tiết thuật toán xây dựng mô hình MMR bán thực nghiệm kết hợp ANN và FIS:

  • Dữ liệu thực nghiệm thu được từ 6 chương trình thử nghiệm trên half-cell $\text{Ni-8YSZ}$ ở 700 °C, 750 °C và 800 °C.
  • Thuật toán trích xuất tốc độ tiêu thụ $\text{CH}4$ chuẩn ($r{\text{CH}_4}^\circ$) bằng hàm mũ công suất (power-law), sau đó đưa vào mạng $\text{ANN}N$ ($N=5$) có hàm kích hoạt sigmoid dạng đối xứng để tính các tỷ số tốc độ $q{\text{CH}4}^T$ và $q{\text{H}_2}^T$.
  • Hệ suy luận mờ FIS sử dụng các hàm thuộc loại tam giác/hình thang nội suy chính xác tốc độ phản ứng tại các giá trị nhiệt độ trung gian.
  • Mô hình được kiểm chứng độc lập trên thiết bị phản ứng reforming phẳng sử dụng chất xúc tác cấu trúc giấy (PSC) có hệ thống đo nhiệt độ phân bố bằng camera hồng ngoại. Biên dạng nhiệt độ tính toán bằng mô hình MMR-CFD hoàn toàn trùng khớp với biên dạng thực nghiệm ở tốc độ không gian GHSV $2880\ \text{h}^{-1}$.

Chương 4: Mô hình hóa và mô phỏng DIR-SOFC vận hành bằng khí sinh học (Modeling and simulation of a DIR-SOFC operating with biogas)

Chương này tích hợp mô hình MMR vào mô hình 3D-CFD hoàn chỉnh cho tế bào phẳng ASC và đề xuất chiến lược thẩm định 3 bước:

  1. Bước 1 (Non-DIR fitting): Hiệu chỉnh các thông số điện hóa và tổn hao thuần trở trên đường cong $i\text{-}V$ đo với nhiên liệu $\text{H}_2$ ($3%\ \text{H}_2\text{O}$).
  2. Bước 2 (DIR fitting): Hiệu chỉnh hệ số bề mặt hoạt tính xúc tác $k_{\text{cat}}$ trên đường cong $i\text{-}V$ đo với hỗn hợp biogas mô phỏng ($\text{CH}_4/\text{CO}_2 = 1$).
  3. Bước 3 (Arbitrary condition check): Kiểm tra độ chính xác của mô hình ở điều kiện vận hành bất kỳ mà không điều chỉnh thêm bất kỳ tham số nào.

Kết quả mô phỏng trên tế bào ASC $14 \times 14\ \text{mm}^2$ (vận hành với $40\ \text{mL}\cdot\text{min}^{-1}$ biogas $\text{CH}_4/\text{CO}_2=1$) chỉ ra:

  • Phản ứng reforming diễn ra tập trung tại $20%$ chiều dài đầu kênh nạp, khiến nhiệt độ tại đây sụt giảm sâu (tạo điểm cực tiểu nhiệt độ).
  • Khi mật độ dòng tăng lên $2\ \text{A}\cdot\text{cm}^{-2}$, nhiệt tỏa ra từ phản ứng điện hóa ($Q_{\text{elec}}$) bù trừ cho nhiệt thu của reforming ($Q_{\text{chem}}$), nâng nhiệt độ cực đại của tế bào lên $812\ ^\circ\text{C}$.
  • Ứng suất kéo cơ-nhiệt cực đại ($\sigma_{T,\max}$) xuất hiện chính xác tại vị trí có gradient nhiệt độ cao nhất gần cửa nạp.
  • So sánh đối chứng chứng minh rằng cách tiếp cận hoàn nguyên song song của Meng Ni sai lệch lớn so với thực nghiệm (dự báo quá mức tốc độ phản ứng và sai lệch đường cong $i\text{-}V$), trong khi mô hình MMR của luận án tái hiện chính xác dữ liệu thực nghiệm.

Chương 5: Các khái niệm DIR cải tiến cho SOFC vận hành bằng khí sinh học (Advanced DIR concepts for SOFCs operating with biogas)

Chương 5 khảo sát hai giải pháp kỹ thuật nhằm tối ưu hóa phân bố nhiệt và nâng cao công suất của tế bào $20 \times 50\ \text{mm}^2$ ở 800 °C với lưu lượng $200\ \text{mL}\cdot\text{min}^{-1}$ biogas ($\text{CH}_4/\text{CO}_2 = 1$):

  • Khái niệm I (Concept-I - Mặt nạ cản khí Anode / Anode gas-barrier mask): Phủ một lớp chắn khí không thấm tại vùng đầu anode nhằm kiểm soát dòng khuếch tán nhiên liệu theo chiều dọc dòng chảy.
    • Kết quả trên cả hai loại tế bào ASC-A ($d_a=950\ \mu\text{m}$) và ASC-B ($d_a=200\ \mu\text{m}$): Giảm khoảng $20%$ ứng suất nhiệt cực đại ($\sigma_{T,\max}$), loại bỏ nguy cơ nứt vỡ màng điện phân, trong khi mật độ công suất cực đại ($P_{\max}$) chỉ giảm nhẹ $\sim 8%$.
  • Khái niệm II (Concept-II - Tích hợp bộ Reforming PSC trong kênh nạp / In-cell PSC reformer): Đặt tấm xúc tác cấu trúc giấy PSC vào kênh dẫn nhiên liệu phía trên anode.
    • Với ASC-A ($d_a=950\ \mu\text{m}$): $P_{\max}$ tăng từ $1,08\ \text{W}\cdot\text{cm}^{-2}$ (cấu hình chuẩn Normal) lên $1,45\ \text{W}\cdot\text{cm}^{-2}$.
    • Với ASC-B ($d_a=200\ \mu\text{m}$): $P_{\max}$ tăng mạnh từ $0,95\ \text{W}\cdot\text{cm}^{-2}$ lên $1,45\ \text{W}\cdot\text{cm}^{-2}$.
    • Cơ chế: Đế anode mỏng của ASC-B giúp sản phẩm $\text{H}2\text{O}$ sinh ra từ phản ứng điện hóa dễ dàng thoát ngược ra kênh nạp, ngăn ngừa sự tích tụ áp suất riêng phần $p{\text{H}_2\text{O}}$ tại lớp chức năng anode (AFL) ở mật độ dòng cao, qua đó duy trì sức điện động Nernst cao dọc theo chiều dài tế bào.

Chương 6: Kết luận và định hướng nghiên cứu tiếp theo (Conclusions and outlook for future work)

Tổng kết các kết quả đạt được về mặt phương pháp luận mô hình hóa, quy trình thẩm định mô phỏng số và các khuyến nghị thiết kế tế bào đơn. Đề xuất các nội dung mở rộng cho các nghiên cứu tiếp theo trong tương lai.

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

Đóng góp về mặt lý luận và khoa học

  1. Xây dựng mô hình động học MMR phi tuyến: Luận án đã giải quyết hạn chế tồn tại nhiều năm trong việc mô phỏng DIR-SOFC chạy biogas bằng cách xây dựng mô hình bán thực nghiệm dựa trên mạng nơ-ron ANN kết hợp hệ suy luận mờ FIS. Mô hình phản ánh chính xác ảnh hưởng đồng thời và cạnh tranh của $\text{CO}_2$ và $\text{H}_2\text{O}$ lên tâm xúc tác Ni.
  2. Chiến lược thẩm định mô hình 3 bước nghiêm ngặt: Đưa ra phương pháp luận chuẩn hóa để thẩm định mô hình CFD đa vật lý cho SOFC, phân tách rõ ràng việc cân chỉnh tham số truyền điện tích/thuần trở (bước non-DIR), tham số diện tích xúc tác hoạt tính (bước DIR) và kiểm chứng độc lập ở điều kiện tùy ý.

Đóng góp về mặt kỹ thuật và thực tiễn

  1. Lượng hóa tương quan giữa độ dày Anode và cơ chế thoát nước: Chỉ ra rằng đối với tế bào anode mỏng (ASC-B, $200\ \mu\text{m}$), việc áp dụng bộ xúc tác PSC tích hợp (Concept-II) mang lại hiệu quả vượt trội trong việc giải phóng $\text{H}_2\text{O}$ khỏi AFL, hạn chế sụt giảm thế Nernst ở mật độ dòng cao.
  2. Giải pháp kiểm soát ứng suất nhiệt bằng màng cản khí: Chứng minh giải pháp phủ mặt nạ cản khí (Concept-I) giúp giảm $20%$ ứng suất kéo cơ-nhiệt cực đại trên màng điện phân, cung cấp giải pháp thiết kế kéo dài tuổi thọ cơ học cho các hệ thống DIR-SOFC.
Cấu hình khảo sát Độ dày đế Anode ($d_a$) Mật độ công suất đỉnh ($P_{\max}$) Mức giảm ứng suất nhiệt cực đại ($\Delta \sigma_{T,\max}$)
Normal ASC-A $950\ \mu\text{m}$ $1,08\ \text{W}\cdot\text{cm}^{-2}$ Cơ sở so sánh ($0%$)
Concept-I ASC-A $950\ \mu\text{m}$ $\sim 0,99\ \text{W}\cdot\text{cm}^{-2}$ (giảm $8%$) Giảm $\sim 20%$
Concept-II ASC-A $950\ \mu\text{m}$ $1,45\ \text{W}\cdot\text{cm}^{-2}$ (tăng $34%$) Tăng (do tập trung nhiệt reforming)
Normal ASC-B $200\ \mu\text{m}$ $0,95\ \text{W}\cdot\text{cm}^{-2}$ Cơ sở so sánh ($0%$)
Concept-I ASC-B $200\ \mu\text{m}$ $\sim 0,87\ \text{W}\cdot\text{cm}^{-2}$ (giảm $8%$) Giảm $\sim 20%$
Concept-II ASC-B $200\ \mu\text{m}$ $1,45\ \text{W}\cdot\text{cm}^{-2}$ (tăng $52,6%$) Tăng (cần tối ưu hóa phân bố)

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

Hạn chế được ghi nhận từ văn bản

  • Mô hình ANN/FIS được xây dựng dựa trên tập dữ liệu thực nghiệm trong phạm vi nhiệt độ $700\text{–}800\ ^\circ\text{C}$ và dải lưu lượng/nồng độ xác định; việc ngoại suy ra ngoài dải tham số này đòi hỏi phải bổ sung dữ liệu huấn luyện mới.
  • Mô hình 3D-CFD trong luận án tập trung vào trạng thái dừng (steady-state), chưa mô phỏng động học quá trình suy giảm hiệu năng do ngộ độc lưu huỳnh $\text{H}_2\text{S}$ theo thời gian thực hoặc tích tụ cặn carbon động học ở điều kiện quá tải dòng.
  • Bức xạ nhiệt giữa các bề mặt kênh dẫn khí và điện cực được đơn giản hóa trong một số miền tính toán để tiết kiệm tài nguyên vi tính.

Hướng nghiên cứu tiếp theo

  • Mở rộng mô hình đa vật lý 3D-CFD cho quy mô bó pin (stack) nhiều tầng tế bào và tích hợp các kênh phân phối dòng khí phức tạp.
  • Phát triển mô hình động học quá độ (transient model) để phân tích ứng suất nhiệt và cơ chế phản ứng trong các chu trình khởi động/tắt máy (thermal cycling) và thay đổi tải đột ngột.
  • Nghiên cứu kết hợp cấu trúc xúc tác PSC phân cấp chức năng (functionally-graded PSC) với màng cản khí để vừa triệt tiêu ứng suất nhiệt vừa đạt mật độ công suất cực đại.

Giá trị tham khảo

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

  • Nghiên cứu sinh và học viên sau đại học: Thuộc các chuyên ngành Kỹ thuật Hóa học, Kỹ thuật Năng lượng Hydrogen, Khoa học Vật liệu và Kỹ thuật Nhiệt; có thể tham khảo phương pháp kết hợp trí tuệ nhân tạo (ANN/FIS) vào mô hình hóa động học phản ứng hóa học phức tạp.
  • Giảng viên và nhà nghiên cứu pin nhiên liệu: Cung cấp nguồn tài liệu chi tiết về phương pháp đo đạc phổ tổng trở điện hóa EIS, phân tích cơ chế phản ứng trên anode cermet Ni-YSZ/Ni-ScSZ và quy trình thẩm định mô hình mô phỏng số 3D.
  • Kỹ sư thiết kế và phát triển hệ thống SOFC: Sử dụng các số liệu về suy giảm ứng suất nhiệt, thông số cấu trúc xốp và các thiết kế cải tiến (Concept-I, Concept-II) để tối ưu hóa hình học tế bào và kênh dẫn khí cho các tổ máy phát điện phân tán chạy biogas quy mô trang trại hoặc nhà máy xử lý rác.

Các nội dung có giá trị tham khảo cao trong tài liệu

  • Bảng tổng hợp chi tiết các phương trình bảo toàn khối lượng, động lượng Darcy-Brinkman, khuếch tán đa cấu tử Stefan-Maxwell và phương trình động học điện hóa Butler-Volmer (Chương 4).
  • Bộ tham số thực nghiệm về độ xốp, độ xoắn, độ thấm và diện tích hoạt tính riêng của các lớp chức năng anode/cathode (Bảng 4.2 – 4.4 trong văn bản).
  • Quy trình đo đạc động học reforming trên đĩa half-cell và sơ đồ ghép nối thuật toán ANN-FIS (Chương 3 và Phụ lục B, C).

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

1. Tại sao cách tiếp cận hoàn nguyên song song truyền thống (Parallel-reforming approach) lại gây sai lệch khi mô phỏng DIR-SOFC chạy biogas?

Cách tiếp cận truyền thống tính toán tốc độ phản ứng tổng bằng cách cộng độc lập tốc độ hoàn nguyên khô ($\text{CH}_4\text{-}\text{CO}_2$) và hoàn nguyên ướt ($\text{CH}_4\text{-}\text{H}_2\text{O}$). Tuy nhiên, trong thực tế vận hành DIR với biogas, $\text{CO}_2$ và $\text{H}_2\text{O}$ cùng hiện diện và cạnh tranh trực tiếp trên cùng các tâm xúc tác kim loại Ni. Cách cộng tuyến tính bỏ qua hiệu ứng tương tác đồng thời này, dẫn đến việc ước tính sai tốc độ tiêu thụ $\text{CH}_4$ cục bộ, làm sai lệch phân bố nhiệt độ và không thể dự báo chính xác đường đặc trưng dòng-áp ($i\text{-}V$) của tế bào pin.

2. Mô hình MMR dựa trên ANN/FIS hoạt động theo nguyên lý nào?

Mô hình sử dụng mạng nơ-ron nhân tạo ANN 5 đầu vào (tương ứng với nồng độ phần mol của $\text{CH}_4, \text{CO}_2, \text{H}_2\text{O}, \text{H}2, \text{CO}$) để học các tỷ số tốc độ phi tuyến $q{\text{CH}4}^T$ và $q{\text{H}_2}^T$ từ dữ liệu thực nghiệm tại 3 mức nhiệt độ cơ sở (700 °C, 750 °C, 800 °C). Khi cần tính toán ở một nhiệt độ bất kỳ trong dải khảo sát, Hệ suy luận mờ (FIS) sẽ thực hiện nội suy mờ các giá trị tốc độ chuẩn để đưa ra tốc độ tiêu thụ và sinh khí thực tế của từng cấu tử.

3. Khái niệm cải tiến Concept-I (Anode gas-barrier mask) mang lại lợi ích gì cho độ bền cơ học của tế bào SOFC?

Concept-I sử dụng một lớp mặt nạ cản khí đặt tại đầu nạp của anode nhằm hạn chế bớt lưu lượng khí nhiên liệu thâm nhập quá nhanh vào thể tích xốp ở cửa vào. Nhờ đó, phản ứng hoàn nguyên methane thu nhiệt mạnh bị kìm hãm bớt ở vùng đầu kênh, giúp phân bố nhiệt độ dọc tế bào trở nên đồng đều hơn. Mô phỏng số xác nhận Concept-I giúp giảm khoảng $20%$ ứng suất kéo cơ-nhiệt cực đại ($\sigma_{T,\max}$) trên màng điện phân, giảm thiểu nguy cơ nứt vỡ điện phân trong khi công suất đỉnh chỉ giảm khoảng $8%$.

4. Tại sao bộ reforming tích hợp PSC (Concept-II) lại đạt hiệu quả vượt trội hơn trên tế bào có đế anode mỏng (ASC-B)?

Trên tế bào ASC-B ($d_a = 200\ \mu\text{m}$), chiều dày lớp anode xốp mỏng hơn nhiều so với ASC-A ($d_a = 950\ \mu\text{m}$), tạo điều kiện cho hơi nước ($\text{H}_2\text{O}$) sinh ra từ phản ứng oxy hóa điện hóa tại lớp chức năng AFL nhanh chóng thoát ra ngoài kênh dẫn nhiên liệu. Việc này ngăn ngừa sự gia tăng áp suất riêng phần của $\text{H}_2\text{O}$ tại AFL khi tế bào vận hành ở mật độ dòng điện cao, duy trì sức điện động Nernst ở mức cao và giúp công suất đỉnh đạt tới $1,45\ \text{W}\cdot\text{cm}^{-2}$ (tăng $52,6%$ so với cấu hình chuẩn).

5. Quy trình thẩm định mô hình ba bước (Three-step model validation) gồm những giai đoạn nào?

Quy trình gồm:

  • Bước 1: Tinh chỉnh các tham số điện hóa và điện trở thuần dựa trên dữ liệu đường cong $i\text{-}V$ đo ở chế độ không hoàn nguyên (chạy bằng $\text{H}_2$ ẩm).
  • Bước 2: Tinh chỉnh hệ số hiệu chỉnh diện tích bề mặt xúc tác hoạt tính ($k_{\text{cat}}$) dựa trên dữ liệu thực nghiệm ở chế độ DIR (chạy bằng hỗn hợp $\text{CH}_4\text{-}\text{CO}_2$).
  • Bước 3: Kiểm tra tính hợp lệ của toàn bộ mô hình đã thiết lập bằng cách đối chiếu với kết quả vận hành thực tế ở một điều kiện hoạt động tùy ý mà không thay đổi bất kỳ tham số nào.

Kết luận

Luận án đã phát triển thành công mô hình động học đa hoàn nguyên methane (MMR) mới dựa trên mạng nơ-ron nhân tạo và hệ suy luận mờ (ANN/FIS), giải quyết triệt để sự thiếu chính xác của các cách tiếp cận mô hình hóa hoàn nguyên song song truyền thống. Bằng việc tích hợp mô hình MMR vào công cụ mô phỏng 3D-CFD đa vật lý kết hợp quy trình thẩm định thực nghiệm ba bước chặt chẽ, nghiên cứu đã làm sáng tỏ mối quan hệ tương hỗ giữa động học phản ứng hóa học, quá trình truyền chất, phân bố gradient nhiệt độ và ứng suất cơ-nhiệt trong pin SOFC nạp trực tiếp biogas. Hai khái niệm thiết kế cải tiến (mặt nạ cản khí Concept-I và bộ reforming PSC Concept-II) cùng các chỉ dẫn kỹ thuật về độ dày chất nền anode cung cấp cơ sở khoa học và công cụ tính toán định lượng quan trọng cho việc thiết kế các hệ thống pin nhiên liệu oxit rắn vận hành bằng biogas hiệu quả và bền bỉ.