Tổng quan về luận án

Sự cấp thiết của việc chuyển dịch cơ cấu năng lượng và giảm thiểu phát thải độc hại từ động cơ đốt trong đang là thách thức trung tâm của kỹ thuật cơ khí động lực toàn cầu. Nguồn năng lượng hóa thạch ngày càng cạn kiệt với trữ lượng ước tính chỉ đủ đáp ứng nhu cầu khai thác trong khoảng 50 năm, trong khi việc đốt cháy nhiên liệu khoáng phát thải khối lượng lớn các tác nhân gây ô nhiễm nghiêm trọng như hydrocarbon ($HC$), carbon monoxide ($CO$), oxide của nitrogen ($NO_x$), muội than ($soot$) và bụi hạt ($particulate\ matter - PM$). Trong bối cảnh Việt Nam thực hiện "Đề án phát triển nhiên liệu sinh học đến năm 2015, tầm nhìn đến năm 2025" theo Quyết định số 177/2007/QĐ-TTg và lộ trình phối trộn theo Quyết định số 53/2012/QĐ-TTg, ngành công nghiệp ethanol sinh học nội địa đã xây dựng 7 nhà máy với tổng công suất thiết kế lên tới $600.000\ m^3/\text{năm}$. Tuy nhiên, sản lượng ethanol tiêu thụ thực tế để pha chế xăng sinh học E5 RON 92 chỉ đạt khoảng $177.900\ m^3/\text{năm}$ (tương đương 44% công suất vận hành của 4 nhà máy chủ lực gồm Tùng Lâm, Bình Phước, Dung Quất), dẫn đến tình trạng dư thừa nghiêm trọng nguồn cồn sinh học sản xuất từ sắn. Ngược lại, nguồn cung diesel sinh học (biodiesel) trong nước lại khan hiếm do thiếu hụt vùng nguyên liệu dầu thực vật tập trung quy mô công nghiệp. Do đó, việc nghiên cứu ứng dụng ethanol để thay thế một phần diesel khoáng cho đội xe thương mại tải trọng nhẹ và trung bình là hướng đi mang tính đột phá khoa học và thực tiễn cấp bách.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định: Việc pha trộn trực tiếp ethanol vào dầu diesel gặp phải các rào cản vật lý - hóa học căn bản bao gồm độ bền tách pha kém ở nhiệt độ thấp ($T < 10^\circ C$), trị số cetane thấp (ethanol có $CN \approx 6-8$ so với diesel $CN \approx 48-51.5$), độ nhớt động học quá thấp ($1.2\ cP$ ở $20^\circ C$ so với $2.95\ cP$ của diesel) gây nguy cơ rò rỉ và mòn cặp đôi piston - xylanh bơm cao áp, cùng nhiệt ẩn hóa hơi cao ($854\ kJ/kg$ so với $301\ kJ/kg$) làm suy giảm nhiệt độ môi chất cuối quá trình nén và kéo dài thời gian cháy trễ. Các công trình quốc tế đã khảo sát hỗn hợp nhị phân diesel-ethanol (DE) hoặc diesel-biodiesel (DB), nhưng các nghiên cứu về hệ nhiên liệu tam phân (ternary blend) diesel-ethanol-biodiesel (DEB) sử dụng phụ gia cải thiện là biodiesel gốc dầu cọ (Palm Oil Methyl Ester - PME) trên các dòng động cơ diesel cơ khí đang lưu hành tại các quốc gia đang phát triển còn rất hạn chế, thiếu sự gắn kết chặt chẽ giữa mô phỏng nhiệt động lực học 1D, khảo sát quang học chùm tia phun tốc độ cao và đánh giá thực nghiệm phát thải hạt bụi chuyên sâu.

Luận án xác lập 3 câu hỏi nghiên cứu (Research Questions - RQ) và 3 giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Sự biến thiên tỷ lệ ethanol ($5 - 30%$) và biodiesel ($5%$) tác động như thế nào đến tính chất lý hóa, cấu trúc chùm tia phun và quá trình bay hơi tạo hỗn hợp trong buồng cháy động cơ diesel?
    • H1: Hàm lượng $5%$ biodiesel đóng vai trò chất đồng dung môi (co-solvent) giúp ổn định pha hỗn hợp, đồng thời bù đắp sự suy giảm độ nhớt và trị số cetane do ethanol gây ra mà không đòi hỏi thay đổi kết cấu hệ thống cấp nhiên liệu.
  • RQ2: Quy luật diễn biến quá trình cháy, áp suất chỉ thị cực đại ($P_{max}$), tốc độ tỏa nhiệt ($ROHR$) và cơ chế hình thành phát thải độc hại ($CO, HC, NO_x, soot, PN$) thay đổi ra sao dưới tác động đồng thời của hiệu ứng nhiệt ẩn hóa hơi cao và hiệu ứng giàu oxy nội tại?
    • H2: Tỷ lệ hòa trộn ethanol tối ưu sẽ cắt giảm đột phá phát thải muội than ($soot$) và nồng độ số lượng hạt bụi ($PN$) nhờ sự xuất hiện của oxy liên kết trong phân tử $C_2H_5OH$ ($34.73%$) và biodiesel ($10.9%$), trong khi duy trì tính năng kinh tế - kỹ thuật ở mức tương đương diesel khoáng.
  • RQ3: Mô hình mô phỏng nhiệt động lực học trên phần mềm chuyên dụng AVL Boost có khả năng dự báo chính xác các chỉ tiêu công tác của động cơ khi sử dụng nhiên liệu tam phân DEB hay không?
    • H3: Việc hiệu chỉnh mô hình cháy và mô hình truyền nhiệt theo dữ liệu thực nghiệm sẽ đạt độ hội tụ cao với sai số dự báo dưới $5%$ đối với công suất, mô-men và suất tiêu hao nhiên liệu.

Khung lý thuyết (theoretical framework) của nghiên cứu được xây dựng dựa trên lý thuyết chu trình nhiệt động cấp nhiệt hỗn hợp Sabathé - Seiliger, mô hình cháy kiểm soát bằng hòa trộn (Mixing Controlled Combustion - MCC) kết hợp mô hình Wiebe hai pha, lý thuyết phân rã tán sương tia phun Rayleigh - Taylor và Kelvin - Helmholtz, cơ chế động học Zeldovich mở rộng cho phản ứng tạo $NO_x$, cùng mô hình tạo muội than Hiroyasu kết hợp oxy hóa Nagle - Strickland - Constable.

Đóng góp đột phá của luận án nằm ở việc lượng hóa chính xác tỷ lệ phối trộn tối ưu DE10B5 ($85%$ diesel khoáng + $10%$ ethanol sinh học từ sắn + $5%$ biodiesel từ dầu cọ theo thể tích). Nhiên liệu DE10B5 cho phép thay thế trực tiếp $10%$ diesel hóa thạch, giúp giảm phát thải muội than và số lượng hạt mịn tới $30 - 50%$ trên toàn dải tải, trong khi mô-men xoắn chỉ suy giảm không đáng kể ($< 2.5%$) và suất tiêu hao nhiên liệu riêng tăng nhẹ ($3 - 5%$). Phạm vi nghiên cứu được thực hiện trên động cơ diesel Hyundai D4BB (4 xylanh, dung tích $2.607\ cm^3$, hệ thống phun nhiên liệu gián tiếp sử dụng bơm cao áp cơ khí phân phối VE) trong điều kiện phòng thí nghiệm tiêu chuẩn ở chế độ đường đặc tính ngoài ($1000 - 4000\ v/ph$) và đường đặc tính tải tại tốc độ định mức $2000\ v/ph$.

Literature Review và Positioning

Tổng quan y văn quốc tế ghi nhận ba dòng nghiên cứu chính trong lĩnh vực ứng dụng nhiên liệu sinh học cho động cơ diesel:

Dòng nghiên cứu hỗn hợp nhị phân Diesel - Ethanol (DE): Giakoumis et al. (2011, 2013) khi khảo sát hỗn hợp DE5 và DE10 trên động cơ tăng áp 6 xylanh phun trực tiếp đã chứng minh độ khói giảm mạnh khi tăng tỷ lệ ethanol, $NO_x$ ít biến động, $CO$ giảm nhưng phát thải $HC$ và suất tiêu hao nhiên liệu có xu hướng tăng nhẹ. E. Ajav et al. (1999) thử nghiệm hỗn hợp DE5 đến DE20 trên động cơ một xylanh TV110 cho thấy công suất động cơ giảm tuyến tính theo tỷ lệ ethanol từ $0.46%$ đến $2.98%$. Venkata Subbaiah et al. (2010) và Bang-Quan He et al. (2003) phân tích sâu hơn về sự sụt giảm áp suất phun và chỉ ra rằng với DE10 và DE30, hàm lượng bồ hóng giảm đáng kể nhưng $NO_x$ có xu hướng gia tăng ở các chế độ tải nặng do sự xuất hiện của oxy phân tử trong cồn.

Dòng nghiên cứu hỗn hợp Diesel - Biodiesel (DB): Ekrem Buyukkaya (2010) thử nghiệm các tỷ lệ B5, B20, B70, B100 gốc dầu hạt cải (Rapeseed Methyl Ester - RME) trên động cơ 6 xylanh tăng áp ($164\ kW$), xác nhận công suất ở mẫu B5 hoàn toàn tương đương diesel, nhưng ở B20, B70 và B100 công suất giảm lần lượt $6\ kW, 8\ kW$ và $10\ kW$; phát thải $CO, HC$ và độ khói giảm tới $45 - 60%$, trong khi $NO_x$ tăng từ $6%$ đến $12%$. B. Tesfa et al. (2010) và Hansen et al. (2005) lý giải sự gia tăng tiêu hao nhiên liệu ($10 - 15%$ với B50, B100) xuất phát từ nhiệt trị thấp hơn ($13.3%$) cùng độ nhớt và khối lượng riêng cao của biodiesel. Sahoo PK et al. (2009) so sánh ba nguồn biodiesel từ Jatropha, Karanja và Polanga, khẳng định phát thải bụi hạt PM giảm từ $28.57%$ đến $69.48%$ nhưng $NO_x$ luôn tăng tương ứng. Raheman & Phadatare (2004) ghi nhận mức giảm $CO$ từ $75%$ đến $94%$ với biodiesel gốc Karanja. Agarwal et al. (2007) phân tích đặc tính cháy chỉ ra rằng trị số cetane cao của biodiesel làm rút ngắn giai đoạn cháy trễ, dịch chuyển thời điểm bắt đầu cháy sớm hơn và hạ thấp áp suất đỉnh so với diesel khoáng.

Dòng nghiên cứu hệ tam phân Diesel - Ethanol - Biodiesel (DEB): Al-Hassan et al. (2012) chứng minh hỗn hợp nhị phân DE5 - DE20 phân tách pha sau 2 đến 80 giờ, nhưng khi bổ sung $10%$ biodiesel, thời gian bền pha của DE10B10 và DE20B10 kéo dài lên 1 đến 9 ngày, riêng DE5B10 hoàn toàn không tách pha. Qiang Fang et al. (2013) nghiên cứu trên động cơ Common Rail 4 xylanh với biodiesel từ dầu ăn phế thải (Waste Cooking Oil Methyl Ester - WME), ghi nhận tốc độ giải phóng nhiệt cực đại tăng khi tăng tỷ lệ ethanol do kéo dài thời gian cháy trễ. X. Shi et al. (2006) thử nghiệm DE3B12 và DE4B16 (biodiesel dầu nành - SME) trên động cơ Euro II, ghi nhận PM giảm tới $48%$, song $NO_x$ tăng tối đa $32%$. T. Krishnaswamy et al. (2013) với DE10B10 gốc Jatropha và Czerwinski et al. (2007) với DE30B15 ghi nhận $NO_x$ và PM đồng thời giảm ở chế độ toàn tải, nhưng $CO$ và $HC$ tăng tại vùng tải nhỏ do nhiệt độ buồng cháy bị dập tắt bởi ethanol.

                           KHÔNG GIAN NGHIÊN CỨU Y VĂN QUỐC TẾ

Trong y văn tồn tại hai cuộc tranh luận học thuật sâu sắc:

  1. Tranh luận về phát thải $NO_x$: Một nhóm tác giả (Shi, 2006; Venkata Subbaiah, 2010) cho rằng hàm lượng oxy cao trong hỗn hợp DEB thúc đẩy hình thành $NO_x$ nhiệt theo cơ chế Zeldovich; trong khi nhóm tác giả khác (Czerwinski, 2007; Nadir Yilmaz, 2014; Su Han Park, 2012) chứng minh nhiệt ẩn hóa hơi lớn của ethanol hấp thụ mạnh nhiệt lượng buồng cháy, kéo tụt nhiệt độ đỉnh ngọn lửa và triệt tiêu sự gia tăng của $NO_x$.
  2. Tranh luận về áp suất cực đại và tốc độ tăng áp: Nghiên cứu của Qiang Fang (2013) chỉ ra $P_{max}$ của DEB thấp hơn diesel ở tải trung bình do trễ bắt cháy, nhưng tăng vọt ở chế độ toàn tải do lượng hòa khí tích lũy trong pha cháy nhanh bùng phát mạnh mẽ.

Định vị của luận án: Trong khi các công trình quốc tế chủ yếu tập trung vào động cơ Common Rail hiện đại hoặc động cơ nghiên cứu 1 xylanh với nguồn biodiesel đắt tiền (Jatropha, hạt cải, đậu nành), luận án này giải quyết trực tiếp bài toán thực chứng trên động cơ đa xylanh có hệ thống phun cơ khí truyền thống chiếm đa số tại Việt Nam. Nghiên cứu sử dụng nguồn cồn sắn khan ($> 99.5%$) và biodiesel dầu cọ PME nội địa, tích hợp chuỗi bằng chứng từ mô phỏng nhiệt động 1D trên AVL Boost, giải mã cơ chế tán sương bằng camera quang học tốc độ cao 3000 FPS đến kiểm nghiệm phát thải hạt bụi chuyên sâu theo quy chuẩn quốc tế.

Đó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 lý thuyết quá trình cháy nén động cơ đốt trong (Taylor, Heywood) đối với hệ nhiên liệu đa thành phần có sự phân cực và khác biệt pha nhiệt động sâu sắc. Khi cồn ethanol (chất phân cực, chứa nhóm $-OH$) được phối trộn vào diesel khoáng (chuỗi hydrocarbon không phân cực) dưới sự hiện diện của biodiesel (các alkyl ester axit béo có cấu trúc lưỡng tính), biodiesel đóng vai trò chất hoạt động bề mặt làm giảm sức căng liên bề mặt, ngăn chặn quá trình tách pha lỏng - lỏng ở mức độ phân tử.

Nghiên cứu làm sâu sắc thêm lý thuyết động học hòa trộn và phản ứng cháy thông qua việc lượng hóa tác động tương hỗ giữa sự suy giảm trị số cetane và hiệu ứng làm giàu oxy (oxygen enrichment). Hàm lượng oxy liên kết trong phân tử $C_2H_5OH$ ($34.73%$) và biodiesel ($10.9%$) tham gia trực tiếp vào phản ứng chuỗi gốc tự do tại vùng lõi tia phun nghèo oxy, làm thay đổi căn bản cơ chế nhiệt phân hydrocarbon thành tiền chất muội than (polycyclic aromatic hydrocarbons - PAH).

Luận án xác lập và kiểm chứng công thức thực nghiệm tính toán trị số cetane ($CN_{ebdg}$) và độ nhớt động học ($\eta_{kev}$) cho hệ tam phân diesel-ethanol-biodiesel:

$$CN_{ebdg} = CN_d - 0.59 \cdot n_e + 0.55 \cdot n_{bd}$$

Trong đó $CN_d$ là trị số cetane của diesel gốc, $n_e$ và $n_{bd}$ lần lượt là phần trăm thể tích của ethanol và biodiesel. Phương trình này chứng minh trên phương diện định lượng rằng mỗi $1%$ biodiesel bổ sung có khả năng bù đắp tương đương cho sự suy giảm cetane gây ra bởi $0.93%$ ethanol, tạo cơ sở lý thuyết cho việc tối ưu hóa tỷ lệ pha chế.

Đồng thời, hàm lượng oxy tổng thể trong nhiên liệu hỗn hợp ($M_o$, tính bằng gam) được lượng hóa chính xác qua biểu thức:

$$M_o = 0.34 \cdot \rho_e \cdot V_e + 0.11 \cdot \rho_b \cdot V_b$$

với $\rho_e, \rho_b$ và $V_e, V_b$ là khối lượng riêng và thể tích của các cấu tử ethanol và biodiesel.

       MÔ HÌNH KHUNG PHÂN TÍCH TÍCH HỢP QUÁ TRÌNH NHIỆT ĐỘNG VÀ PHÁT THẢI

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

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết chuyên sâu:

  1. Lý thuyết nhiệt động học chu trình công tác Sabathé: Áp dụng định luật thứ nhất và thứ hai nhiệt động lực học trong hệ thống hở có biến thiên thành phần hóa học và tỷ nhiệt phụ thuộc nhiệt độ môi chất. Quá trình cháy được mô hình hóa thông qua mô hình tỏa nhiệt đa vùng MCC trên nền tảng phần mềm AVL Boost.
  2. Lý thuyết cơ học chất lỏng và tán xạ tia phun: Phân tích động lực học dòng tia hai pha lỏng - khí qua lỗ phun vòi cơ khí, mô tả quá trình phân rã nguyên phát (primary breakup) và thứ phát (secondary breakup) chịu tác động bởi sự thay đổi mật độ riêng ($\rho$), độ nhớt ($\mu$) và sức căng bề mặt ($\sigma$).
  3. Động học hóa học phát thải: Ứng dụng cơ chế nhiệt Zeldovich cải biên (gồm chuỗi 3 phản ứng thuận nghịch của các gốc $N, NO, O, OH$) để tính toán $NO_x$, kết hợp mô hình động học bán thực nghiệm hai bước của Hiroyasu - Magnussen để mô tả tốc độ tạo thành và oxy hóa muội than:

$$\frac{dm_{soot}}{dt} = \frac{dm_{form}}{dt} - \frac{dm_{oxid}}{dt}$$

Điều kiện biên (boundary conditions): Khung phân tích được xác lập trong giới hạn nhiệt độ nhiên liệu $T \ge 10^\circ C$, áp suất phun định mức của bơm phân phối cơ khí ($120 - 250\ bar$), tỷ lệ ethanol $\le 30%$ và biodiesel $\le 5%$ về thể tích nhằm đảm bảo độ tin cậy cơ học của động cơ không qua cải tạo.

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi hệ hình thực chứng luận (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), tiếp cận đối tượng thông qua chuỗi thực nghiệm khách quan có kiểm soát, kết hợp mô phỏng số học để giải thích cơ chế nội tại.

Thiết kế nghiên cứu đa cấp độ (multi-level design) được triển khai tuần tự:

  • Cấp độ vi mô: Xác định các thuộc tính hóa lý cơ bản của các mẫu nhiên liệu (Diesel, DE5, DE10, DE15, DE20, DE5B5, DE10B5) gồm tỷ trọng ở $20^\circ C$ ($0.8279 - 0.8437$), nhiệt trị thấp ($39.62 - 42.15\ MJ/kg$), độ nhớt động học ($1.802 - 2.435\ cP$), hàm lượng oxy tự do ($1.75 - 6.95%$) và nhiệt ẩn hóa hơi ($319 - 407\ kJ/kg$).
  • Cấp độ trung mô: Khảo sát cấu trúc tia phun không bay hơi và bay hơi thông qua hệ thống buồng áp suất quang học đồng bộ máy quay siêu tốc độ nét cao.
  • Cấp độ vĩ mô: Thử nghiệm toàn diện trên băng thử động cơ tĩnh nhằm thu thập các dữ liệu động lực học, tiêu hao nhiên liệu, áp suất trong xylanh theo góc quay trục khuỷu ($P-\varphi$), phát thải khí ô nhiễm và phổ phân bố hạt bụi.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu kết hợp giữa mô hình số và thực nghiệm kiểm chứng được chuẩn hóa theo các tiêu chuẩn quốc tế nghiêm ngặt:

  1. Quy trình tạo mẫu và kiểm tra độ bền pha: Mẫu nhiên liệu được phối chế bằng thiết bị khuấy cơ học tốc độ cao trong 30 phút ở nhiệt độ phòng. Quá trình kiểm tra độ ổn định quang học và hiện tượng phân lớp được quan sát liên tục trong buồng ổn nhiệt.
  2. Hệ thống quang học chụp tia phun: Vòi phun nguyên bản của động cơ D4BB được gắn trên thiết bị cân chỉnh áp suất mở kim phun ($120\ bar$). Chùm tia phun được rọi sáng bởi nguồn sáng liên tục công suất cao và ghi hình bằng máy quay siêu tốc Photron Fastcam SA3 với tốc độ 3000 khung hình/giây ($3000\ FPS$). Các thông số hình học tia phun như chiều dài tán sương ($L$) và góc nón tia phun ($\theta$) được xử lý qua phần mềm phân tích ảnh số chuyên dụng.
  3. Băng thử động lực học và thiết bị phân tích khí xả:
    • Động cơ D4BB được kết nối trực tiếp với phanh điện động lực cao APA 100 có khả năng điều khiển chính xác mô-men tải và tốc độ vòng quay.
    • Suất tiêu hao nhiên liệu được ghi nhận tự động theo thời gian thực bằng cân đo nhiên liệu trọng lượng AVL Fuel Balance 733S.
    • Thành phần khí thải được dẫn qua hệ thống ống gia nhiệt ($190^\circ C$) vào tủ phân tích khí xả AVL CEBII: đo $CO, CO_2$ bằng cảm biến hồng ngoại không tán sắc (NDIR), đo $HC$ bằng đầu dò ion hóa ngọn lửa gia nhiệt (HFID), đo $NO_x$ bằng detector phát quang hóa học gia nhiệt (HCLD), và xác định hệ số dư lượng không khí $\lambda$.
    • Độ đen của khói được xác định bằng thiết bị đo độ khói lọc giấy AVL 415S Smoke Meter (đơn vị FSN - Filter Smoke Number).
    • Nồng độ số lượng hạt bụi ($PN$) và phổ kích thước hạt được đo bằng hệ thống lấy mẫu pha loãng hai cấp kết hợp buồng bốc hơi nhiệt hóa và thiết bị đếm hạt ngưng tụ (Condensation Particle Counter - CPC) theo quy chuẩn PMP (Particle Measurement Programme).

Data và phân tích

Thông số kỹ thuật của đối tượng thử nghiệm động cơ Hyundai D4BB:

  • Kiểu loại: Động cơ diesel 4 kỳ, 4 xylanh thẳng hàng, làm mát bằng nước.
  • Đường kính xylanh $\times$ Hành trình piston: $91.1\ mm \times 100\ mm$.
  • Dung tích công tác: $2.607\ cm^3$.
  • Tỷ số nén: $\varepsilon = 22.0$.
  • Công suất định mức catalog: $59\ kW$ ($80\ HP$) tại $4000\ v/ph$.
  • Mô-men xoắn cực đại catalog: $167\ N.m$ tại $2200\ v/ph$.
  • Hệ thống cung cấp nhiên liệu: Bơm cao áp phân phối Bosch VE, điều tốc cơ năng, vòi phun kín nhiều lỗ.

Dữ liệu thực nghiệm được thu thập tại các điểm vận hành ổn định:

  • Đường đặc tính ngoài: Tải $100%$ vị trí bàn đạp ga, quét dải tốc độ từ $1000\ v/ph$ đến $4000\ v/ph$ với bước nhảy $500\ v/ph$.
  • Đường đặc tính tải: Giữ cố định tốc độ vòng tua tại $2000\ v/ph$ (vùng mô-men xoắn cực đại và dải hoạt động phổ biến nhất của xe tải nhẹ), biến thiên mức tải ở $25%, 50%, 75%$ và $100%$ tải định mức.

Kiểm định độ vững (Robustness checks): Mô hình mô phỏng 1D trên AVL Boost được hiệu chuẩn lặp (iterative calibration) dựa trên số liệu thực nghiệm của diesel khoáng chuẩn. Kết quả đối chuẩn cho thấy sai lệch giữa mô phỏng (MP) và thực nghiệm (TN) đối với công suất và suất tiêu hao nhiên liệu luôn nằm trong ngưỡng $\pm 3.5%$, đối với phát thải $CO, NO_x, soot$ nằm trong dải sai số chấp nhận được $\pm 8 - 12%$, khẳng định độ tin cậy cao của công cụ dự báo số học.

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

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

  1. Biến đổi tích cực trong động lực học tia phun: Kết quả ghi hình tốc độ cao 3000 FPS cho thấy nhiên liệu hỗn hợp DE10B5 có độ nhớt ($2.283\ cP$) và sức căng bề mặt thấp hơn diesel gốc ($2.95\ cP$), giúp quá trình xé tơi giọt nhiên liệu diễn ra sớm hơn, góc nón tia phun ($\theta$) mở rộng hơn từ $1.5^\circ$ đến $3.2^\circ$ so với diesel khoáng. Mặc dù tầm với tia phun ($L$) giảm nhẹ, diện tích bề mặt tiếp xúc của chùm hạt sương với không khí nén tăng lên rõ rệt, thúc đẩy tốc độ bay hơi hòa trộn trong buồng cháy lốc xoáy.
  2. Duy trì tính năng động lực học và biến thiên suất tiêu hao: Trên đường đặc tính ngoài, mô-men xoắn của động cơ khi sử dụng DE10B5 chỉ giảm bình quân $1.8 - 2.4%$ so với diesel khoáng (ví dụ tại $2200\ v/ph$, mô-men đạt xấp xỉ $163.5\ N.m$ so với $167\ N.m$ của diesel). Suất tiêu hao nhiên liệu có ích ($g_e$) tăng trung bình từ $3.2%$ đến $4.8%$ do nhiệt trị của DE10B5 ($40.1\ MJ/kg$) thấp hơn so với diesel ($42.95\ MJ/kg$). Tuy nhiên, hiệu suất nhiệt có ích ($\eta_e$) của động cơ lại tăng từ $1.1%$ đến $2.3%$ nhờ quá trình cháy diễn ra hoàn thiện hơn.
  3. Đột phá trong việc cắt giảm phát thải bồ hóng (Soot) và bụi hạt (PN): Phát thải độ khói (đo bằng AVL 415S) giảm sâu từ $35%$ đến $52%$ trên toàn bộ dải tải khi sử dụng DE10B5. Tại chế độ $100%$ tải ở $2000\ v/ph$, chỉ số FSN giảm từ $2.15$ xuống còn $1.08$. Nồng độ số lượng hạt bụi ($PN$) trong dải kích thước siêu mịn ($Dp < 100\ nm$) giảm hơn $40%$. Cơ chế giải thích: Sự phân bố oxy nội tại ($3.36%$ trọng lượng) trực tiếp tại vùng lõi giàu nhiên liệu đã ngăn chặn quá trình polyme hóa axetylen thành vòng thơm, đồng thời kéo dài thời gian oxy hóa muội than trong giai đoạn cháy chính và cháy rớt.
  4. Quy luật phát thải $CO, HC$ và $NO_x$:
    • Phát thải $CO$ giảm mạnh từ $15%$ đến $28%$ ở các chế độ tải trung bình và cao nhờ sự có mặt của oxy tự do thúc đẩy chuyển hóa $CO \rightarrow CO_2$.
    • Phát thải $HC$ tăng nhẹ ($5 - 12%$) ở vùng phụ tải thấp ($25%$ tải) do hiệu ứng làm mát cục bộ của ethanol (nhiệt ẩn hóa hơi cao $854\ kJ/kg$) dẫn đến hiện tượng dập tắt màng lửa sát vách buồng cháy (flame quenching), nhưng giảm dần khi tải tăng cao.
    • Phát thải $NO_x$ biến thiên không đáng kể, duy trì ở mức tương đương diesel khoáng (dao động trong biên độ $-2%$ đến $+4%$). Hiện tượng này là kết quả của sự triệt tiêu lẫn nhau giữa hai cơ chế đối nghịch: hiệu ứng giảm nhiệt độ đỉnh ngọn lửa do ethanol bay hơi thu nhiệt và hiệu ứng tăng nồng độ oxy cục bộ thúc đẩy phản ứng Zeldovich từ biodiesel.
       SO SÁNH CÁC CHỈ TIÊU KỸ THUẬT VÀ PHÁT THẢI (DE10B5 SO VỚI DIESEL GỐC)

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Công trình chứng minh tính hiệu quả của mô hình cân bằng nhiệt động - hóa học trong hệ nhiên liệu tam phân; cung cấp bộ dữ liệu thực nghiệm chuẩn mực về tương quan giữa đặc tính phun tán sương, diễn biến áp suất xylanh và phổ phát thải hạt trên động cơ diesel phun gián tiếp.
  • Về mặt phương pháp luận: Thiết lập quy trình nghiên cứu tích hợp từ mô phỏng nhiệt động AVL Boost 1D, khảo sát quang học chùm tia chất lỏng đến phân tích phát thải khí - hạt thực tế, có thể chuyển giao áp dụng cho các nghiên cứu nhiên liệu thay thế khác (như butanol, methanol, dầu thực vật hydro hóa - HVO).
  • Về mặt thực tiễn và công nghiệp: Khẳng định tỷ lệ DE10B5 hoàn toàn tương thích với hệ thống cấp nhiên liệu và động cơ diesel cơ khí đang lưu hành mà không cần thay đổi bất kỳ kết cấu phần cứng nào, không gây mòn kẹt bơm cao áp nhờ đặc tính bôi trơn cao của biodiesel (chỉ số HFRR đạt chuẩn).
  • Về mặt chính sách và kinh tế - xã hội: Đưa ra giải pháp trực tiếp hấp thụ lượng cồn sắn dư thừa ($> 400.000\ m^3/\text{năm}$) cho ngành nông nghiệp và các nhà máy ethanol sinh học tại Việt Nam, giảm gánh nặng nhập khẩu nhiên liệu hóa thạch, đồng thời góp phần cải thiện rõ rệt chất lượng không khí đô thị thông qua việc cắt giảm hơn $40%$ phát thải hạt mịn PM2.5 từ phương tiện vận tải thương mại.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn nghiên cứu (limitations):

  1. Giới hạn về điều kiện nhiệt độ vận hành: Các thử nghiệm được thực hiện ở nhiệt độ phòng thí nghiệm tiêu chuẩn ($20 - 30^\circ C$). Hiện tượng phân lớp pha lỏng của hỗn hợp DEB khi nhiệt độ môi trường giảm sâu xuống dưới $10^\circ C$ trong mùa đông tại vùng núi phía Bắc chưa được khảo sát kèm các phụ gia ức chế đông tụ chuyên biệt.
  2. Giới hạn về chủng loại đối tượng thử nghiệm: Đề tài tập trung chuyên sâu trên động cơ diesel Hyundai D4BB sử dụng bơm cao áp cơ khí phân phối VE. Các kết luận chưa thể khái quát hóa hoàn toàn cho các thế hệ động cơ diesel hiện đại sử dụng hệ thống phun nhiên liệu tích áp Common Rail áp suất siêu cao ($> 2000\ bar$), kim phun điện từ áp điện đa điểm (Piezo injector) và các bộ xử lý xúc tác sau khí xả phức tạp (DOC, DPF, SCR).
  3. Giới hạn về thời gian đánh giá độ bền dài hạn: Nghiên cứu chủ yếu đánh giá các chỉ tiêu tức thời theo chu trình thử nghiệm ổn định. Chưa có điều kiện thực hiện phép thử độ bền mòn động cơ gia tốc liên tục ($500 - 1000\ giờ$) để kiểm chứng mức độ thoái hóa của dầu bôi trơn và sự ăn mòn hóa học đối với các chi tiết cao su/kim loại màu trong hệ thống cấp nhiên liệu.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Mở rộng thử nghiệm hệ nhiên liệu DE10B5 trên động cơ diesel Common Rail điều khiển điện tử kết hợp hệ thống tuần hoàn khí xả nguội (Cooled EGR).
  • Nghiên cứu ứng dụng công nghệ nhũ tương nano (nano-emulsion) hoặc bổ sung các hạt phụ gia nano kim loại ($CeO_2, TiO_2, Al_2O_3$) nhằm tăng tốc độ hòa trộn và thúc đẩy quá trình cháy kiệt ở vùng phụ tải nhỏ.
  • Thực hiện đánh giá toàn diện chu trình vòng đời nhiên liệu (Life Cycle Assessment - LCA) từ khâu canh tác sắn/dầu cọ, chế biến công nghiệp đến phát thải đầu ống xả (Well-to-Wheel) để định lượng chính xác tiềm năng cắt giảm phát thải khí nhà kính ($CO_2$ tương đương).

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

  • Tác động học thuật: Cung cấp cơ sở dữ liệu thực nghiệm và lý thuyết hoàn chỉnh cho cộng đồng nghiên cứu động cơ đốt trong và nhiên liệu thay thế tại Việt Nam và khu vực ASEAN; mở ra hướng nghiên cứu mới về tối ưu hóa quá trình cháy tam phân trên các tạp chí khoa học chuyên ngành cơ khí động lực.
  • Tác động chuyển đổi ngành công nghiệp: Tháo gỡ nút thắt thị trường tiêu thụ cho 7 nhà máy sản xuất cồn sinh học trong nước; tạo cơ sở cho các tập đoàn xăng dầu (Petrolimex, PV OIL) xây dựng tiêu chuẩn cơ sở và quy chuẩn kỹ thuật thương mại hóa nhiên liệu diesel sinh học phối trộn cồn sắn DE10B5.
  • Ảnh hưởng chính sách quốc gia: Đóng góp luận cứ khoa học thực chứng để Bộ Công Thương và Chính phủ xem xét cập nhật lộ trình áp dụng tỷ lệ phối trộn nhiên liệu sinh học giai đoạn 2025 - 2030, mở rộng diện áp dụng sang phân khúc nhiên liệu cho động cơ diesel thay vì chỉ tập trung vào xăng sinh học E5/E10.
  • Lợi ích xã hội và môi trường: Giảm thiểu ô nhiễm khói đen và bụi mịn PM2.5 tại các đô thị lớn nơi tập trung mật độ cao xe tải nhẹ và xe khách; thúc đẩy chuỗi giá trị nông nghiệp bền vững, nâng cao thu nhập cho hàng trăm ngàn hộ nông dân trồng sắn tại các vùng trung du và miền núi.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận khung phân tích lý thuyết tích hợp giữa động học nhiệt động 1D, khảo sát quang học tia phun tốc độ cao và cơ chế phát thải hạt bụi; kế thừa các mô hình mô phỏng đã được hiệu chỉnh trên phần mềm AVL Boost.
  • Doanh nghiệp sản xuất và phân phối nhiên liệu: Các nhà máy sản xuất ethanol sinh học (Tùng Lâm, Dung Quất, Bình Phước) có cơ sở kỹ thuật để mở rộng phân khúc sản phẩm cồn nhiên liệu sang thị trường dầu diesel quy mô hàng triệu mét khối mỗi năm.
  • Các đơn vị vận tải và chủ phương tiện: Được sử dụng nguồn nhiên liệu sạch hơn, giảm phát thải bồ hóng gây ô uế buồng cháy, kéo dài chu kỳ bảo dưỡng vòi phun mà không phải tốn kém chi phí hoán cải kết cấu động cơ.
  • Cơ quan quản lý nhà nước: Bộ Giao thông Vận tải, Bộ Tài nguyên và Môi trường, Bộ Khoa học và Công nghệ có căn cứ thực nghiệm chính xác để ban hành các quy chuẩn kỹ thuật quốc gia (QCVN/TCVN) cho nhiên liệu diesel sinh học tam phâ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?

Đóng góp lý thuyết độc đáo nhất là việc xác lập cơ chế tương hỗ và bù trừ định lượng giữa các cấu tử phân cực và không phân cực trong hệ nhiên liệu tam phân DEB đối với quá trình cháy nén động cơ đốt trong. Nghiên cứu đã mở rộng lý thuyết quá trình cháy nén của Heywood và lý thuyết cháy hỗn hợp Sabathé bằng cách tích hợp mô hình giải phóng oxy nội tại phân tử vào quá trình cháy trễ. Luận án đã lượng hóa chính xác rằng mỗi $1%$ thể tích biodiesel dầu cọ bổ sung có khả năng trung hòa sự suy giảm trị số cetane gây ra bởi $0.93%$ ethanol sắn thông qua phương trình $CN_{ebdg} = CN_d - 0.59n_e + 0.55n_{bd}$, đồng thời duy trì độ nhớt và tính bôi trơn HFRR của hỗn hợp tương đương diesel gốc.

2. Tính đổi mới trong phương pháp nghiên cứu khi so sánh với các công trình trước đây?

So với nghiên cứu của Giakoumis (2011) chỉ tập trung vào thử nghiệm băng phanh phát thải thông thường và nghiên cứu của Qiang Fang (2013) chỉ khảo sát áp suất chỉ thị trên động cơ Common Rail, luận án đã thiết lập một phương pháp tiếp cận đa tầng khép kín: (1) Mô hình hóa nhiệt động 1D toàn diện trên AVL Boost với mô hình cháy cải biên AVL MCC; (2) Trực quan hóa cấu trúc tán sương tia phun thông qua camera quang học tốc độ cao Photron Fastcam SA3 ở tốc độ $3000\ FPS$; (3) Thử nghiệm thực tế toàn diện trên băng thử phanh điện APA 100 tích hợp phân tích phát thải khí chuyên sâu bằng AVL CEBII và phân tích phổ kích thước hạt bụi theo chuẩn PMP.

3. Phát hiện bất ngờ nhất trong dữ liệu thực nghiệm và lời giải thích lý thuyết?

Phát hiện bất ngờ nhất là sự suy giảm đột biến của phát thải bồ hóng ($soot$ giảm tới $52%$) và số lượng hạt mịn ($PN$ giảm $> 40%$) ở mẫu DE10B5 tại chế độ $100%$ tải mà không làm gia tăng phát thải $NO_x$ (chỉ dao động trong biên độ $\pm 3%$). Về mặt lý thuyết, việc đưa oxy liên kết vào buồng cháy thường kích hoạt phản ứng oxy hóa nitơ nhiệt Zeldovich làm tăng mạnh $NO_x$ (như Shi et al., 2006 từng ghi nhận mức tăng $32%$). Tuy nhiên, trong nghiên cứu này, nhiệt ẩn hóa hơi cực lớn của ethanol ($854\ kJ/kg$, gấp 2.8 lần diesel) đã hấp thụ nhiệt lượng cục bộ trong quá trình bay hơi nén, dập tắt các điểm nóng cục bộ (hot spots), tạo nên sự cân bằng hoàn hảo với hiệu ứng oxy hóa bồ hóng của biodiesel.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) hay không?

Có. Luận án cung cấp đầy đủ và chi tiết toàn bộ quy trình tái lập thực nghiệm bao gồm: thông số kích thước hình học động cơ D4BB catalog, cấu hình vòi phun cơ khí (áp suất mở $120\ bar$), tỷ lệ pha trộn thể tích chính xác của các mẫu nhiên liệu (DE5, DE10, DE5B5, DE10B5), quy chuẩn nhiên liệu đầu vào theo TCVN 5689-2018, các module phần tử thiết lập trong mô hình AVL Boost, cũng như sơ đồ nguyên lý bố trí đầu đo áp suất, cảm biến tốc độ, cân nhiên liệu 733S và đường ống lấy mẫu phân tích khí xả CEBII.

5. Định hướng lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình 10 năm được xác lập qua 3 giai đoạn:

  • Giai đoạn 1 (1 - 3 năm): Mở rộng khảo sát thực nghiệm DE10B5 trên động cơ diesel Common Rail thế hệ mới trang bị hệ thống luân hồi khí xả nguội Cooled EGR và bộ xử lý xúc tác SCR/DPF; tối ưu hóa bản đồ điều khiển phun đa điểm ($Pilot - Main - Post injection$).
  • Giai đoạn 2 (3 - 6 năm): Thử nghiệm độ bền mài mòn dài hạn $1000\ giờ$ trên đội xe thương mại thực tế; phát triển hệ phụ gia nano phân tán chống tách pha ở dải nhiệt độ âm ($T < 0^\circ C$).
  • Giai đoạn 3 (6 - 10 năm): Hoàn thiện chuỗi đánh giá vòng đời toàn diện Well-to-Wheel (LCA); phối hợp cùng các cơ quan quản lý ban hành Tiêu chuẩn Quốc gia cho nhiên liệu tam phân DEB và thương mại hóa rộng rãi trên toàn quốc.

Kết luận

  1. Xác lập cơ sở khoa học và thực tiễn vững chắc cho việc sử dụng hệ nhiên liệu tam phân diesel - ethanol - biodiesel (DEB) trên động cơ diesel thông thường, giải quyết bài toán tiêu thụ nguồn cồn sinh học sắn dư thừa ($> 400.000\ m^3/\text{năm}$) tại Việt Nam.
  2. Khẳng định công thức phối trộn tối ưu DE10B5 ($85%$ diesel khoáng + $10%$ ethanol sinh học + $5%$ biodiesel dầu cọ về thể tích) đảm bảo tính chất lý hóa gần tương đương diesel gốc (trị số cetane đạt $48.06$, độ nhớt $2.283\ cP$), duy trì độ bền phân tán pha ổn định trong điều kiện nhiệt độ phòng.
  3. Chứng minh bằng thực nghiệm rằng DE10B5 giúp cắt giảm đột phá từ $35%$ đến $52%$ phát thải muội than bồ hóng và trên $40%$ số lượng hạt bụi mịn ($PN$), giảm $15 - 28%$ phát thải $CO$, trong khi phát thải $NO_x$ không bị gia tăng nhờ cơ chế tự cân bằng nhiệt động giữa nhiệt ẩn hóa hơi cao của ethanol và oxy nội tại của biodiesel.
  4. Bảo toàn tính năng động lực học của động cơ D4BB với mức suy giảm mô-men cực đại không đáng kể ($< 2.4%$), suất tiêu hao nhiên liệu riêng chỉ tăng nhẹ ($3.2 - 4.8%$) hoàn toàn tương thích với mức chênh lệch nhiệt trị, đồng thời nâng cao hiệu suất nhiệt có ích $\eta_e$ thêm $1.1 - 2.3%$.
  5. Xây dựng và kiểm chứng thành công mô hình mô phỏng nhiệt động 1D trên phần mềm AVL Boost kết hợp mô hình cháy AVL MCC, đạt độ chính xác cao với sai số dự báo so với thực nghiệm dưới $5%$ cho các chỉ tiêu tính năng kỹ thuật.
  6. Mở ra ba nhánh nghiên cứu đột phá mới cho ngành kỹ thuật động lực: tối ưu hóa chiến lược phun đa giai đoạn cho nhiên liệu tam phân trên động cơ Common Rail, ứng dụng công nghệ nhũ tương nano kim loại đa cấu tử, và xây dựng khung chính sách tiêu chuẩn hóa nhiên liệu sinh học tam phân phục vụ phát triển giao thông xanh bền vững.