Tổng quan nghiên cứu

Theo báo cáo của Hiệp hội Kho lạnh Quốc tế, tổng dung tích kho lạnh toàn cầu đạt khoảng 458 triệu m³, trong đó hệ sinh thái lưu trữ lạnh phục vụ chuỗi cung ứng thực phẩm chiếm trên 310 triệu m³. Tại Việt Nam, tiềm năng xuất khẩu nông thủy sản đạt từ 10 đến 15 tỷ USD, song kim ngạch thực tế giai đoạn 2009–2010 mới đạt khoảng 4,3 tỷ USD, tương đương 15% dung lượng thị trường. Nguyên nhân cốt lõi bắt nguồn từ việc công nghệ lạnh đông chưa đồng bộ, thời gian vận chuyển sơ bộ kéo dài từ 6 đến 48 giờ và thiếu các tiêu chuẩn tính toán nhiệt động học chính xác cho từng loại thực phẩm bản địa. Tình trạng này khiến tỷ lệ hao hụt khối lượng và suất tiêu hao năng lượng trong nước cao hơn đáng kể so với các quốc gia trong khu vực.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Nhiệt – Lạnh của học viên Võ Đình Hiệp, dưới sự hướng dẫn khoa học của TS. Nguyễn Việt Dũng tại Trường Đại học Bách khoa Hà Nội, tập trung giải quyết bài toán cấp thiết: Nghiên cứu ảnh hưởng các thông số chính tới quá trình cấp đông thịt bò trong thiết bị đông gió bằng mô phỏng CFD. Mục tiêu trọng tâm của đề tài là xây dựng mô hình toán học giải đồng thời trường dòng chảy và quá trình truyền nhiệt chuyển pha, phát triển bộ công cụ mô phỏng số ANSYS Fluent kết hợp hàm người dùng định nghĩa, và thực nghiệm đối chứng trên hệ thống lạnh hai cấp nén.

Nghiên cứu được triển khai trong phạm vi thời gian từ tháng 10 năm 2010 đến tháng 8 năm 2012 tại Viện Khoa học và Công nghệ Nhiệt - Lạnh. Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp bộ dữ liệu chuẩn xác giúp rút ngắn 15% đến 25% thời gian cấp đông, đồng thời giảm 10% đến 18% chi phí điện năng tiêu thụ cho các hệ thống đông gió công nghiệp.

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 phương trình vi phân dẫn nhiệt phi tuyến tổng quát của Stefan (1889) kết hợp định luật dẫn nhiệt Fourier trong môi trường biến đổi pha. Quá trình lưu động của không khí lạnh được mô tả bằng hệ phương trình Navier-Stokes khép kín qua mô hình dòng rối hai phương trình k-epsilon ($k$-$\varepsilon$). Bên cạnh đó, luận văn kế thừa và phân tích giới hạn của các mô hình giải tích kinh điển như phương trình Plank, phương pháp ba giai đoạn của Phạm Quang Tuấn, công thức Nagaoka và mô hình thực nghiệm của Cleland - Earle.

Mô hình nghiên cứu tích hợp bốn khái niệm then chốt:

  • Nhiệt dung riêng hiệu dụng ($c_e$): Đại lượng đặc trưng cho sự biến đổi nhiệt lượng bao gồm cả nhiệt hiện và nhiệt ẩn đông đặc ($L_0 = 333,6\text{ kJ/kg}$) phân bố liên tục theo nhiệt độ.
  • Tỷ lệ đóng băng ($\omega$): Hàm số phi tuyến phụ thuộc vào nhiệt độ đóng băng ban đầu ($t_f \approx -1,5^\circ\text{C}$) và hàm lượng nước liên kết.
  • Mật độ dòng nhiệt ($q$): Đại lượng phản ánh cường độ tỏa nhiệt từ bề mặt thực phẩm ra môi trường lưu chất ($W/m^2$).
  • Tiêu chuẩn đồng dạng Biot ($Bi$) và Stefan ($Ste$): Các chỉ số không thứ nguyên xác định tỷ lệ nhiệt trở dẫn nhiệt so với nhiệt trở tỏa nhiệt đối lưu và tỷ số giữa nhiệt hiện với nhiệt ẩn chuyển pha.

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

Nghiên cứu kết hợp chặt chẽ giữa phương pháp mô phỏng số động lực học dòng chảy (CFD) và phương pháp thực nghiệm vật lý:

  • Cỡ mẫu và phương pháp chọn mẫu: Thực nghiệm sử dụng 12 mẫu thịt bò nạc tươi đồng nhất, được cắt gọt chuẩn hóa dạng hình hộp chữ nhật kích thước $100\text{ mm} \times 100\text{ mm} \times 50\text{ mm}$. Phương pháp chọn mẫu ngẫu nhiên có kiểm soát được áp dụng nhằm đảm bảo độ ẩm ban đầu đạt khoảng 74,5% và hàm lượng chất béo dưới 3,5%.
  • Nguồn dữ liệu và phân tích: Quá trình mô phỏng sử dụng phần mềm Gambit để tạo lưới phân tố phi cấu trúc với hơn 180.000 phần tử, kết hợp thuật toán SIMPLE trong ANSYS Fluent. Các đặc tính nhiệt vật lý phi tuyến phụ thuộc nhiệt độ (khối lượng riêng $\rho$, hệ số dẫn nhiệt $\lambda$, nhiệt dung riêng $c$) được lập trình bằng ngôn ngữ C thông qua hàm UDF (User Defined Functions). Dữ liệu thực nghiệm được thu thập tự động thông qua bộ ghi dữ liệu HIOKI 9768 Smart Suite Utility Pro, đầu đo nhiệt độ cặp nhiệt ngâm sâu tại tâm và bề mặt mẫu thịt, cùng thiết bị đo vận tốc gió đa điểm TSI Model 8346.
  • Lý do lựa chọn phương pháp: Phương pháp số kết hợp CFD là giải pháp tối ưu duy nhất cho phép giải quyết bài toán truyền nhiệt phi tuyến đa chiều có chuyển pha, khắc phục triệt để các sai số từ 20% đến 40% của các công thức giải tích kinh nghiệm do giả định hệ số dẫn nhiệt không đổi.
  • Thời gian nghiên cứu: Toàn bộ quy trình mô phỏng và thực nghiệm được tiến hành liên tục trong 22 tháng (10/2010 – 08/2012).

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

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

  1. Quy luật tác động của nhiệt độ buồng cấp đông: Khi hạ nhiệt độ môi trường buồng đông từ $-25^\circ\text{C}$ xuống $-40^\circ\text{C}$ ở vận tốc gió cố định 2,5 m/s, thời gian hạ nhiệt độ tâm thịt bò từ $+15^\circ\text{C}$ xuống $-18^\circ\text{C}$ giảm mạnh từ 215 phút xuống còn 118 phút, tương đương mức rút ngắn 45,1%.
  2. Ngưỡng hiệu quả của vận tốc gió: Tăng tốc độ gió từ 1,5 m/s lên 3,5 m/s giúp rút ngắn thời gian làm lạnh đông 28,4%. Tuy nhiên, khi tiếp tục nâng vận tốc từ 3,5 m/s lên 5,0 m/s, thời gian cấp đông chỉ giảm thêm 7,2%, trong khi công suất tiêu thụ của quạt gió tăng xấp xỉ 2,9 lần do áp lực động học tỷ lệ với bình phương vận tốc.
  3. Đặc tính biến thiên mật độ dòng nhiệt: Mật độ dòng nhiệt ($q$) đạt giá trị cực đại khoảng $1.820\text{ W/m}^2$ trong giai đoạn đầu của quá trình chuyển pha (vùng nhiệt độ từ $-1,5^\circ\text{C}$ đến $-5^\circ\text{C}$), sau đó giảm dần đều về mức $310\text{ W/m}^2$ khi toàn bộ thể tích mẫu thịt đã đóng băng hoàn toàn.
  4. Độ chính xác và độ tin cậy của mô hình: Sự so sánh giữa đường cong nhiệt độ mô phỏng CFD và dữ liệu đo đạc thực nghiệm cho thấy độ lệch cực đại tại tâm khối thịt không vượt quá 6,2%, và sai số trung bình toàn chu kỳ duy trì ở mức 3,8%.

Thảo luận kết quả

Diễn biến nhiệt độ trong mẫu thịt thể hiện rõ 3 giai đoạn nhiệt động học: làm lạnh sơ bộ giải phóng nhiệt hiện, giai đoạn dừng nhiệt giải phóng nhiệt ẩn kết tinh nước đá (chiếm hơn 62% tổng thời gian làm lạnh), và giai đoạn quá lạnh thể rắn hạ nhiệt độ xuống $-18^\circ\text{C}$. Trong thực tế, các dữ liệu này được thể hiện trực quan thông qua đồ thị đường biểu diễn quan hệ giữa nhiệt độ tâm – bề mặt theo thời gian và biểu đồ phân bố đường dòng vận tốc không khí bao quanh vật thể.

So với mô hình giải tích của Plank vốn xem bề mặt đóng băng phẳng chuyển động tịnh tiến đều, mô hình CFD phản ánh chính xác hiệu ứng biên dòng chảy và sự không đồng đều của hệ số trao đổi nhiệt đối lưu ($\alpha$) trên các mặt đón gió và khuất gió (dao động từ 25 đến $68\text{ W/m}^2\text{K}$). Kết quả này chứng minh rằng việc duy trì dải nhiệt độ buồng đông $-35^\circ\text{C}$ cùng vận tốc gió 3,0 m/s mang lại điểm cân bằng tối ưu nhất giữa chất lượng sản phẩm và chi phí năng lượng.

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

  1. Chuẩn hóa thông số vận hành thiết bị đông gió: Thiết lập nhiệt độ buồng cấp đông ở mức $-35^\circ\text{C}$ và dải tốc độ gió bề mặt từ 2,5 m/s đến 3,2 m/s. Giải pháp này giúp cắt giảm 14% điện năng tiêu hao so với chế độ vận hành cực trị $-40^\circ\text{C}$, áp dụng ngay trong quý 1 cho các kỹ sư vận hành tại nhà máy chế biến thực phẩm.
  2. Tái cấu trúc hình học và khí động học buồng đông: Điều chỉnh khoảng cách giữa các khay chứa thịt tối thiểu 60 mm và lắp đặt tấm hướng dòng không khí để triệt tiêu các vùng quẩn khí, nâng độ đồng đều vận tốc gió lên trên 92%. Mục tiêu hoàn thành cải tạo trong 6 tháng do các đơn vị tư vấn thiết kế cơ điện lạnh chủ trì.
  3. Ứng dụng biến tần điều khiển lưu lượng gió động: Tích hợp hệ thống biến tần cho quạt dàn lạnh nhằm tự động giảm 30% lưu lượng gió sau khi tâm sản phẩm vượt qua điểm kết tinh nhiệt ẩn (dưới $-6^\circ\text{C}$), tiết kiệm thêm 18% điện năng quạt trong lộ trình triển khai 9 tháng bởi phòng kỹ thuật bảo trì.
  4. Ban hành quy chuẩn công nghệ cấp đông thịt đỏ: Xây dựng bộ hướng dẫn kỹ thuật cấp đông tiêu chuẩn áp dụng các hàm nhiệt vật lý phi tuyến nhằm kiểm soát độ hao hụt khối lượng dưới 1,2%, do Hiệp hội Chế biến Thực phẩm phối hợp cùng các viện nghiên cứu thực hiện trong vòng 12 tháng.

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

  1. Kỹ sư thiết kế và vận hành hệ thống lạnh công nghiệp: Nắm vững phương pháp tính toán phụ tải nhiệt biến đổi theo thời gian thực để lựa chọn công suất máy nén và dàn bay hơi chính xác, tránh lãng phí 15% đến 20% chi phí đầu tư ban đầu.
  2. Cán bộ kỹ thuật và quản lý chất lượng tại doanh nghiệp thực phẩm: Ứng dụng chế độ cấp đông tối ưu để bảo vệ cấu trúc tế bào thịt bò, hạn chế sự phá hủy của tinh thể đá lớn, duy trì tỷ lệ giữ nước trên 98% sau rã đông.
  3. Giảng viên, nghiên cứu sinh và học viên cao học ngành Kỹ thuật Nhiệt: Sử dụng luận văn làm tài liệu tham khảo chuyên sâu về kỹ thuật lập trình UDF trong ANSYS Fluent và phương pháp mô hình hóa nhiệt động học hệ đa thành phần có chuyển pha.
  4. Chuyên viên tư vấn tiết kiệm năng lượng và kiểm toán năng lượng: Khai thác các chỉ số định lượng về dòng nhiệt và tổn thất khí động học để xây dựng giải pháp nâng cao hiệu suất năng lượng cho các kho lạnh quy mô trên 10.000 tấn.

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

Mô phỏng CFD mang lại ưu thế gì vượt trội so với các công thức tính toán truyền thống như Plank?

Phương pháp CFD giải quyết trực tiếp hệ phương trình vi phân truyền nhiệt kết hợp dòng chảy rối 3 chiều với các thông số nhiệt vật lý phụ thuộc nhiệt độ. Ngược lại, công thức Plank coi hệ số dẫn nhiệt là hằng số và bỏ qua nhiệt hiện, dẫn đến sai số thực tế có thể lên tới 25% đến 40%.

Vì sao không nên tăng vận tốc gió vượt quá ngưỡng 3,5 m/s trong buồng cấp đông gió?

Khi vận tốc gió tăng từ 3,5 m/s lên 5,0 m/s, thời gian cấp đông chỉ rút ngắn thêm khoảng 7,2% do nhiệt trở dẫn nhiệt bên trong miếng thịt trở thành yếu tố chi phối chính. Trong khi đó, công suất tiêu thụ của động cơ quạt tăng theo hàm bậc ba, gây tiêu hao năng lượng không hiệu quả.

Điểm đóng băng ban đầu của thịt bò là bao nhiêu và có ý nghĩa gì trong tính toán?

Điểm đóng băng ban đầu của thịt bò nằm trong khoảng $-1,5^\circ\text{C}$ đến $-1,8^\circ\text{C}$. Đây là mốc nhiệt độ bắt đầu xuất hiện các mầm tinh thể đá và giải phóng lượng nhiệt ẩn kết tinh $333,6\text{ kJ/kg}$, là thông số nền tảng để thiết lập hàm nhiệt dung riêng hiệu dụng trong mô hình số.

Sai số giữa kết quả mô phỏng số CFD và đo đạc thực nghiệm trong nghiên cứu là bao nhiêu?

Độ sai lệch tương đối giữa mô hình ANSYS Fluent nhúng UDF và kết quả đo thực tế bằng hệ thống HIOKI 9768 dao động trong khoảng 3,8% đến 6,2%. Sai số này hoàn toàn nằm trong giới hạn cho phép của các tiêu chuẩn nghiên cứu nhiệt lạnh quốc tế.

Kết quả nghiên cứu có thể mở rộng áp dụng cho các loại thực phẩm khác không?

Mô hình toán và cấu trúc mã nguồn UDF hoàn toàn có thể tái sử dụng cho thịt lợn, thịt gia cầm hoặc thủy hải sản. Người dùng chỉ cần cập nhật lại các tỷ lệ thành phần hóa học ban đầu gồm hàm lượng nước ($x_w$), protein ($x_p$), chất béo ($x_f$) và nhiệt độ đóng băng ban đầu tương ứng.

Kết luận

  • Xây dựng thành công mô hình toán học toàn diện và các module hàm UDF bằng ngôn ngữ C mô tả chính xác sự phụ thuộc phi tuyến của tính chất nhiệt vật lý thịt bò vào nhiệt độ.
  • Xác định định lượng ảnh hưởng của nhiệt độ buồng cấp đông (từ $-25^\circ\text{C}$ đến $-40^\circ\text{C}$) và vận tốc gió (từ 1,5 m/s đến 5,0 m/s) đến động học quá trình đóng băng và phân bố mật độ dòng nhiệt.
  • Đạt độ chính xác thực nghiệm cao với sai số mô phỏng chỉ từ 3,8% đến 6,2% trên hệ thống thiết bị đo chuẩn HIOKI và TSI.
  • Đề xuất dải thông số công nghệ tối ưu tại nhiệt độ buồng $-35^\circ\text{C}$ và tốc độ gió 3,0 m/s, giúp tiết kiệm từ 12% đến 18% năng lượng vận hành.
  • Định hình khung phương pháp luận hiện đại mở đường cho việc tối ưu hóa thiết kế buồng cấp đông công nghiệp cho nhiều chủng loại nông sản xuất khẩu khác nhau.

Kế hoạch tiếp theo bao gồm việc mở rộng mô hình cho các dạng bao gói hình học phức tạp trong 6 tháng tới và tích hợp thuật toán học máy dự báo phụ tải lạnh theo thời gian thực trong vòng 12 tháng. Các nhà máy chế biến và đơn vị thiết kế cơ điện lạnh nên áp dụng ngay bộ thông số tối ưu này vào dây chuyền sản xuất để nâng cao năng lực cạnh tranh và chất lượng sản phẩm xuất khẩu.