Tổng quan về luận án

Sự cạn kiệt của nguồn năng lượng hóa thạch cùng các tác động tiêu cực từ biến đổi khí hậu toàn cầu do phát thải khí nhà kính ($\text{CO}_2, \text{SO}_x, \text{NO}_x$) đặt ra yêu cầu cấp bách về việc phát triển các nguồn năng lượng tái tạo. Trong số các dạng năng lượng sạch, sinh khối (biomass) là nguồn năng lượng duy nhất chứa carbon có khả năng chuyển hóa trực tiếp thành nhiên liệu lỏng, rắn và khí. Theo Tổ chức Năng lượng Quốc tế (IEA), sinh khối chiếm khoảng 63% tổng nguồn năng lượng tái tạo và đáp ứng khoảng 10% nhu cầu năng lượng toàn cầu. Tại Việt Nam, trữ lượng phụ phẩm nông lâm nghiệp ước tính đạt 45 triệu tấn/năm (tương đương 15 triệu tấn dầu quy đổi), chủ yếu là bột gỗ, mùn cưa và bã mía.

Công nghệ nhiệt phân nhanh (fast pyrolysis) nổi lên như một giải pháp tiên phong so với quá trình đốt cháy truyền thống và khí hóa (gasification). Quá trình nhiệt phân nhanh diễn ra trong điều kiện không có oxy, tốc độ gia nhiệt cực cao ($>1000^\circ\text{C/s}$), nhiệt độ phản ứng trung bình ($450 - 550^\circ\text{C}$) và thời gian lưu của pha hơi cực ngắn ($\tau < 2\text{ s}$), cho phép tối ưu hóa hiệu suất thu hồi nhiên liệu lỏng sinh học (bio-oil) lên đến $60 - 70\text{ wt}%$. Dầu sinh học có mật độ năng lượng cao ($17 - 19\text{ MJ/kg}$), thuận tiện cho việc lưu trữ, vận chuyển và thay thế nhiên liệu đốt lò công nghiệp hoặc nâng cấp thành nhiên liệu động cơ.

Mặc dù công nghệ nhiệt phân nhanh trên thế giới đã đạt được những bước tiến đáng kể (Bridgwater, 2012; Mohan et al., 2006), việc ứng dụng tại các quốc gia đang phát triển vẫn đối mặt với research gap sâu sắc:

  • Thiếu hụt mô hình toán học giải tích mô tả chính xác trường nhiệt độ phi đối xứng và động học truyền nhiệt nội hạt đối với các loại sinh khối nhiệt đới đặc thù.
  • Thiếu các dữ liệu động học thực nghiệm xác định năng lượng hoạt hóa ($E_a$) và hằng số tốc độ phản ứng ($k_i$) của bột gỗ keo và bã mía Việt Nam trong môi trường lò tầng sôi bọt (bubbling fluidized bed).
  • Sự thiếu đồng bộ trong tính toán thiết kế tích hợp hệ thống từ lò phản ứng tầng sôi, thiết bị tách hạt rắn cyclone đến hệ thống ngưng tụ thu gom sương dầu (aerosol).

Luận án của NCS. Phạm Duy Vũ (2017) với đề tài "Nghiên cứu quá trình nhiệt phân biomass sản xuất nhiên liệu sinh học" thuộc chuyên ngành Kỹ thuật nhiệt (Mã số: 62.52.01.15) tại Đại học Đà Nẵng, dưới sự hướng dẫn khoa học của PGS. Hoàng Dương Hùng và PGS. Trần Văn Vang, đã giải quyết toàn diện các khoảng trống nghiên cứu này thông qua hệ thống câu hỏi và giả thuyết khoa học:

  • Research Question 1 (RQ1): Quy luật phân bố trường nhiệt độ phi dừng bên trong hạt sinh khối dạng hình cầu phụ thuộc như thế nào vào chuẩn số Biot ($Bi$), chuẩn số Fourier ($Fo$) và hệ số trao đổi nhiệt phức hợp ($\alpha$) trong lò tầng sôi?
  • Research Question 2 (RQ2): Cấu trúc khí động học, sự phân bố nồng độ thể tích của cát truyền nhiệt và hạt sinh khối trong lò tầng sôi biến thiên ra sao theo vận tốc khí nitơ ($N_2$)?
  • Research Question 3 (RQ3): Các thông số công nghệ (nhiệt độ phản ứng, kích thước hạt, lưu lượng khí mang) ảnh hưởng như thế nào đến hiệu suất thu hồi và thành phần hóa học của dầu sinh học?
  • Research Question 4 (RQ4): Giá trị các thông số động học phản ứng phân hủy cạnh tranh bậc một ($E_{a1}, E_{a2}, E_{a3}$ và $A_1, A_2, A_3$) của bột gỗ và bã mía trong lò tầng sôi là bao nhiêu?

Hệ thống giả thuyết khoa học được xác lập:

  • Hypothesis 1 (H1): Khi bán kính hạt sinh khối $R \le 0.25\text{ mm}$ ($d \le 0.5\text{ mm}$), gradient nhiệt độ nội hạt được giảm thiểu tối đa, thời gian hoàn thành nhiệt phân $\tau_p \le 1.2\text{ s}$, tối ưu hóa quá trình tạo pha hơi và ức chế phản ứng tạo cốc.
  • Hypothesis 2 (H2): Nhiệt độ phản ứng trong khoảng $480 - 500^\circ\text{C}$ với thời gian lưu pha hơi $<2\text{ s}$ sẽ triệt tiêu phản ứng bẻ gãy thứ cấp (secondary cracking), cho hiệu suất thu hồi dầu sinh học cực đại vượt trên $60\text{ wt}%$.
  • Hypothesis 3 (H3): Mô hình động học 3 phản ứng song song cạnh tranh bậc một phản ánh chính xác cơ chế phân hủy nhiệt của sinh khối trong môi trường tầng sôi thực tế.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp của Định luật truyền nhiệt Fourier phi dừng vi mô, Lý thuyết thủy động lực học môi trường xốp tầng sôi (Eulerian-Eulerian Granular Flow Theory) và Mô hình động học phân hủy nhiệt Shafizadeh-Bradbury. Luận án đã thiết kế, chế tạo và vận hành thành công hệ thống nhiệt phân nhanh lò tầng sôi năng suất $500\text{ g/h}$, xác lập hiệu suất thu hồi dầu sinh học đạt $64.8%$ đối với bột gỗ và $61.2%$ đối với bã mía, đóng góp cơ sở dữ liệu thực nghiệm chuẩn xác cho ngành kỹ thuật nhiệt và năng lượng tái tạo.

Literature Review và Positioning

Nghiên cứu về nhiệt phân nhanh sinh khối trên thế giới phát triển qua ba dòng tư tưởng học thuật chính:

  1. Dòng nghiên cứu công nghệ và tối ưu hóa thiết bị phản ứng: Tiên phong bởi Scott và Piskorz (1984) tại Đại học Waterloo với mô hình tầng sôi nông, sau đó được Bridgwater (2012) hệ thống hóa tại Đại học Aston (Vương quốc Anh). Dòng nghiên cứu này tập trung chứng minh lò tầng sôi bọt và tầng sôi tuần hoàn là giải pháp tối ưu nhờ khả năng truyền nhiệt đối lưu - dẫn nhiệt vượt trội từ môi chất cát sang hạt sinh khối với hệ số trao đổi nhiệt $\alpha$ đạt $300 - 600\text{ W/m}^2\text{K}$.
  2. Dòng nghiên cứu mô hình hóa truyền nhiệt và động học phân hủy: Dẫn đầu bởi Di Blasi (2008), Bellais et al. (2003) và Kersten et al. (2005). Các tác giả tập trung vào cơ chế giới hạn truyền nhiệt nội hạt (internal heat transfer limitation) và sự tương tác giữa số Biot ($Bi = \alpha R / \lambda$) với số Pyrolysis ($Py = \lambda / (\rho C_p k R^2)$), chỉ ra rằng kích thước hạt chi phối trực tiếp đến phân bố sản phẩm lỏng - rắn - khí.
  3. Dòng nghiên cứu hóa học dầu sinh học và xúc tác nâng cấp: Mohan, Pittman và Steele (2006), Zhang et al. (2007), Lappas et al. (2002) đi sâu vào phân tích cấu trúc oligomer phức tạp của bio-oil có nguồn gốc từ sự phân hủy cellulose, hemicellulose và lignin, đồng thời phát triển các giải pháp ngưng tụ phân đoạn và xúc tác zeolite (H-ZSM-5).

Trong y văn tồn tại hai cuộc tranh luận học thuật then chốt:

  • Tranh luận 1: Khống chế động học (Kinetic control) đối đầu Khống chế truyền nhiệt (Heat transfer control). Trường phái Di Blasi (2008) cho rằng trong hầu hết các điều kiện thực tế, hạt sinh khối luôn chịu sự cản trở nhiệt nội tại lớn ($Bi > 1$), tạo ra sự chênh lệch pha phản ứng giữa tâm và vỏ hạt. Ngược lại, trường phái Bridgwater (2012) giả định rằng nếu nghiền hạt đủ nhỏ ($d < 0.5\text{ mm}$), có thể coi hạt đẳng nhiệt ($Bi \to 0$) và quá trình hoàn toàn do động học hóa học kiểm soát. Luận án của Phạm Duy Vũ đã giải quyết mâu thuẫn này bằng cách thiết lập nghiệm giải tích phương trình truyền nhiệt vi phân, xác định chính xác bán kính giới hạn $R = 0.25\text{ mm}$ là điểm chuyển tiếp điều kiện biên mà tại đó gradient nhiệt độ nội hạt không còn gây ảnh hưởng tiêu cực đến hiệu suất bốc hơi.
  • Tranh luận 2: Cơ chế bẻ gãy sơ cấp (Primary decomposition) đối đầu Bẻ gãy thứ cấp (Secondary cracking). Shafizadeh & Chin (1977) cùng Bradbury et al. (1979) nhận định dầu sinh học là sản phẩm trực tiếp của quá trình bẻ gãy liên kết glycosidic và ete sơ cấp. Tuy nhiên, Liden et al. (1988) nhấn mạnh rằng nếu thời gian lưu của pha hơi trong vùng nhiệt độ cao vượt quá $2\text{ s}$, các phân tử hơi hữu cơ sẽ trải qua phản ứng đồng thể thứ cấp tạo ra khí không ngưng và cốc thứ cấp. Luận án đã định vị và chứng minh bằng thực nghiệm rằng việc duy trì lưu lượng khí nitơ ở mức $1.2 - 1.4\text{ m}^3\text{/h}$ đảm bảo thời gian lưu $\tau \approx 1.0 - 1.5\text{ s}$, bảo tồn tối đa các hợp chất thơm và đường khan.

So sánh với các nghiên cứu quốc tế tiêu biểu:

  • Nghiên cứu của Garcia-Perez et al. (2008) trên gỗ mềm Bắc Mỹ bằng lò chân không chỉ đạt hiệu suất dầu $55\text{ wt}%$, do tốc độ gia nhiệt thấp thúc đẩy tạo cốc.
  • Công trình của Heo et al. (2010) trên bã mía bằng lò tầng sôi tầng cố định tại Hàn Quốc đạt hiệu suất dầu $52.4\text{ wt}%$, thấp hơn mức $61.2\text{ wt}%$ của luận án do hiệu quả ngưng tụ aerosol chưa tối ưu.

Luận án đã định vị xuất sắc đóng góp học thuật: Lần đầu tiên tại Việt Nam, một nghiên cứu đã khép kín từ mô hình hóa truyền nhiệt giải tích vi mô, mô phỏng số khí động học CFD đa pha 3D trung mô, đến chế tạo hệ thống thực nghiệm tầng sôi hoàn chỉnh và xác lập bộ hằng số động học $E_a, A$ cho hai nguồn sinh khối chủ lực nhiệt đới.

Đó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 phương trình dẫn nhiệt phi dừng Fourier sang bài toán biên loại 3 có điều kiện biên phức hợp đối với hạt sinh khối hình cầu:

$$\frac{\partial t(r, \tau)}{\partial \tau} = a \left( \frac{\partial^2 t(r, \tau)}{\partial r^2} + \frac{2}{r} \frac{\partial t(r, \tau)}{\partial r} \right)$$

với các điều kiện đơn trị:

  • Điều kiện ban đầu: $t(r, 0) = t_0 = \text{const}$ khi $\tau = 0$.
  • Điều kiện biên đối xứng tại tâm: $\left. \frac{\partial t}{\partial r} \right|_{r=0} = 0$.
  • Điều kiện biên loại 3 trên bề mặt hạt ($r = R$):

$$-\lambda \left. \frac{\partial t(r, \tau)}{\partial r} \right|_{r=R} = \alpha \left[ t(R, \tau) - t_f \right]$$

Bằng phương pháp tách biến Fourier, luận án đã thiết lập công thức nghiệm không thứ nguyên giải tích mô tả trường nhiệt độ trong hạt:

$$\theta(\xi, Fo) = \frac{t(r, \tau) - t_f}{t_0 - t_f} = \sum_{n=1}^{\infty} A_n \frac{\sin(\mu_n \xi)}{\mu_n \xi} \exp(-\mu_n^2 Fo)$$

trong đó $\xi = r/R$ là bán kính tương đối, $Fo = a\tau/R^2$ là chuẩn số Fourier, $Bi = \alpha R/\lambda$ là chuẩn số Biot, và $\mu_n$ là các nghiệm dương của phương trình đặc trưng siêu việt:

$$\tan \mu_n = \frac{\mu_n}{1 - Bi}$$

với biên độ dao động sóng nhiệt $A_n$:

$$A_n = \frac{2 (\sin \mu_n - \mu_n \cos \mu_n)}{\mu_n - \sin \mu_n \cos \mu_n}$$

Từ nghiệm giải tích này, luận án rút ra công thức xác định thời gian nhiệt phân hoàn toàn $\tau_p$ của hạt sinh khối bán kính $R$:

$$\tau_p = \frac{R^2}{a \mu_1^2} \ln \left[ \frac{A_1}{\theta_p} \frac{\sin \mu_1}{\mu_1} \right]$$

Mô hình lý thuyết này tạo ra bước chuyển dịch quan niệm (paradigm shift): Chứng minh rằng việc gia tăng hệ số trao đổi nhiệt $\alpha$ chỉ có ý nghĩa rút ngắn thời gian phản ứng khi $Bi < 1$; khi $Bi > 5$, tốc độ gia nhiệt bị chặn hoàn toàn bởi độ dẫn nhiệt nội tại $\lambda$ của sinh khối ($\lambda_{\text{gỗ}} \approx 0.12 - 0.15\text{ W/mK}$, $\lambda_{\text{bã mía}} \approx 0.08 - 0.11\text{ W/mK}$).

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

Khung phân tích của luận án là sự tích hợp đồng bộ giữa 3 cấu phần lý thuyết:

  1. Lý thuyết truyền nhiệt vi mô (Micro-scale heat transfer): Xác định động thái tăng nhiệt của từng lớp vỏ sinh khối từ bề mặt vào tâm, chỉ ra rằng nhiệt độ bắt đầu phân hủy nhiệt của hemicellulose và cellulose diễn ra tại $t_p \ge 250^\circ\text{C}$ và hoàn tất ở $t_f = 500^\circ\text{C}$.
  2. Lý thuyết thủy động lực học đa pha tầng sôi (Meso-scale hydrodynamics): Mô hình hóa tương tác giữa pha khí ($N_2$), pha rắn trơ (cát thạch anh đóng vai trò chất mang nhiệt, kích thước $0.2 - 0.3\text{ mm}$, khối lượng riêng $\rho = 2650\text{ kg/m}^3$) và pha hạt sinh khối ($\rho \approx 400 - 600\text{ kg/m}^3$). Sử dụng mô hình Euler-Euler kết hợp Thuyết động học chất hạt (Kinetic Theory of Granular Flow - KTGF) để mô phỏng trường vận tốc và trở lực lớp sôi.
  3. Mô hình động học phản ứng phân hủy cạnh tranh bậc một (Macro-scale kinetics): Phân chia sinh khối ($S$) thành 3 nhánh sản phẩm song song độc lập:

$$\begin{aligned} S &\xrightarrow{k_1} \text{Dầu sinh học (Bio-oil)} \ S &\xrightarrow{k_2} \text{Khí không ngưng (Gas)} \ S &\xrightarrow{k_3} \text{Cốc sinh học (Biochar)} \end{aligned}$$

Hằng số tốc độ phản ứng $k_i$ tuân theo phương trình Arrhenius:

$$k_i = A_i \exp\left( -\frac{E_{a,i}}{R_g T} \right) \quad (i = 1, 2, 3)$$

Khung phân tích xác lập rõ các điều kiện biên: Áp suất phản ứng khí quyển $P = 1\text{ atm}$, môi trường trơ $N_2$ hoàn toàn không oxy, hình học hạt được quy đổi về hình cầu đẳng hướng tương đương, độ ẩm ban đầu của sinh khối $W \le 10%$.

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 lý nghiên cứu thực chứng và kỹ thuật thực nghiệm phân tích (Positivism & Empirical-Analytical Engineering Approach). Thiết kế nghiên cứu đa tầng bậc (Multi-level scale design) được tổ chức chặt chẽ:

  • Cấp độ vi mô: Mô phỏng số trên MATLAB để giải hệ phương trình vi phân đạo hàm riêng của trường nhiệt độ $t(r, \tau)$ và khảo sát ảnh hưởng của số $Bi$, $Fo$ đến thời gian phản ứng.
  • Cấp độ trung mô: Sử dụng phần mềm ANSYS Fluent để mô phỏng khí động học 3D trong lò phản ứng hình trụ đường kính trong $D = 50\text{ mm}$, chiều cao $H = 600\text{ mm}$, xác định vận tốc tạo lớp sôi tối ưu $\omega_s = 0.15 - 0.25\text{ m/s}$ và vị trí cấp liệu $h_{\text{cl}} = 150\text{ mm}$ phía trên ghi phân phối khí để tối đa hóa sự tiếp xúc pha.
  • Cấp độ vĩ mô: Chế tạo mô hình thực nghiệm lò tầng sôi bọt gia nhiệt điện trở ngoài công suất $3.5\text{ kW}$, tích hợp trục vít cấp liệu kín khí điều khiển bằng biến tần, cyclone thu hồi cốc hiệu suất cao và hệ thống ngưng tụ 2 cấp (bình ngưng ống chùm giải nhiệt nước và bình ngưng thu hồi aerosol).

Quy trình nghiên cứu rigorous

Quy trình chuẩn bị mẫu và thực nghiệm tuân thủ các tiêu chuẩn quốc tế nghiêm ngặt:

  • Chuẩn bị mẫu sinh khối: Bột gỗ keo lai (Acacia mangium) lấy từ xưởng chế biến gỗ tại Đà Nẵng và bã mía (Saccharum officinarum) lấy từ Nhà máy đường Quảng Ngãi. Mẫu được phơi khô tự nhiên, sấy ở $105^\circ\text{C}$ trong tủ sấy chân không đến khi độ ẩm đạt $W < 8%$, sau đó nghiền và sàng rây thành 4 dải kích thước hạt chính xác: $<0.25\text{ mm}$, $0.25 - 0.5\text{ mm}$, $0.5 - 1.0\text{ mm}$, và $1.0 - 2.0\text{ mm}$.
  • Hệ thống phân tích hóa lý và nhiệt:
    • Phân tích nhiệt khối lượng TGA/DTG trên thiết bị Netzsch STA 409 PC/PG trong dải nhiệt độ từ $30^\circ\text{C}$ đến $800^\circ\text{C}$ với tốc độ gia nhiệt $10, 20, 40^\circ\text{C/min}$ dưới dòng khí $N_2$ tinh khiết $99.999%$ ($60\text{ ml/min}$).
    • Phân tích thành phần nguyên tố (C, H, O, N, S) bằng máy Flash 2000 CHNS/O Analyzer (Thermo Scientific) theo tiêu chuẩn ASTM D5291.
    • Phân tích thành phần hóa học của dầu sinh học bằng hệ thống Sắc ký khí khối phổ GC-MS (Agilent 7890A GC kết hợp 5975C MSD, cột mao quản HP-5MS kích thước $30\text{ m} \times 0.25\text{ mm} \times 0.25\text{ }\mu\text{m}$).
    • Xác định nhiệt trị cấp cao (HHV) bằng bom nhiệt lượng tự động theo tiêu chuẩn ASTM D240.
    • Xác định độ nhớt động học bằng nhớt kế Cannon-Fenske ở $40^\circ\text{C}$ theo tiêu chuẩn ASTM D445.
    • Xác định hàm lượng nước trong bio-oil bằng phương pháp chuẩn độ điện lượng Karl Fischer theo tiêu chuẩn ASTM E203.
    • Xác định độ pH bằng máy đo điện cực pH Metrohm 827.
  • Quy trình kiểm soát sai số (Triangulation & Reliability): Mỗi chế độ công nghệ được lặp lại độc lập ít nhất 3 lần. Sai số cân bằng khối lượng (mass balance closure) toàn hệ thống đạt mức cao $95.5 - 98.2%$. Các kết quả đo đạc thống kê có độ lệch chuẩn $SD < 2.5%$.

Data và phân tích

Dữ liệu đặc tính hóa lý cơ bản của hai loại nguyên liệu sinh khối được xác định qua phân tích thực nghiệm:

  • Bột gỗ: Độ ẩm $W = 7.2%$; Chất bốc $V^k = 78.5%$; Carbon cố định $C^k = 13.8%$; Tro $A^k = 0.5%$. Thành phần nguyên tố: $\text{C} = 48.2%$, $\text{H} = 6.1%$, $\text{O} = 45.1%$, $\text{N} = 0.5%$, $\text{S} = 0.1%$.
  • Bã mía: Độ ẩm $W = 8.1%$; Chất bốc $V^k = 76.2%$; Carbon cố định $C^k = 13.2%$; Tro $A^k = 2.5%$. Thành phần nguyên tố: $\text{C} = 45.8%$, $\text{H} = 5.9%$, $\text{O} = 47.6%$, $\text{N} = 0.6%$, $\text{S} = 0.1%$.

Dữ liệu phân tích nhiệt TGA chỉ ra quá trình phân hủy nhiệt diễn ra qua 3 giai đoạn:

  1. Giai đoạn 1 ($30 - 150^\circ\text{C}$): Bay hơi ẩm tự do và ẩm liên kết.
  2. Giai đoạn 2 ($200 - 400^\circ\text{C}$): Vùng phản ứng chính, hemicellulose phân hủy mạnh ở $220 - 315^\circ\text{C}$, cellulose phân hủy ở $315 - 390^\circ\text{C}$ với tốc độ mất khối lượng cực đại tại $345^\circ\text{C}$ đối với bột gỗ và $335^\circ\text{C}$ đối với bã mía.
  3. Giai đoạn 3 ($400 - 600^\circ\text{C}$): Lignin phân hủy chậm và kéo dài, hình thành cốc rắn.

Dữ liệu thực nghiệm động học được xử lý hồi quy phi tuyến trên phần mềm để xác định hệ số góc và điểm cắt của đường biểu diễn $\ln(k_i)$ theo $1/T$, từ đó trích xuất các giá trị $E_{a,i}$ và $A_i$ với hệ số tương quan $R^2 > 0.985$.

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

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

Luận án đã đạt được 5 phát hiện then chốt mang tính đột phá khoa học:

  1. Xác lập nhiệt độ phản ứng tối ưu $T_{\text{opt}} = 480 - 500^\circ\text{C}$ cho hiệu suất thu hồi dầu sinh học cực đại:

    • Đối với bột gỗ: Hiệu suất thu hồi dầu sinh học đạt cực đại $64.8\text{ wt}%$ tại $500^\circ\text{C}$ (trong đó tỷ lệ cốc là $16.5\text{ wt}%$, khí không ngưng là $18.7\text{ wt}%$). Khi nhiệt độ $<400^\circ\text{C}$, hiệu suất dầu chỉ đạt $42.3\text{ wt}%$ do nhiệt lượng không đủ phá vỡ cấu trúc cellulose và lignin (cốc chiếm $38.2\text{ wt}%$). Khi nhiệt độ tăng lên $>550^\circ\text{C}$, phản ứng bẻ gãy thứ cấp pha hơi chiếm ưu thế, hiệu suất dầu giảm xuống $51.2\text{ wt}%$ trong khi khí không ngưng tăng vọt lên $33.5\text{ wt}%$.
    • Đối với bã mía: Hiệu suất dầu sinh học đạt đỉnh $61.2\text{ wt}%$ tại $480^\circ\text{C}$ (cốc chiếm $17.8\text{ wt}%$, khí chiếm $21.0\text{ wt}%$). Bã mía đạt đỉnh nhiệt phân sớm hơn bột gỗ $20^\circ\text{C}$ do hàm lượng hemicellulose kém bền nhiệt cao hơn.

  2. Chứng minh quy luật kích thước hạt tới hạn $d \le 0.5\text{ mm}$ ($R \le 0.25\text{ mm}$):

    • Kết quả mô phỏng và thực nghiệm khẳng định với hạt kích thước $d < 0.25\text{ mm}$ và $0.25 - 0.5\text{ mm}$, thời gian gia nhiệt tâm hạt đạt $500^\circ\text{C}$ chỉ mất tương ứng $0.4\text{ s}$ và $0.9\text{ s}$.
    • Khi kích thước hạt tăng lên $1.0 - 2.0\text{ mm}$, hiệu ứng cản trở truyền nhiệt nội hạt làm thời gian nhiệt phân kéo dài lên $>3.5\text{ s}$, dẫn đến hiện tượng cốc hóa bề mặt (char crusting), làm suy giảm hiệu suất thu hồi dầu sinh học xuống dưới $48.5\text{ wt}%$ và tăng lượng cốc lên $26.8\text{ wt}%$.

  3. Xác lập bộ thông số động học phản ứng nhiệt phân nhanh chính xác trong lò tầng sôi:

    • Bột gỗ:
      • Nhánh tạo Bio-oil: $E_{a1} = 72.45\text{ kJ/mol}$, $A_1 = 2.45 \times 10^4\text{ s}^{-1}$ ($R^2 = 0.991$).
      • Nhánh tạo Khí không ngưng: $E_{a2} = 88.62\text{ kJ/mol}$, $A_2 = 8.12 \times 10^5\text{ s}^{-1}$ ($R^2 = 0.988$).
      • Nhánh tạo Cốc: $E_{a3} = 61.20\text{ kJ/mol}$, $A_3 = 1.18 \times 10^3\text{ s}^{-1}$ ($R^2 = 0.986$).
    • Bã mía:
      • Nhánh tạo Bio-oil: $E_{a1} = 68.12\text{ kJ/mol}$, $A_1 = 1.87 \times 10^4\text{ s}^{-1}$ ($R^2 = 0.993$).
      • Nhánh tạo Khí không ngưng: $E_{a2} = 83.54\text{ kJ/mol}$, $A_2 = 5.40 \times 10^5\text{ s}^{-1}$ ($R^2 = 0.989$).
      • Nhánh tạo Cốc: $E_{a3} = 57.85\text{ kJ/mol}$, $A_3 = 0.95 \times 10^3\text{ s}^{-1}$ ($R^2 = 0.987$).
    • Kết quả chỉ ra rằng năng lượng hoạt hóa nhánh tạo dầu sinh học luôn nhỏ hơn nhánh tạo khí ($E_{a1} < E_{a2}$), chứng minh về mặt nhiệt động học rằng ở vùng nhiệt độ trung bình ($480 - 500^\circ\text{C}$), tốc độ phản ứng tạo lỏng luôn vượt trội so với phản ứng tạo khí.

  4. Định lượng đặc tính hóa lý và thành phần phân tử của Bio-oil:

    • Dầu sinh học thô từ bột gỗ có khối lượng riêng $\rho = 1.18\text{ g/cm}^3$, độ nhớt động học $\nu = 22.5\text{ cSt}$ ($40^\circ\text{C}$), giá trị $\text{pH} = 2.85$, hàm lượng nước $W = 19.5\text{ wt}%$, nhiệt trị cấp cao $\text{HHV} = 18.6\text{ MJ/kg}$ (đạt gần $45%$ so với dầu FO khoáng $\sim 41.8\text{ MJ/kg}$).
    • Phân tích GC-MS nhận diện được 48 hợp chất hữu cơ chính trong dầu sinh học bột gỗ: Dẫn xuất Phenol (Guaiacol, 4-methylguaiacol, Syringol, Eugenol) chiếm $41.2%$ diện tích pic; Axit hữu cơ (chủ yếu là Axit axetic) chiếm $14.5%$; Hợp chất Furan (Furfural, 5-HMF) chiếm $9.8%$; Đường khan (Levoglucosan) chiếm $11.6%$; nhóm Hydrocacbon và các hợp chất khác chiếm $22.9%$. Sự hiện diện của nồng độ cao các hợp chất phenolic mở ra tiềm năng ứng dụng sản xuất keo dán phenolic sinh học (phenolic resins) và hóa chất tinh khiết.

  5. Giải mã thành phần khí không ngưng và cơ chế tái sử dụng năng lượng:

    • Khí không ngưng sinh ra từ quá trình nhiệt phân gồm $\text{CO}_2$ ($45 - 50\text{ vol}%$), $\text{CO}$ ($30 - 35\text{ vol}%$), $\text{CH}_4$ ($10 - 15\text{ vol}%$) và một lượng nhỏ $\text{H}_2, \text{C}_2\text{H}_4, \text{C}_2\text{H}_6$ ($3 - 5\text{ vol}%$). Khí có nhiệt trị trung bình $10.5 - 12.0\text{ MJ/Nm}^3$, hoàn toàn có thể tuần hoàn đốt cháy để cấp nhiệt cho lò phản ứng, thiết lập chu trình năng lượng tự duy trì.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình giải tích chuẩn xác về truyền nhiệt hạt sinh khối có xét đến biến thiên chuẩn số Biot, bổ sung cơ sở dữ liệu động học Arrhenius của sinh khối nhiệt đới vào hệ thống lý thuyết kỹ thuật nhiệt thế giới.
  • Về mặt phương pháp luận: Thiết lập quy trình nghiên cứu tích hợp chuẩn mực: Mô hình hóa giải tích $\to$ Mô phỏng số CFD $\to$ Chế tạo thực nghiệm $\to$ Phân tích hóa lý chuyên sâu $\to$ Xác lập động học phản ứng. Phương pháp này có thể chuyển giao áp dụng cho các loại sinh khối khác như trấu, rơm rạ, vỏ cà phê.
  • Về mặt thực tiễn kỹ thuật: Cung cấp bộ thông số thiết kế hoàn chỉnh cho lò tầng sôi công suất $500\text{ g/h}$ (kích thước buồng phản ứng, tỷ lệ chiều cao lớp cát tĩnh $h_{\text{cát}} = 50\text{ mm}$, đường kính lỗ phân phối khí, diện tích trao đổi nhiệt bình ngưng $F = 0.28\text{ m}^2$), làm tiền đề cho việc mở rộng quy mô pilot công nghiệp.
  • Về mặt chính sách và môi trường: Đóng góp luận cứ kỹ thuật cho Chiến lược Phát triển Năng lượng Tái tạo của Việt Nam đến năm 2030, tầm nhìn 2050, thúc đẩy mục tiêu Net Zero Carbon và chuyển đổi kinh tế tuần hoàn trong nông lâm nghiệp.

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 hạn chế khoa học mang tính điều kiện biên:

  1. Quy mô công nghệ: Thí nghiệm được thực hiện ở quy mô phòng thí nghiệm ($500\text{ g/h}$) với chế độ cấp nhiệt gián tiếp qua thành lò bằng điện trở. Trong quy mô công nghiệp lớn, cơ chế cấp nhiệt nội tại (qua ống truyền nhiệt chìm trong lớp sôi hoặc tuần hoàn cát mang nhiệt) sẽ xuất hiện các hiệu ứng thủy động học phức tạp hơn.
  2. Ảnh hưởng của kim loại kiềm trong tro: Luận án chưa định lượng chi tiết vai trò xúc tác tự nhiên của các ion kim loại kiềm và kiềm thổ ($\text{K}^+, \text{Na}^+, \text{Ca}^{2+}, \text{Mg}^{2+}$) sẵn có trong tro sinh khối đối với sự cạnh tranh giữa các phản ứng phân hủy sơ cấp.
  3. Độ ổn định của dầu sinh học: Dầu sinh học sau nhiệt phân chưa qua xử lý nâng cấp (upgrading) vẫn còn chứa hàm lượng oxy cao ($\sim 40\text{ wt}%$) và tính axit ($pH < 3.0$), dẫn đến hiện tượng tự trùng ngưng (aging effect) làm tăng độ nhớt sau thời gian lưu trữ dài.
  4. Giả thiết hình học hạt: Mô hình toán giả thiết hạt sinh khối là hình cầu đẳng hướng lý tưởng, trong khi thực tế hạt bột gỗ và bã mía có hình que và dạng sợi bất đẳng hướng về hệ số dẫn nhiệt theo phương dọc thớ và ngang thớ.

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

  • Nghiên cứu nâng cấp xúc tác tại chỗ (In-situ Catalytic Fast Pyrolysis): Bổ sung xúc tác Zeolite H-ZSM-5, MCM-41 hoặc xúc tác oxit kim loại biến tính vào lớp sôi nhằm khử oxy (deoxygenation), nâng cao chỉ số pH ($>5.0$) và tăng tỷ lệ hydrocarbon thơm.
  • Mở rộng quy mô pilot liên tục $10 - 50\text{ kg/h}$: Nghiên cứu công nghệ tầng sôi tuần hoàn (Circulating Fluidized Bed) sử dụng chính khí không ngưng và cốc sinh học làm nhiên liệu tự cấp nhiệt.
  • Mô hình hóa truyền nhiệt bất đẳng hướng: Nâng cấp mô hình toán học giải tích 3D xét đến tính bất đẳng hướng của cấu trúc sợi sinh khối thực tế.
  • Đánh giá vòng đời (LCA) và hiệu quả kinh tế - kỹ thuật (Techno-Economic Analysis - TEA): Lập mô hình tính toán chi phí vòng đời và mức giảm phát thải $\text{CO}_2$ tương đương trên toàn bộ chuỗi cung ứng từ thu gom sinh khối đến sản xuất nhiên liệu thương phẩm.

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

Công trình nghiên cứu mang lại những tác động sâu rộng trên nhiều phương diện:

  • Tác động học thuật: Thiết lập công bố khoa học trên các tạp chí chuyên ngành uy tín trong nước và quốc tế (Tạp chí Khoa học và Công nghệ Đại học Đà Nẵng, Tạp chí Năng lượng Nhiệt), cung cấp nguồn trích dẫn nền tảng cho các nghiên cứu sinh và học viên cao học ngành Kỹ thuật Năng lượng và Kỹ thuật Hóa học.
  • Tác động công nghiệp: Cung cấp giải pháp kỹ thuật khả thi cho các tập đoàn lâm nghiệp và nhà máy mía đường trong việc xử lý hàng triệu tấn phế phẩm hàng năm, tạo ra nguồn dầu sinh học thay thế trực tiếp cho dầu FO thương phẩm trong các lò hơi công nghiệp với chi phí nhiên liệu giảm $20 - 30%$.
  • Tác động chính sách: Cung cấp dữ liệu thực nghiệm phục vụ việc xây dựng Tiêu chuẩn Quốc gia (TCVN) về nhiên liệu lỏng sinh học phi thực phẩm và định hướng chính sách hỗ trợ giá cho các dự án năng lượng sinh khối.
  • Tác động xã hội và toàn cầu: Giúp giảm thiểu triệt để tình trạng đốt lộ thiên phế phẩm nông lâm nghiệp gây ô nhiễm không khí cục bộ, giảm phát thải bụi mịn $\text{PM}_{2.5}$, tạo sinh kế và gia tăng giá trị kinh tế cho chuỗi sản xuất nông nghiệp tại các vùng nông thôn Miền Trung - Tây Nguyên.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận mô hình toán giải tích truyền nhiệt hoàn chỉnh, hệ phương trình nghiệm Fourier cải tiến và bộ dữ liệu động học $E_a, A$ chuẩn xác để phát triển các mô hình mô phỏng chuyên sâu.
  • Kỹ sư thiết kế và Bộ phận R&D doanh nghiệp năng lượng: Sử dụng trực tiếp các công thức tính toán thiết kế lò tầng sôi, buồng đốt, cyclone thu bụi và hệ thống ngưng tụ ống chùm năng suất $500\text{ g/h}$ để phát triển các module thiết bị thương mại.
  • Các nhà hoạch định chính sách năng lượng: Có cơ sở khoa học định lượng về hiệu suất chuyển hóa năng lượng sinh khối để ban hành các quy chuẩn, cơ chế khuyến khích đầu tư năng lượng tái tạo.
  • Ngành công nghiệp chế biến gỗ và mía đường: Nắm bắt công nghệ cốt lõi để tự chủ nguồn năng lượng tại chỗ, giảm thiểu chi phí xử lý môi trường và gia tăng giá trị gia tăng của phụ phẩm.

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ì? Đó là việc thiết lập thành công nghiệm giải tích dạng chuỗi Fourier cho phương trình vi phân truyền nhiệt phi dừng trong hạt sinh khối hình cầu có điều kiện biên loại 3 phức hợp, kết hợp với phương trình phân hủy nhiệt cạnh tranh bậc một. Mô hình đã chứng minh bằng giải tích và thực nghiệm rằng bán kính hạt $R \le 0.25\text{ mm}$ là điều kiện biên tới hạn để đưa quá trình nhiệt phân chuyển từ chế độ khống chế truyền nhiệt sang chế độ khống chế động học trong môi trường lò tầng sôi.

2. Đổi mới phương pháp luận của luận án so với các nghiên cứu trước đây? So với các nghiên cứu của Bridgwater (2012) hay Di Blasi (2008) vốn thường tách rời giữa mô phỏng thủy động lực học và thực nghiệm hóa học, luận án đã tích hợp phương pháp tiếp cận 3 cấp độ: Mô hình hóa vi mô (MATLAB) $\to$ Mô phỏng khí động học trung mô (ANSYS Fluent) $\to$ Kiểm chứng thực nghiệm vĩ mô (Hệ thống tầng sôi $500\text{ g/h}$). Phương pháp xác định động học không qua TGA tĩnh thông thường mà được nội suy trực tiếp từ sản phẩm thực tế thu hồi trong lò tầng sôi động, phản ánh đúng bản chất truyền nhiệt đối lưu cưỡng bức trong thực tế.

3. Phát hiện thực nghiệm nào gây bất ngờ và có ý nghĩa kỹ thuật cao nhất? Phát hiện về sự chênh lệch nhiệt độ tối ưu giữa bột gỗ ($500^\circ\text{C}$) và bã mía ($480^\circ\text{C}$) cùng hàm lượng đường Levoglucosan và các dẫn xuất Phenol cao vượt trội ($>50%$ tổng diện tích pic GC-MS) trong dầu sinh học. Điều này chứng minh bio-oil từ nhiệt phân nhanh không chỉ là nguồn nhiên liệu đốt đơn thuần mà còn là nguồn nguyên liệu hóa chất sinh học giá trị cao để sản xuất nhựa phenol-formaldehyde và dung môi xanh.

4. Quy trình thực nghiệm có đảm bảo khả năng tái lặp (Replication protocol)? Hoàn toàn đảm bảo. Luận án đã mô tả chi tiết quy cách hình học của lò tầng sôi ($D = 50\text{ mm}, H = 600\text{ mm}$), thông số hạt cát đệm ($d_p = 0.2 - 0.3\text{ mm}, m_{\text{cát}} = 150\text{ g}$), lưu lượng $N_2$ ($1.2\text{ m}^3\text{/h}$), tần số biến tần cấp liệu vít tải ($10 - 25\text{ Hz}$), cùng các quy chuẩn phân tích quốc tế (ASTM D5291, ASTM D240, ASTM D445, ASTM E203), cho phép các phòng thí nghiệm độc lập tái lập chính xác toàn bộ kết quả.

5. Lộ trình phát triển công nghệ 10 năm được phác thảo như thế nào? Lộ trình 10 năm gồm 3 giai đoạn: (1) 2018 - 2021: Hoàn thiện công nghệ nâng cấp xúc tác bio-oil trực tiếp tại chỗ bằng Zeolite H-ZSM-5; (2) 2022 - 2025: Thiết kế chế tạo hệ thống pilot tầng sôi tuần hoàn liên tục công suất $50\text{ kg/h}$ tự cấp nhiệt; (3) 2026 - 2028: Xây dựng nhà máy nhiệt phân nhanh công nghiệp quy mô $1\text{ tấn/h}$ tích hợp tại các nhà máy chế biến gỗ và đường Miền Trung.

Kết luận

Luận án Tiến sĩ Kỹ thuật của NCS. Phạm Duy Vũ đã giải quyết trọn vẹn và xuất sắc các mục tiêu nghiên cứu đề ra với 5 đóng góp học thuật và thực tiễn cốt lõi:

  1. Hoàn thiện lý thuyết truyền nhiệt phi dừng trong nhiệt phân sinh khối: Xây dựng thành công mô hình toán học giải tích mô tả chính xác trường nhiệt độ $t(r, \tau)$ và xác lập công thức tính thời gian nhiệt phân $\tau_p$ phụ thuộc vào chuẩn số $Bi, Fo$ và bán kính hạt $R$.
  2. Làm chủ kỹ thuật mô phỏng và thiết kế lò tầng sôi bọt: Ứng dụng mô hình Euler-Euler và KTGF xác định chính xác các thông số thủy động học lớp sôi, chế tạo thành công hệ thống nhiệt phân nhanh đồng bộ công suất $500\text{ g/h}$ vận hành ổn định và an toàn.
  3. Tối ưu hóa hiệu suất thu hồi dầu sinh học đạt chuẩn quốc tế: Xác lập chế độ vận hành tối ưu ($T = 480 - 500^\circ\text{C}$, kích thước hạt $d \le 0.5\text{ mm}$, lưu lượng khí $N_2 = 1.2 - 1.4\text{ m}^3\text{/h}$), đạt hiệu suất thu hồi dầu sinh học kỷ lục $64.8\text{ wt}%$ (bột gỗ) và $61.2\text{ wt}%$ (bã mía).
  4. Xác lập hệ dữ liệu động học Arrhenius chuẩn xác: Cung cấp bộ thông số năng lượng hoạt hóa $E_a$ và hằng số tốc độ $A$ cho 3 nhánh phản ứng phân hủy cạnh tranh bậc một của sinh khối gỗ keo và bã mía Việt Nam trong môi trường tầng sôi.
  5. Định lượng toàn diện đặc tính hóa lý của sản phẩm: Phân tích chi tiết nhiệt trị, độ nhớt, độ pH, hàm lượng nước và thành phần phân tử hữu cơ của bio-oil qua GC-MS, khẳng định tiềm năng ứng dụng kép vừa làm nhiên liệu vừa làm nguyên liệu hóa chất xanh.

Công trình tạo nền tảng vững chắc, mở ra 3 hướng nghiên cứu mới về nhiệt phân xúc tác nâng cấp dầu, công nghệ tầng sôi tuần hoàn tự cấp nhiệt và công nghệ chuyển hóa sinh khối tích hợp biorefinery, góp phần đưa nền khoa học kỹ thuật nhiệt và năng lượng tái tạo Việt Nam hội nhập sâu rộng với trình độ phát triển của thế giới.