Tổng quan về luận án
Nghiên cứu của nghiên cứu sinh Nguyễn Đức Hạnh với đề tài "Nghiên cứu ảnh hưởng của hỗn hợp nhiên liệu dầu thực vật - diesel đến phun nhiên liệu, tạo hỗn hợp, cháy và tính năng của động cơ diesel tàu thủy" (Luận án Tiến sĩ Kỹ thuật, Chuyên ngành Khai thác và bảo trì tàu thủy, Mã số: 95.116, Trường Đại học Hàng hải Việt Nam, năm 2020) được thực hiện trong bối cảnh ngành hàng hải quốc tế đối mặt với những tiêu chuẩn môi trường nghiêm ngặt chưa từng có. Công ước quốc tế MARPOL 73/78 (Phụ lục VI) của Tổ chức Hàng hải Quốc tế (IMO) đặt ra các quy định bắt buộc về giới hạn phát thải ôxít nitơ ($NO_x$ Tier I, II, III) và giới hạn hàm lượng lưu huỳnh trong nhiên liệu hàng hải không vượt quá 0,5% áp dụng trên phạm vi toàn cầu từ ngày 01/01/2020. Đối với đội tàu vận tải biển Việt Nam có tổng trọng tải trên 8 triệu tấn DWT, hàng năm tiêu thụ khoảng 400 triệu tấn nhiên liệu hóa thạch, việc chuyển đổi sang các nguồn năng lượng thân thiện với môi trường là yêu cầu cấp bách sống còn nhằm duy trì quyền cập cảng quốc tế và tối ưu hóa hiệu quả khai thác kinh tế.
Khoảng trống nghiên cứu (Research Gap) then chốt mà luận án nhận diện và tập trung giải quyết là sự thiếu hụt cơ sở khoa học mang tính hệ thống về động lực học tia phun, quá trình phân rã chất lỏng, động học hòa trộn - cháy và phương pháp hiệu chỉnh thông số vận hành cho hệ thống phun nhiên liệu (HTPNL) kiểu cơ khí truyền thống trên động cơ diesel tàu thủy cỡ vừa và lớn khi sử dụng trực tiếp dầu thực vật (Straight Vegetable Oil - SVO) pha trộn với dầu diesel (DO). Phần lớn các công trình trước đây chỉ khảo sát trên động cơ ô tô cao tốc hoặc nghiên cứu nhiên liệu diesel sinh học (Biodiesel/FAME) - vốn đòi hỏi quy trình este hóa phức tạp, chi phí sản xuất cao và sử dụng nhiều hóa chất xúc tác ($NaOH$) gây ô nhiễm thứ cấp, đồng thời có thời gian bảo quản ngắn (dưới 1 năm), không phù hợp với đặc thù hành trình dài ngày của tàu biển.
Luận án xác lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học cụ thể:
- Câu hỏi nghiên cứu 1 (RQ1): Sự khác biệt về các đặc tính hóa lý (độ nhớt động học, khối lượng riêng, sức căng bề mặt, trị số Cetan) của hỗn hợp dầu cọ - diesel (PO) ảnh hưởng như thế nào đến đặc tính vĩ mô (chiều dài tia phun $S$, chiều dài phân rã sơ cấp $L_b$, góc nón tia phun $\theta_s$) và vi mô (đường kính trung bình hạt Sauter $SMD$) của chùm tia phun?
- Câu hỏi nghiên cứu 2 (RQ2): Quy luật biến thiên của thời điểm bắt đầu phun, thời gian cháy trễ ($\tau_i$), diễn biến áp suất cháy trong xy lanh ($p_z$) và trường nhiệt độ - vận tốc cháy diễn ra như thế nào trong không gian buồng đốt động cơ chính tàu thủy Hanshin 6LU32?
- Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để thiết lập một mô hình toán học tối ưu hóa đa mục tiêu nhằm đồng thời giảm thiểu suất tiêu hao nhiên liệu ($g_e$) và nồng độ phát thải $NO_x$ thông qua việc hiệu chỉnh góc phun sớm (GPS) và áp suất phun (ASP)?
Hệ thống giả thuyết nghiên cứu:
- Giả thuyết 1 (H1): Khi tỉ lệ dầu cọ tăng từ 10% đến 30% (PO10, PO20, PO30), độ nhớt cao và sức căng bề mặt lớn sẽ làm tăng chiều dài phân rã sơ cấp $L_b$, tăng độ xâm nhập tia phun $S$, nhưng làm giảm góc tán sương $\theta_s$ và tăng đường kính hạt $SMD$, cản trở quá trình xé tơi nhiên liệu nếu giữ nguyên kết cấu HTPNL ban đầu.
- Giả thuyết 2 (H2): Trị số Cetan cao hơn của dầu cọ (dầu cọ nguyên chất đạt $CN = 52,92$ so với $42,89$ của DO) sẽ bù trừ một phần độ nhớt bằng cách rút ngắn thời gian cháy trễ $\tau_i$, dịch chuyển thời điểm bắt đầu cháy (SOC) về sớm hơn, dẫn đến sự gia tăng áp suất cháy cực đại $p_z$ ở các góc quay trục khuỷu (GQTK) trước điểm chết trên (ĐCT).
- Giả thuyết 3 (H3): Tồn tại một miền giá trị tối ưu cục bộ của bộ ba thông số (tỉ lệ phối trộn dầu cọ, góc phun sớm, áp suất nâng kim phun) thỏa mãn đồng thời tiêu chuẩn phát thải $NO_x$ theo MARPOL 73/78 và duy trì suất tiêu hao nhiên liệu có ích $g_e$ ở mức kinh tế kỹ thuật tối ưu.
Khung lý thuyết của công trình tích hợp lý thuyết tia phun chất lỏng trong môi trường nén ép (mô hình Dent và Hiroyasu), lý thuyết động lực học lưu chất tính toán (CFD) dựa trên phương trình Navier-Stokes và mô hình rối $k-\epsilon$, cùng lý thuyết quy hoạch thực nghiệm Bề mặt đáp ứng (Response Surface Methodology - RSM). Luận án được thực hiện trong khuôn khổ Đề tài NCKH & PTCN cấp Quốc gia mã số ĐT.11/NLSH thuộc “Đề án phát triển NLSH đến năm 2015, tầm nhìn đến năm 2025” do PGS. Đặng Văn Uy chủ nhiệm. Đối tượng thực nghiệm trực tiếp là động cơ diesel chính tàu thủy Hanshin 6LU32 tại Trung tâm thí nghiệm hệ động lực tàu thủy - Trường Đại học Hàng hải Việt Nam, thử nghiệm ở hai chế độ phụ tải 400 kW (40% tải) và 600 kW (60% tải), tốc độ định mức 273 vòng/phút, khảo sát toàn diện 5 loại nhiên liệu: DO, PO10, PO20, PO30 và PO100 có pha phụ gia "Nano Fuel Booster" tỉ lệ 1/800.
+----------------------------------------------------------------------------------------------------+
| KHUNG PHÂN TÍCH TỔNG QUAN |
+----------------------------------------------------------------------------------------------------+
| [Đặc tính nhiên liệu PO] --> [Động lực học tia phun] --> [Quá trình cháy] --> [Chỉ tiêu] |
| - Độ nhớt: 2,6 -> 6,45 cSt - Tăng chiều dài tia S - Rút ngắn cháy - Suất tiêu |
| - Tỉ trọng: 850 -> 866,8 kg/m3 - Tăng phân rã sơ cấp Lb trễ \tau_i hao ge |
| - Cetan: 42,89 -> 52,11 - Tăng kích thước hạt SMD - Dịch chuyển áp - Phát thải |
| - Sức căng: 20 -> 24,2 mN/m - Giảm góc nón \theta_s suất đỉnh p_z NOx (IMO) |
+----------------------------------------------------------------------------------------------------+
|
v
+-------------------------------------------------------+
| TỐI ƯU HÓA ĐA MỤC TIÊU (RSM - BỀ MẶT ĐÁP ỨNG) |
| Hiệu chỉnh: Góc phun sớm (GPS) + Áp suất phun (ASP) |
+-------------------------------------------------------+
Literature Review và Positioning
Tổng quan y văn quốc tế và trong nước về ứng dụng nhiên liệu sinh học cho động cơ diesel phân hóa thành hai trường phái kỹ thuật chính:
Trường phái thứ nhất tập trung vào công nghệ biến tính este hóa dầu thực vật thành Biodiesel (FAME). Các công trình của Ayhan Demirbas (2007), Nguyễn Hoàng Vũ (2009 với đề tài ĐT.12/NLSH cho xe quân sự), Viện Hóa học Công nghiệp Việt Nam (2009 với B5 từ mỡ cá), Phan Đắc Yến (2015) và Phan Văn Quân (2015) đã chứng minh tính tương thích tốt của Biodiesel với hệ thống nhiên liệu tiêu chuẩn. Tuy nhiên, rào cản lớn nhất của trường phái này trong vận tải biển là chi phí sản xuất cao, tiêu hao nhiều năng lượng và hóa chất xúc tác kiềm, cùng nguy cơ suy giảm chất lượng nhanh chóng do hiện tượng oxy hóa khi lưu trữ trong các két nhiên liệu tàu thủy hoạt động dài ngày trên biển.
Trường phái thứ hai hướng đến giải pháp sử dụng trực tiếp dầu thực vật (SVO) hoặc pha trộn cơ học với dầu diesel. Điển hình là các phát triển của Hãng Elsbett (Đức) trên động cơ ô tô chuyên dụng, nghiên cứu của Nguyễn Thạch (2009) về bộ chuyển đổi SVO cho máy phát điện dưới 100 HP, công trình của Trần Thế Nam (2015) về thiết bị hòa trộn liên tục trên động cơ 6LU32, Phùng Minh Lộc (2014) khảo sát dầu dừa trên tàu cá cỡ nhỏ, và Mai Đức Nghĩa (2014) tối ưu hóa góc phun sớm cho hỗn hợp B15 trên động cơ cao tốc.
+--------------------------------------------------------+
| CÁC NHÁNH NGHIÊN CỨU CHÍNH |
+--------------------------------------------------------+
/ \
/ \
[Trường phái Biodiesel (FAME)] [Trường phái Dầu thực vật trực tiếp (SVO)]
- Demirbas (2007) - Elsbett (Đức)
- Nguyễn Hoàng Vũ (2009) - Nguyễn Thạch (2009)
- Phan Đắc Yến (2015) - Trần Thế Nam (2015)
- Phan Văn Quân (2015) - Nguyễn Đức Hạnh (2020 - Luận án)
* Hạn chế: Giá thành cao, dễ biến chất * Ưu điểm: Giá rẻ, tận dụng hệ thống hâm FO
Trong dòng nghiên cứu SVO, tồn tại các cuộc tranh luận học thuật sâu sắc về tương quan giữa chất lượng hòa trộn, tính kinh tế và phát thải khí xả:
- Quan điểm của Herchel Thaddeus và cộng sự (2001) cùng Kalam & Masjuki / A. (2003) tại Đại học Malaya chỉ ra rằng dầu cọ nguyên chất và hỗn hợp dầu cọ - diesel làm giảm phát thải bồ hóng đáng kể nhờ cấu trúc phân tử giàu oxy (khoảng 16% khối lượng oxy) và tỉ lệ C/H thấp. Tuy nhiên, các tác giả này nhấn mạnh quá trình tạo hỗn hợp bị suy giảm nghiêm trọng do độ nhớt cao, dẫn đến tăng suất tiêu hao nhiên liệu và suy giảm công suất nếu không can thiệp điều chỉnh thông số phun.
- Ngược lại, nhóm nghiên cứu của Venkataramana Reddy và cộng sự (2010) với dầu Honge (H20), Siddalingappa (2011) với dầu Karanja (K15), Purushothaman và cộng sự (2009) với tinh dầu cam (OSP30), và Sapit Azwan và cộng sự (2011) chứng minh rằng việc tăng áp suất phun từ 180-200 bar lên 210-235 bar và tăng góc phun sớm thêm 4-5° GQTK có thể phục hồi hoàn toàn hiệu suất nhiệt tương đương dầu DO, đồng thời cắt giảm phát thải HC và CO.
So sánh với các nghiên cứu quốc tế tiêu biểu:
- So với nghiên cứu của Tamilporai và Mohd. Shabir (2008) trên động cơ 55 kW dùng dầu Jatropha (chỉ khảo sát mô phỏng diễn biến áp suất mà chưa phân tích cấu trúc vi mô tia phun 10 lỗ), luận án của NCS. Nguyễn Đức Hạnh đã giải quyết triệt để bài toán dòng chảy nhiều pha qua vòi phun và cấu trúc hình học tia phun thực tế.
- So với công trình của Naga Sarada và cộng sự (2010) về dầu bông (chỉ tối ưu đơn biến áp suất phun trên động cơ tĩnh tại), luận án đã nâng cấp phương pháp luận lên quy hoạch thực nghiệm đa biến trên động cơ đẩy chính tàu thủy công suất thực tế 600 kW.
Định vị học thuật của luận án: Luận án đã lấp đầy khoảng trống nghiên cứu bằng cách kết hợp cơ chế gia nhiệt tận dụng hệ thống hâm dầu nặng (HFO/FO) có sẵn trên tàu biển (để hạ độ nhớt dầu cọ từ 40,24 cSt ở 40°C xuống mức tương thích) với việc xác lập phương pháp hiệu chỉnh đồng thời góc phun sớm và áp suất phun, tạo bước tiến vững chắc về lý thuyết và thực tiễn cho việc thay thế nhiên liệu hóa thạch trong ngành vận tải biển.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án đã đóng góp những luận điểm khoa học quan trọng vào lý thuyết động cơ đốt trong và cơ học chất lưu:
- Mở rộng lý thuyết tia phun nhiên liệu của Dent (1971) và Hiroyasu (1980): Thiết lập mô hình giải tích mô tả quá trình phân rã chất lỏng phi tiêu chuẩn có độ nhớt và sức căng bề mặt biến thiên. Luận án chứng minh rằng đối với nhiên liệu PO, chiều dài phân rã sơ cấp $L_b$ tăng tỉ lệ thuận với độ nhớt động học $\mu_l$ và sức căng bề mặt $\sigma$:
$$L_b = C \cdot d_0 \cdot \left(\frac{\rho_l}{\rho_g}\right)^{0.5} \cdot Re_l^a \cdot We_l^b$$
Đồng thời, đường kính trung bình Sauter ($SMD$) tăng theo hàm số của hệ số Reynolds ($Re_l$) và Weber ($We_l$), giải thích chính xác hiện tượng hạt sương nhiên liệu có kích thước lớn hơn và góc nón tia phun $\theta_s$ bị thu hẹp khi nồng độ dầu cọ tăng.
- Bổ sung quy luật động học cháy trễ cho nhiên liệu sinh học gốc este/triglyceride: Luận án làm rõ tác động đối nghịch giữa hai yếu tố: sự chậm trễ trong quá trình bay hơi vật lý (do độ nhớt và nhiệt độ sôi cao) và sự gia tốc phản ứng hóa học tự bốc cháy (do trị số Cetan của dầu cọ $CN = 52,92$ cao hơn DO $CN = 42,89$). Kết quả chỉ ra rằng yếu tố động học hóa học chiếm ưu thế, làm rút ngắn thời gian cháy trễ tổng thể $\tau_i$ từ 0,15 đến 0,35 ms tùy theo chế độ tải.
- Chuyển dịch mô hình đánh giá chu trình công tác: Chuyển từ quan điểm hiệu chỉnh đơn biến truyền thống sang lý thuyết cân bằng động giữa động năng tia phun, phân bố năng lượng rối và tốc độ tỏa nhiệt cục bộ nhằm khống chế cơ chế hình thành $NO_x$ theo thuyết Zeldovich nhiệt độ cao mà không làm suy giảm hiệu suất kinh tế có ích $\eta_e$.
+------------------------------------------------------------------------------------------------+
| KHUNG MÔ HÌNH LÝ THUYẾT MỞ RỘNG |
+------------------------------------------------------------------------------------------------+
| 1. ĐỘNG LỰC HỌC TIA PHUN (Dent - Hiroyasu mở rộng) |
| - Phân rã sơ cấp: L_b = f(d_0, \rho_l/\rho_g, Re_l, We_l) |
| - Kích thước hạt: SMD = f(\mu_l, \sigma, \Delta p) |
| |
| 2. ĐỘNG HỌC CHÁY VÀ PHÁT THẢI (Zeldovich - Arrhenius) |
| - Cháy trễ: \tau_i = f(p_c, T_c, E_A, CN_{mix}) --> Rút ngắn thời gian cháy trễ |
| - Nhiệt độ cháy cục bộ T_max & Nồng độ O2 nội tại (16%) --> Kiểm soát phát thải NOx |
| |
| 3. TỐI ƯU HÓA BỀ MẶT ĐÁP ỨNG (Response Surface Methodology) |
| - Hàm mục tiêu kép: Min(g_e) & Min(NOx) = f(GPS, ASP, % Dầu cọ) |
+------------------------------------------------------------------------------------------------+
Khung phân tích độc đáo
Khung phân tích của luận án là sự hợp nhất của ba nền tảng lý thuyết:
- Lý thuyết thủy động lực học chất lưu (Hydrodynamics): Vận dụng hệ phương trình Navier-Stokes trung bình Reynolds (RANS) kết hợp mô hình rối hai phương trình $k-\epsilon$ ($C_\mu = 0,09; C_\epsilon = 1,92$) để mô phỏng dòng chảy đa pha qua chùm 10 lỗ phun có đường kính $d_0$ của vòi phun kín động cơ Hanshin 6LU32.
- Lý thuyết truyền chất và phản ứng cháy (Combustion Kinetics): Mô hình hóa quá trình oxy hóa hỗn hợp đa thành phần trong buồng đốt thực tế thông qua phần mềm Ansys Fluent, phân tích sự phân bố năng lượng động năng rối, đường đồng mức vận tốc và trường nhiệt độ theo từng độ GQTK.
- Lý thuyết tối ưu hóa bề mặt đáp ứng (Response Surface Methodology - RSM): Xây dựng hàm mục tiêu hai biến đáp ứng ${g_e, NO_x}$ phụ thuộc vào ba biến điều khiển độc lập: Tỉ lệ pha trộn dầu cọ ($X_1 \in [10%, 30%]$), Góc phun sớm ($X_2 \in [14^\circ, 22^\circ]$ trước ĐCT) và Áp suất nâng kim phun ($X_3 \in [180, 240\text{ bar}]$).
Điều kiện biên xác lập chặt chẽ: Áp suất không khí nạp $p_n$, nhiệt độ nạp $T_n$, tỷ số nén $\epsilon$, hình học đỉnh piston và nắp xi lanh được số hóa chính xác bằng phần mềm AutoCAD Mechanical 2017 trước khi chuyển đổi sang lưới tính toán (Gambit/Grid).
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu tuân thủ triết lý thực chứng (Positivism), kết hợp chặt chẽ giữa tính toán giải tích lý thuyết, mô phỏng số học CFD và kiểm chứng thực nghiệm trên đối tượng quy mô công nghiệp (động cơ Hanshin 6LU32). Thiết kế nghiên cứu đa cấp độ được thiết lập tuần tự:
- Cấp độ vi mô: Đo lường các thuộc tính hóa lý cơ bản của 5 mẫu nhiên liệu (DO, PO10, PO20, PO30, PO100).
- Cấp độ trung mô (Mesoscopic): Mô phỏng 3D trường phân bố áp suất, vận tốc dòng chảy qua 10 lỗ vòi phun và không gian buồng đốt bằng Ansys Fluent.
- Cấp độ vĩ mô: Thử nghiệm tải thực tế trên băng thử động lực học, đo áp suất chỉ thị trong xy lanh, chụp ảnh quang học quá trình tạo hỗn hợp - cháy và đo phổ phát thải khí xả.
+----------------------------------------------------------------------------------------------------+
| THIẾT KẾ NGHIÊN CỨU ĐA CẤP ĐỘ |
+----------------------------------------------------------------------------------------------------+
| [Cấp độ Vi mô] --> Phân tích thuộc tính hóa lý 5 mẫu nhiên liệu (ASTM standard) |
| | |
| [Cấp độ Trung mô] --> Mô phỏng CFD 3D trường dòng 10 lỗ vòi phun & buồng đốt (Ansys Fluent) |
| | |
| [Cấp độ Vĩ mô] --> Thử nghiệm thực tế trên động cơ Hanshin 6LU32 (400 kW & 600 kW) |
| | |
| [Tối ưu hóa] --> Quy hoạch thực nghiệm bề mặt đáp ứng RSM (Box-Behnken matrix) |
+----------------------------------------------------------------------------------------------------+
Quy trình nghiên cứu rigorous
Quy trình nghiên cứu bao gồm các giao thức kỹ thuật chuẩn hóa:
- Chuẩn hóa nhiên liệu: Các mẫu nhiên liệu PO được phân tích tại Trung tâm Đào tạo và Tư vấn KH&CN Bảo vệ Môi trường Thủy theo Quy chuẩn Kỹ thuật Quốc gia QCVN: 2015/BKH&CN và tiêu chuẩn ASTM quốc tế: Khối lượng riêng (ASTM D1298), Độ nhớt động học (ASTM D445), Trị số Cetan (ASTM D613), Nhiệt trị (ASTM D240), Sức căng bề mặt tại 40°C và Áp suất hơi bão hòa tại 25°C.
- Mô phỏng động lực học: Lưới tính toán 3D được thiết lập trong Gambit với mật độ lưới dày đặc tại vùng lỗ vòi phun và bề mặt đỉnh piston, kiểm tra tính độc lập của lưới nhằm đảm bảo sai số tính toán dưới 2%.
- Hệ thống đo lường thực nghiệm: Bệ thử động cơ diesel Hanshin 6LU32 kết nối với phanh thủy lực chuyên dụng; cảm biến áp suất áp điện Kistler lắp tại xy lanh số 1 để đo áp suất cháy cực đại $p_z$; thiết bị nội soi quang học Visio Scope ghi hình ngọn lửa buồng đốt qua cửa sổ thạch anh; thiết bị phân tích khí xả chuyên dụng đo nồng độ $NO_x$ (ppm và quy đổi g/kWh).
+-------------------------------------------------------+
| SƠ ĐỒ HỆ THỐNG THỬ NGHIỆM THỰC TẾ |
+-------------------------------------------------------+
|
+---------------------------+---------------------------+
| |
v v
[Động cơ Hanshin 6LU32] [Hệ thống đo lường]
- Công suất: 400 kW / 600 kW - Cảm biến áp suất Kistler (p_z)
- Vòng quay: 273 v/p - Nội soi Visio Scope (-11,7° -> 22°)
- Vòi phun kín 10 lỗ - Phân tích khí xả NOx
- Két hòa trộn & hâm dầu - Phanh thủy lực đo công suất
- Tam giác hóa dữ liệu (Triangulation): Kết hợp đồng thời ba nguồn dữ liệu: kết quả giải tích từ công thức thực nghiệm, trường thông số mô phỏng CFD từ Ansys Fluent và số liệu thực nghiệm trên động cơ.
Data và phân tích
Số liệu phân tích thuộc tính nhiên liệu (Trích xuất từ Bảng 1.5 của tài liệu gốc):
| STT |
Tên chỉ tiêu |
Đơn vị |
DO |
PO10 |
PO20 |
PO30 |
PO100 |
| 1 |
Khối lượng riêng ở 15°C |
$\text{kg/m}^3$ |
850,0 |
853,8 |
859,9 |
866,8 |
922,5 |
| 2 |
Độ nhớt động học ở 40°C |
$\text{cSt}$ |
2,60 |
3,42 |
5,31 |
6,45 |
40,24 |
| 3 |
Trị số Cetan (CN) |
- |
42,89 |
50,13 |
50,91 |
52,11 |
52,92 |
| 4 |
Nhiệt độ chớp cháy cốc kín |
°C |
72 |
73 |
75 |
77 |
135 |
| 5 |
Nhiệt trị |
$\text{MJ/kg}$ |
43,40 |
39,72 |
39,55 |
38,69 |
37,11 |
| 6 |
Sức căng bề mặt ở 40°C |
$\text{mN/m}$ |
20,0 |
21,4 |
22,8 |
24,2 |
34,0 |
| 7 |
Áp suất hơi bão hòa ở 25°C |
$\text{MPa}$ |
0,0477 |
0,0456 |
0,0432 |
0,0407 |
0,0500 |
Xử lý số liệu và độ tin cậy: Sử dụng phần mềm Tecplot để trực quan hóa trường vector vận tốc và biểu đồ áp suất theo góc quay trục khuỷu. Phân tích phương sai (ANOVA) trong quy hoạch thực nghiệm được ứng dụng để kiểm định mức ý nghĩa thống kê của các hệ số hồi quy ($p < 0,05$), xây dựng biểu đồ Pareto chuẩn hóa và phân tích số dư nhằm bảo đảm mô hình tối ưu không bị chệch.
Phát hiện đột phá và implications
Những phát hiện then chốt
- Quy luật biến thiên cấu trúc tia phun: Khi tỷ lệ dầu cọ tăng từ PO10 đến PO30, chiều dài phân rã sơ cấp $L_b$ tăng từ 8,5% đến 24,3% so với dầu DO; chiều dài tia phun $S$ tăng thêm do động lượng hạt chất lỏng lớn (khối lượng riêng đạt 866,8 $\text{kg/m}^3$ ở PO30), nhưng góc nón tia phun $\theta_s$ suy giảm rõ rệt. Đường kính trung bình hạt Sauter ($SMD$) tăng từ 18,2 $\mu\text{m}$ (DO) lên mức trên 24,5 $\mu\text{m}$ (PO30), khẳng định độ tơi sương bị suy giảm nếu không tăng áp suất phun.
- Hiện tượng rút ngắn thời gian cháy trễ ($\tau_i$): Trái với quan niệm thông thường cho rằng dầu thực vật khó cháy sẽ làm kéo dài pha cháy trễ, dữ liệu thực nghiệm trên động cơ 6LU32 ở cả hai chế độ 400 kW và 600 kW chứng minh $\tau_i$ của PO10, PO20, PO30 đều ngắn hơn so với dầu DO từ 0,5° đến 1,8° GQTK. Nguyên nhân vật lý là do trị số Cetan vượt trội của dầu cọ ($CN = 50,13 - 52,11$ so với $42,89$ của DO) thúc đẩy phản ứng oxy hóa tự kích hoạt nhiệt độ cao diễn ra sớm hơn.
- Diễn biến áp suất cháy cực đại ($p_z$) và thị giác hóa buồng đốt: Hình ảnh nội soi Visio Scope tại các mốc $-11,7^\circ, -7^\circ, -3^\circ, +5^\circ, +11^\circ, +22^\circ$ GQTK cho thấy ngọn lửa bùng phát sớm hơn ở vùng gần đỉnh ĐCT. Đồ thị áp suất chỉ thị xác nhận đỉnh $p_z$ của hỗn hợp PO xuất hiện sớm hơn 1 - 2° GQTK so với DO khi giữ nguyên góc phun sớm định mức ban đầu, gây ra hiện tượng tăng áp suất cản ngược trong hành trình nén cuối của piston.
- Nghịch lý phát thải $NO_x$ và tiêu hao nhiên liệu ($g_e$): Do nhiệt trị của PO thấp hơn DO (PO30 đạt 38,69 MJ/kg so với 43,40 MJ/kg của DO), suất tiêu thụ nhiên liệu $g_e$ thực tế tăng từ 4,2% đến 9,8% để duy trì cùng mức công suất kéo phanh. Nồng độ phát thải $NO_x$ có xu hướng gia tăng cục bộ do hàm lượng oxy nội tại trong phân tử dầu cọ (chiếm 16% khối lượng) thúc đẩy phản ứng nhiệt Zeldovich tại vùng hòa khí đồng nhất.
- Hiệu chỉnh tối ưu hóa đa mục tiêu bằng Bề mặt đáp ứng (RSM): Luận án đã xác định thành công bộ thông số vận hành tối ưu cho động cơ 6LU32: Sử dụng hỗn hợp PO20 kết hợp với việc giảm góc phun sớm (GPS) lùi lại 1,5° - 2,0° GQTK và tăng áp suất nâng kim phun lên mức 215 - 220 bar. Giải pháp này giúp phục hồi chất lượng phun sương hạt nhỏ ($SMD \approx 19,5\ \mu\text{m}$), hạ đỉnh nhiệt độ cục bộ trong buồng đốt, qua đó giảm nồng độ $NO_x$ xuống dưới ngưỡng giới hạn của IMO MARPOL Annex VI Tier II mà vẫn duy trì suất tiêu hao nhiên liệu ở mức kinh tế kỹ thuật chấp nhận được.
+----------------------------------------------------------------------------------------------------+
| TÓM TẮT PHÁT HIỆN THỰC NGHIỆM |
+----------------------------------------------------------------------------------------------------+
| Chỉ tiêu Nhiên liệu DO Nhiên liệu PO20 (chưa hiệu chỉnh) PO20 (đã hiệu chỉnh tối ưu)
| --------------------------------------------------------------------------------------------------
| Độ nhớt động học (40°C) 2,60 cSt 5,31 cSt 5,31 cSt (hâm nóng ~40°C)
| Trị số Cetan 42,89 50,91 50,91
| Góc phun sớm (GPS) Định mức Định mức (quá sớm với PO) Lùi 1,5° - 2,0° GQTK
| Áp suất phun (ASP) Định mức Định mức (hạt sương to) Tăng lên 215 - 220 bar
| Thời gian cháy trễ (\tau_i) Chuẩn Rút ngắn (bốc cháy sớm) Tối ưu hóa gần ĐCT
| Đường kính hạt (SMD) 18,2 \mu m 22,8 \mu m 19,5 \mu m
| Phát thải NOx Tiêu chuẩn IMO Tăng cục bộ Đạt chuẩn MARPOL Tier II
| Suất tiêu hao (g_e) Chuẩn (100%) Tăng (+6,5%) Kiểm soát tốt (+3,2%)
+----------------------------------------------------------------------------------------------------+
Implications đa chiều
- Về mặt lý thuyết: Cung cấp mô hình định lượng liên kết giữa tính chất nhiệt động của nhiên liệu tái tạo với quá trình phát triển tia phun và động lực học buồng đốt kín, đóng góp luận cứ khoa học cho ngành động cơ nhiệt tàu thủy.
- Về mặt phương pháp luận: Thiết lập quy trình chuẩn hóa kết hợp giữa phần mềm mô phỏng CFD 3D và thuật toán quy hoạch thực nghiệm Bề mặt đáp ứng (RSM) để giải bài toán cải tiến hệ thống cơ khí cũ sang sử dụng nhiên liệu mới.
- Về mặt ứng dụng thực tiễn: Cho phép các chủ tàu, công ty vận tải biển tại Việt Nam tận dụng ngay đội tàu hiện hữu và hệ thống két hâm sấy FO sẵn có để chuyển sang sử dụng nhiên liệu sinh học PO20 mà không cần đầu tư thay mới động cơ tốn kém hàng triệu USD.
- Về mặt chính sách: Cung cấp cơ sở khoa học xác đáng cho Bộ Giao thông Vận tải và Bộ Khoa học & Công nghệ trong việc xây dựng lộ trình quy chuẩn hóa việc phối trộn nhiên liệu sinh học cho vận tải thủy nội địa và ven biển theo Quyết định 177/2007/QĐ-TTg và Quyết định 53/2012/QĐ-TTg của Thủ tướng Chính phủ.
Limitations và Future Research
Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:
- Phạm vi thử nghiệm động cơ: Nghiên cứu thực nghiệm tập trung chuyên sâu trên dòng động cơ diesel 4 kỳ trung/thấp tốc Hanshin 6LU32; chưa mở rộng kiểm chứng trên các dòng động cơ 2 kỳ siêu lớn (Crosshead marine engines) lai chân vịt trực tiếp vốn có chu trình công tác và cơ chế quét khí khác biệt.
- Giới hạn dải tải thử nghiệm: Dữ liệu thực nghiệm mới giới hạn ở hai mức tải đại diện là 40% (400 kW) và 60% (600 kW) tại bệ thử phòng thí nghiệm, chưa khảo sát đầy đủ dải tải 85% - 100% (Continuous Maximum Rating - MCR) khi tàu hành hải trong điều kiện bão gió khắc nghiệt.
- Độ bền dài hạn và suy thoái chi tiết: Luận án chưa có điều kiện theo dõi thử nghiệm liên tục trong 3.000 - 5.000 giờ hoạt động để đánh giá mức độ mài mòn bơm cao áp, hiện tượng lắng cặn carbon tại kim phun, và tốc độ suy thoái của dầu bôi trơn các-te.
Chương trình nghiên cứu tương lai cần tập trung vào các định hướng:
- Mở rộng nghiên cứu sang hệ thống phun nhiên liệu điện tử Common Rail Diesel Injection (CRDi) điều khiển áp suất và thời điểm phun biến thiên linh hoạt theo tải.
- Nghiên cứu độ bền ăn mòn hóa học và tích tụ muội than trên xéc-măng và bề mặt nắp xy lanh sau chu kỳ vận hành dài hạn 10.000 giờ.
- Khảo sát các nguồn dầu thực vật thế hệ mới phi thực phẩm như dầu Jatropha, vi tảo hoặc dầu mỡ cá phế thải pha trộn đa hợp phần.
- Phát triển hệ thống điều khiển tự động tỷ lệ hòa trộn nhiên liệu liên tục theo thời gian thực (Real-time Adaptive Fuel Blending System) tích hợp cảm biến đo $NO_x$ trực tiếp tại đường xả tàu thủy.
Tác động và ảnh hưởng
Nghiên cứu mang lại những giá trị tác động sâu rộng:
- Tác động học thuật: Mở ra hướng nghiên cứu liên ngành chuyên sâu giữa Kỹ thuật cơ khí động lực, Thủy khí động lực học và Hóa học năng lượng sinh học. Dự kiến tạo tiền đề cho nhiều ấn phẩm khoa học trên các tạp chí quốc tế chuyên ngành (như Applied Energy, Fuel, Journal of Marine Science and Technology).
- Chuyển đổi công nghiệp hàng hải: Cung cấp giải pháp kỹ thuật có tính khả thi cực cao cho hơn 1.500 tàu biển Việt Nam đang sử dụng động cơ cơ khí thế hệ cũ, giúp các doanh nghiệp vận tải biển tiết kiệm 15 - 20% chi phí chuyển đổi công nghệ so với việc lắp đặt tháp lọc khí thải (Scrubber) hoặc hoán cải sang động cơ LNG đắt đỏ.
- Tác động chính sách và xã hội: Hỗ trợ hiện thực hóa Chiến lược Phát triển Bền vững Kinh tế Biển Việt Nam đến năm 2030, tầm nhìn 2045, giảm thiểu hàng trăm ngàn tấn phát thải $SO_x, CO_2$ và khói đen mỗi năm tại các vùng duyên hải và cảng biển quốc gia.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận nghiên cứu kết hợp giữa mô phỏng số học Ansys Fluent và kỹ thuật chẩn đoán quang học Visio Scope trong buồng đốt kín thực tế.
- Kỹ sư khai thác và Máy trưởng tàu biển: Nắm bắt quy trình vận hành và thông số kỹ thuật chuẩn mực để hiệu chỉnh vòi phun, bơm cao áp và hệ thống hâm dầu khi tiếp nhận nhiên liệu sinh học PO tại các cảng.
- Doanh nghiệp vận tải biển và Đóng tàu: Sở hữu cẩm nang giải pháp kỹ thuật chi phí thấp nhằm thỏa mãn ngay các tiêu chuẩn MARPOL Phụ lục VI mà không làm gián đoạn kế hoạch khai thác thương mại của đội tàu.
- Cơ quan Quản lý Nhà nước (Cục Hàng hải Việt Nam, Cục Đăng kiểm Việt Nam): Có cơ sở dữ liệu thực nghiệm chuẩn xác để ban hành các quy chuẩn kỹ thuật quốc gia về đăng kiểm an toàn kỹ thuật cho tàu biển sử dụng nhiên liệu sinh học hỗn hợp.
Câu hỏi chuyên sâu
1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?
Luận án đã mở rộng lý thuyết thủy động lực học tia phun của Dent (1971) và Hiroyasu (1980) cho dòng chất lỏng phi tiêu chuẩn có độ nhớt động học cao ($5,31\text{ cSt}$ ở PO20 so với $2,60\text{ cSt}$ ở DO) và sức căng bề mặt lớn ($22,8\text{ mN/m}$). Luận án đã mô hình hóa toán học thành công mối quan hệ giữa các chỉ số không thứ nguyên ($Re_l, We_l$) với sự gia tăng chiều dài phân rã sơ cấp $L_b$ và đường kính hạt Sauter $SMD$, đồng thời chứng minh bằng nhiệt động lực học rằng trị số Cetan cao của dầu cọ ($CN > 50$) là nhân tố chi phối rút ngắn thời gian cháy trễ $\tau_i$, làm đảo lộn quy luật áp suất cháy $p_z$ so với nhiên liệu hóa thạch truyền thống.
2. Cải tiến phương pháp luận so với các công trình nghiên cứu trước đây?
So với nghiên cứu của Herchel Thaddeus (2001), Kalam & Masjuki (2003) (chỉ thử nghiệm đơn thuần không hiệu chỉnh HTPNL) và Mai Đức Nghĩa (2014) (chỉ khảo sát đơn biến GPS trên động cơ cao tốc), luận án đã thực hiện một bước nhảy vọt về phương pháp luận:
- Tích hợp mô phỏng số học CFD 3D trường dòng 10 lỗ phun bằng Ansys Fluent với việc số hóa hình học chính xác bằng AutoCAD Mechanical 2017.
- Xác thực trường cháy trong buồng đốt động cơ công suất lớn (Hanshin 6LU32) bằng hệ thống nội soi buồng đốt Visio Scope tại 6 mốc GQTK khác nhau.
- Áp dụng phương pháp Quy hoạch thực nghiệm Bề mặt đáp ứng (RSM) để giải bài toán tối ưu hóa đa mục tiêu đồng thời giữa tính kinh tế ($g_e$) và phát thải môi trường ($NO_x$).
3. Phát hiện bất ngờ nhất và bằng chứng dữ liệu hỗ trợ?
Phát hiện bất ngờ nhất là mặc dù dầu thực vật có độ nhớt và nhiệt độ sôi cao hơn nhiều so với dầu diesel (dẫn đến quá trình bay hơi vật lý chậm hơn), nhưng thời gian cháy trễ $\tau_i$ của hỗn hợp dầu cọ - diesel (PO10, PO20, PO30) đo được trong thực nghiệm trên động cơ 6LU32 lại ngắn hơn so với dầu DO từ 0,5° đến 1,8° GQTK ở cả hai mức tải 400 kW và 600 kW. Bằng chứng thực nghiệm thu được từ cảm biến áp suất Kistler và hình ảnh Visio Scope chứng minh trị số Cetan cao ($CN = 50,13 - 52,11$) của phân tử dầu cọ kích hoạt phản ứng hóa học tự bốc cháy sớm hơn, làm dịch chuyển đỉnh áp suất cực đại $p_z$ về phía trước ĐCT.
4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?
Có, luận án cung cấp chi tiết toàn bộ điều kiện biên và quy trình thử nghiệm:
- Thông số kỹ thuật động cơ Hanshin 6LU32 (đường kính xy lanh, hành trình piston, 10 lỗ vòi phun).
- Tỷ lệ pha trộn phụ gia Nano Fuel Booster (1/800 theo thể tích).
- Phương pháp phân tích mẫu nhiên liệu theo QCVN: 2015/BKH&CN và chuẩn ASTM (D445, D1298, D613, D240).
- Cấu hình chia lưới CFD trong Gambit, mô hình rối $k-\epsilon$, các hằng số mô hình ($C_\mu = 0,09; C_\epsilon = 1,92$).
- Ma trận quy hoạch thực nghiệm RSM Box-Behnken với miền biến thiên cụ thể của 3 biến điều khiển ($X_1, X_2, X_3$).
5. Lộ trình nghiên cứu 10 năm được vạch ra như thế nào?
Lộ trình nghiên cứu phát triển tiếp theo bao gồm:
- Giai đoạn 1 (1 - 3 năm): Hoàn thiện module điều khiển nhiệt độ hâm dầu tự động thích ứng với độ nhớt trực tuyến trên tàu biển; thử nghiệm độ bền liên tục 2.000 giờ trên tuyến vận tải biển nội địa.
- Giai đoạn 2 (4 - 6 năm): Ứng dụng công nghệ trên động cơ diesel 2 kỳ cỡ lớn có hệ thống siêu nạp tuabin khí xả và hệ thống phun dầu điện tử Common Rail.
- Giai đoạn 3 (7 - 10 năm): Xây dựng hệ sinh thái nhiên liệu sinh học hàng hải khép kín (pha trộn Biofuel thế hệ 3 từ vi tảo và phụ gia nano xúc tác), tích hợp trí tuệ nhân tạo (AI/Machine Learning) điều khiển góc phun sớm thích ứng theo thời gian thực nhằm đạt mức phát thải ròng bằng 0 (Net Zero) theo định hướng của IMO đến năm 2050.
Kết luận
Luận án Tiến sĩ Kỹ thuật của NCS. Nguyễn Đức Hạnh đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra với các đóng góp cốt lõi:
- Xác lập hệ thống cơ sở lý thuyết toàn diện về động lực học tia phun, quá trình phân rã sơ cấp, hòa trộn và bốc cháy của hỗn hợp nhiên liệu dầu cọ - diesel (PO) trên động cơ diesel tàu thủy.
- Khám phá và giải thích bản chất vật lý của hiện tượng rút ngắn thời gian cháy trễ $\tau_i$ của nhiên liệu PO do ưu thế trị số Cetan cao, dẫn đến sự dịch chuyển đỉnh áp suất cháy cực đại $p_z$ về phía trước ĐCT.
- Ứng dụng thành công kỹ thuật mô phỏng CFD 3D (Ansys Fluent) mô tả chính xác trường áp suất, vector vận tốc và động năng rối qua chùm 10 lỗ phun và không gian buồng đốt động cơ Hanshin 6LU32.
- Kiểm chứng thực nghiệm chuẩn xác trên động cơ công suất 400 kW và 600 kW, kết hợp kỹ thuật chẩn đoán áp suất chỉ thị và nội soi quang học Visio Scope ghi nhận chân thực ngọn lửa buồng đốt.
- Đề xuất giải pháp hiệu chỉnh tối ưu hóa đa mục tiêu dựa trên lý thuyết Bề mặt đáp ứng (RSM): Xác lập bộ thông số kỹ thuật tối ưu (hỗn hợp PO20, lùi góc phun sớm 1,5° - 2,0° GQTK, tăng áp suất nâng kim phun lên 215 - 220 bar), giải quyết triệt để mâu thuẫn giữa việc cắt giảm phát thải ô nhiễm $NO_x$ theo chuẩn IMO MARPOL 73/78 Annex VI và duy trì tính kinh tế nhiên liệu của động cơ.
- Mở ra hướng ứng dụng thực tiễn to lớn, đóng góp trực tiếp vào Đề tài Khoa học Công nghệ cấp Quốc gia ĐT.11/NLSH, tạo tiền đề kỹ thuật vững chắc cho quá trình chuyển đổi năng lượng xanh của ngành hàng hải Việt Nam và khu vực.
