Một bài toán cơ học chất lưu: Tính thời gian chất lỏng chảy khỏi bình qua lỗ nhỏ

Trong cơ học chất lưu, bài toán xác định thời gian chất lỏng chảy khỏi bình qua lỗ nhỏ là một ví dụ kinh điển, xuất hiện từ các giáo trình vật lý phổ thông cho đến các tài liệu chuyên sâu về thủy động lực học. Mặc dù mô hình đơn giản nhất sử dụng phương trình Torricelli cho phép xác định nhanh tốc độ chảy của chất lỏng, nhưng mô hình này chỉ đúng khi chất lỏng được coi là lý tưởng và bỏ qua mọi lực cản.

Đối với chất lỏng thực, dòng chảy luôn chịu ảnh hưởng của độ nhớt và các tổn thất năng lượng do ma sát. Khi lỗ thoát có chiều dài đáng kể hoặc chất lỏng chuyển động trong một ống hẹp, việc sử dụng phương trình Torricelli có thể dẫn đến sai số đáng kể trong việc tính thời gian chảy của chất lỏng. Khi đó, cần sử dụng các mô hình đầy đủ hơn như công thức Darcy–Weisbach hoặc định luật Poiseuille để mô tả chính xác hơn chuyển động của dòng chảy.

Bài viết này trình bày lời giải cho một bài toán cơ học chất lưu quen thuộc bằng ba mô hình có mức độ gần đúng khác nhau. Qua đó, người đọc có thể thấy rõ phạm vi áp dụng của từng mô hình cũng như ảnh hưởng của độ nhớt và lực cản đến thời gian nước chảy khỏi bình.

Bài toán

Một bình chứa nước có tiết diện ngang không đổi S, trên thành bên gần đáy có một lỗ nhỏ để nước chảy ra ngoài. Ban đầu mực nước trong bình cao h₁ so với tâm lỗ. Sau một khoảng thời gian, mực nước giảm xuống còn h₂ (với h₁ > h₂).

Hãy xác định khoảng thời gian để mực nước giảm từ h₁ xuống h₂ trong các trường hợp sau:

1. Xem nước là chất lỏng lý tưởng, bỏ qua mọi tổn thất năng lượng. Tốc độ nước chảy qua lỗ được xác định theo phương trình Torricelli.
2. Lỗ thoát có dạng một ống trụ ngắn dài L, đường kính D. Tính đến tổn thất cột áp do ma sát trong ống bằng công thức Darcy–Weisbach với hệ số cản thủy lực f.
3. Giả sử nước chảy tầng trong ống ngắn nói trên. Biết nước có độ nhớt động lực học η. Hãy xác định thời gian chảy bằng cách sử dụng định luật Poiseuille.

So sánh các kết quả thu được và cho biết phạm vi áp dụng của từng mô hình.

I. Trường hợp lý tưởng: Phương trình Torricelli

Phương trình Torricelli được suy ra từ phương trình Bernoulli cho chất lỏng lý tưởng: $$ \frac{p}{\rho}+\frac{v^{2}}{2}+gh=\text{const} $$ Trong đó $h$ là chiều cao của cột chất lỏng phía trên lỗ, $p$ là áp suất, $\rho$ là khối lượng riêng của chất lỏng. Nếu lớp trên cùng đứng yên và áp suất đồng đều theo độ cao, ta có: $$ v=\sqrt{2gh} $$ Đây chính là phương trình Torricelli.

Giả sử $S$ là diện tích bề mặt chất lỏng trong bình, $s$ là diện tích lỗ, $h$ là chiều cao hiện tại của cột chất lỏng (Hình 1), $h_1$ là chiều cao ban đầu của cột chất lỏng, $h_2$ là chiều cao lúc sau. Từ điều kiện liên tục của dòng chảy, ta có: $$ S\frac{dh}{dt}=-s\sqrt{2gh} $$ hoặc: $$ dt=-\frac{S}{s\sqrt{2g}}\frac{dh}{\sqrt{h}} $$ Lấy tích phân, ta nhận được thời gian thay đổi chiều cao của cột chất lỏng: $$ t_{1}=\frac{S}{s}\sqrt{\frac{2}{g}}(\sqrt{h_{1}}-\sqrt{h_{2}}) $$

Hình 1: Sơ đồ bình chứa chất lỏng có diện tích mặt thoáng $S$ và diện tích lỗ $s$ ở độ cao $h$.

II. Tính đến tổn thất: Phương trình Bernoulli với tổn thất cột áp

Trong thực tế, dòng chảy bị cản trở bởi ma sát trong. Khi đó, trong phương trình Bernoulli sẽ xuất hiện thêm một số hạng: $$ \frac{p}{\rho}+\frac{v^{2}}{2}+gh=h_{\text{tt}} + \text{const} $$ Trong trường hợp chất lỏng chảy qua một ống ngang, phần tổn thất cột áp được xác định bởi công thức thực nghiệm Darcy-Weisbach: $$ h_{\text{tt}}=f\frac{L}{D}\frac{v^{2}}{2g} $$ Ở đây $f$ là hệ số cản thủy lực, $L$ là chiều dài ống, $D$ là đường kính của ống (Hình 2). Hệ số $f$ phụ thuộc vào nhiều nguyên nhân, và trước hết là phụ thuộc vào chế độ chảy: chảy tầng hay chảy rối.

Hình 2: Chất lỏng chảy qua lỗ có dạng ống ngang với chiều dài $L$ và đường kính $D$.

Ví dụ, đối với dòng chảy tầng trong các ống có thành cứng, hệ số này có thể được xác định thông qua số Reynolds $Re$: $$ f=\frac{64}{Re} $$ Do đó, tốc độ chảy sẽ bị giảm xuống: $$ v=\sqrt{\frac{2gh}{1+\frac{fL}{D}}} $$ Trong trường hợp của chúng ta, thời gian thay đổi chiều cao của cột chất lỏng trong bình sẽ bằng: $$ t_{2}=\frac{S}{s}\sqrt{2\frac{D+fL}{gD}}(\sqrt{h_{1}}-\sqrt{h_{2}}) $$ Phương pháp này mang lại kết quả thực tế hơn đối với các tính toán kỹ thuật, vì nó đã tính đến lực cản trong ống và ảnh hưởng từ hình dạng hình học của ống.

III. Độ nhớt: Biên dạng vận tốc Poiseuille

Đối với dòng chảy tầng của chất lỏng bên trong một ống hẹp (Hình 3), vận tốc được phân bố dọc theo bán kính theo một quy luật parabol, được Poiseuille thiết lập thông qua thực nghiệm: $$ v(r)=\frac{\Delta p}{4\eta L}(R^{2}-r^{2}) $$ Ở đây $\Delta p = \rho gh$ là độ chênh lệch áp suất, $\eta$ là độ nhớt động lực học của chất lỏng, $L$ là chiều dài ống, $R$ là bán kính của ống.

Hình 3: Biên dạng phân bố vận tốc của dòng chảy tầng bên trong ống hẹp theo định luật Poiseuille.

Bằng cách lấy tích phân sự phân bố vận tốc trên toàn bộ tiết diện, ta thu được lưu lượng thể tích của chất lỏng: $$ Q=\frac{\pi\Delta p R^{4}}{8\eta L} $$ Sau khi thay $\Delta p = \rho gh$, phương trình mô tả sự thay đổi chiều cao của chất lỏng theo thời gian có dạng: $$ S\frac{dh}{dt}=-\frac{\pi\rho g R^{4}h}{8\eta L} $$ Từ đó suy ra: $$ dt=-\frac{8\eta LS}{\pi\rho g R^{4}}\frac{dh}{h} $$ Lấy tích phân theo chiều cao từ $h_1$ đến $h_2$, ta tìm được thời gian chảy: $$ t_{3}=\frac{8\eta LS}{\pi\rho g R^{4}}\ln\frac{h_{1}}{h_{2}} $$ Cách tiếp cận này mô tả chính xác nhất sự chảy tầng chậm, đặc biệt trong trường hợp các chất lỏng có độ nhớt, khi mà lực cản đối với dòng chảy bên trong lỗ đóng vai trò chi phối các hiệu ứng khác.

IV. Kết luận

Nghiên cứu đã chỉ ra rằng thời gian chảy của chất lỏng phụ thuộc rất nhiều vào mô hình mô tả được lựa chọn. Phương trình Torricelli, tuy đơn giản và trực quan, nhưng chỉ có thể áp dụng cho chất lỏng lý tưởng không có độ nhớt và không mất mát năng lượng. Trong điều kiện thực tế, nó cho thời gian chảy ngắn hơn so với thực tế.

Việc tính toán tổn thất theo phương trình Darcy-Weisbach cho phép mô tả tình huống một cách thực tế hơn, đặc biệt nếu lỗ có chiều dài và lực cản đáng kể. Phương pháp này rất tiện lợi khi cần một sự xấp xỉ nhanh có tính đến các yếu tố cản trở chính.

Định luật Poiseuille mang lại một cái nhìn toàn diện nhất về dòng chảy tầng, chỉ ra rằng vận tốc trên tiết diện được phân bố theo quy luật parabol. Cách tiếp cận này đóng vai trò quan trọng trong các bài toán mà độ nhớt là yếu tố then chốt — ví dụ như khi chất lỏng chảy chậm trong các kênh hẹp.

Tóm lại, có thể nói rằng mỗi mô hình đều có phạm vi áp dụng riêng của nó, và việc lựa chọn đúng đắn phụ thuộc vào mức độ chi tiết cần thiết để mô tả một quá trình vật lý.

📚 Bạn nên đọc

• Lại bàn về đường cong đoản thời (Brachistochrone)
  Bài toán kinh điển của Johann Bernoulli về thời gian chuyển động nhỏ nhất, liên quan đến tư duy mô hình hóa trong cơ học và nguyên lý biến phân.
• Các chuyên đề Vật lý
  Hệ thống các chuyên đề Vật lý từ cơ bản đến nâng cao dành cho học sinh yêu thích Vật lý và giáo viên giảng dạy.
• Các bài toán trộn khí lý tưởng
  Một ví dụ về phương pháp xây dựng mô hình vật lý và giải quyết bài toán bằng các định luật cơ bản.