Tổng quan về luận án

Ngành thủy sản đóng vai trò là một trong những trụ cột kinh tế biển quan trọng nhất của Việt Nam, quốc gia sở hữu 3.670 km bờ biển và xếp hạng 32 trong số 156 quốc gia có biển. Theo thống kê từ Hiệp hội Chế biến và Xuất khẩu Thủy sản Việt Nam (VASEP), "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" [1]. Tuy nhiên, công nghệ chế biến tôm khô nội địa hiện tại chủ yếu phụ thuộc vào phương thức phơi nắng truyền thống hoặc lò sấy than, củi và điện trở gia nhiệt gián tiếp [2, 3]. Các phương pháp sơ khai này dẫn đến thời gian sấy kéo dài (10–12 giờ), làm biến tính sâu hệ protein cơ thịt, suy giảm sắc tố astaxanthin tự nhiên, thất thoát chất ngấm ra hòa tan và không đáp ứng các tiêu chuẩn nghiêm ngặt về vi sinh và vệ sinh an toàn thực phẩm của các thị trường xuất khẩu cao cấp [2, 3].

+-------------------------------------------------------------------------+
|                        BỐI CẢNH KHOA HỌC & THỰC TIỄN                    |
|  Xuất khẩu tôm 2019: 3,38 tỷ USD | Tôm thẻ chân trắng: 2,36 tỷ USD (70%)|
|  Hạn chế: Phơi nắng / Lò thủ công -> Biến tính protein, mất Astaxanthin |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
|                        KHOẢNG TRỐNG NGHIÊN CỨU                          |
|  - Thiếu mô hình liên kết Truyền nhiệt - Truyền chất - Áp suất (TNTC-P) |
|  - Bỏ qua gradient áp suất mao quản (Luikov, Brooker, Niamnuya)         |
|  - Chưa tích hợp độ co rút phi tuyến tính DeffST & tính dị hướng cơ thịt|
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
|                         GIẢI PHÁP ĐỘT PHÁ CỦA ĐỀ TÀI                    |
|  Mô hình hóa đa vật lý HP-IR: Bơm nhiệt tách ẩm + Bức xạ Hồng ngoại     |
|  Đo áp suất buồng cơ thịt bằng đầu dò vi kim y tế tích hợp COMSOL       |
+-------------------------------------------------------------------------+

Khoảng trống học thuật (research gap) then chốt nằm ở sự thiếu vắng một cơ sở lý thuyết mô hình hóa toàn diện về quá trình truyền nhiệt và truyền chất (TNTC) đa pha trong cấu trúc vật liệu xốp mao dẫn sinh học của tôm thẻ chân trắng (Penaeus vannamei). Phần lớn các công trình trước đây chỉ tiếp cận đơn thuần theo phương trình dẫn nhiệt Fourier một chiều hoặc khuếch tán ẩm đẳng nhiệt Fickian, hoàn toàn triệt tiêu biến số áp suất nội tại bên trong các vi mao quản $(\nabla P = 0)$ và bỏ qua ảnh hưởng của hiện tượng co rút tế bào cơ thịt theo thời gian thực [44, 45, 48, 49]. Nhằm giải quyết triệt để bài toán này, luận án tập trung giải quyết các câu hỏi nghiên cứu và giả thuyết sau:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế tương tác động học giữa dòng nhiệt bức xạ hồng ngoại bước sóng trung $(\lambda = 2,7 \div 25\ \mu\text{m})$ và tác nhân sấy không khí khô nhiệt độ thấp từ bơm nhiệt đối với trường nhiệt độ, trường áp suất và trường ẩm bên trong con tôm thẻ diễn ra như thế nào?

  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để xây dựng và giải hệ phương trình vi phân đạo hàm riêng liên kết phi tuyến mô tả đồng thời động học nhiệt Fourier, khuếch tán Fickian và dòng đối lưu Darcy có xét đến biến dạng co rút $D_{effST}$ và tính dị hướng cấu trúc sợi cơ thịt?

  • Câu hỏi nghiên cứu 3 (RQ3): Thông số vận hành tối ưu đa mục tiêu (nhiệt độ $t$, vận tốc gió $V$, khoảng cách bức xạ $H$, công suất hồng ngoại $IP$) của hệ thống sấy bơm nhiệt kết hợp hồng ngoại (HP-IR) là gì để cực đại hóa tỷ lệ hoàn nguyên phục hồi (RR), bảo toàn hàm lượng protein, giữ màu sắc astaxanthin và tối thiểu hóa suất tiêu hao năng lượng riêng (SEC)?

  • Giả thuyết 1 (H1): Sự kết hợp giữa gia nhiệt bức xạ hồng ngoại đâm xuyên tâm $(\delta = 9\text{ mm})$ và làm mát bề mặt bằng không khí khô đối lưu từ bơm nhiệt sẽ thiết lập các gradient nhiệt độ $(\nabla T)$ và gradient áp suất $(\nabla P)$ cùng chiều từ trong tâm ra ngoài bề mặt, tạo lực đẩy đối lưu mao quản vượt trội giúp giảm trên 50% thời gian sấy so với bơm nhiệt đơn thuần.

  • Giả thuyết 2 (H2): Việc tích hợp hệ số khuếch tán ẩm hiệu dụng có xét đến độ co rút thể tích $D_{effST}$ và vận tốc dòng ẩm đối lưu Darcy $u$ vào mô hình phần tử hữu hạn (FEM) trên phần mềm COMSOL Multiphysics sẽ tăng độ chính xác dự báo động học sấy với sai số $R^2 > 0,95$.

Khung lý thuyết của luận án kế thừa và phát triển từ Lý thuyết truyền nhiệt - truyền chất trong vật liệu xốp mao dẫn của A.V. Luikov [47], định luật chuyển động dòng lưu chất qua môi trường xốp Darcy, định luật Fick mở rộng, định luật bức xạ Planck-Wien-Lambert và lý thuyết nhiệt động học chất khô Pakowski [104]. Phạm vi nghiên cứu thực nghiệm được chuẩn hóa trên mẫu tôm thẻ chân trắng nuôi tại các tỉnh duyên hải miền Trung (Quảng Ngãi, Bình Định, Phú Yên, Khánh Hòa) với kích thước hình học chính xác: chiều dài thân $L = 70,8\text{ mm}$, đường kính đốt đầu $d_1 = 14,5\text{ mm}$, đường kính đốt đuôi $d_2 = 6,6\text{ mm}$, khối lượng cơ thịt ăn được chiếm $60,00%$, độ ẩm ban đầu $M_0 = 76,63%$ (tương đương $2,333\text{ kg/kg vlk}$ hay $41.264,97\text{ mol/m}^3$), hàm lượng protein chiếm $19,27%$ chất tươi ($70 \div 80%$ chất khô) và lipid chiếm $0,92%$ [2].

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về kỹ thuật sấy thực phẩm phân hóa thành ba dòng nghiên cứu chính:

                                  Y VĂN QUỐC TẾ & TRONG NƯỚC
   DÒNG BƠM NHIỆT (HP)                DÒNG HỒNG NGOẠI (IR)                  DÒNG LAI GHÉP (HP-IR)
• Strommen (2002), Taşeri (2018):    • Rosenthal (1992),                   • Aktas (2016):
  Tiết kiệm 60-80% năng lượng          Krishnamurthy (2008):                 Rút ngắn 69% thời gian,
• Dirk (2004), Wang & Chang:           $\lambda_{H_2O} = 2,7-25\ \mu\text{m}$ giảm 43,2% năng lượng
  Nhiệt độ thấp, giữ protein           Cơ chế cộng hưởng dao động          • Yun Deng (2014):
• Võ Duy Mạnh, Phạm Văn Tùy:         • Ratti (2011), Supawan (2008):         Tăng Deff mực ống
  Ứng dụng nông sản Việt Nam           Tăng tốc độ, giữ màu đỏ tôm         • Trần Đại Tiến, Lê Như Chính:
• Hạn chế: Thời gian sấy dài (11h)   • Hạn chế: Nóng cục bộ, cong vênh       Thử nghiệm mực & lá hẹ
  1. Dòng nghiên cứu sấy bơm nhiệt (Heat Pump - HP): Các công trình của Strommen [12], Levent Taşeri [13], Dirk et al. [10], Wang & Chang [11], Claussen et al. [15] chỉ ra rằng sấy bơm nhiệt tiêu thụ năng lượng ít hơn $60% \div 80%$ so với sấy không khí nóng đối lưu tiêu chuẩn nhờ chu trình lạnh tuần hoàn thu hồi ẩn nhiệt ngưng tụ. Tại Việt Nam, Võ Duy Mạnh và Lê Chí Hiệp [4], Phạm Văn Tùy [5] chứng minh hiệu quả duy trì màu sắc và hoạt tính sinh học của cà rốt và kẹo dẻo khi sấy ở $35 \div 40^\circ\text{C}$. Tuy nhiên, hạn chế lớn nhất của sấy HP đối với thủy sản là tốc độ bốc hơi bề mặt bị giới hạn bởi quá trình khuếch tán nội, dẫn đến thời gian sấy kéo dài tới $10 \div 11\text{ giờ}$ [6, 7], gây ức chế năng suất công nghiệp.
  2. Dòng nghiên cứu sấy bức xạ hồng ngoại (Infrared - IR): Nghiên cứu của Rosenthal [19], Krishnamurthy et al. [21], Matsuura [24] làm rõ cơ chế truyền năng lượng điện từ bước sóng $\lambda = 2,7 \div 25\ \mu\text{m}$, tác động trực tiếp vào các liên kết cộng hóa trị của phân tử nước $(\nu_1, \nu_2, \nu_3$ với tỷ lệ cường độ hấp thụ $0,07 : 1,47 : 1,00)$, sinh nhiệt ma sát nội tại tức thời. Ratti [27] và Supawan [25] chứng minh sấy IR rút ngắn thời gian làm khô $2 \div 3\text{ lần}$ và duy trì màu đỏ tôm tốt hơn ở $70^\circ\text{C}$ so với không khí nóng. Dẫu vậy, Lê Văn Hoàng [23] và Nguyễn Thị Bích Thủy [26] nhấn mạnh nguy cơ sấy IR độc lập dễ gây cháy cục bộ bề mặt, nứt nẻ và làm biến tính sợi myosin/actin do mật độ thông lượng nhiệt cao.
  3. Dòng nghiên cứu công nghệ sấy lai ghép bơm nhiệt - hồng ngoại (HP-IR): Mustafa Aktas et al. [29, 30] công bố rằng hệ thống HP-IR cắt giảm $69%$ thời gian sấy và giảm $43,2%$ năng lượng tiêu thụ toàn phần. Yun Deng et al. [31, 32] phát hiện cường độ IR $500\text{ W}$ kết hợp bơm nhiệt $50^\circ\text{C}$ nâng cao hệ số khuếch tán ẩm $D_{eff}$, tăng khả năng phục hồi nước và bảo toàn nồng độ axit amin của mực ống. Xiaoyong Song [34, 35] khẳng định HP-IR ức chế phản ứng sẫm màu phi enzym Maillard.
Chỉ số / Tiêu chí so sánh Sấy phơi nắng thủ công [2, 3] Sấy đối lưu gió nóng (HAD) [44, 45] Sấy bơm nhiệt đơn (HP) [6, 7, 12] Sấy bức xạ hồng ngoại (IR) [25, 27] Sấy bơm nhiệt - hồng ngoại (HP-IR) [Luận án]
Thời gian sấy trung bình $24 \div 36\text{ h}$ (phụ thuộc thời tiết) $6 \div 8\text{ h}$ $(t = 80 \div 120^\circ\text{C})$ $10 \div 11\text{ h}$ $(t = 35 \div 45^\circ\text{C})$ $3 \div 4\text{ h}$ $(t = 70 \div 90^\circ\text{C})$ $3,5 \div 4,5\text{ h}$ $(t = 45^\circ\text{C}, H = 35\text{ cm})$
Mức tiêu hao năng lượng 0 (nhưng rủi ro hư hỏng cao) Rất cao (thất thoát nhiệt khí thải) Thấp (giảm $60 \div 80%$ so với HAD) Trung bình (hiệu suất bức xạ $39%$) Rất thấp (giảm $43,2 \div 59,8%$ năng lượng)
Biến tính Protein & Sắc tố Biến tính nặng, ôi hóa lipid Biến tính nhiệt, mất liên kết sợi cơ Giữ nguyên vẹn nhưng dễ nhiễm khuẩn chậm Co cứng bề mặt, biến dạng nhiệt cao Bảo toàn tối đa protein, astaxanthin đỏ tươi
Cơ chế truyền ẩm chủ đạo Khuếch tán Fickian chậm chạp Khuếch tán Fickian do nhiệt độ cao Khuếch tán Fickian do chênh lệch ẩm Bốc hơi nhanh kèm tăng nhiệt cục bộ Khuếch tán Fickian + Đối lưu Darcy ($\nabla T, \nabla P$ cùng chiều)
Mô hình hóa toán học Không áp dụng Fick 2D/3D (bỏ qua $\nabla P$ và co rút) Fick 1D/2D thuần túy Dẫn nhiệt Fourier có nguồn nhiệt thể tích Hệ phương trình đa vật lý FEM (Fourier - Fick - Darcy - Shrinkage)

Về mặt lý thuyết truyền khối, tồn tại một cuộc tranh luận học thuật sâu sắc: Hệ phương trình kinh điển của Luikov [47] đã xác lập mối quan hệ tương hỗ giữa ba trường thế: thế ẩm $M$, thế nhiệt $T$ và thế áp suất $P$:

$$\begin{cases} \frac{\partial M}{\partial \tau} = K_{11} \nabla^2 M + K_{12} \nabla^2 T + K_{13} \nabla^2 P \ \frac{\partial T}{\partial \tau} = K_{21} \nabla^2 M + K_{22} \nabla^2 T + K_{23} \nabla^2 P \ \frac{\partial P}{\partial \tau} = K_{31} \nabla^2 M + K_{32} \nabla^2 T + K_{33} \nabla^2 P \end{cases}$$

Tuy nhiên, các công trình quốc tế tiêu biểu của Chalida Niamnuya et al. [44] (Thái Lan) khi mô phỏng sấy tầng sôi tôm thẻ trên COMSOL 3.3, Worachard Chawanasporn et al. [45] sử dụng phần mềm FEMLAB, cũng như Brooker [48] và Trần Văn Phú [49] đều quy ước đơn giản hóa mô hình bằng cách triệt tiêu toàn bộ ảnh hưởng của gradient áp suất $(K_{13} = K_{23} = K_{33} = 0)$ vì cho rằng áp suất trong lòng vật liệu không đáng kể.

Luận án này định vị sự đột phá bằng việc chứng minh luận điểm của Mingheng Shi et al. [85] và Soner Celen et al. [86]: "Cường độ của dòng dịch chuyển ẩm đối lưu lớn hơn nhiều lần dòng dịch chuyển ẩm dưới tác động của lực mao dẫn và khuếch tán". Khi có nguồn gia nhiệt bức xạ thể tích đâm xuyên, nước tự do trong mao quản lập tức chuyển pha sinh hơi, gia tăng áp suất bão hòa cục bộ $\Delta P$, tạo nên động lực dòng chảy Darcy mãnh liệt từ trong tâm ra bề mặt. Luận án đã lấp đầy khoảng trống tồn tại nhiều thập kỷ bằng cách giải quyết tường minh bài toán áp suất mao quản trong cấu trúc thịt tôm thẻ.

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

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

Luận án mở rộng thuyết truyền nhiệt truyền chất kinh điển thông qua việc xây dựng hệ phương trình vi phân liên kết đa vật lý 3 chiều không ổn định, dung hợp ba quy luật cơ bản:

                            HỆ THỐNG PHƯƠNG TRÌNH LIÊN KẾT ĐA VẬT LÝ
  1. Phương trình vi phân phân bố nhiệt độ (Định luật Fourier mở rộng):

$$\rho C_{eq} \frac{\partial T}{\partial \tau} = \nabla \cdot (\lambda_{TB} \nabla T) + q_{v,IR}$$

Trong đó, dung lượng nhiệt tương đương và hệ số dẫn nhiệt hiệu dụng được tích hợp từ các pha vật chất: $C_{eq} = \theta_p C_p + (1 - \theta_p) \rho_d C_d$; $\lambda_{TB} = \frac{\lambda_{//} + \lambda_{\bot}}{2}$ phản ánh tính chất dị hướng của sợi cơ thịt tôm; $q_{v,IR} = \frac{Q_{IR}}{V_{shr}}\ [\text{W/m}^3]$ là mật độ nguồn nhiệt thể tích hấp thụ từ bức xạ hồng ngoại có xét đến độ co rút thể tích tức thời $V_{shr}$.

  1. Phương trình vi phân phân bố hàm lượng ẩm (Định luật Fick II kết hợp Định luật Darcy):

$$\frac{\partial M}{\partial \tau} = \nabla \cdot (D_{effST} \nabla M) - \nabla \cdot (u M)$$

Dòng ẩm tổng cộng bên trong cấu trúc mao quản tôm là sự chồng chập của dòng khuếch tán phân tử $J_{diff} = -D_{effST} \nabla M$ và dòng đối lưu áp suất $J_{conv} = u M$.

  1. Phương trình dòng ẩm đối lưu theo gradient áp suất (Định luật Darcy):

$$u = -\frac{K}{\mu_{eff}} \nabla P = -\left( \frac{K K_{r,l}}{\mu_l} \nabla P_l + \frac{K K_{r,v}}{\mu_v} \nabla P_v \right)$$

Định luật này thể hiện "sự chênh lệch nhiệt độ giữa tâm và bề mặt ngoài của VLS lớn sẽ xuất hiện dòng ẩm đối lưu di chuyển trong các mao quản của VLS do chênh lệch áp suất theo định luật Darcy và dòng ẩm khuếch tán di chuyển do chênh lệch nhiệt độ theo định luật Fick", tạo nên bước chuyển biến mô hình (paradigm shift) từ mô hình khuếch tán thụ động sang mô hình động lực học đa pha chủ động.

                          CƠ CHẾ ĐỘNG LỰC HỌC TRUYỀN ẨM HP-IR
        
            BỨC XẠ HỒNG NGOẠI (IR)             KHÔNG KHÍ KHÔ BƠM NHIỆT (HP)
             [Bước sóng 2,7-25 µm]               [Nhiệt độ thấp 35-45°C]
             Đâm xuyên tâm (δ=9mm)               Làm mát & Tách ẩm bề mặt
              Nhiệt độ tâm T_core                 Nhiệt độ bề mặt T_surf
             (T_core > T_surf)                   (Ẩm bề mặt bốc hơi nhanh)

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

Khung phân tích của luận án tích hợp bốn lý thuyết cốt lõi trong một hệ quy chiếu đồng nhất:

  • Lý thuyết bức xạ nhiệt Lambert - Wien: Xác định chiều dài bước sóng phát xạ cực đại của đèn hồng ngoại tại nhiệt độ bề mặt $t_{IR} = 160^\circ\text{C}$ ($T_{IR} = 433\text{ K}$) theo định luật Wien: $\lambda_{\max} = \frac{2898}{T_{IR}} = 6,69\ \mu\text{m}$. Bức xạ này nằm chuẩn xác trong vùng hấp thụ cực đại của nước sinh học $(\lambda = 2,7 \div 25\ \mu\text{m})$. Độ sâu thâm nhập bức xạ đạt $\delta = 0,009\text{ m}$ ($9\text{ mm}$), vượt qua bán kính trung bình thân tôm ($R_{TB} = 7,25\text{ mm}$), đảm bảo toàn bộ thể tích tôm tiếp nhận năng lượng đồng đều theo định luật hấp thụ Lambert: $Q_{IR} = q_0 \exp(-\alpha \delta)$.
  • Lý thuyết nhiệt vật lý thành phần Ashrae - Paul Singh: Xác định chính xác các đại lượng nhiệt vật lý thông qua phân tích thành phần hóa học cơ bản gồm nước $(w_w)$, protein $(w_p)$, lipid $(w_l)$, tro $(w_a)$, carbohydrate $(w_c)$:
    • Dẫn nhiệt song song: $\lambda_{//} = \sum \lambda_i \theta_{pi}$
    • Dẫn nhiệt vuông góc: $\frac{1}{\lambda_{\bot}} = \sum \frac{\theta_{pi}}{\lambda_i}$
    • Nhiệt dung riêng: $C = \sum C_i w_i$
    • Khối lượng riêng: $\rho = \frac{1}{\sum (w_i / \rho_i)}$
  • Lý thuyết động học co rút và biến tính mao quản: Hiệu chỉnh hệ số khuếch tán ẩm hiệu dụng phụ thuộc đồng thời vào nhiệt độ cục bộ và sự co ngót diện tích bề mặt $A_p(\tau)$ theo mô hình Chandan Kumar [117]:

$$D_{effST} = \frac{D_0 \exp\left(-\frac{E_a}{R T}\right) + D_{effF} \left(\frac{R_{TB,\tau}}{R_{TB,0}}\right)^2}{2}$$

Điều kiện biên (boundary conditions) được thiết lập chặt chẽ cho ba bề mặt hình học của thân tôm:

  • Biên truyền nhiệt: $-k \nabla T \cdot \vec{n} = q_0 + \alpha_{dl}(T_s - T_a) + \varepsilon \sigma (T_s^4 - T_{surr}^4) + \alpha_m \rho_w r (M_s - M_e)$
  • Biên truyền chất: $-D_{effST} \nabla M \cdot \vec{n} = \alpha_m (M_s - M_e)$

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý nghiên cứu thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), đi từ bản chất cấu trúc vi mô đến lượng hóa thực nghiệm vĩ mô. Thiết kế nghiên cứu đa tầng (Multi-level Multiphysics Design) bao gồm:

  1. Tầng mô hình hóa giải tích và số trị: Lập trình giải hệ phương trình vi phân phi tuyến liên kết trên bộ giải phần tử hữu hạn COMSOL Multiphysics 4.x.
  2. Tầng thực nghiệm định chuẩn: Chế tạo hệ thống máy sấy bơm nhiệt kết hợp hồng ngoại quy mô pilot và cảm biến đo đạc áp suất nội vi mô.
  3. Tầng tối ưu hóa đa mục tiêu: Áp dụng phương pháp bề mặt đáp ứng (Response Surface Methodology - RSM) dựa trên quy hoạch thực nghiệm trực giao tâm xoay (CCRD).
                      THIẾT KẾ NGHIÊN CỨU ĐA TẦNG TOÀN DIỆN

Quy trình nghiên cứu rigorous

  • Chuẩn bị mẫu sinh học: Tôm thẻ chân trắng tươi sống (Penaeus vannamei) thu gom tại đầm nuôi đạt kích cỡ thương phẩm đồng nhất $60 \div 70\text{ con/kg}$. Quy trình xử lý sơ chế nghiêm ngặt: rửa bằng nước đá lạnh $4^\circ\text{C}$, vặt đầu, bóc vỏ, rút tim gan và ruột nội tạng, chỉ giữ lại cơ thịt nguyên vẹn.
  • Xác định thông số hình học: Sử dụng thước kẹp điện tử Mitutoyo (độ chính xác $\pm 0,01\text{ mm}$) đo đạc đa điểm trên 100 mẫu thử, số hóa hình học đốt thân 2D và 3D chuyển nạp vào môi trường CAD/COMSOL.
  • Thiết bị đo áp suất mao quản sáng tạo: Chế tạo đầu dò đo áp suất bên trong lòng cơ thịt tôm bằng cách cải tiến kim tiêm truyền dịch y tế vô trùng (đường kính ngoài $0,45\text{ mm}$) gắn kết nối trực tiếp với cảm biến áp suất vi sai điện tử piezo có độ nhạy cao. Đầu kim được cố định chính xác vào tâm đốt thân thứ 3 của tôm, truyền tín hiệu áp suất $\Delta P$ theo mili-giây về bộ thu thập dữ liệu (Data Acquisition System).
  • Hệ thống sấy HP-IR: Hệ thống kín tuần hoàn không khí, dàn bay hơi tách ẩm, máy nén lạnh R134a, dàn ngưng nhiệt sưởi ấm, kết hợp cụm bóng đèn hồng ngoại công suất bức xạ điều biến $(300 \div 1000\text{ W})$, biến tần điều khiển quạt gió vận tốc $0,5 \div 4,0\text{ m/s}$.

Data và phân tích

                        CÁC CÔNG THỨC NHIỆT VẬT LÝ ĐẦU VÀO ĐƯỢC LẬP TRÌNH TRÊN COMSOL

Dữ liệu mô phỏng được rời rạc hóa bằng lưới tam giác phi tuyến (Free Triangular Mesh) với hơn 45.000 phần tử biên và phần tử khối. Thuật toán bước thời gian ngầm định (Implicit Time-stepping Solver - BDF) được áp dụng với bước nhảy $\Delta \tau = 1\text{ s}$. Kiểm định độ tương thích giữa mô phỏng và thực nghiệm thông qua các chỉ số thống kê chuẩn mực: Hệ số xác định $R^2$, Sai số bình phương trung bình gốc (RMSE) và Sai số phần trăm tuyệt đối trung bình (MAPE):

$$\text{RMSE} = \sqrt{\frac{1}{N} \sum_{i=1}^N (M_{exp,i} - M_{pre,i})^2} \le 0,023\text{ kg/kg vlk}; \quad R^2 \ge 0,985$$

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

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

                               CÁC PHÁT HIỆN THEN CHỐT TỪ THỰC NGHIỆM
  1. Sự cộng hưởng đâm xuyên của bước sóng hồng ngoại tối ưu: Tại nhiệt độ bóng phát $160^\circ\text{C}$, bước sóng đỉnh $\lambda_{\max} = 6,69\ \mu\text{m}$ xuyên sâu $9\text{ mm}$, bao phủ $100%$ thể tích hình học thân tôm. Năng lượng bức xạ kích hoạt dao động kéo giãn và uốn khúc của liên kết phân tử nước $\text{H-O-H}$ ở tần số $\nu_2$, chuyển hóa tức thời thành nội nhiệt tại tâm mà không phụ thuộc vào độ dẫn nhiệt kém của lớp chất béo bề mặt.
  2. Xác lập sự tồn tại của xung áp suất mao quản $(\Delta P)$: Thực nghiệm bằng đầu dò vi kim y tế ghi nhận áp suất trong lõi tôm tăng vọt lên $1,35 \div 1,68\text{ bar}$ trong $30 \div 45\text{ phút}$ đầu của giai đoạn sấy đẳng tốc, sau đó giảm tiệm cận về áp suất khí quyển khi hàm lượng ẩm hạ xuống dưới mức tới hạn $M_k = 0,8\text{ kg/kg vlk}$. Dòng Darcy đối lưu do gradient áp suất đóng góp tới $42%$ tổng lưu lượng ẩm thoát ra ở nửa đầu chu trình sấy.
  3. Hiệu ứng hiệp đồng gia tốc thoát ẩm: Nhờ $\nabla T$, $\nabla P$ và $\nabla M$ cùng hướng từ tâm ra bề mặt, thời gian sấy tôm thẻ giảm từ $11\text{ giờ}$ (khi sấy bơm nhiệt thuần túy) xuống còn $3,8\text{ giờ}$ trên hệ thống HP-IR, tiết kiệm $52,1%$ thời gian và giảm $59,8%$ năng lượng điện tiêu thụ riêng.
  4. Bảo tồn cấu trúc sinh hóa và phẩm chất cảm quan vượt trội: Do bề mặt tôm liên tục được làm mát bởi tác nhân sấy nhiệt độ thấp $(40 \div 45^\circ\text{C})$ từ bơm nhiệt, nhiệt độ tối đa của cơ thịt không vượt quá ngưỡng biến tính nhiệt cục bộ $(55^\circ\text{C})$. Phân tích sinh hóa cho thấy sản phẩm giữ được $94,2%$ hàm lượng astaxanthin, tỷ lệ phục hồi hút nước (RR) đạt $88,6%$ (cao hơn $24,5%$ so với sấy gió nóng đối lưu) và tổ chức tế bào cơ thịt có độ xốp mao quản $\varepsilon_p = 0,42$ đồng nhất.

Implications đa chiều

  • Về mặt học thuật: Mô hình toán Fourier - Fick - Darcy đa vật lý kèm bộ thông số nhiệt vật lý giải tích cung cấp công cụ mô phỏng chuẩn xác cho ngành kỹ thuật nhiệt và cơ học lưu chất trong vật liệu xốp sinh học.
  • Về mặt công nghệ chế tạo: Cung cấp thông số kỹ thuật cốt lõi (công suất nhiệt, bề mặt bức xạ, lưu lượng gió, tỷ lệ bypass không khí) để các kỹ sư thiết kế chế tạo buồng sấy HP-IR quy mô công nghiệp tự động hóa.
  • Về mặt chuỗi giá trị thủy sản: Cho phép các nhà máy chế biến thủy sản tại Đồng bằng sông Cửu Long và Duyên hải miền Trung chuyển đổi từ xuất khẩu nguyên liệu thô đông lạnh sang chế biến sâu tôm khô cao cấp xuất khẩu đạt tiêu chuẩn HACCP và ISO 22000.

Limitations và Future Research

Mặc dù đạt được những đóng góp đột phá, luận án thẳng thắn thừa nhận 3 hạn chế học thuật:

  1. Giả thiết hình học bán cứng: Quá trình mô phỏng sử dụng hệ số khuếch tán ẩm có xét co rút $D_{effST}$ trên lưới phần tử cố định, chưa áp dụng phương pháp lưới chuyển động biến dạng động lực học phi tuyến (Arbitrary Lagrangian-Eulerian - ALE) để quan sát sự méo mó hình học 3D thực tế của các đốt tôm theo thời gian thực.
  2. Quy mô lớp vật liệu mỏng đơn tầng: Các tính toán thực nghiệm giới hạn ở chế độ sấy một lớp (single layer) trên lưới đỡ nhằm bảo đảm tính đồng nhất bức xạ, chưa giải quyết trọn vẹn hiện tượng che khuất quang học và cản trở khí động học khi sấy lớp dày (deep-bed drying) hoặc sấy tầng sôi liên tục.
  3. Giới hạn biến thiên thành phần sinh học: Các công thức nhiệt vật lý được tính toán trên giá trị trung bình mẫu tôm thương phẩm khỏe mạnh, chưa bao quát sự biến động tỷ lệ mỡ/protein theo mùa vụ sinh sản hoặc độ mặn của nguồn nước nuôi.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Phát triển mô hình mô phỏng biến dạng cơ học đàn hồi - nhớt - dẻo 3D (Viscoelastic Shrinkage FEM) kết hợp giải thuật ALE trên COMSOL Multiphysics.
  • Nghiên cứu cơ chế điều khiển tối ưu thời gian thực (Real-time Adaptive Control) sử dụng thuật toán trí tuệ nhân tạo (AI) điều biến bước sóng đèn hồng ngoại LED (NIR/MIR) thích ứng theo sự biến đổi độ ẩm bề mặt.
  • Mở rộng ứng dụng công nghệ HP-IR đa năng sang các đối tượng thủy hải sản cao cấp khác như hải sâm, bào ngư, mực ống và cá ngừ đại dương.

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

  • Tác động học thuật: Dự báo công trình sẽ thu hút các trích dẫn trong các tạp chí Q1 chuyên ngành kỹ thuật nhiệt và kỹ thuật thực phẩm quốc tế như Journal of Food Engineering, Drying Technology, và Applied Thermal Engineering.
  • Tác động kinh tế ngành: Công nghệ HP-IR rút ngắn thời gian sấy và tối ưu hóa năng suất giúp giảm giá thành chế biến tôm khô xuất khẩu từ $18 \div 22%$, nâng cao biên lợi nhuận của doanh nghiệp chế biến thủy sản.
  • Tác động môi trường và xã hội: Việc loại bỏ hoàn toàn lò đốt than đá, củi giúp cắt giảm $100%$ khí thải độc hại $\text{CO}_2, \text{SO}_x, \text{NO}_x$ và bụi mịn tại các làng nghề chế biến hải sản truyền thống, đồng thời đảm bảo tạo ra sản phẩm sạch, bảo vệ sức khỏe người tiêu dùng toàn cầu.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận hệ phương trình đa vật lý toàn diện và bộ số liệu nhiệt vật lý thực nghiệm tôm thẻ chân trắng Việt Nam được chuẩn hóa.
  • Kỹ sư R&D và Nhà sản xuất thiết bị: Sở hữu công thức thiết kế chính xác, giải thuật mô phỏng và bộ thông số tối ưu để thương mại hóa máy sấy HP-IR thế hệ mới.
  • Doanh nghiệp và Hợp tác xã chế biến thủy sản: Nâng cao chất lượng sản phẩm tôm khô đạt tiêu chuẩn xuất khẩu sang các thị trường khó tính (EU, Hoa Kỳ, Nhật Bản), gia tăng giá trị chuỗi cung ứng tôm Việt Nam.
  • Nhà hoạch định chính sách: Căn cứ khoa học chuẩn xác phục vụ việc xây dựng chiến lược phát triển kinh tế biển xanh và chính sách hỗ trợ chuyển đổi năng lượng trong công nghiệp chế biến sau thu hoạch.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo và mang tính nền tảng nhất của luận án là gì?

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 của A.V. Luikov bằng việc thiết lập phương trình liên kết đa vật lý 3 chiều tích hợp đồng thời trường nhiệt Fourier (có nguồn nhiệt bức xạ thể tích $q_{v,IR}$), trường ẩm Fickian (có hệ số co rút $D_{effST}$) và đặc biệt là trường dòng ẩm đối lưu Darcy chịu sự chi phối của gradient áp suất mao quản sinh ra do bốc hơi nội vi thể $(\nabla P)$. Đây là công trình đầu tiên định lượng và giải thành công biến số áp suất nội tôm thẻ trong công nghệ sấy HP-IR.

2. Sự đổi mới về phương pháp luận so sánh với các công trình quốc tế tiền nhiệm?

So với nghiên cứu của Chalida Niamnuya et al. [44] (COMSOL 3.3 sấy tầng sôi) và Worachard Chawanasporn et al. [45] (FEMLAB sấy khí nóng) vốn áp dụng phương pháp khuếch tán đẳng áp $(\nabla P = 0)$ trên hình học bán cầu đơn giản hóa, luận án này:

  • Ứng dụng giải phẫu học hình học 3 chiều thực tế của thân tôm theo từng đốt phân đoạn.
  • Sáng tạo kỹ thuật đo áp suất nội vi mô bằng đầu dò vi kim y tế kết nối cảm biến áp điện.
  • Tích hợp tính chất nhiệt vật lý dị hướng $(\lambda_{//}, \lambda_{\bot})$ theo phương thớ cơ thịt.

3. Phát hiện nào bất ngờ nhất đi ngược lại nhận thức kỹ thuật truyền thống?

Trong kỹ thuật sấy thông thường, gia nhiệt nhiệt độ cao là điều kiện bắt buộc để tăng tốc độ thoát ẩm. Tuy nhiên, luận án chứng minh rằng trong công nghệ HP-IR, việc hạ thấp nhiệt độ tác nhân sấy bơm nhiệt xuống $40 \div 45^\circ\text{C}$ kết hợp chiếu xạ hồng ngoại kiểm soát lại tạo ra tốc độ bốc hơi nước nhanh hơn và hiệu quả năng lượng cao hơn sấy ở $70 \div 80^\circ\text{C}$. Điều này là do nhiệt độ bề mặt thấp giúp duy trì độ mở của các vi mao quản, ngăn ngừa hoàn toàn hiện tượng tạo màng cứng (case hardening) gây tắc nghẽn đường thoát ẩm.

4. Luận án có cung cấp đầy đủ quy trình và thông số để tái lập thực nghiệm không?

Luận án cung cấp bộ quy trình chuẩn hóa $100%$ có khả năng tái lập:

  • Phương trình toán học, điều kiện biên, điều kiện ban đầu, và bảng tham số nhiệt vật lý giải tích hóa theo nhiệt độ và thành phần hóa học $(w_w, w_p, w_l, w_a, w_c)$.
  • Kích thước hình học mẫu tôm tiêu chuẩn và thông số cài đặt bộ giải trên COMSOL Multiphysics.
  • Thông số chế độ vận hành tối ưu đa mục tiêu: Nhiệt độ sấy $t = 45^\circ\text{C}$, Vận tốc gió $V = 2,5\text{ m/s}$, Khoảng cách bức xạ hồng ngoại $H = 35 \div 41\text{ cm}$, Độ ẩm không khí cấp $\varphi = 20 \div 40%$.

5. Định hướng lộ trình nghiên cứu 10 năm tiếp theo (2026–2035) từ nền tảng luận án?

  • Giai đoạn 2026–2028: Phát triển mô hình biến dạng phi tuyến ALE 3D và thử nghiệm sấy HP-IR liên tục dạng băng tải rung cho cụm nhà máy công nghiệp.
  • Giai đoạn 2029–2032: Tích hợp nguồn bức xạ hồng ngoại LED đa bước sóng (Multi-spectral IR) điều khiển thông minh bằng thuật toán học sâu (Deep Learning) nhận diện ẩm trực tiếp qua ảnh nhiệt phổ (Hyperspectral Imaging).
  • Giai đoạn 2033–2035: Tự động hóa toàn diện quy trình chế biến không phát thải (Zero-emission Post-harvest Drying) tích hợp hệ thống điện mặt trời áp mái và bơm nhiệt môi chất lạnh tự nhiên $\text{CO}_2$ (R744).

Kết luận

  1. Xây dựng thành công mô hình toán học đa vật lý hoàn chỉnh: Hợp nhất chuẩn xác các định luật Fourier, Fick và Darcy, giải quyết triệt để bài toán liên kết truyền nhiệt, truyền chất và động lực học áp suất mao quản có xét đến hiện tượng co rút thể tích trong tôm thẻ chân trắng sấy HP-IR.
  2. Đột phá phương pháp đo đạc thực nghiệm: Ứng dụng thành công kỹ thuật cảm biến đầu dò vi kim y tế để đo đạc trực tiếp áp suất nội vi mô trong lòng vật liệu sinh học, cung cấp bằng chứng thực nghiệm vững chắc chứng minh sự tồn tại của dòng ẩm đối lưu Darcy trong pha sấy đẳng tốc.
  3. Chuẩn hóa bộ cơ sở dữ liệu nhiệt vật lý chuyên sâu: Thiết lập hệ thống công thức giải tích xác định chính xác hệ số dẫn nhiệt dị hướng, nhiệt dung riêng, khối lượng riêng và hệ số khuếch tán ẩm theo nhiệt độ và thành phần sinh hóa của tôm thẻ chân trắng Việt Nam.
  4. Xác lập chế độ công nghệ sấy tối ưu mang tính thực tiễn cao: Chứng minh chế độ sấy HP-IR tối ưu đa mục tiêu ($t = 45^\circ\text{C}, V = 2,5\text{ m/s}, H = 35 \div 41\text{ cm}$) giúp rút ngắn thời gian sấy xuống $3,5 \div 4,5\text{ giờ}$, tiết kiệm hơn $50%$ năng lượng tiêu thụ, nâng cao tỷ lệ hoàn nguyên lên $88,6%$ và bảo toàn phẩm chất màu sắc astaxanthin đỏ tươi đặc trưng.
  5. Mở ra hướng nghiên cứu và ứng dụng công nghiệp sâu rộng: Đặt nền móng vững chắc cho việc chuyển giao công nghệ, thiết kế chế tạo hệ thống máy sấy HP-IR quy mô công nghiệp, góp phần hiện đại hóa ngành chế biến thủy sản sau thu hoạch và nâng cao năng lực cạnh tranh quốc gia của thương hiệu tôm Việt Nam trên trường quốc tế.