Tổng quan luận án

Vấn đề bảo quản sau thu hoạch đối với quả vải thiều tươi (Litchi chinensis Sonn) tại Việt Nam mang tính cấp thiết xuất phát từ thực trạng sản xuất và đặc tính sinh lý của nông sản. Theo số liệu thống kê năm 2016, tổng sản lượng vải cả nước đạt 312.556,3 tấn, tập trung chủ yếu tại tỉnh Bắc Giang (khoảng 142.000 tấn với diện tích 30.000 ha riêng tại huyện Lục Ngạn), Hải Dương (khoảng 25.000 tấn) và Quảng Ninh (khoảng 8.000 tấn). Quả vải có đặc tính hô hấp thường biến (non-climacteric), hàm lượng nước cao (80 ÷ 84%), vỏ mỏng nhưng thô ráp, dễ mất nước nhanh (tỷ lệ hao hụt khối lượng ở nhiệt độ môi trường trong 24 giờ sau thu hoạch đạt 3 ÷ 5%) dẫn đến hiện tượng nâu hóa vỏ và suy giảm chất lượng trong vòng 48 giờ. Mùa vụ thu hoạch diễn ra ngắn (khoảng 30 ÷ 45 ngày từ giữa tháng 6 đến cuối tháng 7), gây áp lực tiêu thụ và vận chuyển. Tỷ lệ tổn thất sau thu hoạch của quả vải tại Việt Nam lên tới 25 ÷ 30% do thiếu quy trình công nghệ lạnh đồng bộ, chủ yếu xuất khẩu thô hoặc tiêu thụ nội địa. Để thực hiện Quyết định số 1895/QĐ-TTg ngày 17/12/2012 của Thủ tướng Chính phủ và Nghị quyết số 48/NQ-CP ngày 23/09/2009 nhằm giảm tổn thất rau quả sau thu hoạch xuống dưới 10%, việc nghiên cứu các thông số công nghệ lạnh cho quả vải là yêu cầu bắt buộc.

Khoảng trống nghiên cứu được luận án chỉ ra là: trước thời điểm nghiên cứu, tại Việt Nam chưa có công trình nghiên cứu hệ thống và đồng bộ về đặc tính nhiệt vật lý hiệu dụng, quá trình truyền nhiệt dẫn nhiệt không ổn định, cường độ hô hấp phụ thuộc nhiệt độ và cơ chế tỏa ẩm hao hụt tự nhiên của quả vải. Các nghiên cứu quốc tế trước đó thường áp đặt thông số nhiệt vật lý của loại quả khác (như quả cherry) hoặc chỉ tiến hành thí nghiệm đơn lẻ, thiếu mô hình toán học tổng quát gắn với chuỗi công nghệ lạnh hoàn chỉnh.

Mục tiêu chung của luận án là xây dựng cơ sở khoa học cho việc phát triển công nghệ chế biến và bảo quản lạnh quả vải quy mô tập trung. Bốn mục tiêu cụ thể được đánh số bao gồm:

  1. Xây dựng mô hình tính toán tính chất nhiệt vật lý hiệu dụng của quả vải.
  2. Xây dựng mô hình mô phỏng quá trình làm lạnh quả vải.
  3. Xây dựng mô hình hô hấp của quả vải phụ thuộc nhiệt độ.
  4. Xây dựng mô hình xác định độ hao hụt tự nhiên của quả vải trong quá trình bảo quản quy mô tập trung.

Đối tượng nghiên cứu là quả vải thiều tươi thu hoạch tại huyện Lục Ngạn, tỉnh Bắc Giang, đạt độ chín 3 (khối lượng trung bình 15 ÷ 25 g, đường kính 25 ÷ 35 mm, chiều cao 26 ÷ 40 mm, hàm lượng nước trong cùi chiếm 80 ÷ 84%).

Phạm vi nghiên cứu bao gồm: các tính chất nhiệt vật lý (nhiệt dung riêng đẳng áp, hệ số dẫn nhiệt, khối lượng riêng, enthalpy), quá trình truyền nhiệt và tỏa ẩm, cường độ hô hấp và hao hụt khối lượng xảy ra trong các khâu làm lạnh sơ bộ (bằng nước đá và không khí đối lưu), vận chuyển và bảo quản lạnh. Dữ liệu thực nghiệm được thu thập liên tiếp trong 3 mùa vụ (2014, 2015, 2016).

Tổng quan tài liệu và vị trí của luận án

Luận án đã phân loại và điểm luận có hệ thống các công trình nghiên cứu trong và ngoài nước về công nghệ bảo quản và truyền nhiệt ẩm nông sản:

Về phương pháp bảo quản hóa học: Luận án phân tích phương pháp xử lý khí $SO_2$ và $H_2SO_3$ nhằm ức chế enzym polyphenoloxidase chống nâu hóa vỏ (được nêu qua các tài liệu của ngành), tuy nhiên phương pháp này để lại dư lượng độc hại nên bị cấm; phương pháp sử dụng màng bao phủ chitosan kết hợp axit hữu cơ citric/tartaric của Jacques Joas (2005); phương pháp xử lý bằng hợp chất 1-methylcyclopropene (1-MCP) nồng độ 300 ÷ 1.000 ppm của Blankenship và Dole (2003), Mahajan (2009).

Về phương pháp bảo quản vật lý:

  • Công nghệ khí kiểm soát (CA): Nghiên cứu của Datta và cộng sự (1963), Kader (1993), Cameren (1995) chỉ ra môi trường nhiệt độ 7°C, nồng độ $O_2 \le 5%$ và $CO_2 = 3 \div 5%$ làm giảm tốc độ tổn thất vitamin C và chất rắn hòa tan, nhưng chi phí thiết bị cao và quả dễ hỏng nhanh khi đưa ra môi trường thường.
  • Công nghệ bao gói khí điều biến (MAP): Các công trình của Chen và cộng sự (1982), Kader (1994) chứng minh nồng độ $O_2 \approx 5%$ và $CO_2 = 3 \div 5%$ kéo dài thời gian bảo quản quả vải 7 ÷ 8 ngày ở nhiệt độ thường và 25 ÷ 30 ngày ở nhiệt độ thấp.
  • Bảo quản lân cận điểm băng: Nghiên cứu của Shen Jiang và cộng sự (2011) hạ nhiệt độ theo các bước (5°C $\rightarrow$ 3°C $\rightarrow$ 0°C $\rightarrow$ -0,5°C $\rightarrow$ -1,2°C; độ ẩm 87%), kế thừa nguyên lý của I. Trumac từ cuối thập niên 1970.
  • Bảo quản trong điện từ trường yếu: Khống chế chuyển động ion $H^+$ và phân tử nước trong quá trình hô hấp.
  • Chế độ bảo quản lạnh tiêu chuẩn: Khan và cộng sự (2012), Paull và Chen (2014), Kader (2000b), Cao Văn Hùng (2006) thống nhất dải nhiệt độ bảo quản tối ưu cho quả vải là 2 ÷ 5°C, độ ẩm tương đối 90 ÷ 95%.

Về quá trình làm lạnh sơ bộ và mô phỏng truyền nhiệt:

  • Brosnan (2001) khẳng định làm lạnh sơ bộ trong vòng 4 giờ sau thu hoạch là khâu quyết định.
  • Ketsa và Leelawatana (1992) thực nghiệm hạ nhiệt độ vải từ 26°C xuống 5°C trong 18 phút bằng nước lạnh 3 ÷ 5°C.
  • Olesen và cộng sự (2004) làm lạnh tâm quả từ 28°C xuống 5°C mất 10 phút bằng nước 2°C, nhưng mất 3 giờ nếu dùng không khí 2°C; mô hình mô phỏng của nhóm tác giả này tồn tại sai số do sử dụng thông số nhiệt vật lý của quả cherry thay cho quả vải.
  • Yu Shen Liang và cộng sự (2013) làm lạnh vải bằng nước 5 ± 2°C trong 20 phút để nhiệt độ tâm đạt 7,1°C.
  • Holcroft và cộng sự (1996) đối chứng làm lạnh bằng nước (18 phút) và không khí đối lưu tốc độ 2 m/s (71 phút).
  • Bryant (2012) định lượng tỷ lệ mất nước khi làm lạnh bằng không khí đối lưu nhẹ (0,38%), đối lưu cưỡng bức (0,20%) và nước lạnh (0%).

Về tính chất nhiệt vật lý và hao hụt khối lượng:

  • Mô hình nhiệt dung riêng ($C_p$): Siebel (1892), Charm (1978), Heldman và Singh (1981), Choi và Okos (1983, 1986), Chen (1985).
  • Mô hình hệ số dẫn nhiệt ($\lambda$): Novikov (1977), Sweat (1974, 1986), Murakami và Okos (1989), Dincer (1995, 1997), Carson (2006), Wang Handong (2006).
  • Xác định hệ số trao đổi nhiệt ($\alpha$): Mariangela Amendola và cộng sự (2009).
  • Mô hình hao hụt tự nhiên và truyền ẩm: Định luật Fick, định luật Dalton, Arce và Sweat (1980), Becker và Fricke (1996), Markarian và cộng sự (2006), phương pháp tia quá trình $\varepsilon$ của Jadan.

Vị trí của luận án: Luận án giải quyết khoảng trống về việc thiếu hụt mô hình tính toán nhiệt vật lý chuyên biệt cho quả vải nhiều lớp (vỏ, cùi, hạt); phát triển phương pháp số giải bài toán dẫn nhiệt không ổn định kết hợp thuật toán tối ưu gián tiếp xác định hệ số trao đổi nhiệt đối lưu; và thiết lập mô hình hao hụt tự nhiên theo hệ nhiệt động tựa cân bằng đẳng áp.

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

Cơ sở lý thuyết và khung phân tích

  • Lý thuyết dẫn nhiệt không ổn định: Xem xét quả vải là một khối cầu đẳng hướng tương đương, phương trình vi phân dẫn nhiệt 1 chiều trong hệ tọa độ cầu có dạng:

$$\rho(t) C(t) \frac{\partial t(r, \tau)}{\partial \tau} = \frac{1}{r^2} \frac{\partial}{\partial r} \left( \lambda(t) r^2 \frac{\partial t(r, \tau)}{\partial r} \right) + q(r, \tau)$$

Trong đó $q(r, \tau)$ là nguồn nhiệt nội sinh do hô hấp; điều kiện biên loại 3 (Robin) tại bề mặt bán kính $R$:

$$-\lambda \left. \frac{\partial t}{\partial r} \right|{r=R} = \alpha (t_b - t\infty)$$

Điều kiện đối xứng tại tâm quả: $\left. \frac{\partial t}{\partial r} \right|_{r=0} = 0$.

  • Mô hình tính chất nhiệt vật lý hiệu dụng: Coi quả vải là hệ nhiệt động cân bằng đa thành phần không chuyển pha (gồm vỏ, cùi, hạt với các thành phần nước, carbohydrate, protein, chất béo, tro).
    • Nhiệt dung riêng hiệu dụng đẳng áp kế thừa mô hình Choi và Okos (1986): $C_p = \sum C_{pi} X_i$.
    • Hệ số dẫn nhiệt hiệu dụng xây dựng theo mô hình song song của Murakami và Okos (1989): $\lambda = \sum V_i \lambda_i$.
    • Enthalpy hiệu dụng: $i = \sum \int C_{pi} X_i dt$.
  • Động học hô hấp: Tốc độ phản ứng sinh hóa của quá trình oxy hóa đường glucoza ($C_6H_{12}O_6 + 6O_2 \rightarrow 6CO_2 + 6H_2O + Q$) được mô hình hóa theo phương trình Arrhenius: $R_{CO2} = A \cdot \exp\left(-\frac{E_a}{R_\mu T}\right)$.
  • Cơ sở tính toán hao hụt tự nhiên: Dựa trên phân tích trạng thái nhiệt ẩm của không khí đi qua khối quả trong không gian bảo quản, thiết lập tia quá trình $\varepsilon = \frac{\Delta i}{\Delta d}$ trên đồ thị $i - d$ của không khí ẩm, xác định cân bằng ẩm giữa quả vải và môi trường bảo quản.

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

  • Phương pháp số: Sử dụng phương pháp thể tích hữu hạn (Finite Volume Method - FVM) để rời rạc hóa phương trình vi phân dẫn nhiệt từng lớp; giải hệ phương trình đại số tuyến tính bằng thuật toán ma trận ba đường chéo (TDMA).
  • Phương pháp nghịch đảo/gián tiếp xác định tham số: Ứng dụng thuật toán tối ưu hóa lát cắt vàng (Golden Section Search) để cực tiểu hóa hàm mục tiêu sai lệch nhiệt độ giữa mô phỏng và thực nghiệm, từ đó xác định chính xác hệ số trao đổi nhiệt đối lưu $\alpha$ và tọa độ thực $r$ của que đo nhiệt độ.
  • Phương pháp thực nghiệm:
    • Thu hái quả vải Lục Ngạn đúng độ chín 3, phân loại, đo kích thước hình học bằng thước kẹp điện tử để quy đổi đường kính tương đương ($D_{tđ}$).
    • Đo trường nhiệt độ làm lạnh sơ bộ trong bể nước đá và buồng không khí đối lưu cưỡng bức/tự nhiên bằng hệ thống cặp nhiệt điện chuẩn.
    • Đo nồng độ khí $O_2$ và $CO_2$ bằng máy phân tích khí Model 6600 trong bình kín đẳng nhiệt ở 5 mức nhiệt độ để xác định cường độ hô hấp.
    • Xác định khối lượng và độ hao hụt tự nhiên định kỳ bằng cân phân tích điện tử độ chính xác cao.
  • Nguồn tư liệu/dữ liệu: Các chuỗi số liệu thực nghiệm thu thập độc lập từ 3 vụ vải liên tiếp (2014, 2015, 2016) tại Phòng thí nghiệm Viện Khoa học & Công nghệ Nhiệt - Lạnh (Đại học Bách khoa Hà Nội), phối hợp đo kiểm chuẩn tại Viện Vệ sinh Dịch tễ Trung ương và Viện Cơ điện Nông nghiệp và Công nghệ sau thu hoạch.

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương 1 tổng hợp thực trạng ngành rau quả và chuỗi giá trị quả vải Việt Nam. Phân tích thành phần dinh dưỡng và đặc tính cơ bản của quả vải Lục Ngạn được tổng hợp qua các bảng số liệu:

Thành phần trong 100 g cùi vải Giá trị Chỉ tiêu vật lý - hóa học quả vải Lục Ngạn độ chín 3 Giá trị
Hàm lượng nước 81,9 ÷ 84,83 % Khối lượng trung bình 21,20 g (19,76 ± 2,48 g)
Năng lượng 63 ÷ 64 kcal Thể tích trung bình 18,52 ± 1,2 ml
Protein 0,68 ÷ 1,0 g Khối lượng riêng 1070 ± 44,72 $\text{kg}\cdot\text{m}^{-3}$
Chất béo 0,3 ÷ 0,58 g Đường kính trung bình 33,86 mm
Carbohydrate 13,31 ÷ 16,4 g Chiều cao trung bình 34,10 mm
Chất xơ 0,23 ÷ 0,4 g Hàm lượng đường tổng 18,00 %
Tro 0,37 ÷ 0,5 g Hàm lượng chất khô hòa tan 18,52 %
Vitamin C 24 ÷ 60 mg Hàm lượng axit 0,13 %

Chương này phân tích hạn chế của các phương pháp đo trực tiếp nhiệt dung riêng và hệ số dẫn nhiệt đối với quả tươi (do cấu trúc dị tính, kích thước lớp mỏng và sự sinh nhiệt hô hấp làm sai lệch dòng nhiệt ổn định), từ đó khẳng định tính đúng đắn của việc phát triển mô hình tính toán nhiệt vật lý lý thuyết kết hợp phương pháp số.

Chương 2: Nội dung và phương pháp nghiên cứu

Chương 2 trình bày chi tiết quy trình thu hái vải tại Lục Ngạn, quy cách phân loại quả đạt độ chín 3, phương pháp xác định đường kính tương đương $D_{tđ} = \sqrt[3]{D_h \cdot D_w \cdot D_d}$. Thiết kế và lắp đặt hệ thống thực nghiệm gồm:

  • Thiết bị làm lạnh sơ bộ bằng nước đá có tuần hoàn cưỡng bức.
  • Thiết bị làm lạnh bằng luồng không khí đối lưu có điều chỉnh vận tốc ($0 \div 2\text{ m}\cdot\text{s}^{-1}$) và nhiệt độ.
  • Hệ thống tủ bảo quản thử nghiệm có kiểm soát nhiệt độ (dải $0 \div 10^\circ\text{C}$, sai số $\pm 0,5^\circ\text{C}$) và độ ẩm (đạt $90 \div 95%$).
  • Thiết lập quy trình phân tích sai số thí nghiệm, đánh giá độ ổn định nhiệt ẩm môi trường không tải và có tải trước khi thực hiện các phép đo chính thức.

Chương 3: Mô phỏng quá trình làm lạnh

Chương 3 phát triển các phương trình tính toán thông số nhiệt vật lý cho từng lớp (vỏ, cùi, hạt) và toàn bộ quả vải:

  • Khối lượng riêng hiệu dụng $\rho$, nhiệt dung riêng đẳng áp $C_p(t)$, hệ số dẫn nhiệt $\lambda(t)$ và enthalpy $i(t)$.
  • Xây dựng sơ đồ thuật toán và chương trình tính toán phân bố trường nhiệt độ theo bán kính và thời gian.
  • Kiểm chứng độ tin cậy của mô hình số bằng cách đối chiếu với dữ liệu thực nghiệm của Dincer (1995) đối với vật thể hình cầu và mô hình của Olesen (2004).
  • Kết quả mô phỏng quá trình làm lạnh sơ bộ: Làm lạnh quả vải có $D_{tđ} = 32\text{ mm}$ từ $31,31^\circ\text{C}$ xuống $7,05^\circ\text{C}$ bằng nước đá $4,28^\circ\text{C}$ (tốc độ $0,2\text{ m}\cdot\text{s}^{-1}$) diễn ra nhanh gấp nhiều lần so với làm lạnh bằng không khí đối lưu $4,27^\circ\text{C}$ (tốc độ $0,5\text{ m}\cdot\text{s}^{-1}$).
  • Ứng dụng phương pháp tối ưu phiếm hàm sai số xác định gián tiếp hệ số trao đổi nhiệt đối lưu:
    • Quá trình làm lạnh trong nước đá: hiệu chỉnh hệ số $\alpha$ và vị trí que đo thực tế đạt độ tương thích cao với thực nghiệm.
    • Quá trình làm lạnh trong không khí đối lưu tự nhiên: xác định chính xác giá trị $\alpha$ tại bề mặt quả, giải quyết sai số vị trí que đo khi cắm vào quả vải thực nghiệm.

Chương 4: Xây dựng mô hình xác định độ hao hụt tự nhiên và thông số công nghệ lạnh cho quả vải

Chương 4 hoàn thiện các mô hình công nghệ phục vụ thiết kế:

  • Xây dựng mô hình cường độ hô hấp $R_{CO2}$ theo nhiệt độ từ số liệu đo thực nghiệm 3 vụ (2014 ÷ 2016) ở 5 mức nhiệt độ ($0^\circ\text{C}, 5^\circ\text{C}, 10^\circ\text{C}, 15^\circ\text{C}, 20^\circ\text{C}$), so sánh và kiểm chứng với mô hình của Paull.
  • Xây dựng phần mềm chuyên dụng tính toán nhiệt tỏa do hô hấp $q_{resp}$ của quả vải.
  • Xây dựng mô hình xác định độ hao hụt khối lượng tự nhiên do bay hơi nước dựa trên phương trình cân bằng ẩm và tia quá trình $\varepsilon$ của không khí trong dải nhiệt độ $0 \div 25^\circ\text{C}$.
  • Kiểm chứng mô hình độ hao hụt tự nhiên qua thí nghiệm bảo quản thực tế tại nhiệt độ $+4,5^\circ\text{C}$, độ ẩm $92%$ RH.
  • Xác định các thông số định mức công nghệ lạnh:
    • Suất tiêu hao đá cây cho quá trình làm lạnh sơ bộ tại vùng nguyên liệu.
    • Suất tiêu hao đá bổ sung trong quá trình đóng thùng cách nhiệt vận chuyển đường dài.
    • Năng suất lạnh cần thiết của thiết bị tái lạnh tại kho bảo quản tập trung.
  • Đề xuất 3 mô hình kho lạnh bảo quản quả vải quy mô công nghiệp: kho lạnh có dàn quạt thổi cưỡng bức trực tiếp, kho lạnh có hệ thống phân phối không khí kiểu kênh dẫn khí, và kho lạnh có hệ thống phân phối không khí kiểu ống gió đục lỗ đồng đều.

Chương 5: Kết luận và kiến nghị

Tổng kết toàn bộ các kết quả đạt được về mặt khoa học lý thuyết và giá trị thực tiễn; đề xuất các hướng nghiên cứu mở rộng trong giai đoạn tiếp theo.

Kết quả và những đóng góp mới

Đóng góp mới về mặt khoa học lý luận

  1. Xây dựng mô hình mô phỏng quá trình làm lạnh của quả vải trên cơ sở hệ phương trình vi phân dẫn nhiệt không ổn định viết cho hệ tọa độ cầu 1 chiều với điều kiện biên loại 3, áp dụng cho đối tượng không đồng chất nhưng đẳng hướng. Phát triển phương pháp tính toán tính chất nhiệt vật lý hiệu dụng ($C_p, \lambda, \rho, i$) bằng cách xem quả vải là hệ nhiệt động cân bằng đa thành phần không biến pha dựa trên thành phần của các lớp vỏ, cùi và hạt.
  2. Xây dựng mô hình cường độ hô hấp của quả vải Lục Ngạn phụ thuộc vào nhiệt độ theo phương trình Arrhenius, dựa trên việc xử lý dữ liệu thực nghiệm của hơn 20 mẻ thí nghiệm trong 3 mùa vụ liên tiếp (2014 ÷ 2016) ở 5 chế độ nhiệt độ.
  3. Phát triển thành công mô hình xác định độ hao hụt tự nhiên của quả vải trong quá trình bảo quản dựa trên nguyên lý hệ nhiệt động tựa cân bằng đẳng áp và phương pháp tia quá trình $\varepsilon$ trên đồ thị nhiệt ẩm $i - d$.
  4. Thiết lập phương pháp gián tiếp kết hợp mô hình số và thuật toán tối ưu hóa lát cắt vàng để xác định đồng thời hệ số trao đổi nhiệt đối lưu bề mặt $\alpha$ và vị trí thực tế của đầu đo nhiệt độ bên trong quả vải.

Đóng góp mới về mặt thực tiễn và ứng dụng

  1. Xây dựng phần mềm tính toán nhiệt tỏa do hô hấp của quả vải theo nhiệt độ, làm công cụ tính toán tải nhiệt thiết kế hệ thống lạnh.
  2. Xác định các chỉ tiêu định mức công nghệ lạnh trong chuỗi cung ứng: suất tiêu hao đá cây làm lạnh sơ bộ, lượng đá duy trì khi vận chuyển bằng thùng cách nhiệt, năng suất lạnh tái lạnh khi nhập kho.
  3. Đề xuất các yêu cầu kỹ thuật và thông số công nghệ thiết kế các mẫu kho lạnh bảo quản quả vải quy mô tập trung (chế độ tối ưu: nhiệt độ $2 \div 5^\circ\text{C}$, độ ẩm $90 \div 95%$, thời gian bảo quản đạt từ 30 ngày trở lên đảm bảo tiêu chuẩn thương phẩm).

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án:

  • Mô hình mô phỏng truyền nhiệt giả thiết quả vải có dạng hình cầu đối xứng lý tưởng, chưa bao quát đầy đủ sự sai lệch hình học bất đối xứng của các quả vải thực tế.
  • Các khảo sát thực nghiệm tập trung chuyên biệt vào giống vải thiều Lục Ngạn ở độ chín 3, chưa mở rộng sang các độ chín khác hoặc các giống vải khác nhau (như vải tu hú, vải lai chín sớm).
  • Quá trình khuếch tán ẩm trong mô hình hao hụt tự nhiên được tiếp cận theo phương pháp vĩ mô qua trạng thái không khí kho bảo quản, chưa tích hợp hệ phương trình vi phân truyền ẩm 3 chiều bên trong cấu trúc mao dẫn của vỏ quả.

Hướng nghiên cứu tiếp theo được tác giả đề xuất:

  • Phát triển mô hình mô phỏng truyền nhiệt và truyền chất đồng thời 3 chiều có xét đến hình dạng thực phi tiêu chuẩn của quả vải bằng kỹ thuật quét 3D và phương pháp phần tử hữu hạn (FEM).
  • Nghiên cứu sâu hơn về biến đổi vi sinh vật, hoạt tính enzym và các chỉ tiêu hóa sinh chi tiết trong điều kiện bảo quản kết hợp công nghệ lạnh với bao gói khí điều biến quy mô công nghiệp.
  • Hoàn thiện gói phần mềm tự động hóa tính toán thiết kế kho lạnh chuyên dụng cho chuỗi bảo quản nông sản nhiệt đới.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu trong các lĩnh vực:

  • Giảng viên, nghiên cứu sinh và học viên ngành Kỹ thuật Lạnh, Công nghệ Nhiệt, Máy và Thiết bị Năng lượng: Tham khảo phương pháp xây dựng mô hình số giải bài toán truyền nhiệt không ổn định, kỹ thuật xác định tính chất nhiệt vật lý hiệu dụng của thực phẩm dị tính, và thuật toán tối ưu hóa gián tiếp tham số truyền nhiệt.
  • Các nhà nghiên cứu công nghệ sau thu hoạch và chế biến nông sản: Sử dụng dữ liệu thực nghiệm về cường độ hô hấp, nhiệt hô hấp và động học mất nước của quả vải để thiết kế quy trình bảo quản tối ưu.
  • Kỹ sư thiết kế, doanh nghiệp chế biến và cơ quan quản lý nông nghiệp: Ứng dụng các thông số định mức suất tiêu hao đá, tải nhiệt tái lạnh và sơ đồ công nghệ phân phối không khí trong kho lạnh để đầu tư xây dựng cơ sở hạ tầng chuỗi lạnh tại các vùng nguyên liệu tập trung (như Bắc Giang, Hải Dương).

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

1. Tại sao tác giả không sử dụng phương pháp đo thực nghiệm trực tiếp để xác định nhiệt dung riêng ($C_p$) và hệ số dẫn nhiệt ($\lambda$) của quả vải?
Trả lời: Do quả vải có cấu trúc dị tính phức tạp gồm 3 lớp riêng biệt (vỏ mỏng $0,8 \div 1,5\text{ mm}$, cùi mọng nước chứa $>80%$ ẩm, hạt cứng), không đồng chất và đẳng hướng. Đo trực tiếp bằng phương pháp đĩa nóng hay que thăm nhiệt đòi hỏi mẫu khô, đồng nhất và thời gian đo kéo dài hàng chục giờ; khi đó nhiệt tỏa do hô hấp và sự mất nước sẽ làm sai lệch hoàn toàn kết quả đo. Vì vậy, tác giả sử dụng phương pháp mô hình hóa tính toán bán thực nghiệm theo tỷ lệ thành phần hóa học kết hợp mô hình song song của Murakami & Okos.

2. Phương pháp gián tiếp xác định hệ số trao đổi nhiệt đối lưu ($\alpha$) và vị trí que đo nhiệt độ được thực hiện như thế nào?
Trả lời: Tác giả kết hợp mô hình mô phỏng FVM với dữ liệu thực nghiệm đo nhiệt độ theo thời gian. Bằng cách thiết lập hàm mục tiêu phiếm hàm $f$ đánh giá độ lệch chuẩn giữa nhiệt độ mô phỏng và nhiệt độ đo thực tế, thuật toán lát cắt vàng được dùng để tối ưu hóa đồng thời giá trị hệ số $\alpha$ và vị trí bán kính $r$ thực tế của que đo trong quả vải, khắc phục sai số định vị cơ học khi cắm que đo vào tâm quả.

3. Mối quan hệ giữa cường độ hô hấp của quả vải và nhiệt độ được mô hình hóa theo dạng toán học nào?
Trả lời: Cường độ hô hấp ($R_{CO2}$) được mô hình hóa theo phương trình động học tốc độ phản ứng Arrhenius:

$$R_{CO2} = A \cdot \exp\left(-\frac{E_a}{R_\mu T}\right)$$

Dựa trên dữ liệu thực nghiệm ở 5 mức nhiệt độ ($0^\circ\text{C}, 5^\circ\text{C}, 10^\circ\text{C}, 15^\circ\text{C}, 20^\circ\text{C}$) qua 3 mùa vụ liên tiếp (2014 ÷ 2016), mô hình phản ánh sự sụt giảm phi tuyến của lượng $CO_2$ sinh ra khi nhiệt độ bảo quản hạ thấp.

4. Sự khác biệt về hiệu quả làm lạnh sơ bộ giữa môi trường nước đá và không khí đối lưu là gì?
Trả lời: Thí nghiệm và mô phỏng chỉ ra làm lạnh sơ bộ bằng nước đá ($4,28^\circ\text{C}$, tốc độ $0,2\text{ m}\cdot\text{s}^{-1}$) hạ nhiệt độ tâm quả từ $31,31^\circ\text{C}$ xuống $7,05^\circ\text{C}$ chỉ mất khoảng $15 \div 20$ phút và tỷ lệ mất nước bằng 0%. Trong khi đó, làm lạnh bằng không khí đối lưu ($4,27^\circ\text{C}$, tốc độ $0,5 \div 2\text{ m}\cdot\text{s}^{-1}$) cần thời gian kéo dài từ $1 \div 3$ giờ và làm mất $0,20 \div 0,38%$ khối lượng do bay hơi nước.

5. Các thông số công nghệ tối ưu được luận án khuyến nghị để bảo quản lạnh quả vải tươi là gì?
Trả lời: Chế độ bảo quản tối ưu trong kho lạnh là nhiệt độ duy trì trong dải $2 \div 5^\circ\text{C}$, độ ẩm tương đối $\phi = 90 \div 95%$. Quả vải sau thu hoạch cần được làm lạnh sơ bộ bằng nước đá trong vòng 4 giờ đầu, sau đó bảo quản trong kho lạnh có hệ thống phân phối gió cưỡng bức đồng đều (kiểu kênh hoặc ống đục lỗ) để kéo dài thời gian thương phẩm trên 30 ngày.

Kết luận

Luận án của tác giả Lê Xuân Tuấn đã giải quyết trọn vẹn cơ sở khoa học và tính toán công nghệ cho chuỗi bảo quản lạnh quả vải thiều tươi quy mô tập trung. Bằng việc kết hợp chặt chẽ giữa mô hình hóa nhiệt động lực học, mô phỏng số truyền nhiệt truyền chất và chuỗi thực nghiệm 3 mùa vụ, công trình đã xây dựng thành công các mô hình xác định tính chất nhiệt vật lý hiệu dụng, động học hô hấp và độ hao hụt tự nhiên. Các kết quả tính toán định mức suất tiêu hao đá, tải nhiệt tái lạnh và thiết kế kho bảo quản cung cấp luận cứ kỹ thuật phục vụ trực tiếp ngành cơ điện lạnh và công nghệ sau thu hoạch tại Việt Nam.