Tổng quan nghiên cứu

Sự gia tăng nhanh chóng của số lượng phương tiện cơ giới trên toàn cầu đang đặt ra những thách thức nghiêm trọng về cạn kiệt tài nguyên dầu mỏ và ô nhiễm môi trường. Theo các số liệu khảo sát năng lượng quốc tế, trữ lượng dầu mỏ thế giới ước tính còn khoảng 150 tỷ tấn, với mức tiêu thụ vượt mốc 3,6 tỷ tấn mỗi năm, nguy cơ cạn kiệt nguồn nhiên liệu hóa thạch có thể diễn ra trong vòng 40 đến 50 năm tới nếu không có các phát hiện mới. Động cơ diesel hiện chiếm hơn 50% tổng số lượng động cơ đốt trong trên thế giới nhờ hiệu suất nhiệt cao; tuy nhiên, khí xả của động cơ diesel lại chứa nhiều thành phần độc hại, đặc biệt là oxit nitơ và chất ô nhiễm dạng hạt gây nguy hại cho sức khỏe con người.

Để giải quyết tình trạng này, xu hướng đa dạng hóa năng lượng bằng cồn sinh học như butanol, ethanol và methanol phối trộn với dầu diesel khoáng đang được ứng dụng rộng rãi. Mặc dù các loại cồn chứa hàm lượng oxy cao giúp quá trình cháy sạch hơn, đặc tính bay hơi nhanh và thời gian cháy trễ kéo dài lại làm tăng tốc độ tỏa nhiệt tức thời, tác động trực tiếp đến phụ tải nhiệt của các chi tiết bao quanh buồng cháy. Bề mặt gương xi lanh bắt buộc phải duy trì dưới ngưỡng 180 độ C để ngăn chặn hiện tượng phá hủy màng dầu bôi trơn và mài mòn kết cấu.

Mục tiêu cốt lõi của đề tài là xây dựng mô hình mô phỏng chu trình nhiệt động và đánh giá định lượng ảnh hưởng của ba loại nhiên liệu phối trộn 10% cồn sinh học gồm B10 (pha 10% butanol), E10 (pha 10% ethanol) và M10 (pha 10% methanol) với 90% diesel khoáng (D100) đến trạng thái nhiệt của ống lót xi lanh trên động cơ V12 ở chế độ công suất định mức. Nghiên cứu được thực hiện tại Trường Đại học Kỹ thuật Công nghiệp - Đại học Thái Nguyên, mang ý nghĩa khoa học và thực tiễn sâu sắc khi cung cấp cơ sở dữ liệu nhiệt động học tin cậy, làm chủ công nghệ mô phỏng và hỗ trợ tối ưu hóa các thông số vận hành động cơ lưỡng nhiên liệu tại Việt Nam.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng nhiệt động lực học động cơ đốt trong và lý thuyết truyền nhiệt phức hợp bao gồm dẫn nhiệt, đối lưu và bức xạ nhiệt. Trong buồng cháy động cơ diesel, nhiệt độ môi chất công tác biến thiên từ khoảng 310 K đến 420 K ở cuối kỳ nạp và đạt đỉnh từ 1750 K đến 2800 K trong quá trình cháy, tạo ra dòng nhiệt truyền qua thành vách buồng cháy có thể lên tới 10 MW/m2. Truyền nhiệt bức xạ từ các hạt muội than trong ngọn lửa diesel chiếm từ 20% đến 35% tổng lượng nhiệt trao đổi, cao gấp khoảng 5 lần so với bức xạ từ pha khí.

Để giải quyết bài toán nhiệt động phức tạp, mô hình cháy đa vùng cho tia phun trực tiếp được áp dụng, phân chia tia phun thành 5 vùng theo phương hướng kính và tối đa 80 vùng dọc theo trục tia để tính toán chính xác sự thâm nhập của không khí, độ bay hơi của hạt nhiên liệu và động học tỏa nhiệt. Quá trình trao đổi nhiệt đối lưu giữa môi chất và thành buồng cháy được mô hình hóa theo phương trình thực nghiệm Hohenberg với các hằng số hiệu chỉnh đặc trưng (C1 bằng 130 và C2 bằng 1,4), sử dụng chiều dài đặc trưng là đường kính khối cầu tương đương với thể tích xi lanh tức thời. Các khái niệm cốt lõi được chuẩn hóa bao gồm: ống lót xi lanh, áp suất chỉ thị trung bình, suất tiêu hao nhiên liệu có ích, hệ số truyền nhiệt đối lưu và gradient nhiệt độ thành vách.

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

Đối tượng nghiên cứu là động cơ diesel cao tốc V12 4 kỳ, 12 xi lanh bố trí chữ V, không tăng áp, buồng cháy thống nhất với đường kính xi lanh 150 mm, hành trình pít tông 180 mm đối với dãy bên trái và 186,7 mm đối với dãy bên phải, chiều dài thanh truyền 320 mm, tỷ số nén 15 và góc phun sớm nhiên liệu từ 30 đến 33 độ trước điểm chết trên. Hệ thống cấp nhiên liệu sử dụng vòi phun 7 lỗ với đường kính mỗi lỗ 0,25 mm, lượng nhiên liệu cung cấp ở chế độ định mức đạt 138,90 mg trên một chu trình.

Cỡ mẫu tính toán bao trùm toàn bộ chu trình công tác 720 độ góc quay trục khuỷu với bước tích phân chia nhỏ từ 0,1 đến 0,5 độ trên phần mềm mô phỏng nhiệt động một chiều chuyên dụng của hãng Gamma Technologies. Mô hình hình học ống lót xi lanh được chia lưới phần tử hữu hạn với hơn 120.000 phần tử bậc cao trên phần mềm ANSYS. Phương pháp phân tích số một chiều kết hợp phần tử hữu hạn ba chiều được lựa chọn vì tính chính xác cao, giúp mô phỏng chi tiết các điều kiện biên nhiệt phức tạp và tiết kiệm tối đa chi phí so với phương pháp thử nghiệm phá hủy trên động cơ công suất lớn.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng chu trình công tác trên động cơ V12 với các loại nhiên liệu D100, B10, E10 và M10 đã chỉ ra những chuyển biến rõ rệt về đặc tính nhiệt và chỉ tiêu kinh tế kỹ thuật:

Thứ nhất, nhiệt độ môi chất công tác cực đại trong xi lanh khi sử dụng nhiên liệu hỗn hợp cồn - diesel có xu hướng tăng nhẹ so với diesel nguyên chất D100, đạt giá trị trong khoảng 2150 K đến 2350 K. Nguyên nhân do hàm lượng oxy có sẵn trong phân tử cồn (methanol chứa khoảng 50% oxy, ethanol khoảng 34,8% oxy và butanol khoảng 21,6% oxy) thúc đẩy quá trình oxy hóa triệt để hơn trong giai đoạn cháy khuếch tán.

Thứ hai, hệ số truyền nhiệt tức thời từ khí cháy tới thành vách buồng cháy tăng từ 3% đến 6,8% đối với nhiên liệu E10 và M10 so với D100 tại thời điểm lân cận điểm chết trên của kỳ cháy giãn nở. Nhiên liệu B10 thể hiện sự tương đồng nhất với D100 khi mức tăng hệ số truyền nhiệt chỉ dao động dưới 2,5%.

Thứ ba, trường nhiệt độ trên ống lót xi lanh phân bố không đồng đều theo chiều cao thân xi lanh. Nhiệt độ cao nhất tập trung tại vùng gờ phía trên tiếp giáp nắp máy, dao động trong khoảng từ 145 độ C đến 168 độ C đối với cả 4 loại nhiên liệu, vẫn duy trì dưới ngưỡng giới hạn nguy hiểm 180 độ C. Vùng thân dưới tiếp xúc trực tiếp với áo nước làm mát có nhiệt độ ổn định từ 80 độ C đến 95 độ C.

Thứ tư, suất tiêu hao nhiên liệu có ích khi sử dụng M10 và E10 tăng từ 2,2% đến 4,8% so với D100 do nhiệt trị của cồn thấp hơn diesel khoáng, trong khi công suất và áp suất chỉ thị trung bình của động cơ hầu như không bị suy giảm đáng kể.

Thảo luận kết quả

Diễn biến biến thiên nhiệt độ và hệ số truyền nhiệt được minh chứng rõ nét qua đồ thị đường cong nhiệt động theo góc quay trục khuỷu và các bảng tổng hợp chỉ tiêu công tác. Đồ thị cho thấy các đỉnh nhiệt độ của E10 và M10 xuất hiện sớm hơn và dốc hơn so với D100, phản ánh tốc độ tỏa nhiệt ban đầu cao do thời gian cháy trễ kéo dài và nhiệt ẩn bay hơi lớn của cồn làm tăng lượng hòa khí hòa trộn trước khi bốc cháy.

Dữ liệu trường nhiệt độ của ống lót xi lanh được trình bày trực quan qua các sơ đồ phân bố đẳng nhiệt dạng khối không gian. Gradient nhiệt độ lớn nhất xuất hiện theo phương hướng kính tại khu vực 1/3 phía trên hành trình pít tông, tạo nên ứng suất nhiệt cục bộ. Tuy nhiên, nhờ cấu trúc áo nước làm mát tuần hoàn hiệu quả, mức nhiệt tối đa 168 độ C của M10 vẫn đảm bảo độ bền nhiệt của vật liệu gang hợp kim và bảo toàn màng dầu bôi trơn. Kết quả này hoàn toàn tương thích với các công trình nghiên cứu của các nhà khoa học trong và ngoài nước về phụ tải nhiệt động cơ đốt trong, đồng thời khẳng định tính an toàn cơ học khi chuyển đổi sang nhiên liệu lưỡng nhiên liệu cồn sinh học tỷ lệ 10%.

Đề xuất và khuyến nghị

Nhằm tối ưu hóa hiệu quả vận hành và kiểm soát trạng thái nhiệt của ống lót xi lanh khi sử dụng nhiên liệu cồn sinh học, nghiên cứu đưa ra 4 nhóm giải pháp cụ thể:

Thứ nhất, ưu tiên lựa chọn nhiên liệu butanol sinh học B10 cho các dòng động cơ diesel công suất lớn. Do butanol có chuỗi carbon dài hơn, trị số cetane cao hơn và khả năng hòa trộn đồng nhất với diesel mà không cần chất trợ tan, B10 giúp động cơ duy trì tải nhiệt ổn định tương đương diesel truyền thống, hướng tới mục tiêu giữ vững tuổi thọ động cơ trên 15.000 giờ làm việc.

Thứ tư, hiệu chỉnh góc phun sớm nhiên liệu lùi lại từ 1 đến 3 độ góc quay trục khuỷu khi chuyển sang sử dụng E10 và M10. Biện pháp này giúp hạ thấp áp suất cực đại và giảm nhiệt độ đỉnh của khí cháy khoảng 50 K đến 80 K, từ đó giảm ứng suất nhiệt lên bề mặt ống lót xi lanh và giảm nồng độ phát thải oxit nitơ trong khí xả. Kế hoạch thực hiện nên được áp dụng ngay trong quá trình bảo dưỡng định kỳ tại các xưởng kỹ thuật.

Thứ ba, nâng cấp hệ thống làm mát động cơ thông qua việc duy trì nhiệt độ nước làm mát đầu vào ổn định trong khoảng 75 độ C đến 85 độ C, kết hợp tăng lưu lượng tuần hoàn qua khu vực đỉnh xi lanh khoảng 5% đến 8%. Giải pháp này do các kỹ sư vận hành thực hiện nhằm kiểm soát gradient nhiệt độ ngang thành vách dưới mức 60 độ C, ngăn ngừa nguy cơ rạn nứt vì mỏi nhiệt.

Thứ tư, các cơ quan chức năng thuộc ngành năng lượng và giao thông vận tải cần phối hợp với các doanh nghiệp sản xuất ban hành tiêu chuẩn kỹ thuật quốc gia về phụ gia chống phân tách pha và phụ gia chống ăn mòn cho hỗn hợp nhiên liệu cồn - diesel trong lộ trình từ nay đến năm 2030, tạo hành lang pháp lý vững chắc cho việc thương mại hóa nhiên liệu sinh học.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị thực tiễn cho nhiều nhóm chuyên môn:

Nhóm 1: Các kỹ sư nghiên cứu và phát triển tại các nhà máy chế tạo động cơ, cơ sở thiết kế và cải tạo phương tiện cơ giới. Luận văn cung cấp quy trình hoàn chỉnh về việc thiết lập mô hình tính toán nhiệt động lực học trên phần mềm chuyên dụng và mô phỏng ứng suất nhiệt phần tử hữu hạn.

Nhóm 2: Các cán bộ kỹ thuật quản lý và vận hành đội xe tải nặng, máy công trình, trạm phát điện công nghiệp và phương tiện thủy. Nghiên cứu giúp nắm vững các thông số kỹ thuật, cách thức điều chỉnh hệ thống nhiên liệu và làm mát khi sử dụng nhiên liệu sinh học pha cồn an toàn.

Nhóm 3: Giảng viên, nghiên cứu sinh và sinh viên chuyên ngành Kỹ thuật Cơ khí Động lực, Kỹ thuật Ô tô, Máy năng lượng tại các trường đại học kỹ thuật. Luận văn là tài liệu tham khảo chuyên sâu về mô hình hóa chu trình cháy đa vùng và lý thuyết truyền nhiệt trong không gian buồng cháy.

Nhóm 4: Các nhà hoạch định chính sách môi trường và an ninh năng lượng tại các cơ quan quản lý nhà nước. Dữ liệu từ luận văn cung cấp luận cứ khoa học thực tế để xây dựng các kịch bản giảm phát thải khí nhà kính và đa dạng hóa nguồn nhiên liệu quốc gia.

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

Tại sao nhiệt độ bề mặt gương xi lanh bắt buộc phải được khống chế dưới 180 độ C? Nhiệt độ bề mặt gương xi lanh vượt quá 180 độ C sẽ khiến màng dầu nhờn bôi trơn bị phân hủy nhiệt, giảm độ nhớt và mất khả năng bôi trơn thủy động. Điều này dẫn đến ma sát khô giữa séc-măng và thành xi lanh, gây ra hiện tượng bó kẹt pít tông, mài mòn nghiêm trọng và phá hủy kết cấu động cơ.

Nhiên liệu butanol B10 có điểm gì ưu việt hơn so với ethanol E10 và methanol M10? Butanol có mạch 4 carbon nên sở hữu trị số cetane cao hơn, nhiệt trị lớn hơn và tính hòa tan trong dầu diesel vượt trội mà không bị tách pha khi gặp ẩm. Ngoài ra, nhiệt ẩn bay hơi của butanol thấp hơn ethanol và methanol, giúp chu trình cháy diễn ra êm dịu hơn và duy trì phụ tải nhiệt ổn định tương đương diesel khoáng.

Mô hình truyền nhiệt Hohenberg có ưu điểm gì so với mô hình Woschni truyền thống? Mô hình Woschni thường dự báo quá cao dòng nhiệt cực đại trong quá trình cháy và dự báo quá thấp trong quá trình nén và thải đối với động cơ diesel cao tốc. Mô hình Hohenberg hiệu chỉnh công thức dựa trên chiều dài đặc trưng của thể tích hình cầu tức thời, mang lại kết quả sát với thực nghiệm đo đạc dòng nhiệt thực tế.

Hiệu suất và công suất của động cơ V12 thay đổi như thế nào khi sử dụng nhiên liệu hỗn hợp cồn 10%? Công suất định mức của động cơ hầu như được giữ nguyên do hệ thống vòi phun cung cấp đủ khối lượng chu trình. Tuy nhiên, suất tiêu hao nhiên liệu có ích tăng nhẹ từ 2% đến 4,8% khi dùng E10 và M10 để bù đắp phần mật độ năng lượng thấp hơn của cồn, đổi lại lượng khói đen và bồ hóng giảm đáng kể.

Việc kết hợp giữa phần mềm GT-Power và ANSYS mang lại lợi ích gì cho bài toán nhiệt động cơ? GT-Power đóng vai trò giải chu trình nhiệt động một chiều để xuất ra quy luật áp suất, nhiệt độ môi chất và hệ số truyền nhiệt theo thời gian. Các thông số này được sử dụng làm điều kiện biên chính xác cho ANSYS giải bài toán phần tử hữu hạn ba chiều, giúp xác định chuẩn xác trường nhiệt độ cục bộ của chi tiết máy.

Kết luận

  • Luận văn đã xây dựng thành công mô hình chu trình nhiệt động lực học trên phần mềm GT-Power kết hợp mô hình cháy DI-Jet và mô hình truyền nhiệt Hohenberg cho động cơ diesel V12.
  • Kết quả mô phỏng đã lượng hóa chi tiết sự biến thiên của nhiệt độ môi chất và hệ số truyền nhiệt đối lưu khi sử dụng ba loại nhiên liệu sinh học B10, E10 và M10 so với diesel nguyên chất D100.
  • Nghiên cứu đã tính toán trường nhiệt độ ba chiều của ống lót xi lanh trên phần mềm ANSYS, chứng minh nhiệt độ đỉnh của thành vách đạt mức 145 đến 168 độ C, hoàn toàn nằm trong giới hạn an toàn nhiệt của vật liệu.
  • Nhiên liệu butanol B10 được xác định là giải pháp thay thế tối ưu nhất về độ ổn định nhiệt động lực học và tính tương thích với kết cấu động cơ hiện hành.
  • Nghiên cứu mở ra hướng phát triển thử nghiệm thực nghiệm trên băng thử công suất và mở rộng khảo sát ứng suất nhiệt cho cụm nắp xi lanh và pít tông trong các giai đoạn tiếp theo.

Các kết quả khoa học đạt được từ công trình là tiền đề quan trọng thúc đẩy ứng dụng nhiên liệu sinh học cồn - diesel trên các phương tiện vận tải cơ giới tại Việt Nam, góp phần giảm thiểu ô nhiễm môi trường và phát triển năng lượng bền vững.