Reynolds và miền chuyển tiếp: khi nào hệ số ma sát trở nên bất định

Một con số Reynolds có thể làm thay đổi cả cách tính
Một tuyến dầu giả định cho kết quả tổn thất nhỏ khi dùng công thức dòng tầng. Sau khi thay độ nhớt bằng giá trị ở nhiệt độ khác, bảng tính vẫn tự động áp dụng công thức cũ. Kết quả trông hợp lý vì độ nhớt giảm thì áp suất giảm, nhưng có thể người tính đã đi qua miền mà giả thiết dòng tầng không còn được bảo đảm. Bài toán vì vậy bắt đầu bằng kiểm tra chế độ chảy, trước khi tối ưu từng chữ số của áp suất.

Reynolds biểu thị tỷ lệ đặc trưng giữa tác động quán tính và nhớt. NASA trình bày vai trò của nó như một tham số tương tự của dòng chảy. [1] Trong bài này, chiều dài đặc trưng là đường kính trong của ống tròn; vận tốc là trung bình trên tiết diện. Không dùng đường kính danh nghĩa hoặc vận tốc ở tâm dòng thay cho hai đại lượng đó.

Công thức, đơn vị và các giả thiết
Với một chất lỏng Newton, một pha, có khối lượng riêng và độ nhớt không đổi trên đoạn xét:

v = 4Q/(πD²)
Re = ρvD/μ = vD/ν = 4ρQ/(πμD)

Q là lưu lượng thể tích tại điều kiện trong ống, m³/s; D là đường kính trong, m; v là vận tốc trung bình, m/s; ρ là khối lượng riêng, kg/m³; μ là độ nhớt động lực, Pa·s; ν = μ/ρ là độ nhớt động học, m²/s. Re không có thứ nguyên. Đơn vị mPa·s phải đổi sang Pa·s bằng hệ số 0,001; cSt phải đổi sang m²/s bằng hệ số 10⁻⁶.

Hai cách ghi độ nhớt không thể hoán đổi trực tiếp. Nếu một phiếu tính chất cho ν = 20 cSt, điều đó không có nghĩa μ = 0,020 Pa·s cho mọi chất lỏng. Với ρ = 900 kg/m³, μ tương ứng là 0,018 Pa·s. Sai khác nhỏ trong thao tác này có thể đáng chú ý khi trạng thái đang gần miền chuyển tiếp.

Đối với dòng tầng phát triển đầy đủ trong ống tròn, mô hình của MathWorks sử dụng hằng số ma sát dòng tầng 64. [2] Viết theo hệ số Darcy, fD = 64/Re và Δp = fD(L/D)ρv²/2. Thế Re vào rồi rút gọn thu được Δp = 32μLv/D² = 128μLQ/(πD⁴). L là chiều dài đoạn thẳng, m; Δp tính bằng Pa. Bài viết chỉ xét ma sát đoạn thẳng nằm ngang, chưa cộng tổn thất phụ kiện hoặc chênh cao.

Ví dụ cơ sở: kiểm tra bằng hai đường tính
Chọn Q = 0,001 m³/s, tương đương 3,6 m³/h; D = 0,050 m; L = 120 m; ρ = 900 kg/m³ và μ = 0,020 Pa·s. Đây là bộ số giả định, không phải dữ liệu của một sản phẩm dầu cụ thể. Diện tích trong ống bằng 0,00196350 m², nên v = 0,509296 m/s.

Re bằng 1.145,92 và fD = 0,0558505. Nhánh tính Darcy cho Δp = 0,0558505 × (120/0,050) × 900 × 0,509296²/2 = 15.645,57 Pa, hay 15,65 kPa. Nhánh tính đã rút gọn cho 128 × 0,020 × 120 × 0,001/(π × 0,050⁴), cũng bằng 15.645,57 Pa.

Sự khớp nhau kiểm tra phép thế và đơn vị, không xác nhận mọi giả thiết của hiện trường. Nếu một phần của tuyến là đoạn vào chưa phát triển, hoặc có thay đổi tiết diện đáng kể, mô hình đoạn thẳng phát triển đầy đủ cần được xem lại. Cũng không dùng sự khớp của hai công thức tương đương như hai bằng chứng vật lý độc lập: cả hai cùng dựa trên một mô hình dòng tầng.

Đáng chú ý, ρ biến mất khỏi công thức Δp dòng tầng sau rút gọn khi Q, μ và hình học cố định. Điều này không nói khối lượng riêng luôn không quan trọng. Nó vẫn tham gia Re, trọng lượng cột chất lỏng và việc đổi giữa μ với ν. Nếu giữ ν thay vì giữ μ, thay đổi ρ cũng làm thay đổi μ và tổn thất tính được.

Giảm độ nhớt: khi phép ngoại suy cần dừng lại
Giữ nguyên Q và hình học, giảm μ xuống 0,010 Pa·s làm Re tăng gấp đôi, thành 2.291,83. Nếu vẫn áp dụng nhánh dòng tầng, fD bằng 0,0279253 và Δp bằng 7.822,78 Pa. Đây là kết quả có điều kiện của nhánh công thức, chưa phải một dự đoán đã được xác nhận cho dòng thực ở Reynolds này.

Tài liệu COMSOL phân biệt các tương quan dành cho dòng tầng, rối và cách nối giữa chúng; đồng thời lưu ý một cách nối công thức không nhất thiết cho ước lượng chính xác trong vùng chuyển tiếp. [3] Vì vậy, không nên coi một đường cong trơn do phần mềm vẽ ra là bằng chứng rằng trạng thái vật lý đã được xác định duy nhất.

Để lượng hóa ảnh hưởng của lựa chọn mô hình, thử hai hệ số Darcy giả định fD = 0,028 và 0,045 tại cùng vận tốc. Chúng cho Δp lần lượt 7,84 và 12,61 kPa, chênh 4,76 kPa. Hai hệ số này chỉ là hai kịch bản khảo sát, không phải cận bảo đảm hoặc khuyến nghị lựa chọn. Kết quả giúp đặt câu hỏi thực tế: phần bất định do mô hình có lớn hơn phần thay đổi đang cố tối ưu hay không?

Nếu giảm tiếp μ xuống 0,005 Pa·s, Re lên 4.583,66. Công thức 64/Re vẫn trả về một con số, nhưng việc tính được không làm nó phù hợp. Cần chọn tương quan với phạm vi áp dụng và độ nhám phù hợp, rồi tính lại. Đối với chất lưu phi Newton, mô hình độ nhớt còn phụ thuộc quan hệ ứng suất và tốc độ biến dạng; lấy một giá trị μ duy nhất rồi dùng mọi công thức Newton là một giả thiết bổ sung đáng kiểm tra. [3]

Độ nhạy không chỉ nằm ở lưu lượng
Ở trạng thái μ = 0,010 Pa·s, sai khác độ nhớt ±10% cho Re khoảng 2.083,48 đến 2.546,48. Khoảng này không đối xứng quanh 2.291,83 vì Re tỷ lệ nghịch với μ. Khi dữ liệu độ nhớt thiếu chắc chắn, in Re đến hai chữ số thập phân không phản ánh mức hiểu biết thật. Nên giữ khoảng đầu vào cùng các trường hợp tính thay vì chỉ một giá trị trung tâm.

Quay lại trạng thái cơ sở dòng tầng với μ = 0,020 Pa·s, giảm đường kính 10% từ 50 xuống 45 mm làm Δp tăng theo (50/45)⁴, đạt 23,85 kPa, tăng 52,42%. Tăng đường kính lên 55 mm cho 10,69 kPa, giảm 31,70%. Re cũng thay đổi theo nghịch đảo D ở Q cố định, nên sau mỗi thay đổi đường kính vẫn phải kiểm tra chế độ chảy.

Nếu thay đổi Q trong miền dòng tầng và các tính chất không đổi, Δp tỷ lệ với Q. Giữ fD cố định để áp dụng luật bình phương lưu lượng cho trường hợp này sẽ tạo sai mô hình, bởi chính fD giảm khi Re tăng. Mối liên hệ giữa Q, Re và fD phải được cập nhật cùng nhau; không chọn riêng quy luật nào cho kết quả thuận tiện nhất.

Một hồ sơ tính có thể kiểm tra lại
Đầu vào nên ghi nguồn và nhiệt độ của μ, ρ; đường kính trong thực dùng; lưu lượng tại điều kiện nào; chiều dài và phạm vi đoạn thẳng. Đầu ra nên có v, Re, tương quan ma sát, phạm vi áp dụng và Δp. Nếu khảo sát chuyển tiếp, cần phân biệt dữ liệu đo, lựa chọn mô hình và các hệ số giả định dùng để thử độ nhạy.

Giá trị cần giữ lại từ ví dụ là cách ra quyết định: kiểm tra chế độ trước, tính tổn thất trong phạm vi phù hợp, rồi xem bất định nào đang chi phối. Cách đó giúp người đọc biết khi nào một kết quả đủ dùng cho so sánh sơ bộ và khi nào cần tính chất chất lưu, khảo sát hình học hoặc mô hình chi tiết hơn.

Nguồn tham khảo
[1] Similarity Parameters — NASA Glenn Research Center. Truy cập ngày 09/11/2026.
[2] Pipe (TL) — Pipe that transports fluid in thermal liquid networks — MathWorks. Truy cập ngày 09/11/2026.
[3] Flow Equations — Pipe Flow Module — COMSOL. Truy cập ngày 09/11/2026.

Rà soát nội dung: 09/11/2026.
 
Top