Cách tính tốc độ dòng chảy của băng nhỏ giọt và tổn thất áp suất

Jun 09, 2026

Để lại lời nhắn

Băng nhỏ giọttốc độ dòng chảymất áp suấtlà hai con số quyết định xem bạn cóhệ thống tưới nhỏ giọtcung cấp nước đồng đều hoặc để lại những vết khô ở đầu đuôi. Khi xác định kích thước của một hệ thống, bạn cần biết chính xác lượng nước mà mỗi bộ phát cung cấp và lượng áp suất bạn bị mất trong quá trình chạy.

 

Tại sao tính toán lưu lượng băng nhỏ giọt của bạn có thể sai?

Ba lỗi xuất hiện lặp đi lặp lại trong thiết kế thủy lực băng nhỏ giọt:

1. Sử dụng phương trình Hazen-Williams với C=150.Hệ số đó được hiệu chỉnh cho nguồn điện PVC cứng.Băng nhỏ giọt có thành-mỏngvới các đường dẫn dòng chảy mê cung liên tục có hệ số ma sát cao hơn đáng kể. Nghiên cứu được công bố tạiNước(MDPI) đã thử nghiệm hai băng nhỏ giọt có thành mỏng-thương mại và nhận thấy hệ số Blasius phải bằngmột=0.3225 đến 0,3442, không phải tiêu chuẩn 0,3164 được sử dụng cho ống cứng trơn. Sử dụng giá trị trong sách giáo khoa sẽ đánh giá thấp tổn thất do ma sát tới 8%.

2. Bỏ qua hệ số khử Christiansen.Một băng nhỏ giọt có hàng chục hoặc hàng trăm lỗ thoát nước. Nước rời khỏi đường ống ở mỗi đầu phát, do đó tốc độ dòng chảy giảm dọc theo chiều dài. Nếu bạn tính tổn thất do ma sát như thể toàn bộ dòng khí vào di chuyển dọc theo chiều dài, thì bạn sẽ đánh giá quá cao theo hệ số 2–3. Hệ số F{4}} Christiansen khắc phục điều này.

3. Sử dụng tốc độ dòng danh nghĩa mà không tính đến sự thay đổi áp suất.Bộ phát băng nhỏ giọt (loại bù áp suất-không áp suất) tuân theo q=k × H^x. Bộ phát 1,38 L/h được định mức ở cột áp 10 m sẽ chỉ cung cấp khoảng 1,07 L/h ở cột áp 6 m - giảm 22%. Dòng "định mức" chỉ áp dụng ở một áp suất cụ thể.

 

Ⅰ. Làm thế nào để tính toán tốc độ dòng phát của băng nhỏ giọt ở bất kỳ áp suất nào?

Bộ phát băng nhỏ giọt không bù áp suất-áp suất{1}} tuân theo phương trình phóng điện của bộ phát:q = k × H^x

Biểu tượng Nghĩa Đơn vị
q Tốc độ dòng phát L/h
k Hệ số phóng điện (xác định bằng hình học bộ phát) -
H Đầu áp suất vận hành m nước
x Số mũ phát (chỉ báo chế độ dòng chảy) -

số mũxcho bạn biết dòng chảy nhạy cảm như thế nào với sự thay đổi áp suất:

giá trị x Chế độ dòng chảy Nó có nghĩa là gì
0.0–0.2 Bù áp- Dòng chảy hầu như không thay đổi theo áp suất
0.4–0.6 Hỗn loạn (hầu hết các băng nhỏ giọt) Dòng chảy thay đổi khoảng √H
0.7–1.0 Đường dẫn nhiều lớp hoặc đường dẫn dài Dòng chảy rất nhạy cảm với áp suất-

 

Bộ phát-phẳng nhất vàbăng nhỏ giọt mê cungrơi vào vùng hỗn loạn vớix ≈ 0.47–0.57. Một nghiên cứu về sáu băng nhỏ giọt thương mại cho thấy giá trị x trung bình là 0,486.Vì mục đích ước tính khi giá trị k và x của nhà sản xuất không được công bố,x = 0.5là giá trị mặc định hợp lý cho-bộ phát băng nhỏ giọt dòng chảy hỗn loạn và k có thể được tính toán ngược-từ tốc độ dòng danh nghĩa ở áp suất định mức.

 

Ví dụ hoạt động 1: Lưu lượng giảm bao nhiêu ở áp suất thấp hơn?

Được cho:

băng nhỏ giọt phát phẳng, tốc độ dòng danh nghĩa: 1,38 L/h ở 0,1 MPa (đầu ≈10,2 m)

Giả sử x=0.5 (bộ phát hỗn loạn)

Bước 1 Quay lại-tính k:

k = q / H^x = 1.38 / 10.2^0.5 = 1.38 / 3.194 = 0.432

Bước 2 Tính toán lưu lượng ở áp suất 0,06 MPa (cột áp ≈6,1 m):

q = 0.432 × 6.1^0.5 = 0.432 × 2.470 = 1.07 L/h

Đó là mộtgiảm 22%từ mức 1,38 L/h định mức - chỉ khi vận hành ở 60% áp suất định mức.

 

Bạn sẽ nhận được tốc độ dòng chảy nào ở các áp suất khác nhau?

Sử dụng cùng một phương pháp tính toán ngược{0}}(x=0.5) cho thông số kỹ thuật của băng nhỏ giọt bộ phát phẳng:

Dòng danh nghĩa @

đầu 10m

k (ước tính) đầu 4m đầu 6m đầu 8m đầu 10m đầu 12m đầu 15m
0.8 L/h 0.253 0.51 0.62 0.71 0.80 0.88 0.98
1.1 L/h 0.348 0.70 0.85 0.98 1.10 1.20 1.35
1.38 L/h 0.436 0.87 1.07 1.23 1.38 1.51 1.69
2.0 L/h 0.632 1.26 1.55 1.79 2.00 2.19 2.45
3.0 L/h 0.949 1.90 2.32 2.68 3.00 3.29 3.67

Lưu ý: giá trị k được ước tính từ thông số kỹ thuật danh nghĩa giả sử x=0.5. Giá trị thực tế có thể khác nhau từ ±5–10% tùy thuộc vào hình dạng bộ phát. Khi có sẵn, hãy luôn sử dụng hệ số k và x do nhà sản xuất-công bố.

info-1297-864

 

Ⅱ. Làm thế nào để tính toán tổn thất ma sát trong các mặt bên của băng nhỏ giọt?

Phương trình Darcy-Weisbach là tiêu chuẩn để tính toán tổn thất áp lực ma sát trong đường ống:hf=f × (L/D) × (v²/2g)

Biểu tượng Nghĩa Đơn vị
hf Mất đầu ma sát m
f Darcy-hệ số ma sát Weisbach không thứ nguyên
L Chiều dài ống m
D Đường kính trong m
v Vận tốc dòng chảy m/s
g Gia tốc trọng trường (9,81) m/s²

 

⒈ Hệ số ma sát cho-Băng dính có thành mỏng

Để có dòng chảy rối trơn tru trong các ống nhựa có đường kính-nhỏ (4.000 < Re < 100.000), hệ số ma sát được tính bằng phương trình loại Blasius-:f=a / Re^0,25

trong đó Re=vD/υ (số Reynold) và υ=độ nhớt động học của nước (1,01 × 10⁻⁶ m2/s ở 20 độ ).

Hệ số a phụ thuộc vào loại ống:

Loại ống/băng một giá trị Nguồn
Ống cứng trơn tiêu chuẩn 0.3164 Blasius (bản gốc)
Ống PE có đường kính{0}}nhỏ (12–25mm) 0.300–0.302 Bagarello và cộng sự; Frizzone và cộng sự.
Turbo Tape (mê cung liên tục) 0.3442 Reti và cộng sự.
Băng nhỏ giọt bạc (mê cung liên tục) 0.3225 Reti và cộng sự.
Băng nhỏ giọt-bộ phát phẳng (ước tính) 0.32–0.34 Dự toán kỹ thuật

Mê cung liên tục được hàn bên trong băng nhỏ giọt có thành{0}}mỏng làm tăng ma sát vượt xa những gì mà công thức đường ống trơn-dự đoán. Bạn nên sử dụng=0.33 làm giá trị trung bình thận trọng cho băng nhỏ giọt-bộ phát phẳng khi không có dữ liệu thử nghiệm cụ thể.

 

⒉ Christiansen F{0}}Hệ số cho nhiều cửa hàng

Băng nhỏ giọt ở bên không phải là một ống đơn giản, nó có các đầu ra cách đều nhau để xả dòng chảy dọc theo chiều dài. Hệ số khử Christiansen giải thích cho điều này:hf_actual=F × hf_full_flow.Đối với bất kỳ bên nào có hơn ~ 20 bộ phát, F ≈ 0,35 là giá trị an toàn.

Số lượng ổ cắm (N) F
1 0.500
5 0.381
10 0.364
20 0.352
50 0.350
100+ 0.350

 

Ví dụ thực hiện 2: Tính toán tổn thất ma sát đầy đủ

Được cho:

Băng nhỏ giọt phát phẳng 16mm

Độ dày của tường: 0,2mm; đường kính trong ước tính: 15,6mm (0,0156 m)

Tốc độ dòng phát: 1,38 L/h ở đầu 10m

Khoảng cách phát: 30cm (0,3m)

Chiều dài bên: 150m

Áp suất đầu vào: 0,1 MPa (đầu 10,2 m)

Địa hình: bằng phẳng (độ dốc 0%)

Nhiệt độ nước: 20 độ

Bước 1: Tổng số nguồn phát:

N = 150 / 0.3 = 500 máy phát

Bước 2: Tổng tốc độ dòng chảy ngang (giả sử tất cả các bộ phát ở lưu lượng định mức):

Tổng Q_= 500 × 1.38=690 L/h =0.000192 m³/s

Trong thực tế, tốc độ dòng chảy giảm dọc theo chiều ngang khi áp suất giảm. Việc sử dụng tốc độ dòng chảy vào là biện pháp thận trọng và tiêu chuẩn cho thiết kế ban đầu.

Bước 3: Vận tốc dòng chảy ở cửa vào:

v = 4Q / (πD²) = 4 × 0.000192 / (π × 0.0156²) = 1.00 m/s

Bước 4: Số Reynolds:

Re=vD/υ=1.00 × 0,0156 / (1,01 × 10⁻⁶) =15,446

Đây là phạm vi hỗn loạn trơn tru (4.000 < Re < 100.000), do đó áp dụng phương trình Blasius đã sửa đổi.

Bước 5: Hệ số ma sát (a=0.33 đối với băng phát- phẳng):

f = 0.33 / 15446^0.25 = 0.33 / 11.16 = 0.0296

Bước 6: Tổn thất ma sát dòng chảy-toàn bộ (không hiệu chỉnh đầu ra):

hf_raw=0.0296 × (150 / 0,0156) × (1,00² / 19,62)=0.0296 × 9615 × 0.0510 =14.50 m

Bước 7: Áp dụng hệ số F{1}} Christiansen (N=500, F=0.35):

hf_actual=0.35 × 14.50 =5.08 m0,050 MPa

Bước 8: Áp lực ở phần đuôi:

P_đuôi=10.2 - 5.08 =5.12 m0,050 MPa

Phán quyết:Áp suất-cuối đuôi 0,050 MPa vừa phải với áp suất vận hành tối thiểu được khuyến nghị dành cho băng nhỏ giọt bộ phát phẳng (0,05 MPa) [3]. Ở độ cao 150m, mặt bên này đạt giới hạn thiết kế. Bất kỳ tổn thất bổ sung nào từ phụ kiện, bộ lọc hoặc độ cao sẽ đẩy phần đuôi xuống dưới thông số kỹ thuật.

Điều gì thay đổi ở 120m?Chạy tính toán tương tự cho 120m:

  • N = 400, Q = 0.000154 m³/s
  • hf_actual =3.25 m(0,032 MPa)
  • P_tail=10.2 - 3.25=6.95 m (0,068 MPa) → lề thoải mái

 

Ⅲ. Khi nào nên sử dụng Hazen-Williams để giảm ma sát trên băng nhỏ giọt?

Phương trình Hazen-Williams đơn giản hơn và được sử dụng rộng rãi trong thiết kế tưới tiêu:hf=10.67 × L × Q^1,852 / (C^1,852 × D^4,87)

Biểu tượng Nghĩa Đơn vị
hf Mất đầu m
L Chiều dài ống m
Q Tốc độ dòng chảy L/s
C Hazen-Hệ số độ nhám Williams không thứ nguyên
D Đường kính trong m

Đối với băng nhỏ giọt polyetylen, giá trị C trong tài liệu nằm trong khoảng từ 130 đến 150. Tiện ích mở rộng UF/IFAS sử dụng C=130 cho các đường đa cạnh ¾-inch trong tính toán tưới nhỏ giọt.

 

Darcy-Weisbach vs. Hazen-Williams: Công thức giảm ma sát nào chính xác hơn cho băng nhỏ giọt?

Sử dụng các thông số tương tự như Ví dụ 2 (băng 16mm, 1,38 L/h, khoảng cách 30cm, 150m, lưu lượng vào 0,192 L/s, D=0.0156 m):

Phương pháp hf (m) Áp suất đuôi (MPa) nhà tdentle amtdle
Darcy-Weisbach (một=0.33) 5.08 0.050 Đường cơ sở
Hazen-Williams (C=150) 4.35 0.057 −14,4% (đánh giá thấp)
Hazen-Williams (C=140) 4.80 0.053 −5.5%
Hazen-Williams (C=130) 5.36 0.047 +5.5% (đánh giá quá cao)

Mua mang về:Hazen-Williams với C=140–145 gần đúng với kết quả Darcy-Weisbach trong phạm vi ±5% cho kịch bản này. C=150 quá lạc quan. C=130 đưa ra ước tính thận trọng. Đối với thiết kế cuối cùng, hãy luôn xác minh với Darcy{10}}Weisbach bằng hệ số Blasius đã sửa đổi.

 

Ⅳ. Bạn có thể chạy băng nhỏ giọt bên trong bao lâu?

Hạn chế thiết kế chính cho các mặt bên của băng nhỏ giọt làsự thay đổi dòng chảy- chênh lệch giữa tốc độ dòng phát cao nhất và thấp nhất trên một bên không được vượt quá 10% (theo ISO và tiêu chuẩn quốc gia Trung Quốc GB/T 50485).

  • Đối với các nguồn phát hỗn loạn có x ≈ 0,5, sự thay đổi dòng chảy 10% tương ứng với sự thay đổi áp suất khoảng 20% ​​(vì Δq/q ≈ x × ΔH/H). Điều này có nghĩa là:Biến đổi áp suất cho phép=±10% đầu vào
  • Đối với địa hình bằng phẳng, toàn bộ sự thay đổi áp suất xuất phát từ tổn thất ma sát, vì vậy:hf_allowable ≈ 0,20 × H_inlet

 

Chiều dài chạy băng nhỏ giọt tối đa theo tốc độ dòng chảy và khoảng cách là bao nhiêu?

Bảng sau đây trình bày khoảng thời gian chạy tối đa ước tính cho băng nhỏ giọt bộ phát phẳng trên địa hình bằng phẳng, giả sử biến đổi dòng chảy 10% (biến thiên áp suất 20%) và áp suất đầu vào ở cột nước 10m. Được tính toán bằng cách sử dụng Darcy-Weisbach với=0.33 và Christiansen F=0.35.

Băng 16mm (ID ước tính: 15,6mm):

Dòng phát Khoảng cách Độ dài chạy tối đa Số lượng máy phát
0.8 L/h 20cm 254m 1270
0.8 L/h 30cm 327m 1090
1.38 L/h 20cm 135m 675
1.38 L/h 30cm 174m 580
2.0 L/h 20cm 93m 465
2.0 L/h 30cm 120m 400

Băng 22mm (ID ước tính: 21,4mm):

Dòng phát Khoảng cách Độ dài chạy tối đa Số lượng máy phát
0.8 L/h 20cm 468m 2340
0.8 L/h 30cm 603m 2010
1.38 L/h 20cm 249m 1245
1.38 L/h 30cm 321m 1070
2.0 L/h 20cm 171m 855
2.0 L/h 30cm 220m 733

Xác minh:Các giá trị này nhất quán với dữ liệu về thời lượng chạy tối đa-do nhà sản xuất công bố cho các sản phẩm băng nhỏ giọt tương đương. Ví dụ: Băng nhỏ giọt Dripmax Silver (khoảng cách 16 mm, 0,4 L/h, 30 cm) liệt kê 371m ở mức biến thiên lưu lượng 10% và đầu vào 1,0 bar. Giá trị tính toán của chúng tôi cho tốc độ dòng chảy thấp hơn (0,8 L/h so với{10}} L/h) ở cùng đường kính sẽ ngắn hơn, điều này được mong đợi vì tốc độ dòng chảy cao hơn trên mỗi bộ phát với khoảng cách gần hơn sẽ tạo ra nhiều ma sát hơn.

Ghi chú:Tất cả các giá trị giả định địa hình bằng phẳng. Xem phần tiếp theo để điều chỉnh độ dốc.

Ⅴ. Độ dốc ảnh hưởng đến áp suất băng nhỏ giọt như thế nào?

Những thay đổi về độ cao sẽ cộng hoặc trừ áp suất sẵn có tại mỗi điểm dọc theo mặt bên:ΔH_độ cao=± Δz

trong đó Δz là sự thay đổi độ cao (dương khi lên dốc, âm khi xuống dốc). Sự thay đổi áp suất tính bằng MPa trên 10m thay đổi độ cao là:ΔP=0.098 MPa trên độ cao 10m

Hoặc tương đương:Độ cao 1m=0.0098 MPa=0.1 bar ≈ 1,42 PSI

 

Tác động thực tế đến thời lượng chạy

Độ dốc Thay đổi áp suất trên mỗi chiều dài 100m Ảnh hưởng đến thời lượng chạy tối đa
Lên dốc 0,5% −0,0049 MPa Giảm độ dài tối đa xuống ~15–20%
Lên dốc 1% −0,0098 MPa Giảm độ dài tối đa xuống ~30–40%
Phẳng 0 Sử dụng độ dài tối đa được tính toán
Xuống dốc 0,5% +0.0049 MPa Tăng độ dài tối đa thêm ~15–20%
Xuống dốc 1% +0.0098 MPa Tăng độ dài tối đa lên ~30–40%

 

Ví dụ hoạt động 3: Điều gì xảy ra với áp lực lên độ dốc 1%?

Được cho:Băng 16mm, bộ phát 1,38 L/h, khoảng cách 30cm, chiều ngang 150m, đầu vào 0,1 MPa

Tình trạng Mất ma sát Thay đổi độ cao Thay đổi áp suất ròng Áp lực đuôi Phán quyết
Phẳng 0,050 MPa 0 −0,050 MPa 0,050 MPa Ở giới hạn
Lên dốc 1% 0,050 MPa +0.015 MPa −0,065 MPa 0,035 MPa Thất bại
Xuống dốc 1% 0,050 MPa −0,015 MPa −0,035 MPa 0,065 MPa Vượt qua với lề

Trên độ dốc lên 1%, cùng độ dốc 150m đó giảm xuống 0,035 MPa ở phần đuôi - thấp hơn nhiều so với mức tối thiểu 0,05 MPa. Bạn cần rút ngắn chiều ngang xuống còn khoảng 100m hoặc chuyển sang băng 22mm.

Khi xuống dốc 1%, độ cao tăng lên sẽ bù đắp một phần tổn thất ma sát và áp suất ở đuôi là 0,065 MPa. Bạn có thể mở rộng phạm vi bên này lên khoảng 200m trước khi đạt đến giới hạn áp suất.

 

Câu hỏi thường gặp: 5 sai lầm thường gặp trong thiết kế thủy lực băng nhỏ giọt

Tại sao bạn không nên tin tưởng vào tốc độ dòng danh nghĩa trên bảng thông số kỹ thuật

+

-

Lưu lượng định mức trên bảng thông số kỹ thuật của sản phẩm áp dụng ở đúng một áp suất. Bộ phát 1,38 L/h ở cột nước 10m chỉ cung cấp 1,07 L/h ở cột nước 6m. Nếu thiết kế của bạn giả định 1,38 L/h ở mọi nơi, bạn sẽ đánh giá quá cao lượng nước cung cấp lên tới 22% ở phần cuối.

Sửa chữa:Luôn tính toán lưu lượng thực tế ở áp suất cuối-đuôi bằng cách sử dụng q=k × H^x.

Điều gì xảy ra khi bạn bỏ qua yếu tố F{0}} Christiansen

+

-

Một chiều ngang 150m với 500 bộ phát có tổn thất ma sát chỉ bằng 35% so với những gì bạn tính toán với giả định toàn bộ dòng chảy dọc theo toàn bộ chiều dài. Việc bỏ qua hệ số F-sẽ đánh giá quá cao tổn thất do ma sát ~3×, điều này có thể dẫn đến việc bạn sử dụng các đường ống quá khổ một cách không cần thiết - hoặc tệ hơn, mang lại cho bạn cảm giác tự tin sai lầm vì bạn cho rằng tổn thất là rất lớn và bạn "đã tính đến điều đó".

Sửa chữa:Áp dụng F=0.35 cho bất kỳ bên nào có hơn 20 bộ phát.

Tại sao Hệ số Blasius Tiêu chuẩn (a=0.3164) lại sai đối với băng nhỏ giọt

+

-

Băng nhỏ giọt có thành mỏng-có mê cung liên tục hoặc bộ phát phẳng có độ ma sát cao hơn ống cứng trơn. Nghiên cứu đã công bố cho thấy=0.3225–0,3442 đối với băng có thành-mỏng có các mê cung liên tục [1]. Sử dụng 0,3164 đánh giá thấp tổn thất ma sát từ 2–8%.

Sửa chữa:Sử dụng=0.33 cho băng nhỏ giọt-bộ phát phẳng khi không có dữ liệu thử nghiệm cụ thể.

Tại sao chỉ riêng tổn hao ma sát không nói lên toàn bộ câu chuyện

+

-

Ma sát chỉ là một thành phần của sự thay đổi áp suất dọc theo một bên. Độ cao có thể cộng hoặc trừ nhiều như vậy. Trên địa hình đồi núi, việc bỏ qua độ cao có thể dẫn đến lỗi hệ thống ở các điểm cao hoặc ngập lụt ở các điểm thấp.

Sửa chữa:Tổng thay đổi áp suất=tổn thất ma sát ± thay đổi độ cao. Luôn bao gồm cả hai.

Tại sao C=150 lại quá lạc quan đối với băng nhỏ giọt

+

-

C=150 thích hợp cho nguồn điện PVC mới, trơn. Điều này quá lạc quan đối với các mặt bên của băng nhỏ giọt, có bộ phát bên trong và các mặt cắt ngang (trong trường hợp-băng có thành mỏng)-biến dạng dưới áp lực. Việc sử dụng C=150 đánh giá thấp tổn thất do ma sát từ 10–15% so với Darcy-Weisbach với hệ số Blasius đã được hiệu chỉnh.

Sửa chữa:Sử dụng C=130 để ước tính H-W thận trọng hoặc tốt hơn nữa là sử dụng Darcy-Weisbach.

 

Tham khảo nhanh: Tóm tắt các công thức chính

Những gì bạn cần Công thức Thông số chính
Dòng phát ở mọi áp suất q = k × H^x k từ thông số kỹ thuật danh nghĩa; x ≈ 0,5 đối với nguồn phát nhiễu loạn
Tổn thất do ma sát (Darcy-Weisbach) hf=f × (L/D) × (v²/2g) × F f=a/Re^0,25; a ≈ 0,33; F ≈ 0,35
Tổn hao do ma sát (Hazen-Williams) hf=10.67 × L × Q^1,852 / (C^1,852 × D^4,87) C=130–140 cho băng nhỏ giọt
Thay đổi áp suất độ cao ΔP=±0,0098 MPa trên 1m độ cao +lên dốc, −xuống dốc
Ma sát cho phép đối với sự thay đổi dòng chảy 10% hf_allowable ≈ 0,20 × H_inlet Giả sử x ≈ 0,5

 

Tài liệu tham khảo

1. Reti, C. và cộng sự. "Mất đầu trong-Băng nhỏ giọt có tường mỏng với mê cung liên tục."Nước(MDPI), 2019. PMC6925943

2. "内镶贴片式滴头流道结构参数对水力性能影响的试验研究." 节水灌溉, 2023. Liên kết

3. Zazueta, FS "Những cân nhắc về thủy lực cho hệ thống tưới vi sinh cho cây có múi." Phần mở rộng UF/IFAS, Ấn bản CH156. liên kết

4.Dữ liệu kỹ thuật băng nhỏ giọt Dripmax Silver. liên kết

5.Rivulis T-Đặt tên sản phẩm băng nhỏ giọt và tính toán lưu lượng. liên kết

6.Bagarello, V. và cộng sự. "Nghiên cứu thực nghiệm về khả năng cản dòng chảy trong ống nhựa có đường kính{2}}nhỏ."Tạp chí Kỹ thuật Thủy lợi và Thoát nước, 1997.