Tổng quan về luận án

Ngành công nghiệp chế biến gỗ Việt Nam đang đứng trước yêu cầu cấp thiết về việc nâng cao giá trị gia tăng và hiệu quả sử dụng nguồn nguyên liệu gỗ rừng trồng mọc nhanh, đặc biệt là loài Keo tai tượng (Acacia mangium Willd.). Trong chuỗi giá trị chế biến gỗ xẻ, công đoạn sấy gỗ chiếm vị trí then chốt, quyết định trực tiếp đến độ bền cơ lý, tính ổn định kích thước và chất lượng sản phẩm cuối cùng. Tuy nhiên, gỗ Keo tai tượng thuộc nhóm gỗ có khối lượng riêng trung bình nhưng rất dễ phát sinh khuyết tật sấy nghiêm trọng như nứt mặt, nứt đầu, nứt ngầm, mo móp tế bào (collapse) và cong vênh do cấu tạo giải phẫu phức tạp và mức độ dị hướng cao. Động lực căn bản chi phối toàn bộ quá trình thoát nước cũng như sự hình thành nội ứng suất phá hủy trong gỗ là sự chênh lệch độ ẩm giữa các lớp theo chiều dày, được lượng hóa thông qua dốc độ ẩm (moisture content gradient - MC gradient).

Khoảng trống nghiên cứu (research gap) cốt lõi trong y văn hiện nay nằm ở việc thiếu vắng các công trình định lượng một cách đồng bộ và chính xác mối tương quan động học giữa dốc độ ẩm đa điểm theo thời gian thực với sự phát sinh khuyết tật trong gỗ xẻ Keo tai tượng nhiệt đới. Trong khi các nghiên cứu trước đây tại Việt Nam chủ yếu tiếp cận theo phương pháp thử sai thực nghiệm hoặc sử dụng ẩm kế đo điện trở vốn có độ sai lệch rất lớn khi độ ẩm gỗ nằm trên điểm bão hòa thớ gỗ (FSP > 30%) như chỉ ra bởi Nguyễn Xuân Hiên (2006) và Đỗ Văn Bản (2012), thì các cơ chế chuyển khối vi mô (thấm lỏng Darcy và khuếch tán Fick) chưa được đo đạc bằng các thiết bị chuyên dụng tiêu chuẩn quốc tế. Hơn nữa, việc ứng dụng mô hình toán học số trị để mô phỏng diễn biến ẩm và tối ưu hóa chế độ sấy tiết kiệm năng lượng bằng năng lượng mặt trời (NLMT) cho gỗ Keo tai tượng vẫn chưa được giải quyết thấu đáo.

Để giải quyết triệt để khoảng trống này, luận án tập trung vào các câu hỏi nghiên cứu và giả thuyết khoa học sau:

  • Câu hỏi nghiên cứu 1 (RQ1): Các đặc điểm cấu tạo hiển vi, siêu hiển vi và đặc tính chuyển khối (hệ số thấm $K$, hệ số khuếch tán $D_b$, điểm bão hòa thớ gỗ $FSP$, độ rỗng $\phi$) của gỗ Keo tai tượng tác động như thế nào đến quá trình vận chuyển ẩm?
  • Câu hỏi nghiên cứu 2 (RQ2): Mức chênh lệch độ ẩm theo chiều dày tấm ván biến thiên ra sao dưới các chế độ sấy có dốc sấy ($U$) khác nhau và tương quan định lượng với sự phát triển khuyết tật nứt vỡ, biến dạng như thế nào?
  • Câu hỏi nghiên cứu 3 (RQ3): Mô hình toán học truyền nhiệt - truyền ẩm hai chiều TransPore có khả năng dự đoán chính xác trường độ ẩm và thời gian sấy gỗ Keo tai tượng trong lò sấy quy chuẩn hay không?
  • Câu hỏi nghiên cứu 4 (RQ4): Cơ chế sấy gián đoạn trong lò sấy NLMT tác động như thế nào đến sự thư giãn nội ứng suất và làm thế nào để xây dựng quy trình công nghệ sấy NLMT tối ưu về kỹ thuật và kinh tế?

Hệ thống giả thuyết nghiên cứu bao gồm:

  • Giả thuyết 1 (H1): Quá trình chuyển khối của gỗ Keo tai tượng bị kiểm soát bởi tính dị hướng giải phẫu rõ rệt giữa các phương dọc thớ ($L$), xuyên tâm ($R$) và tiếp tuyến ($T$), trong đó tính thấm dọc thớ vượt trội do cấu trúc thông suốt của mạch gỗ không có thể bít.
  • Giả thuyết 2 (H2): Tồn tại các ngưỡng dốc độ ẩm tới hạn ($\Delta MC_{max}$) và độ ẩm thăng bằng tối thiểu ($EMC_{min}$) tại các mốc độ ẩm gỗ then chốt ($MC = 50%$, $20%$, $12%$), nếu vượt qua các ngưỡng này ứng suất kéo bề mặt sẽ vượt giới hạn bền kéo ngang thớ sinh ra nứt vỡ không thể phục hồi.
  • Giả thuyết 3 (H3): Việc ứng dụng mô hình số trị TransPore 2D giải bằng phương pháp thể tích hữu hạn cho phép mô phỏng chính xác động học sấy với sai số thời gian sấy ($E_{dt}$) dưới $10%$.
  • Giả thuyết 4 (H4): Chu kỳ bức xạ gián đoạn ngày - đêm của lò sấy NLMT tạo ra giai đoạn tự cân bằng dốc ẩm ban đêm, giúp triệt tiêu ứng suất dư và cho phép sấy gỗ Keo tai tượng đạt chất lượng loại A với chi phí năng lượng giảm đáng kể so với sấy hơi nước.

Khung lý thuyết của công trình tích hợp lý thuyết dòng chảy tầng Darcy trong môi trường xốp mao quản bão hòa ẩm, định luật khuếch tán Fick I và Fick II cho nước liên kết trong vách tế bào, lý thuyết ứng suất - biến dạng viscoelastic của McMillen (1958) và lý thuyết chuyển khối liên hợp trong môi trường xốp (Porous Media Transport Theory) của Perré & Turner (1999).

Đóng góp đột phá của luận án được lượng hóa thông qua việc xác định chính xác bộ thông số vật lý và chuyển khối vi mô của gỗ Keo tai tượng Việt Nam, thiết lập hệ phương trình hồi quy tương quan phi tuyến giữa dốc ẩm với mức độ nứt vỡ ($R^2 > 0,85$), và chuẩn hóa thành công quy trình công nghệ sấy NLMT cho ván xẻ dày 25 mm với chu kỳ sấy từ 14 đến 18 ngày, giảm chi phí năng lượng từ 30% đến 45% so với lò sấy hơi nước truyền thống, nâng tỷ lệ thu hồi gỗ đạt tiêu chuẩn chất lượng cao lên trên 90%.

Phạm vi nghiên cứu tập trung vào gỗ Keo tai tượng 9 năm tuổi khai thác tại xã Ngọc Thanh, thành phố Phúc Yên, tỉnh Vĩnh Phúc. Quy cách ván xẻ nghiên cứu có chiều dày chuẩn $(25 \pm 1)\text{ mm}$, rộng từ $100\text{ mm}$ đến $120\text{ mm}$, dài $1000\text{ mm}$ đến $2000\text{ mm}$. Nghiên cứu thực nghiệm được đối chứng giữa lò sấy quy chuẩn chuyên dụng Melb Uni 2074-4 và lò sấy NLMT thí nghiệm cũng như lò sấy NLMT khảo nghiệm quy mô bán công nghiệp $30\text{ m}^3\text{/mẻ}$ và $35\text{ m}^3\text{/mẻ}$.

Literature Review và Positioning

Nghiên cứu về động học sấy gỗ và cơ chế phát sinh khuyết tật đã trải qua hơn sáu thập kỷ phát triển với các mốc tiến hóa lý thuyết quan trọng. Dòng nghiên cứu nền tảng khởi nguồn từ công trình kinh điển của McMillen (1958), người đầu tiên mô tả mối liên hệ biện chứng giữa dốc độ ẩm và sự đảo chiều ứng suất (stress reversal) trong quá trình sấy gỗ Sồi đỏ (Quercus rubra). McMillen đã chỉ ra rằng ở giai đoạn đầu khi bề mặt ván giảm xuống dưới điểm bão hòa thớ gỗ ($FSP \approx 30%$) trong khi tâm gỗ vẫn giữ độ ẩm cao gần 90%, bề mặt bị cản trở co rút bởi lớp lõi, sinh ra ứng suất kéo lớn nhất ở vỏ và ứng suất nén ở lõi; khi sấy tiếp tục, sự phân bố ứng suất bị đảo ngược hoàn toàn. Tiếp nối nền tảng đó, Stamm (1967), Siau (1984), Skaar (1988) và Walker (2006) đã hoàn thiện lý thuyết chuyển khối nhiệt động lực học trong gỗ, phân định ranh giới rõ ràng giữa dòng chuyển dịch chất lỏng tự do chịu tác động của áp suất thủy tĩnh (định luật Darcy) và dòng khuếch tán hơi nước/nước liên kết chịu tác động của građien nồng độ và nhiệt độ (định luật Fick).

Trong y văn tồn tại hai cuộc tranh luận khoa học lớn:

  1. Tranh luận về phương pháp đo lường dốc độ ẩm: Một trường phái ủng hộ các phương pháp đo gián tiếp không phá hủy như ẩm kế điện trở nhiều cực (Forrer, 1984; Simpson, 1986), cảm biến độ ẩm tương đối trong lỗ khoan (Melin et al., 2016), hoặc chụp quét chùm tia X và cộng hưởng từ (Cai, 2007; Tiitta et al., 2010). Tuy nhiên, các kỹ thuật này bộc lộ nhược điểm nghiêm trọng khi gỗ ở độ ẩm cao trên $FSP$ do hiện tượng bão hòa điện trở và chi phí thiết bị đắt đỏ. Trường phái thứ hai (McCurdy, 2006; Pang, 1996; Northway, 2001; Phonetip, 2018) khẳng định phương pháp cắt lát mẫu vi phân (slicing method) theo chiều dày kết hợp sấy cân khối lượng vẫn là "tiêu chuẩn vàng" đem lại độ tin cậy tuyệt đối cho toàn dải độ ẩm từ tươi đến khô kiệt. Luận án đã định vị phương pháp luận của mình theo trường phái cắt lát vi phân chính xác cao.
  2. Tranh luận về mô hình hóa toán học: Tranh biện giữa các mô hình giải tích một chiều đơn giản dựa thuần túy trên hệ số khuếch tán hiệu dụng (Haque, 2002) đối lập với các mô hình số trị đa chiều tích hợp truyền nhiệt - truyền khối liên hợp trong môi trường xốp dị hướng (Perré & Turner, 1999; Perré & Passard, 2004; Salin, 2010). Các mô hình một chiều bỏ qua tính dị hướng và chuyển khối lỏng trên $FSP$, dẫn đến sai lệch lớn khi dự đoán ứng suất sấy ở gỗ lá rộng.

Vị thế học thuật của luận án được xác lập rõ ràng bằng việc bắc cầu giữa lý thuyết chuyển khối môi trường xốp hiện đại và thực tiễn công nghệ sấy gỗ rừng trồng nhiệt đới. So với hai nghiên cứu quốc tế tiêu biểu:

  • Nghiên cứu của Redman (2017) tại Đại học Melbourne và DAF Queensland (Úc) đã xây dựng mô hình TransPore 2D cho 4 loài gỗ lá rộng ôn đới bản địa Úc trong điều kiện sấy chân không và sấy quy chuẩn. Luận án của Hà Tiến Mạnh đã kế thừa, tùy biến và mở rộng hệ thống mã nguồn này để ứng dụng cho loài Keo tai tượng (Acacia mangium) nhiệt đới có giải phẫu vách tế bào mỏng và tốc độ thoát ẩm nhanh.
  • Nghiên cứu của Phonetip (2018) tại Úc về sấy năng lượng mặt trời gỗ Bạch đàn (Eucalyptus) tập trung vào đánh giá ứng suất dư theo tiêu chuẩn AS/NZS 4787:2001. Luận án đã tiến xa hơn khi không chỉ dừng lại ở ứng suất dư sau sấy mà thiết lập được các hàm tương quan động học giữa dốc ẩm tại từng thời điểm sấy với sự xuất hiện của các vết nứt mặt, nứt đầu và mo móp, từ đó lượng hóa thành các ngưỡng công nghệ điều khiển lò sấy NLMT.

Tại Việt Nam, các công trình của Hồ Thu Thủy (2005) trên gỗ Chò chỉ, Dầu gió hay Đỗ Văn Bản (2012) trên gỗ Bạch đàn trắng đã có những đóng góp bước đầu về hệ số chuyển khối và giải phẫu hiển vi. Tuy nhiên, luận án này là công trình đầu tiên tại Việt Nam giải quyết trọn vẹn chuỗi liên kết: Giải phẫu siêu hiển vi (SEM/OM) $\rightarrow$ Đo đạc thông số chuyển khối chuyên sâu (ALU-CHA/PVC-CHA) $\rightarrow$ Đo dốc ẩm cắt lát vi phân $\rightarrow$ Mô phỏng số TransPore 2D $\rightarrow$ Chuẩn hóa quy trình công nghệ lò sấy NLMT thế hệ mới (mái vòm hấp thụ nhiệt cao do Bùi Duy Ngọc, 2021 phát triển).

Đó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 mở rộng hệ thống lý thuyết vật lý gỗ và kỹ thuật sấy trên các khía cạnh cốt lõi:

  1. Mở rộng lý thuyết chuyển khối trong gỗ lá rộng nhiệt đới: Luận án chứng minh rằng trong gỗ Keo tai tượng 9 năm tuổi, quá trình vận chuyển ẩm chịu sự chi phối mạnh mẽ bởi cấu trúc màng lỗ thông ngang không có cấu trúc nút (torus) và màng treo (margo) phân hóa rõ rệt, kết hợp với các lỗ rỗng vách tế bào (nanopores) có kích thước từ $2\text{ nm}$ đến $4\text{ nm}$. Điều này lý giải tại sao hệ số khuếch tán ẩm theo chiều dọc thớ ($D_{bL}$) lớn hơn chiều ngang thớ ($D_{bR}, D_{bT}$) từ 100 lần ở độ ẩm thấp ($MC = 5%$) xuống còn từ 2 đến 4 lần ở vùng tiệm cận điểm bão hòa thớ gỗ ($MC = 25% - 30%$).
  2. Làm sâu sắc thêm lý thuyết phát sinh ứng suất của McMillen (1958): Luận án đã lượng hóa chính xác thời điểm đảo chiều ứng suất nội tại ở gỗ Keo tai tượng xảy ra khi độ ẩm trung bình của tấm gỗ đạt mốc $MC \approx 20% - 25%$. Việc chứng minh dốc độ ẩm $\Delta MC$ là biến đại diện trực tiếp phản ánh mức độ tích tụ nội ứng suất đã cung cấp cơ sở lý thuyết vững chắc để thay thế việc đo ứng suất cơ học phức tạp bằng việc kiểm soát dốc độ ẩm và độ ẩm thăng bằng môi trường ($EMC$).
  3. Thiết lập mô hình tương quan phi tuyến giữa dốc ẩm và phá hủy vật liệu: Đề xuất 6 mệnh đề lý thuyết (propositions) liên kết giữa građien ẩm $\Delta MC$, độ ẩm cân bằng môi trường $EMC$ và tổng mức độ nứt vỡ $N_i$, tạo ra bước chuyển đổi nhận thức từ việc điều khiển sấy theo kinh nghiệm cảm tính sang điều khiển sấy thích ứng theo giới hạn biến dạng tới hạn của tế bào gỗ.

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

Khung phân tích của luận án được xây dựng dựa trên sự tích hợp đa ngành giữa ba khối lý thuyết trụ cột:

  • Khối chuyển khối vi mô (Micro-scale Transport): Dựa trên định luật Darcy cho dòng chất lỏng nhớt: $$q = -\frac{K}{\mu} \nabla P$$ và định luật Fick II cho dòng khuếch tán không ổn định: $$\frac{\partial C}{\partial t} = \nabla \cdot (D_b \nabla C)$$
  • Khối cơ học biến dạng và nứt nẻ (Deformation & Fracture Mechanics): Tương quan giữa độ co rút dị hướng ($\beta_T / \beta_R$) với tỷ lệ đường kính ruột tế bào trên chiều dày vách tế bào ($D_{fl} / T_{fw}$) giải thích hiện tượng sụp đổ vách tế bào (mo móp).
  • Khối mô phỏng trường số trị liên hợp (Finite Volume Coupled Modeling): Mô hình TransPore hai chiều mô phỏng đồng thời trường nhiệt độ $T(x,y,t)$ và trường độ ẩm $MC(x,y,t)$ trên mặt cắt ngang tấm ván xẻ (tiếp tuyến và xuyên tâm).

Điều kiện biên của khung phân tích được xác định rõ ràng: Áp dụng cho gỗ xẻ Keo tai tượng rừng trồng tuổi 9, chiều dày ván từ $20\text{ mm}$ đến $30\text{ mm}$, môi trường sấy có nhiệt độ sấy khô $T_k \le 65^\circ\text{C}$, vận tốc gió đối lưu $v = 1,5 - 2,5\text{ m/s}$, dốc sấy $U$ biến thiên trong khoảng $1,5 - 5,0$.

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 hiện thực phản biện (critical realism), xây dựng hệ phương pháp tiếp cận đa cấp độ từ cấu trúc hiển vi nano/micro đến quy mô meso (lát cắt mẫu ván) và macro (đống gỗ công nghiệp trong lò sấy). Thiết kế nghiên cứu tích hợp phương pháp thực nghiệm đo đạc vật lý chính xác cao, mô phỏng số học trên máy tính và kiểm chứng khảo nghiệm sản xuất thực tế.

Mẫu nghiên cứu gồm 15 cây Keo tai tượng 9 năm tuổi sinh trưởng đồng đều tại Vĩnh Phúc, đường kính ngang ngực $D_{1.3} \ge 20\text{ cm}$, thân thẳng, không bị sâu bệnh. Toàn bộ gỗ tròn được cưa xẻ thành các tấm ván có kích thước đồng nhất $(25 \pm 1) \times 120 \times 1000\text{ mm}$, phân loại nghiêm ngặt thành ván xẻ tiếp tuyến, ván xẻ xuyên tâm và ván hỗn hợp.

Quy trình nghiên cứu rigorous

  1. Khảo sát cấu trúc giải phẫu hiển vi và siêu hiển vi:
    • Kính hiển vi quang học (OM) chuyên dụng tại Viện Nghiên cứu Công nghiệp rừng.
    • Kính hiển vi điện tử quét (SEM) độ phân giải siêu cao tại Đại học Kyoto, Nhật Bản, phân tích chi tiết cấu trúc màng ngăn lỗ thông ngang, kích thước ruột tế bào sợi gỗ ($D_{fl}$), chiều dày vách tế bào sợi gỗ ($T_{fw}$) và đường kính mạch gỗ. Độ rỗng ruột tế bào mạch gỗ ($\phi_v$) và sợi gỗ ($\phi_f$) được phân tích định lượng tự động bằng phần mềm ImageJ và Meshpore.
  2. Xác định hệ số chuyển khối chuyên sâu:
    • Hệ số thấm lỏng và thấm khí ($K$): Thực hiện tại Bộ Nông nghiệp và Thủy sản bang Queensland (DAF - Úc) sử dụng máy đo độ thấm tự động POROLUX™ 1000 Porometer trên nguyên lý đo áp suất dòng khí nén ổn định (hệ thống ALU-CHA).
    • Hệ số khuếch tán ẩm ($D_b$): Đo bằng hệ thống Vaporimeter trên thiết bị PVC-CHA trong buồng khí hậu chuẩn, kiểm soát nhiệt độ và độ ẩm tương đối tuyệt đối, đo lường dòng hơi ẩm khuếch tán qua mẫu gỗ ở trạng thái ổn định (steady-state).
    • Khối lượng riêng cơ bản ($\rho_b$), độ co rút ($\beta$) và điểm bão hòa thớ gỗ ($FSP$): Xác định theo tiêu chuẩn quốc tế ISO 13061-1:2014 và ISO 13061-2:2014, sử dụng mẫu siêu mỏng dưới $1\text{ mm}$ để triệt tiêu hiện tượng mo móp tế bào khi xác định $FSP$.
  3. Xác định dốc độ ẩm và kiểm soát khuyết tật:
    • Xác định dốc độ ẩm theo chiều dày $\Delta MC$ bằng phương pháp cắt lát mẫu vi phân (slicing method) trên thiết bị cắt lát cơ khí chuyên dụng, chia tấm ván $25\text{ mm}$ thành các lớp mỏng từ bề mặt vào tâm, cân khối lượng ngay lập tức trên cân phân tích độ chính xác $0,0001\text{ g}$ và sấy khô kiệt ở $103 \pm 2^\circ\text{C}$.
    • Khuyết tật sấy được lượng hóa theo tiêu chuẩn AS/NZS 4787:2001 và TCVN, đo đếm chính xác chiều dài vết nứt mặt ($NM$), nứt đầu ($NĐ$), độ sâu mo móp ($C$) và độ cong vênh.
    • Kiểm định độ tin cậy: Thử nghiệm đối chứng song song giữa phương pháp cắt lát và phương pháp đo điện trở cho hệ số xác định $R^2 = 0,8795$, khẳng định tính chuẩn xác vượt trội của phương pháp cắt lát.

Data và phân tích

Dữ liệu thực nghiệm được xử lý bằng các công cụ toán thống kê và số trị tiên tiến:

  • Sử dụng mô hình TransPore 2D lập trình trên môi trường MATLAB để giải hệ phương trình đạo hàm riêng phi tuyến bảo toàn khối lượng và năng lượng trong môi trường xốp bằng phương pháp thể tích hữu hạn (Finite Volume Method) với lưới tam giác phi cấu trúc.
  • Phân tích hồi quy phi tuyến đa biến và phân tích phương sai (ANOVA) bằng phần mềm SPSS và MATLAB để thiết lập các hàm tương quan giữa $U$, $\Delta MC$, $EMC$ và mức độ nứt vỡ $N_i$.
  • Đánh giá sai số mô hình thông qua chỉ số sai số toàn phương trung bình gốc (RMSE) và sai số thời gian sấy: $$E_{dt} = \frac{|t_s - t_e|}{t_e} \times 100%$$ trong đó $t_s$ là thời gian dự đoán của mô hình và $t_e$ là thời gian sấy thực nghiệm.

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

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

  1. Bộ thông số cấu tạo vi mô và chuyển khối đặc trưng của Keo tai tượng 9 tuổi:
    • Khối lượng riêng cơ bản trung bình đạt $\rho_b = 0,486 - 0,532\text{ g/cm}^3$; Điểm bão hòa thớ gỗ xác định bằng phương pháp co rút tuyến tính đạt $FSP = 28,6%$.
    • Chiều dày vách tế bào sợi gỗ trung bình $T_{fw} = 1,85\text{ }\mu\text{m}$, đường kính ruột tế bào sợi gỗ $D_{fl} = 14,2\text{ }\mu\text{m}$, tỷ lệ $D_{fl} / T_{fw} \approx 7,67$ cho thấy gỗ có nguy cơ sụp đổ vách tế bào (mo móp) ở mức trung bình cao khi sấy nhiệt độ cao ở giai đoạn ẩm tự do.
    • Hệ số thấm khí dọc thớ ($K_L$) đạt giá trị trung bình $1,42 \times 10^{-12}\text{ m}^2$, cao hơn hệ số thấm xuyên tâm ($K_R \approx 2,15 \times 10^{-16}\text{ m}^2$) và tiếp tuyến ($K_T \approx 1,08 \times 10^{-16}\text{ m}^2$) từ 4 đến 5 bậc độ lớn ($p < 0,001$).
    • Hệ số khuếch tán ẩm ở trạng thái ổn định theo chiều dọc thớ $D_{bL} = 8,45 \times 10^{-10}\text{ m}^2\text{/s}$, trong khi theo chiều xuyên tâm $D_{bR} = 1,12 \times 10^{-10}\text{ m}^2\text{/s}$ và tiếp tuyến $D_{bT} = 0,95 \times 10^{-10}\text{ m}^2\text{/s}$.
  2. Quy luật hình thành dốc độ ẩm và động học khuyết tật dưới các dốc sấy khác nhau:
    • Ở mẻ sấy cứng ($U = 4,4 - 5,0$), dốc độ ẩm giữa tâm và mặt ván ($\Delta MC$) tăng vọt lên mức đỉnh điểm $38,5%$ chỉ sau 48 giờ sấy; ứng suất kéo bề mặt vượt giới hạn bền kéo ngang thớ làm xuất hiện nứt mặt sớm với tổng chiều dài nứt mặt khả kiến lên tới $1250\text{ mm/m}^2$, nứt đầu xuất hiện trên $85%$ số thanh mẫu tiếp tuyến.
    • Ở mẻ sấy mềm ($U = 2,0 - 2,5$), dốc độ ẩm $\Delta MC$ được duy trì ổn định dưới mức $15% - 18%$ trong suốt quá trình rút ẩm tự do; mức độ nứt mặt giảm $92%$, nứt đầu giảm $88%$, gỗ hầu như không bị mo móp.
  3. Xác lập hàm tương quan và các ngưỡng công nghệ tới hạn:
    • Xây dựng thành công các phương trình tương quan đa biến dạng hàm mũ và bậc hai giữa tổng mức độ nứt vỡ ($N_i$) với $\Delta MC$ và $EMC$ tại ba mốc kiểm soát độ ẩm chiến lược:
      • Tại mốc $MC = 50%$: $N_i = 0,042(\Delta MC)^2 - 0,31(\Delta MC) + 1,15$ ($R^2 = 0,892$).
      • Tại mốc $MC = 20%$: $N_i = 0,028(\Delta MC)^2 + 0,15(\Delta MC) - 0,45$ ($R^2 = 0,864$).
      • Tại mốc $MC = 12%$: $N_i = 0,019(\Delta MC)^2 + 0,08(\Delta MC) + 0,12$ ($R^2 = 0,851$).
    • Lượng hóa chính xác ngưỡng dốc độ ẩm tối đa $\Delta MC_{max} \le 16,2%$ và độ ẩm thăng bằng môi trường tối thiểu $EMC_{min} \ge 11,8%$ ở giai đoạn $MC > FSP$ để đảm bảo tỷ lệ nứt vỡ bằng $0%$.
  4. Kiểm chứng thành công mô hình TransPore 2D trên gỗ Keo tai tượng:
    • Kết quả mô phỏng số khớp nối chặt chẽ với dữ liệu thực nghiệm sấy trong lò quy chuẩn: đường cong giảm ẩm lý thuyết và thực nghiệm có độ tương đồng cao với sai số thời gian sấy $E_{dt} = 6,84%$ đối với ván xẻ xuyên tâm và $E_{dt} = 7,42%$ đối với ván xẻ tiếp tuyến ($RMSE < 1,8%$).
  5. Cơ chế vượt trội của công nghệ sấy NLMT đối với chất lượng gỗ Keo tai tượng:
    • Trong lò sấy NLMT, vào ban đêm khi bức xạ mặt trời bằng không, nhiệt độ lò giảm tự nhiên và độ ẩm tương đối ($RH$) tăng lên, tạo ra chu kỳ "hồi ẩm bề mặt" (conditioning cycle) tự động. Quá trình này giúp ẩm từ tâm tiếp tục khuếch tán ra ngoài làm phẳng đường cong dốc ẩm, triệt tiêu sự tích tụ ứng suất dư nguy hiểm, ngăn chặn hoàn toàn hiện tượng nứt ngầm và chai cứng bề mặt mà không cần cấp hơi nhiệt cưỡng bức tốn kém.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp cơ sở dữ liệu khoa học đầy đủ nhất từ trước đến nay về các hằng số chuyển khối vi mô của Keo tai tượng nhiệt đới, đóng góp vào thư viện vật liệu cho các mô hình truyền nhiệt - truyền khối quốc tế.
  • Về mặt phương pháp luận: Khẳng định tính khả thi của việc kết hợp giữa phân tích vi cấu trúc giải phẫu hiển vi với phương pháp cắt lát vi sai và mô hình hóa số trị hai chiều, tạo hình mẫu phương pháp luận cho nghiên cứu sấy các loài gỗ rừng trồng nhiệt đới khác.
  • Về mặt thực tiễn sản xuất: Cung cấp bảng chế độ sấy NLMT chuẩn hóa cho ván xẻ Keo tai tượng dày 25 mm với 5 giai đoạn điều khiển dốc sấy linh hoạt ($U = 2,0 - 3,2$), rút ngắn thời gian sấy xuống còn 14 - 18 ngày (nhanh hơn phơi tự nhiên 2,5 - 3 lần), đạt độ ẩm cuối cùng $MC_c = 10 \pm 1%$, tỷ lệ khuyết tật loại A đạt trên $92%$.
  • Về mặt chính sách và môi trường: Đẩy mạnh lộ trình chuyển đổi xanh trong ngành chế biến gỗ Việt Nam, giảm phát thải khí nhà kính $CO_2$ nhờ thay thế nồi hơi đốt nhiên liệu hóa thạch/củi trấu bằng năng lượng mặt trời sạch, phù hợp với cam kết Net Zero vào năm 2050 của Chính phủ.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận 4 hạn chế khoa học:

  1. Giả định về hệ số đẳng nhiệt nhả ẩm: Mô hình số học TransPore 2D trong luận án giả định hệ số đẳng nhiệt nhả ẩm bằng 1 (bỏ qua khoảng trễ hấp phụ $W$ giữa đường cong hút ẩm và nhả ẩm do thiếu thiết bị đo sorption isotherm động học liên tục ở các dải nhiệt độ cao).
  2. Không gian mô phỏng 2 chiều: Mô hình mới dừng lại ở không gian 2D trên mặt cắt ngang thanh ván ($x, y$), chưa tích hợp dòng chuyển dịch ẩm dọc thớ đa chiều 3D dọc theo chiều dài thanh gỗ ($z$), vốn có ý nghĩa tại khu vực đầu mút ván xẻ.
  3. Phạm vi lập địa và tuổi gỗ: Mẫu thí nghiệm thu thập tập trung tại một vùng lập địa (Vĩnh Phúc) ở độ tuổi 9 năm; chưa khảo sát biến động chuyển khối theo các cấp tuổi khác nhau (5, 7, 12 năm) hoặc giữa các lập địa có sự khác biệt lớn về khí hậu và thổ nhưỡng (Tây Nguyên, Nam Bộ).
  4. Sự phụ thuộc thời tiết của lò sấy NLMT: Lò sấy NLMT thí nghiệm chưa tích hợp bộ điều khiển thông minh tự động dự báo thời tiết (AI-based predictive control) để bù nhiệt tối ưu trong các đợt mưa kéo dài.

Chương trình nghiên cứu tương lai (Future Research Agenda) cần tập trung vào:

  • Mở rộng mô hình TransPore sang không gian 3 chiều (3D finite volume) kết hợp cơ cấu liên hợp phân tích ứng suất hữu hạn (FEA) để dự báo hình thái nứt vỡ ba chiều theo thời gian thực.
  • Xác định đường đẳng nhiệt hấp phụ - giải hấp phụ ($sorption\text{ }isotherms$) thực nghiệm ở các cấp nhiệt độ từ $30^\circ\text{C}$ đến $80^\circ\text{C}$ cho gỗ Keo tai tượng.
  • Nghiên cứu tích hợp hệ thống sấy NLMT lai ghép (Hybrid Solar Kiln) sử dụng bơm nhiệt (Heat Pump) hoặc phụ trợ sinh khối điều khiển tự động qua IoT.
  • Khảo sát biến thiên tính chất chuyển khối và dốc ẩm của các dòng vô tính Keo lai (Acacia mangium $\times$ Acacia auriculiformis) và Keo lá tràm.

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

Luận án tạo ra những tác động lan tỏa sâu rộng:

  • Tác động học thuật: Mở ra hướng nghiên cứu liên ngành chuyên sâu giữa giải phẫu học thực vật, vật lý nhiệt - ẩm và mô phỏng số trị cơ cấu xốp tại Việt Nam; cung cấp tài liệu tham khảo chuẩn mực cho các chương trình đào tạo sau đại học ngành Kỹ thuật chế biến Lâm sản.
  • Tác động công nghiệp: Kết quả nghiên cứu đã được chuyển giao và ứng dụng trực tiếp tại các cơ sở chế biến gỗ quy mô vừa và nhỏ tại Vĩnh Phúc, Phú Thọ và Đồng Nai. Việc chuyển đổi sang lò sấy NLMT mái vòm giúp doanh nghiệp tiết kiệm từ $180.000$ đến $250.000\text{ VNĐ/m}^3$ chi phí năng lượng sấy so với lò hơi truyền thống, giảm tỷ lệ hao hụt gỗ do nứt vỡ từ $12%$ xuống dưới $3%$.
  • Tác động môi trường và xã hội: Góp phần hiện thực hóa chiến lược sản xuất sạch hơn, cải thiện môi trường làm việc cho công nhân khi không phải vận hành lò hơi áp lực cao và giảm thiểu khói bụi phát thải ra cộng đồng dân cư lân cận.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận nghiên cứu thực nghiệm chuẩn mực quốc tế kết hợp công cụ mô phỏng số trị hiện đại trong khoa học gỗ.
  • Các nhà khoa học và giảng viên ngành Lâm sản: Sử dụng cơ sở dữ liệu về hằng số vật lý và chuyển khối vi mô của gỗ Keo tai tượng để giảng dạy và phát triển các mô hình lý thuyết chuyên sâu.
  • Kỹ sư công nghệ và doanh nghiệp chế biến gỗ: Ứng dụng trực tiếp bảng chế độ sấy và quy trình công nghệ sấy NLMT để tối ưu hóa chi phí sản xuất, nâng cao chất lượng gỗ xẻ xuất khẩu.
  • Cơ quan quản lý và hoạch định chính sách (Bộ Nông nghiệp và PTNT, Bộ KH&CN): Có thêm luận cứ khoa học để xây dựng tiêu chuẩn kỹ thuật quốc gia về sấy gỗ rừng trồng và ban hành chính sách khuyến khích ứng dụng năng lượng tái tạo trong chế biến nông 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ì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc lượng hóa và tích hợp thành công lý thuyết dốc ẩm vi phân với lý thuyết đảo chiều ứng suất của McMillen (1958) trên đối tượng gỗ lá rộng nhiệt đới mọc nhanh. Luận án đã làm sáng tỏ cơ chế vi mô: sự co rút bất đẳng hướng và dốc ẩm vượt ngưỡng ($\Delta MC > 16,2%$) là nguồn gốc trực tiếp kích hoạt ứng suất kéo vượt quá ứng suất phá hủy ngang thớ của vách tế bào sợi gỗ mỏng ($T_{fw} = 1,85\text{ }\mu\text{m}$), đồng thời chứng minh quá trình sấy gián đoạn NLMT tạo ra sự hồi phục biến dạng từ biến (viscoelastic creep recovery), làm triệt tiêu nội ứng suất dư mà không cần giai đoạn xử lý nhiệt ẩm cuối mẻ sấy.

2. Tính đổi mới về phương pháp luận so với ít nhất hai nghiên cứu quốc tế trước đây?

So với nghiên cứu của Haque (2002) vốn sử dụng mô hình khuếch tán 1D đơn giản bỏ qua pha thấm lỏng, luận án đã ứng dụng mô hình số TransPore 2D đầy đủ giải bằng thể tích hữu hạn, tích hợp đồng thời dòng thấm Darcy và dòng khuếch tán Fick. So với nghiên cứu của Yuniarti (2015) và Phonetip (2018) tại Úc vốn chỉ đo ứng suất dư sau sấy bằng phương pháp cắt ngón tay (prong test) theo tiêu chuẩn AS/NZS 4787:2001, luận án đã đổi mới bằng việc theo dõi động học dốc độ ẩm đa điểm liên tục từ lúc gỗ tươi ($MC > 100%$) bằng kỹ thuật cắt lát vi phân chính xác kết hợp đo thông số chuyển khối trên các hệ thống chuyên dụng POROLUX 1000 và PVC-CHA.

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ờ và thú vị nhất là ở mẻ sấy cứng ($U = 4,4 - 5,0$), mặc dù các vết nứt mặt phát triển rất mạnh ở giai đoạn đầu ($MC$ từ $100%$ xuống $50%$), nhưng khi độ ẩm trung bình xuống dưới $15%$, các vết nứt mặt này có xu hướng tự khép miệng lại (tự đóng kín) do lớp lõi bên trong bắt đầu co rút mạnh và kéo lớp bề mặt vào trạng thái chịu nén. Tuy nhiên, vết nứt khép miệng này tạo ra "khuyết tật ẩn", sẽ tái nứt vỡ nghiêm trọng khi xẻ gia công hoặc thay đổi độ ẩm môi trường sử dụng. Điều này khẳng định không thể đánh giá chất lượng sấy thuần túy bằng mắt thường sau khi sấy xong mà phải kiểm soát dốc ẩm trong suốt quá trình sấy.

4. Luận án có cung cấp quy trình nhân bản thí nghiệm (Replication Protocol) hoàn chỉnh không?

Có. Luận án cung cấp đầy đủ và chi tiết giao thức nhân bản thực nghiệm bao gồm: tiêu chuẩn chọn cây mẫu, kích thước và vị trí lấy mẫu lõi khoan đo thấm/khuếch tán, bản vẽ chi tiết thiết bị cắt lát mẫu vi phân, toàn bộ thông số cài đặt điều kiện biên (nhiệt độ khô $T_k$, nhiệt độ ướt $T_u$, vận tốc gió $v$, độ dốc sấy $U$, độ ẩm thăng bằng $EMC$) cho phần mềm MATLAB chạy mô hình TransPore 2D, cùng bảng hướng dẫn quy đổi thông số điều khiển quạt gió và cửa thoát ẩm theo từng mốc độ ẩm trong lò sấy NLMT công nghiệp.

5. Chương trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Chương trình 10 năm tới tập trung vào 3 trọng tâm:

  1. Xây dựng bản đồ chuyển khối số hóa (Digital Mass Transfer Atlas) cho 10 loài gỗ rừng trồng chủ lực của Việt Nam.
  2. Phát triển phần mềm thương mại mô phỏng và điều khiển thông minh lò sấy gỗ thời gian thực ứng dụng trí tuệ nhân tạo (AI-Kiln Control) tích hợp dữ liệu cảm biến IoT dốc ẩm.
  3. Hoàn thiện công nghệ lò sấy NLMT module hóa tích hợp pin lưu trữ năng lượng và bơm nhiệt công nghiệp phục vụ chế biến gỗ sâu tại các vùng nguyên liệu quy mô tập trung.

Kết luận

  1. Xác định hoàn chỉnh và chính xác bộ cơ sở dữ liệu giải phẫu hiển vi và hằng số chuyển khối của gỗ Keo tai tượng 9 năm tuổi: $\rho_b = 0,486\text{ g/cm}^3$, $FSP = 28,6%$, $K_L = 1,42 \times 10^{-12}\text{ m}^2$, $K_R = 2,15 \times 10^{-16}\text{ m}^2$, $K_T = 1,08 \times 10^{-16}\text{ m}^2$, $D_{bL} = 8,45 \times 10^{-10}\text{ m}^2\text{/s}$, $D_{bR} = 1,12 \times 10^{-10}\text{ m}^2\text{/s}$, $D_{bT} = 0,95 \times 10^{-10}\text{ m}^2\text{/s}$.
  2. Thiết lập định lượng bản chất mối quan hệ giữa dốc độ ẩm theo chiều dày và sự hình thành khuyết tật sấy; chứng minh dốc ẩm vượt ngưỡng $\Delta MC > 16,2%$ ở giai đoạn $MC > FSP$ là nguyên nhân cốt lõi gây nứt mặt và nứt đầu ván xẻ.
  3. Xây dựng thành công hệ phương trình hồi quy phi tuyến xác thực tương quan giữa dốc độ ẩm, độ ẩm thăng bằng môi trường với mức độ nứt vỡ tại các mốc kiểm soát độ ẩm chiến lược ($MC = 50%$, $20%$, $12%$).
  4. Ứng dụng và kiểm chứng thành công mô hình toán học TransPore 2D giải bằng phương pháp thể tích hữu hạn trong mô phỏng sấy quy chuẩn gỗ Keo tai tượng với sai số thời gian sấy nhỏ ($E_{dt} < 7,5%$).
  5. Khám phá và luận giải khoa học cơ chế tự cân bằng ẩm và thư giãn ứng suất ban đêm của công nghệ sấy NLMT, làm sáng tỏ ưu thế vượt trội của sấy gián đoạn đối với gỗ lá rộng nhiệt đới.
  6. Đề xuất quy trình công nghệ sấy gỗ Keo tai tượng dày 25 mm bằng lò sấy NLMT kiểu mái vòm hấp thụ nhiệt cao, cho phép rút ngắn thời gian sấy xuống 14 - 18 ngày, nâng tỷ lệ gỗ đạt chất lượng cao lên trên $90%$ và giảm $30% - 45%$ chi phí năng lượng so với sấy hơi nước, mang lại hiệu quả kinh tế - kỹ thuật và môi trường vượt bậc cho ngành chế biến gỗ Việt Nam.