Tổng quan về luận án

Ngành công nghiệp chế biến lâm sản Việt Nam duy trì vị thế mũi nhọn trong nền kinh tế với giá trị xuất khẩu gỗ và lâm sản đạt 15,87 tỉ USD, xuất siêu đạt 12,94 tỉ USD, đóng góp trên 30% tổng kim ngạch xuất khẩu nông lâm thủy sản toàn quốc. Trong bối cảnh nguồn cung nguyên liệu gỗ cứng nhiệt đới có giá trị kinh tế cao như Căm xe (Xylia xylocarpa) phụ thuộc lớn vào nhập khẩu với nhu cầu gia tăng nhanh chóng, bài toán công nghệ đặt ra là phải tối ưu hóa quy trình gia công, trong đó sấy gỗ là khâu then chốt quyết định chất lượng thành phẩm nội ngoại thất. Gỗ Căm xe thuộc nhóm gỗ rất nặng, kết cấu dị hướng phức tạp, vách tế bào dày, đặc biệt chứa hàm lượng cao chất chiết xuất, gôm nhựa và tinh thể phân bố dày đặc trong hệ thống ống mạch ($7,5 - 12,5\text{ lỗ/mm}^2$) và tia gỗ (chiếm 5 - 10% thể tích thân cây). Khi áp dụng các phương pháp sấy đối lưu cưỡng bức truyền thống ở nhiệt độ cao, gôm nhựa bị hóa lỏng, dịch chuyển và kết tủa làm bít tắc các mao dẫn thoát ẩm, làm gia tăng cục bộ gradient ẩm độ và ứng suất sấy, dẫn đến tỷ lệ khuyết tật nứt tét, co rút và cong vênh rất cao.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ sự thiếu vắng các dữ liệu nhiệt vật lý thực nghiệm có tính hệ thống của gỗ Căm xe theo hàm độ ẩm, đồng thời chưa có mô hình toán học giải tích số mô tả chính xác quá trình truyền nhiệt và truyền ẩm đa chiều trong môi trường chân không kết hợp bức xạ nhiệt hồng ngoại. Nghiên cứu của Bùi Thị Thiên Kim (2023) tại Trường Đại học Nông Lâm TP. Hồ Chí Minh đã giải quyết triệt để vấn đề này thông qua hai câu hỏi nghiên cứu và hai giả thuyết khoa học tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Các thông số nhiệt vật lý cơ bản (khối lượng riêng, nhiệt dung riêng, hệ số dẫn nhiệt dọc và ngang thớ, điểm bão hòa thớ gỗ, độ ẩm cân bằng) của gỗ Căm xe biến thiên như thế nào theo hàm độ ẩm $W \in [10%, 40%]$?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật phân bố động học của trường nhiệt và trường ẩm bên trong thanh gỗ Căm xe khi sấy chân không bức xạ hồng ngoại diễn ra như thế nào và đâu là chế độ công nghệ tối ưu để tối thiểu hóa thời gian sấy lẫn tỷ lệ khuyết tật?
  • Giả thuyết nghiên cứu 1 (H1): Quá trình truyền nhiệt và chuyển khối trong gỗ Căm xe dưới điều kiện chân không bức xạ hồng ngoại tuân theo hệ phương trình đạo hàm riêng liên kết phi tuyến tính mô tả dòng dẫn nhiệt Fourier, bức xạ bề mặt Stefan-Boltzmann và cơ chế khuếch tán ẩm do gradient áp suất và gradient độ ẩm.
  • Giả thuyết nghiên cứu 2 (H2): Hạ áp suất môi trường xuống mức áp suất chân không thấp kết hợp gia nhiệt thể tích bằng bức xạ hồng ngoại dải sóng ngắn $\lambda_{max} = 2,5 - 3,5\text{ }\mu\text{m}$ sẽ hạ điểm sôi của nước tự do xuống dưới nhiệt độ nóng chảy của gôm nhựa, tạo động lực thoát ẩm mãnh liệt mà không làm phá vỡ cấu trúc cơ lý của gỗ.

Khung lý thuyết của luận án được xây dựng dựa trên lý thuyết nhiệt động lực học môi trường xốp dị hướng của Luikov, định luật dẫn nhiệt Fourier, định luật khuếch tán Fick, định luật bức xạ Planck và giải pháp giải số sai phân - phần tử hữu hạn (Finite Element Method - FEM). Đóng góp đột phá của công trình được lượng hóa cụ thể: xây dựng thành công bộ 4 phương trình tương quan nhiệt vật lý với hệ số xác định $R^2 \ge 0,9532$; mô phỏng số học trường nhiệt - ẩm trên phần mềm COMSOL Multiphysics với sai số trung bình so với thực nghiệm dưới 5%; và xác lập bộ thông số công nghệ sấy chân không tối ưu đa mục tiêu ($T_s = 58,5 - 58,9^\circ\text{C}$; $p = 0,2\text{ bar}$; $P_{hn} = 625 - 641\text{ W/m}^2$) trên 4 quy cách kích thước phôi gỗ từ 20 mm đến 50 mm chiều dày. Phạm vi nghiên cứu được thực hiện trên nguyên liệu gỗ Căm xe tự nhiên tròn có đường kính $0,5 - 0,7\text{ m}$, chiều dài $3 - 7\text{ m}$, độ ẩm ban đầu đạt giá trị trung bình $W_{woIN} = 40,36%$.

Literature Review và Positioning

Lịch sử phát triển công nghệ sấy chân không vật liệu gỗ ghi nhận các cột mốc quan trọng từ sáng chế buồng sấy chân không tiếp xúc đầu tiên của Vincenzo Pagnozzi (1965), quy trình sấy trong dầu sôi chân không Royal của Bror Hager (1971), đến các công trình sấy nén ép chân không của Haygreen (1986) và Ressel (1994, 1999). Bước sang thế kỷ 21, các nghiên cứu quốc tế chuyển dịch mạnh mẽ sang các hệ thống sấy lai (hybrid drying) kết hợp chân không với các nguồn năng lượng trường như sóng cao tần, vi sóng và bức xạ điện từ nhằm khắc phục nhược điểm truyền nhiệt kém trong môi trường áp suất thấp.

Trong bức tranh tổng quan học thuật, tồn tại cuộc tranh luận sâu sắc giữa hai trường phái lý thuyết về động lực thoát ẩm trong môi trường chân không:

  • Trường phái 1 (Dòng mao dẫn dọc thớ chiếm ưu thế): Zhangjing Chen (1997) khi khảo sát động học sấy chân không ở nhiệt độ $60^\circ\text{C}$ và áp suất $18\text{ mmHg}$ đã lập luận rằng tốc độ sấy không phụ thuộc đáng kể vào chiều dày ván mà bị chi phối áp đảo bởi chiều dài thanh gỗ, do các mao quản mạch gỗ dọc thớ tạo thành đường thoát ẩm chính khi áp suất hơi nước nội tại gia tăng.
  • Trường phái 2 (Gradient ẩm đa chiều và trở lực cấu trúc ngang thớ): Ngược lại, Lin Zhang, Stavros Avramidis và Savvas G. Hatzikiriakos (1997), kết hợp với các nghiên cứu mô hình hóa thế nước 2D của Maurice Defo, Alain Cloutier và Yves Fortin (2000) trên gỗ Thích (Acer saccharum), chứng minh rằng sự thoát ẩm theo phương xuyên tâm và tiếp tuyến đóng vai trò quyết định trong việc hình thành ứng suất sấy dư thừa và khuyết tật nứt bề mặt, đòi hỏi phải mô tả chính xác cả hai phương dẫn nhiệt ($k_l, k_r$) và truyền ẩm ($D$).

Khi định vị với các công bố quốc tế tiêu biểu:

  • Nghiên cứu của Thitinan Sattho và Ram Yamsaengsung (2005) trên gỗ Cao su (Hevea brasiliensis) ứng dụng sấy chân không kết hợp hơi quá nhiệt ở áp suất $190 - 210\text{ mmHg}$ và dải nhiệt $60 - 80^\circ\text{C}$ đạt thời gian sấy 24 giờ từ độ ẩm $80%$ về $12%$. Tuy nhiên, phương pháp sử dụng hơi quá nhiệt đòi hỏi hệ thống nồi hơi áp lực cao và khó kiểm soát sự đồng đều màu sắc đối với các loại gỗ chứa tannin nhạy cảm nhiệt.
  • Công trình của Xian Jun Li, Bi Guang Zhang và Wen Jun Li (2008) về sấy chân không vi sóng gỗ Thông Masson (Pinus massoniana) với công suất $115\text{ kW/m}^3$ tại $0,06\text{ MPa}$ đã phân định rõ ba giai đoạn động học sấy, nhưng gặp trở ngại lớn về nguy cơ tạo điểm quá nhiệt cục bộ (hotspots) do sự phân bố trường điện từ không đồng đều trong các khối gỗ dày.
  • Nghiên cứu của Guler và Dilek (2020) khi sấy chân không cao tần cho gỗ Sồi (Quercus petraea) và Óc chó (Juglans regia) chỉ ra rằng tốc độ thoát ẩm quá nhanh ở nhiệt độ cao làm suy giảm độ bền uốn tĩnh từ 17% đến 20% do biến tính cấu trúc sợi gỗ.

Luận án của Bùi Thị Thiên Kim đã định vị xuất sắc khoảng trống công nghệ bằng việc lựa chọn phương thức gia nhiệt bức xạ hồng ngoại kết hợp buồng hút chân không liên tục. Bức xạ hồng ngoại với bước sóng cực đại $\lambda_{max} = 2,5 - 3,5\text{ }\mu\text{m}$ (xác lập qua định luật Wien $\lambda_{max} = 2886/T$) kích thích trực tiếp dao động cộng hưởng của liên kết $\text{O}-\text{H}$ trong phân tử nước tự do, biến nước thành "vật đen" hấp thụ năng lượng tối đa trong khi chất gỗ hữu cơ tương đối trong suốt với bức xạ. Cơ chế này tạo nên sự gia nhiệt thể tích êm dịu, duy trì nhiệt độ sấy thấp ($58,5 - 58,9^\circ\text{C}$) ở áp suất $0,2\text{ bar}$, triệt tiêu hiện tượng hóa lỏng bít tắc của gôm nhựa Căm xe, giúp rút ngắn thời gian sấy mà vẫn bảo toàn hoàn hảo cơ tính và màu sắc tự nhiên của gỗ.

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

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

Luận án đóng góp quan trọng vào lý thuyết khoa học gỗ và kỹ thuật sấy thông qua việc mở rộng hệ phương trình truyền nhiệt - truyền chất Luikov cho đối tượng gỗ cứng nhiệt đới có hàm lượng trích ly cao. Nghiên cứu đã thực nghiệm và xác lập hệ thống các phương trình hàm số thực nghiệm biểu diễn chính xác mối quan hệ phi tuyến giữa các đặc tính nhiệt vật lý của gỗ Căm xe với độ ẩm $W \in [10%, 40%]$:

  1. Khối lượng riêng thể tích: $$\rho_{wo} = 0,0068 \cdot W + 0,7875 \quad (R^2 = 0,9532)$$ (trong đó $W$ là độ ẩm %, $\rho_{wo}$ có đơn vị $\text{g/cm}^3$)

  2. Nhiệt dung riêng: $$C_{wop} = -0,058 \cdot W + 1538,2 \quad (R^2 = 0,9998)$$ (trong đó $C_{wop}$ có đơn vị $\text{J/kg}\cdot\text{K}$)

  3. Hệ số dẫn nhiệt dọc thớ ($k_l$) và ngang thớ ($k_r$): $$k_l = -0,0001 \cdot W + 0,1995 \quad (R^2 = 0,9776)$$ $$k_r = 0,0002 \cdot W + 0,2083 \quad (R^2 = 0,9920)$$ (trong đó $k_l, k_r$ có đơn vị $\text{W/m}\cdot\text{K}$)

  4. Các điểm giới hạn ẩm học đặc trưng:

    • Điểm bão hòa thớ gỗ trung bình: $W_{woFSP} = 20,556%$ (dao động trong khoảng $18,14 - 23,07%$).
    • Độ ẩm thăng bằng trong điều kiện khí hậu TP. Hồ Chí Minh ($T = 32^\circ\text{C}$, $\varphi = 75%$): $W_{woEQ} = 14%$.
    • Độ ẩm ban đầu tự nhiên của nguyên liệu: $W_{woIN} = 40,36%$ (dao động trong khoảng $39,21 - 41,63%$).

Sự phát hiện ra xu hướng suy giảm nhẹ của hệ số dẫn nhiệt dọc thớ $k_l$ theo độ ẩm và sự gia tăng của hệ số dẫn nhiệt ngang thớ $k_r$ là một đóng góp lý thuyết mới, phản ánh sự biến đổi cấu trúc pha nước - không khí - vách tế bào trong hệ thống mao dẫn phức tạp của loài Xylia xylocarpa.

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

Khung phân tích của luận án là sự hợp nhất đa ngành giữa ba trụ cột lý thuyết:

  1. Nhiệt động học pha chân không: Ứng dụng phương trình trạng thái và nguyên lý Clapeyron-Clausius để xác định áp suất chân không vận hành ($p = 0,2\text{ bar}$ tương đương khoảng $150\text{ mmHg}$), tại đó nhiệt độ sôi của nước hạ xuống mức xấp xỉ $60^\circ\text{C}$, tạo động lực chênh lệch áp suất hơi nước bão hòa giữa tâm gỗ và bề mặt mà không làm cháy gỗ.
  2. Quang nhiệt học bức xạ hồng ngoại: Tích hợp định luật bức xạ Stefan-Boltzmann và mô hình hấp thụ bức xạ chọn lọc của liên kết $\text{O}-\text{H}$ trong nước, tính toán thông lượng nhiệt bức xạ $q_R = \varepsilon_{eff} \sigma (T_R^4 - T_s^4)$ truyền trực tiếp vào lớp bề mặt và tán xạ vào các lớp tế bào lân cận.
  3. Cơ học môi trường liên tục xốp: Mô hình hóa phương trình vi phân truyền nhiệt và ẩm liên kết hai chiều (2D) và ba chiều (3D): $$\rho C_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + m_v \Delta h_{evap}$$ $$\frac{\partial W}{\partial t} = \nabla \cdot (D \nabla W) - m_v$$

Điều kiện biên (boundary conditions) được thiết lập chặt chẽ với điều kiện biên loại 2 (thông lượng nhiệt bức xạ từ nguồn đèn hồng ngoại công suất $P_{hn}$) và điều kiện biên loại 3 (trao đổi ẩm bề mặt với môi trường chân không thông qua hệ số truyền khối $h_m$). Giới hạn biên điều kiện áp dụng rõ ràng cho các dải quy cách phôi gỗ xẻ có chiều dày từ 20 mm đến 50 mm, chiều rộng từ 50 mm đến 150 mm và chiều dài danh định 500 mm.

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ủ thế giới quan thực chứng (positivism) kết hợp nguyên lý kiểm chứng thực nghiệm chặt chẽ (empirical validation). Thiết kế nghiên cứu được thực hiện theo quy trình đa tầng khép kín:

[Khảo sát đặc tính thực nghiệm] 
[Mô hình hóa lý thuyết] 
[Giải số & Mô phỏng FEM] 
[Thực nghiệm kiểm chứng] 
[Quy hoạch thực nghiệm RSM] 
[Tối ưu hóa đa mục tiêu]

Quy trình nghiên cứu rigorous

  • Thu thập mẫu và chuẩn bị nguyên liệu: Căm xe (Xylia xylocarpa) nhóm II được tuyển chọn từ các lóng gỗ tròn chuẩn đường kính $0,5 - 0,7\text{ m}$, không bị khuyết tật tự nhiên (mục mục, sâu đục, nứt tâm). Mẫu được xẻ thành 4 nhóm kích thước đại diện cho công nghiệp đồ gỗ: $20 \times 50 \times 500\text{ mm}$, $50 \times 50 \times 500\text{ mm}$, $20 \times 150 \times 500\text{ mm}$, và $50 \times 150 \times 500\text{ mm}$.
  • Quy trình đo nhiệt vật lý:
    • Khối lượng riêng xác định theo phương pháp thể tích - khối lượng chuẩn trên mẫu $20 \times 20 \times 30\text{ mm}$ ở các dải ẩm từ sấy khô tuyệt đối đến bão hòa nước.
    • Nhiệt dung riêng xác định bằng bình nhiệt lượng kế đoạn nhiệt chất lượng cao, cân bằng nhiệt với môi trường nước chuẩn.
    • Hệ số dẫn nhiệt dọc thớ ($k_l$) và ngang thớ ($k_r$) đo bằng thiết bị truyền nhiệt bản phẳng chuyên dụng trong buồng cách nhiệt có đầu dò nhiệt điện trở PT100 độ nhạy cao.
  • Hệ thống thiết bị sấy chân không thực nghiệm: Buồng sấy chân không hình trụ chế tạo bằng thép không gỉ chịu áp lực cao, trang bị bơm hút chân không vòng nước/dầu đạt độ chân không tuyệt đối tới $0,01\text{ bar}$, giàn đèn bức xạ hồng ngoại công suất điều khiển biến tần linh hoạt, hệ thống cảm biến đo nhiệt độ đa điểm gắn sâu vào tâm và bề mặt thanh gỗ, hệ thống cân khối lượng điện tử online tích hợp bộ ghi dữ liệu tự động (Data Logger).

Data và phân tích

  • Kỹ thuật mô phỏng số học: Giải hệ phương trình vi phân đạo hàm riêng liên kết bằng phương pháp phần tử hữu hạn (FEM) trên nền tảng phần mềm COMSOL Multiphysics. Miền tính toán được chia lưới tự do với các phần tử biên cực mịn (boundary layer mesh) nhằm nắm bắt chính xác gradient nhiệt và ẩm cực lớn tại bề mặt tiếp xúc bức xạ.
  • Phương pháp quy hoạch thực nghiệm: Ứng dụng phương pháp quy hoạch thực nghiệm trực giao bậc 1 và bậc 2 (Response Surface Methodology - RSM) với 2 yếu tố đầu vào chính:
    • $X_1$: Nhiệt độ sấy bề mặt ($T_s \in [50, 70]^\circ\text{C}$)
    • $X_2$: Cường độ bức xạ hồng ngoại ($P_{hn} \in [500, 800]\text{ W/m}^2$)
    • Áp suất chân không được khống chế tương ứng tại $p = 0,2\text{ bar}$.
  • Thông số đầu ra (Hàm mục tiêu):
    • $Y_{tg}$: Thời gian sấy từ độ ẩm ban đầu $W_{woIN} \approx 40%$ về độ ẩm cuối $W_{end} = 10 - 12%$ (giờ).
    • $Y_{kt}$: Tỷ lệ thể tích gỗ khuyết tật sau sấy (%) tính theo tổng diện tích và thể tích các vết nứt mặt, nứt đầu, cong vênh theo tiêu chuẩn kiểm tra chất lượng lâm sản.
  • Kiểm định thống kê: Toàn bộ các mô hình hồi quy bậc 1 và bậc 2 đều được đánh giá ý nghĩa thống kê thông qua phân tích phương sai (ANOVA), kiểm định Fisher ($F$-test) và giá trị $p < 0,05$, hệ số tương quan $R^2$ đạt từ 0,92 đến 0,98.

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

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

  1. Quy luật phân bố nhiệt - ẩm không đồng đều: Kết quả mô phỏng FEM và thực nghiệm chứng minh rằng trong giai đoạn đầu (0 - 15 phút), nhiệt độ bề mặt thanh gỗ nhận bức xạ hồng ngoại tăng nhanh lên mức $58 - 60^\circ\text{C}$, trong khi tâm gỗ có độ trễ nhiệt nhất định. Tuy nhiên, nhờ môi trường chân không $0,2\text{ bar}$, nước tự do tại tâm nhanh chóng đạt điểm sôi thấp, tạo áp suất nội tại đẩy dòng ẩm di chuyển hướng tâm ra ngoài theo chiều ngược với dòng nhiệt dẫn vào, làm giảm đáng kể hiện tượng "cháy vỏ" (case hardening).
  2. Khống chế sai số mô hình dưới 5%: Sự so sánh giữa đường cong phân bố nhiệt độ tâm và đường cong giảm ẩm thực nghiệm với kết quả tính toán từ mô hình toán giải bằng COMSOL Multiphysics cho thấy sự trùng khớp cao về biên dạng và xu thế động học, với sai số trung bình lớn nhất trên tất cả các quy cách mẫu đều dưới 5%.
  3. Mối quan hệ phi tuyến giữa chế độ công nghệ và khuyết tật gỗ: Phân tích biểu đồ 3D mặt đáp ứng và đường đồng mức (contour plots) chỉ ra rằng khi tăng cường độ bức xạ $P_{hn} > 700\text{ W/m}^2$ và nhiệt độ $T_s > 65^\circ\text{C}$, thời gian sấy $Y_{tg}$ giảm mạnh nhưng tỷ lệ khuyết tật $Y_{kt}$ tăng vọt phi tuyến do gradient ẩm quá lớn vượt quá giới hạn ứng suất bền kéo ngang thớ của gỗ Căm xe.
  4. Bộ thông số công nghệ sấy tối ưu đa mục tiêu: Luận án đã giải bài toán tối ưu hàm tổng hợp ($Y_{sum}$) cho từng dạng kích thước mẫu, xác định chính xác các thông số công nghệ sấy chân không hồng ngoại phù hợp nhất:
    • Quy cách $20 \times 50 \times 500\text{ mm}$: $T_s = 58,5^\circ\text{C}$; $p = 0,2\text{ bar}$; $P_{hn} = 625\text{ W/m}^2$.
    • Quy cách $50 \times 50 \times 500\text{ mm}$: $T_s = 58,7^\circ\text{C}$; $p = 0,2\text{ bar}$; $P_{hn} = 633\text{ W/m}^2$.
    • Quy cách $20 \times 150 \times 500\text{ mm}$: $T_s = 58,6^\circ\text{C}$; $p = 0,2\text{ bar}$; $P_{hn} = 630\text{ W/m}^2$.
    • Quy cách $50 \times 150 \times 500\text{ mm}$: $T_s = 58,9^\circ\text{C}$; $p = 0,2\text{ bar}$; $P_{hn} = 641\text{ W/m}^2$.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp cơ sở khoa học hoàn chỉnh về cơ chế động học truyền nhiệt và thoát ẩm của gỗ cứng nhiệt đới nhiều dầu nhựa trong môi trường chân không bức xạ hồng ngoại, lấp đầy khoảng trống dữ liệu nhiệt vật lý thực nghiệm của loài Căm xe trong cơ sở dữ liệu lâm sản thế giới.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn mực kết hợp giữa mô hình hóa giải tích số học (COMSOL FEM) và quy hoạch thực nghiệm đa mục tiêu (RSM-ANOVA), có thể chuyển giao và áp dụng trực tiếp cho việc nghiên cứu chế độ sấy các loài gỗ quý khó sấy khác (gỗ Hương, Gõ đỏ, Cẩm lai, Mun).
  • Về mặt ứng dụng thực tiễn: Rút ngắn thời gian sấy gỗ Căm xe xuống nhiều lần so với sấy đối lưu thông thường (từ hàng tháng xuống còn vài chục giờ tùy chiều dày phôi), giảm tỷ lệ khuyết tật nứt nẻ xuống dưới 3 - 5%, tiết kiệm đáng kể chi phí năng lượng và diện tích mặt bằng nhà xưởng.
  • Về mặt chính sách và công nghiệp: Cung cấp luận cứ khoa học giúp các nhà hoạch định chính sách ngành chế biến gỗ xây dựng quy chuẩn kỹ thuật quốc gia về công nghệ sấy tiết kiệm năng lượng, khuyến khích các doanh nghiệp đầu tư đổi mới công nghệ từ sấy lò hơi truyền thống sang công nghệ chân không hồng ngoại tự động hóa cao.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, công trình nghiên cứu vẫn tồn tại một số giới hạn học thuật cần được mở rộng:

  1. Giới hạn quy mô thực nghiệm: Nghiên cứu được thực hiện trên hệ thống thiết bị sấy chân không quy mô phòng thí nghiệm (pilot) với chiều dài mẫu danh định 500 mm; các hiệu ứng truyền khối dọc trục trên phôi gỗ dài công nghiệp ($2,0 - 3,0\text{ m}$) cần được hiệu chỉnh thêm về độ suy giảm bức xạ theo khoảng cách.
  2. Thiếu vắng đo lường ứng suất thực thời (real-time stress): Luận án đánh giá tỷ lệ khuyết tật sau khi kết thúc quá trình sấy mà chưa tích hợp cảm biến đo trực tiếp sự biến thiên ứng suất kéo/nén trong vách tế bào gỗ trong suốt thời gian thực của mẻ sấy.
  3. Chưa khảo sát đa dạng biến thiên xuất xứ gỗ: Nguồn mẫu Căm xe nghiên cứu chủ yếu đại diện cho nhóm gỗ nhập khẩu phổ biến, chưa bao quát toàn bộ sự sai khác về cấu tạo giải phẫu và hàm lượng gôm nhựa theo từng tiểu vùng lập địa sinh thái khác nhau.

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

  • Phát triển mô hình mô phỏng 3D mở rộng có tích hợp trường ứng suất cơ học đàn - nhớt - dẻo (viscoelastic stress modeling) để dự báo chính xác thời điểm xuất hiện vết nứt vi mô.
  • Nghiên cứu thiết kế hệ thống sấy chân không hồng ngoại liên tục quy mô công nghiệp ($10 - 20\text{ m}^3/\text{mẻ}$) với giàn đèn hồng ngoại điều khiển độc lập theo từng vùng nhiệt độ (multi-zone PID control).
  • Ứng dụng trí tuệ nhân tạo (AI/Machine Learning) kết hợp thuật toán tối ưu hóa bầy đàn (PSO/GA) để tự động điều chỉnh thông số sấy online theo tín hiệu cảm biến độ ẩm không phá hủy.

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

  • Tác động học thuật: Dự báo công trình sẽ trở thành tài liệu tham khảo cốt lõi cho các nghiên cứu sinh, học viên cao học chuyên ngành Kỹ thuật Chế biến Lâm sản, Kỹ thuật Nhiệt và Khoa học Vật liệu Màng xốp. Bộ dữ liệu nhiệt vật lý có tiềm năng trích dẫn cao trong các tạp chí uy tín thuộc danh mục Scopus/Web of Science trong lĩnh vực Drying Technology, Wood Science and Technology, BioResources.
  • Chuyển đổi công nghiệp: Thúc đẩy cuộc cách mạng công nghệ trong ngành sản xuất ván sàn và đồ gỗ xuất khẩu cao cấp tại Việt Nam. Giúp các doanh nghiệp chế biến gỗ giải quyết dứt điểm tình trạng tồn đọng vốn lưu động do chu kỳ sấy đối lưu quá dài, nâng cao tỷ lệ lợi dụng gỗ nguyên liệu thêm 10 - 15%.
  • Lợi ích xã hội và môi trường: Công nghệ sấy chân không hồng ngoại sử dụng 100% năng lượng điện sạch, loại bỏ hoàn toàn khí thải $\text{CO}_2, \text{SO}_x$ và tro xỉ từ các lò hơi đốt củi/than truyền thống, đóng góp trực tiếp vào cam kết Net Zero và lộ trình phát triển bền vững của ngành công nghiệp chế biến gỗ Việt Nam.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận nghiên cứu kết hợp giữa mô hình hóa vật lý - toán học, giải số phần tử hữu hạn và quy hoạch thực nghiệm tối ưu hóa; thừa hưởng bộ thông số nhiệt vật lý chuẩn của gỗ Căm xe.
  • Kỹ sư R&D và Nhà chế tạo thiết bị: Sở hữu cơ sở tính toán nhiệt, tính chọn công suất giàn đèn hồng ngoại, công suất bơm hút chân không và quy trình công nghệ 8 bước để thiết kế chế tạo máy sấy chân không thương mại.
  • Doanh nghiệp chế biến gỗ: Áp dụng trực tiếp chế độ sấy chuẩn ($T_s, p, P_{hn}$) cho 4 quy cách phôi ván sàn Căm xe, giúp nâng cao chất lượng sản phẩm xuất khẩu, giảm thiểu khiếu nại về nứt tét hay cong vênh từ khách hàng quốc tế.
  • Cơ quan quản lý nhà nước: Sử dụng các số liệu khoa học để xây dựng định mức kinh tế kỹ thuật, tiêu chuẩn tiêu hao năng lượng trong công nghiệp chế biến lâm sản.

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ì? Đó là việc mở rộng hệ phương trình truyền nhiệt - truyền chất Luikov cho vật liệu gỗ cứng nhiệt đới có tỷ trọng cao và giàu gôm nhựa (Xylia xylocarpa), thông qua việc xác lập đồng bộ 4 phương trình tương quan nhiệt vật lý ($\rho_{wo}, C_{wop}, k_l, k_r$) phụ thuộc trực tiếp vào độ ẩm $W \in [10%, 40%]$ với độ tin cậy vượt trội ($R^2 \ge 0,9532$).

  2. Điểm mới về mặt phương pháp luận so với các nghiên cứu tiền nhiệm là gì? So với mô hình sấy chu kỳ đối lưu của Assouad và Jomaa (2003) hay mô hình thế nước tiếp xúc của Maurice Defo et al. (2000), luận án đã tích hợp thành công mô hình bức xạ hồng ngoại bước sóng ngắn $\lambda_{max} = 2,5 - 3,5\text{ }\mu\text{m}$ vào hệ phương trình giải số phần tử hữu hạn trên COMSOL Multiphysics, cho phép mô phỏng chính xác động học phân bố nhiệt - ẩm với sai số dưới 5% so với thực nghiệm.

  3. Phát hiện thực nghiệm nào gây bất ngờ nhất trong nghiên cứu? Đó là phát hiện về quy luật biến thiên trái chiều của hệ số dẫn nhiệt: trong khi hệ số dẫn nhiệt ngang thớ $k_r$ tăng thuận theo độ ẩm ($k_r = 0,0002 \cdot W + 0,2083$), thì hệ số dẫn nhiệt dọc thớ $k_l$ lại có xu hướng giảm nhẹ khi độ ẩm tăng ($k_l = -0,0001 \cdot W + 0,1995$). Điều này được giải thích bởi cấu trúc vách tế bào sợi gỗ dày đặc và sự hiện diện của dịch bào gôm nhựa trong các tế bào mô mềm quanh mạch.

  4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không? Luận án cung cấp chi tiết toàn bộ điều kiện biên, thông số vật liệu, lưu đồ thuật toán giải số, mã lệnh thiết lập trên COMSOL, sơ đồ nguyên lý buồng sấy và ma trận thực nghiệm trực giao bậc 2 cùng các bảng phân tích ANOVA, đảm bảo khả năng tái lập thực nghiệm và mô phỏng với độ chính xác tuyệt đối.

  5. Chương trình nghiên cứu 10 năm tới mở ra từ luận án này là gì? Xây dựng hệ sinh thái công nghệ sấy chân không thông minh (Smart Vacuum Drying 4.0) bao gồm: thư viện số hóa đặc tính nhiệt vật lý 50 loài gỗ nhiệt đới, phần mềm chuyên dụng dự báo nứt gỗ dựa trên AI và chế tạo buồng sấy chân không hồng ngoại liên tục quy mô công nghiệp lớn.

Kết luận

Luận án tiến sĩ kỹ thuật của NCS. Bùi Thị Thiên Kim là một công trình khoa học mẫu mực, kết hợp nhuần nhuyễn giữa lý thuyết cơ bản sâu sắc, kỹ thuật mô phỏng số tiên tiến và thực nghiệm công nghệ chính xác. Sáu đóng góp cụ thể của luận án bao gồm:

  1. Xác định chính xác các thông số nhiệt vật lý cơ bản và hàm số tương quan thực nghiệm theo độ ẩm cho gỗ Căm xe (Xylia xylocarpa).
  2. Xây dựng thành công mô hình vật lý và mô hình toán học mô tả bản chất truyền nhiệt và truyền ẩm liên kết trong môi trường chân không bức xạ hồng ngoại.
  3. Ứng dụng phương pháp phần tử hữu hạn trên phần mềm COMSOL Multiphysics mô phỏng trực quan diễn biến phân bố trường nhiệt và ẩm với sai số thực nghiệm $< 5%$.
  4. Thiết lập hệ phương trình hồi quy bậc hai mô tả ảnh hưởng của nhiệt độ và cường độ bức xạ hồng ngoại đến thời gian sấy và tỷ lệ khuyết tật bằng phương pháp quy hoạch thực nghiệm RSM.
  5. Xác lập bộ thông số công nghệ sấy tối ưu đa mục tiêu ($T_s = 58,5 - 58,9^\circ\text{C}$; $p = 0,2\text{ bar}$; $P_{hn} = 625 - 641\text{ W/m}^2$) cho các dải quy cách chiều dày 20 - 50 mm.
  6. Đề xuất quy trình công nghệ sấy chân không hồng ngoại hoàn chỉnh, mở ra hướng ứng dụng công nghiệp đột phá nhằm nâng cao giá trị gia tăng cho ngành chế biến lâm sản Việt Nam trên trường quốc tế.