Lấy thể tích chia lưu lượng đầu có thể thiếu thời gian
Khi một bồn hở xả nhờ trọng lực, cột áp giảm cùng mức chất lỏng. Một phép tính dùng lưu lượng lúc đầu cho cả quá trình thường cho thời gian ngắn hơn mô hình biến thiên. Tuy nhiên, thay nó bằng một công thức “xả cạn” cũng chưa đủ nếu phía sau lỗ còn đường ống hoặc có dòng cấp vào bồn.
Bài này xét chất lỏng không nén được, nhiệt độ không đổi, bồn thông khí đủ và đầu ra là tia tự do ra cùng áp suất khí quyển. Diện tích mặt thoáng lớn hơn nhiều diện tích lỗ, nên động năng của mặt thoáng được bỏ qua. Không có bơm hoặc đường ống dài phía sau lỗ trong trường hợp cơ sở. Đây là mô hình tính thời gian, không phải hướng dẫn mở xả một thiết bị đang chứa lưu chất công nghệ.
Ghép lưu lượng tức thời với cân bằng thể tích
Đặt h là độ cao mặt chất lỏng trên tâm lỗ, m; a là diện tích hình học lỗ, m²; C_d là hệ số lưu lượng không thứ nguyên. Quan hệ dùng trong ví dụ:
Q_out=C_d a√(2gh)=k√h; k=C_d a√(2g).
Q_out có đơn vị m³/s, k có đơn vị m^(5/2)/s và g=9,81 m/s². Virtual Labs phân biệt hệ số lưu lượng với hệ số vận tốc: tia có thể vừa co hẹp vừa giảm vận tốc so với trường hợp lý tưởng. [1] Vì vậy, C_d không chỉ là hệ số giảm tốc độ tại một tiết diện tia bất kỳ.
Với tiết diện ngang bồn A_t không đổi và không có dòng vào, cân bằng là A_t dh/dt=−k√h. Tích phân giữa h_i và h_f cho:
t_f=2A_t(√h_i−√h_f)/k.
√h(t)=√h_i−kt/(2A_t).
Mức h không giảm tuyến tính theo thời gian, nhưng căn mức giảm tuyến tính trong mô hình này. Công thức chỉ dùng trong khoảng h còn thỏa giả thiết lỗ ngập và C_d có thể coi không đổi. Không tiếp tục bình phương vế phải sau khi nó đã âm để tạo ra một mực nước tăng lại.
Một bồn trụ và lỗ tròn giả định
Cho đường kính trong bồn 2,00 m, A_t=π m²; đường kính lỗ d=0,050 m, a=0,001963495 m²; C_d=0,620. Mức ban đầu h_i=2,40 m, mức kết thúc h_f=0,20 m. Các kích thước, hệ số và mức đều là dữ liệu giả định, không lấy từ một thiết bị thực.
Tính được k=0,005392263 m^(5/2)/s và thời gian hạ mức t_f=1.284,051 s, tương đương 21,4009 phút. Lưu lượng đầu là 0,008353658 m³/s; cuối là 0,002411493 m³/s. Thể tích đã xả bằng A_t(h_i−h_f)=6,911504 m³.
Nếu dùng lưu lượng đầu để chia thể tích, thời gian chỉ là 827,363 s. Sai khác không đến từ hệ số tổn thất khác nhau: cả hai phép dùng cùng C_d. Nó đến từ việc một phép giữ Q ở giá trị cao nhất, còn phép kia cập nhật cột áp liên tục. Trong đúng mô hình bồn trụ có diện tích và C_d không đổi, không có dòng vào, căn bậc hai của h giảm tuyến tính theo thời gian. Vì Q tỷ lệ với căn bậc hai đó, Q cũng giảm tuyến tính; lưu lượng trung bình theo thời gian bằng chính (Q_đầu + Q_cuối)/2, khoảng 0,005382575 m³/s. Đẳng thức riêng này không còn được bảo đảm khi đổi hình học bồn hoặc thêm dòng vào.
Một phép kiểm thuận tiện là thay thời gian vừa tìm vào quan hệ căn mức. Giá trị trở về √0,20, xác nhận rằng cùng mốc thời gian và cùng cao độ đã được sử dụng. Nếu lấy h đo từ đáy bồn trong khi tâm lỗ cao hơn đáy, sai lệch mốc sẽ xuất hiện ngay ở phép kiểm này.
Độ nhạy: đường kính và C_d có vai trò khác nhau
Giữ mọi đầu vào khác, t tỷ lệ nghịch C_d và nghịch d². Nếu C_d thực chỉ 0,560, thời gian tăng thành 1.421,628 s; nếu là 0,680, thời gian giảm còn 1.170,753 s. Với lỗ giảm từ 50 xuống 40 mm và C_d vẫn 0,620, thời gian thành 2.006,330 s, tăng 56,25%.
Đây là các kịch bản thay đổi từng tham số riêng. Trong thiết bị thật, đổi đường kính hoặc hình dạng mép lỗ có thể làm C_d thay đổi theo. Không dùng quan hệ d⁻² để chứng minh một phương án mới chỉ cần đổi diện tích. Nghiên cứu của Joye và Barrett đối chiếu nhiều cấu hình bồn–ống và chỉ ra tầm quan trọng của việc mô tả đúng tổn thất cùng hình học. [2]
Nếu một thử nghiệm giả định với chính bồn và lỗ này mất 900 s để xuống từ 2,00 tới 0,80 m, đảo công thức cho C_d=0,417237. Con số ấy là hệ số hiệu dụng phù hợp với cặp mức và thời gian đã đo. Nó chưa chứng minh mép lỗ có một tính chất cố định như vậy: sai số mức, một đoạn ống bị bỏ sót hoặc thông khí kém đều có thể bị hấp thụ vào tham số hiệu dụng.
Khi dòng cấp vẫn còn
Xét độc lập trường hợp một dòng cấp không đổi Q_in=0,001000 m³/s tiếp tục vào bồn. Khi đó A_t dh/dt=Q_in−k√h. Mức cân bằng theo mô hình bằng h_eq=(Q_in/k)²=0,0343920 m. Từ mức cao hơn, bồn tiến về mức ấy; nó không thể tự hạ tới 0,020 m khi dòng cấp giữ nguyên.
Đặt b=Q_in/k, có đơn vị √m. Với h_i>h_f>h_eq, tích phân cho:
t_f=(2A_t/k){√h_i−√h_f+b ln[(√h_i−b)/(√h_f−b)]}.
Tỷ số trong logarithm không có đơn vị. Với hai mức 2,40 và 0,20 m, thời gian thành 1.640,721 s. Khi h_f tiến sát h_eq từ phía trên, mẫu bên trong logarithm tiến về không và thời gian tăng không giới hạn. Đây là khác biệt vật lý giữa “còn một ít dòng vào” với việc chỉ cộng một lượng thể tích cố định vào lúc đầu.
Công thức này không còn dùng được nếu Q_in thay đổi theo thời gian hoặc theo mức. Khi đó cần giải phương trình vi phân với lịch sử dòng vào. Có thể kiểm lời giải số bằng cân bằng tích lũy: thể tích giảm bằng tích phân dòng ra trừ tích phân dòng vào.
Hình dáng bồn thay đổi tích phân
Nếu tiết diện ngang phụ thuộc mức, phương trình tổng quát là A_t(h)dh/dt=−k√h. Với một bồn hình nón lý tưởng, đỉnh dưới tại h=0, bán kính miệng R=1,00 m và chiều cao H=2,40 m, A_t(h)=πR²h²/H². Bỏ qua hình học rất nhỏ quanh lỗ, thời gian giữa hai mức là:
t_f=2πR²(h_i^(5/2)−h_f^(5/2))/(5H²k).
Lấy h_i=2,40 và h_f=0,20 m, kết quả 360,307 s. Không diễn giải rằng đáy nón luôn xả nhanh hơn bồn trụ cùng lượng chứa: ở hai mức này bồn nón chỉ xả 2,511820 m³, ít hơn đáng kể 6,911504 m³ của ví dụ trụ. So sánh thời gian cần đi kèm thể tích và cơ sở hình học.
Giới hạn gần cuối quá trình
Khi mức xuống gần lỗ, cột áp tại tâm không còn đại diện tốt cho toàn diện tích nếu một phần lỗ mất ngập. Autodesk mô tả việc chuyển mô hình sang hành vi tràn khi mặt nước thấp hơn đỉnh lỗ. [3] Vì vậy, thời gian tới mức nhỏ nên được tách khỏi phần đường cong đã dùng công thức lỗ ngập.
Nếu đầu ra chìm trong một bể khác, phải dùng chênh cột áp hữu hiệu và xét sự đổi mức bên nhận. Nếu nối ống, cần mô hình ma sát, phụ kiện, cao độ đầu ống và chế độ dòng phù hợp. Bài học thực dụng là xác định đúng hai biên áp suất, mốc h và diện tích chứa trước khi chọn công thức. Khi dữ liệu đo được dùng để hiệu chỉnh C_d, hãy giữ lại đường cong mức–thời gian để kiểm xem một tham số duy nhất có giải thích được toàn khoảng xả hay không.
Nguồn tham khảo
[1] Orifices Experiment — Theory — Virtual Labs, International Institute of Information Technology Hyderabad. Truy cập ngày 09/11/2026.
[2] The Tank Drainage Problem Revisited: Do These Equations Actually Work? — Donald D. Joye và Branden C. Barrett; The Canadian Journal of Chemical Engineering 81, October 2003; bản lưu Đại học São Paulo. Truy cập ngày 09/11/2026.
[3] Orifice — InfoDrainage User Guide — Autodesk. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Khi một bồn hở xả nhờ trọng lực, cột áp giảm cùng mức chất lỏng. Một phép tính dùng lưu lượng lúc đầu cho cả quá trình thường cho thời gian ngắn hơn mô hình biến thiên. Tuy nhiên, thay nó bằng một công thức “xả cạn” cũng chưa đủ nếu phía sau lỗ còn đường ống hoặc có dòng cấp vào bồn.
Bài này xét chất lỏng không nén được, nhiệt độ không đổi, bồn thông khí đủ và đầu ra là tia tự do ra cùng áp suất khí quyển. Diện tích mặt thoáng lớn hơn nhiều diện tích lỗ, nên động năng của mặt thoáng được bỏ qua. Không có bơm hoặc đường ống dài phía sau lỗ trong trường hợp cơ sở. Đây là mô hình tính thời gian, không phải hướng dẫn mở xả một thiết bị đang chứa lưu chất công nghệ.
Ghép lưu lượng tức thời với cân bằng thể tích
Đặt h là độ cao mặt chất lỏng trên tâm lỗ, m; a là diện tích hình học lỗ, m²; C_d là hệ số lưu lượng không thứ nguyên. Quan hệ dùng trong ví dụ:
Q_out=C_d a√(2gh)=k√h; k=C_d a√(2g).
Q_out có đơn vị m³/s, k có đơn vị m^(5/2)/s và g=9,81 m/s². Virtual Labs phân biệt hệ số lưu lượng với hệ số vận tốc: tia có thể vừa co hẹp vừa giảm vận tốc so với trường hợp lý tưởng. [1] Vì vậy, C_d không chỉ là hệ số giảm tốc độ tại một tiết diện tia bất kỳ.
Với tiết diện ngang bồn A_t không đổi và không có dòng vào, cân bằng là A_t dh/dt=−k√h. Tích phân giữa h_i và h_f cho:
t_f=2A_t(√h_i−√h_f)/k.
√h(t)=√h_i−kt/(2A_t).
Mức h không giảm tuyến tính theo thời gian, nhưng căn mức giảm tuyến tính trong mô hình này. Công thức chỉ dùng trong khoảng h còn thỏa giả thiết lỗ ngập và C_d có thể coi không đổi. Không tiếp tục bình phương vế phải sau khi nó đã âm để tạo ra một mực nước tăng lại.
Một bồn trụ và lỗ tròn giả định
Cho đường kính trong bồn 2,00 m, A_t=π m²; đường kính lỗ d=0,050 m, a=0,001963495 m²; C_d=0,620. Mức ban đầu h_i=2,40 m, mức kết thúc h_f=0,20 m. Các kích thước, hệ số và mức đều là dữ liệu giả định, không lấy từ một thiết bị thực.
Tính được k=0,005392263 m^(5/2)/s và thời gian hạ mức t_f=1.284,051 s, tương đương 21,4009 phút. Lưu lượng đầu là 0,008353658 m³/s; cuối là 0,002411493 m³/s. Thể tích đã xả bằng A_t(h_i−h_f)=6,911504 m³.
Nếu dùng lưu lượng đầu để chia thể tích, thời gian chỉ là 827,363 s. Sai khác không đến từ hệ số tổn thất khác nhau: cả hai phép dùng cùng C_d. Nó đến từ việc một phép giữ Q ở giá trị cao nhất, còn phép kia cập nhật cột áp liên tục. Trong đúng mô hình bồn trụ có diện tích và C_d không đổi, không có dòng vào, căn bậc hai của h giảm tuyến tính theo thời gian. Vì Q tỷ lệ với căn bậc hai đó, Q cũng giảm tuyến tính; lưu lượng trung bình theo thời gian bằng chính (Q_đầu + Q_cuối)/2, khoảng 0,005382575 m³/s. Đẳng thức riêng này không còn được bảo đảm khi đổi hình học bồn hoặc thêm dòng vào.
Một phép kiểm thuận tiện là thay thời gian vừa tìm vào quan hệ căn mức. Giá trị trở về √0,20, xác nhận rằng cùng mốc thời gian và cùng cao độ đã được sử dụng. Nếu lấy h đo từ đáy bồn trong khi tâm lỗ cao hơn đáy, sai lệch mốc sẽ xuất hiện ngay ở phép kiểm này.
Độ nhạy: đường kính và C_d có vai trò khác nhau
Giữ mọi đầu vào khác, t tỷ lệ nghịch C_d và nghịch d². Nếu C_d thực chỉ 0,560, thời gian tăng thành 1.421,628 s; nếu là 0,680, thời gian giảm còn 1.170,753 s. Với lỗ giảm từ 50 xuống 40 mm và C_d vẫn 0,620, thời gian thành 2.006,330 s, tăng 56,25%.
Đây là các kịch bản thay đổi từng tham số riêng. Trong thiết bị thật, đổi đường kính hoặc hình dạng mép lỗ có thể làm C_d thay đổi theo. Không dùng quan hệ d⁻² để chứng minh một phương án mới chỉ cần đổi diện tích. Nghiên cứu của Joye và Barrett đối chiếu nhiều cấu hình bồn–ống và chỉ ra tầm quan trọng của việc mô tả đúng tổn thất cùng hình học. [2]
Nếu một thử nghiệm giả định với chính bồn và lỗ này mất 900 s để xuống từ 2,00 tới 0,80 m, đảo công thức cho C_d=0,417237. Con số ấy là hệ số hiệu dụng phù hợp với cặp mức và thời gian đã đo. Nó chưa chứng minh mép lỗ có một tính chất cố định như vậy: sai số mức, một đoạn ống bị bỏ sót hoặc thông khí kém đều có thể bị hấp thụ vào tham số hiệu dụng.
Khi dòng cấp vẫn còn
Xét độc lập trường hợp một dòng cấp không đổi Q_in=0,001000 m³/s tiếp tục vào bồn. Khi đó A_t dh/dt=Q_in−k√h. Mức cân bằng theo mô hình bằng h_eq=(Q_in/k)²=0,0343920 m. Từ mức cao hơn, bồn tiến về mức ấy; nó không thể tự hạ tới 0,020 m khi dòng cấp giữ nguyên.
Đặt b=Q_in/k, có đơn vị √m. Với h_i>h_f>h_eq, tích phân cho:
t_f=(2A_t/k){√h_i−√h_f+b ln[(√h_i−b)/(√h_f−b)]}.
Tỷ số trong logarithm không có đơn vị. Với hai mức 2,40 và 0,20 m, thời gian thành 1.640,721 s. Khi h_f tiến sát h_eq từ phía trên, mẫu bên trong logarithm tiến về không và thời gian tăng không giới hạn. Đây là khác biệt vật lý giữa “còn một ít dòng vào” với việc chỉ cộng một lượng thể tích cố định vào lúc đầu.
Công thức này không còn dùng được nếu Q_in thay đổi theo thời gian hoặc theo mức. Khi đó cần giải phương trình vi phân với lịch sử dòng vào. Có thể kiểm lời giải số bằng cân bằng tích lũy: thể tích giảm bằng tích phân dòng ra trừ tích phân dòng vào.
Hình dáng bồn thay đổi tích phân
Nếu tiết diện ngang phụ thuộc mức, phương trình tổng quát là A_t(h)dh/dt=−k√h. Với một bồn hình nón lý tưởng, đỉnh dưới tại h=0, bán kính miệng R=1,00 m và chiều cao H=2,40 m, A_t(h)=πR²h²/H². Bỏ qua hình học rất nhỏ quanh lỗ, thời gian giữa hai mức là:
t_f=2πR²(h_i^(5/2)−h_f^(5/2))/(5H²k).
Lấy h_i=2,40 và h_f=0,20 m, kết quả 360,307 s. Không diễn giải rằng đáy nón luôn xả nhanh hơn bồn trụ cùng lượng chứa: ở hai mức này bồn nón chỉ xả 2,511820 m³, ít hơn đáng kể 6,911504 m³ của ví dụ trụ. So sánh thời gian cần đi kèm thể tích và cơ sở hình học.
Giới hạn gần cuối quá trình
Khi mức xuống gần lỗ, cột áp tại tâm không còn đại diện tốt cho toàn diện tích nếu một phần lỗ mất ngập. Autodesk mô tả việc chuyển mô hình sang hành vi tràn khi mặt nước thấp hơn đỉnh lỗ. [3] Vì vậy, thời gian tới mức nhỏ nên được tách khỏi phần đường cong đã dùng công thức lỗ ngập.
Nếu đầu ra chìm trong một bể khác, phải dùng chênh cột áp hữu hiệu và xét sự đổi mức bên nhận. Nếu nối ống, cần mô hình ma sát, phụ kiện, cao độ đầu ống và chế độ dòng phù hợp. Bài học thực dụng là xác định đúng hai biên áp suất, mốc h và diện tích chứa trước khi chọn công thức. Khi dữ liệu đo được dùng để hiệu chỉnh C_d, hãy giữ lại đường cong mức–thời gian để kiểm xem một tham số duy nhất có giải thích được toàn khoảng xả hay không.
Nguồn tham khảo
[1] Orifices Experiment — Theory — Virtual Labs, International Institute of Information Technology Hyderabad. Truy cập ngày 09/11/2026.
[2] The Tank Drainage Problem Revisited: Do These Equations Actually Work? — Donald D. Joye và Branden C. Barrett; The Canadian Journal of Chemical Engineering 81, October 2003; bản lưu Đại học São Paulo. Truy cập ngày 09/11/2026.
[3] Orifice — InfoDrainage User Guide — Autodesk. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
