Tổng quan luận án

Luận án tiến sĩ kỹ thuật của tác giả Võ Tuyển với đề tài "Nghiên cứu, ứng dụng dòng phun rối xoáy trong hệ thống thiết bị tưới phun", thuộc chuyên ngành Công nghệ chế tạo máy (mã số: 62 52 04 01), được bảo vệ năm 2013 tại Trường Đại học Bách khoa – Đại học Quốc gia Thành phố Hồ Chí Minh dưới sự hướng dẫn khoa học của PGS. Nguyễn Thanh Nam và PGS. Trần Thị Hồng.

Tính cấp thiết và khoảng trống nghiên cứu

Kỹ thuật tưới phun là phương thức cung cấp nước dưới dạng hạt mưa hoặc hạt sương, đóng vai trò quan trọng trong việc thỏa mãn nhu cầu sinh lý của cây trồng, tiết kiệm nước và năng lượng. Tuy nhiên, hai dòng thiết bị tưới phun thông dụng thời điểm nghiên cứu đều bộc lộ hạn chế kỹ thuật:

  • Đầu phun tự xoay kiểu búa đập có bán kính tưới lớn (6 đến 20 m) nhưng yêu cầu áp lực làm việc cao (trên 2,5 bar) và độ đồng đều tưới phun không cao, hạt nước thô dễ gây tổn thương cây trồng mềm yếu.
  • Đầu phun cố định tỏa tròn có hạt mưa đồng đều hơn nhưng bán kính phun hẹp (dưới 5 m), phụ thuộc vào núm chặn cơ học.

Khoảng trống nghiên cứu được tác giả xác định: Các nghiên cứu về dòng chảy rối xoáy trước đây hầu như chỉ tập trung vào lĩnh vực buồng đốt công nghiệp hoặc thiết bị trao đổi nhiệt nhằm ổn định ngọn lửa và tăng cường hòa trộn nhiên liệu - không khí. Việc nghiên cứu một cách có hệ thống về dòng phun rối xoáy ứng dụng trong hệ thống thiết bị tưới phun nông nghiệp vẫn còn bỏ ngỏ. Đặc biệt, việc mô hình hóa toán học có tích hợp quá trình trao đổi chất (nồng độ pha lỏng) để mô tả dòng hai pha đồng nhất, cũng như việc chế tạo đầu phun tạo xoáy bằng rãnh xoáy cố định để tối ưu hóa thủy lực tưới tiêu chưa được giải quyết đầy đủ.

Mục tiêu nghiên cứu

Mục đích bao trùm của luận án là nghiên cứu bản chất thủy lực và ứng dụng dòng phun rối xoáy vào hệ thống thiết bị tưới phun nhằm mở rộng bán kính phun, tăng độ đồng đều tưới, giảm áp lực làm việc và bảo đảm tiết kiệm năng lượng.

Năm nhiệm vụ nghiên cứu cụ thể được đánh số theo văn bản luận án gồm:

  1. Xây dựng mô hình toán cho dòng phun rối xoáy ứng dụng trong kỹ thuật tưới phun và xây dựng thuật toán giải hệ phương trình toán bằng phương pháp khối hữu hạn.
  2. Mô phỏng số các thông số vật lý của dòng phun rối xoáy và xác định mối quan hệ giữa các thông số vật lý, thông số hình học của dòng phun với hệ số cường độ xoáy.
  3. Thiết kế, chế tạo đầu phun tạo xoáy và xây dựng mô hình thực nghiệm đo đạc các thông số kích thước hình học và thông số vật lý của dòng phun.
  4. Nghiên cứu, phân tích tác động của hiệu ứng xoáy tới các thông số kỹ thuật của dòng phun và bằng phương pháp quy hoạch thực nghiệm, xác định các thông số tối ưu của quá trình tưới phun gồm đường kính vòi phun, hệ số xoáy và lưu lượng tưới.
  5. Ứng dụng tính toán, thiết kế và lắp đặt hệ thống tưới phun mưa sử dụng đầu phun tạo xoáy cho một mô hình cụ thể và tiến hành so sánh kết quả ứng dụng thực nghiệm hệ thống tưới phun tại ba đơn vị khác nhau khi sử dụng đầu phun tạo xoáy.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Dòng phun rối xoáy trong hệ thống thiết bị tưới phun và tác động của hiệu ứng xoáy tới trường dòng chảy.
  • Phạm vi nghiên cứu:
    • Khảo sát dòng phun rối xoáy chuyển động dừng, phun từ miệng vòi tròn vào môi trường không khí tĩnh.
    • Về mặt không gian mô phỏng: Giới hạn ở giai đoạn đầu của dòng phun từ miệng vòi ($x = 0$) đến độ cao cực đại ($x = H$) trước khi các giọt nước rơi tự do.
    • Về mặt thực nghiệm và ứng dụng: Thiết kế, chế tạo vòi phun với đường kính vòi từ 2,0 mm đến 4,5 mm; khảo sát tại phòng thí nghiệm và lắp đặt đối chứng thực tế tại 3 cơ sở sản xuất, đào tạo thuộc khu vực Thành phố Hồ Chí Minh, Đồng Nai và Bình Dương.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước đây được điểm lại

Trong luận án, tác giả tổng hợp tình hình nghiên cứu qua ba nhóm chuyên đề:

  1. Về kỹ thuật tưới phun:

    • Trên thế giới: Sterling David (thập niên 1960) nghiên cứu so sánh kỹ thuật tưới phun với các kỹ thuật tưới tiết kiệm nước khác tại Mỹ. Giữa thập niên 1970, các nhà khoa học tại Israel và Australia bắt đầu chế tạo vòi phun lưu lượng nhỏ. Don Olson (cuối thập niên 1970) đóng góp quan trọng trong việc hoàn thiện đầu phun mưa lưu lượng nhỏ. Dữ liệu của Tổ chức Lương thực và Nông nghiệp Liên hiệp quốc (FAO) ghi nhận trên 1.000 ha cây trồng áp dụng công nghệ này.
    • Tại Việt Nam: Bắt đầu nghiên cứu ứng dụng từ năm 1993. Năm 1994, Trường Đại học Thủy lợi triển khai thử nghiệm hệ thống phun mưa tại Phủ Quỳ, Nghệ An. Tiếp đó là các thử nghiệm đầu phun sương bằng nhựa (Micro Sprinkle) mô phỏng mẫu từ Israel. Bài báo "Vòi phun và công nghệ tưới phun mưa" (2007) khảo sát góc độ của vòi phun nhưng dừng lại ở dòng phun không xoáy.

  2. Về mô hình dòng chảy rối:

    • Boussinesq (1887) thiết lập giả thiết mối liên hệ tuyến tính giữa ứng suất rối và tensor tốc độ biến dạng thông qua độ nhớt rối động lực $\mu_t$ đẳng hướng.
    • Reynolds (1894) phân tách vận tốc tức thời thành thành phần trung bình thời gian và thành phần mạch động.
    • Prandtl (1925) đề xuất mô hình chiều dài xáo trộn $l_m$; tiếp đó năm 1945 đưa ra mô hình một phương trình cho động năng rối $k$.
    • Launder và Jones thiết lập mô hình hai phương trình $k-\varepsilon$ chính tắc, trở thành công cụ phổ biến cho dòng chảy rối một pha và hai pha đồng nhất.
    • Tại Việt Nam: Mô hình chiều dài xáo trộn được áp dụng giải bài toán lan truyền chất trên sông bằng lược đồ sai phân ADI; mô hình độ nhớt rối mô phỏng lan truyền khói bụi địa hình phức tạp; mô hình Boussinesq nghiên cứu sóng tàu ven bờ; mô hình $k-\varepsilon$ khảo sát hải lưu bồi lắng đáy biển, khí động lực học quanh cánh máy bay, tia phun trong buồng đốt động cơ đốt trong, dòng chảy qua dàn ống phun mưa, và hoạt động của ống phun hỏa tiễn trong điều kiện quá giãn nở.

  3. Về dòng phun rối xoáy:

    • Peckham và Atkinson (1957), Harvey (1962) nghiên cứu dòng xoáy sau cánh và xác định cường độ xoáy qua tỷ số vận tốc tiếp tuyến và vận tốc dọc trục.
    • McManus, Legner và David (1994) khảo sát hiệu quả vòi phun xoáy dao động ngắt quãng.
    • Chen và Driscoll (1998) nghiên cứu dòng xoáy mạnh ($S = 1,0$) hình thành vùng hồi lưu làm giảm chiều dài ngọn lửa 5 lần.
    • Lin (1998) phân tích mô hình tia đồng trục xoáy dựa trên hệ phương trình trung bình Reynolds (RANS).
    • A. Syred (1984), Jiri Brada (2001), Vaidyanathan Sankaran và Suresh Menon (2002) thiết lập mô hình toán dòng xoáy biểu diễn qua hàm dòng $\psi$ trong tọa độ trụ.
    • Jian Zhang và Chengkai Zhu (2009) đưa thông số nồng độ (ACM) vào mô phỏng dòng phun rối và cháy trong buồng đốt xoáy.
    • Tại Việt Nam: Các tác giả khảo sát dòng hai pha rối xoáy không đẳng nhiệt trong buồng đốt công nghiệp, thiết bị đầu đốt tạo xoáy và máy lạnh hiệu ứng xoáy Rank.
Tác giả / Nhóm tác giả Năm Nội dung trọng tâm Phương pháp / Mô hình
Boussinesq 1887 Khái niệm ứng suất rối, giả thiết độ nhớt rối Giải tích lý thuyết
Reynolds 1894 Phân tách đại lượng trung bình và mạch động Thống kê dòng rối
Prandtl 1925, 1945 Chiều dài xáo trộn $l_m$, mô hình một phương trình $k$ Bán thực nghiệm
Launder & Jones 1972 Mô hình hai phương trình vi phân $k-\varepsilon$ Số trị, RANS
Peckham, Harvey 1957, 1962 Khảo sát dòng xoáy sau cánh, tỷ số vận tốc $w/u$ Thực nghiệm thủy động
Chen & Driscoll 1998 Vùng hồi lưu ở hệ số xoáy mạnh $S = 1,0$ Thực nghiệm buồng đốt
Syred, Brada, Menon 1984, 2002 Thiết lập mô hình toán dòng xoáy dùng hàm dòng $\psi$ Phương trình vi phân
Schraiber 1980 Phương trình trao đổi chất của pha phân tán Vận tải nồng độ hạt

Khoảng trống luận án lựa chọn giải quyết

Luận án kết nối hai nhánh lý thuyết: cơ học dòng chảy rối xoáy và kỹ thuật tưới phun nông nghiệp. Thay vì chỉ áp dụng hiệu ứng xoáy trong buồng đốt ở nhiệt độ cao, tác giả đưa hiệu ứng này vào vòi phun thủy lực nước lạnh. Để giải quyết việc thiếu hụt mô hình toán học phù hợp cho tưới phun, tác giả bổ sung phương trình vận tải nồng độ hạt nước theo Schraiber vào hệ phương trình RANS và mô hình $k-\varepsilon$, tạo nên một hệ 7 phương trình khép kín giải bằng phương pháp khối hữu hạn (FVM).

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết và khung phân tích

Dòng phun nước từ đầu phun vào không khí cuốn theo không khí xung quanh, đồng thời các giọt nước phân tán vào môi trường tạo thành dòng hai pha nước - không khí. Do các bọt khí và giọt lỏng hòa trộn đồng đều, không có sự chênh lệch lớn về vận tốc và áp suất giữa hai pha, tác giả thiết lập giả thiết xem dòng phun là dòng hai pha đồng nhất, chuyển động dừng ($\partial/\partial t = 0$) và có tính chất đối xứng trục.

Hệ tọa độ trụ $(x, r, \theta)$ được sử dụng với $x$ là trục dọc theo hướng tia phun, $r$ là phương hướng kính, và $\theta$ là góc vòng quanh trục $x$. Vì dòng đối xứng trục, các đạo hàm theo $\theta$ triệt tiêu ($\partial/\partial\theta = 0$).

Hệ phương trình vi phân đạo hàm riêng gồm 7 phương trình khép kín:

  1. Phương trình liên tục: $$\frac{\partial (\rho u)}{\partial x} + \frac{1}{r} \frac{\partial (r \rho v)}{\partial r} = 0$$
  2. Phương trình cân bằng động lượng theo phương dọc trục ($x$): Xét ảnh hưởng của trọng lực ($F_x = -g$ hoặc $+g$ tùy chiều đặt vòi), áp suất $p$ và độ nhớt rối $\nu_t$.
  3. Phương trình cân bằng động lượng theo phương hướng kính ($r$): Chứa thành phần lực ly tâm do vận tốc tiếp tuyến $w^2/r$.
  4. Phương trình cân bằng động lượng theo phương tiếp tuyến ($\theta$): Biểu diễn sự biến thiên của vận tốc xoáy $w$.
  5. Phương trình trao đổi chất (nồng độ khối lượng pha nước $C$): Dựa trên mô hình Schraiber và định luật Fick, có xét đến vận tốc rơi tự do của hạt nước $u_s$: $$u \frac{\partial C}{\partial x} + v \frac{\partial C}{\partial r} = \frac{\partial}{\partial x}\left(\frac{\nu_t}{Sc}\frac{\partial C}{\partial x}\right) + \frac{1}{r}\frac{\partial}{\partial r}\left(r \frac{\nu_t}{Sc}\frac{\partial C}{\partial r}\right) - u_s \frac{\partial C}{\partial x}$$ Trong đó, số Schmidt chọn $Sc = 0,7$; đường kính hạt mưa trung bình chọn $d = 0,115\text{ mm}$; tỷ trọng $s = \rho_w/\rho_g = 1000 / 1,228 \approx 814,3$; vận tốc rơi Stokes $u_s = \frac{(s - 1) g d^2}{18 \nu}$.
  6. Phương trình động năng rối ($k$): Vận tải động năng rối với số hạng nguồn sinh rối $P_k$ và tiêu tán $\varepsilon$.
  7. Phương trình tốc độ tiêu tán động năng rối ($\varepsilon$): Đóng kín hệ phương trình với độ nhớt rối động học $\nu_t = C_\mu \frac{k^2}{\varepsilon}$.

Các hằng số mô hình chuẩn được sử dụng: $$C_\mu = 0,09; \quad \sigma_k = 1,00; \quad \sigma_\varepsilon = 1,30; \quad C_1 = 1,44; \quad C_2 = 1,92$$

Hệ số cường độ xoáy $S$ đặc trưng cho dòng xoáy được tính gần đúng theo công thức của A. K. Gupta: $$S = \frac{2}{3} \frac{w}{u} = \frac{2}{3} \tan\alpha$$ Với $\alpha$ là góc nghiêng của rãnh tạo xoáy trong đầu phun. Dòng được phân chia thành xoáy yếu khi $S < 0,6$ và xoáy mạnh khi $S > 0,6$ (ngưỡng tới hạn $S_{th} = 0,6$).

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

  • Phương pháp mô phỏng số: Đưa hệ phương trình 7 ẩn ($u, v, w, p, C, k, \varepsilon$) về dạng không thứ nguyên tổng quát. Sử dụng phương pháp khối hữu hạn (FVM) trên lưới cấu trúc 2D, lập trình thuật toán lặp giải hệ phương trình sai phân trên phần mềm Matlab 2008a.
  • Phương pháp thực nghiệm:
    • Thiết kế và chế tạo các đầu phun xoáy bằng phương pháp gia công cơ khí chính xác với các góc nghiêng rãnh $\alpha$ khác nhau tạo ra các hệ số xoáy $S = 0; 0,2; 0,4; 0,7; 1,2; 2,0$ và các đường kính miệng vòi $d = 3,0\text{ mm}; 3,5\text{ mm}; 4,0\text{ mm}; 4,5\text{ mm}$.
    • Thiết bị đo đạc: Bơm ly tâm cung cấp nước công suất $N = 370\text{ W}$, tốc độ quay $n = 2900\text{ vg/ph}$, lưu lượng danh định $Q = 25\text{ l/ph}$; đồng hồ đo lưu lượng dải đo 1 - 100 l/ph (độ chính xác 0,1 - 0,2%); đồng hồ đo áp suất dải đo 0,1 - 5 bar (độ chính xác 0,1%); các cốc hứng nước tiêu chuẩn đường kính miệng 70 mm, chiều cao 90 mm bố trí theo mạng lưới tỏa tròn để xác định cường độ mưa và độ đồng đều tưới ($CU$); thước dây đo bán kính phun $R$ và chiều cao dòng phun $H$.
  • Phương pháp quy hoạch thực nghiệm: Sử dụng ma trận quy hoạch thực nghiệm trực giao bậc một và bậc hai nhằm thiết lập các phương trình hồi quy mô tả tương tác giữa 3 thông số vào ($X_1$: đường kính vòi phun $d$; $X_2$: hệ số xoáy $S$; $X_3$: áp suất/lưu lượng $p, Q$) tới 3 hàm mục tiêu ($Y_1$: bán kính tưới phun $R$; $Y_2$: độ đồng đều tưới $CU$; $Y_3$: chi phí điện năng tiêu thụ). Kiểm định phương sai bằng tiêu chuẩn thống kê (ANOVA) để loại bỏ các hệ số không có nghĩa và giải bài toán tối ưu đa mục tiêu.

Nội dung chính theo từng chương

Chương 1: Tổng quan, mục tiêu và nhiệm vụ nghiên cứu

Chương này tập trung trình bày bức tranh tổng quát về công nghệ tưới phun trên thế giới và tại Việt Nam, phân loại các hệ thống tưới phun (cố định, di động, bán di động). Tác giả so sánh cấu tạo và đặc tính làm việc của hai loại đầu phun phổ biến: loại tự xoay búa đập (bán kính tưới rộng 6 - 20 m nhưng cần áp lực cao $> 2,5\text{ bar}$, phân bố kém đồng đều) và loại tỏa tròn có núm chặn (áp lực thấp 1,2 - 2 bar, hạt tương đối đồng đều nhưng bán kính $< 5\text{ m}$).

Phần tiếp theo trình bày bản chất vật lý của chuyển động rối, quá trình hình thành dòng rối xoáy từ rãnh dẫn hướng, ống tiếp tuyến hoặc đầu phun tự quay. Tác giả định nghĩa hệ số cường độ xoáy $S$, phân tích sự khác nhau giữa dòng xoáy yếu ($S < 0,6$) không tạo dòng tuần hoàn ngược và dòng xoáy mạnh ($S > 0,6$) hình thành vùng áp suất âm và dòng hồi lưu trung tâm ở miệng phun. Trên cơ sở tổng quan tài liệu trong và ngoài nước, tác giả xác định 4 nội dung nghiên cứu lớn và hệ thống các phương pháp thực hiện từ lý thuyết, mô phỏng số, chế tạo đến đo đạc thực nghiệm.

Chương 2: Nghiên cứu mô phỏng số dòng phun rối xoáy

Chương 2 là phần trọng tâm về mặt lý thuyết mô hình hóa. Tác giả thiết lập hệ phương trình Navier-Stokes trung bình hóa Reynolds (RANS) kết hợp phương trình bảo toàn khối lượng, phương trình nồng độ Schraiber và mô hình rối hai phương trình $k-\varepsilon$ trong hệ tọa độ trụ đối xứng 2D. Hệ phương trình được chuẩn hóa không thứ nguyên và chuyển về dạng phương trình truyền vận tổng quát cho biến $\Phi$: $$\frac{\partial (u \Phi)}{\partial x} + \frac{1}{r}\frac{\partial (r v \Phi)}{\partial r} = \frac{\partial}{\partial x}\left(\Gamma_\Phi \frac{\partial \Phi}{\partial x}\right) + \frac{1}{r}\frac{\partial}{\partial r}\left(r \Gamma_\Phi \frac{\partial \Phi}{\partial r}\right) + S_\Phi$$

Các điều kiện biên được xác lập chặt chẽ:

  • Tại miệng phun ($x = 0$): Giá trị $u, w, p$ cho trước theo chế độ làm việc;
  • Tại biên ngoài dòng phun ($r = R$): $u = v = w = p = C = k = \varepsilon = 0$;
  • Tại độ cao giới hạn ($x = H$): Toàn bộ các biến triệt tiêu về 0 trước khi hạt nước chuyển sang trạng thái rơi tự do;
  • Tại trục đối xứng ($r = 0$): $v = 0$, $w = 0$, đạo hàm riêng theo bán kính của các biến khác bằng không ($\partial u/\partial r = \partial p/\partial r = \partial C/\partial r = \partial k/\partial r = \partial \varepsilon/\partial r = 0$).

Tác giả áp dụng phương pháp khối hữu hạn (FVM), chia lưới miền tính toán và lập trình giải thuật toán trên Matlab 2008a. Kết quả mô phỏng số cung cấp trường phân bố chi tiết của:

  • Đường dòng và hàm dòng $\psi$ tại $S = 0$ và $S = 0,7$;
  • Vận tốc dọc trục $u/u_0$, vận tốc hướng kính $v/v_0$, vận tốc tiếp tuyến $w/w_0$;
  • Áp suất tĩnh $p/p_0$;
  • Động năng rối $k/k_0$ và tỷ lệ tiêu tán rối $\varepsilon/\varepsilon_0$;
  • Nồng độ pha lỏng $C/C_0$.

Cuối chương, tác giả so sánh kết quả mô phỏng số của mình với các dữ liệu thực nghiệm đã công bố của các tác giả đi trước (như A. K. Gupta, Syred) về sự suy giảm vận tốc trục $u_m$, phân bố vận tốc tiếp tuyến $w_m$, phân bố áp suất trên trục $p_m$, và nồng độ tại tiết diện ngang $x/d = 4$, cho thấy sự tương thích tốt về quy luật biến thiên thủy lực.

Chương 3: Nghiên cứu thực nghiệm dòng phun rối xoáy trong kỹ thuật tưới phun

Chương 3 trình bày các bước chế tạo mô hình thực nghiệm và kiểm chứng tính toán số. Đầu phun tạo xoáy được chế tạo với các chi tiết rãnh nghiêng góc $\alpha$ khác nhau nhằm thay đổi hệ số xoáy $S$.

Tác giả tiến hành đo đạc các thông số hình học và thủy lực ở các mức áp suất $p = 1,2\text{ bar}; 2,0\text{ bar}; 2,2\text{ bar}$ với đường kính vòi phun $d = 3,0\text{ mm}; 3,5\text{ mm}; 4,0\text{ mm}; 4,5\text{ mm}$. Các dữ liệu thực nghiệm thu được gồm:

  • Góc mở dòng phun $\beta$, bán kính tưới $R$, lưu lượng qua vòi $Q$;
  • Phân bố cường độ mưa $i$ (mm/h) và chỉ số độ đồng đều tưới $CU$ (%);
  • Độ thô hạt mưa $K$.

Kiểm chứng giữa mô hình số và thực nghiệm cho thấy mức độ sai lệch nhỏ đối với góc phun, bán kính dòng phun và lưu lượng. Kết quả xác định mối quan hệ giữa hệ số xoáy $S$ với các chỉ tiêu tưới:

  • Khi tăng $S$ từ 0 lên 1,2, góc phun mở rộng, diện tích phủ mưa tăng lên rõ rệt, hạt mưa mịn hơn (độ thô giảm), độ đồng đều $CU$ tăng;
  • Tuy nhiên nếu $S$ quá lớn ($S > 1,2$), lực ly tâm làm tia nước tản rộng ngay sát miệng phun khiến chiều cao và tầm phun xa bị suy giảm nhanh.

Áp dụng phương pháp quy hoạch thực nghiệm đa yếu tố:

  • Xây dựng ma trận thực nghiệm bậc một và bậc hai cho 3 thông số ảnh hưởng: $X_1$ (đường kính lỗ vòi), $X_2$ (hệ số xoáy), $X_3$ (lưu lượng/áp lực).
  • Phân tích phương sai ANOVA cho 3 hàm mục tiêu: Bán kính tưới $Y_1$, độ đồng đều tưới $Y_2$, và chi phí năng lượng $Y_3$.
  • Loại bỏ các hệ số hồi quy không có ý nghĩa thống kê, xây dựng các mặt đáp 3D ($X_1-X_2$, $X_1-X_3$, $X_2-X_3$).
  • Giải bài toán tối ưu tìm ra bộ thông số kỹ thuật hợp lý bảo đảm bán kính tưới lớn nhất, độ đồng đều cao nhất trong khi mức tiêu hao điện năng là thấp nhất.

Chương 4: Ứng dụng thiết kế hệ thống tưới phun

Chương 4 chuyển hóa các kết quả nghiên cứu lý thuyết và thực nghiệm vào bài toán thiết kế hệ thống tưới phun hoàn chỉnh cho đối tượng cây trồng cụ thể.

Nội dung thiết kế thủy lực bao gồm:

  1. Xác định nhu cầu tưới và chế độ tưới dựa trên tốc độ thấm của từng loại đất và nhu cầu sinh lý cực đại của cây trồng;
  2. Lựa chọn dạng bố trí đầu phun (hình vuông hoặc tam giác) căn cứ vào bán kính tưới tối ưu $R$ xác định ở Chương 3;
  3. Tính toán lưu lượng trên các đoạn ống chính, ống nhánh và chọn đường kính ống theo vận tốc dòng chảy kinh tế $v_{KT}$;
  4. Xác định tổn thất cột áp ma sát dọc đường và tổn thất cục bộ trên đường ống cơ bản $h_{wd}$;
  5. Tính toán chiều cao hút, chiều cao đẩy và áp lực làm việc để chọn bơm ly tâm phù hợp;
  6. Thống kê vật tư, phụ kiện và van điều tiết.

Tác giả triển khai khảo sát, lắp đặt đối chứng đầu phun xoáy tại 3 đơn vị:

  • Ứng dụng 1: Mô hình tưới thực nghiệm rau màu và hoa tại Trường Đại học Công nghiệp Thực phẩm TP. Hồ Chí Minh.
  • Ứng dụng 2: Khu nhân giống mía thuộc Công ty Cổ phần Mía Đường La Ngà (Đồng Nai), so sánh thực nghiệm với đầu phun xoay ngẫu lực T92 và đầu phun 098.
  • Ứng dụng 3: Vườn ươm cây lâm nghiệp tại Trường Cao đẳng nghề Công nghệ và Nông Lâm Nam Bộ (Bình Dương), sử dụng bồn cấp nước dung tích 5.000 lít, so sánh trực tiếp với đầu phun tỏa tròn 096.

Kết quả đối chứng tại cả 3 mô hình khẳng định: Đầu phun xoáy hoạt động ổn định ở dải áp suất thấp (1,2 đến 2,0 bar), tạo hạt mưa nhỏ mịn không làm dập nát lá mầm, cải thiện tiểu khí hậu (hạ nhiệt độ, tăng ẩm) và có độ đồng đều phân bố lượng nước cao hơn các đầu phun đối chứng cùng kích cỡ lỗ vòi.

Chương 5: Kết luận chung và hướng phát triển của đề tài

Chương cuối tóm lược các kết quả chính đã hoàn thành của đề tài qua ba nhóm nội dung: lý thuyết mô hình hóa, thực nghiệm chế tạo vòi phun và ứng dụng triển khai thực tế. Đồng thời, tác giả vạch ra các định hướng mở rộng nghiên cứu tiếp theo nhằm hoàn thiện giải pháp kỹ thuật.

Kết quả và những đóng góp mới

Đóng góp mới về mặt lý luận khoa học

  1. Xây dựng hoàn chỉnh mô hình toán học cho dòng phun rối xoáy hai pha đồng nhất ứng dụng trong thiết bị tưới phun, khép kín hệ phương trình động lượng RANS và mô hình rối hai phương trình $k-\varepsilon$ bằng phương trình trao đổi chất (vận tải nồng độ hạt nước theo Schraiber có xét vận tốc rơi Stokes).
  2. Xây dựng thuật toán và chương trình giải số hệ 7 phương trình vi phân phi tuyến dạng không thứ nguyên trong hệ tọa độ trụ bằng phương pháp khối hữu hạn (FVM) trên nền Matlab.
  3. Mô phỏng và làm rõ bản chất thủy lực của dòng phun xoáy: Khi hệ số xoáy vượt qua giá trị tới hạn ($S > 0,6$), gradient áp suất ngược hình thành vùng tuần hoàn áp suất âm cục bộ ở miệng vòi, làm thay đổi cấu trúc vận tốc trục và gia tăng độ phân tán hạt nước theo phương hướng kính.

Đóng góp mới về mặt thực tiễn

  1. Thiết kế và chế tạo thành công mẫu đầu phun tạo xoáy sử dụng rãnh xoáy cố định bên trong vòi phun. Kết cấu này không sử dụng chi tiết chuyển động cơ học (như cơ cấu búa đập hay trục quay ngẫu lực), khắc phục triệt để hiện tượng kẹt cát, mài mòn cơ khí, kéo dài tuổi thọ thiết bị.
  2. Thông qua quy hoạch thực nghiệm đa yếu tố, xác lập được quy luật tác động định lượng của đường kính vòi phun $d$, hệ số xoáy $S$ và áp suất làm việc $p$ tới bán kính tưới $R$, độ đồng đều tưới $CU$, độ thô hạt mưa $K$ và chi phí điện năng.
  3. Xác định được dải thông số làm việc hiệu quả của đầu phun xoáy cho tưới cây trồng mềm yếu: Hệ số xoáy $S$ trong khoảng $0,7$ đến $1,2$; áp suất làm việc từ 1,2 đến 2,0 bar (thấp hơn đáng kể so với mức $> 2,5\text{ bar}$ của đầu phun xoay truyền thống), giúp tiết kiệm điện năng cho trạm bơm.
  4. Đưa ra quy trình tính toán thiết kế thủy lực hoàn chỉnh và kiểm chứng ứng dụng thành công tại 3 cơ sở thực tế (Trường ĐH Công nghiệp Thực phẩm TP.HCM, Công ty CP Mía Đường La Ngà, Trường CĐ Nghề Công nghệ và Nông Lâm Nam Bộ).

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu

Từ phạm vi và các giả thiết được nêu trong văn bản luận án, công trình tồn tại các giới hạn kỹ thuật:

  • Mô hình toán học chỉ xem xét dòng chảy ở trạng thái dừng ($\partial/\partial t = 0$), chưa mô tả được các dao động quá độ khi đóng mở van hoặc dao động áp lực trên đường ống.
  • Coi dòng hai pha nước - không khí là dòng đồng nhất và đẳng hướng, bỏ qua hiện tượng trượt vận tốc giữa hạt nước và không khí ở các vùng biên xa miệng vòi.
  • Miền mô phỏng số chỉ giới hạn trong giai đoạn đầu của dòng phun từ miệng vòi tới độ cao cực đại ($x \le H$), chưa tích hợp bài toán quỹ đạo đường đạn của giọt nước trong giai đoạn rơi tự do chịu ảnh hưởng của các điều kiện khí tượng tự nhiên (vận tốc gió ngoài trời, nhiệt độ bốc hơi).

Hướng nghiên cứu tiếp

Văn bản luận án đề xuất các định hướng phát triển tiếp theo:

  1. Nghiên cứu mở rộng mô hình toán học cho dòng phun rối xoáy ở trạng thái không dừng (unsteady flows) và dòng hai pha không đồng nhất có xét tương tác pha chi tiết.
  2. Khảo sát thực nghiệm ảnh hưởng của gió tự nhiên và nhiệt độ môi trường đến sự phân tán, bay hơi của hạt sương khi tưới trên đồng ruộng quy mô lớn.
  3. Hoàn thiện công nghệ khuôn mẫu chính xác để chế tạo hàng loạt đầu phun tạo xoáy bằng vật liệu polyme kỹ thuật cao, giảm giá thành sản xuất so với gia công kim loại.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên ngành phục vụ trực tiếp cho các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên ngành Cơ kỹ thuật, Thủy khí động lực học, Máy thủy khí: Luận án cung cấp chi tiết hệ phương trình toán học RANS kết hợp mô hình $k-\varepsilon$ và phương trình vận chuyển nồng độ hạt phân tán; phương pháp rời rạc hóa phương trình vi phân từng phần bằng FVM và lưu đồ thuật toán giải trên Matlab (được trình bày ở Chương 2 và Phụ lục 1).
  • Kỹ sư cơ khí chế tạo máy và thiết bị nông nghiệp: Chương 3 và Phụ lục cung cấp toàn bộ kích thước hình học, góc nghiêng rãnh xoáy $\alpha$, kết cấu đầu phun tạo xoáy cố định, cùng các số liệu thực nghiệm đo đạc lưu lượng, áp suất và phân bố cỡ hạt mưa.
  • Kỹ sư thủy lợi, cán bộ kỹ thuật thiết kế hệ thống tưới tiêu: Nội dung Chương 4 cung cấp quy trình tính toán thủy lực mạng lưới đường ống, phương pháp xác định tổn thất áp lực $h_{wd}$, phương pháp bố trí đầu phun xoáy và chọn công suất bơm ly tâm phù hợp với yêu cầu tưới tiết kiệm nước trong nhà kính, vườn ươm cây giống, và vùng canh tác cây công nghiệp.

Câu hỏi thường gặp

1. Điểm khác biệt cơ bản về nguyên lý làm việc giữa đầu phun tạo xoáy trong luận án và các đầu phun truyền thống là gì?

Các đầu phun truyền thống tạo dòng phun mưa nhờ cơ cấu đòn gánh búa đập làm quay thân vòi (đầu phun tự xoay) hoặc dùng núm chặn hình nón đặt phía trước miệng ra (đầu phun tỏa tròn). Đầu phun của luận án tạo hiệu ứng xoáy thông qua các rãnh dẫn hướng cố định có góc nghiêng $\alpha$ bên trong thân vòi, làm lưu chất xuất hiện vận tốc tiếp tuyến $w$ ngay tại miệng ra mà không cần bất kỳ bộ phận cơ khí chuyển động nào.

2. Mô hình toán học trong luận án đã bổ sung phương trình gì để mô tả dòng hai pha trong kỹ thuật tưới phun?

Tác giả đã bổ sung phương trình trao đổi chất (phương trình bảo toàn nồng độ khối lượng của pha nước $C$) dựa trên mô hình của Schraiber và định luật khuếch tán Fick, có tính đến số Schmidt rối ($Sc = 0,7$) và vận tốc rơi tự do Stokes của hạt nước ($u_s$). Phương trình này kết hợp với phương trình liên tục, 3 phương trình động lượng RANS và 2 phương trình mô hình rối $k-\varepsilon$ tạo thành hệ 7 phương trình khép kín.

3. Ngưỡng phân chia mức độ xoáy trong luận án được xác định như thế nào?

Cường độ xoáy được đặc trưng bởi hệ số xoáy $S \approx \frac{2}{3}\tan\alpha$. Luận án xác định hệ số xoáy tới hạn là $S_{th} = 0,6$. Khi $S < 0,6$ dòng thuộc trạng thái xoáy yếu, gradient áp suất ngược không đủ để tạo dòng chảy ngược dọc trục. Khi $S > 0,6$ dòng thuộc trạng thái xoáy mạnh, xuất hiện vùng tuần hoàn áp suất âm gần miệng vòi hút ngược dòng lưu chất, giúp tia phun mở rộng góc tưới nhanh chóng.

4. Quy hoạch thực nghiệm trong luận án khảo sát những thông số và hàm mục tiêu nào?

Nghiên cứu sử dụng quy hoạch thực nghiệm đa yếu tố khảo sát 3 thông số vào: Đường kính lỗ vòi $X_1$ ($d$), hệ số xoáy $X_2$ ($S$), và lưu lượng/áp lực $X_3$ ($Q/p$). Ba hàm mục tiêu đầu ra được đánh giá và tối ưu hóa gồm: Bán kính tưới phun $Y_1$ ($R$), độ đồng đều khi tưới $Y_2$ ($CU$), và chi phí điện năng tiêu thụ $Y_3$.

5. Kết quả ứng dụng thực tế của đầu phun xoáy tại 3 cơ sở đối chứng ghi nhận những ưu điểm gì?

Tại Trường ĐH Công nghiệp Thực phẩm TP.HCM, Công ty CP Mía Đường La Ngà và Trường CĐ Nghề Công nghệ và Nông Lâm Nam Bộ, đầu phun xoáy làm việc ở áp lực thấp (1,2 - 2,0 bar), tạo hạt mưa nhỏ mịn phù hợp cho rau màu, mầm mía và cây giống lâm nghiệp; không làm dập nát bề mặt lá, phân bố lượng nước đồng đều và hỗ trợ cải tạo vi khí hậu mát mẻ cho khu vực canh tác.

Kết luận

Luận án của tác giả Võ Tuyển đã giải quyết trọn vẹn việc nghiên cứu lý thuyết, mô phỏng số, chế tạo thử nghiệm và ứng dụng thực tiễn dòng phun rối xoáy trong hệ thống thiết bị tưới phun nông nghiệp. Công trình đã khép kín mô hình toán học 7 phương trình RANS kết hợp mô hình rối $k-\varepsilon$ và phương trình vận tải nồng độ hạt nước, giải thành công bằng phương pháp khối hữu hạn trên nền Matlab 2008a. Việc chế tạo và tối ưu hóa đầu phun có rãnh tạo xoáy cố định đã chứng minh được hiệu quả kỹ thuật rõ rệt qua việc mở rộng vùng phủ mưa, tăng độ đồng đều và giảm áp lực làm việc, mở ra hướng ứng dụng hiệu quả cho kỹ thuật tưới tiết kiệm nước tại Việt Nam.