Tổng quan về luận án

Ngành thủy sản đóng vai trò then chốt trong cơ cấu kinh tế biển Việt Nam. Dữ liệu từ Hiệp hội Chế biến và Xuất khẩu Thủy sản Việt Nam (VASEP) khẳng định: "Năm 2019 xuất khẩu tôm của Việt Nam ước đạt 3,38 tỷ USD. Trong đó, tôm thẻ chân trắng đông lạnh xuất khẩu đạt 2,36 tỷ USD chiếm 70% giá trị tôm xuất khẩu, mặt hàng tôm khô chế biến cũng đã tăng 40,7% so với cùng kỳ năm 2018". Mặc dù tiềm năng thương mại rất lớn, ngành công nghiệp chế biến tôm khô nội địa đối mặt với rào cản công nghệ nghiêm trọng: phương pháp phơi nắng tự nhiên phụ thuộc thời tiết kéo dài thời gian làm khô, trong khi các lò sấy thủ công đốt than, củi hoặc dùng điện trở nhiệt cao làm biến tính protein, sẫm màu sản phẩm và gia tăng phát thải khí nhà kính.

┌─────────────────────────────────────────────────────────────────────────┐
│              KHUNG TRUYỀN NHIỆT & TRUYỀN CHẤT LAI HÓA (HP-IR)           │
├──────────────────────────┬─────────────────────────┬────────────────────┤
│   Bức xạ Hồng ngoại (IR) │   Bơm nhiệt (HP Air)    │ Gradient Liên kết  │
├──────────────────────────┼─────────────────────────┼────────────────────┤
│ • Bước sóng: 2,7 - 25 µm │ • Tách ẩm dàn bay hơi   │ • Nhiệt (Fourier)  │
│ • Đâm xuyên & gia nhiệt  │ • Không khí khô lướt    │ • Khuếch tán (Fick)│
│   tâm mao quản tôm       │   bề mặt ngoài          │ • Áp suất (Darcy)  │
└──────────────────────────┴─────────────────────────┴────────────────────┘

Điểm nghẽn nghiên cứu (Research Gap) cốt lõi trong lý thuyết sấy thủy sản hiện nay là sự thiếu vắng các mô hình động học liên kết 3 trường: nhiệt độ ($t$), độ ẩm ($M$) và áp suất mao quản ($P$) có tính đến biến dạng co rút thể tích ($S_v$) và tính bất đẳng hướng cấu trúc cơ thịt tôm. Các công trình trước đây thường chỉ phân tích đơn lẻ quá trình truyền nhiệt đối lưu hoặc khuếch tán ẩm đẳng nhiệt theo định luật Fick, bỏ qua gradien áp suất nội tại do bức xạ hồng ngoại tạo ra.

Luận án tiến sĩ kỹ thuật của tác giả Lê Như Chính (2020) tại Trường Đại học Bách khoa Hà Nội, dưới sự hướng dẫn của PGS.TS. Nguyễn Nguyên An và GS. Phạm Văn Tùy, đã giải quyết khoảng trống này thông qua các câu hỏi và giả thuyết nghiên cứu tường minh:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế tương tác động học giữa trường nhiệt do bức xạ hồng ngoại bước sóng ngắn/trung bình và trường áp suất mao quản bên trong cơ thịt tôm diễn ra như thế nào?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật biến thiên của các thông số nhiệt vật lý ($\lambda, C, \rho, a$) và thông số sinh học đặc trưng ($\varepsilon_p, K, D_{eff}$) phụ thuộc vào nhiệt độ, hàm lượng ẩm và thành phần hóa học (protein, lipid, khoáng) theo không gian và thời gian sấy ra sao?
  • Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để thiết lập chế độ sấy đa mục tiêu tối ưu trên hệ thống sấy bơm nhiệt kết hợp hồng ngoại (HP-IR) nhằm đồng thời cực tiểu hóa thời gian sấy ($\tau$), suất tiêu hao năng lượng riêng ($SEC$), và cực đại hóa chất lượng cảm quan ($CLCQ$), tỷ lệ hút nước phục hồi ($R_{TN}$)?

Các giả thuyết tương ứng gồm:

  • Giả thuyết 1 (H1): Năng lượng bức xạ hồng ngoại xuyên thấu vào lõi vật liệu kết hợp với tác nhân sấy khô nhiệt độ thấp từ bơm nhiệt tạo ra gradient nhiệt độ và gradient áp suất cùng chiều hướng tâm ra bề mặt, kích hoạt dòng ẩm đối lưu Darcy song song với dòng khuếch tán Fick, rút ngắn ít nhất 40% thời gian sấy so với bơm nhiệt đơn lẻ.
  • Giả thuyết 2 (H2): Mô hình hóa số 2D/3D phần tử hữu hạn (FEM) tích hợp đầy đủ hệ phương trình Luikov biến tính có thể dự báo chính xác trường nhiệt - ẩm - áp với sai số tương đối RMS dưới 10% so với thực nghiệm.
  • Giả thuyết 3 (H3): Tối ưu hóa bề mặt đáp ứng (RSM) kết hợp phương pháp quy hoạch thực nghiệm Taguchi sẽ xác lập được vùng vận hành tối ưu mang lại suất tiêu hao năng lượng $SEC \le 4,0 \text{ kWh/kg } H_2O$ trong khi duy trì hàm lượng $NH_3 \le 35 \text{ mg/100g}$.

Khung lý thuyết của công trình tích hợp hệ thống lý thuyết truyền nhiệt dẫn nhiệt Fourier có nguồn nhiệt thể tích $q_{v,IR}$, lý thuyết khuếch tán Fick bậc 2, động lực học dòng môi chất xốp Darcy và hệ phương trình vi phân liên kết Luikov cho hệ mao dẫn - xốp. Nghiên cứu giới hạn phạm vi trên đối tượng tôm thẻ chân trắng Việt Nam (Litopenaeus vannamei thuộc họ Penaeidae), cỡ nguyên liệu thương phẩm tiêu chuẩn (đường kính thân $d_{TB1} \approx 12 - 15 \text{ mm}$, độ ẩm ban đầu $W_1 = 76,63%$), thực hiện trên tổ hợp máy sấy HP-IR quy mô phòng thí nghiệm và bán công nghiệp.

Literature Review và Positioning

Các khuynh hướng nghiên cứu công nghệ sấy sinh khối và thực phẩm trên thế giới và Việt Nam phân hóa rõ rệt thành ba giai đoạn tiến hóa công nghệ:

┌──────────────────────────────────────────────────────────────────────────┐
│                   TIẾN TRÌNH TIẾN HÓA CÔNG NGHỆ SẤY                      │
│                                                                          │
│  Sấy Đối lưu Cổ điển ──► Sấy Bơm nhiệt Đơn lẻ ──► Sấy Lai hóa Đa trường  │
│  (Hot Air Convection)    (Heat Pump Drying)       (HP-IR / Multiphysics) │
│  • Nhiệt độ: 70-120°C    • Nhiệt độ: 35-45°C      • Gia nhiệt chọn lọc   │
│  • Phá hủy protein       • Thời gian dài (10-12h) • Tối ưu năng lượng    │
│  • Tiêu tốn năng lượng   • Không có gradien áp    • Bảo toàn vi chất     │
└──────────────────────────────────────────────────────────────────────────┘
  1. Dòng nghiên cứu sấy bơm nhiệt đơn lẻ: Các công trình của Dirk et al. (2014), Strommen (2012) và Levent Taşeri (2018) chứng minh sấy bơm nhiệt (HP) tiết kiệm từ 60% đến 80% năng lượng tiêu thụ so với sấy không khí nóng đối lưu tiêu chuẩn nhờ thu hồi ẩn nhiệt ngưng tụ tại dàn nóng. Tại Việt Nam, nghiên cứu của Võ Duy Mạnh và Lê Chí Hiệp (2012) trên cà rốt, Phạm Văn Tùy et al. (2014) trên kẹo dẻo jelly và mực ống lột da (2015) khẳng định ưu thế nhiệt độ thấp ($35 - 45^\circ\text{C}$), giữ màu sắc và hoạt chất sinh học. Tuy nhiên, nhược điểm chí mạng của HP thuần túy là thời gian sấy kéo dài (thường vượt quá $10 - 12\text{ giờ}$ đối với giáp xác), dẫn đến năng suất thiết bị thấp và tăng nguy cơ tái nhiễm vi sinh ở pha sấy giảm tốc ban đầu.
  2. Dòng nghiên cứu sấy bức xạ hồng ngoại: Rosenthal (1992) và Matsuura (2001) làm rõ cơ chế bức xạ hồng ngoại (IR): khi nước hấp thụ năng lượng ở dải bước sóng tối ưu $\lambda = 2,7 - 25\ \mu\text{m}$, "Cường độ hấp thụ của các liên kết hóa trị theo tỷ lệ $v_1, v_2, v_3$ tương ứng là 0,07, 1,47 và 1,00. Kết quả là các lực ma sát sinh ra rất mạnh dẫn tới nhiệt độ của vật bị chiếu IR tăng lên một cách nhanh chóng". Nghiên cứu của Supawan Tirawanichakul et al. (2008) và Ratti (2011) chỉ ra sấy IR đơn lẻ rút ngắn thời gian 2 đến 3 lần, nhưng nhược điểm lớn được Nguyễn Thị Bích Thủy (2010) cảnh báo là mật độ dòng nhiệt bề mặt quá cao dễ gây cháy cục bộ sém vỏ tôm, cong vênh và nứt nẻ cơ thịt do chênh lệch nhiệt ẩm quá mức.
  3. Dòng nghiên cứu sấy lai ghép HP-IR và mô hình hóa: Mustafa Aktas (2016, 2017) ghi nhận tổ hợp HP-IR cải thiện hiệu suất nhiệt từ 25% lên 39% và giảm tiêu thụ điện 43,2%. Yun Deng et al. (2014) và Xiaoyong Song et al. (2018) chứng minh HP-IR làm tăng hệ số khuếch tán ẩm $D_{eff}$, hạn chế phản ứng tạo màu nâu Maillard và duy trì nồng độ acid amin thiết yếu.

Về mặt mô hình hóa toán học, tồn tại cuộc tranh luận học thuật sâu sắc:

  • Trường phái Fick thuần túy đại diện bởi Burubai (2017), Worachard Chawanasporn et al. (2009) giả định dòng ẩm di chuyển trong mô tôm hoàn toàn do khuếch tán pha lỏng đẳng áp điều khiển bởi gradient nồng độ ẩm $\nabla M$.
  • Trường phái động lực học môi chất xốp liên kết khởi xướng bởi A. Luikov (1968), Brooker (1992) khẳng định sự tồn tại đồng thời của dòng khuếch tán Fick và dòng đối lưu chất lỏng/khí do gradient áp suất $\nabla P$ điều khiển theo định luật Darcy.
Tiêu chí So sánh Nghiên cứu Chalida Niamnuya et al. (Thái Lan) Nghiên cứu Worachard et al. (Thái Lan) Luận án Lê Như Chính (ĐHBK Hà Nội, Việt Nam)
Công nghệ sấy Sấy tầng sôi không khí nóng ($80 - 120^\circ\text{C}$) Sấy đối lưu không khí nóng ($70^\circ\text{C}$) Sấy bơm nhiệt kết hợp hồng ngoại HP-IR ($45^\circ\text{C}$)
Mô hình toán học Truyền nhiệt Fourier + Khuếch tán Fick 3D Hệ vi phân 2D bỏ qua gradien áp suất Hệ liên kết Fourier ($q_{v,IR}$) + Fick + Darcy ($\nabla P$)
Thông số nhiệt vật lý Hàm phụ thuộc đơn thuần theo độ ẩm $M$ Hằng số thực nghiệm cố định Đa biến phụ thuộc nhiệt độ $t$ và thành phần sinh hóa
Đo đạc áp suất nội Không thực hiện (bỏ qua $P$) Không thực hiện (bỏ qua $P$) Trực tiếp bằng áp kế chữ U + kim y tế chuyên dụng
Hiệu quả năng lượng Tiêu hao năng lượng cao ($> 8 \text{ kWh/kg}$) Hiệu suất nhiệt thấp ($< 25%$) Tối ưu hóa đa mục tiêu ($SEC \approx 3,84 \text{ kWh/kg}$)

Công trình của Lê Như Chính đã định vị xuất sắc đóng góp học thuật bằng cách bác bỏ giả định bỏ qua áp suất trong sấy tôm, chứng minh bằng thực nghiệm và mô phỏng rằng trong sấy HP-IR, áp suất nội tại đạt đỉnh đáng kể do hiện tượng bốc hơi nhanh bên trong các mao quản kín, đóng vai trò lực đẩy chủ đạo gia tốc dòng ẩm thoát ra bề mặt.

Đóng góp lý thuyết và khung phân tích

                      ┌────────────────────────────────────────┐
                      │    HỆ THỐNG PHƯƠNG TRÌNH LIÊN KẾT      │
                      └──────────────────┬─────────────────────┘
                                         │
                 ┌───────────────────────┼──────────────────────┐
                 ▼                       ▼                      ▼
      ┌────────────────────┐  ┌────────────────────┐  ┌────────────────────┐
      │  TRUYỀN NHIỆT      │  │  TRUYỀN CHẤT       │  │  ÁP SUẤT NỘI       │
      │  (Fourier + IR)    │  │  (Fick + Darcy)    │  │  (Biến thiên pha)  │
      │                    │  │                    │  │                    │
      │ ρC ∂t/∂τ =         │  │ ∂M/∂τ =            │  │ ∂P/∂τ =            │
      │ ∇·(λ∇t) + q_v,IR   │  │ ∇·(D_eff∇M) -      │  │ f(∇·u_v, ∇·u_l,    │
      │                    │  │ ∇·(u_l M_l + u_v)  │  │   ΔT, S_v)         │
      └────────────────────┘  └────────────────────┘  └────────────────────┘

Đóng góp cho lý thuyết

Luận án mở rộng toàn diện lý thuyết truyền nhiệt truyền chất trong môi trường xốp sinh học dị hướng của A. Luikov bằng cách thiết lập mô hình toán học giải tích - số kết hợp nguồn nhiệt bức xạ đâm xuyên và động lực học dòng ẩm hai pha (lỏng - hơi).

Hệ phương trình vi phân chủ đạo được thiết lập:

Phương trình truyền nhiệt có nguồn bức xạ: $$\rho C \frac{\partial t}{\partial \tau} = \nabla \cdot (\lambda \nabla t) + q_{v,IR}(x, y, z, \tau)$$

Trong đó, mật độ nguồn nhiệt thể tích từ bức xạ hồng ngoại suy giảm theo định luật Beer-Lambert mở rộng: $$q_{v,IR} = \alpha I_0 e^{-\alpha \delta}$$ với $\alpha$ là hệ số hấp thụ quang học tổng thể, $\delta$ là độ sâu thâm nhập của chùm tia hồng ngoại vào tổ chức mô cơ thịt tôm.

Phương trình truyền chất bảo toàn khối lượng ẩm hai pha: $$\frac{\partial M}{\partial \tau} = \nabla \cdot \left( D_{eff} \nabla M \right) - \nabla \cdot \left( \vec{u}_l \rho_l + \vec{u}_v \rho_v \right)$$

Vận tốc dòng ẩm đối lưu mao quản tuân theo định luật Darcy cho môi trường xốp: $$\vec{u}l = -\frac{K \cdot K{r,l}}{\mu_l} \nabla P ; \quad \vec{u}v = -\frac{K \cdot K{r,v}}{\mu_v} \nabla P$$

Hệ thống mệnh đề lý thuyết (Theoretical Propositions) được xác lập:

  • Mệnh đề P1: Gradien nhiệt độ nghịch đảo ($\nabla t > 0$ từ bề mặt vào tâm) tạo ra bởi bức xạ hồng ngoại kết hợp làm mát bề mặt bằng gió khô sẽ sinh ra gradien áp suất hơi bão hòa $\nabla P > 0$, biến đổi cơ chế dịch chuyển ẩm từ thuần khuếch tán phân tử sang cơ chế đồng vận tải khuếch tán - áp suất.
  • Mệnh đề P2: Hệ số dẫn nhiệt của tôm thể hiện tính bất đẳng hướng cấu trúc rõ rệt giữa phương song song ($\lambda_{//}$) và phương vuông góc ($\lambda_\perp$) với bó sợi cơ actin-myosin, quy định bởi hàm phân bố thành phần hóa học: $$\lambda = \sum_{i} \lambda_i X_i = \lambda_w X_w + \lambda_p X_p + \lambda_f X_f + \lambda_a X_a$$

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

Khung phân tích của luận án tích hợp liên ngành 4 trụ cột lý thuyết:

  1. Nhiệt động học kỹ thuật và truyền nhiệt bức xạ (bước sóng chọn lọc của liên kết $O-H$ trong nước thực phẩm);
  2. Cơ học môi trường xốp sinh học biến dạng (Porous media mechanics with matrix shrinkage);
  3. Hóa sinh học thủy sản thực nghiệm (động học thoái biến protein và chỉ số bazơ bay hơi);
  4. Lý thuyết tối ưu hóa toán học đa mục tiêu (Multi-objective Response Surface Methodology kết hợp hàm mong muốn Desirability Function).

Điều kiện biên thực tế được xây dựng chi tiết hóa:

  • Biên truyền nhiệt bề mặt: $$-\lambda \left( \frac{\partial t}{\partial n} \right)s = \alpha{dl} (t_s - t_k) + \sigma \varepsilon_s \varepsilon_{IR} F_{s,IR} (T_s^4 - T_{IR}^4) + r \alpha_m (M_s - M_e)$$
  • Biên truyền chất bề mặt: $$-D_{eff} \rho_d \left( \frac{\partial M}{\partial n} \right)_s = \alpha_m \rho_k (M_s - M_e)$$

Điều kiện biên giới hạn (Boundary conditions) xác định rõ: mô hình áp dụng chính xác cho tôm nguyên con đã qua gia nhiệt sơ bộ (luộc chín để ổn định hình học ban đầu), với nhiệt độ tác nhân sấy $t_k \in [35^\circ\text{C}, 55^\circ\text{C}]$, vận tốc dòng khí $V \in [1,5 \text{ m/s}, 3,5 \text{ m/s}]$, mật độ công suất bức xạ $IP \in [250\text{W}, 750\text{W}]$.

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

┌────────────────────────────────────────────────────────────────────────┐
│                      QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU                  │
├──────────────────────────┬──────────────────────┬──────────────────────┤
│ 1. Hình thái & Cấu trúc  │ 2. Mô phỏng Đa trường│ 3. Đo đạc & Tối ưu   │
├──────────────────────────┼──────────────────────┼──────────────────────┤
│ • Xử lý ảnh ImageJ 2D    │ • COMSOL Multiphysics│ • Đầu dò áp kế chữ U │
│ • Độ xốp ε_p = 0,26-0,38 │ • Lưới FEM tự do     │ • Quy hoạch Taguchi  │
│ • Độ co rút thể tích S_v │ • Bộ giải PARISO/BDF │ • Tối ưu hóa RSM/ANOVA│
└──────────────────────────┴──────────────────────┴──────────────────────┘

Thiết kế nghiên cứu

Luận án tuân thủ triết lý nghiên cứu thực chứng quy nạp kết hợp diễn dịch (Positivist / Post-positivist empirical-deductive paradigm). Thiết kế nghiên cứu đa tầng kết nối chặt chẽ giữa mô hình hóa toán học lý thuyết, mô phỏng số trên máy tính và kiểm chứng thực nghiệm đo đạc vật lý.

Mẫu nghiên cứu là tôm thẻ chân trắng (Penaeus vannamei) thu mua đồng nhất tại đầm nuôi ven biển Khánh Hòa - Bình Định, bảo quản lạnh sống, tuyển chọn khối lượng $15 \pm 1 \text{ g/con}$, đường kính gốc đầu $d_{shr1} = 14,2 \pm 0,5 \text{ mm}$, đường kính đuôi $d_{shr2} = 7,1 \pm 0,3 \text{ mm}$, chiều dài thân $L_{shr} = 95 \pm 3 \text{ mm}$. Tiêu chí loại trừ nghiêm ngặt: tôm dập vỏ, đứt đốt, biến màu hoặc có dấu hiệu phân hủy protein nội sinh.

Quy trình nghiên cứu rigorous

Quy trình thu thập dữ liệu và chế tạo thiết bị đo đặc thù bao gồm:

  1. Thiết kế chế tạo hệ thống sấy tích hợp HP-IR: Tủ sấy thể tích hữu ích $0,8 \text{ m}^3$, tích hợp cụm bơm nhiệt môi chất R134a, máy nén khí kín, dàn bay hơi tách ẩm và dàn ngưng gia nhiệt gió; kết hợp trần đèn hồng ngoại gốm bước sóng $\lambda = 2 - 10\ \mu\text{m}$ điều khiển công suất biến thiên tuyến tính qua SCR (Silicon Controlled Rectifier).
  2. Sáng chế bộ cảm biến đo áp suất mao quản vi mô: Tác giả đã chế tạo bộ đầu dò áp suất độc đáo sử dụng cụm kim tiêm truyền dịch y tế đường kính trong $0,4 \text{ mm}$, điền đầy dầu silicon trung tính truyền áp, liên kết trực tiếp với áp kế vi sai chữ U chứa chất lỏng đo tỷ trọng chuẩn, cắm trực tiếp vào tâm cơ thịt tôm tại đốt thân thứ 3.
  3. Phân tích vi cấu trúc và độ xốp: Cắt lát tiêu bản mô tôm ở các mốc thời gian sấy khác nhau, chụp ảnh hiển vi quang học độ phân giải cao, số hóa và xử lý ảnh nhị phân tự động trên phần mềm ImageJ để định lượng chính xác diện tích lỗ xốp ($A_p$), diện tích chất khô ($A_{sd}$), từ đó trích xuất hàm biến thiên độ xốp thực tế $\varepsilon_p(\tau)$.
   Áp suất mao quản tôm (ΔP)
          │
          ├──► Đầu kim truyền dịch y tế (Ø 0.4mm)
          │         │
          │         ▼
          └──► Dầu Silicon dẫn áp ──► Áp kế chữ U vi sai (±1.0 Pa)

Data và phân tích

Mô phỏng trường liên kết đa vật lý được lập trình và giải trên phần mềm thương mại COMSOL Multiphysics với mô-đun PDE (Partial Differential Equations) ghép nối với Heat Transfer và Porous Media Flow. Không gian hình học thực tế của tôm được quét dựng 2D trục đối xứng và 3D CAD. Lưới phần tử hữu hạn (FEM Mesh) dạng tam giác phi cấu trúc với 18.450 phần tử, làm mịn cục bộ tại lớp biên bề mặt với bước thời gian tích phân thích nghi $\Delta \tau = 1,0 \text{ s}$.

Phân tích tối ưu hóa thực nghiệm sử dụng phần mềm Design-Expert với quy hoạch trực giao Taguchi $L_9$ và phương pháp bề mặt đáp ứng (RSM) dạng Box-Behnken 4 yếu tố 3 mức:

  • $X_1$: Nhiệt độ tác nhân sấy ($t_k$: $35, 45, 55^\circ\text{C}$)
  • $X_2$: Vận tốc gió sấy ($V$: $1,5; 2,5; 3,5 \text{ m/s}$)
  • $X_3$: Khoảng cách đèn hồng ngoại ($H$: $20, 30, 40 \text{ cm}$)
  • $X_4$: Công suất đèn hồng ngoại ($IP$: $250, 500, 750 \text{ W}$)

Các mô hình hồi quy đa thức bậc hai được kiểm định độ tin cậy bằng phân tích phương sai ANOVA, hệ số tương quan $R^2 \ge 0,965$, hệ số biến thiên $CV < 4,5%$, và giá trị p-value kiểm định sai số thiếu mức khớp ($Lack\ of\ Fit$) $p > 0,05$.

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

               ┌───────────────────────────────────────────────┐
               │    SO SÁNH CÁC CHỈ TIÊU KỸ THUẬT QUAN TRỌNG   │
               └───────────────────────┬───────────────────────┘
                                       │
                 ┌─────────────────────┴─────────────────────┐
                 ▼                                           ▼
      ┌─────────────────────┐                     ┌─────────────────────┐
      │  THỜI GIAN SẤY (h)  │                     │  SEC (kWh/kg H2O)   │
      ├─────────────────────┤                     ├─────────────────────┤
      │ HP-IR:   ██ 5,2h    │                     │ HP-IR:   ██ 3,84    │
      │ HP thuần:████ 10,2h │                     │ HP thuần:█████ 6,85 │
      │ Đối lưu: ██████14,5h│                     │ Đối lưu: ████████12, │
      └─────────────────────┘                     └─────────────────────┘

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

  • Phát hiện 1 (Động thái áp suất nội và dòng ẩm cưỡng bức): Lần đầu tiên ghi nhận và định lượng thành công sự xuất hiện của xung áp suất dư bên trong tôm sấy HP-IR. Trong 60 phút đầu, nhiệt độ tâm tôm tăng nhanh từ $28^\circ\text{C}$ lên $46,5^\circ\text{C}$ do bức xạ hồng ngoại đâm xuyên, làm tăng áp suất mao quản lên cực đại $\Delta P_{max} = 1420 \text{ Pa}$ tại $\tau = 45 \text{ phút}$. Gradien áp suất này tạo ra dòng ẩm đối lưu cưỡng bức chiếm tới $38,5%$ tổng lưu lượng ẩm thoát ra trong giai đoạn đẳng tốc, bác bỏ quan niệm cũ cho rằng áp suất trong sấy tôm luôn bằng 0.
  • Phát hiện 2 (Bất đẳng hướng nhiệt vật lý): Hệ số dẫn nhiệt tôm thẻ biến thiên phi tuyến theo chiều dài bước sấy: ở trạng thái tươi luộc, $\lambda_{//} = 0,542 \text{ W/m.K}$ (dọc thớ cơ) cao hơn $16,8%$ so với $\lambda_\perp = 0,464 \text{ W/m.K}$ (vuông góc thớ cơ). Khi hàm lượng ẩm giảm từ $76,63%$ về $18,0%$, nhiệt dung riêng $C$ giảm từ $3,65 \text{ kJ/kg.K}$ xuống $1,82 \text{ kJ/kg.K}$, trong khi khối lượng riêng thực tế $\rho$ tăng từ $1052 \text{ kg/m}^3$ lên $1245 \text{ kg/m}^3$ do sự co rút mạnh của khung collagen.
  • Phát hiện 3 (Động lực học độ xốp và co rút): Phân tích hình ảnh ImageJ chứng minh độ xốp của tôm giảm nhanh trong pha đầu từ $\varepsilon_{p0} = 0,345$ xuống $\varepsilon_{p} = 0,182$ do cơ chế co rút mao quản bề mặt, sau đó phục hồi nhẹ lên $\varepsilon_{p} = 0,224$ ở giai đoạn cuối do các vi bọt khí giãn nở dưới tác dụng nhiệt IR. Độ co rút thể tích tổng thể đạt $S_v = 44,2%$.
  • Phát hiện 4 (Bộ thông số tối ưu hóa thực nghiệm đa mục tiêu): Giải bài toán tối ưu bề mặt đáp ứng RSM xác định được chế độ vận hành tối ưu chuẩn xác:
    • Nhiệt độ không khí sấy: $t_k = 45,0^\circ\text{C}$
    • Vận tốc dòng tác nhân sấy: $V = 2,50 \text{ m/s}$
    • Khoảng cách giàn đèn hồng ngoại: $H = 41,0 \text{ cm}$
    • Công suất bức xạ hồng ngoại: $IP = 500 \text{ W}$ (tương ứng độ dày lớp sấy $D_{VLS} = 1 \text{ lớp}$)
    • Thời gian sấy thực tế: $\tau_{TN} = 5,20 \text{ giờ}$ (rút ngắn $49,0%$ so với bơm nhiệt đơn thuần $10,2 \text{ giờ}$, và giảm $64,1%$ so với đối lưu gió nóng $14,5 \text{ giờ}$).
    • Suất tiêu hao năng lượng: $SEC_{TN} = 3,84 \text{ kWh/kg } H_2O$ (tiết kiệm $43,9%$ điện năng so với máy sấy bơm nhiệt thông thường).
  • Phát hiện 5 (Bảo toàn chất lượng dinh dưỡng và cảm quan): Sản phẩm tôm khô sấy ở chế độ tối ưu đạt điểm chất lượng cảm quan $CLCQ = 18,65/20$ điểm (màu đỏ cam sáng đồng đều, vỏ bóng, thịt dai chắc). Hàm lượng $NH_3$ chỉ ở mức $26,4 \text{ mg/100g}$ (thấp hơn nhiều ngưỡng giới hạn $45 \text{ mg/100g}$), tỷ lệ hút nước phục hồi đạt $R_{TN} = 2,48 \text{ g vla/g vlk}$, hoạt độ nước đạt ngưỡng an toàn tuyệt đối $a_w = 0,582$, tổng số vi sinh vật hiếu khí $< 10^3 \text{ CFU/g}$, không phát hiện Salmonella và E. coli.

Implications đa chiều

┌────────────────────────────────────────────────────────────────────────┐
│                        HỆ GIÁ TRỊ TÁC ĐỘNG                             │
├──────────────────────────┬──────────────────────┬──────────────────────┤
│ LÝ THUYẾT & PHƯƠNG PHÁP  │ CÔNG NGHIỆP THỰC TẾ  │ CHÍNH SÁCH & XÃ HỘI  │
├──────────────────────────┼──────────────────────┼──────────────────────┤
│ • Bổ sung gradient áp    │ • Tiết kiệm 43,9% SEC│ • Chuỗi giá trị tôm  │
│   suất Darcy vào mô hình │ • Rút ngắn 49% t_sấy │ • Kinh tế carbon thấp│
│ • Cảm biến áp lực kim y tế│ • Giữ protein, vi sinh│ • Tiêu chuẩn xuất khẩu│
└──────────────────────────┴──────────────────────┴──────────────────────┘
  • Học thuật: Cung cấp cơ sở dữ liệu hoàn chỉnh về bộ thông số nhiệt vật lý động học và mở ra phương pháp luận tính toán mô phỏng liên kết đa vật lý (Multiphysics) chính xác cho các đối tượng sinh học có hình dạng phức tạp và cấu trúc bất đẳng hướng.
  • Kỹ thuật công nghệ: Cung cấp giải pháp chế tạo thiết bị sấy lai hóa hiệu suất cao, làm chủ công nghệ thiết kế module điều khiển tự động tỷ lệ công suất nhiệt/hồng ngoại theo thời gian sấy thực tế.
  • Kinh tế - Xã hội: Giúp các doanh nghiệp chế biến thủy sản nâng cao tỷ suất lợi nhuận thông qua giảm chi phí điện năng, loại bỏ phụ thuộc mùa vụ và thời tiết, đồng thời đáp ứng các tiêu chuẩn vệ sinh an toàn thực phẩm khắt khe của các thị trường xuất khẩu cao cấp như EU, Hoa Kỳ, Nhật Bản.

Limitations và Future Research

  1. Ràng buộc về hình học mô hình: Mô hình số 2D và 3D trong COMSOL vẫn xem tiết diện thân tôm có dạng hình học elip đồng nhất trơn nhẵn, chưa mô phỏng chi tiết độ nhám bề mặt vỏ kitin và sự phân tách cục bộ giữa lớp vỏ ngoài với các bó cơ bên trong khi co ngót.
  2. Giới hạn quy mô mẻ sấy: Các số liệu thực nghiệm được thu thập trên lớp vật liệu sấy mỏng (1 lớp tôm xếp khay rải đều). Động thái khí động học và truyền nhiệt truyền chất trong lớp tôm dày xếp chồng nhiều tầng trong các buồng sấy công nghiệp lớn ($> 500 \text{ kg/mẻ}$) chưa được khảo sát đầy đủ.
  3. Bước sóng bức xạ hồng ngoại: Nghiên cứu sử dụng giàn đèn hồng ngoại trung bình phát xạ dải liên tục; chưa thử nghiệm chuyên sâu các loại đèn phát xạ hồng ngoại xung (Pulsed-IR) hoặc đơn sắc có bước sóng điều biến chính xác theo phổ hấp thụ từng phút của tôm.

Chương trình nghiên cứu mở rộng trong tương lai bao gồm:

  • Phát triển thuật toán học máy (Machine Learning / Deep Neural Networks) kết hợp mạng cảm biến IoT để giám sát trực tuyến trường độ ẩm bề mặt và điều khiển thích nghi công suất đèn IR theo thời gian thực (Adaptive Real-time Control).
  • Mở rộng mô hình đa pha tích hợp biến dạng phi tuyến đàn hồi - dẻo (Viscoelastic deformation) mô tả chính xác vết nứt vi mô trong sợi cơ tôm.
  • Nghiên cứu ứng dụng công nghệ sấy HP-IR cho các đối tượng hải sản giá trị cao khác như hải sâm, bào ngư, vi cá và mực một nắng.

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

  • Tác động học thuật: Luận án tạo ra bước tiến khoa học quan trọng trong chuyên ngành Kỹ thuật Nhiệt và Công nghệ Thực phẩm, mở ra hướng nghiên cứu kết hợp mô hình số phần tử hữu hạn với kỹ thuật chẩn đoán áp suất vi mô. Kết quả công trình đã công bố trên các tạp chí khoa học chuyên ngành uy tín trong nước và kỷ yếu hội nghị quốc tế.
  • Chuyển đổi công nghiệp thủy sản: Giải pháp sấy HP-IR mang lại tiềm năng thay thế hoàn toàn các lò sấy thủ công gây ô nhiễm môi trường tại các làng nghề truyền thống dọc duyên hải miền Trung và Đồng bằng sông Cửu Long, nâng tỷ lệ thu hồi tôm loại 1 lên trên $92%$.
  • Hiệu quả năng lượng và môi trường: Với mức tiết kiệm điện $43,9%$ và cắt giảm thời gian sấy phân nửa, công nghệ này đóng góp trực tiếp vào cam kết quốc gia về giảm phát thải khí nhà kính ròng bằng 0 (Net Zero) trong ngành chế biến nông lâm thủy sản.

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

┌────────────────────────────────────────────────────────────────────────┐
│                        MA TRẬN HƯỞNG LỢI TRỰC TIẾP                     │
├──────────────────────────┬─────────────────────────────────────────────┤
│ 1. NCS & Học giả         │ Kế thừa khung toán học Luikov - Fick - Darcy│
│                          │ Bộ tham số nhiệt vật lý tôm thẻ chân trắng  │
├──────────────────────────┼─────────────────────────────────────────────┤
│ 2. Doanh nghiệp R&D      │ Bản vẽ nguyên lý & module tủ sấy HP-IR      │
│                          │ Bộ thông số sấy tối ưu: 45°C, 2.5m/s, 500W  │
├──────────────────────────┼─────────────────────────────────────────────┤
│ 3. Nhà máy Chế biến      │ Tiết kiệm 43,9% điện năng tiêu thụ (SEC)    │
│                          │ Nâng cao CLCQ (18,65/20), hạ NH3 < 35mg/100g│
├──────────────────────────┼─────────────────────────────────────────────┤
│ 4. Nhà hoạch định        │ Cơ sở khoa học ban hành quy chuẩn kỹ thuật  │
│                          │ Thúc đẩy chuỗi giá trị tôm khô xuất khẩu    │
└──────────────────────────┴─────────────────────────────────────────────┘
  • Nghiên cứu sinh và Giảng viên chuyên ngành Nhiệt - Năng lượng - Thực phẩm: Kế thừa khung toán học mô phỏng liên kết đa vật lý và phương pháp đo kiểm động học sấy hiện đại.
  • Kỹ sư R&D thiết kế thiết bị: Tiếp cận bản thiết kế hệ thống sấy lai ghép điều khiển chính xác các thông số nhiệt độ, độ ẩm, trường bức xạ và lưu lượng gió.
  • Doanh nghiệp xuất khẩu thủy sản: Sở hữu quy trình công nghệ chế biến tôm khô đạt tiêu chuẩn kiểm nghiệm vi sinh và cảm quan quốc tế, giảm giá thành sản xuất.
  • Cơ quan quản lý nhà nước: Có cơ sở dữ liệu định lượng phục vụ công tác xây dựng tiêu chuẩn kỹ thuật quốc gia về hiệu quả năng lượng cho máy sấy công nghiệ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?

Đóng góp lý thuyết độc đáo nhất là việc hoàn thiện hệ phương trình vi phân liên kết 3 trường: nhiệt độ ($t$), độ ẩm ($M$) và áp suất ($P$) dựa trên việc mở rộng hệ phương trình Luikov cho vật liệu xốp sinh học. Khác với các nghiên cứu truyền thống chỉ thừa nhận dòng ẩm khuếch tán Fick, luận án đã tích hợp thành công dòng ẩm đối lưu theo định luật Darcy bị thúc đẩy bởi gradient áp suất sinh ra từ sự gia nhiệt sâu của bức xạ hồng ngoại.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của Chalida Niamnuya et al. (2008) và Worachard et al. (2009) vốn bỏ qua áp suất mao quản và coi các thông số nhiệt vật lý là hằng số hoặc chỉ phụ thuộc độ ẩm, luận án của Lê Như Chính đã:

  • Xây dựng mô hình nhiệt vật lý biến thiên liên tục theo nhiệt độ và tỷ lệ 4 thành phần hóa học chính (nước, protein, lipid, tro).
  • Sáng chế phương pháp thực nghiệm dùng kim tiêm y tế nối áp kế chữ U để trực tiếp đo đạc áp suất nội tôm khi sấy.
  • Tích hợp hiện tượng co rút thể tích $S_v(\tau)$ và biến thiên độ xốp qua xử lý ảnh ImageJ vào mô phỏng FEM trên COMSOL Multiphysics.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là sự xuất hiện của đỉnh áp suất dư $\Delta P = 1420 \text{ Pa}$ tại lõi thân tôm chỉ sau 45 phút sấy HP-IR. Chính đỉnh áp suất này tạo ra gradient áp suất cực lớn hướng từ tâm ra bề mặt, giúp đẩy nước tự do trong các mao quản ra ngoài với tốc độ cao mà không làm nhiệt độ bề mặt tôm bị tăng vọt, nhờ luồng không khí khô $45^\circ\text{C}$ của bơm nhiệt liên tục bốc hơi làm mát bề mặt.

4. Luận án có cung cấp giao thức tái lập thí nghiệm (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp giao thức tái lập chi tiết ở mức độ cao: từ kích thước mẫu tôm chuẩn hóa, thời gian luộc sơ bộ, sơ đồ đấu nối cảm biến, thông số kỹ thuật bơm nhiệt R134a, dải bước sóng đèn gốm hồng ngoại, ma trận thí nghiệm Taguchi/RSM cho đến toàn bộ hệ phương trình toán học kèm điều kiện biên vi phân và lưu đồ thuật toán giải trên COMSOL Multiphysics.

5. Lộ trình phát triển nghiên cứu 10 năm tiếp theo được định hình ra sao?

Lộ trình 10 năm tập trung vào 3 trọng tâm:

  1. Công nghiệp hóa thiết bị sấy HP-IR quy mô liên tục (Continuous conveyor HP-IR dryer) công suất nhiều tấn/ngày.
  2. Tích hợp trí tuệ nhân tạo (AI/Computer Vision) để nhận diện màu sắc, độ co rút trực tiếp và điều khiển tối ưu hóa thích nghi công suất bức xạ hồng ngoại theo thời gian thực.
  3. Thương mại hóa quy trình sấy lạnh lai ghép cho toàn bộ chuỗi sản phẩm thủy sản cao cấp xuất khẩu của Việt Nam.

Kết luận

Công trình luận án tiến sĩ của tác giả Lê Như Chính đã tạo lập một dấu ấn khoa học xuất sắc với 6 đóng góp cốt lõi:

  1. Thiết lập hoàn chỉnh mô hình toán học giải tích - số mô tả quá trình truyền nhiệt truyền chất liên kết đa trường ($t - M - P$) trong con tôm thẻ chân trắng sấy HP-IR, giải quyết triệt để bài toán tích hợp nguồn nhiệt hồng ngoại thể tích $q_{v,IR}$ và dòng đối lưu áp suất Darcy.
  2. Xây dựng hệ thống công thức giải tích xác định chính xác các thông số nhiệt vật lý ($\lambda, C, \rho, a$) bất đẳng hướng theo nhiệt độ và thành phần hóa sinh học của nguyên liệu thủy sản Việt Nam.
  3. Phát triển thành công giải pháp kỹ thuật đo áp suất mao quản vi mô trong quá trình sấy bằng cảm biến kim y tế kết hợp áp kế chữ U vi sai.
  4. Mô phỏng thành công động thái biến thiên không - thời gian của trường nhiệt, trường ẩm và trường áp suất trên phần mềm COMSOL Multiphysics kết hợp số hóa vi cấu trúc độ xốp bằng ImageJ.
  5. Xác lập chế độ sấy đa mục tiêu tối ưu trên hệ thống HP-IR: $t_k = 45^\circ\text{C}$, $V = 2,5 \text{ m/s}$, $H = 41 \text{ cm}$, $IP = 500 \text{ W}$, giúp rút ngắn $49%$ thời gian sấy và tiết kiệm $43,9%$ điện năng tiêu thụ riêng ($SEC = 3,84 \text{ kWh/kg } H_2O$).
  6. Chứng minh sự vượt trội về chất lượng sản phẩm tôm khô thành phẩm: đạt điểm cảm quan $18,65/20$, bảo toàn màu sắc tự nhiên, duy trì độ dai cơ thịt, khống chế $NH_3 \le 26,4 \text{ mg/100g}$ và đảm bảo các tiêu chuẩn an toàn vi sinh xuất khẩu khắt khe.

Công trình không chỉ nâng tầm lý thuyết sấy môi trường xốp sinh học mà còn mở ra bước ngoặt hiện đại hóa công nghệ cho ngành chế biến thủy hải sản Việt Nam trên thị trường quốc tế.