Tổng quan về luận án

Nghiên cứu của tác giả Hà Tiến Mạnh tại Viện Khoa học Lâm nghiệp Việt Nam, dưới sự hướng dẫn khoa học của GS. Phạm Văn Chương và TS. Bùi Duy Ngọc, tập trung giải quyết bài toán cốt lõi trong công nghệ chế biến gỗ: "Ảnh hưởng của chênh lệch độ ẩm trong gỗ xẻ Keo tai tượng (Acacia mangium Willd.) đến khuyết tật khi sấy, ứng dụng cho lò sấy năng lượng mặt trời". Trong bối cảnh diện tích rừng trồng tại Việt Nam tăng trưởng mạnh mẽ hơn ba thập kỷ qua, Keo tai tượng (Acacia mangium Willd.) đã trở thành loài cây trồng chủ lực cung cấp nguyên liệu cho ngành công nghiệp chế biến và xuất khẩu gỗ. Tuy nhiên, do thuộc nhóm gỗ III theo phân loại của Hồ Xuân Các (1999) – 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 như nứt mặt, nứt đầu, mo móp (cell collapse) và cong vênh – quá trình sấy công nghiệp đối mặt với thách thức lớn về tỷ lệ phế phẩm và chi phí năng lượng.

Tính tiên phong của luận án thể hiện ở việc lần đầu tiên tại Việt Nam thiết lập mối quan hệ định lượng giữa cấu tạo hiển vi, các thông số chuyển khối cơ bản (hệ số thấm $K$, hệ số khuếch tán $D_b$, điểm bão hòa thớ gỗ $\text{FSP}$) với mức chênh lệch độ ẩm vi phân theo chiều dày (MC gradient) và mức độ khuyết tật, từ đó xây dựng mô hình toán học mô phỏng 2D và tối ưu hóa quy trình sấy năng lượng mặt trời (NLMT).

Khoảng trống nghiên cứu (research gap) mà luận án giải quyết xuất phát từ thực tế các nghiên cứu sấy gỗ trong nước trước đây (Nguyễn Xuân Hiên, 2006; Đỗ Văn Bản, 2012) chủ yếu sử dụng máy đo điện trở để ước lượng chênh lệch ẩm giữa bề mặt và tâm. Phương pháp này bộc lộ sai số rất lớn khi độ ẩm gỗ nằm trên điểm bão hòa thớ gỗ ($\text{MC} > \text{FSP}$). Đồng thời, các nghiên cứu sấy NLMT tại Việt Nam (Trần Đức Sinh, 2009; Hứa Thị Huần và Hồ Xuân Các, 2012; Bùi Duy Ngọc, 2021) đa phần tập trung vào thiết kế cơ khí và thiết bị mà thiếu cơ sở dữ liệu động học chuyển khối nội tại để giải thích bản chất phát sinh khuyết tật.

Luận án thiết lập 4 câu hỏi nghiên cứu (RQ) và 2 giả thuyết khoa học ($H$) chính:

  • RQ1: Các đặc điểm cấu tạo giải phẫu và thông số chuyển khối ($K, D_b, \rho_b, \text{FSP}$) của gỗ Keo tai tượng chi phối động học vận chuyển ẩm như thế nào?
  • RQ2: Quy luật biến thiên của mức chênh lệch độ ẩm theo chiều dày (MC gradient) tác động ra sao đến sự hình thành các khuyết tật nứt vỡ, mo móp và biến dạng?
  • RQ3: Mô hình số TransPore hai chiều có khả năng mô phỏng chính xác trường độ ẩm và nhiệt độ trong lò sấy quy chuẩn hay không?
  • RQ4: Cơ chế sấy gián đoạn bằng NLMT điều tiết dốc sấy $U$ và độ ẩm thăng bằng ($\text{EMC}$) như thế nào để triệt tiêu khuyết tật và tối ưu chi phí?
  • $H_1$: Sự gia tăng dốc ẩm vượt ngưỡng tới hạn sẽ tạo ra chênh lệch co rút dị hướng, làm nội ứng suất kéo bề mặt vượt quá cường độ chịu kéo ngang thớ, trực tiếp kích hoạt hiện tượng nứt vỡ.
  • $H_2$: Pha hồi phục ẩm ban đêm trong chu trình sấy NLMT gián đoạn giúp tự cân bằng dốc ẩm, giải tỏa ứng suất dư và làm giảm khuyết tật sấy so với phương pháp sấy cưỡng bức liên tục.

Khung lý thuyết của nghiên cứu tích hợp cơ học vật liệu xốp dị hướng, định luật thấm Darcy điều chỉnh, định luật khuếch tán Fick và lý thuyết động học nội ứng suất sấy của McMillen (1958). Phạm vi nghiên cứu thực nghiệm được thực hiện trên gỗ Keo tai tượng 9 năm tuổi khai thác tại Phúc Yên, Vĩnh Phúc; bao gồm 200 thanh ($1000 \times 80 \times 27\text{ mm}$) để kiểm tra cấu tạo - vật lý, 60 thanh ($1322 \times 80 \times 25\text{ mm}$) có độ ẩm ban đầu $\text{MC} > 140%$ cho các mẻ sấy quy chuẩn đo dốc ẩm bằng kỹ thuật cắt lát, 200 thanh ($1000 \times 80 \times 25\text{ mm}$) cho mẻ sấy thí nghiệm NLMT, và kiểm chứng sản xuất tại lò sấy NLMT công suất $30\text{ m}^3/\text{mẻ}$ so sánh với lò hơi nước $35\text{ m}^3/\text{mẻ}$.

Literature Review và Positioning

Nghiên cứu về sấy gỗ trên thế giới đã phát triển qua nhiều giai đoạn với các trường phái lý thuyết chuyên sâu:

                  ┌─────────────────────────────────────────────────────────┐
                  │                 LÝ THUYẾT VẬN CHUYỂN ẨM                 │
                  │   Stamm (1967), Siau (1984), Skaar (1988), Walker (2006) │
                  └────────────────────────────┬────────────────────────────┘
                                               │
                       ┌───────────────────────┴───────────────────────┐
                       ▼                                               ▼
     ┌───────────────────────────────────┐           ┌───────────────────────────────────┐
     │         CƠ CHẾ THẤM DÒNG          │           │       CƠ CHẾ KHUẾCH TÁN ẨM        │
     │      (MC > FSP, Áp suất ΔP)       │           │      (MC < FSP, Nồng độ ΔC)       │
     │  Darcy; Perré (1987); Rousset     │           │  Fick I & II; Siau (1984);        │
     │  (2004) ALU-CHA / Porolux 1000    │           │  Perré & Agoua (2002) PVC-CHA     │
     └─────────────────┬─────────────────┘           └─────────────────┬─────────────────┘
                       │                                               │
                       └───────────────────────┬───────────────────────┘
                                               ▼
                  ┌─────────────────────────────────────────────────────────┐
                  │              MÔ HÌNH HÓA VÀ NỘI ỨNG SUẤT                │
                  │   McMillen (1958) Stress Inversion; Perré & Turner     │
                  │   (1999) TransPore 2D; Carr (2013); Redman (2017) FEA   │
                  └────────────────────────────┬────────────────────────────┘
                                               │
                                               ▼
                  ┌─────────────────────────────────────────────────────────┐
                  │              VỊ TRÍ ĐÓNG GÓP CỦA LUẬN ÁN                │
                  │   Định lượng MC Gradient bằng cắt lát vi phân           │
                  │   (R²=0.8795 vs Điện trở) + TransPore 2D Acacia mangium │
                  │   + Tối ưu hóa dốc sấy NLMT gián đoạn (U=3.2)           │
                  └─────────────────────────────────────────────────────────┘

Dòng nghiên cứu về cơ chế vận chuyển ẩm và chuyển khối: Stamm (1967), Siau (1984), Skaar (1988) và Walker (2006) đã chuẩn hóa lý thuyết vận chuyển ẩm trong vật liệu gỗ lá rộng. Nước tự do trong lòng tế bào di chuyển chủ yếu nhờ cơ chế thấm dưới tác động của chênh lệch áp suất khi $\text{MC} > \text{FSP}$. Ngược lại, nước liên kết trong vách tế bào và hơi nước trong ruột tế bào di chuyển qua cơ chế khuếch tán dưới tác động của chênh lệch nồng độ khi $\text{MC} < \text{FSP}$. Perré (1987) và Rousset et al. (2004) đã phát triển hệ thống ALU-CHA xác định hệ số thấm khí và lỏng, trong khi Perré & Agoua (2002) thiết lập hệ thống PVC-CHA để xác định hệ số khuếch tán ở trạng thái ổn định.

Dòng nghiên cứu về động lực học nội ứng suất và mô hình hóa: McMillen (1958) đã đặt nền móng khi chứng minh mối quan hệ giữa dốc ẩm và sự đảo chiều ứng suất sấy (stress inversion) ở gỗ Sồi đỏ (Quercus spp.). Giai đoạn đầu, bề mặt gỗ nhanh chóng xuống dưới $\text{FSP}$ và co rút, chịu ứng suất kéo trong khi tâm chịu nén; khi tâm gỗ đạt $\text{FSP}$, ứng suất đảo chiều, bề mặt chuyển sang chịu nén và tâm chịu kéo, dẫn đến nứt ngầm nếu vượt giới hạn bền. Về mô hình số, Perré & Turner (1999) đã công bố mô hình TransPore 2D mô phỏng dòng chuyển ẩm và nhiệt trong vật liệu xốp, sau đó được Carr et al. (2013) nâng cấp thuật toán bước thời gian phi tuyến và Redman (2017) lập trình trên Matlab cho các loài gỗ lá rộng của Úc.

Tranh luận khoa học cốt lõi: Tồn tại hai luồng quan điểm đối lập trong mô hình hóa sấy gỗ. Trường phái truyền thống (Salin, 2010) ủng hộ mô hình thuần khuếch tán (diffusion-based models) vì sự đơn giản trong tính toán. Ngược lại, trường phái cơ học môi trường xốp (Perré, 2007; Redman, 2017) khẳng định mô hình thuần khuếch tán chỉ chính xác dưới $\text{FSP}$ và bắt buộc phải kết hợp phương trình thấm Darcy để mô tả chính xác quá trình thoát nước tự do ở độ ẩm cao. Về công nghệ, các nghiên cứu sấy cưỡng bức (Hồ Thu Thủy, 2005; Nguyễn Cảnh Mão, 2015) ưu tiên tốc độ sấy nhưng tiêu tốn năng lượng cao, trong khi các công trình sấy NLMT (Plumptre, 1979; Read et al., 1974; Phonetip, 2018) chứng minh tính ưu việt về chất lượng nhờ cơ chế hồi phục ứng suất tự nhiên.

Định vị của luận án: Luận án đã kết nối khoảng trống giữa vật lý gỗ cơ bản và công nghệ sấy ứng dụng tại Việt Nam. Bằng việc so sánh đối sánh với hai nghiên cứu quốc tế điển hình – công trình của Jusoh et al. (2014) tại Úc về cấu tạo giải phẫu và độ co rút của Acacia mangium, và nghiên cứu của Phonetip (2018) tại Đại học Melbourne về mô phỏng ứng suất sấy NLMT trên gỗ Bạch đàn – luận án đã mở rộng mô hình TransPore 2D cho gỗ Keo tai tượng rừng trồng nhiệt đới và chuẩn hóa công nghệ sấy NLMT với dốc sấy $U = 3,2$.

Đó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 và kiểm chứng thực nghiệm hệ thống lý thuyết sấy gỗ quốc tế trên đối tượng gỗ rừng trồng mọc nhanh nhiệt đới thông qua các đóng góp cụ thể:

  1. Kiểm chứng và mở rộng lý thuyết ứng suất sấy của McMillen (1958): Thay thế việc ước lượng ứng suất định tính bằng các hàm tương quan định lượng giữa mức chênh lệch độ ẩm vi phân theo chiều dày ($\Delta\text{MC}$), độ ẩm thăng bằng ($\text{EMC}$) và tỷ lệ khuyết tật nứt vỡ ($N_i, NM_i, ND_i$).
  2. Bổ sung cơ sở dữ liệu chuyển khối vào mô hình TransPore (Perré & Turner, 1999): Xác định toàn diện bộ hằng số vật lý thực nghiệm cho gỗ Keo tai tượng 9 tuổi gồm hệ số thấm lỏng/khí ($K_L, K_R, K_T$), hệ số khuếch tán hơi nước ($D_{bL}, D_{bR}, D_{bT}$), độ rỗng ruột tế bào ($\phi_t, \phi_v, \phi_f$) và điểm $\text{FSP}$ thực tế.
  3. Thiết lập mệnh đề chuyển pha và sụp đổ tế bào: Chứng minh tỷ lệ giữa đường kính ruột tế bào và bề dày vách sợi gỗ ($D_{fl}/T_{fw}$) là chỉ số hình thái học quyết định mức độ nhạy cảm của gỗ đối với khuyết tật mo móp trong giai đoạn sấy ban đầu trên $\text{FSP}$.

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

Khung phân tích của luận án tích hợp đa chiều giữa 3 lý thuyết nền tảng:

  • Lý thuyết dòng chảy mao dẫn Darcy điều chỉnh cho môi trường dị hướng ($K_L \gg K_R, K_T$).
  • Lý thuyết khuếch tán phân tử Fick I & II cho pha nước liên kết trong vách tế bào.
  • Lý thuyết biến dạng đàn - nhớt dẻo (viscoelasticity) giải thích sự tích tụ và giải tỏa nội ứng suất.
┌─────────────────────────────────────────────────────────────────────────────────────────────┐
│                            KHUNG PHÂN TÍCH CHUYỂN KHỐI VÀ KHUYẾT TẬT                        │
├───────────────────────────────┬──────────────────────────────┬──────────────────────────────┤
│    ĐẶC TÍNH CẤU TRÚC VÀ CƠ LÝ  │    ĐỘNG HỌC VẬN CHUYỂN ẨM    │    ỨNG SUẤT VÀ KHUYẾT TẬT    │
├───────────────────────────────┼──────────────────────────────┼──────────────────────────────┤
│ • Giải phẫu: OM & SEM         │ • Thấm khí/lỏng: Porolux 1000│ • Dốc ẩm vi phân: Cắt lát    │
│ • Độ rỗng: ImageJ (φt, φv, φf)│ • Khuếch tán: PVC-CHA        │ • Khuyết tật: Nứt mặt, nứt   │
│ • Dfl/Tfw ratio: Mo móp       │ • Mô hình hóa: TransPore 2D  │   đầu, mo móp, cong vênh     │
│ • Co rút dị hướng, FSP, ρb    │ • Điều tiết môi trường: Tk,  │ • Tối ưu hóa: Ngưỡng dốc sấy │
│                               │   Tư, RH, EMC                │   U, chế độ sấy NLMT         │
└───────────────────────────────┴──────────────────────────────┴──────────────────────────────┘

Điều kiện biên của khung phân tích được xác định rõ: gỗ xẻ Keo tai tượng 9 năm tuổi, chiều dày danh nghĩa $25 \pm 1\text{ mm}$, nhiệt độ sấy môi trường từ $30^\circ\text{C}$ đến $65^\circ\text{C}$, giả định hệ số đẳng nhiệt trao đổi ẩm bằng 1.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng (positivism) kết hợp chặt chẽ giữa phân tích giải phẫu vi mô, đo đạc vật lý thực nghiệm độ chính xác cao và mô phỏng số học đa cấp độ. Thiết kế nghiên cứu phân tầng gồm ba cấp độ:

  • Cấp độ vi mô (Micro-scale): Xác định đặc điểm hình thái tế bào, vách thứ cấp $S_2$, vi sợi cellulose, cấu trúc lỗ thông ngang có vành, màng không nút trên kính hiển vi quang học (OM) và kính hiển vi điện tử quét (SEM).
  • Cấp độ trung mô (Meso-scale): Xác định các thông số chuyển khối động lực học ($K, D_b$) và thông số co rút tế bào qua hệ thống buồng đo chuyên dụng.
  • Cấp độ vĩ mô (Macro-scale): Thử nghiệm sấy quy chuẩn, sấy NLMT và ứng dụng sản xuất quy mô công nghiệp.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu tuân thủ các tiêu chuẩn quốc tế nghiêm ngặt:

  • Chuẩn bị mẫu và đo đạc chuyển khối: Mẫu xác định hệ số thấm khí và lỏng được khoan rút lõi theo ba hướng giải phẫu: dọc thớ (L), xuyên tâm (R), tiếp tuyến (T) và đo trên thiết bị tự động Porolux™ 1000 Porometer (nguyên lý ALU-CHA). Hệ số khuếch tán ẩm ($D_b$) được đo bằng phương pháp buồng đo thăng bằng ổn định trên hệ thống PVC-CHA vaporimeter đặt trong tủ khí hậu kiểm soát nhiệt độ và độ ẩm tuyệt đối.
  • Phương pháp cắt lát vi phân xác định dốc ẩm: Khắc phục hoàn toàn nhược điểm của ẩm kế điện trở. Mẫu thử ($1322 \times 80 \times 25\text{ mm}$) được bọc kín và cắt lát theo chiều dày tấm gỗ bằng thiết bị cắt lát chuyên dụng ngay trong quá trình sấy, sau đó cân sấy khô kiệt theo tiêu chuẩn ISO 4470:1981 và TCVN 13701-1:2023. Tác giả chỉ ra:

    "Nghiên cứu thăm dò cho thấy hệ số tương quan về kết quả đo của phương pháp cắt lát so với phương pháp điện trở là rất cao ($R^2 = 0,8795$) khi xác định mức chênh lệch độ ẩm."

  • Đo lường hình thái giải phẫu: Sử dụng phần mềm ImageJ phân tích ảnh chụp SEM (tại Đại học Kyoto, Nhật Bản) để đo diện tích ruột mạch gỗ ($A_j$), diện tích sợi gỗ ($A_i$), từ đó tính toán độ rỗng sợi gỗ ($\phi_f$), độ rỗng mạch ($\phi_v$) và tổng độ rỗng ($\phi_t$).
┌─────────────────────────────────────────────────────────────────────────────────────────────┐
│                          SƠ ĐỒ QUY TRÌNH THỰC NGHIỆM VÀ ĐỐI CHỨNG                           │
├─────────────────────────────────────────────────────────────────────────────────────────────┤
│  Gỗ Keo tai tượng 9 tuổi (Vĩnh Phúc) ──► Cưa xẻ thanh tiêu chuẩn (1322 x 80 x 25 mm)       │
│                                      │                                                      │
│  ┌───────────────────────────────────┴───────────────────────────────────┐                  │
│  ▼                                                                       ▼                  │
│  [PHÒNG THÍ NGHIỆM DAF - ÚC & KYOTO]                     [THỰC NGHIỆM SẤY QUY CHUẨN]        │
│  • SEM / OM ──► ImageJ ──► Độ rỗng φt                    • Mẻ sấy cứng (U = 4,4 - 5,0)      │
│  • Porolux 1000 ──► Hệ số thấm K (L, R, T)               • Mẻ sấy mềm (U = 2,0 - 2,5)       │
│  • PVC-CHA ──► Khuếch tán Db (L, R, T)                   • Cắt lát vi phân xác định ΔMC     │
│  • Co rút màng mỏng ──► FSP, ρb                          • Đo khuyết tật: NM, ND, Mo móp    │
│  └───────────────────────────────────┬───────────────────────────────────┘                  │
│                                      │                                                      │
│                                      ▼                                                      │
│  [MÔ HÌNH TOÁN TRANSPORE 2D (MATLAB)] ◄─── Nhập dữ liệu chuyển khối                         │
│  • Mô phỏng phân bố MC(x,y,t) và T(x,y,t)                                                   │
│  • Đối chứng mô hình vs thực nghiệm (Sai số thời gian Edt, RMSE)                             │
│                                      │                                                      │
│                                      ▼                                                      │
│  [TỐI ƯU HÓA VÀ KHẢO NGHIỆM LÒ SẤY NĂNG LƯỢNG MẶT TRỜI]                                     │
│  • Thử nghiệm các dốc sấy U (U = 3,2 tối ưu)                                                │
│  • Kiểm chứng sản xuất 30 m³ NLMT vs 35 m³ lò hơi nước ──► Tiết kiệm chi phí 31,3%           │
└─────────────────────────────────────────────────────────────────────────────────────────────┘

Data và phân tích

Dữ liệu thực nghiệm được xử lý bằng phân tích phương sai (ANOVA) và hồi quy phi tuyến trên phần mềm SPSS và Matlab:

  • Dữ liệu cấu tạo và chuyển khối cung cấp trực tiếp vào hệ thống phương trình vi phân đạo hàm riêng của mô hình TransPore 2D (Perré & Turner, 1999; Redman, 2017). Mô hình giải phương trình bảo toàn khối lượng (nước lỏng, hơi nước, không khí) và bảo toàn năng lượng (entanpi) bằng phương pháp thể tích hữu hạn (Finite Volume Method - FVM) trên lưới phần tử tam giác trong miền chữ nhật hai chiều.
  • Thuật toán bước nhảy thời gian phi tuyến theo hàm mũ (exponential time-stepping) được áp dụng để kiểm soát sai số hội tụ, đảm bảo chỉ số sai số toàn phương trung bình gốc ($\text{RMSE}$) và sai số dự đoán thời gian sấy ($E_{dt}$) ở mức tối thiểu.

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

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

  1. Tính dị hướng cực đoan trong đặc tính chuyển khối: Kết quả đo đạc trên hệ thống Porolux 1000 và PVC-CHA cho thấy sự chênh lệch vượt trội giữa các chiều giải phẫu. Hệ số thấm dọc thớ ($K_L$) cao hơn hàng trăm đến hàng nghìn lần so với hệ số thấm xuyên tâm ($K_R$) và tiếp tuyến ($K_T$). Luận án nhấn mạnh trích dẫn lý thuyết:

    "Tỷ lệ lớn hơn của hệ số khuếch tán dọc thớ so với ngang thớ từ 100 lần ở độ ẩm 5 % đến từ 2 đến 4 lần ở độ ẩm từ 25 % đến 30 %." Đồng thời, cấu tạo sợi gỗ Keo tai tượng có tỷ lệ đường kính ruột trên bề dày vách ($D_{fl}/T_{fw}$) cao, vách sợi mỏng ($1,72 - 2,03\text{ }\mu\text{m}$), giải thích nguyên nhân gỗ bị mo móp mạnh khi sấy ở nhiệt độ cao ở giai đoạn $\text{MC} > \text{FSP}$.

  2. Quy luật động học dốc ẩm và ngưỡng khuyết tật nứt vỡ: So sánh mẻ sấy cứng ($U = 4,4 - 5,0$) và mẻ sấy mềm ($U = 2,0 - 2,5$) đã chứng minh: dốc ẩm đạt cực đại ở thời điểm độ ẩm trung bình đạt khoảng $50%$. Ở mẻ sấy cứng, tốc độ thoát ẩm bề mặt quá nhanh tạo ra dốc ẩm vượt ngưỡng chịu đựng, khiến chiều dài vết nứt mặt và nứt đầu tăng vọt. Thiết lập thành công hệ phương trình tương quan giữa mức độ nứt vỡ ($N_i$), mức chênh lệch độ ẩm ($\Delta\text{MC}$) và $\text{EMC}$ tại các mốc độ ẩm $50%$, $20%$ và $12%$.

    Xác định chính xác các ngưỡng dốc ẩm tối đa và $\text{EMC}$ tối thiểu để kiểm soát khuyết tật:

    • Để đạt mức độ nứt vỡ $0%$: Dốc ẩm $\Delta\text{MC}$ tại thời điểm $\text{MC} = 50%$ không được vượt quá $28,5%$, và $\text{EMC}$ môi trường phải duy trì $\ge 13,8%$.
    • Khi $\text{MC} = 20%$, ngưỡng $\Delta\text{MC}$ cho phép là $\le 11,2%$ với $\text{EMC} \ge 9,5%$.

  3. Kiểm chứng độ chính xác của mô hình số TransPore 2D: Mô phỏng trường độ ẩm và nhiệt độ trên mặt cắt ngang ván xuyên tâm và ván tiếp tuyến cho thấy sự tương thích rất cao với số liệu thực nghiệm. Đường cong giảm ẩm lý thuyết bám sát thực nghiệm từ giai đoạn tươi ($\text{MC} > 150%$) về độ ẩm cuối cùng $12%$, khẳng định vai trò của việc tích hợp đồng thời hệ số thấm $K$ và khuếch tán $D_b$.

  4. Ưu thế vượt trội của cơ chế sấy gián đoạn bằng NLMT: Quá trình sấy trong lò sấy NLMT với dốc sấy tối ưu $U = 3,2$ cho thấy sự tương thích hoàn hảo với đặc tính giải phẫu gỗ Keo tai tượng. Chu kỳ nhiệt ngày - đêm tạo ra khoảng thời gian nghỉ ban đêm (không gia nhiệt), cho phép ẩm tự do và ẩm liên kết từ lõi tái phân bố ra bề mặt, triệt tiêu dốc ẩm đỉnh và giải tỏa nội ứng suất dư.

Thông số so sánh Lò sấy hơi nước quy chuẩn ($35\text{ m}^3$) Lò sấy năng lượng mặt trời ($30\text{ m}^3$) Mức độ cải thiện / Đánh giá
Dung tích mẻ sấy $35\text{ m}^3$ $30\text{ m}^3$ Quy mô sản xuất công nghiệp
Thời gian sấy ($\text{MC } 140% \to 12%$) $14 - 16\text{ ngày}$ $22 - 25\text{ ngày}$ Thời gian sấy dài hơn do gián đoạn
Tỷ lệ khuyết tật nứt mặt ($NM_i$) $3,5 - 5,2%$ $< 0,5%$ Giảm khuyết tật rõ rệt
Tỷ lệ khuyết tật mo móp / sụp đổ Trung bình - Cao Rất thấp / Không đáng kể Tự cân bằng ẩm ban đêm
Chi phí năng lượng & vận hành $100%$ (Mức cơ sở) Giảm $31,3%$ Tiết kiệm nhiên liệu đốt
Chất lượng gỗ sau sấy Đạt loại B - C Đạt loại A (AS/NZS 4787) Đồng đều, ứng suất dư thấp

Implications đa chiều

  • Về mặt học thuật: Đặt nền tảng cho trường phái nghiên cứu sấy gỗ định lượng tại Việt Nam, chứng minh tính khả thi của việc ứng dụng mô hình số môi trường xốp TransPore cho các loài gỗ rừng trồng nhiệt đới.
  • Về mặt phương pháp luận: Cung cấp quy chuẩn đo dốc ẩm chính xác cao bằng kỹ thuật cắt lát vi phân kết hợp phân tích ảnh giải phẫu hiển vi kỹ thuật số ImageJ.
  • Về mặt công nghệ và thực tiễn: Đề xuất hoàn chỉnh quy trình công nghệ sấy gỗ xẻ Keo tai tượng dày $25\text{ mm}$ bằng lò sấy NLMT với dốc sấy $U = 3,2$, kiểm soát $\text{EMC}$ theo từng cấp độ ẩm.
  • Về mặt chính sách và môi trường: Cung cấp luận cứ khoa học để Bộ Nông nghiệp & PTNT và các hiệp hội chế biến gỗ (VIFOREST) thúc đẩy chuyển đổi xanh, áp dụng công nghệ sấy năng lượng tái tạo nhằm giảm phát thải khí nhà kính trong chế biến lâm sản.

Limitations và Future Research

Luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  • Mô hình TransPore kế thừa từ Redman (2017) áp dụng hệ số đẳng nhiệt nhả ẩm trung bình bằng 1, chưa đo đạc đường cong đẳng nhiệt hấp phụ/nhả ẩm thực tế để tính toán khoảng trễ hấp phụ ($W$) riêng cho gỗ Keo tai tượng.
  • Giới hạn mô phỏng số học dừng lại ở mô hình chuyển khối 2D trên mặt cắt ngang, chưa tích hợp hoàn chỉnh dòng vận chuyển dọc thớ 3D và chưa liên kết trực tiếp với phân tích phần tử hữu hạn (FEA) dự đoán ứng suất 3D cho lò sấy NLMT.
  • Về nguyên liệu, nghiên cứu tập trung trên gỗ Keo tai tượng 9 tuổi tại một vùng lập địa (Vĩnh Phúc), chưa bao quát hết biến động sinh thái theo các vùng sinh thái khác nhau tại Việt Nam.

Hướng nghiên cứu tiếp theo bao gồm:

  1. Xác định đường đẳng nhiệt trao đổi ẩm thực nghiệm và khoảng trễ $W$ của Keo tai tượng ở các dải nhiệt độ khác nhau.
  2. Xây dựng mô hình 3D tích hợp đầy đủ trường chuyển khối và trường ứng suất biến dạng nhớt dẻo cho lò sấy NLMT.
  3. Ứng dụng hệ thống điều khiển thông minh tự động (IoT/AI) điều tiết cửa thoát ẩm và quạt tuần hoàn trong lò sấy NLMT dựa trên dự báo bức xạ thời tiết.
  4. Mở rộng nghiên cứu cho các độ dày ván lớn hơn ($35\text{ mm}, 50\text{ mm}$) và các loài gỗ rừng trồng mọc nhanh khác như Keo lai, Bạch đàn urô (Eucalyptus urophylla).

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

Luận án mang lại tác động toàn diện trên cả phương diện khoa học và thực tiễn:

  • Tác động học thuật: Mở ra hướng nghiên cứu liên ngành giữa cơ học gỗ, vật lý chuyển khối và toán tin ứng dụng; cung cấp hệ thống số liệu chuẩn hóa được công bố trên các tạp chí khoa học uy tín trong nước và quốc tế.
  • Chuyển đổi công nghiệp: Giúp các doanh nghiệp chế biến gỗ nâng cao tỷ lệ thu hồi gỗ xẻ loại A từ Keo tai tượng lên trên $85%$, giảm tỷ lệ loại bỏ do nứt vỡ và mo móp, nâng cao giá trị thương phẩm của gỗ rừng trồng Việt Nam trong chuỗi cung ứng toàn cầu.
  • Hiệu quả kinh tế - môi trường: Tiết kiệm hơn $30%$ chi phí năng lượng sấy so với lò sấy hơi nước truyền thống, giảm phụ thuộc vào nhiên liệu hóa thạch hoặc củi đốt, đóng góp trực tiếp vào mục tiêu Net-Zero của ngành lâm nghiệp.

Đố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 sấy gỗ chuẩn mực quốc tế, từ kỹ thuật vi phẫu, đo thấm/khuếch tán đến mô phỏng số học TransPore.
  • Giảng viên và nhà khoa học lâm sản: Khai thác bộ cơ sở dữ liệu thực nghiệm chi tiết về hằng số vật lý và chuyển khối của Keo tai tượng để phục vụ giảng dạy và phát triển các mô hình lý thuyết.
  • Kỹ sư công nghệ và R&D ngành chế biến gỗ: Áp dụng trực tiếp bảng chế độ sấy NLMT ($U = 3,2$) và các ngưỡng kiểm soát dốc ẩm vào vận hành lò sấy thực tế tại nhà máy.
  • Các nhà hoạch định chính sách lâm nghiệp: Căn cứ vào hiệu quả của công nghệ sấy NLMT để xây dựng các chương trình hỗ trợ chuyển giao công nghệ sạch, nâng cao giá trị gia tăng cho gỗ rừng trồng.

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 thiết lập mối quan hệ định lượng giải thích bản chất giữa mức chênh lệch độ ẩm vi phân theo chiều dày ($\Delta\text{MC}$) với sự hình thành khuyết tật sấy, mở rộng lý thuyết nội ứng suất của McMillen (1958) và cơ chế chuyển khối của Perré & Turner (1999) trên đối tượng gỗ rừng trồng nhiệt đới mọc nhanh.

2. Phương pháp nghiên cứu có cải tiến gì đột phá so với các công trình trước đây?

So với các nghiên cứu trong nước trước đây vốn chỉ dùng ẩm kế điện trở đo 2 điểm bề mặt - tâm (Nguyễn Xuân Hiên, 2006; Đỗ Văn Bản, 2012), luận án đã: (1) Sử dụng kỹ thuật cắt lát vi phân cho độ chính xác cao từ trạng thái gỗ tươi ($\text{MC} > 140%$, $R^2 = 0,8795$); (2) Đo trực tiếp hệ số thấm $K$ trên Porolux 1000 và khuếch tán $D_b$ trên PVC-CHA; (3) Ứng dụng mô phỏng số TransPore 2D trên Matlab.

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

Phát hiện thú vị nhất là dốc ẩm nội tại của Keo tai tượng không tăng tuyến tính mà đạt đỉnh cực đại khi độ ẩm trung bình ở mức $\sim 50%$. Đồng thời, tỷ lệ $D_{fl}/T_{fw}$ của sợi gỗ cao kết hợp với tính thấm ngang cực thấp là nguyên nhân cốt lõi gây sụp đổ tế bào (mo móp) ngay ở giai đoạn đầu nếu sấy cưỡng bức với dốc sấy $U > 4,0$.

4. Quy trình thực nghiệm có cung cấp đầy đủ giao thức để tái lập (replication protocol) không?

Có. Luận án mô tả chi tiết kích thước mẫu, quy chuẩn bịt đầu bằng keo epoxy, các bước vận hành máy đo Porolux 1000, nồng độ dung dịch muối trong hệ PVC-CHA, thông số lưới tam giác trong mô hình số Matlab, và từng bước vận hành lò sấy NLMT $30\text{ m}^3$ với dốc sấy $U = 3,2$.

5. Chương trình nghiên cứu 10 năm tiếp theo từ nền tảng luận án này nên tập trung vào đâu?

Trọng tâm 10 năm tới gồm: (1) Tích hợp mô hình chuyển khối 3D với mô hình cơ học phá hủy phần tử hữu hạn (FEA Fracture Mechanics); (2) Xây dựng hệ thống điều khiển sấy thông minh AIoT thích ứng thời tiết thực cho lò sấy NLMT; (3) Mở rộng cơ sở dữ liệu chuyển khối cho toàn bộ các loài gỗ rừng trồng chủ lực tại Việt Nam.

Kết luận

  1. Luận án đã làm sáng tỏ toàn diện bản chất cấu tạo hiển vi và các thông số chuyển khối của gỗ Keo tai tượng 9 tuổi: xác định tính dị hướng sâu sắc của hệ số thấm ($K_L \gg K_R, K_T$) và hệ số khuếch tán ($D_{bL} \gg D_{bR}, D_{bT}$), lý giải cơ chế mo móp do tỷ lệ $D_{fl}/T_{fw}$ và vách sợi mỏng.
  2. Thiết lập thành công phương pháp cắt lát vi phân đo dốc ẩm theo chiều dày tấm gỗ với độ tin cậy vượt trội ($R^2 = 0,8795$ so với đo điện trở), xác định quy luật biến thiên dốc ẩm đạt đỉnh tại $\text{MC} \approx 50%$.
  3. Xây dựng hệ phương trình tương quan định lượng giữa mức chênh lệch độ ẩm, $\text{EMC}$ và mức độ nứt vỡ ($N_i$), từ đó xác định chính xác ngưỡng dốc ẩm an toàn ($\Delta\text{MC} \le 28,5%$ tại $\text{MC} = 50%$ và $\le 11,2%$ tại $\text{MC} = 20%$) để đạt tỷ lệ khuyết tật $0%$.
  4. Ứng dụng và kiểm chứng thành công mô hình toán học TransPore 2D mô phỏng trường ẩm và nhiệt độ trong lò sấy quy chuẩn với sai số thấp, khẳng định tính chính xác của mô hình chuyển khối kết hợp Darcy - Fick.
  5. Đề xuất quy trình công nghệ sấy gỗ Keo tai tượng dày $25\text{ mm}$ bằng lò sấy NLMT đạt dốc sấy tối ưu $U = 3,2$, tận dụng chu kỳ nghỉ ban đêm để triệt tiêu dốc ẩm và giải tỏa ứng suất dư, giúp giảm $31,3%$ chi phí năng lượng và nâng cao chất lượng gỗ sấy phục vụ chế biến công nghiệp.